Effect of Geometrical Parameters on the Coefficient of Performance of the Ranque-Hilsch Vortex Tube ()
1. Introduction
The vortex tube (VT) is a device with a simple geometry without moving mechanical parts [1] [2] [3] [4] . It enables the separation of hot and cold air vortices from the inlet air, as it is supplied with compressed air that flows tangentially through the inlet nozzles [5] [6] [7] . The Vortex tube was first discovered in 1933 by Georges J. Ranque [8] , and in 1947 Rudolf Hilsch improved and modified the design of Ranque [9] . The vortex tube consists of many parts [10] [11] including one or more inlet nozzles, a vortex chamber, a control valve that is located at the hot outlet, a cold orifice and a working tube [12] . When compressed air is injected tangentially through the inlet nozzles into the vortex chamber a strong rotational flow (vortex) is formed at the periphery of the tube wall [13] . The vortex propagates till the end of the tube where some air leaves the tube via the hot outlet. However, adjusting the control valve at the hot outlet forces the rest of the air to reverse its direction and exit from the cold outlet along the centerline of the tube. It is noted that two vortices occur inside the tube, one at the periphery of the tube which exits at the hot outlet at a temperature much higher than the inlet temperature and the other one at the core of the tube with opposite flow direction which exits at the cold orifice with a temperature much lower than the inlet temperature.
The VT is used in many applications such as cooling of airborne electronic components, cooling of gas samples, and cooling of soldered parts including spot welding and ultrasonic welding. It is also used in the separation of air into nitrogen and oxygen rich fluid stream [14] [15] . The VT has many advantages [16] [17] [18] [19] because of its low cost, light weight, reliability, compactness and free maintenance as it has no moving parts. It has an adjustable temperature range and cools without using any refrigerant or electricity.
The results presented in this paper show that the COP of the VT can be enhanced by adjusting the cold orifice to tube diameter ratio (
) and the length to the tube diameter ratio while depending on the adjusted cold mass fraction. This opens a new door to research on enhancing the performance of the VT based on the requirements of the application.
2. Energy Separation inside the RHVT
An energy separation occurs in the vortex tube between the cold and hot streams which causes the temperature changes between the inlet and exit streams. Many researchers suggested various theories to explain the energy separation inside the vortex tube. For example, Ahlborn [20] [21] suggested that the RHVT is a refrigeration device that can be modeled and analyzed as a classic thermodynamic refrigeration cycle with significant temperature splitting, refrigerant and coolant loops, expansion and compression loops and heat exchangers. A secondary circulation is formed in the tube and acts as a refrigerant to transfer the energy from the cold stream to the hot stream. Some researchers [1] [22] [23] [24] reported that the energy separation inside the RHVT is due to work transfer caused by a torque produced by viscous shear between inner and outer stream. Kurosaka [25] attempted to explain the energy separation inside the VT by acoustic streaming. He proposed that the energy separation is due to damping of acoustic streaming.
The effect of friction and turbulence in the VT has always been a major concern of many researchers and was considered by some researchers as the reason for energy separation. For example, using a simple vortex tube model with CFX software and κ-ε turbulence model, Kazantseva et al. [26] attributed the energy separation inside the VT to the energy and gas dynamic interaction of vortices inside the VT. Behara et al. [27] modeled the RHVT using STAR-CD software and concluded that the energy separation inside the VT exists due to the shear work and heat transfer between the hot and cold streams. Bovand et al. [28] studied a 3D model of the VT and concluded that a large portion of the energy transfer inside the RHVT is due to the tangential viscous shear. Aljuwayehl et al. [1] explained the operation inside the VT by the existence of viscous shear acting on a control rotating surface separating the cold and hot flow regions. A torque is produced which causes work transfer between the cold and hot flows. El May et al. [29] reported the formation of a great amount of tangential shear stress of the flow near the wall of the tube, and proposed that the tangential shear work can be considered the most important mechanism of energy separation in the tube. Abdelghany and Kandil [30] presented an operational theory for the vortex tube such the mechanism inside the tube is in analogy to the gas refrigeration cycle (reversed Brayton cycle) where the inner vortex acts as a turbine and that outer vortex acts as a compressor. They [30] concluded that a power transfer exist between bother vortices in which the inner vortex loses energy and the outer vortex gains energy based on the transfer of angular momentum between both vortices.
Some researchers noted the existence of secondary circulations in the vortex tube. They attributed the energy separation to the secondary circulations. For example, Behara et al. [31] confirmed the existence of secondary circulations inside the RHVT at small cold orifice to tube diameter ratio (
) and that such circulations become weak as the
ratio increases and they disappear at
ratios higher than 0.583. They used a CFD model of the RHVT and showed that the total energy transfer was reduced by 21% due to the existence of secondary circulation. Therefore, they concluded that it can be a performance degrading mechanism. Such conclusion was confirmed by Farouk [32] [33] using a large eddy simulation (LES) technique on a CFD model of the VT which showed that the secondary circulation resulted in enhancing the mixing between cold and hot flow which increased the temperature of air exiting at cold outlet. Also, Bovand et al. [28] confirmed that the secondary circulations are performance degrading mechanism. Recently, Kandil and Abdelghany [34] investigated the existence of secondary circulations in the vortex tube at different cold orifice to tube diameter ratio (
) cases. It was confirmed [34] that operating at high cold mass fraction doesn’t have any effect on the formation of secondary circulations. Nevertheless, when operating at low
ratio secondary circulations are formed near the cold outlet but only at low cold mass fraction. They concluded that the secondary circulations have almost no effect in degrading the performance of the vortex tube.
3. Performance Optimization of the RHVT
The performance of the RHVT is affected by geometrical properties such as the length to tube diameter ratio (
) and cold orifice to tube diameter ratio (
). The length to tube diameter ratio
is a very important parameter that attracted the interest of many researches s such as [1] [2] [35] [36] [37] [38] [39] . For example, Bramo and Pourmahmoud [2] [35] [37] conducted a study on the
ratio effect on the performance of the VT at cold mass fraction of 0.3 (at the maximum cooling effect) and concluded that at
ratio of 9.3 the VT has the highest total temperature change. Almost all the other
ratios have the same performance and that the high performance at
ratio of 9.3 is due to the stagnation point along the centerline of the tube. Patel and Bartaria [38] showed that as
increases
increases to reach its maximum value at
ratio of 30 after which the performance of the VT decreases. They reported that the performance of the tube in terms of maximum temperature difference between the hot and cold outlets,
, increases with increasing
ratio up to a critical value after which the performance degrades. Aydin and Baki [40] confirmed that there is an optimum value for the
ratio for the vortex tube. Maurya and Bhavsar [36] concluded that as the
ratio increases
increases to reach an asymptotic state, and to obtain the maximum energy separation effect the
ratio should be greater than 10. Aljuwayhel et al. [1] studied the
ratio effect on the performance at cold mass fraction of 0.3, near the maximum refrigeration point, and reported that as the length increases the performance of the tube increases up to a certain length of tube after which the effect of the
ratio on the performance becomes negligible. Their results showed that there is a critical length of the VT over which the majority of energy transfer takes place and exceeding this length doesn’t affect the energy separation.
Kandil and Abdelghany [34] studied the effect of (
) ratio on the performance of the vortex tube at different cold mass fractions of 0.34, 0.4, 0.428, 0.5 and 0.72.
The effects of varying
on the performance of the VT was investigated by Maurya and Bhavsar [36] by fixing the length and diameter of the tube at 120 mm and 10 mm, respectively, and changing the diameter of the cold orifice. The best performance was achieved at the lowest cold temperature (
), which is the temperature of air exiting at the cold outlet The cold temperature,
, decreased with increasing
till a critical diameter where the
then
increased significantly. They suggested that the reduction of pressure gradient due to increasing the cold orifice diameter (
) is the reason behind the rise of the cold exit temperature after reaching the critical
. Andalibi et al. [41] showed that maximum temperature separation occurs at
ratio of 0.45 and that operating at higher ratios will cause backflow at the cold outlet which will degrade the tube’s performance. Nimbalkar and Muller [3] studied the effect of the
ratio on the energy separation inside the tube and showed that there is an optimum cold orifice diameter and operating at a lower diameter will result in mixing the cold stream with the hot one and operating at a higher diameter will result in mixing the hot stream with the cold one.
Recently, Kandil and Abdelghany [34] studied the effect of the
ratio on the performance of vortex tube in terms of the cold temperature drop and the hot temperature rise. Their results [34] showed that all the
ratio cases exhibit similar behavior where the cold temperature drop increases when decreasing the cold mass fraction until a certain value of the cold mass fraction after which the cold temperature drop decreases significantly. They explained the reason behind the decrease of the cold temperature drop after certain cold mass fraction by the occurrence of back flow at the cold outlet
The results presented in this paper support and build on the research of Kandil and Abdelghany [34] where the changes in the geometrical parameters of the VT can be adjusted to enhance one parameter which is the COP.
4. Numerical Model
An Axisymmetric model of the RHVT was designed and simulated using ANSYS Fluent® software. The model was validated using the experimental measurements of Skye et al. [42] . ANSYS Fluent software solves the following compressible turbulent flow conservation of mass, momentum and energy and state equation [43] (ANSYS Fluent Theory Guide):
Continuity equation in tensor notation:
(1)
Momentum equations in tensor notations (all three components):
(2)
Where
(3)
Energy Equation:
(4)
Equation of State (ideal gas law)
(5)
where
is the Kronecker delta function (
if
and
if
);
denotes viscous stress tensor; H represents the total enthalpy; T is the static temperature, p is the pressure,
is the density and R is the universal gas constant at standard conditions.
4.1. Turbulence Model
The results of the CFD model of the RHVT using the κ-ε Turbulence model showed the best agreement with the experimental results of Skye et al. [42] . Therefore, the κ-ε Turbulence model was selected to simulate the turbulence inside the VT in this study.
The geometry used in the current CFD model is based on the ExairTM 708 slpm vortex tube used by Skye et al. [42] as shown in Figure 1. The current model shown in Figure 2 is an axisymmetric model of the tube with matching geometrical properties to the VT used by Skye et al. [42] , where the working tube length is 106 mm, tube diameter is 11.4 mm and cold exit diameter is 6.2 mm.
The results of the CFD model of the VT became mesh independent when the number of elements reached 17,000 elements. Due to the complexity and sensitivity of the flow inside the VT, the time step of the solution had to be very small
Figure 1. Schematic of vortex tube used in the experiment [42] .
Figure 2. Current CFD Model of the RHVT.
to allow converge such that the computation time of one simulation ranged from 12 to 48 hours using 8 processors on a 3.64 GHz HP workstation Z800 with 16 GB Ram.
4.2. Boundary Conditions
The experimental data from Skye et al. [42] provided some of the boundary conditions for the inlet nozzles and both outlets while the rest of the data were acquired from different previous CFD models [2] [35] which were implemented based on the current configuration. The boundary conditions are summarized in Table 1.
4.3. Numerical Model Validation
The current CFD model was validated by comparing its results with the experimental measurements of Skye et al. [42] using the cold temperature drop (
), and the hot temperature rise (
) defined as the difference in temperatures between the outlet hot air and the inlet air.
Figure 3 shows the comparison of the cold temperature drop of the current model with the experimental measurements and the computational results of Skye et al. [42] .
The results in Figure 3 show that the cold temperature drop increases as the cold mass fraction (
), which is the ratio between the cold air outlet mass flow rate divided by the inlet mass flow rate, decreases up to a maximum
of 40.4 K at a cold mass fraction of 0.34 after which
decreases to reach a temperature of 31.7 K at
of 0.2099. The
reaches a minimum value of 20 K at
of 0.87 which shows that as the cold mass fraction increases to a maximum value the cold temperature drop decreases to its minimum value. It is noted that the current model has better agreement with the experimental measurements than that of Skye et al. [42] model.
Figure 4 shows the comparison between the hot temperature rise results of the current study with the experimental measurements and the computational results of Skye et al. [42] .
The results show that the hot temperature rise increases as the cold mass fraction increases to reach a maximum value of 86.7 K at cold mass fraction of 0.87. Increasing the cold mass fraction further will be impractical due to the very low
Table 1. Boundary conditions used in the current study.
Figure 3. Comparison of the cold temperature drop results with the experimental measurements and computational results of Skye et al. [42] .
Figure 4. Comparison of the hot temperature rise results with the experimental measurements and computational results of Skye et al. [42] .
flow rate of air from the hot outlet. The results show also that at very low cold mass fraction, the hot temperature rise decreases till it reaches a minimum value of 11.4 K at cold mass fraction of 0.2099. The results show also that when operating at cold mass fractions from 0.2 to 0.45 the increase in temperature rise is slow while as the cold mass fraction increases, the increase in temperature rise becomes rapid especially at very high cold mass fractions. The comparison shows that the current model has better agreement with the experimental measurements than that of Skye et al. [42] model.
5. Performance Optimization of the RHVT
5.1. Effects of the Length to Tube Diameter Ratio on the Total Coefficient of the Vortex Tube (COPtotal)
The performance of the VT was determined by Kandil and Abdelghany [34] in terms of total temperature change. When studying the performance of the RHVT, not only the hot temperature rise or the cold temperature drop should considered but also the cold mass fraction, the pressure at the inlet and the pressure at the cold outlet. The VT can be used for cooling as it may be used for heating therefore a coefficient of performance (COP) for the two operations are defined.
The coefficient of performance for the RHVT when it is used as a refrigerator is the cooling power
divided by the input power
. The input power for any refrigeration system is the input to the compressor. But for a RHVT system a compressed air source is used which makes it difficult to define the input power. Fulton [44] proposed that the work used to compress the air through a reversible isothermal compression process from the exhaust pressure up to the input pressure the input work is expressed as follows:
(6)
where
(7)
where
is the specific gas constant,
is the input pressure and
is the cold exit pressure
Therefore the COP of the RHVT as a cooler is expressed as
(8)
The COP for the RHVT when it is used as a heat pump is defined as the heating power
divided by the input power
. The input power used by the system takes the same from as the one used for the cooling RHVT
Therefore, the COP of the RHVT as a heat pump is expressed as
(9)
The total performance of the RHVT can be determined using the summation of the COP as a cooler and as a heat pump which is expressed as:
(10)
Figure 5 shows the relationship between the
and different
ratios at cold mass fractions of 0.34, 0.4 and 0.43. The results in Figure 5 at cold mass fraction of 0.34 show that at
ratio of 2.8 the
is 0.3309 and as the
increases the
increases to reach a maximum value of 0.3746 at
of 7.02 with a difference of 0.0437 from the
at the lowest
ratio of 2.8. After which any further increase in the
ratio will result in a decrease in the
and the performance of the VT to reach a value of 0.3601 at
ratio of 20.2.
The performance of the VT in terms of
at different
ratios was studied at different values of cold mass fraction. The results in Figure 5 and Figure 6 at cold mass fraction of 0.4 show that the
at
ratios above 3.5 decreased significantly changing the trend of the curve. Increasing μ to 0.439 decreased the
even more such that the relation between the
and the
ratio at
in Figure 6 shows a completely different pattern than that at
shown in Figure 5 where the lowest
ratio of 2.8 has the highest
of 0.4936. As the
ratio increases the
decreases to reach an asymptotic value of 0.476 at
ratio of 20.2.
These results show that the
at different
ratios is dependent upon the cold mass fraction, where at low cold mass fraction, as the
increases the
increases to reach a maximum value at a critical
ratio. Any further increase in the
ratio will result in decreasing the
. While as the cold mass fraction increases the difference between the maximum and minimum
is reduced from 0.04 at
to 0.0086 at
. Any further increase in the cold mass fraction will change the trend of the curve
Figure 5. Total COP of the VT at different L/D ratios at cold mass fractions of 0.34, 0.4 and 0.43.
Figure 6. Total COP of the VT at different L/D ratios at cold mass fractions of 0.5, 0.64 and 0.72.
between the
and different
ratios such that the lowest
ratio has the highest
and as the
ratio increases the
decreases to reach asymptotic value at critical
ratio.
In Figure 6, it is shown that the
at a cold mass fraction of 0.64 is higher than that at any other cold mass fraction. Such results show that the best
value will be achieved when operating the tube at a cold mass fraction of 0.64 at any length to tube diameter ratio.
Such results don’t follow the same pattern of the results of total temperature of Kandil and Abdelghany [34] where the best performance of the VT in terms of total temperature is when operating at the maximum possible cold mass fraction. This can be interpreted that the
depend not only on the cold and hot temperature changes as the total temperature change but also on the cold mass fraction and the pressure at the inlet and cold outlet.
5.2. Effects of the Cold Orifice to Tube Diameter Ratio on the Performance of the Vortex Tube in Terms of COP
The cold orifice diameter is an important parameter that affects the performance of the VT. Therefore, several trials were done to optimize the performance of the tube of a diameter of 11.4 mm by varying the cold orifice diameter to pipe diameter ratio,
.
5.2.1. Highest COP Cold (COPc)
When studying the performance of the RHVT taking into consideration not just the highest temperature drop but also the cold mass fraction, the pressure at the inlet and the pressure at the cold outlet, a new dimensionless parameter had to be defined to combine the effect of all the previously mentioned parameters on the performance of the vortex tube. This dimensionless parameter is the
which was defined by Equation (8). Figure 7 shows the relationship between the
and the cold mass fraction.
It is shown in Figure 7 that all
ratio cases follow a similar pattern whereas the cold mass fraction increases the
increases till a critical cold mass fraction which is almost the same in all cases of value of 0.64 after which the
decreases with any further increase in cold mass fraction. The results show that the
ratio of 0.6316 has the highest
of 0.319 at a cold mass fraction of 0.64.
It is shown in Figure 7 that at cold mass fraction of 0.3 the
of the
ratios above 0.5 drops to be lower than the
of the
ratios below 0.5 but the pattern starts when the highest
ratio of 0.6316 drops at cold mass fraction of 0.4. These results show that decreasing the cold mass fraction beyond a threshold, especially for high
ratios, causes a drop in
, resulting in an inversely proportional relationship beyond this threshold.
5.2.2. Highest COP Hot (COPh)
When studying the performance of the RHVT taking into consideration not just the highest temperature rise but also the hot mass fraction, the pressure at the inlet and the pressure at the cold outlet, a new dimensionless parameter had to be defined to combine the effect of all the previously mentioned parameters on the performance of the vortex tube. This dimensionless parameter is the
which was defined by Equation (9). Figure 8 shows the relationship between the
and the cold mass fraction.
It is shown in Figure 8 that all
ratio cases follow a similar pattern whereas cold mass fraction increases the
increases up to a critical cold
Figure 7. Cold COP of the RHVT at different dc/D ratio.
Figure 8. Hot COP of the RHVT at different dc/D ratio.
mass fraction which is almost the same value of 0.64 in all cases after which any further increase in the cold mass fraction will result in decreasing the
. The highest
of different
ratio cases existed at this exact cold mass fraction as stated in the previous Section 5.2.1. The
ratio of 0.6316 has the highest
of 0.2145 at cold mass fraction of 0.64.
5.2.3. Highest COP Total (COPtotal)
From the previous sections, it was noted that the highest
and
were at a cold mass fraction of 0.64 for all
ratio cases The
was defined in Equation (10) as the summation of the
and
, therefore when studying the performance of the RHVT in terms of the
, the maximum
for all the different
ratio cases should exist at cold mass fraction of 0.64 which is confirmed in Figure 9 which shows the relationship between the
and the cold mass fraction. It is shown in Figure 9 that the
ratio of 0.6361 has the highest
of 0.534 at cold mass fraction of 0.64.
It is noted from the results in Section 5.1 and Section 5.2 that the best performance of the vortex tube in terms of
occurs when operating the vortex tube at a cold mass fraction of 0.64 regardless of the length to tube diameter ratio or the cold orifice to tube diameter ratio.
For Future research, these results can be more investigated in terms of flow analysis using 3D modeling of the RHVT to gain more insights about what actually occurs at the cold mass fraction of 0.64. Building on these results, the flow visualization at different
ratios at this cold mass fraction will also be very beneficial.
Figure 9. Total COP of the RHVT at different dc/D ratio.
6. Conclusions
The results of the current model are in good agreement with the available experimental measurements. Many geometrical parameters affect the performance of the tube. In this study, the length to tube diameter ratio effect in terms of
was investigated and the results showed that its effect on the performance of the tube is different at different cold mass fractions where the maximum
at any length to tube diameter ratio occurs at cold mass fraction of 0.64.
Also the effect of the cold orifice to tube diameter ratio on
,
and
was investigated and the results showed that all the studied VTs at different
ratios have the maximum
and the maximum
at a cold mass fraction of 0.64. The highest
ratio of 0.6361 has the highest
and
.
It can be concluded that the best performance of the studied vortex tube at any
or
ratio occurs at a cold mass fraction of 0.64. Future investigation of the physical properties of the tube at such cold mass fraction and its relation to other operating points of the vortex tube should be considered.