1. Introduction
The reduction of fossil fuel consumption is a priority for the industry stimulated by the worldwide policies to reach a climate-neutral scenario in 2050 [
1]. In these circumstances, natural gas (NG) suppliers are reconsidering the procedure to accommodate the gas from the distribution grid to the consumption points’ requirements. This adaptation is carried out in a pressure reduction station where the NG is throttled to liberate the excess pressure. However, based on the Joule–Thomson effect, the NG temperature plunges, leading to the unsuitable hydrates formation [
2,
3]. To avoid this issue, the NG is traditionally preheated before the expansion by means of a boiler, which consumes part of the inlet gas. Thus, efforts have been focused either on the reduction of the NG boiler consumption [
3] or a replacement of this device by, e.g., solar collectors [
4], geothermal exchangers [
5] or turbo-expanders [
6]. In this context, using Ranque–Hilsch vortex tubes (VT) seems a promising alternative since energy input is not required to increase the gas temperature before the expansion. This device is able to generate two separated streams at different temperatures when a compressed gas flow is injected tangentially. The cold stream exits from the central outlet nearby the inlet, and the hot stream is exhausted from the exit at the other side of the tube.
Figure 1 represents the flow distribution inside a dual path or counterflow VT. Although single path VT can also be used, some researchers report that this configuration is less efficient than the counterflow one [
7]. The main advantage of a VT is its simplicity of design and lack of moving parts. Thus, power input is not required and its maintenance is greatly simplified. Based on these advantages, the VT can be considered an excellent choice for heating and cooling applications, e.g., heating, ventilating and air-conditioned (HVAC), cooling of electronics, gas liquefaction or gas heating [
8,
9].
The purpose of the present research is to demonstrate the feasibility of substituting the hot water boiler in an actual pressure reduction station with two vortex tubes. The VTs used are specially designed for working with typical pressure gradients found in NG distribution plants and incorporate a novel system to avoid condensation in the cold inlet. This system diverts part of the flow of the hot exit towards the cold one. Due to this novelty, one of the tubes is first tested in laboratory conditions to assess its performance. This is done using nitrogen as a working fluid. The testing is complemented with computational fluid dynamics (CFD) simulations of the same tube using methane. This allows us to estimate the performance of the tube when installed in the actual station. Finally, the field test under real conditions of demand is performed using two tubes to cover higher demands of NG.
The article is divided into four main sections. The introduction briefly reviews previous experimental and CFD studies on VT. Secondly, the methods and materials section details the laboratory setup, design of experiments, and experimental methods for testing the small VT, as well as the methodology for CFD simulation. The results section is divided into laboratory tests, CFD simulations and field tests, in which the obtained results are discussed and compared between them. Finally, the conclusion section summarizes and highlights the obtained results.
1.1. Review of Experimental Works of Vortex Tubes
The VT has been investigated experimentally to assess its performance under varying conditions of geometrical and thermodynamic nature. Thus, in [
10,
11], the temperature separation as a function of the pressure ratio and flow of natural gas was studied. It was concluded that the VT could be used for liquefying and purifying NG. A more comprehensive research was carried out in [
12] encompassing the effect of geometrical parameters, i.e., diameter and length of the main tube, the diameter of the outlet orifice, the shape of the entrance nozzle; and thermo-physical parameters, inlet gas pressure, type of gas, cold gas/mass ratio and moisture of inlet gas. It was found that the efficiency increases with the ratio length to diameter of the tube, while the diameter of the cold outlet presents an optimum value with respect to efficiency and temperature separation. The temperatures’ separation was found to decrease with increasing cold mass fraction (only fractions higher than
were studied) and increasing number of nozzles. Instead, they are improved by increasing the pressure ratio.
The conditions inside the tube were investigated in [
13]. They measured the pressure, temperature and velocity inside a vortex tube with nitrogen as a working fluid. The inlet condition was changed by rounding the entrance. It was found that this influences the velocities inside the tube and can improve the temperature separation. Similarly, ref. [
14] experimentally measured velocities and pressures inside a VT. They studied the influence of geometrical parameters such as inlet nozzles, cold exit, hot exit and length of the tube. The authors found that a forced vortex is formed at the injection point, which gradually transforms into a free vortex towards the hot exit. The free vortex in the central region transfers energy outwards. The same authors carried out an energy analysis of the VT in [
14]. Their interpretation is that the temperature drop at the cold exit is due to the pressure gradient, while the temperature rise at the hot exit is due to partial stagnation of the flow. The effect of the number of nozzles and their aspect ratio was studied by [
15]. It was also found that increasing number and aspect ratio of nozzle leads to larger mixing zones, which reduces the difference between hot/cold exit temperatures.
An optimization and sensibility analysis were performed in [
16], considering the effect of inlet air pressure, hot tube length, hot tube internal diameter, orifice diameter, and nozzle diameter, on hot-outlet air temperature. The Taguchi method was used to optimize the response. It was observed that all parameters studied have a significant impact on hot temperature. Another study on optimization can be found in [
17]. They experimentally studied the influence of the number of nozzles, cold orifice diameter and inlet pressure on the cooling performance of a VT. The results were analyzed using the response surface methodology. They show that reducing the number of nozzles increases the cold temperature difference for all pressures studied. Instead, the cold exit diameter presents an optimum value which maximizes that difference. The sensitivity analysis proved that the latter parameter was the most influential of the three analyzed.
Regarding the connection of vortex tubes, [
18] studied the performance of two VT in series and in parallel. It was found that the serial setup yields a higher cold temperature difference, while the parallel gives rise to higher hot temperature separation. As with single tubes, it was found that the performance increases with the cold mass fraction and pressure ratio. The present research uses a parallel connection between two tubes in the field tests.
Investigations regarding the use of VT in pressure reduction stations have also been performed recently. Thus, refs. [
2,
19] assess the feasibility of using a VT to reduce the expense of fossil fuels in the preheating of NG, necessary before distributing the gas to consumers. Both studies conclude that this is an interesting solution allowing us to recover energy and simplifying the station maintenance.
The literature on VT is mounting, ranging from design and operation to basic physics. Some reviews can be found in [
7,
20,
21].
1.2. Review of Simulation of Vortex Tubes
Simulations and models aimed at understanding the physics of vortex tubes have been reviewed by [
20]. Ref. [
22] employed CFX software to numerically simulate the compressible flow and energy separation of the VT used in Bruun’s experiment [
23]. They developed an axisymmetric model and solved the mass, momentum, and energy conservation equations, neglecting the gravity contribution and employing the
-
model for turbulence. Ref. [
24] conducted CFD simulations for different types and numbers of nozzles, evaluating the swirl, axial and radial velocity components, as well as flow patterns including secondary circulation flow. Optimum cold end diameter and length/diameter ratio were obtained with CFD and validated with experiments. They used the RNG (Renormalization group)
-
model for turbulence. Ref. [
25] performed a comprehensive 2D axisymmetric CFD model of a real VT and compared it with real experiment measurements. In general, they obtained good agreement between model and experiments. They explored a wide range of cold fractions and temperature separation along it. Ref. [
26] made a comparison of different RANS (Reynolds Averaged Navier–Stokes) turbulence models in predicting the temperature separation in VT. As [
25], they studied temperature separation against cold fraction, and compared against experimental studies. A 2D axisymmetric model was used. They concluded that standard
-
gives the best agreement with experiments in terms of temperature separation. Ref. [
27] present a 3D CFD model to study the effect of cold end diameter in temperature separation of VT. They found an optimal cold end diameter for a given VT and claim that CFD is a powerful tool for obtaining optimum parameters in the design of VT. Ref. [
28] investigated experimentally and through a 3D CFD model the effect of a conical valve length, as well as the inlet pressure and nozzle intakes. They obtained an optimum cone length for maximum efficiency. Ref. [
29] explored the effect of changing the divergence angle of a flared VT and compared it with a straight one, by simulating a 3D CFD model. They applied a periodic boundary condition to reduce the domain to a
sector of the VT.
Regarding the influence of the working fluid on the performance of VT, and specifically with methane or NG there are fewer studies in literature. Ref. [
12] performed experiments with three different gases inside a VT, namely oxygen, helium, and air. They found that, among other parameters, the specific heat ratio of gas influenced the temperature separation. Specifically, they found that the higher the specific heat ratio, the higher the cold temperature separation. Ref. [
30] found a good correlation between experiments and a mathematical similarity relation with the same three working gases and believed that the better energy separation of helium compared with oxygen and air is due to the much lower molecular weight of helium. Ref. [
31] investigated with a 2D compressible CFD model the effect of working gases on energy separation effects. They tested helium, air, nitrogen, oxygen, carbon dioxide, ammonia, and water. They confirmed the results from [
12,
30] by stating that cold temperature separation was favored by increasing the specific heat capacity ratio and by reducing molecular weight. Ref. [
32] performed a 3D CFD simulation, together with an analytical model, comparing the performance of a VT with air and methane as working fluids. They found that both hot and cold temperature separation were lower for methane. In [
33], a 3D simulation based on the geometry of [
25] is validated and then used to explore the energy distribution inside the tube and the role of the Reynolds stresses. Ref. [
34] also performed a 3D CFD simulation of a six nozzle VT, with methane and air as working fluids, arriving at the same conclusion as [
32]: methane simulations result in a slightly lower temperature separation for both hot and cold ends.
An important conclusion from the CFD simulation survey is that a 2D axisymmetric model using a compressible fluid and the - model for turbulence is adequate to reproduce the separation of temperatures at the exit of the VT and the pressure ratio. For its balance of simplicity and accuracy, this will be the simulation setup of choice for this work.
2. Materials and Methods
This section explains in detail the materials and procedures followed to perform the experimental tests and simulation studies. The results are revealed in the next section.
2.1. Laboratory Tests
The experimental tests were done with a VT whose nominal flow is 74 Nm
/h (GR1, 600 mm in length and 31 mm in diameter). A schematic of the experimental setup is shown in
Figure 2. In the field test, an additional VT with double nominal flow (GR2, 710 mm in length and 40 mm in diameter) will be used. This was not tested in the laboratory due to the difficulty to achieve enough gas supply. However, both VTs have the same design and are expected to behave analogously under the right conditions of pressure ratio and flow rate. The objectives of the tests were:
Verify that the flow division corresponds to the design, i.e., cold.
Verify the temperature increase in the hot side and the temperature decrease in the cold side.
Study the dependence between the flow rate and the pressure ratio .
Study the performance of the vortex tube at different ratios, to check the right operation of the system under different conditions.
The tests were carried out under these general conditions:
The working gas was nitrogen.
The working pressure in the inlet manifold ranged from 3 to 9 barg.
The pressure in the outlet manifold (of the testing installation) was set to 1 barg.
The estimated stabilization time of the installation was 2 min. An additional 8 min were added for data collection and verification. In two cases, due to nitrogen supply problems, the tests were shortened by two minutes.
The inlet pressure range is limited from below by the minimum pressure ratio required by the VT to achieve a good separation of hot/cold streams. According to the manufacturer, this ratio needs to be higher than 2. The maximum pressure is limited by the supply system, which consists of a rack of pressurized cylinders. The pressure at the outlet is chosen to be similar to that found in the actual plant in which field tests will be performed. The duration of each test was chosen considering the stabilization time of the VT and the availability of gas at the required pressure.
The experimental tests performed are summarized in
Table 1. The time recorded is the total time, i.e., including the stabilization of the VT.
During a test, the pressure and temperature can be measured by a PLC (Programmable Logic Controller) system, while the flow has to be obtained as the difference between the initial and final count in a flow meter. In
Table 2, we list the instruments used in the laboratory and field tests. The measurements are annotated at a frequency of 1 min. For analysis, the average of the measured values during the test and the final value reached have been used.
2.2. CFD Simulations
The purpose of the simulations is to obtain a prediction of the effect of switching the working gas of the VT from nitrogen to methane. The first step is to develop a model that reproduces the observed behavior with nitrogen, using the experimental tests as a guide. Afterwards, the nitrogen is replaced by methane.
Since the detailed geometry of the VT is not known, there is no advantage in modeling an inaccurate 3D model, given that, as the literature consulted shows, a 2D axisymmetric model can replicate the behavior of a real VT. For that reason, a 2D axisymmetric model was developed. The gas inlet has been modeled as an annular inlet, following the approach of [
25]. The cold and hot outlets are modeled as cylinders, with the hot side partially covered by a cylindrical valve.
Figure 3 shows a recreation of the resulting 3D VT, and a detail of the hot outlet valve.
The flow of gas enters the tube with a swirl momentum due to the angle
. This parameter recreates the orientation of nozzles in real tubes. The swirl is a result of the radial and tangential components of the velocity, while the axial component is set to zero. The cold and hot outlets have been lengthened to allow a developed outlet flow profile, which helps to stabilize the numerical solution. The resulting 2D geometry is depicted in
Figure 4, though only half of it had to be simulated thanks to symmetry respect to the axis. “L” stands for length, “D” for diameter, and subscripts “t”, “h”, “c”, “i”, “hv” stand for “Vortex tube”, “hot outlet”, “cold outlet”, “inlet” and “hot valve”, respectively. The dimensions of GR1 in the simulation are gathered in
Table 3.
Since the detailed geometry is unknown, the procedure to replicate the tube is the following:
- (1)
Develop a 2D model with three free geometric parameters. These parameters are:
- (a)
Size of the hot valve (
in
Figure 4).
- (b)
- (c)
Angle of inclination of the inlet (
in
Figure 5).
- (2)
Perform a parametric sweep of those parameters by CFD to obtain the values that replicate the results provided by the manufacturer, in terms of temperature separation and pressure ratio.
- (3)
Once the right values for the parameters are obtained, perform simulations to replicate the laboratory tests.
Before the parametric sweep, it is necessary to ensure that the used mesh is good enough to trust its results. For that purpose, a mesh sensitivity analysis was performed. Since the above-mentioned geometric parameters are unknown, generic values were used,
mm,
mm,
. As the monitored variables to assess the convergence of the mesh, temperature separation in cold and hot outlets were selected, since they are the most important performance indicators of the VT. The results for the analysis are shown graphically in
Figure 6. In
Table 4, the relative differences between the results of a mesh and its previous one are shown. Given the low variation for values higher than ∼90 k cells, the selected mesh is the one with 95,064 elements. The mesh comprises quadrilateral cells with a few triangular cells. The maximum cell size is
mm while at the inlet it decreased to
mm, and at the exits it is
mm in the hot and
mm in the cold. The average aspect ratio is very close to 1 indicating that the mesh has a high quality.
Both nitrogen and methane are considered as compressible ideal gases, with temperature-dependent thermophysical properties. Navier–Stokes equations for a 2D axisymmetric problem with swirl were solved, along with energy equation and the two equations of the turbulence model. Turbulence has been modeled through the
-
standard model, as suggested by many other researchers, such as [
25,
26,
27,
28]. An
insensitive near-wall modeling was employed, the so-called Menter–Lechner formulation. The viscous heating option was enabled since, for compressible flows, this factor can be important. According to [
35], this effect is important when the Brinkman number
approaches or exceeds the unit. In the performed simulations, a value of
∼3 was found, which justifies the use of the viscous heating effect.
For solving the set of equations, a coupled pressure–velocity coupling scheme was used, with a pseudo-time method, using the conservative length scale method, and a second order upwind scheme was used for all the equation discretizations. Furthermore, the equations were solved using a double-precision solver. Convergence was assessed by monitoring the residuals, the average temperature at both outlets, as well as the pressure at the inlet and outlets. When the temperature and pressure showed steady conditions, the simulation was considered to be converged.
Regarding the boundary conditions, walls are considered adiabatic and non-slipping, since they are well insulated from the environment; the mass flow inlet is fixed, as well as the outlets, for ensuring the cold ratio of
provided by the manufacturer as an operating condition. These boundary conditions, along with their numerical values, are collected in
Table 5. For the inlet turbulent boundary conditions, the hydraulic diameter and a turbulence intensity of
were set.
With the above boundary conditions, an optimization study was performed to find the three values of the geometric parameters, which result in the pressure ratio and temperature separation found in the laboratory for nitrogen gas, reported in the results section.
After replicating the experiments with nitrogen, the working fluid is changed to methane in the simulations to obtain predictions of temperature separation when VT is used in the pilot natural gas plant. The same boundary conditions were used (
Table 5; note that, while normalized volumetric flow is the same, mass flow changes to adapt to the change in density between nitrogen and methane).
4. Field Tests
The last step of the research was to implement the studied vortex tubes in a real environment, under typical natural gas regasification station conditions. The station where the two VT (see
Section 2.1) have been tested (see schematic in
Figure 11 top) has a single liquefied natural gas (LNG) tank and a series of atmospheric heat exchangers, where the regasification occurs. Finally, the gas is heated in a boiler, odorized and injected into the grid. The purpose of the vortex tubes is to substitute the gas heating system, thus avoiding water-heating costs and maintenance, and reducing the pressure of the re-gasified natural gas to that of distribution and supply. Hence, the vortex tubes are placed (
Figure 11 bottom) after the first set of atmospheric heat exchangers—this way they can receive the pre-heated gas (at near atmospheric temperatures) and separate it in the hot and cold flows. The hot flow goes directly to the outlet manifold, while the cold one passes through another atmospheric heat exchanger, with the purpose of raising its temperature and thus obtaining a greater final temperature. Then, the cold flow becomes mixed with the hot one in the outlet manifold. Finally, the gas is odorized and injected into the grid.
The gas inlet pressure here was 4 barg, which is the outlet pressure of the LNG tank, and the outlet pressure is around 1 barg. Actually, the outlet pressure varies slightly depending on the demand of gas from the grid. The gas pipes and the atmospheric heat exchangers are both exposed to the ambient temperature, which during the tests ranged from 0 C to 15 C. The main objective of these tests was to assess whether the vortex tubes can replace the current heating system and compare their performance with the observations made in the laboratory and the simulations.
The results of the above operations under various VT flow and pressure parameters during over a month (34 days) period are presented in
Table 10. Note that the first 3 rows of the
Table 10 correspond to performance of the whole facility, while the other 10 rows are divided into two blocks of five for each vortex tube.
The presented data are averaged over time. Due to the dynamic changes in the demand of gas, over a day there are periods in which only one tube works, or both or none. In practice, however, the case in which both VT are required at the same time never arose. For this reason, the flows shown do not add up to give the total of the facility. Actually, one can deduce that most of the time it is the bigger tube (GR2) that is working. Due to the average, the peak values reached in the temperatures are somewhat higher than those reported in the
Table 10, but those peaks are short-lived and, as such, do not represent any long-term behavior.
Period 1 and Period 2 were dedicated to the fine tuning of the facility and the values reported are not representative. In Period 3, an adequate adjustment of the elements of the facility was reached. Effectively, the pressure ratio is above 2 for both tubes, both have a hot fraction of and the flow is very close to the nominal. During this period, it is achieved the maximum average gain in temperature, C. During Period 4, unfortunately, an error in the manipulation of the hot exit valve of GR1 caused an increase in this fraction, up to . This immediately causes a drop in the performance of the tube, especially noticeable in the hot exit temperature. This, in turn, produces a decrease in the overall gain of temperature of 1 C. Instead, GR2 continues to perform well in periods 4 to 7, even though the pressure ratio is barely 2, i.e., in the limit of validity by design.
Regarding the maximum temperature differences (see
Figure 12), it has been observed that the largest ones are registered when the inlet flow rate is close to the nominal flow. This is especially noteworthy in vortex GR1, where the maximum temperature differences are the highest at 73 Nm
/h, reaching the largest hot temperature difference at
C, and the largest low temperature difference at
C, with pressure ratios always above 2. It is expected, according to laboratory results, that when operating with higher pressure ratios larger temperature differences should be reached. In the case of vortex GR2, the largest temperature differences are
C for the hot side and
C for the cold side, for an inlet flow of
Nm
/h and a pressure ratio above 2. A larger temperature difference should also be expected in this case for higher pressure ratios, as for the GR1 case. Additionally, it is also worth noting that the vortex should always be operated with pressure ratios (
) above 2, as this also determines the amount of energy that can be transferred to the gas in the form of heat. The higher this pressure ratio is, the larger temperature differences should be observed. Summarizing, the largest temperature differences should be expected when the inlet flow rate is close to the nominal flow and the pressure ratio is as high as possible.
These results are different in comparison with those obtained in the laboratory, where the mean temperature differences are around 25
C (hot side) and 5
C (cold side). The reduction in temperature gain was predicted by our CFD simulations (see
Table 9). Some authors [
31] have reported the following observations regarding the temperature separation:
When the working gas is switched from nitrogen to natural gas (), the molecular weight is decreased (from to ) and the ratio is decreased (from to at 300 K and 1 atmosphere). Clearly, the matter of the dependence of the exit temperatures with respect to the properties of the working gas deserves more investigation.
It can also be noted that vortex GR2 raises the total temperature differences slightly less than GR1. This is due to the average flow being below nominal and the pressure ratio just above 2. These two parameters are key to the optimum performance of the tubes. Should the conditions in the plant allow for a higher inlet pressure, both parameters would increase, thus allowing it to reach greater gas temperature increments, as was observed in the laboratory tests.
As for the whole facility, the maximum obtained temperature differences are shown in
Figure 13, again showing a clear dependence on the pressure ratio and flow rate. In this case, the final temperatures are lower than those observed in the laboratory due to the reasons explained above. These temperatures are around 7
C at their best, with mean values of around
C. Thanks to these data, the dependence on the pressure ratio and flow ratio can also be confirmed.
4.1. Thermodynamic Analysis
The enthalpies at the inlets and outlets of each vortex and for the whole vortex facility are reported in
Table 11. The relative differences are computed as
for the VTs and as
for the whole vortex facility. The mismatch in enthalpy between inlets and outlets is very small, almost negligible. For the field tests, the pipes of the facility were insulated, which explains the reduction in losses with respect to the laboratory results. The thermal efficiency (
) of each vortex improves greatly compared to the laboratory observations. This improvement has two sources. For one side, the insulation of the pipes since, as explained before, the probes cannot be placed right at the inlet/outlets of the VTs. On the other side, the definition of
depends on the gas through the ratio of heat capacities
. Since the laboratory and field tests have been made with different gases, this quantity is not directly comparable.
The results for the whole facility are analogous, showing negligible losses in energy. Although the figures are very small, it is worth noting that the losses for the whole facility are systematically larger than those for each VT. This is not only consistent, but strongly indicates that the losses come from the connecting pipes rather than from the tubes.
The relative entropy production is in the range 40–50% for both VTs and for the whole facility, with some values above 50% in the latter. These values are lower than those measured in the laboratory at similar pressure ratio, indicating that the latter were affected by heat losses that have been avoided in the field thanks to the insulation.
Note that the isentropic efficiency is not computed for the whole facility because in this case there is no cold exit, as required by the definition.
4.2. Reduction of Emissions
From the results of period 3 we estimated the non-emitted CO. This period was chosen because the VT worked with the best adjustment of all elements of the facility. The calculation was made as follows:
Considering that the test could not be made during the coldest months of the year, the final reduction in emissions will be higher. This is a significant reduction for a proof of concept as performed in this work.
5. Summary and Conclusions
The purpose of this work has been to assess the possibility of substituting the boiler used in a NG distribution station by a system of two VTs. This substitution allows avoiding gas consumption and simplifies the maintenance of the station. Since the used VTs incorporate a system to prevent frost in the cold exit, the first part of the work consisted on the characterization of one of the tubes under controlled laboratory conditions using nitrogen as the working gas. Secondly, the tubes were tested for several weeks under real conditions in a NG distribution plant in which the pressure ratio changes dynamically according to the demand of the consumers. The field test were made with the same facility and and measurement equipment than the laboratory tests.
In the laboratory tests, it was observed that the temperature of the hot exit rises with increasing pressure ratio (or flow), while the temperature of the cold exit remains essentially unchanged. The latter is due to the mentioned system to prevent frost. The rise of temperature in the hot exit starts at , increases fast up to and somewhat more slowly afterwards. The maximum temperature rise is above 40 C at the highest-pressure ratio that could be tested.
The role of CFD simulations has been to estimate the changes in temperature jumps when the working gas is switched from nitrogen to NG. Since the design of the tube is not accessible, and a faithful simulation is not possible, a workaround has been used. This simple approach allows us to determine that the temperature rises expected with NG are around lower than those observed with nitrogen for a given pressure ratio.
The main conclusions that can be drawn from the laboratory and field tests are:
The overall average gain in temperature is 4 C. The peak gains are about 6 C to 7 C.
A in CO emissions has been achieved.
The optimal hot flow fraction is .
The minimum pressure ratio is , but higher is recommended.
The use of vortex tubes is widespread in HVAC applications but is still limited in other industrial environments. In this work, we report a novel application in the natural gas industry. Moreover, the tubes used are specially designed to work under low pressure ratios available in NG distribution plants. For these reasons, we have developed a complete workflow that runs from the concept of the project on paper to its actual implementation.
Our work confirms the validity of vortex tubes for gas reheating in NG distribution plants. The gain of 4 C may seem modest, but this can be improved increasing the pressure at the exit of the LNG tank. Besides, it is also possible to act on a number of aspects such as facility design, VTs design, their price, and the regulations concerning LNG plants, at least in Spain, where the test was carried out. Combining these factors, it is possible to achieve important gains in temperature with VTs. The gain in temperature combines with the reduction of emissions, which has also been proven, and the savings in energy and maintenance costs, to make the concept presented a very attractive alternative to traditional solutions. We have also confirmed the usefulness of CFD simulation to model the behavior of the tubes. This technique can be used for predicting their performance under changing conditions, optimizing a tube for different applications or designing tubes for specific uses.