257 1 Numerical Simulation of the Effect of Heat Conductivity in a 4 K Regenerator M. Xu, Q. Bao, A. Tsuchiya Technology Research Center, Sumitomo Heavy Industries, Ltd. Nishitokyo-city, Tokyo, Japan 188-8585 ABSTRACT In a 4K cryocooler, the performance of the second stage is strongly dependent on the efficiency of the second stage regenerator. In order to improve the efficiency of the second stage regenerator, the effect of the heat conductivity of the regenerator material is analyzed by numerical simulation. The temperature profile in the second stage regenerator with an artificially increased heat conductivity of the regenerator material at the warm end is analyzed. Based on the simulation results, a novel and simple way, which is to replace the regenerator material at the warm end with a material having a higher heat conductivity, is proposed. By increasing the heat conductivity to 1x10 6 times that of lead, the P-V power is increased by 0.7 W at the first stage and 0.15 W at the second stage. Correspondingly, the cooling capacity is increased by 0.4 W at the first stage and 0.03 W at the second stage. The simulation results are reported in this paper. INTRODUCTION In a 4K cryocooler, the performance of the second stage is strongly dependent on the real gas effect and the efficiency of the second stage regenerator. The real gas effect can be estimated by the method suggested by de Waele and Xu et al. 1,2 Then, the regenerator loss at the second stage can be estimated by the enthalpy flow excluding the real gas effect. 3,4 In order to improve the efficiency of a 4K pulse tube cryocooler, the extra enthalpy flow caused by the real gas effect was utilized by adding a heat exchanger to the second stage regenerator. 5 However, it is difficult to extract cooling from the second stage of a 4K GM cryocooler because there is a large thermal conducting resistance caused by the clearance between the cylinder and the displacer, and the low thermal conducting displacer wall. Xu, et al. reported that the performance of a 4K GM cryocooler can be improved by shifting the temperature profile at the warm end of the second stage regenerator to a higher level. 6 A novel and simple way to shift the temperature profile, which is to replace the regenerator material at the warm end of the second stage regenerator with a material having a lower heat capacity, was proposed and the effect was confirmed with a simulation method. 6 In this paper, the temperature profile of the second stage regenerator with respect to the heat conductivity of the regenerator material and the mass variation in the void volume of the second stage regenerator are analyzed. Another novel and simple way for changing the temperature profile in the second stage regenerator is proposed. The effect of the novel way is analyzed by the simulation method we developed. 4 Cryocoolers 18, edited by S.D. Miller and R.G. Ross, Jr. International Cryocooler Conference, Inc., Boulder, CO, 2014 277
258 Regenerator & Recuperator Investigations 278 REGENERATOR & RECUPERATOR INVESTIGATIONS 2 Figure 1. Schematic diagram of a second stage regenerator SIMULATION We reported a simulation method and confirmed the accuracy of the simulation results by comparing them with the results measured in a 1W 4K GM cryocooler. 4 The simulation method has been verified to be an effective way to improve the efficiency of a 4K GM cryocooler. 3-4, 6 Figure 1 shows a schematic of the second stage regenerator, which is divided into three parts. Part I is at the warm end, Part II is at the middle, and Part III is at the cold end. The heat conductivities of the regenerator materials in the three parts are called ê w,1, ê w,2 and ê w,3, respectively. During the simulation, the heat conductivity of the regenerator material is assumed to be equal to the normal value of lead in Part II, while the heat conductivity in Part I is assumed to be equal to or larger than that of lead. For the material in Part I, all other properties, such as density and heat capacity, are assumed to be the same as those of lead. In the simulation model, the two-stage G-M cryocooler is divided into many elements and the state in each element can be calculated by solving the basic equations. 6 The energy conservation equations of gas flow and regenerator material are listed as follows: The energy conservation equation of gas flow is, where Ù is the heat transfer coefficient and A is the heat transfer area. The energy conservation equation of the regenerator material is, where V w is the volume of material, c w is the volumetric specific heat capacity of the material, A w is the cross-section area of the material, and ê w is the heat conductivity of the material. Since the regenerator is filled with spheres, the equivalent heat conductivity of the regenerator packing bed is much smaller than that of a solid material. The equivalent heat conductivity of the regenerator material can be calculated by the following correlation, (1) (2) where ê g is the heat conductivity of helium gas. (3) SIMULATION RESULTS AND DISCUSSIONS Simulations have been performed using the parameters of a RDK-408D2 1W 4K two-stage GM cryocooler manufactured by Sumitomo Heavy Industries, Ltd. (SHI). For typical calculations, the cold-head is operated at 1.0 Hz, the high and low pressures are 2.25 MPa and 0.76 MPa, respectively. The inner diameters of the first and second stage cylinder are 82 mm and 35 mm, respectively. The diameter of the second stage regenerator is 30 mm. The normalized length of Part I, II and III are about one-third for each. The second stage regenerator is filled with lead and HoCu 2 spheres. The porosity of the second stage regenerator is assumed to be 0.3 for all calculations. The average sphere diameter is 0.44 mm for lead and 0.33 mm for HoCu 2.
EFFECT OF HEAT CONDUCTIVITY IN 4K REGENERATOR 279 259 3 Figure 2. Temperature profile in the second stage regenerator. Temperature Profile in the Second Stage Regenerators Figure 2 shows the simulation results of the temperature profile in the second stage regenerator. The heat conductivity of regenerator material in Part I is assumed to be 1x10 2, 1x10 4 and 1x10 6 times that of lead. As shown in Fig. 2, when the heat conductivity in Part I is assumed to be larger than that of lead, the temperatures in the second stage regenerator increase. Accordingly, the average temperature in the second stage regenerator increases. Pressure Variation in Expansion Volumes Figures 3a and 3b show the pressure variation in the first and the second stage expansion volumes, respectively. As shown in the figures, when the heat conductivity in Part I is assumed to be 1x10 6 times that of lead, the pressures in the expansion volumes goes down slightly more slowly after approaching the high pressure and then rises up slightly more slowly after approaching the low pressure. See Reference 6 for more details about the mechanism. As shown in Fig. 2, the tempera- Figure 3a. Pressure variation in the first-stage expansion volume.
260 Regenerator & Recuperator Investigations 280 REGENERATOR & RECUPERATOR INVESTIGATIONS 4 Figure 3b. Pressure variation in the second-stage expansion volume. tures slightly increase in the region from 8 K to 30 K, where the density change caused by pressure change is large. Accordingly, the average temperature in the second stage regenerator increases and a smaller amount of gas is needed to fill the void volume of the second stage regenerator. Therefore, the pressure in the cold head goes down slightly more slowly. Similarly, since the average temperature in the second stage regenerator is increased and a smaller amount of gas is needed to fill the void volume of the second stage regenerator, the pressure in the cold head rises up slightly more slowly. The pressure change caused by an increased heat conductivity is relatively small compared to that caused by a reduced heat capacity reported in Ref. 6. The reason can be analyzed from Equation (3). The equivalent heat conductivity of the regenerator packing bed is strongly dependent on the heat conductivity of the helium gas. In contrast, the impact on the equivalent heat conductivity caused by the heat conductivity of the regenerator material is relatively small. P-V Power and Cooling Capacity Figure 4 shows the P-V diagram of the first and the second stage expansion volumes. As shown in Fig.4, when the heat conductivity in Part I is assumed to be 1x10 6 times that of lead, the areas of the first and second-stage P-V diagrams increase slightly. The reason is that the pressures in the expansion volumes go down slightly more slowly after approaching the high pressure and rise up slightly more slowly after approaching the low pressure. Therefore, the phase shift between the pressure and the expansion volume is improved, and the P-V power is slightly increased. The simulation results of the cooling capacity and the P-V power when the heat conductivity of the regenerator material in Part I of the second stage regenerator is increased are shown in Table 1. The enthalpy flux from the first and the second stage regenerator is also shown in Table 1. As shown in Table 1, when the heat conductivity is increased, the P-V power and the cooling capacity at both stages slightly increase. By increasing the heat conductivity to 1x10 6 times that of lead, the P-V power is increased by 0.7 W at the first stage and 0.15 W at the second stage. The enthalpy flux from the regenerator is increased by 0.1 W at the first stage, but increased 0.12 W at the second stage. Accordingly, the cooling capacity is increased by 0.4 W at the first stage and 0.03 W at the second stage. Apparently, the cooling capacity of the first stage is strongly dependent on the P-V power, while the cooling capacity of the second stage is mainly dependent on the regenerator efficiency. Although the effect is relatively small, the tendency is the same as that with a reduced heat capacity at the warm end of the second stage regenerator.
EFFECT OF HEAT CONDUCTIVITY IN 4K REGENERATOR 281 261 5 Figure 4. P-V diagram of the expansion volumes. (a) first stage; (b) second stage. Table1. Simulation Results with Increased Heat Conductivity in Part I of the Second Stage Regenerator.
262 282 REGENERATOR Regenerator & RECUPERATOR & Recuperator INVESTIGATIONS Investigations CONCLUSIONS A novel and simple way to shift the temperature profile is proposed, that is to replace the regenerator material at the warm end of the second stage regenerator with a material having a higher heat conductivity. According to the simulation results, by increasing the heat conductivity to 1x10 6 times that of lead, the P-V power is increased by 0.7 W at the first stage and 0.15 W at the second stage. The cooling capacity is increased by 0.4 W at the first stage and 0.03 W at the second stage. 6 REFERENCES 1. Xu, M.Y., De Waele, A.T.A.M. and Ju, Y.L., A Pulse Tube Refrigerator below 2K, Cryogenics, vol. 39 (1999), pp. 865-869. 2. De Waele, A.T.A.M., Xu, M.Y. and Ju, Y.L., Nonideal-gas Effect in Regenerators, Cryogenics, vol. 39 (1999), pp. 847-851. 3. Xu, M. and Morie, T., Development of High-efficiency 4K GM Cryocoolers, Proceeding of the 24th International Cryogenic Engineering Conference and International Cryogenic Material Conference, Cryogenics and Superconductivity Society of Japan, Tokyo (2012), pp. 403-406. 4. Xu, M. and Morie, T., Numerical Simulation of 4K GM Cryocooler, Cryocoolers 17, ICC Press, Boulder, CO (2012), pp. 253-259. 5. Zhu, S., Ichikawa, M. and Nogawa, M., et al., 4K Pulse Tube Refrigerator and Excess Cooling Power, Advances in Cryogenic Engineering, Vol 47B, American Institute of Physics, New York (2002), pp. 633-640. 6. Xu, M. and Morie, T., Numerical Simulation of the Second Stage Regenerator in a 4K GM Cryocooler, Advances in Cryogenic Engineering, Vol 59B, American Institute of Physics, New York (2014), pp. 633-640.