Presented at the 16th International Cryocooler Conference, held May 17-20, 2008 in Atlanta, Georgia.A phasor diagram for a pulse tube refrigerator (PTR) is a vectorial representation of mass flow rate, pressure, and temperature at different locations as a function of time. With the help of a phasor diagram, the operation of different types of pulse tube refrigerators can be well understood. Phasor analysis based on these diagrams gives an idea regarding the underlying complex phenomena of the PTR. In the present work, a simplified model has been presented based on the assumption that there is no phase difference between temperature and pressure throughout the working space. The phasor analysis is extended to a two-stage Orifice Pulse Tube Refrigerator (OPTR) and to a Double Inlet PTR (DIPTR). The important contribution of the work is that it highlights the condition for which the DIPTR will work better than the OPTR

Atrey, M. D.

Mohanta. L.

Scholarly Materials And Research @ Georgia Tech

Phasor Analysis of Pulse Tube Refrigerator L. Mohanta, M.D. Atrey Department of Mechanical Engineering Indian Institute of Technology Bombay, Powai Mumbai-400076, India ABSTRACT A phasor diagram for a pulse tube refrigerator (PTR) is a vectorial representation of mass flow rate, pressure, and temperature at different locations as a function of time. With the help of a phasor diagram, the operation of different types of pulse tube refrigerators can be well understood.  Phasor analysis based on these diagrams gives an idea regarding the underlying complex phenomena of the PTR. In the present work, a simplified model has been presented based on the assumption that there is no phase difference between temperature and pressure throughout the working space. The phasor analysis is extended to a two-stage Orifice Pulse Tube Refrigerator (OPTR) and to a Double Inlet PTR (DIPTR).  The important contribution of the work is that it highlights the condition for which the DIPTR will work better than the OPTR.INTRODUCTIONIn recent years, pulse tube refrigerator (PTR) technology has undergone a lot of improvement. PTRs are now widely used for various applications, e.g. MRI, NMR, space technology, cooling of SQUIDS etc. This is due to the absence of cold moving part in PTRS, which improves their vibration characteristics and reliability. The PTR works in a cyclic manner which consists of a) Compression, b) Constant Pressure Regenerative cooling, c) Expansion and d) Constant pressure regenerative heating [1]. Different analyses have been presented in order to analyse the underlying phenomena taking place in the PTR. Various numerical models have been proposed that give quite good result [2, 3].  However, for a DIPTR or for a two-stage PTR, these analyses become more mathematical in nature and therefore are not easy to understand.  Also, they do not help us to imagine the physical processes that are taking place in the PTR. In view of this, the present work highlights the working of different types of PTR using phasor analysis, which is a vectorial representation of the working of the PTR.   Basic analysis of the PTR using phasor diagrams was introduced by Radebaugh [4], while Kittel et al. [5] carried out similar analyses using Taylor series approximation.  Hofmann and Pan [1] explained the working of different types of PTR with phasor diagrams and analyzed them using an electrical analogy and complex geometry.  A phasor diagram for a Double Inlet PTR was suggested by Chokhawala et al. [6].  In the present work, a detailed phasor analysis for a two-stage PTR, both OPTR and DIPTR, is presented based on the law of conservation of mass. The phasor analysis helps to understand the importance of phase difference between mass flow rate at the cold end and the pressure pulse in the pulse tube. The refrigerating effect for different types of PTR strongly depends on the 299300 PULSE TUBE ANALYSIS AND EXPERIMENTAL MEASUREMENTSphase shift arrangement and also on the phase difference. Due to the complexities of the flowinvolved in the two-stage DIPTR and OPTR, the phasor diagram put forth in the present work will prove to be very helpful to understand the importance of these parameters. The present work also introduces the concept and importance of phase lead of the double inlet flow entry in the DIPTR. The analysis also comes out with a necessary condition on the value of the phase leadangle in order to make a DIPTR function better in terms of refrigerating effect, as compared to an OPTR ANALYSIS FOR OPTR Assumptionsa) The gas obeys the ideal gas law. b) There is no phase difference between the pressure and the mass flow rate at the hot endthrough the orifice of the pulse tube. c) There is no phase difference between the pressure and the temperature throughout the working space of the PTR. d) The variation of pressure, temperature and mass flow at a given location is sinusoidal. e) The pulse tube behaves adiabatically during the operation. Pressure and Temperature Variations In the pulse tube, pressure and temperature variations are given as: 0 1 cos( )P P P t  (1) 0 1 cos( )T T T t  (2) Applying the adiabatic condition to the pulse tube,10 0T PT P (3) Using Taylor series expansion and neglecting higher order terms [5],  10 01TT P1P  (4) Mass low in the Pulse Tube  Considering the whole pulse tube as a control volume and applying the law of conservation of mass, ptcd Vv vdt h (5) Assuming that the pulse tube is at a mean temperature, density can be taken asmptPRTptcmptd PVm mRT dth  (6) ptcmptV dPm mRT dth  (7) Mass flow rate at the hot end, , through orifice is directly proportional to pressure difference [5].  hm1 cosP P t  (8) Mass flow rate at hot end is given as [5], F1 1 coshhcTm C PTt (9)Combining Eq. (9) and Eq. (7), 1 1 1cos cos( / 2)pthcc mptVTm C P t P tT RT (10) From Eq (10), it is clear that mass flow rate at the cold end consists of two components. The first term in the right hand side is the contribution resulting from the mass flow rate at the hot end through the orifice which is in phase with the pressure. The second term represents the additional amount of mass entering the pulse tube due to pressurization, which is leading the earlier term by 90o.  Fig.1 shows the equivalent phasor representation of Eq. (10). The phase angle ï represents the phase difference between the pressure pulse and T cm T. The numerical analysis of the first law of thermodynamics applied to the PTR also shows a similar phase difference in the mass flow rate at different locations of the pulse tube [7]. Enthalpy flow Due to the pulsating nature of the flow in the pulse tube, time-averaged enthalpy is considered for energy balance for a cycle.  Fig. 2 shows enthalpy flow in the pulse tube. Applying the first law of thermodynamics to the cold end heat exchanger of the PTR as shown in Fig. 2 gives, Q rH H  (11) where Q is the refrigerating effect,  is the time-averaged enthalpy from the cold end heat exchanger of the pulse tube and  is the time-averaged enthalpy flow from regenerator to cold end heat exchanger. HrHFor perfect regeneration, =0r     HHence the refrigerating effect obtained will be, 0QpcCH m Tdt  (12) Substituting the expression for mass flow rate from Eq. (10) and using Eq. (3),  10( )( cos )pc cCH m T T t dt  (13) Figure 1.  Phasor diagram showing mass flow rate at hot end and cold end of pulse tube. 301PHASOR ANALYSIS OF PULSE TUBE REFRIGERATORFigure 2.   Enthalpy flow in the pulse tube. The integration of the second term in Eq. (10) becomes zero as this term is 900 out of phase. This shows there is no contribution to the refrigerating effect by the mass flow due to pressurization.  Hence Eq. (13) reduces to1 1 10cos * ( cos )p hccC TH C P t T TTt dt (14)On further integration, 21105phC PHPT C (15)This can be expressed in terms of mass flow rate at the cold end, .  From Fig.1 m  and can be related as follows. cm hcmcoshcmm And 101cos2c cPQ H RT mP (16) This result is in agreement with Radebaugh [4] and Kittel [5].  Eqs. (15), (16) show that thesmaller is the phase difference, ï , and larger is the mass flow rate, c , and the more refrigerating effect will be obtained.  In other words, the product mmcosc  has to be as large as possible. In the case of BPTR, the mass flow rate at the hot end ( h ) is zero, which implies ï is equalto 90m0.  Hence, theoretically, the refrigerating effect should be equal to zero for this configu-ration. However, due to heat transfer from gas to the wall of the pulse tube, ï decreases slightlybelow 900.  The numerical value of ï could also be indicative of the losses in the system. As the loss increases, ï  increases and that will lead to decrease in the refrigerating effect.PHASOR DIAGRAM FOR TWO STAGE OPTRThe phasor diagram for a two-stage OPTR can be drawn in a similar fashion as that of the single-stage OPTR. For the two-stage PTR, both the pulse tubes are considered separately. Fig. 3 shows the phasor diagram for both stages of the pulse tube. Fig. 3(a) gives the phasor diagram for the second stage pulse tube of the two-stage OPTR where  and  represent the phasors for mass flow rate at the second stage hot and cold ends, respectively.  It can be seen 2hm 2cmFigure 3.  Phasor diagram for pulse tube (a) Second stage (b) First stage.302 PULSE TUBE ANALYSIS AND EXPERIMENTAL MEASUREMENTSfrom Fig. 3(a) that theRe 2g hm , phasor for mass flow rate at the hot end of the second stage regenerator, can be obtained by adding the phasor for regenerator-2 loss, (Re 21Re 2gm gVPR T), to the phasor for mass flow rate at the cold end of the pulse tube-2,    Fig. 3(b) gives the phasor diagram for the first stage of the two stage OPTR.  The flow coming out of regenerator-1 gets divided into two streams; one stream goes to the first stage pulse tube and the second one goes through regenerator-2. Therefore, the mass flow rate through the cold end of the regenerator-1 (2cmRe 1g cm ) can be obtained by adding the phasors  and 1cmRe 2g hm as given in Eq. (17). R e 1 R e 2 1g c g h cm m m (17)Redefining  as given in Eq. (10) for the first stage, 1cm11 1 1 11 1cos( ) cos( / 2)pthcc mptVTm C P t P tT RT (18) As shown in Fig. (3),Re 2g hm  can be expressed as, Re 22 Re 21 1 1 12 2 Re 2cos cos( / 2) ( / 2)g hpt ghc mpt m gmV VTC P t P t P tT RT R T (19) Fig. 4 shows the combined phasor diagram for the two-stage OPTR including the losses inthe compressor. The effect of void volume in the heat exchangers is neglected. The refrigeratingeffect (Q) can be calculated as earlier using Eq. (15) or (16) separately for both the stages.ANALYSIS FOR DIPTRThe analysis for the OPTR provides the background for the analysis for the DIPTR. A DIPTR differs from an OPTR by a connection from the high pressure line to the hot end of the pulse tube as shown in Fig. 5. All the assumptions are the same as those for the OPTR. In addition, a phase lead angle,  is considered between the flow through double inlet valve and theFigure 4.  Phasor diagram corresponding to the two-stage OPTR. pressure pulse in the PTR. 303PHASOR ANALYSIS OF PULSE TUBE REFRIGERATOR304 PULSE TUBE ANALYSIS AND EXPERIMENTAL MEASUREMENTSFigure 5.  Schematic diagram of Double Inlet Pulse Tube Refrigerator (DIPTR).Mass Flow in the Pulse Tube (DIPTR) Considering the whole pulse tube as a control volume and applying the law of conservation of mass,  ptc od Vv v vdt DI (20) ptc o Dmptd PVm m mRT dtII (21) h o Dm m m  (22) where is the net mass flow rate leaving the hot end of the pulse tube. hmptcmptd PVm mRT dth  (23) Thus, o  and m DIm are directly proportional to the pressure difference across the valve. As per the assumption, DI  leads the pressure by phase angle . Substituting  in Eq. (23) in terms of pressure we get, m hm1 1 1 2 1cos cos( / 2) cos( )pth hcc mpt cVT Tm C P t P t C P tT RT T (24) Figs. 6 (a) and (b) show the equivalent representation of the expression given in Eq. (24) in terms of phasors.  It can be seen from these figures that h  depends on o andm m DIm  which are controlled by the respective valve openings. Fig. 6 (a) shows the phasor diagram when the DI phase lead angle  is less than 900, while Fig. 6(b) shows the phasor diagram when  is greaterthan 900.  It can be seen from these figures that coscm  is more than when  is more than 90°.omEnthalpy low  Substituting the mass flow rate expression from Eq. (24) for calculating refrigerating effect,  1 1 1 2 1 10 0cos *( cos ) cos( )*( cos )p h hc cc cC T TH C P t T T t dt C P t T TT Tt dt (25)Figure 6.  Phasor diagram for DIPTR when phase lead (a) less than 90o (b) more than 90oF2 21 11 20 0cos5 5p pDIPTR h hC CP PQ H T C T CP P (26) Eq. (26) shows that the refrigerating effect for the DIPTR depends on . The first term is the same as that for the OPTR, while the second term represents the contribution due to the flow through the DI valve.From Fig. 6 we can rewrite the Eq. (26) as, 11 cos2DIPTR c coPQ RT mP(27)Eq. (27) indicates the fact that for maximum QDIPTR, cos  has to be as large as possible. cmPHASOR DIAGRAM FOR TWO STAGE DIPTRThe phasor diagram can be drawn for a two-stage DIPTR by considering two pulse tubes separately and finally adding them together as done in the earlier case for the two-stage OPTR. Fig. 7 shows the phasor diagram for the two-stage DIPTR.  The mass flow rate throughregenerator-1 (Re 1g cm ) is divided through pulse tube-1 and regenerator-2 which is given below,Re 1 Re 2 1g c g hm m mc(28)The mass flow rate at the cold end of the pulse tube-1 m can be redefined from Eq (24),1cwhileRe 2g hm is given as, 2 12Re 22 R1 1 1 12 2 Re 2cos( )cos cos( / 2) ( / 2)hcg hpt ghc mpt m gTm C P tTV VTC P t P t P tT RT R Te 2 (29)The refrigerating effect for both the stages can be calculated separately using Eq. (26). Figure 7.  Phasor diagram corresponding to the two-stage DIPTR. 305PHASOR ANALYSIS OF PULSE TUBE REFRIGERATORRESULTS AND DISCUSSIONDouble Inlet Phase Lead Angle, From Eq. (26) and the phasor diagrams given in Fig.6, it can be seen that QDIPTR depends on the value of  when compared with QOPTR .  It can be observed from Eq. (26) that  has to be more than 90° to achieve higher QDIPTR compared to QOPTR, for the same opening of orifice. The value of  can be varied in practice by varying the opening of the double inlet valve. From Fig. 6 it can be noted that the projection of c on the pressure phasor, that ism coscm , increases withincrease in . When the double-inlet valve is closed, DI is equal to zero, and the expression for the QDIPTR given by Eq. (26) reduces to Eq. (15), which corresponds to an OPTR. As the valueof  increase, the value of  decreases, and the refrigerating effect increases. This fact isexplained in detail in the next section.mPhasor Analysis for an Existing Single-Stage PTR The above analysis is now applied to an existing PTR [8]. The flow coefficient for orificeand DIPTR valves are taken as given below. C1=6E-10kg/sec.pa [1, 9] C2=1E-10kg/sec.pa [1, 9] Fig. 8 shows the variation of refrigerating effect (Q) with pressure ratio for differentconfigurations of the PTR. Based on arguments made above, it can be seen that QDIPTR is less as compared to QOPTR when  is 45º (less than 90º) and more when it is 135º (more than 90º). Also, it can be seen from these figures that an increase in pressure ratio increases the refrigeratingeffect.The variation of h  and c with  is shown in the Fig. 9. At  equal to zero, although cm isquite large, Qm mDIPTR is minimum due to the fact that  is also large for this case. It can be observed from Eq. (27) that, QDIPTR depends on the product of c and com s . This can also be confirmedfrom Fig. 6 that  is maximum for  equal to zero, amounting to minimum QDIPTR.Fig.10 shows the variation of QDIPTR and DIPTR with  and QOPTR is plotted on the graph for comparison.  As  increases from zero onwards, c  decreases up to a certain value, but mcoscm  increases continuously due to a decrease in . As a result, the refrigerating effect QDIPTRincreases. At  equal to 90°, the value of coscm  for the DIPTR is the same with or without the presence of DIm . In this case, QDIPTR with DIm  is the same as QOPTR ( DIm =0). This can be seen Figure 8.  Variation of refrigerating effect with pressure ratio for different double inlet entry. 306 PULSE TUBE ANALYSIS AND EXPERIMENTAL MEASUREMENTSfrom Fig. 6. The presence of DI  results in less gas passing through the regenerator, which causes c to be smaller and formm DIPTR to be reduced.  In the absence of DI , the value of c increases along with an increase in value of , which causes more losses in reality.  Fig. 9 shows that, as the value of  increases in DIPTR, the phase difference  reduces up to a certain minimum value and then increases. As discussed earlier, QmmDIPTR will be more comparable to QOPTR when  is greater than 90°.  From Figs. 9 and 10, it can be seen that QDIPTR is maximumwhen  is 180°. CONCLUSION  The present work highlights the fundamental relationship between mass flow rates at the hot and cold ends of the pulse tube and also the mass flow rate through the orifice and DI valve.  Mass flow rate and enthalpy are expressed in terms of pressure, temperature, and harmonic terms. The work concludes that the phase difference,  and , are important parameters in estimating mass flow rates and refrigerating effect. The phasor diagrams have been presented for two-stage OPTR and DIPTR from which the concept of phase difference  and DI phase lead angle  can be very well understood.  It has been found that the refrigerating effect is directly proportional to the mass flow rate at the cold end and the amplitude of dynamic pressure. Also, Figure 9.  Variation of mass flow at hot and cold end of the pulse tube with phase lead angle. Figure 10. Variation of phase angle and refrigerating effect of DIPTR with phase lead angle. 307PHASOR ANALYSIS OF PULSE TUBE REFRIGERATORfor the  DIPTR, the refrigerating effect will be more than for the OPTR when the DI phase lead angle  is greater than 90º.   The phase diagrams presented for two-stage OPTR and DIPTR systems provide an insight into the physical phenomena taking place in these PTRs. REFERENCES 1.  Hofmann A, Pan H, “Phase shifting in pulse tube refrigerators,” Cryogenics (1999), 39 (6), pp. 529-537.2.  Zhu S, Matsubara Y, “Numerical method of inertance tube pulse tube refrigerator,” Cryogenics(2004), 44 (9), pp. 649-660. 3.  Wang C, Wu PY and Chen ZQ, “Numerical analysis of double-inlet pulse tube refrigerator,” Cryogenics (1993), 33 (5), pp. 526-530. 4.  Radebaugh R. “A Review of pulse tube refrigeration,” Advances in Cryogenic Engineering, 35 (1990), pp. 1191-1205. 5.  Kittel P, Kashani A, Lee J M and Roach P R, “General Pulse Tube Theory,” Cryogenics (1996) 36 (10), pp. 849-857.  6.  Chokhawala M D, Desai K P, Barik H B, Narayankhedkar KG, Phasor analysis for double inlet pulse tube cryocooler,” Advances in Cryogenic Engineering , 45 (2000), pp. 159-165. 7.  He Y L, Huang J, Zhao C F, Liu Y W, “First and second law analysis of pulse tube refrigerator,” Applied Thermal Engineering, 26 (2006), pp. 2301–2307. 8.  Gawli BS, Design and development of miniature pulse tube nitrogen liquefier, PhD thesis, IIT Bombay, India. 9.  Boer PCT, “Performance of the inertance pulse tube,” Cryogenics (2002), 42 (3-4), pp. 209-221.  308 PULSE TUBE ANALYSIS AND EXPERIMENTAL MEASUREMENTS

Phasor Analysis of Pulse Tube Refrigerator

https://smartech.gatech.edu/bitstream/1853/38790/1/038.pdf

Phasor Analysis of Pulse Tube Refrigerator

Abstract

Similar works

Full text

Available Versions

Scholarly Materials And Research @ Georgia Tech