Next Article in Journal
Valorization of Corn Cobs for Xylitol and Bioethanol Production through Column Reactor Process
Previous Article in Journal
Torque Ripple Minimization of Variable Reluctance Motor Using Reinforcement Dual NNs Learning Architecture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reducing Energy Consumption and CO2 Emissions in Natural Gas Preheating Stations Using Vortex Tubes

by
Jaime Guerrero
1,*,
Antonio Alcaide-Moreno
1,
Ana González-Espinosa
1,
Roberto Arévalo
1,2,
Lev Tunkel
3,
María Dolores Storch de Gracia
4,5,6 and
Eduardo García-Rosales
4
1
CIRCE-Technology Center, Ranillas Av. 3D, 1st Floor, 50018 Zaragoza, Spain
2
STEAM, Universidad Europea de Valencia, Paseo de la Alameda, 7, 46010 Valencia, Spain
3
Universal Vortex Inc., 410 Princeton Hightstown Rd., Princeton Junction, West Windsor, NJ 08550, USA
4
REDEXIS, Ranillas Av. 1D, 2nd Floor, 50018 Zaragoza, Spain
5
Department of Organizational Engineering, Business Administration and Statistics, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
6
Grupo de Investigación en Organizaciones Sostenibles (GIOS), Universidad Politécnica de Madrid, 28006 Madrid, Spain
*
Author to whom correspondence should be addressed.
Energies 2023, 16(13), 4840; https://doi.org/10.3390/en16134840
Submission received: 12 May 2023 / Revised: 6 June 2023 / Accepted: 16 June 2023 / Published: 21 June 2023
(This article belongs to the Topic Advanced Technologies and Methods in the Energy System)

Abstract

:
This work proposes an innovative method for adjusting the natural gas from the grid to the consumer pipeline requirements in a full-scale pressure reduction station. The use of two counterflow vortex tubes instead of the traditional boiler to preheat the gas before throttling is demonstrated as a powerful alternative. Thus, a reduction of fossil fuel consumption is reached, which amounts to 7.1 % less CO 2 emitted. To ensure the optimal configuration, the vortex tube was thoroughly characterized in laboratory facilities using nitrogen as the working fluid. Various operating conditions were tested to determine the most efficient setup. Computational Fluid Dynamics (CFD) simulations were conducted with nitrogen to validate the behavior of the vortex tube. Subsequently, the working fluid was switched to methane to assess the performance differences between the two gases. Finally, the vortex tubes were deployed at a full-scale installation and tested under real consumption demand. The results obtained from this study offer promising insights into the practical implementation of the proposed method for adjusting the natural gas flow, highlighting its potential for reducing fossil fuel consumption and minimizing CO 2 emissions. Further improvements and refinements can be made based on these findings.

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 0.4 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 60 ° 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 3 /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., 60 % 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 P i n / P o u t .
  • 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 ( L h v in Figure 4).
(b)
Size of the inlet ( L i in Figure 4).
(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, L h v = 75 mm, L i = 2.5 mm, θ = 5 . 5 ° . 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 0.5 mm while at the inlet it decreased to 0.073 mm, and at the exits it is 0.15 mm in the hot and 0.122 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 y + 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 B r = μ U e 2 k Δ T approaches or exceeds the unit. In the performed simulations, a value of B r ∼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 60 % 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 5 % 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).

3. Results and Discussion

3.1. Laboratory Tests

The results obtained are presented: (i) by averaging the temperature jumps during the duration of the experiment (average separation) and (ii) the separation at the end of the test, without averaging (final separation).
The temperature separations are reported in an absolute value. That is, the hot temperature separation represents the increment measured in the hot exit with respect to the inlet. Similarly, the cold separation is the decrement measured in the cold exit respect to the inlet.
We have estimated the uncertainty in our measurements by analyzing the results of repeated experiments under the same conditions. The relative errors are in the range 5–12%. However, uncertainties do not affect equally all experiments. Those carried out with large pressure ratio and large flow rate show less uncertainty than experiments with small pressure ratio and flow rate. As can be seen in the following sections, when the operating conditions are close to their lower end, the behavior of the vortex is not stable and exit temperatures fluctuate. The standard deviations expressed as a percentage with respect to the average lie in a similar range to the relative error and show the same trend.

3.1.1. Temperature Separation Measurements

Figure 7a shows the temperature separation in GR1 as a function of the pressure ratio. On the hot side, there is a monotonic increase in temperature, which reaches around 45 ° C at P i n / P o u t = 4.9 . On the cold side, there is no clear trend. All temperature drops are between 0 ° C and 10 ° C. This is due to the by-pass implemented by the vortex manufacturer to bring part of the flow from the hot outlet to the cold outlet to avoid condensation.
Figure 7b shows the same information as a function of the gas flow rate. The gas mass flow rate increases linearly with the pressure ratio (see Figure 8); therefore, the temperature separation follows the same trend as with the pressure ratio.
The results of this section confirm that the temperature separation increases with the pressure ratio. Compared with traditional designs that do not include a by-pass from the hot exit to the cold one, the hot temperatures reached are around 10 ° C lower. On the other side, the cold temperatures are around zero instead of below, which can be advantageous when mixing the two currents in the exit collector. The cold exit temperature will rise in the field test by means of an atmospheric exchanger. An important observation is that the VT needs to reach its nominal flow rate to start to increase the temperature jump. This can be checked in Figure 7b, where the hot temperature rise decays slightly for flow rates below ≈75 Nm 3 /h.

3.1.2. Thermodynamic Analysis

In this section, the VT is analyzed in terms of its efficiency. To account for the energy losses we compute the enthalpy at the inlet of the tube, H i n , and the enthalpies at the cold and hot exits, H h o t and H c o l d , reported in Table 6. The difference is given as a relative difference (RD) respect to the inlet:
R D = 100 × H i n H c o l d H h o t H i n .
The enthalpies increase linearly with the pressure ratio (or the flow rate); however, there is no systematic trend in the relative difference. These differences are consistently close to 3 % for most experiments, which is a satisfactorily low value. It must be considered that these losses account for the effect of the pipes, given that the measurements cannot be made right at the inlets/outlets of the tube. So, the results can be improved using insulation in the pipes.
The production of entropy ( Δ S ) is given as the relative increase respect to the inlet. Entropy increases notably as it always does in gas expansion, but in this case the effect is enhanced due to the high pressures in the inlet. Correspondingly, the production of entropy increases with the pressure ratio.
To compare these results with those of other groups, we compute the isentropic efficiency, η i s e n , as defined in [36]. This is the ratio of enthalpy change of the inlet and cold outlet to the enthalpy change of the isentropic expansion of the process. Assuming that the gas is ideal, one obtains
η i s e n = T i n l e t T c o l d T i n l e t 1 P o u t / P i n γ 1 γ
where γ = C p / C v is the ratio of heat capacities of the gas. The results are given in Table 6 as a percentage. They are close to 6 % for most cases without any systematic trend. The efficiency reported in [36] changes with the cold mass fraction, showing a maximum between 60–70%. For a cold mass fraction of 60 % , the design value for GR1, their efficiency is ≈12%, indicating that the tube used in this research is less efficient.
From this analysis, it can be concluded that the vortex tube presents fairly low losses of energy, and an efficiency in the range of similar tubes studied in the literature, albeit somewhat lower. The losses in energy could be compensated by using insulation in the connecting pipes.

3.2. CFD Simulations

The results from the parametric study, in which inlet length, inlet angle and hot valve length were swept to reproduce laboratory tests with nitrogen gas (see Table 7), are gathered in Table 8. The internal system for avoiding condensation in the cold stream needs to be accounted for in the CFD simulations. This was done by substituting the complex real system by its thermal effect on cold outlet walls. Thus, a positive heat flux boundary condition has been set, so the stream in the cold side is heated up to the temperature separation obtained in the experiments. Since the real heat flux to the cold stream is not known, it is added as another optimization parameter. The criterion used to decide which is the best set of parameters is to minimize the relative sum of the absolute differences between calculated and target temperatures (cold and hot sides). The values of the parameters were L i = 2.2 mm, θ = 80 ° , L h v = 6 mm and the heat input was 212 W.
The maximum difference between target and calculated values is 11.4 % . The low pressure ratio is due to the point of measurement, which is different in simulations and experiments. In the former, the pressure ratio is computed between inlet and cold outlet; in the real vortex tubes, the pressure is measured before gas enters the nozzles. It is reasonable to think that a significant pressure gradient takes place between nozzles inlet and tube inlet.
A characteristic temperature profile inside the VT is shown in Figure 9a. Taking advantage of the capabilities of the simulation, we show the result of removing the heating system (Figure 9b), i.e., the result that would be obtained with a conventional VT tube, for comparison. The effect of the heating device can be seen in the cold outlet, more marked in the corners. The cold plume (up to yellow color) penetrates slightly more into the tube than in the non-heated case. But beyond this point, no appreciable differences in the profile are discernible. This result is interesting because it shows that the new design (diverting part of the hot flow to the cold side) does not significantly alter the conditions inside the tube.
It is not possible to validate the inner tube temperature distribution, since there are no measurements of it at any point. Nevertheless, the obtained profiles perfectly match in shape and qualitative features with all the consulted literature.
For the simulations with methane as a working gas, the results, along with the relative differences with respect to the nitrogen simulations, are collected in Table 9.
All the relevant variables take lower values than those from nitrogen. As mentioned in the introduction, reference [32] also obtained a lower value of temperature separation when comparing air (mostly nitrogen) and methane. The hot temperature separation is a 30 % lower than that obtained without heating, but cold temperature separation is very near to zero, significantly differing from the 4 ° C obtained with nitrogen. The pressure ratio to obtain the 60–40 (cold-hot) fraction ratios and with the operating mass flow is about 15–20% lower than those needed by nitrogen.
The temperature profiles inside the tube are depicted in Figure 10, along with a profile of a non-heated case for comparison. The temperature gradients are smoother in the case of methane compared to nitrogen. However, the profiles of both working gases share similar qualitative features. The region with the lowest temperatures is the same in both cases, in a half oval-shaped region around could outlet. And for both gases, the differences between the non-heated and the heated cases are very small and confined to the region near the cold outlet.

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 40 % and the flow is very close to the nominal. During this period, it is achieved the maximum average gain in temperature, 4.16 ° 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 57 % . 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 3 /h, reaching the largest hot temperature difference at 15.35 ° C, and the largest low temperature difference at 14.76 ° 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 13.86   ° C for the hot side and 12.14   ° C for the cold side, for an inlet flow of 145.38 Nm 3 /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 ( P i n / P o u t ) 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:
  • It increases upon decreasing the molecular weight of the gas.
  • It increases upon increasing the C p / C v ratio of the gas.
When the working gas is switched from nitrogen to natural gas ( CH 4 ), the molecular weight is decreased (from 28.02 to 16.04 ) and the C p / C v ratio is decreased (from 1.4 to 1.299 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 5.5   ° 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 R D = 100 × ( H i n H c o l d H h o t ) / H i n for the VTs and as R D = 100 × ( H i n H o u t ) / H i n 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 ( η i s e n ) 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 η i s e n 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 2 . This period was chosen because the VT worked with the best adjustment of all elements of the facility. The calculation was made as follows:
  • The energy gained by the gas between the inlets of the VT and odorization system inlet was calculated by using the following equation for the whole period, with i representing each period sample, m ˙ V T the mass flow rate of gas through VT, and C p , CH 4 the specific heat of natural gas:
    Q ˙ = i = 1 e n d - o f - p e r i o d Δ T i n - o u t m ˙ V T · C p , CH 4
  • The needed CH 4 for heating up the gas is obtained assuming that it is burned in a boiler with an efficiency of η = 0.4 , applying the equation below, with LHV CH 4 as the low heating value of natural gas:
    m ˙ C H 4 , c o m b u s t i o n = Q ˙ L H V CH 4 · η
  • The calculated CH 4 mass is transformed into equivalent CO 2 by using the molecular weight ratio of both gases.
  • The CO 2 mass from the previous step is divided by the total CO 2 emitted before the installation of VT (equivalent to the 0.5% of natural gas flow, according to the facilities’ owner, calculated as in the previous steps, but considering the total flow of natural gas through the installation), therefore obtaining the relative value as follows:
    CO 2 r e d u c t i o n [ % ] = m ˙ CO 2 , f r o m - a v o i d e d - c o m b u s t m ˙ CO 2 , f r o m - N G - c o m b u s t · 100 = 7.07 %
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 P i n / P o u t = 2 , increases fast up to P i n / P o u t = 3.5 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 30 % 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 7.1 % in CO 2 emissions has been achieved.
  • The optimal hot flow fraction is 40 % .
  • The minimum pressure ratio is P i n / P o u t = 2 , 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.

Author Contributions

Conceptualization, A.G.-E., R.A. and E.G.-R.; Methodology, J.G., A.G.-E., R.A. and E.G.-R.; Investigation, J.G., A.A.-M., L.T. and E.G.-R.; Writing—original draft, J.G. and A.A.-M.; Writing—review & editing, J.G., A.A.-M., R.A. and E.G.-R.; Supervision, R.A., M.D.S.d.G. and E.G.-R.; Project administration, E.G.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RED CERVERA H24NEWAGE (CER-20211002) “Desarrollo de tecnologías avanzadas de producción, almacenamiento y distribución de hidrógeno, y su transferencia industrial para la nueva era del hidrógeno en España” that has received funding from the “Centro para el Desarrollo Tecnológico e Industrial” (CDTI) from the Spanish Ministry of Science and Innovation. The APC was funded by the same source.

Data Availability Statement

Data is contained within the article.

Acknowledgments

REDEXIS conceived this project and coordinated its execution. Vortex tubes were designed by UVI. The blueprints of the testing facility were developed by Alberto Frisa (CIRCE) and Ana Isabel Lozano. The facility was built by Ogisa Infraestructuras and the remote control by SICA Ing. During the experiments we had the assistance of Óscar Puyo (CIRCE).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Commission, 2050 Long-Term Strategy. 2018. Available online: https://ec.europa.eu/clima/eu-action/climate-strategies-targets/2050-long-term-strategy_en#commissions-vision (accessed on 9 August 2022).
  2. Danieli, P.; Masi, M.; Lazzaretto, A.; Carraro, G.; Volpato, G. A Smart Energy Recovery System to Avoid Preheating in Gas Grid Pressure Reduction Stations. Energies 2021, 15, 371. [Google Scholar] [CrossRef]
  3. Olfati, M.; Bahiraei, M.; Veysi, F. A novel modification on preheating process of natural gas in pressure reduction stations to improve energy consumption, exergy destruction and CO2 emission: Preheating based on real demand. Energy 2019, 173, 598–609. [Google Scholar] [CrossRef]
  4. Cascio, E.L.; Ma, Z.; Schenone, C. Performance assessment of a novel natural gas pressure reduction station equipped with parabolic trough solar collectors. Renew. Energy 2018, 128, 177–187. [Google Scholar] [CrossRef]
  5. Farzaneh-Gord, M.; Ghezelbash, R.; Arabkoohsar, A.; Pilevari, L.; Machado, L.; Koury, R. Employing geothermal heat exchanger in natural gas pressure drop station in order to decrease fuel consumption. Energy 2015, 83, 164–176. [Google Scholar] [CrossRef]
  6. Arabkoohsar, A.; Farzaneh-Gord, M.; Deymi-Dashtebayaz, M.; Machado, L.; Koury, R. A new design for natural gas pressure reduction points by employing a turbo expander and a solar heating set. Renew. Energy 2015, 81, 239–250. [Google Scholar] [CrossRef]
  7. Yilmaz, M.; Kaya, M.; Karagoz, S.; Erdogan, S. A review on design criteria for vortex tubes. Heat Mass Transf. 2009, 45, 613–632. [Google Scholar] [CrossRef]
  8. Linderstrøm-Lang, C.U. Gas separation in the Ranque-Hilsch vortex tube. J. Heat Mass Transf. 1964, 7, 1195–1206. [Google Scholar] [CrossRef]
  9. Jejurkar, A.N.; Shukla, A.N. An Overview on Vortex Tube Applications. In Proceedings of the 1st National Conference On Recent Research In Engineering & Technology, Mogar, India; 2015. [Google Scholar] [CrossRef]
  10. Poshernev, N.V.; Khodorkov, I.L. Natural-gas tests on a conical vortex tube (CVT) with external cooling. Chem. Pet. Eng. 2004, 40, 212–217. [Google Scholar] [CrossRef]
  11. Poshernev, N.V.; Khodorkov, I.L. Experience from the operation of a conical vortex tube with natural gas. Chem. Pet. Eng. 2003, 39, 602–607. [Google Scholar] [CrossRef]
  12. Saidi, M.; Valipour, M. Experimental modeling of vortex tube refrigerator. Appl. Therm. Eng. 2003, 23, 1971–1980. [Google Scholar] [CrossRef]
  13. Gao, C.M.; Bosschaart, K.J.; Zeegers, J.C.H.; de Waele, A.T.A.M. Experimental study on a simple Ranque–Hilsch vortex tube. Cryogenics 2005, 45, 173–183. [Google Scholar] [CrossRef]
  14. Xue, Y.; Arjomandi, M.; Kelso, R. Experimental study of the flow structure in a counter flow Ranque-Hilsch vortex tube. Int. J. Heat Mass Transf. 2014, 55, 5853–5860. [Google Scholar] [CrossRef]
  15. Avci, M. The effects of nozzle aspect ratio and nozzle number on the performance of the Ranque-Hilsch vortex tube. Appl. Therm. Eng. 2017, 50, 302–308. [Google Scholar] [CrossRef]
  16. Kumar, G.S.; Padmanabhanb, G.; Sarmac, B.D. Optimizing the Temperature of Hot outlet Air of Vortex Tube using Taguchi Method. Procedia Eng. 2014, 97, 828–836. [Google Scholar] [CrossRef] [Green Version]
  17. Bovand, M.; Valipour, M.S.; Dincer, K.; Eiamsa-ard, S. Application of Response Surface Methodology to optimization of a standard Ranque-Hilsch vortex tube refrigerator. Appl. Therm. Eng. 2014, 67, 545–553. [Google Scholar] [CrossRef]
  18. Attalla, M.; Ahmed, H.; Ahmed, M.S.; El-Wafa, A.A. Experimental Investigation for Thermal Performance of Series and Parallel Ranque-Hilsch Vortex Tube Systems. Appl. Therm. Eng. 2017, 123, 327–339. [Google Scholar] [CrossRef]
  19. Ghezelbash, R.; Farzaneh-Gord, M.; Said, M. Performance Assessment of Vortex Tube And Vertical Ground Heat Exchanger in Reducing Fuel Consumption of Conventional Pressure Drop Stations. Appl. Therm. Eng. 2016, 102, 213–226. [Google Scholar] [CrossRef]
  20. Eiamsa-ard, S.; Promvonge, P. Review of Ranque-Hilsch effects in vortex tubes. Renew. Sustain. Energy Rev. 2008, 12, 1822–1842. [Google Scholar] [CrossRef]
  21. Kirmaci, V.; Kaya, H. Effects of Working Fluid, Nozzle Number, Nozzle Material and Connection Type on Thermal Performance of a Ranque-Hilsch Vortex Tube: A review. Int. J. Refrig. 2014, 91, 254–266. [Google Scholar] [CrossRef]
  22. Fröhlingsdorf, W.; Unger, H. Numerical investigations of the compressible flow and the energy separation in the Ranque-Hilsch vortex tube. Int. J. Heat Mass Transf. 2015, 42, 415–422. [Google Scholar] [CrossRef]
  23. Bruun, H.H. Experimental investigation of the energy separation in vortex tubes. J. Mech. Eng. Sci. 1969, 11, 567–582. [Google Scholar] [CrossRef]
  24. Behera, U.; Paul, P.J.; Kasthurirengan, S.; Karunanithi, R.; Ram, S.; Dinesh, K.; Jacob, S. CFD analysis and experimental investigations towards optimizing the parameters of Ranque-Hilsch vortex tube. Int. J. Heat Mass Transf. 1969, 48, 1961–1973. [Google Scholar] [CrossRef] [Green Version]
  25. Skye, H.M.; Nellis, G.; Klein, S.A. Comparison of CFD analysis to empirical data in a commercial vortex tube. Int. J. Refrig. 2006, 29, 71–80. [Google Scholar] [CrossRef]
  26. Dutta, T.; Sinhamahapatra, K.; Bandyopdhyay, S. Comparison of different turbulence models in predicting the temperature separation in a Ranque-Hilsch vortex tube. Int. J. Refrig. 2010, 33, 783–792. [Google Scholar] [CrossRef]
  27. May, O.T.E.; Mokni, I.; Mhiri, H.; Bournot, P. CFD investigation of a vortex tube: Effect of the cold end orifice in the temperature separation mechanism. Sci. Acad. Trans. Renew. Energy Syst. Eng. Technol. SATRESET 2011, 1, 84–89. [Google Scholar]
  28. Rafiee, S.E.; Sadeghiazad, M.M. Three-dimensional and experimental investigation on the effect of cone length of throttle valve on thermal performance of a vortex tube using k-e turbulence model. Appl. Therm. Eng. 2014, 66, 65–74. [Google Scholar] [CrossRef]
  29. Dhillon, A.K.; Bandyopadhyay, S.S. CFD analysis of straight and flared vortex tube. IOP Conf. Ser. Mater. Sci. Eng. 2015, 101, 012067. [Google Scholar] [CrossRef] [Green Version]
  30. Stephan, K.; Lin, S.; Durst, M.; Huang, F.; Seher, D. A similarity relation for energy separation in a vortex tube. Int. J. Heat Mass Transf. 1984, 27, 911–920. [Google Scholar] [CrossRef]
  31. Khazaei, H.; Teymourtash, A.; Malek-Jafarian, M. Effects of gas properties and geometrical parameters on performance of a vortex tube. Sci. Iran. 2012, 19, 454–462. [Google Scholar] [CrossRef] [Green Version]
  32. Ahumada, L.M.; Silvera, A.J.B.; Valencia, K.A.M.; Suárez, J.M. Comparison of an analytical and computational fluid-dynamics models of a commercial Ranque-Hilsch vortex tube operating with Air and Methane. CT F-Cienc. Tecnol. Futuro 2019, 9, 61–71. [Google Scholar] [CrossRef]
  33. Vigneshwar, K.; Balakumar, C.; Disimile, P.J. Analysis of energy separation inside a Ranque-Hilsch vortex tube using RANS based CFD. Results Eng. 2021, 11, 100255. [Google Scholar]
  34. Jafargholinejad, S.; Heydari, N. Simulation of vortex tube using natural gas as working fluid with application in city gas stations. In Proceedings of the 12th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics, Costa de Sol, Spain, 11–13 July 2016; pp. 327–331. [Google Scholar]
  35. Ansys. ANSYS Fluent Theory Guide; Ansys Inc.: Canonsburg, PA, USA, 2020. [Google Scholar]
  36. Farzaneh-Gord, M.; Sadi, M. Improving vortex tube performance based on vortex generator design. Energy 2014, 72, 492–500. [Google Scholar] [CrossRef]
Figure 1. Schematic of the flow separation inside a VT.
Figure 1. Schematic of the flow separation inside a VT.
Energies 16 04840 g001
Figure 2. Schematic of the the experimental setup used to assess GR1 along with the instruments used.
Figure 2. Schematic of the the experimental setup used to assess GR1 along with the instruments used.
Energies 16 04840 g002
Figure 3. Three-dimensional recreation of the hot outlet and the regulation valve of the CFD model.
Figure 3. Three-dimensional recreation of the hot outlet and the regulation valve of the CFD model.
Energies 16 04840 g003
Figure 4. Simulated 2D geometry.
Figure 4. Simulated 2D geometry.
Energies 16 04840 g004
Figure 5. Inlet velocity components.
Figure 5. Inlet velocity components.
Energies 16 04840 g005
Figure 6. Mesh sensitivity analysis results.
Figure 6. Mesh sensitivity analysis results.
Energies 16 04840 g006
Figure 7. Separation of temperatures recorded in the GR1 vortex as a function of pressure ratio (a) and the flow rate (b). The bar indicates averaged values.
Figure 7. Separation of temperatures recorded in the GR1 vortex as a function of pressure ratio (a) and the flow rate (b). The bar indicates averaged values.
Energies 16 04840 g007
Figure 8. Flow rate as a function of pressure ratio in GR1 tube.
Figure 8. Flow rate as a function of pressure ratio in GR1 tube.
Energies 16 04840 g008
Figure 9. Temperature profiles of CFD simulations with nitrogen for replication of experiments. (a) GR1 with heating system. (b) GR1 without heating system, for comparison.
Figure 9. Temperature profiles of CFD simulations with nitrogen for replication of experiments. (a) GR1 with heating system. (b) GR1 without heating system, for comparison.
Energies 16 04840 g009
Figure 10. Temperature profiles of CFD simulations with methane for prediction of temperatures. (a) GR1 with heating system; (b) GR1 without heating system, for comparison.
Figure 10. Temperature profiles of CFD simulations with methane for prediction of temperatures. (a) GR1 with heating system; (b) GR1 without heating system, for comparison.
Energies 16 04840 g010
Figure 11. Diagram of a usual regasification facility (top) and the vortex facility (bottom).
Figure 11. Diagram of a usual regasification facility (top) and the vortex facility (bottom).
Energies 16 04840 g011
Figure 12. Temperature difference obtained as a function of the pressure ratio (left) and the working flow rate (right).
Figure 12. Temperature difference obtained as a function of the pressure ratio (left) and the working flow rate (right).
Energies 16 04840 g012
Figure 13. Temperature increase as a function of the pressure ratio (left) and the working flow rate (right).
Figure 13. Temperature increase as a function of the pressure ratio (left) and the working flow rate (right).
Energies 16 04840 g013
Table 1. Experiments carried out in the laboratory.
Table 1. Experiments carried out in the laboratory.
Test P in  (barg) P out  (barg)Duration (min)
13110
23.5110
34110
44.518
5518
66.518
7918
Table 2. Instrumentation used in the laboratory and field tests.
Table 2. Instrumentation used in the laboratory and field tests.
MeasureInstrumentModelUncertainty
Gas flowRotary gas meter2xG65 DN 500.5–1%
Gas flowRotary gas meter2xG65 DN 800.5–1%
TemperatureThermocoupleNIVELCO Thermocount TNClass A (EN 60751)
PressurePressure transducerGems Sensor 8650.25 ±% FS
Table 3. Relevant fixed geometric dimensions of simulated VT.
Table 3. Relevant fixed geometric dimensions of simulated VT.
L t (mm) L c (mm) L h (mm) D t (mm) D c (mm) L t / D t
60060603112.719.35
Table 4. Mesh sensitivity analysis metrics. Relative variation in temperatures’ separation between successive refinements of the mesh.
Table 4. Mesh sensitivity analysis metrics. Relative variation in temperatures’ separation between successive refinements of the mesh.
Elements Δ T hot ( % ) Δ T cold ( % )
31,488--
46,6862.702.43
65,6091.070.97
87,2790.600.55
95,0640.120.11
108,1060.230.21
117,0330.100.09
140,8350.280.27
165,1930.240.23
193,3520.370.36
Table 5. Applied boundary conditions in CFD simulations.
Table 5. Applied boundary conditions in CFD simulations.
Working GasWall TreatmentInlet Temperature ( ° C)Inlet Mass Flow (Nm 3 /h)Cold Outlet Mass Flow (Nm 3 /h)Hot Outlet Mass Flow (Nm 3 /h)
NitrogenAdiabatic 11 73 43.8 29.2
MethaneAdiabatic 11 73 43.8 29.2
Table 6. Summary of thermodynamic and efficiency analysis for GR1. Enthalpies are given in W.
Table 6. Summary of thermodynamic and efficiency analysis for GR1. Enthalpies are given in W.
Test P in / P out H in H cold H hot RD (%) η isen ( % ) Δ S ( % )
12.127343.712815.814806.55−3.86.0375.2
22.298555.303419.485317.47−2.15.3472.8
32.619807.773809.096016.02−0.215.2276.7
42.8110,527.804286.826548.19−2.91.9578.4
53.0211,676.134741.607270.27−2.97.9179.0
63.7714,860.985826.579482.10−3.06.1081.4
74.8819,309.057423.3712,564.07−3.57.4382.3
Table 7. Mean values of temperature separation from experiments with nitrogen gas.
Table 7. Mean values of temperature separation from experiments with nitrogen gas.
Hot Temperature Separation ( ° C)Cold Temperature Separation ( ° C)Pressure Ratio
294.52.1
Table 8. Obtained values from CFD replication of experiments with nitrogen gas, and comparison with experimental values.
Table 8. Obtained values from CFD replication of experiments with nitrogen gas, and comparison with experimental values.
Hot Temp. Separation ( ° C)Cold Temp. Separation ( ° C)Pressure Ratio
ValueDeviation ( % ) ValueDeviation ( % ) ValueDeviation ( % )
25.8−11.04.74.41.86−11.4
Table 9. Obtained results from CFD simulations with methane, and comparison with nitrogen simulations.
Table 9. Obtained results from CFD simulations with methane, and comparison with nitrogen simulations.
Hot Temp. Separation ( ° C)Cold Temp. Separation ( ° C)Pressure Ratio
ValueDiff. ( % ) ValueDiff. ( % ) ValueDiff. ( % )
16.6−34.90.3−93.61.57−15.6
Table 10. Main data obtained during the field tests. Temperatures are measured in ° C, flow rates in Nm 3 /h, and the hot fraction of flow rate is a % of the total.
Table 10. Main data obtained during the field tests. Temperatures are measured in ° C, flow rates in Nm 3 /h, and the hot fraction of flow rate is a % of the total.
QuantityPeriod 1Period 2Period 3Period 4Period 5Period 6Period 7
Global
T o u t T i n 3.13.94.163.223.33.023.75
Flow rate89.9132.3123.80132.4172.3126.2132.2
P i n / P o u t 1.842.042.092.031.9521.86
GR1
Δ T h o t 11.313.813.606.96.26.37.3
Δ T c o l d 8.910.911.0010.7010.210.512
Flow rate64.175.174.8079.1067.770.476.2
Hot fraction39.939.639.2057.805754.558.8
P i n / P o u t 1.832.052.022.021.9522.14
GR2
Δ T h o t 9.111.211.2011.5010.0111.411.8
Δ T c o l d 8.39.69.609.508.41010
Flow rate73.9112.8119.60128.70136.9130.80126.5
Hot fraction53.94442.7039.00354038.8
P i n / P o u t 1.822.032.362.381.8222
Table 11. Summary of thermodynamic and efficiency analysis for field tests. Enthalpies are given in W, the relative difference, isentropic efficiency η i s e n and relative entropy increase Δ S in %.
Table 11. Summary of thermodynamic and efficiency analysis for field tests. Enthalpies are given in W, the relative difference, isentropic efficiency η i s e n and relative entropy increase Δ S in %.
QuantityPeriod 1Period 2Period 3Period 4Period 5Period 6Period 7
GR1
H i n 7577.79755.59510.28428.47061.85715.76897.6
H h o t 3050.93941.95607.44915.94055.93136.04087.0
H c o l d 4522.75811.93957.93513.13004.82572.42806.7
RD0.10.0−0.60.00.00.10.1
η i s e n 19.820.420.920.121.521.322.3
Δ S 40.945.044.144.043.541.245.8
GR2
H i n 1110.63615.43963.47783.310,492.112,143.013,387.2
H h o t 492.81611.41567.13083.23725.75067.65314.9
H c o l d 616.42002.12389.64691.76745.27066.18056.4
RD0.10.10.20.10.20.10.1
η i s e n 19.219.219.118.319.420.120.1
Δ S 40.544.946.843.943.540.944.5
Global
H i n 8688.313,370.913,473.616,211.817,556.814,965.120,284.8
H o u t 8714.113,425.313,530.216,280.117,624.615,050.120,375.8
RD0.300.410.420.420.390.570.45
Δ S 51.755.844.243.455.940.444.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guerrero, J.; Alcaide-Moreno, A.; González-Espinosa, A.; Arévalo, R.; Tunkel, L.; Storch de Gracia, M.D.; García-Rosales, E. Reducing Energy Consumption and CO2 Emissions in Natural Gas Preheating Stations Using Vortex Tubes. Energies 2023, 16, 4840. https://doi.org/10.3390/en16134840

AMA Style

Guerrero J, Alcaide-Moreno A, González-Espinosa A, Arévalo R, Tunkel L, Storch de Gracia MD, García-Rosales E. Reducing Energy Consumption and CO2 Emissions in Natural Gas Preheating Stations Using Vortex Tubes. Energies. 2023; 16(13):4840. https://doi.org/10.3390/en16134840

Chicago/Turabian Style

Guerrero, Jaime, Antonio Alcaide-Moreno, Ana González-Espinosa, Roberto Arévalo, Lev Tunkel, María Dolores Storch de Gracia, and Eduardo García-Rosales. 2023. "Reducing Energy Consumption and CO2 Emissions in Natural Gas Preheating Stations Using Vortex Tubes" Energies 16, no. 13: 4840. https://doi.org/10.3390/en16134840

APA Style

Guerrero, J., Alcaide-Moreno, A., González-Espinosa, A., Arévalo, R., Tunkel, L., Storch de Gracia, M. D., & García-Rosales, E. (2023). Reducing Energy Consumption and CO2 Emissions in Natural Gas Preheating Stations Using Vortex Tubes. Energies, 16(13), 4840. https://doi.org/10.3390/en16134840

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop