 Open Access
 Total Downloads : 0
 Authors : Pankaj Kumar, Ajay Kumar Gupta, M. Kumar, R. K. Sahoo
 Paper ID : IJERTCONV6IS16018
 Volume & Issue : RDME – 2018 (Volume 06 – Issue 16)
 Published (First Online): 05012019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Numerical Investigation to Predict the Performance of A GM Pulse Tube Refrigerator using of Pressure Wave Profile
Numerical Investigation to Predict the Performance of A GM Pulse Tube Refrigerator using of Pressure Wave Profile
GM Pulse tube Refrigerator
Pankaj Kumar*, Ajay Kumar Gupta, M. Kumar Ph.D. Scholar, Cryogenics Engineering Laboratory, Department of Mechanical Engineering
National Institute of Technology, Rourkela Odisha, India, 769008
R. K. Sahoo
Professor, Cryogenics Engineering Laboratory, Department of Mechanical Engineering National Institute of Technology, Rourkela
Odisha, India
Abstract The development of GM pulse tube over Inertance pulse tube is to come across the requirement of high heat capacity. The performance of the GM pulse tube refrigerator was greatly improved by introducing different types of pressure wave profiles. Here, 3D model of DIPTR is created using ANSYS Design Modeler and Computational Fluid Dynamics (CFD) solution approach is chosen for numerical simulation purpose. The detailed study on cool down behavior at CHX of DIPTR, pressure variation at inlet and also heat transfer has been performed using CFD package FluentÂ®17 software. Four numbers of cases have been chosen in which the pressure User Defined Functions (UDFs) are different and all other parameters remain unchanged. Pressure UDFs are applied at inlet as a boundary condition to define the oscillating motion of piston inside the pistoncylinder arrangement. The four pressure UDFs are of different wave forms such as Sinusoidal, Rectangular and Trapezoidal. The operating frequency for all cases is 2 Hz. After simulation using four pressure UDFs, it has been found that pressure UDF generating trapezoidal pressure wave is more efficient than other pressure UDFs and a temperature of 110 K is obtained using Trapezoidal pressure wave profile.
Keywords Computational Fluid Dynamics (CFD); Cryogenic Temperature; Double Inlet Pulse Tube Refrigerator (DIPTR); Pressure wave profile; Numerical simulation.

INTRODUCTION
Sl.
No.
Component
Radius (m)
Length (m)
Thermal Condition
1.
Transfer Line
5.30E03
11.5002
Isothermal
2.
After Cooler (AFT)
22.00E03
20.00E03
Isothermal
3.
Regenerator
20.00E03
21.00E02
Adiabatic
4.
Cold Heat Exchanger (CHX)
15.00E03
18.0E02
Adiabatic
5.
Pulse tube
15.00E03
25.00E02
Adiabatic
6.
Hot Heat Exchanger (HHX)
22.00E03
20.00E03
Isothermal
7.
Surge Volume
10.0E02
30.00E02
Isothermal
The initial observation on one end of a hole tube called pulse tube with a pressure wave at low frequency reported in early 1960s [1], and known as pulse tube refrigeration. Numerous developments in such field have been reported, however, Gifford McMahan (GM) type pulse tube refrigerator is one of them which has been used in case of high cooling capacity requirement [2, 3, 4].
The argument, need of analyzing a pulse tube is still a question or requirement, is still debatable. Because one hand is to understand the thermoacoustic phenomena and in other hand is establishing the technique to design a pulse tube for a desired requirement. Such intentions provoke us to concentrate the research on establishing a methodology to analyze a GM pulse tube using commercially available CFD
software [5]. The fact, Why, is acknowledged by the audience, however, how with trustworthiness is still a prerequisite.
From aforementioned literature review it can be concluded that 3D analysis of pulse tube is a basic requirement which has not been addressed in detail till now.
In present extensive numerical investigation, the performance of different wave forms such as Sinusoidal, Rectangular and Trapezoidal, has been optimized to attain the lowest cooling temperature at CHX. To start with, the benchmark GMPTR has been examined by applying Sinusoidal wave and a robust CFD technique has been established. Various issues in numerical investigations have been addressed in detail. Moreover, the recommendations have been made to achieve steady state solution with low computational overhead.

SIMULATION MODEL AND BOUNDARY
CONDITIONS
In order to justify the stability and reliability of the proposed numerical technique, a GM pulse tube from literature has been taken which has been experimentally investigated [1].
TABLE 1. DIMENSION AND THERMAL BOUNDARY CONDITIONS OF GM PULSE TUBE (KASTHURIRENGAN S.ET AL.,2001 [2])
The dimension of GM pulse tube under investigation has been mentioned in Table.1. The detailed investigation has been carried out to demonstrate the complexities such as the requirement of desired meshing, appropriate contacts, boundary conditions, solver configurations, and necessity of User Defined Functions (UDF) involved in simulation to achieve desired result. The Fig.1 represents the geometry of the GM pulse tube under investigation. To save the computational time, the symmetric sectional geometry has been simulated.
p v2
E h
2
One thing is very important i.e. all the properties is applicable only for helium as a working fluid. All the above equation is not applying for after cooler, cold heat exchanger (CHX), hot heat exchanger (HHX) and regenerator components. These components are modelled as porous media.
The resistance of gas flow through the porous media is modeled by an additional momentum source term to the standard fluid flow equations. The source term is nothing but it is an additional pressure drop term due to porous matrix. This source term consists of two parts: a viscous loss terms and an inertial loss term. For a simple homogeneous porous media, the source term is [5]
1
Si v C2 2 f
v v
Where, and C2 are the permeability and the inertial resistance factor, respectively. Thus, the mass, momentum and energy equations for the porous media can be expressed as [7]
1 r v v 0
FIG 1. SCHEMATIC GEOMETRY OF THE GM PULSE TUBE UNDER
t f r
f r r f x
INVESTIGATION; (A) 3D GEOMETRY, AND (B) SECTIONAL f v r f vv p . Si
GEOMETRY.
t
f f s s f f
E 1 E . v E p
t
. k 1 k T .v
f s
where is porosity of the porous media; k f and ks
are
thermal conductivity of fluid and solid structure, respectively. The axial thermal conductivity of the solid structure is set to be about 10% of the screen material due to the contact thermal resistance between the screen layers [8].
The following onedimensional empirical equations are considered for the pressure gradient in the porous media [7]:
p
fF
f
2 v2
FIG 2. FINITE ELEMENT DISCRETIZATION OF 3D AXISYMMETRIC MODEL OF AN GMPTR.
x dh
where fF is fanning friction factor based onsteady flow correlation and expressed as:
nd 2 33.6
nd
vd
h h
f h

GOVERNING EQUATIONS
fF
4l
Reh
0.337 , R
4 2 eh
Continuumbased conservation equations can be applied everywhere in the system. The mass, momentum and energy equations solved by fluent are as follows [6]:
Substituting Eq. (10) in Eq. (9) the pressure gradient in porous media becomes
p 33.6nv n 2 v2
0.337 f
f
f r
1 r V
f x
V 0
—— (1)
x 2l
2 2
t r r x
v vv p .
t
f f
where l is the mesh distance; is the opening area ratio of screen; dh is the hydraulic diameter of screen, n is the number of packed screens per length; is the viscosity. Comparing Eq. (11) with Eq. (5) leads to
t f E .v f E p .kf T .v
D 1 33.6n ;C
2
0.337n
Where,
2l 2 2
For a perfect stacking of square mesh screens in which the
weaving causes no inclination of the wires and screen layers
are not separated, and can be calculated as [9]
with porous media using equivalent viscous
1
, x 0.0254 , x 12 / x 2
resistance (1.085e+09 m2), inertial resistance (5750
4x t d m t t t w
where m is mesh per inch; dw is wire diameter of screen in meter. The geometry of porous medium for regenerator is 304SS woven wire mesh of 200 mesh size and for other heat exchangers it is copper wire mesh of 100 mesh size is used for the simulation.

ASSUMPTIONS AND BOUNDARY CONDITIONS
From literature, it can be noticed that the CFD analysis of GM pulse tube still persists as a challenge for researchers due it its computational overhead. The present extensive numerical investigation tries to address such issue in detail. The computational fluid dynamics of the GM pulse tube refrigerator has been investigated in ANSYS/Fluent platform for a given dimension [10]. The finite element discretization has been carried out in ANSYS meshing module. The mid side node has been dropped to reduce the number of nodes. In general, the boundary layer meshing used to be carried out to solve the fluid flow problem. However, in present situation the compressor motion generates a pressure wave of very high amplitude (approximately Â± 5 bar on the operating pressure 16 bar). So, the effect of boundary layer may be neglected to reduce the computational overhead. At the start of the simulation, the temperature in the system has been assumed to be at 293 K. Table.1 shows the required thermal boundary conditions on all components, corresponding to executed experiment [11]. The simulation has been carried out using below said assumptions.

The computational overhead is due to number of elements and nodes. Moreover, the piston simulation
using dynamic meshing adds significant more computational overhead.

In order to reduce the computational time, the piston is simulated using pressure user defined function
(UDF). The symmetric sectional analysis has been carried out to reduce the computational overload for parametric optimization.

The computational overhead is high due to time duration required in transient analysis to attain steady
state, demonstrated in first section. In order to reduce the computational time, the transient analysis has been proposed to initiate from a steady state solution using a known temperature profile, discussed in subsequent section.

The simulation has been carried out using ideal gas
model (Helium) by solving NaiverStrokes continuity, momentum and energy equations including enabling k turbulent model to simulate flow in porous media.

The operating pressure across GM pulse has been taken to be 16 bar and the operating frequency of the compressor piston motion has been taken to be 2 Hz.

In realtime investigation [11], the AFT, regenerator,
CHX, and HHX have been prepared using fine wire meshes (Copper wire of 100 mesh size). In simulation, those components have been modelled
m1) and porosity (0.697), following empirical
models [12].

The conduction equations corresponding to porous media model have been solved by using steel as material apart from CHX, where copper has been taken.


GRID INDEPENDENCE TEST
The grid testing has been carried out in HP Compaq Elite LE1902x computer having 16GB of RAM. The testing has been carried out for different numbers of element and nodes for 2Daxisymmetric model and 3D model separately. The observations have been tabulated below. In order to ensure the accuracy of the computational grid, a details verification of grid independence test for the cold end temperature in the IPTR with five different grid number are checked.
TABLE 2: GRID TESTING OF IPTR (3D ANALYSIS)
Parameter
Casea
Caseb
Casec
Cased
No. of Elements
39878
40148
48306
55248
No. of Nodes
34958
35423
37687
44978
Lowest Temperature
140
125
110
108.5
Sl.
No. (a)
(b)
(c)
(d) Computational
Time 150 175 200 232
From above analyses, it can be inferred that adequate meshing plays a vital role in achieving the desired result. From observations, it can be recommended that the CHX, and HHX should be meshed using triangular or tetrahedron elements for the requirement of adequate heat transfer. The pulse tube and regenerator including transfer lines should be modeled using rectangular or hexahedron elements using sweep method in case of 3D simulation, correspondingly.

RESULTS AND ANALYSIS
The present paper aims for developing a procedure to generate such predictions of the PTR performance for different pressure waveforms. This will help to determine an optimum pressure waveform for a given PTR configuration and will be useful in a significant way for the design of the rotary valve. The performance of the GM PTR system, various wave forms are studied. The performance of the GM PTR system, various wave forms are studied.
Case1. Pressure UDF (Sinusoidal wave form)
Pressure UDF having sinusoidal wave form is as below:
#include "udf.h" DEFINE_PROFILE (unsteady_pressure, thread, position)
{
face_t f;
real t = CURRENT_TIME; begin_f_loop (f, thread)
{
F_PROFILE (f, thread, position) = 101325.0 Ã— (10 + 5:0 Ã— sin (2 Ã— 3:1416 Ã— 2000 Ã— t));
}
end_f_loop (f, thread)
}
This UDF is compiled in the FLUENT software and applied as boundary condition at inlet and then simulation is carried out. The pressure wave generated at inlet of DIPTR which is shown in fig 3.
FIG 3. PRESSURE WAVE GENERATED AT INLET OF DIPTR FOR CASE – 1.
FIG 4. TEMPERATURE DECREMENT BEHAVIOUR OF DIPTR FOR CASE – 1.
It has been observed that by implementing pressure UDF to simulate piston and transient simulation using steady state
solution by sinusoidal wave profile approach, the GM pulse tube in hand has been investigated. From simulation it has been observed that, the corresponding cooling curve for 510 seconds with a time step of 0.001 second has been shown in Fig.4. From Fig.4, it can be noticed that the lowest temperature CHX achieve is very closely to 134 K.
FIG 5. TEMPERATURE CONTOUR OF DIPTR FOR CASE1.
FIG 6. DENSITY CONTOUR OF DIPTR FOR CASE1.
For detailed understanding the contour of temperature and density distribution have been extracted and have been shown in Fig.6 From contour it has been observed that the temperature distribution gradually decreases in regenerator however in pulse tube it is not that gradual.
Case2. Pressure UDF (Rectangular wave form)
Pressure UDF having rectangular wave form is as below:
#include "udf.h" DEFINE_PROFILE (unsteady_pressure, thread, position)
{
face_t f;
real t = CURRENT_TIME; begin_f_loop (f, thread)
{
for (i = 0; i <100; i++)
{
if (t < (p + 0.0025))
F_PROFILE (f, thread, position) = 20; else if (t = = (p + 0.0025)) F_PROFILE (f, thread, position) = 15; else if (t < (p + 0.005))
F_PROFILE (f, thread, position) = 10; else if (t = = (p+ 0.005))
F_PROFILE (f, thread, position) = 15; else
p = p + 0.005;
}
}
end_f_loop (f, thread)
}
The pressure wave generated at inlet of DIPTR, which is shown in Fig.7. When we implemented rectangular pressure wave profile by similar fashion in DIPTR then by simulation it has been observed that, the corresponding cooling curve for 240 seconds with a time step of 0.001 second has been shown in Fig.8. From Fig.8, it can be noticed that approximately after 210 sec. the lowest temperature at CHX is constant which is very closely to 128 K.
FIG 7. PRESSURE WAVE GENERATED AT INLET OF DIPTR FOR CASE 2
FIG 8. TEMPERATURE DECREMENT BEHAVIOUR OF DIPTR FOR CASE 2
FIG 9. TEMPERATURE CONTOUR OF DIPTR FOR CASE2.
FIG 10. DENSITY CONTOUR OF DIPTR FOR CASE2.
Afterward, various parameters over DIPTR, the contour plot of temperature and density have been extracted and are shown in Fig. 9 and Fig.10.
Case3. Pressure UDF (Trapezoidal wave form)
Pressure UDF having trapezoidal wave form is as below:
#include "udf.h" DEFINE_PROFILE (unsteady_pressure, thread, position)
{
face_t f;
real t = CURRENT_TIME; begin_f_loop (f, thread)
{
for (i = 0; i <100; i++)
{
if (t < (p + 0.001))
F_PROFILE (f, thread, position) = 15 + 5000* (tp); else if (t < (p + 0.0015))
F_PROFILE (f, thread, position) = 20; else if (t < (p + 0.0035))
F_PROFILE (f, thread, position) = 20 – 5000*(t0.0015p); else if (t < (p + 0.004))
F_PROFILE (f, thread, position) = 10; else if (t < (p + 0.005))
F_PROFILE (f, thread, position) = 10 + 5000*(t0.004p); else
p = p + 0.005;
}
}
end_f_loop (f, thread)
}
Last, simulation is carried out by trapezoidal pressure wave profile and it gives a better result as compare with other two pressure profile. Fig. 11 shows the trapezoidal wave profile.
FIG 11. PRESSURE WAVE GENERATED AT INLET OF DIPTR for
CASE – 3.
FIG 12. TEMPERATURE DECREMENT BEHAVIOUR OF DIPTR FOR
CASE 3.
FIG 13. TEMPERATURE CONTOUR OF DIPTR FOR CASE – 3.
FIG 14. DENSITY CONTOUR OF DIPTR FOR CASE 3.
As discussed earlier the GM pulse tube has been investigated in a similar manner. From simulation it has been observed that, the computational overhead time reduce by applying trapezoidal wave profile. The corresponding cooling curve for 200 seconds with a time step of 0.001 second has been shown in Fig.12. From Fig.12, it can be noticed that the lowest temperature CHX achieve is very closely to 110 K.

CONCLUSION
In order to observe the difference of refrigeration performance of a DIPTR at different pressure wave profile, a same pressure input has been applied for the simulation by user define function (UDF). Simulation is carried out for DIPTR using the ANSYS FLUENT 15.0 software. All the four cases having different pressure UDFs have been considered and the cold end temperature is determined for all the three cases. From simulation result, temperature obtained at cold end of DIPTR are 134 K, 128 K and 110 K for Sinusoidal, Rectangular and Trapezoidal wave forms respectively.
From the above simulation result, it can be concluded that the temperature obtained using trapezoidal pressure UDF i.e. 110 K is lower than that of the other cases. So trapezoidal wave form of pressure is much better than all other wave forms.

REFERENCES