Numerical method of inertance tube pulse tube refrigerator

A nodal analysis method for simulating inertance tube pulse tube refrigerators is introduced. The energy equation, continuity equation, momentum equation of gas, energy equation of solid are included in this model. Boundary condition can be easily changed to enable the numerical program calculate thermal acoustic engines, inertance tube pulse tube refrigerators, double inlet pulse tube refrigerators, and others. Implicit control volume method is used to solve these equations. In order to increase the calculation speed, the continuity equation is changed to pressure equation with ideal gas assumption, and merged with momentum equation. Then the algebraic equation group from continuity and momentum equation becomes one group. With this numerical method, an example calculation of a large scale inertance tube pulse tube refrigerator is shown.