1. Introduction
In the modern world, the main problem facing all industries is ensuring energy efficiency and environmental protection. This is why, for various industries, it is necessary to develop installations and systems that are energy efficient and do not harm the environment [
1,
2].
Technological processes implemented at below atmospheric pressures [
3,
4,
5] are becoming widely used. One of the new types of devices used in such systems is a liquid–vapor ejector (LVE), which operates on the principle of jet thermal compression. Water that is not heated to saturation is used as the primary stream. The process of boiling occurs when the active flow flows out of the nozzle. Further pumping of passive flow also occurs, which can be the vapor phase of almost any substance or mixture. Water is a cheap and environmentally friendly product to be used as the working substance of a primary flow, and its use in such systems is innovative and promising.
The efficiency of such an apparatus and a vacuum unit depends on the completeness of the mixing of the primary and secondary flows and achieving the required parameters for the mixed flow at the LVE outlet. One of the main indicators is the degree of vapor overproduction. This parameter determines the amount of steam that must be returned to the cycle. It is evaluated as the difference between the amount of steam removed from the separator as a vapor phase and the amount of the ejected passive flow.
This investigation is relevant to increasing the energy efficiency of existing technological systems and creating new, more efficient and environmentally safe schemes. This is because liquid–steam ejectors are a better alternative to existing steam jet ejectors. They do not require the external generation of working steam, which simplifies the process and does not require additional energy consumption.
Simultaneously, in order to evaluate the expected effect of a new apparatus, it is necessary to perform numerical and experimental studies of liquid–vapor ejectors with different flow path geometries. This would enable the implementation of the most rational and constructive solution.
1.1. Literature Review and Problem Statement
Recent European and world studies have been devoted to exploring the operation of two-phase jet devices on different refrigerants or carbon dioxide [
6,
7,
8,
9]. Although they have significant advantages over single-phase devices, they have been studied due to their complex working process [
10,
11]. The main disadvantage of two-phase jet devices is the cost of the working medium (refrigerant) and issues of safety and environmental friendliness.
Recently, an increasing number of studies have suggested the use of water or water vapor, as in steam jet ejectors, as the working fluid of the primary flow [
12,
13,
14]. However, in steam jet ejectors, steam generation takes place in external steam generators. This requires additional energy consumption and reduces the overall efficiency. In addition, water can have negative properties, including the formation of scale, which is inevitable under such working process conditions. This can be simply solved; the authors of [
15] propose using ethers and aqueous solutions that have neutral properties and do not lead to interaction with a secondary flow at the chemical level.
Study of the geometry and design features of the flow path of two-phase jet devices are presented in [
16,
17,
18]. The efficiency of the working process, both in the apparatus itself and in the whole system, directly depending on its parameters. As noted above, the workflow of two-phase jet devices is quite complex, which does not allow for unifying their development and production. Each device is developed for a specific purpose and a certain degree of pressure increase [
19,
20].
In this study, we investigated the influence of the flow path geometry and degree of steam overproduction in the mixed flow at the outlet of the liquid–vapor ejector on the efficiency of the entire operating process.
1.2. The Aim and Objectives of the Research
The aim of this research is to evaluate the influence of the vapor overproduction value in the mixed flow leaving the liquid–vapor ejector with different flow path geometry on the efficiency of the working process of the vacuum unit.
To achieve this aim, the following objectives were formulated:
Perform a numerical study of a liquid–vapor ejector with different geometries of the flow path with exact parameters for the working fluid of the primary flow at the LVE inlet and calculate the degree of vapor overproduction for each of the options;
Conduct an experimental study of a liquid–vapor ejector with different geometries of the flow path with exact parameters of the working fluid of the primary flow at the LVE inlet and calculate the degree of vapor overproduction for each of the options;
Carry out an exergy analysis of a liquid–vapor ejector with different flow path geometries and determine achievable performance indicators of the LVE and its corresponding vacuum unit.
2. Materials and Methods
In the following section, we describe the selection of operating parameters for condition monitoring and the numerical and experimental methods used to study the liquid–vapor ejector as part of a vacuum unit.
2.1. Initial Data and Design Schemes
The initial data for numerical modeling and experimental research were chosen for the operation of a vacuum unit as part of a plant’s vacuum cooling system for biodiesel fuel production [
21,
22].
Water was used as the working substance of the primary flow. At the inlet to the motive nozzle, it had the following parameters: pressure = 0.4 MPa, temperature = 135 °C, and mass flow = 0.4918 kg/s.
For the working substance of the secondary flow, water vapor was assessed, which had a pressure = 0.02 MPa at the inlet to the receiving chamber. The secondary flow temperature corresponded to the saturation temperature at the corresponding pressure. The mass flow rate was variable and depended on the injection ratio for each LVE with the corresponding geometry.
At the LVE outlet, there was a flow of a two-phase vapor droplet fine structure with a mass flow rate of the mixed flow, which is the sum of the flow rates of the primary and secondary flows. The pressure and temperature of the final flow depended on the geometrical and operating parameters of the flow part of the ejector. According to the study’s goals, LVEs with different flow path geometries were considered, which differed in the geometric shape of the mixing chamber and in the presence or absence of a diffuser (
Figure 1).
2.2. Mathematical Model for Calculating the Working Process of a Liquid–Vapor Ejector
To calculate the averaged parameters along the length of the flow path and total characteristics of the LVE, the system of equations of one-dimensional adiabatic motion in the quasi-equilibrium thermodynamic approximation was used for the selected boundaries of the considered flow section:
The equation of state of a thermally metastable vapor-droplet medium
Mass conservation equation (taking into account the phase transition)
Total enthalpy equation (1st law of thermodynamics)
Momentum conservation equation
Entropy production equation (2nd law of thermodynamics)
The equation of contact stresses arising on the wall
where
υ is the specific volume of a particle of the
i-th distribution group [m
3/kg];
t is the temperature of the particle of the
i-th distribution group (deg.);
w is the average particle flow velocity (m/s);
F is the area of the channel through which the particles of the
i-th distribution group pass (m
2);
z is the dimensionless longitudinal coordinate;
χ—the mass moisture content (kg/kg])
h is the enthalpy of a particle of the
i distribution group (J/kg);
is the shear stress on the channel wall (Pa);
P is the perimeter of the channel through which particles of the
i-th distribution group pass (m),
s is the entropy of a particle of the
i-th distribution group (J/(kg
K); the index “
liq” is the parameter in the liquid state; and index “
st” is the parameter in the vapor state.
The system of Equations (1)–(6) includes thermodynamic parameters (e.g., temperature and pressure) that evaluated parameters directly depend on (e.g., entropy, enthalpy, and specific volume). These values are determined from thermodynamic diagrams. Pressure varies in a range of 0.35–0.45 MPa, and the temperature range is 128–143 °C. Operating parameters (e.g., flow velocity, mass flow of secondary flow, and mass vapor content) vary in a very wide range, from subsonic values (inlet section) of the active flow nozzle to supersonic values (outlet section). At the outlet from the ejector, its value is 110–250 m/s.
The flow in the mixing chamber and ejector’s inlet is subcritical. It significantly reduces energy losses when mixing active and passive currents. It also increases the efficiency of the whole process.
The value of the flow rate of the active stream is determined by calculations from the continuity equation. It ranges from 0.45 to 0.51 kg/s. The value of the mass flow rate of the passive flow is also a calculated value. However, it is a result of direct measurement of the volumetric flow rate using a flow meter.
The values of all measured and calculated parameters depend on the chosen design of the ejector (
Figure 1).
Having solved the system of Equations (1)–(6) with the assumptions adopted in work, it is possible to obtain a formula for calculating the mass vapor content in the mixed flow of the finely dispersed structure of a two-phase vapor-droplet, which at the LVE outlet:
where
is the enthalpy of the mixed flow at the LVE outlet [J/kg] and
is the specific heat of vaporization of the mixed flow at the LVE outlet [J/kg].
In the process of ejection of the secondary flow, which will be saturated water vapor, the two flows will be mixed. At the LVE outlet, there is a flow of a two-phase vapor droplet fine structure. The mixed mass flow rate (
) is the sum of the primary (
) and secondary (
) flows:
During the separation process, most of the flow will be condensed and sent back to the operating cycle of the vacuum unit, and the non-condensed part (
) will be removed as a vapor phase from the cycle. To keep the mass flow rate of the primary flow, and thus all cycle parameters constant, it is necessary to constantly return an additional mass of liquid to the cycle:
Substituting (8) into (9) and performing the transformation, we obtain a formula for determining the degree of vapor overproduction in the LVE working process. Its value will determine how much liquid must be returned to the cycle to maintain material balance:
In a dimensionless form, the formula for determining the degree of vapor overproduction will look like
where u is the injection coefficient determined by the formula
2.3. Experimental Bench for Studying the Liquid–Vapor Ejector as Part of a Vacuum Unit
The experimental bench for condition monitoring and LVE studies in vacuum mode (
Figure 2) consists of an ejector, 1 and two heating tanks, 2 and 18, for heating water and auxiliary communications, pipelines, and fittings.
When testing LVE in vacuum mode, the working fluid of the primary flow is supplied from tank 8 with the help of pump 9 to heating tank 2. Shut-off valves 4 and 12 stop the working fluid supply to the heating tank after filling it. The liquid level in the heating tank is controlled by valve 5. Then, the working fluid is heated by heater 7. The pressure in the tank is controlled by pressure gauge 6. Safety valve 3 is used for emergency pressure relief when the permissible value is exceeded. Discharge into the drain is carried out using valve 15.
After heating and reaching the required pressure that is equal to saturation pressure, the working fluid is fed through a tube-in-pipe heat exchanger 17 into the motive nozzle of LVE. The fluid flow is regulated by valve 16. Cooling water for the heat exchanger is supplied from tank 8. Water is circulated through the heat exchanger using pump 9. Shut-off valves 11 and 12 are used to open or close the cooling water supply to the pumps. The cooling water supply is adjusted using the control valve 13.
In the section from heating tank 2 to the motive nozzle, the following parameters are monitored: pressure—by pressure gauge 25, temperatures—by thermocouples with digital multimeters 26 and 27, and the mass flow of cooling water—by mass flow meter 28.
The secondary flow working fluid is supplied from tank 8 with the help of pump 9 to heating tank 18, which is heated by heating element 19. Stop valve 11 serves to stop the working fluid supply to the heating tank after it is filled. The vacuum value in heating tank 18 is provided by liquid ring vacuum pump 14 and is controlled by vacuum gauge 30. A thermocouple measures the temperature in heating tank 18 with digital multimeter 31. Discharge into the drain is carried out using valve 20.
Then, the vapor is fed into the receiving chamber of the secondary LVE flow. Its flow rate is regulated by valve 21. In the section from heating tank 18 to the fitting for supplying it to the receiving chamber of the secondary LVE flow, the following parameters are controlled: vacuum value—with vacuum gauge 33, temperature—with a thermocouple with digital multimeter 34, volume flow—with flow meter 35.
After leaving the diffuser, the mixed flow enters separator 22, where it is separated into vapor and liquid phases. The vapor phase is discharged into the environment through the pipeline, and the liquid phase is drained through the pipelines into tank 8, from which it is pumped out by pump 9 to fill the heating tanks 2 and 18 or cool the working fluid of the primary flow in heat exchanger 17. Shut-off valve 24 works for the liquid phase from the separator. When studying the parameters of the working vapor jet, digital scales 23 were used to measure the momentum (reaction force).
The experimental LVE model (
Figure 3) was made plane-parallel flat and had a flow path width of 7 mm, made of organic glass (polymethyl methacrylate).
The Mastech MT-838 digital multimeter and TR-01A thermocouple used for temperature measurement had an average systematic error of <1.5% in determining the absolute discharge pressure. The maximum error of temperature measurement by Chromel-Copel thermocouples, complete with a class 0.4 multimeter, did not exceed 1 °C. The maximum error in measuring mass costs, according to calibration tests, was 3%. Class 0.4 standard pressure gauges measured static pressure [
23].
As a result of this study, we obtained the distribution of pressures and velocities of the working fluid of primary and secondary flows in the flow part of the LVE. The measurement error did not exceed 5%, with a confidence level of 0.95, which is acceptable for technical experiments. In another work [
24], a similar measurement accuracy was achieved.
The velocity of working fluid at specific points was indirectly determined by measuring the nozzle’s pressure, jet pulsation, and geometrical parameters.
2.4. Exergy Method for Analyzing a Vacuum Unit Based on a Liquid–Vapor Ejector
The authors used the exergy thermodynamic analysis method to assess the effect of vapor overproduction in a liquid–vapor ejector on its efficiency. It was proposed by J. Tsatsaronis in [
25,
26]. It is a modern tool for evaluating the effectiveness of thermomechanical systems. It simplifies the determination of the total loss of work produced and expended due to the irreversibility of processes in the system under consideration. The main criterion for assessing the perfection of the system is the exergy efficiency (
), determined by the formula
For the “liquid–vapor ejector + separator” system (
Figure 4), the formula for the exergy efficiency can be written according to [
27], as:
where
and
represent the exergy at the LVE inlet of the primary flow and secondary flow (saturated vapor), respectively;
and
represent the exergy at the outlet of the liquid and vapor phase separator, respectively;
and
represent the specific exergy at the LVE inlet of the primary flow and secondary flow (saturated vapor), respectively; and
and
represent the specific exergy at the liquid and vapor phase separator outlet, respectively.
After exergy transformations of Equation (14), taking into account formula (11), there is
Similar to [
27], the exergy efficiency for the vacuum unit system can be written as (
Figure 5):
where
is the power of the circulation pump;
is the exergy of the coolant flow in the heat exchanger; and
and
represent the exergy of the coolant flow at the inlet and outlet of the heat exchanger, respectively.
After exergy transformations, there is
where
is the mass flow rate of the coolant in the heat exchanger and
and
represent the specific exergy of the coolant flow at the inlet and outlet of the heat exchanger, respectively.
2.5. Uncertainty Analysis
Table 1 presents all the values used during the study, the ways in which they were evaluated, and what accuracy they were given.
The LVE efficiency, according to the experiment results, is determined from the data of direct measurements and processing of experimental data.
The steam overproduction degree is determined by formula (11), based on direct and indirect measurements of the values presented in
Table 1. Its total relative error is evaluated by the following formula:
where
is the coefficient that depends on the total number of measurements;
,
, and
are the measured or calculated values of pressure, temperature, and mass vapor content in the mixed flow at the outlet of the ejector, respectively;
and
represent the calculated mass flow rates of active and passive flows, respectively;
,
,
,
, and
represent values of absolute errors of measurement or calculation of the corresponding quantities.
According to expression (18), the total relative error of the steam overproduction degree is 1.65%.
3. Results
To verify the reliability of the results experimentally obtained by condition monitoring [
28], the authors of the work performed additional mathematical modeling in the Ansys CFX Academic software package.
The system of Navier–Stokes equations for compressible media was used. The turbulence was considered using the k-ε model. Thus, a joint solution of the Navier–Stokes equations for time-averaged variables and additional equations to determine the pulsation components were used.
The Rayleigh–Plesset equation was also applied to determine the growth dynamics of steam bubbles. The mass transfer between the liquid and vapor phases was considered using the equations of evaporation and condensation.
The nozzle model created in Ansys CFX was axisymmetric. The calculation grid was automatically generated and consisted of approximately 50,000 cells. The mesh was compacted in places where the geometry of the flow part changed. Water, which boils in the diffuser part of the nozzle, was used as the working fluid of the active flow. Water vapor was used as the working fluid of the passive flow. The properties of the working fluids were taken from the Ansys CFX material database, according to the industry standard, IAPWS.
Input and output parameters are given in
Section 2.1. The boundary conditions were defined as “Inlet” with “Volume Flow Rate”, “Wall” without roughness, and “Outlet” with the “Average Static Pressure”. Ambient pressure was considered equal to atmospheric pressure under normal conditions (101,325 Pa).
As a result of the comparison of numerical and experimental data, the average value of the absolute error did not exceed 1.5%.
Figure 6 shows a graph of the pressure distribution along the length of the flow path of the LVE with different geometry, following
Figure 1.
From
Figure 6, it can be seen that LVE with a cylindrical or conical mixing chamber without a diffuser at the outlet will have lower pressures at the level
= 0.1–0.12 MPa. Accordingly, the degree of increase in the pressure of the secondary flow will be decreased, leading to an additional load on the circulation pump, which must be provided with pressure at the inlet to the LVE
= 0.4 MPa.
For LVE with a cylindrical or conical mixing chamber and diffuser, = 0.22–0.25 MPa.
In addition, according to the results of the experimental studies and calculations using formulas (8) and (12), we can determine the injection coefficient for LVE without a diffuser as u = 0.03–0.05 and with a diffuser as u = 0.07–0.12. Thus, LVE with a diffuser pumps out a larger amount of secondary flow per unit of time, that is, an average of 37.5% greater productivity.
Figure 7 shows a graph of the distribution of mass content vapors along the length of the LVE flow path with different geometry, following
Figure 1. It can be seen that in a LVE with a cylindrical or conical mixing chamber without a diffuser, the output will be much larger (
= 0.271) than in a LVE with a diffuser (
= 0.127). Having made calculations according to Formulas (7) and (11), it is possible to determine the degree of vapor overproduction. For LVE without a diffuser, it is
= 1.517–2.273, and for LVE with a diffuser,
= 1.118–1.426.
After carrying out experimental studies of a liquid–vapor ejector without a diffuser, and performing calculations of the steam overproduction degree by formula (10), the value of was greater. Therefore, it is necessary to return more water to the cycle to ensure that the condition is met.
Figure 8 shows a generalized theoretical dependence of the degree of vapor overproduction and injection coefficient on the mass content of vapor in the secondary flow at the LVE inlet. It can be seen that, with an increase in the mass content of vapor in the secondary flow at the LVE inlet, the injection coefficient increases, but the overproduction of vapor in the working process decreases.
The experimental characteristics of the efficiency of the LVE and its corresponding vacuum unit are shown in
Figure 9.
In this paper, we conducted an analysis using numerical and experimental studies on the effect of vapor overproduction in a liquid–vapor ejector on the efficiency of a vacuum unit. Based on our results, several significant advantages could be identified. First, using LVE, we gained the opportunity to implement a fundamentally new energy conversion cycle in which vapor generation occurred within the working process. It did not require additional equipment that would complicate the technological scheme of the installation. Secondly, using water as the working medium of a primary flow increased the efficiency and safety of the proposed unit.
The following features were found. For secondary flow pressure increases < 4, LVE with a cylindrical mixing chamber and diffuser could be used; for > 4, a conical mixing chamber and diffuser could be used. This was due to the time during which the mixing of the working media of primary and secondary flows and the completion of the processes of heat and mass transfer between the phases at different pressures and temperatures occurred.
This study determined how to choose the geometric parameters of the liquid–vapor ejector flow path, based on the initial parameters of the working media of primary and secondary flows. This made it possible to avoid off-design modes of operation of the ejector at its design stage. The experimental study results indicate correctness, making it possible to obtain such performance indicators.
The closest in design to the proposed type of apparatus are vapor jet ejectors with a water vapor working medium. However, they are multi-stage and have a total efficiency of 2–5%. They have low efficiency because, in one degree of the vapor jet ejector, a pressure difference of only two to three times can be created. In a liquid–vapor ejector, this difference is at a level of 8–10.
Limitations in applying this type of ejector are due to the vacuum values at the inlet to the secondary flow nozzle, which can be achieved at the level of 10–15 kPa. To obtain a lower value of this pressure, it is necessary to use fore vacuum pumps (booster or molecular).
The study carried out is a continuation of the research on the working process of liquid–vapor ejectors in the experimental direction described in [
12,
13,
14,
15,
16]. The experimental data obtained in this work make it possible to determine the operating ranges of the liquid–vapor ejector as part of a vacuum unit, in which its efficiency will be maximized. This allows checking the mathematical model developed by the authors and, if necessary, introducing the appropriate correlation coefficients.
Further work could improve LVE workflow efficiency by profiling the diffuser part of the motive nozzle. However, it is necessary to improve the mathematical model that considers different nozzle geometry and conduct additional experimental studies.
4. Conclusions
As a result of our numerical and experimental studies, we found that in an ejector with a cylindrical or conical mixing chamber with a diffuser, the degree of vapor overproduction was
ψ4 = 1.4, and without a diffuser,
ψ4 = 2.3. The presence of a diffuser contributed to a greater degree of completion of the mixing process and obtained a two-phase vapor-droplet fine structure at the outlet of the ejector. The degree of vapor overproduction reached rational values at the level
ψ4= 1.10–1.25. The degree of steam overproduction was evaluated according to Formula (11), in which the values of mass steam content were taken from
Figure 7.
As a result of the numerical study on the influence of mass steam content in a passive flow (
x02) on the steam overproduction in a mixed flow, we reached the following conclusion: at values of
x02 < 0.215, the steam overproduction degree
ψ4 becomes less than 1. This leads to additional condensation of the passive flow steam (
Figure 8).
As a result of the numerical studies, graphical dependence of vapor overproduction on the secondary flow parameters could be obtained. At the design stage, this can determine the rational variant of the liquid–vapor ejector flow path with maximum efficiency indicators.
As a result of the exergy analysis of the liquid–vapor ejector with different flow path geometry, the efficiency indicators of the LVE with a cylindrical and conical chamber without a diffuser were obtained at the level = 30.5–33%, and with a diffuser at = 20–23.5%.