Construction of a Numerical Model for Flow Flash Evaporation with Non-Condensable Gas

: Flash evaporation processes are widely used in petroleum, food, chemical, power, and other industries to separate products or extract heat. The liquid is often entrained by non-condensing gas components. This study develops a multiphase, multicomponent, and pressure-driven phase-change-coupled model to numerically study water ﬂash evaporation with non-condensing CO 2 . The model includes the mass, momentum, energy, volume of ﬂuid (VOF), species transport, turbulence (RNG k-ε ), modiﬁed phase-change Lee, and non-condensing CO 2 release governing equations. The steam generation rate and mechanism for pure water and different concentrations of CO 2 are considered. The results show that the numerical model can accurately predict the ﬂash evaporation process and has high accuracy compared with the experimental data. Both the dissolved and entrained CO 2 that are released can severely disturb the ﬂow ﬁeld, leading to an increase in the steam generation rate. Under a 1–10% volume concentration of dissolved CO 2 and 0.0661–0.1688% mass concentration of entrained CO 2 , the maximum increase ratio of steam generation can reach 20%.


Introduction
Low-pressure flash evaporation is a unique evaporation technique that is different from evaporation from an overheated surface.Steam is primarily produced at the liquid surface because ambient atmospheric pressure is lower than the saturation pressure of the temperature of the liquid [1].This technique offers a high evaporation rate without requiring excessively high liquid temperatures.Flash evaporation is widely used in many industrial fields, such as seawater desalination, chemical engineering, food processing, ocean thermal energy conversion, and waste heat recovery [2,3].
Many factors have a significant effect on the flash evaporation rate, such as the liquid temperature, atmospheric pressure, saturation pressure, initial liquid height, and liquid flow rate [4].Several experimental studies have been conducted on this topic.Wu et al. [5] experimentally studied flash evaporation under low-pressure conditions, with the water temperature ranging from 50.0 to 90.5 • C, initial water height of 400 mm to 550 mm, and superheats ranging from 1.5 K to 17.8 K; the results showed that an increase in water temperature and a decrease in the initial liquid level height can enhance the extent of flash evaporation.Zhang et al. [6] studied water flash evaporation with different values of flow rate (400-1200 L h −1 ), water film concentration (0-10%), and pressure (7.4-47.4kPa).Ni et al. [7] experimentally investigated the effect of different sealing air flows on the flash evaporation rate.The valve opening ratio varied from 0.125 to 1. Chi et al. [8] experimentally studied the thermal characteristics of water flash evaporation from a thin liquid film.The experimental operating conditions were an initial temperature in the range of 60-85 • C, a liquid film height of 3-20 mm, and a superheat in the range of 2-22 K. Lior and Nishiyama [9] experimentally investigated the use of electrolytically generated hydrogen bubbles to improve flash evaporation performance.They found that hydrogen bubbles promoted ebullition in the regions of the flash stage, and thus reduced the non-equilibrium temperature difference by up to 15%.Tian et al. [10] investigated water flash evaporation under a microwave.The temperature ranged from 19 to 50 • C, pressure from 4-12 kPa, flow rate from 10-40 L h −1 , and the microwave power ranged from 0.81 to 1.35 kW.Wang et al. [11] conducted experimental and theoretical studies on water flash evaporation and developed a one-dimensional (1D) mathematical model.The model considered four main physical processes: gas outflow due to depressurization, mild evaporation, flash evaporation, and the updated flashing chamber pressure and gas components.Zhang et al. [12] experimentally studied the static flash evaporation of aqueous NaCl solutions at flash speeds in the range of 4.8 × 10 −4 -2.18 s −1 .
While many experimental studies have explored some physical regulations, there still remain gaps in the understanding.However, owing to the limitations of experimental test methods, our understanding of the internal heat and mass transfer mechanisms remains unclear.In addition, experimental tests are time consuming and costly, and it is difficult to reproduce the experimental results in other working conditions.Similarly, there are problems of errors and assumptions in the calculations and numerical model, making the simulation not always completely accurate, but with the development of computers, numerical methods have attracted increasing attention.Xie et al. [13] used the volume of fluid (VOF) and Lee models to simulate the flash point of ethylene glycol.The results showed that the flash point varied at different measurement points, and a lower flash point could be obtained at the measurement point with a higher vapor volume fraction.The vapor volume fraction and temperature showed strong fluctuations, and the measurement deviations of the flash points were closely related to the fluctuations of the vapor temperature.A mathematical model was developed by Cheng et al. [14] to numerically simulate the droplet and film flash evaporation.The results showed that droplet flash evaporation significantly affected the cooling heat transfer of vacuum flash evaporation.Gao et al. [15] proposed a mathematical model to predict the flash evaporation of LiCl droplets released into a vacuum.They found that the pressure is the core factor affecting the evaporation rate, i.e., the lower the pressure, the more intense the flash.Droplets with a higher initial temperature only strengthened the intensity of flash evaporation at the beginning.The model can accurately predict droplet flash evaporation without radiation heat.Ji et al. [16] conducted a three-dimensional (3D) simulation of the flash evaporation of a non-uniform spray in a saturated vapor environment using the Eulerian-Lagrangian framework in the software Fluent.The spray initially underwent violent evaporation and a significant temperature reduction within the vicinity of the nozzle exit.As the spray develops downstream, evaporation becomes moderate.Smaller droplets were distributed in the central spray region, with a comparatively higher evaporation rate and larger velocity variation.Jin and Low [17] used a single-phase model to study the flow-field distribution in a flash evaporation chamber.Wang et al. [18] established a 1D droplet flash evaporation model by introducing a simple heat balance correlation and a phase-transition model to Fick's law.The effects of ambient pressure, temperature, droplet size, and the fraction of vapor concentration on the flash evaporation and spray cooling rate were presented in detail.
The presence of non-condensable gases in the vapor can have a serious effect on the heat transfer process in the condenser, as it reduces the partial pressure of the vapor around the cooling surface and increases the transfer resistance of the vapor to the cooling surface, which ultimately leads to a reduction in the amount of heat transferred and an increase in the vacuum load of the condensing system.Many scholars have studied the effects of non-condensable gas.Dro ży ński [19] studied the significant influence of non-condensable gases on the steam condensation process; in particular, the presence of non-condensable gases increases the pressure in the equipment and affects the condensation efficiency.Balaji et al. [20] investigated the effect of non-condensable gases in a low-temperature thermal desalination system and showed that each unit mass of non-condensable gas flowing into the condensation chamber resulted in 1.6 times the unit vapor loss from the vacuum pump outlet.Genthner et al. [21] studied the effect of non-condensable gases in multiple-effect distillation systems and showed that the presence of non-condensable gases increased the heat transfer area of a horizontal tube condenser by about 10-30%, while each unit of non-condensable gases led to a loss of 20-30 units of vapor under severe operating conditions.
The aforementioned numerical studies all treated liquid as a single phase.However, in the industry, liquids often include and dissolve a significant amount of non-condensing gases.Therefore, we propose a new coupled numerical model that includes multiphase, multicomponent, phase-change, and turbulent-flow modules that can accurately predict the liquid flash evaporation process with a non-condensing gas.The non-condensing gas was considered to be dissolved and entrained into two components in the liquid.The effects of the two types of non-condensing gases on the flash evaporation mechanism and steam generation rates are discussed in detail.

Physical Model and Boundary Conditions
A horizontal-flow flash chamber with three openings, namely, the fluid inlet, fluid outlet, and steam outlet, was studied, as shown in Figure 1a.The chamber length and width were 1 m and 1.4 m, respectively.A structured grid was used in the chamber, as shown in Figure 1b.Four grid sizes were independently tested, and the grid sizes were 3.0, 2.5, 2.0, and 1.5 mm, corresponding to 154,845, 224,000, 350,000, and 620,712 grids, respectively.A grid size of 2.0 mm was selected because of the computational time and clarity of the results.
Appl.Sci.2023, 13, x FOR PEER REVIEW 3 of 17 efficiency.Balaji et al. [20] investigated the effect of non-condensable gases in a low-temperature thermal desalination system and showed that each unit mass of non-condensable gas flowing into the condensation chamber resulted in 1.6 times the unit vapor loss from the vacuum pump outlet.Genthner et al. [21] studied the effect of non-condensable gases in multiple-effect distillation systems and showed that the presence of non-condensable gases increased the heat transfer area of a horizontal tube condenser by about 10-30%, while each unit of non-condensable gases led to a loss of 20-30 units of vapor under severe operating conditions.The aforementioned numerical studies all treated liquid as a single phase.However, in the industry, liquids often include and dissolve a significant amount of non-condensing gases.Therefore, we propose a new coupled numerical model that includes multiphase, multicomponent, phase-change, and turbulent-flow modules that can accurately predict the liquid flash evaporation process with a non-condensing gas.The non-condensing gas was considered to be dissolved and entrained into two components in the liquid.The effects of the two types of non-condensing gases on the flash evaporation mechanism and steam generation rates are discussed in detail.

Physical Model and Boundary Conditions
A horizontal-flow flash chamber with three openings, namely, the fluid inlet, fluid outlet, and steam outlet, was studied, as shown in Figure 1a.The chamber length and width were 1 m and 1.4 m, respectively.A structured grid was used in the chamber, as shown in Figure 1b.Four grid sizes were independently tested, and the grid sizes were 3.0, 2.5, 2.0, and 1.5 mm, corresponding to 154,845, 224,000, 350,000, and 620,712 grids, respectively.A grid size of 2.0 mm was selected because of the computational time and clarity of the results.

Working Conditions and Boundary Conditions
This study investigated water flow with non-condensable CO2.Two types of noncondensable CO2 were considered, namely, the dissolved and entrained components in water [22,23], and were respectively considered as liquid-phase CO2(aq) and gas-phase CO2(g).Before flash evaporation, the water in the flash chamber had an initial water-level height (h) controlled by the static pressure of the liquid outlet (Pout,l).During the working process, the environmental pressure of the flash chamber (controlled by the pressure of

Working Conditions and Boundary Conditions
This study investigated water flow with non-condensable CO 2 .Two types of noncondensable CO 2 were considered, namely, the dissolved and entrained components in water [22,23], and were respectively considered as liquid-phase CO 2(aq) and gas-phase CO 2(g) .Before flash evaporation, the water in the flash chamber had an initial water-level height (h) controlled by the static pressure of the liquid outlet (P out,l ).During the working process, the environmental pressure of the flash chamber (controlled by the pressure of the vapor outlet P out,v ) was lower than the saturation pressure corresponding to the inlet temperature (T in ), and the pressure difference was the driving force of the flash process.
The inlet conditions were water with a fixed temperature (T in ) and mass flow rate.Different masses of carbon dioxide entered the flash chamber together as a gas or solute in water, depending on the operating conditions.The vapor and liquid outlets are both pressure conditions, where the steam outlet (P out,v ) pressure is equal to 19.9 kPa, and the liquid outlet pressure is equal to 19.9 kPa at the water-height pressure (ρgh).The physical parameters of water and vapor at different temperatures are listed in Table 1.The physical properties of water were assumed to have fixed values corresponding to P out,v , and the vapor was assumed to be an ideal gas for calculating its density under different temperatures and pressures.

Model Construction and Governing Equations
In this study, a multiphase, multicomponent, and pressure-driven Lee phase-changecoupled model was proposed [24].The governing equations and calculation flowchart are shown in Figure 2. The left side of the chart shows the main governing equations, including the transient momentum conservation equations, mass conservation equation, energy conservation equation, phase fraction conservation equation (VOF model), and the non-condensing gas fraction equation.The right side shows the calculated equations for the phase change and released non-condensing CO 2 source term, which is the original model of this work.The phase-change ratio in the conventional Lee model is calculated using the temperature difference (T-T sat ), whereas in this study, T sat is not a fixed value and is calculated using the local pressure, as shown in Equations ( 18) and (19).Two types of dissolved components are released (CO 2(aq) ): supersaturation release and interfacial diffusion release, corresponding to Equations ( 20) and ( 23), respectively.
Appl.Sci.2023, 13, x FOR PEER REVIEW 4 of 17 the vapor outlet Pout,v) was lower than the saturation pressure corresponding to the inlet temperature (Tin), and the pressure difference was the driving force of the flash process.
The inlet conditions were water with a fixed temperature (Tin) and mass flow rate.Different masses of carbon dioxide entered the flash chamber together as a gas or solute in water, depending on the operating conditions.The vapor and liquid outlets are both pressure conditions, where the steam outlet (Pout,v) pressure is equal to 19.9 kPa, and the liquid outlet pressure is equal to 19.9 kPa at the water-height pressure (ρgh).The physical parameters of water and vapor at different temperatures are listed in Table 1.The physical properties of water were assumed to have fixed values corresponding to Pout,v, and the vapor was assumed to be an ideal gas for calculating its density under different temperatures and pressures.

Model Construction and Governing Equations
In this study, a multiphase, multicomponent, and pressure-driven Lee phase-changecoupled model was proposed [24].The governing equations and calculation flowchart are shown in Figure 2. The left side of the chart shows the main governing equations, including the transient momentum conservation equations, mass conservation equation, energy conservation equation, phase fraction conservation equation (VOF model), and the noncondensing gas fraction equation.The right side shows the calculated equations for the phase change and released non-condensing CO2 source term, which is the original model of this work.The phase-change ratio in the conventional Lee model is calculated using the temperature difference (T-Tsat), whereas in this study, Tsat is not a fixed value and is calculated using the local pressure, as shown in Equations ( 18) and (19).Two types of dissolved components are released (CO2(aq)): supersaturation release and interfacial diffusion release, corresponding to Equations ( 20) and ( 23), respectively.

Mass, Momentum, and Energy Equations
Equations ( 1)-( 3) show the conservation of mass, momentum, and energy of the two phases, respectively.The calculated pressure, velocity, and temperature are shared between the two phases, and the relative motion of the two phases can be described by variations in the local phase fraction [25,26].
where t represents time, p is the local pressure, v is the local velocity vector, ρ is the density of the mixture, µ is the dynamic viscosity, λ is the thermal conductivity, T is the local temperature, g is the gravitational acceleration vector, E is the internal energy, and .
S lv is the source term of energy.
The density of the vapor is calculated using the ideal gas equation of state, which is updated before each calculation step, as shown in Equation ( 4), where M W is the molar mass of water, taken as 0.018 kg mol −1 , and R is the ideal gas constant of 8.314 J/(mol•K).
The F CSF is the source term of the volume force used to characterize the surface tension, as shown in Equations ( 5) and ( 6).Its representation was derived from a continuum surface force model [27].The model assumes that the surface tension is constant along the surface and models only the normal force.The magnitude of the pressure drop across the surface is obtained from the magnitude of the curvature of the interface κ and the surface tension coefficient σ, and the calculation is added in the form of the source term of the momentum equation by deriving the divergence theorem.

VOF Equations
By solving the partial differential equation of the phase fraction, the VOF method can obtain the percentage of each phase in different cells and show the overall gas-liquid distribution.The volume-fraction equation of the VOF model is expressed as follows [28]: where α l is the volume fraction of the liquid phase, and .
m lv is the rate of transfer of mass in units of kg/(m 3 s).During the calculations, α l is kept in the range of 0 to 1, and the volume fraction of the gas phase and liquid phase must sum to 1 in a finite volume, which provides the distribution of the liquid and gas phases within a single cell.
Furthermore, the magnitude of the internal energy of the unit that mixed the two phases in the energy equation was calculated using the following equation: where µ, ρ, and λ are the dynamic viscosity, density, and thermal conductivity of the mixture, respectively.These physical properties Φ can be calculated by averaging the volume fraction as shown below [29].

Species Transport Equations
The dissolved CO 2(aq) and entrained CO 2(g) in the water were released and dispersed in the flash chamber.To describe the distribution of CO 2 within the flash chamber, component mass fraction equations were added as shown in the following equations [30]: In the above equations, Y CO2(aq) and Y CO2(g) are the mass fractions of CO 2 in the form of solute in water and gas in vapor, respectively.The left side shows the unsteady-state terms and convection terms of the mass fraction, and the right side shows the diffusion terms and interphase transfer source terms, where the source terms are used to characterize the variation in CO 2 during the release of dissolved gases.In the turbulent state, the diffusion flux J i can be expressed using the following formula: In this formula, the first term on the right side is the diffusion flux that is obtained based on Fick's law, which includes the diffusion caused by both the concentration gradient and turbulence.µ t is the turbulent viscosity, D T,i is the thermal mass diffusion coefficient, and Sc t is the turbulent Schmidt number, which is equal to 0.7.In this numerical model, the diffusion coefficient D ij of CO 2 in water was taken as 4.1 × 10 −9 m 2 s −1 [31], and the diffusion coefficient D ij of CO 2 in vapor was taken as 0.26 × 10 −4 m 2 s −1 [32].

Turbulence Equations
In addition, the RNG k-ε model was selected to simulate the turbulent flow in the flash tank.This model has a better result of turbulence simulation compared to the Standard k-ε and RSM models in the VOF model for stratified flow simulation [33,34].This model calculates the turbulent kinetic energy in the flow by solving the k and ε equations.The control equations are shown in Equations ( 14) and (15).
where k represents the turbulent kinetic energy, ε represents the turbulent dissipation rate, G k represents the generation of turbulence kinetic energy due to the mean velocity gradients, and µ eff is the corrected turbulence viscosity.

Phase-Change Lee Equations
Based on the above model, the values of the source terms .m lv and .
S lv in Equations ( 3) and (7) were further varied to simulate the phase-change process.
The transfer of water from the liquid phase to the gas phase during the phase-change process is described by .m lv in the volume-fraction equation, which ultimately leads to a decrease in α l and an increase in α v .The decrease in sensible heat owing to evaporation is controlled by the energy source term .S lv in the energy equation.In this case, the values of the source terms were obtained from the modified Lee model, as shown in Equations ( 16) and (17). .
where c represents the time relaxation parameter, p w is the partial pressure of local water, and r is the latent heat of evaporation.
As shown in the two equations above, this model first determines whether the cell is undergoing a phase transition, and it then calculates the source terms of the phasetransition cell using the local temperature and static pressure.This approach does not require a distinction between nucleation boiling and interfacial mass transfer and uses the same operation to deal with the boiling phenomenon that occurs in both the phase-interface cell and the cells filled with liquid, which leads to a simple and effective simulation of the full-scale flow boiling phenomenon.
Compared to the original Lee model, the main change in the previously described model is that the driving force of the source term is provided by the local pressure and temperature.This modification is important because a small local pressure change at low ambient pressures can have a significant impact on the saturation temperature in that the distribution of water pressure due to gravity can lead to a non-uniform distribution of the flash vapor rate.

Non-Condensing CO 2 Release Model
The release process was divided into two modes, supersaturated precipitation ( .S i ) and interfacial diffusion ( .S s ), as follows: .
The supersaturated release rate ( .m s ) is caused by the supersaturated solubility, which occurs when the mass fraction of CO 2 in the local liquid phase is greater than the corresponding solubility of the local pressure.This is calculated as follows: .
where c s is the supersaturation release coefficient, Y CO 2(aq) is the dissolved gas mass fraction for a single grid, Henry's constant for carbon dioxide dissolved in water k H is 6 × 10 7 pa [35], and p CO2 is the CO 2 partial pressure obtained using Dalton's law as follows: Appl.Sci.2023, 13, 11638 8 of 17 The interfacial diffusion ( .m id ) of CO 2 is caused by the gas-liquid concentration difference at the interface.The process was activated when the local liquid-phase CO 2 concentration exceeded the local gas-phase CO 2 concentration as follows: .
where c i is the interface release coefficient, and A is the interface area.In this study, a simple geometric reconstruction method is employed to obtain the interface area.For a single grid, only the phase fraction and gradient were considered, and the continuity of the interface was ignored.Four typical gas-liquid distribution types in a single grid are illustrated in Figure 3.
Appl.Sci.2023, 13, x FOR PEER REVIEW 8 of 17 where ci is the interface release coefficient, and A is the interface area.In this study, a simple geometric reconstruction method is employed to obtain the interface area.For a single grid, only the phase fraction and gradient were considered, and the continuity of the interface was ignored.Four typical gas-liquid distribution types in a single grid are illustrated in Figure 3.For the four interfacial types, the specific solution processes of the interface area are as follows: (1) the position of the main body of the liquid phase is obtained by the positive and negative values of the gradient components of the phase fraction, and the coordinate system is established with the grid corner of the main body of the liquid phase as the origin; (2) the normal vector and slope of the interface are obtained by the phase fractional gradient; (3) the interface linear equation intercept is calculated using the phase fraction; and (4) finally, the linear equations of the phase interface and the grid intersection are obtained, and the interface area is solved using the Equation (23).
where X cell and Y cell represent the width of the local grid in the x and y directions, respectively, ∇αv,x and ∇αv,y are the local grid liquid fraction gradients in the x and y directions, respectively, and C x and C y are the intersection points of the interface line with the x-axis and y-axis, respectively.

Numerical Method
The PISO method was used to couple the momentum and mass equations.The gradient was obtained using the least-squares method, and the density, momentum, turbulent kinetic energy, turbulent dissipation rate, and energy equations were in a first-order windward discrete format.The time step used for transient calculations was 5×10 −5 s.The under-relaxation factors were used at values 0.3 (pressure), 0.3 (momentum), 0.5 (turbulent kinetic energy and turbulent dissipation rate), 1 (mass fraction), and 0.8 (energy) [36,37].The convergence criterion was set to a residual of less than 10 −6 (energy equation) and 10 −3 (continuity, momentum, turbulence, and CO2 transport equations).
In addition, Ding et al. [38] showed that numerical oscillations may occur when using the Lee model for phase-transition processes, especially in the presence of excessively large evaporation coefficients.To ensure that the calculation had better stability, this study assumed that the local pressure did not change within 1 ms.This assumption effectively increases the convergence of the calculations.For the four interfacial types, the specific solution processes of the interface area are as follows: (1) the position of the main body of the liquid phase is obtained by the positive and negative values of the gradient components of the phase fraction, and the coordinate system is established with the grid corner of the main body of the liquid phase as the origin; (2) the normal vector and slope of the interface are obtained by the phase fractional gradient; (3) the interface linear equation intercept is calculated using the phase fraction; and (4) finally, the linear equations of the phase interface and the grid intersection are obtained, and the interface area is solved using the Equation (23).
where X cell and Y cell represent the width of the local grid in the x and y directions, respectively, ∇α v,x and ∇α v,y are the local grid liquid fraction gradients in the x and y directions, respectively, and C x and C y are the intersection points of the interface line with the x-axis and y-axis, respectively.

Numerical Method
The PISO method was used to couple the momentum and mass equations.The gradient was obtained using the least-squares method, and the density, momentum, turbulent kinetic energy, turbulent dissipation rate, and energy equations were in a first-order windward discrete format.The time step used for transient calculations was 5 × 10 −5 s.The under-relaxation factors were used at values 0.3 (pressure), 0.3 (momentum), 0.5 (turbulent kinetic energy and turbulent dissipation rate), 1 (mass fraction), and 0.8 (energy) [36,37].The convergence criterion was set to a residual of less than 10 −6 (energy equation) and 10 −3 (continuity, momentum, turbulence, and CO 2 transport equations).
In addition, Ding et al. [38] showed that numerical oscillations may occur when using the Lee model for phase-transition processes, especially in the presence of excessively large evaporation coefficients.To ensure that the calculation had better stability, this study assumed that the local pressure did not change within 1 ms.This assumption effectively increases the convergence of the calculations.

Model Verification
The results obtained from the numerical model constructed in this study were compared with experimental data [39], and the release rates (c s and c id ) were determined based on the experimental data.The experimental test shown in Figure 4a illustrates a cylindrical vessel with an inner diameter of 0.045 m, a height of 0.2 m, and a water height of 0.15 m.The experimental operation used dissolved saturated CO 2 water from 5.4 bar decreasing to three pressure cases (5.31, 5.23, and 5.08 bar).The experimental release rate at 100 s was selected and compared with the numerical results, as shown in Figure 4b, and the maximum error was less than ±10%.The release coefficients of supersaturation CO 2 (c s ) and interface diffusion CO 2 (c id ) were 8.8 × 10 −6 (1/s) and 0.01 (1/m 2 •s), respectively.

Model Verification
The results obtained from the numerical model constructed in this study were compared with experimental data [39], and the release rates (cs and cid) were determined based on the experimental data.The experimental test shown in Figure 4a  In addition, the flash evaporation process without the presence of non-condensable gases was also compared with the experimental results [40].The flash chamber for the flow flash evaporation used in the experimental articles compared had an internal flash space of 480 mm in height, 100 mm in length, and 100 mm in width.Five sets of numerical simulation experiments with different inlet temperatures were carried out under the conditions of an initial water-level height of 120 mm and flow rate of 800 L/h.The comparison results are shown in Figure 5 below, and the steam generation rate obtained from the numerical simulation shows a good accuracy.In addition, the flash evaporation process without the presence of non-condensable gases was also compared with the experimental results [40].The flash chamber for the flow flash evaporation used in the experimental articles compared had an internal flash space of 480 mm in height, 100 mm in length, and 100 mm in width.Five sets of numerical simulation experiments with different inlet temperatures were carried out under the conditions of an initial water-level height of 120 mm and flow rate of 800 L/h.The comparison results are shown in Figure 5 below, and the steam generation rate obtained from the numerical simulation shows a good accuracy.
gases was also compared with the experimental results [40].The flash chamber for the flow flash evaporation used in the experimental articles compared had an internal flash space of 480 mm in height, 100 mm in length, and 100 mm in width.Five sets of numerical simulation experiments with different inlet temperatures were carried out under the conditions of an initial water-level height of 120 mm and flow rate of 800 L/h.The comparison results are shown in Figure 5 below, and the steam generation rate obtained from the numerical simulation shows a good accuracy.

Effect of Entrained CO 2 Concentration on Flash Performance
The 63.5 • C water with four kinds of entrained CO 2 volume concentrations (1 vol%, 4 vol%, 7 vol%, and 10 vol%) was compared.The steam generation rate and CO 2 release rate after flash process stabilization are shown in Figure 6.Compared with pure water, the steam generation rate for water with entrained CO 2 increased by approximately 4.8%, 11.7%, 15.9%, and 20.0%, corresponding to volumes of 1%, 4%, 7%, and 10%, respectively.The increase in the steam generation rate is proportional to the increase in the steam-carrying volume.Compared with pure water, the presence of entrained CO 2 leads to an increase in the steam evaporation rate per step.In addition, the percentage of CO 2 components in the outlet steam of the flash chamber increased significantly with the inlet flow of entrained CO 2 and were approximately 0.13%, 0.50%, 0.84%, and 1.16%, respectively.

Effect of Entrained CO2 Concentration on Flash Performance
The 63.5 °C water with four kinds of entrained CO2 volume concentrations (1 vol%, 4 vol%, 7 vol%, and 10 vol%) was compared.The steam generation rate and CO2 release rate after flash process stabilization are shown in Figure 6.Compared with pure water, the steam generation rate for water with entrained CO2 increased by approximately 4.8%, 11.7%, 15.9%, and 20.0%, corresponding to volumes of 1%, 4%, 7%, and 10%, respectively.The increase in the steam generation rate is proportional to the increase in the steam-carrying volume.Compared with pure water, the presence of entrained CO2 leads to an increase in the steam evaporation rate per step.In addition, the percentage of CO2 components in the outlet steam of the flash chamber increased significantly with the inlet flow of entrained CO2 and were approximately 0.13%, 0.50%, 0.84%, and 1.16%, respectively.The temperature, phase, CO2 concentration, partial pressure, absolute pressure, and local steam production rate distribution are shown in Figure 7, under the 1 vol% entrained CO2 volume concentration.As shown in Figure 7a, there is an area of high temperature above the inlet area, which is caused by non-condensable gas.Under the effect of gravity, the CO2 entrained in the inlet water quickly floats upwards after entering the flash chamber and rapidly drives the water upwards, leading to an increase in water temperature above the inlet region.As can be seen from the distribution of local vapor generation rates in Figure 7, the regions of intense phase change include the low-pressure region at the top and the region near the inlet where carbon dioxide is present in large quantities.The increase in the flash steam generation rate was caused primarily by the high temperature and low water vapor partial pressure around the entrance of the flash chamber.The entrained CO2 gas in the flash chamber had little effect on the total pressure stratification The temperature, phase, CO 2 concentration, partial pressure, absolute pressure, and local steam production rate distribution are shown in Figure 7, under the 1 vol% entrained CO 2 volume concentration.As shown in Figure 7a, there is an area of high temperature above the inlet area, which is caused by non-condensable gas.Under the effect of gravity, the CO 2 entrained in the inlet water quickly floats upwards after entering the flash chamber and rapidly drives the water upwards, leading to an increase in water temperature above the inlet region.As can be seen from the distribution of local vapor generation rates in Figure 7, the regions of intense phase change include the low-pressure region at the top and the region near the inlet where carbon dioxide is present in large quantities.The increase in the flash steam generation rate was caused primarily by the high temperature and low water vapor partial pressure around the entrance of the flash chamber.The entrained CO 2 gas in the flash chamber had little effect on the total pressure stratification phenomenon of the flash tank; however, the local high CO 2 concentration reduced the partial pressure of the steam in the air space near the inlet.Under the combined effect of the low partial pressure of steam and the high temperature at the inlet, there was a severe flash area caused by entrained CO 2 near the inlet area.Owing to the high steam generation rate in this area, the overall flash steam generation rate of the flash chamber also increased.The phase distribution and local CO2 release rate distribution are shown in Figures 8  and 9.Under the condition of a higher inlet entrained CO2 volume fraction, the phasechange process in the intense flash area is slightly changed because the increase in the inlet flow rate of entrained CO2 increased local CO2 concentration; however, the increase in the inlet flow rate also intensified the bubble uplift process, thus narrowing the area of the intense flashing area around the inlet.In general, the more inlet-entrained CO2 shrank the region of the intense flash area, but increasing the intensity of the phase transition, and had little effect on the gas-liquid distribution.The phase distribution and local CO 2 release rate distribution are shown in Figures 8 and 9.Under the condition of a higher inlet entrained CO 2 volume fraction, the phase-change process in the intense flash area is slightly changed because the increase in the inlet flow rate of entrained CO 2 increased local CO 2 concentration; however, the increase in the inlet flow rate also intensified the bubble uplift process, thus narrowing the area of the intense flashing area around the inlet.In general, the more inlet-entrained CO 2 shrank the region of the intense flash area, but increasing the intensity of the phase transition, and had little effect on the gas-liquid distribution.

Effect of Dissolved CO2 Concentration on Flash Performance
Figure 10 shows the steam and CO2 generation rates after stabilization for different dissolved CO2 mass concentrations.Compared with pure water, the steam generation rate increased by 2.1%, 4.8%, 6.2%, and 7.5%, corresponding to dissolved CO2 mass concentrations of 0.0661%, 0.1004%, 0.1348%, and 0.1688% (g/g), respectively.The ratio of the CO2 release rate to the steam generation rate was 0.011%, 0.017%, 0.024%, and 0.030%, respectively.The steam generation rate and CO2 release rate both increased linearly with increasing dissolved CO2 mass concentration.Under different working conditions, the interfacial release CO2 rate was approximately 1.8 times higher than the supersaturated dissolution release CO2 rate.

Effect of Dissolved CO2 Concentration on Flash Performance
Figure 10 shows the steam and CO2 generation rates after stabilization for different dissolved CO2 mass concentrations.Compared with pure water, the steam generation rate increased by 2.1%, 4.8%, 6.2%, and 7.5%, corresponding to dissolved CO2 mass concentrations of 0.0661%, 0.1004%, 0.1348%, and 0.1688% (g/g), respectively.The ratio of the CO2 release rate to the steam generation rate was 0.011%, 0.017%, 0.024%, and 0.030%, respectively.The steam generation rate and CO2 release rate both increased linearly with increasing dissolved CO2 mass concentration.Under different working conditions, the interfacial release CO2 rate was approximately 1.8 times higher than the supersaturated dissolution release CO2 rate.

Effect of Dissolved CO 2 Concentration on Flash Performance
Figure 10 shows the steam and CO 2 generation rates after stabilization for different dissolved CO 2 mass concentrations.Compared with pure water, the steam generation rate increased by 2.1%, 4.8%, 6.2%, and 7.5%, corresponding to dissolved CO 2 mass concentrations of 0.0661%, 0.1004%, 0.1348%, and 0.1688% (g/g), respectively.The ratio of the CO 2 release rate to the steam generation rate was 0.011%, 0.017%, 0.024%, and 0.030%, respectively.The steam generation rate and CO 2 release rate both increased linearly with increasing dissolved CO 2 mass concentration.Under different working conditions, the interfacial release CO 2 rate was approximately 1.8 times higher than the supersaturated dissolution release CO 2 rate.Figure 11 presents the distributions of the CO2 concentration, CO2 interfacial release rate, CO2 supersaturated release rate, and local steam generation rate.Figure 11a shows that the volume fraction of CO2 in gas in the entire flash space is small, which has little effect on the partial pressure of water vapor and did not cause a significant effect on the rate of vapor generation.In addition, as shown in Figure 11b-d, the release rate of CO2 was three to four orders of magnitude lower than the steam generation rate, thus the flow field disturbance by the release of CO2 is not significant.The phase distributions for the different dissolved CO2 mass fractions are shown in Figure 12.The phase distributions for different cases were slightly different owing to the release of CO2.Compared with the pure-water flashing condition, the intense flashing area in the flashing chamber increased, but the change was small.Compared with the flash evaporation process of pure water without CO2, the release of dissolved CO2 causes a low local partial pressure of water vapor in the bottom water body, resulting in denser bubbles of the mixed components of water vapor and CO2 during slow flash evaporation.This causes the steam generation rate to increase slightly and leads to the junction between the severe flash region and slow flash region being not very obvious, especially for a high Figure 11 presents the distributions of the CO 2 concentration, CO 2 interfacial release rate, CO 2 supersaturated release rate, and local steam generation rate.Figure 11a shows that the volume fraction of CO 2 in gas in the entire flash space is small, which has little effect on the partial pressure of water vapor and did not cause a significant effect on the rate of vapor generation.In addition, as shown in Figure 11b-d, the release rate of CO 2 was three to four orders of magnitude lower than the steam generation rate, thus the flow field disturbance by the release of CO 2 is not significant.Figure 11 presents the distributions of the CO2 concentration, CO2 interfacial release rate, CO2 supersaturated release rate, and local steam generation rate.Figure 11a shows that the volume fraction of CO2 in gas in the entire flash space is small, which has little effect on the partial pressure of water vapor and did not cause a significant effect on the rate of vapor generation.In addition, as shown in Figure 11b-d, the release rate of CO2 was three to four orders of magnitude lower than the steam generation rate, thus the flow field disturbance by the release of CO2 is not significant.The phase distributions for the different dissolved CO2 mass fractions are shown in Figure 12.The phase distributions for different cases were slightly different owing to the release of CO2.Compared with the pure-water flashing condition, the intense flashing area in the flashing chamber increased, but the change was small.Compared with the flash evaporation process of pure water without CO2, the release of dissolved CO2 causes a low local partial pressure of water vapor in the bottom water body, resulting in denser bubbles of the mixed components of water vapor and CO2 during slow flash evaporation.This causes the steam generation rate to increase slightly and leads to the junction between the severe flash region and slow flash region being not very obvious, especially for a high The phase distributions for the different dissolved CO 2 mass fractions are shown in Figure 12.The phase distributions for different cases were slightly different owing to the release of CO 2 .Compared with the pure-water flashing condition, the intense flashing area in the flashing chamber increased, but the change was small.Compared with the flash evaporation process of pure water without CO 2 , the release of dissolved CO 2 causes a low local partial pressure of water vapor in the bottom water body, resulting in denser bubbles of the mixed components of water vapor and CO 2 during slow flash evaporation.This causes the steam generation rate to increase slightly and leads to the junction between the severe flash region and slow flash region being not very obvious, especially for a high dissolved CO 2 concentration.The liquid height showed little increase, and a larger liquid-gas interface was produced, which also led to an increase in steam generation.
Appl.Sci.2023, 13, x FOR PEER REVIEW 14 of 17 dissolved CO2 concentration.The liquid height showed little increase, and a larger liquidgas interface was produced, which also led to an increase in steam generation.

Conclusions
In this study, a multiphase, multicomponent, and pressure-driven Lee phase-change model was constructed to numerically study water flash evaporation with non-condensing CO2.The numerical model can effectively predict the flash evaporation process, and an intense flash evaporation area around the inlet and a reduction in the partial pressure of vapor caused by the release of CO2 in the flash chamber were observed using this model, which makes the presence of CO2 enhance the intensity of vapor generation.Overall, the increase ratio of steam generation can reach 4.8-20% and 2.1-7.5% for a 1-10% volume concentration of dissolved CO2 and 0.0661-0.1688%mass concentration of entrained CO2, respectively.

Figure 1 .
Figure 1.Computational domain and grid structure of the flash chamber: (a) computational domain; (b) grid structure.

Figure 1 .
Figure 1.Computational domain and grid structure of the flash chamber: (a) computational domain; (b) grid structure.

Figure 2 .
Figure 2. Flow chart of the numerical model calculation.

Figure 3 .
Figure 3. Four typical gas-liquid distribution types in a single grid: (a) C x < X cell , C y < Y cell ; (b) C x > X cell , C y < Y cell ; (c) C x < X cell , C y > Y cell ; (d) C x > X cell , C y > Y cell .

Figure 4 .
Figure 4. Validation of the numerical model with experimental CO2 release rate data: (a) schematic diagram of the experimental model Re-printed with permission from Ref. [39].2014, Elsevier Inc.; (b) comparison of the experimental and numerical results.

Figure 4 .
Figure 4. Validation of the numerical model with experimental CO 2 release rate data: (a) schematic diagram of the experimental model Re-printed with permission from Ref. [39].2014, Elsevier Inc.; (b) comparison of the experimental and numerical results.

Figure 5 .
Figure 5. Validation of the numerical model with experimental CO 2 release rate data: (a) schematic diagram of the experimental model Re-adapted with permission from Ref. [40].2013, Elsevier Ltd.; (b) comparison of the experimental and numerical results.

17 Figure 5 .
Figure 5. Validation of the numerical model with experimental CO2 release rate data: (a) schematic diagram of the experimental model Re-adapted with permission from Ref. [40].2013, Elsevier Ltd.; (b) comparison of the experimental and numerical results.

Figure 6 .
Figure 6.Variation in the steam generation rate and CO2 release rate of different entrained CO2 volume fractions.

Figure 6 .
Figure 6.Variation in the steam generation rate and CO 2 release rate of different entrained CO 2 volume fractions.

Figure 7 .
Figure 7. Characteristic distribution of the flash chamber for the 1 vol% entrained CO2 volume fraction.

Figure 7 .
Figure 7. Characteristic distribution of the flash chamber for the 1 vol% entrained CO 2 volume fraction.

Figure 9 .
Figure 9. Local CO 2 release rate distribution for different entrained CO 2 volume fractions.

Figure 10 .
Figure 10.Variation in the steam generation rate and CO2 release rate of different dissolved CO2 mass fraction values.

Figure 11 .
Figure 11.Distribution of CO2 release and steam generation.

Figure 10 .
Figure 10.Variation in the steam generation rate and CO 2 release rate of different dissolved CO 2 mass fraction values.

17 Figure 10 .
Figure 10.Variation in the steam generation rate and CO2 release rate of different dissolved CO2 mass fraction values.

Figure 11 .
Figure 11.Distribution of CO2 release and steam generation.

Figure 11 .
Figure 11.Distribution of CO 2 release and steam generation.

Figure 12 .
Figure 12.Phase distribution for different dissolved CO2 mass fractions.

Table 1 .
Physical properties of saturated water and saturated steam at 19.9 kPa.

Table 1 .
Physical properties of saturated water and saturated steam at 19.9 kPa.