Abstract
To predict the refrigeration capacity of a two-stage pulse tube cryocooler operating at less than 4K, a new numerical simulation is developed. Full one-dimensional governing equations for expander are solved using the finite difference method. It considers variable physical properties of real helium gas and regenerator materials. The calculation results show the pressure oscillation in one refrigeration cycle, the equivalent PV diagrams both at the hot end and at the cold end of the pulse tube, and the axial temperature distributions both of the pulse tube and the regenerator.