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.
[1]
A. Shapiro.
The dynamics and thermodynamics of compressible fluid flow.
,
1953
.
[2]
Yoichi Matsubara,et al.
A numerical method of regenerator
,
2004
.
[3]
Y. Matsubara,et al.
Phase Shift Effect of the Long Neck Tube for the Pulse Tube Refrigerator
,
1997
.
[4]
D. E. Gray,et al.
American Institute of Physics Handbook
,
1957
.
[5]
S. Patankar.
Numerical Heat Transfer and Fluid Flow
,
2018,
Lecture Notes in Mechanical Engineering.
[6]
Gregory W. Swift,et al.
Analysis and performance of a large thermoacoustic engine
,
1992
.