I am currently working on 1-D numerical simulation of pulse tube refrigerator. The governing equations are given in attachment.
I have plotted results showing mass flow rate over the cycle (mass flow rate vs pressure for a cycle). As per the cyclic analysis, one should know that mass flow rate should have initial and final state same after a cycle completes or cycle averaged mass flow rate must be zero during that cycle. Any how, I am not getting the same. Please see the results in attachment.
What I am doing is, I am using nozzle equation for orifice mass flow rate calculation which will serve as boundary condition for mass flow rate. Then I will use 1-D continuity equation to get mass flow rate at other faces. (See the attachment)
Please let me know where I am making mistakes.