Next Article in Journal
A Three-Dimensional Elastoplastic Constitutive Model for Geomaterials
Previous Article in Journal
Food Classification and Meal Intake Amount Estimation through Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of Steam–CO2 Mixture Condensation over a Flat Plate Based on the Solubility of CO2

1
Engineering Research Center of BEEE, Ministry of Education, School of Energy and Environment, Southeast University, Nanjing 210096, China
2
Department of Mechanical Engineering, City University of Hong Kong, Hong Kong 999077, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(9), 5747; https://doi.org/10.3390/app13095747
Submission received: 22 March 2023 / Revised: 25 April 2023 / Accepted: 3 May 2023 / Published: 6 May 2023
(This article belongs to the Topic Advanced Heat and Mass Transfer Technologies)

Abstract

:
In order to successfully study the condensation and separation of a steam–CO2 mixture, a boundary layer model was applied to the mixture condensation of steam and CO2 on horizontal and vertical plates. The modified condensation boundary layer model of steam and CO2, given the CO2 solubility in the condensate, was established, numerically solved, and verified with existing experimental data. Different condensation data of steam–air and steam–CO2 mixtures were compared, and the effect of CO2 solubility on the mixed gas condensation was analyzed under multiple pressure conditions (1 atm–10 MPa). The simulation data show that the presence of CO2 will deteriorate the condensation heat transfer, just like air. Given that CO2 is slightly soluble, some CO2 can pass through the gas–liquid interface to enter the condensate film and reduce the accumulated CO2 on the gas–liquid interface, which improves the condensation. However, the solubility of CO2 is only significant under high-pressure conditions, inducing its effects on condensation. A comparison of the condensation coefficients of the steam–CO2 mixture shows the lower impact of CO2 condensation on the horizontal plate compared to that on the vertical plate. For most conditions, the steam–CO2 mixture gas condensation heat transfer coefficient on the vertical plate surface is still larger than that on the horizontal plate surface, and the improvement in the condensation heat transfer coefficient caused by low CO2 solubility (2 or 10%) at 10 MPa on the vertical plate is also larger than that of the horizontal plate.

1. Introduction

The carbon dioxide emission from the industrial sector and power plants has contributed significantly to global climate change for several decades [1]. The rise in the average global temperature is one of the known impacts of emissions of greenhouse gases (mainly CO2) [2]. In 2021, the Intergovernmental Panel on Climate Change (IPCC) alleged that existing technologies, including Carbon Capture and Storage (CCS), are essential to achieve CO2 reduction and climate change mitigation [3]. Climate change mitigation and adaptation measures likely need to include technologies that can reduce carbon emissions at fossil fuel production sites, capturing the carbon before it is released into the air. Water vapor and CO2 are the main components of combustion products (flue gas).
Chen et al. [4] proposed a liquefied natural gas (LNG) and O2 combustion gas–steam mixed cycle (GSMC) power generation scheme in which peak shaving, cold energy storage, and full CO2 capture are all integrated. In GSMC, steam and CO2 mixtures are condensed and separated in condensers and low/high-pressure heaters, with the condensation pressure exceeding 10 MPa in vacuum conditions. For the transformation of CO2 capture in existing coal-fired power plants, Zhao et al. [5] studied the efficient power generation model with a high CO2 capture rate, and the range of its steam–CO2 mixture condensation and separation pressure also ranges from around 30 kPa to 10 Mpa. In a classical chemical cycle combustion process, fuel is oxidized by metal oxides in a fuel reactor. The exhaust stream, primarily comprising CO2 and H2O, can also be purified to mainly CO2 by condensing the H2O [6]. Thus, steam–CO2 mixture condensation may appear in wide areas and state ranges for component separation and CO2 capture [4,5,6].
The condensation of water vapor with air, N2, O2, CH4, H2, He, NH3, or other gases has been widely studied and summarized by scholars. Yi et al. [7] experimentally studied the effect of non-condensable gas (NCG) on the condensation heat transfer of water vapor with a 0.05–0.5 air mass fraction on the isothermal vertical aluminum plate. They also discussed changes in the condensation mode with surface subcooling and air concentration. Kuhn et al. [8] obtained an extensive database of steam–air mixtures and steam–helium mixtures and investigated the local heat transfer of steam in the presence of NCGs in vertical tubes. Rose [9] measured the heat transfer of film condensation on vertical plate surfaces, and the approximate theoretically based equations were obtained in relation to the mass flux of vapor to the condensing surface (condensation rate) to the free-stream and condensate surface conditions. Ilpo et al. [10] conducted condensation air and water vapor experiments on vertical and horizontal plates. Hie et al. [11] conducted a vapor condensation experiment with air as an NCG on a horizontal plate. The above experimental results all show that the existence of a small amount of NCG can greatly affect the condensation heat transfer process. In the laminar flow regime, Diwany et al. [12] found that a 1.5% volume fraction of air will bring about a 50% decrease in the condensation heat transfer coefficient of steam. The condensation features of a non-azeotropic mixture are different from those of steam–NCG mixtures. Morrison et al. [13] analyzed the condensation heat transfer characteristics of a steam/NH3 mixture and considered the effect of gas concentration on condensation.
Sparrow et al. [14,15] and Minkowycz et al. [16] developed a boundary layer theory for the condensation of a steam–air mixture on horizontal and vertical plate surfaces, which is a predictive model that is easy to both build and compute. The boundary layer theory can reveal the condensation mechanism of a steam–NCG mixture with relative clarity, and several studies based on this theory have been carried out [17,18]. In recent decades, with the development of commercial computing software, many scholars have carried out simulation analysis with computational fluid dynamics (CFD) software on fluid flow fields and obtained expected research results. Angelino et al. [19] used the self-similar methodology to describe laminar momentum diffusion, and Boghi et al. [20,21] described laminar passive scalar diffusion with the self-similar methodology. Furthermore, the theory of laminar momentum diffusion is developed, and a wide range of Reynolds and Schmidt studies offered generality to the results, in turn providing theoretical support for the experimental results. Hossein et al. [22] used the finite volume method to simulate the fuel cell processes and introduced innovative designs to improve the performance of cylindrical polymer fuel cells. Moreover, the condensation heat transfer of a steam–NCG mixture could also be calculated using CFD software. Dehbi et al. [23] integrated a model for the wall condensation of a steam–NCG gas mixture in the ANSYS CFD code FLUENT, and the predicted heat transfer rate was found to be generally consistent with experimental results. Li [24] simulated the condensation of water vapor in the presence of NCGs in a vertical cylindrical condenser and proposed a method to address the limitations of the currently used CFD to solve the simulation of flows involving multispecies as well as gas fluids and liquids in separate channels. Choudhury et al. [25] studied the unsteady hydrodynamic-free convective flow of a viscoelastic incompressible fluid past a vertical porous plate in the presence of a variable suction and provided approximate solutions on velocity and temperature fields, shear stress, and the rate of heat transfer using the perturbation technique. Based on a suitable model, CFD simulation can provide additional details on condensation, which may be more accurate, but more modeling and calculation work is also required. Therefore, the boundary layer theory, which only needs to solve the ordinary differential equation (group), still has analytical and engineering significance and can be used as a simple method to study the condensation of a mixture.
Compared with other steam–NCG mixtures, research on the condensation of steam–CO2 mixtures is very limited. Ge et al. [26,27] studied the condensation of water vapor on different surfaces with high and low concentrations of CO2. Heat transfer data were different between the steam–CO2 and steam–air mixture condensation rates in the experiment. Lu et al. [28] investigated steam condensation on a horizontal tube in the presence of CO2 with an experiment and obtained condensation heat transfer coefficients for the sub-atmospheric pressure ranging from 5 kPa to 101 kPa. Takami et al. [29] used Comsol software to simulate the condensation heat transfer and separation of a steam–CO2 and O2 mixture in the condenser. This provides some data that serve as the basis for the design of a volumetric heat exchanger (condenser) that is suitable for steam–CO2 separation in flue gas, thus improving the efficiency of power plants. Lu et al. [30] studied vapor condensation in the presence of non-condensable gases using CFD simulations and calculated the effects of velocity, surface subcooling, and the molar fraction of non-condensable gases on the heat transfer from steam–air or steam–CO2 mixtures. However, there is still a gap in the recent study to address before we can meet the engineering application requirements in steam–CO2 mixture condensation.
Based on current studies, CO2 is usually regarded as an NCG. However, in an engineering application, the solubility of CO2 is not only different from NH3, which is miscible with water in any proportion, but also different from NCGs, such as air, He, etc. CO2 is a slightly soluble gas. Based on the phase equilibrium data of H2O and CO2 [31,32,33,34], under normal pressure or vacuum conditions, the solubility of CO2 is small. However, when the pressure is much higher than atmospheric conditions (such as P > 1 MPa), the solubility of CO2 significantly increases, which may cause some obvious effects on condensation. Studies that focus on these various areas are rarely conducted. Thus, the effect of low CO2 solubility on mixture condensation heat transfer is unclear and worth studying in order to serve related engineering applications [4,5,6].
Therefore, in this paper, boundary layer analysis was carried out specifically for the condensation heat transfer characteristics of steam–CO2 mixtures in a larger range of parameters that fit engineering application conditions. Using the laminar boundary layer model of Sparrow et al. [14,15] and Rose et al. [9], the mixture condensation boundary layer model based on CO2 solubility was further implemented for both horizontal and vertical plates. The condensation heat transfer of steam–CO2 and steam–air mixtures was simulated under normal pressure conditions at the same mass and mole fraction, respectively, for data comparison purposes. Moreover, in the total pressure range of 1 atm–10 MPa, which aims to cover the highest pressure range that is needed for the engineering applications mentioned above [4,5,6], the impact of low CO2 solubility on the condensation heat transfer of steam–CO2 mixtures was analyzed. Through simple laminar boundary layer analysis, the heat and mass transfer characteristics of steam–CO2 mixture condensation can be preliminarily studied, which acts as the basis for the turbulent condensation heat transfer study of a steam–CO2 mixture and can also offer some reference for the design of CO2 separation and capture equipment.

2. Boundary Layer Models of Steam–CO2 Mixture Condensation

2.1. Steam–CO2 Condensation Model Based on CO2 as an NCG

If CO2 is regarded as an NCG, Sparrow’s model [14,15] can be applied for steam–CO2 condensation on horizontal and vertical plate surfaces by introducing the properties of CO2 instead of air. It is abbreviated as an NG model for short in this paper.
Figure 1a,b show a schematic diagram of the mixture condensation boundary layer model. The gas temperature and CO2 mass fraction at the entrance and the infinity of the plate are T and α, respectively. The wall temperature (Tw) is constant and the condensate takes a liquid film on the plate with a thickness of δ. Symbols for the gas–liquid interface are all denoted by subscript i, while subscript x indicates the local parameters. The liquid film covered on the horizontal plate is driven by gas, which has an inlet velocity of U, while the liquid film covered on the vertical plate is driven by gravity, and the gas velocity at infinity is 0 m/s. The equations and boundary conditions of the NG model governed by Sparrow et al. [14,15] are arranged and shown in Table 1.
As for the governing equation of boundary layer theory, Sparrow [14,15] introduced dimensionless intermediate variables. Based on dimensionless simplifications, the NG model with partial differential equations can be reduced to ordinary differential equations, and this transformation makes the equation easier to solve numerically, as shown in Table 2 below.
In addition, the relationship of steam partial pressure PH2O and total pressure P is given in Equation (8) [14] for low-pressure conditions to complete the iterative calculation.
P H 2 O = 1 α i x 1 α i x ( 1 M H 2 O / M CO 2 ) P
Furthermore, in the NG model, in the liquid-phase boundary layer, the condensation mass flux (m) can be represented as:
m = ρ l ( u l d δ d x v l ) = ρ v ( u v v v )
combined with Equations (4a) and (4c) and Equations (5a) and (5c), we can obtain the following:
m = 1 2 ν l U x f ( η i ) = 1 2 ν v U x F ( 0 )
Under the assumption of NCGs, CO2 entering the condensate through the gas–liquid interface is 0:
ρ CO 2 ( u v d δ d x v v ) j CO 2 y = δ = 0
where j is the diffusive mass flux in the boundary layer, and then αix can be given as Equation (5c).

2.2. Model Modification Considering CO2 Solubility in the Condensate

Figure 1c shows the difference in mixture condensation under the NCG and solubility assumptions of CO2. The solubility of CO2 in water (S) has complete data [35,36,37]. Equation (12) gives a method that is used to calculate S in this paper, which is a correlation of CO2 solubility under different temperatures and pressures [35].
P a r = c 1 + c 2 T + c 3 / T + c 4 T 2 + c 5 / T 2 + c 6 P + c 7 P T + c 8 P / T + c 9 P T 2 + c 10 P 2 T + c 11 P 3 S M C O 2 = P P H 2 O exp P a r R CO 2
where Par depends on T and P, and Equation (12) is an empirical formula showing piecewise functions, where T and P are different values, and the coefficients will be different. Moreover, coefficients c1~c11 can be acquired from the study conducted by Mao et al. [35], RCO2 is the fugacity coefficient of CO2, and bar is the unit of pressure in Equation (12).
Based on the solubility data, CO2 is a soluble gas in the pressure range of 101.325 kPa–10 MPa, and the temperature is calculated in the range of 345.15 K–583.2 K; thus, S is in the range of 0.3545–63.12 g/kg(water), which increases with pressure and decreases with temperature. Given the solubility, the low solubility of CO2 in the condensate heat transfer of steam–CO2 cannot be simply ignored.
When considering that CO2 is slightly soluble, the expression of mass transfer at the gas–liquid interface changes, which may affect the solution results of the heat transfer. It is necessary to modify the boundary layer model in Section 2.1. In this paper, the steam–CO2 condensation model based on low CO2 solubility is referred to as the low-solubility condensable gas (CG) model.
In the CG model, the compositions of the condensate are different. The amount of CO2 dissolved into the condensate through the gas–liquid interface, f(αl), is related to the distribution of CO2 in the condensate, αl. Equation (11) in the NG model becomes:
ρ C O 2 ( u v d δ d x v v ) j CO 2 y = δ = f α l
The ionization of CO2 in water [32,38] and the condensation and upstream CO2 dissolution can all influence the αl distribution in the condensate. However, to simplify the analysis, it is assumed that the change in the CO2 percentage in the near-thin film region is small, which is:
α l x = 0 , α l y = 0
Thus, αl is constrained by:
2 α l y 2 = 0 , y = 0 : α l y = 0 , y = δ : α l = S 1 + S
In Equation (15), the form of the conservation equation is the same as Equations (2a) and (3a), but the boundary conditions are different from those of momentum and energy equations. The solution of Equation (15) is shown below:
α l = S 1 + S ( 0 y δ )
This uniform distribution of CO2 at the y direction is obtained under simplifications. However, this is acceptable because the liquid film is thin and the condensate flow rate is small. To form the CO2 distribution in Equation (16), the following can be obtained:
f α l = m S
This expression is similar to the condensation of the steam–NH3 mixture, in which the condensation proportion of ammonia and steam is also determined by the phase equilibrium [39]. Then, in the CG model, Equation (17) can be modified as:
[ ρ CO 2 ( u v d δ d x v v ) j CO 2 ] y = δ = m S
By combining Equations (5a), (5c), and (10), the vapor mass percentage at the gas–liquid interface becomes:
α i x = α + 1 / 2 S c F ( 0 ) φ v ( 0 ) S 1 + 1 / 2 S c F ( 0 ) φ v ( 0 )
When considering the CO2 slight solubility, the situation is different from Equation (7c), and it becomes:
φ v ( 0 ) = 3 S c F ( 0 ) [ φ v ( 0 ) + α S ]
In addition, the solution enthalpy of CO2, hso, is rarely mentioned in the literature, and some of the results on this [40,41] are listed in Table 3. The data in the hso∙S/hfg column show the ratio of CO2 solution enthalpy to the latent heat of steam condensation, which is smaller than 0.5% for the states listed in Table 3. Given that the solution enthalpy data are incomplete and its proportion is very low, hso/hfg is taken as a constant (0.2) for the estimation of solution enthalpy.

2.3. Numerical Method

Through the introduction of the above model, we transformed the partial differential equation, which is not easy to solve into the ordinary differential equation and easier to solve through the intermediate variable. Different from the integral method adopted by Sparrow [14,15], this paper uses a discrete case to solve the equation. This method expands the computational scope of the model and can find the convergence solution quickly. By solving the NG or CG model, the gas–liquid interface data Tix and αix can be obtained.
The solution process is shown in Figure 2, and the fourth-order Runge–Kutta (RK4) and Euler iterative methods are respectively used to calculate and solve the ordinary differential equations on the horizontal and vertical plate surfaces. The process is to assume the physical property of the gas–liquid interface at any position of the plate first, solve the assumed actual value in relation to the existing solution conditions described in the NG model or CG model, complete the cycle, and repeat the iterative calculation until the assumed value and the actual value meet the specified deviation, after which it is considered to solve the actual solution of the gas–liquid interface. For the solution of the physical properties of the mixture, a program in MATLAB R2021a was written that calls the NIST 9.1 [42] database to check the physical parameters to be calculated at the same time in the process of iterative calculation.
After Tix and αix are acquired at the gas–liquid interface, condensation heat transfer features can be calculated, as shown below [15].
The average heat transfer coefficient (HTC) on the horizontal plate:
h ¯ = 1 L 0 L h x d x = 1 L 0 L q ( T T w ) d x = 2 k l η i , l U ν l L T i T w T T w
The average HTC on the vertical plate:
h ¯ = 1 L 0 L h x d x = h f g μ l T T w 4 3 [ C p ( T i T w ) h f g P r ] 3 / 4 ( g L 3 4 ν l 2 ) 1 / 4
where L is the length of the flat plate and hx is the local heat transfer coefficient. TTw is the temperature difference between the bulk and the wall. For saturate inlet conditions, this indicates the subcooling of condensation. Ti is the average temperature along the plate.
For condensation heat transfer reduction, q/q0 (subscript 0 is the data for pure streams) on the horizontal plate is:
q q 0 = k l k l , 0 υ l , 0 υ l η i , 0 η i T i T w T T w
On the vertical plate, this is:
q q 0 = [ C p ( T i T w ) C p , 0 ( T T w ) ] 3 / 4 ( h f g μ l h f g , 0 μ l , 0 ) 1 / 4 ( υ l , 0 υ l )

3. Model Verification

3.1. Comparing the Results with Numerical Data and Existing Experimental Data

The NG model is based on Sparrow’s boundary layer theoretical model, but the difference is that the NCG is CO2 (not air), and the iterative calculation scheme used has been modified. Therefore, it is still meaningful to compare the simulated numerical result with Sparrow’s data result. Table 4 lists some data compared with Sparrow’s calculation results:
In the current literature, experimental data on the condensation of steam–CO2 mixture are scant. To verify the rationality of boundary layer simulation calculation as reported by Othmer [43], these data are selected for data verification purposes in this study. Comparison data are listed in Table 5. A data comparison method by Minkowycz et al. [16] is also used, in which the heat transfer reduction (q/q0) is compared with Othmer’s experimental data and the boundary layer simulation data. Since the working fluid is a steam–air mixture, the NG model is applied with the NGC properties of air. As shown in Table 5, deviations between experimental and simulation data are around 6.08~14.68%.
Experimental data for steam–CO2 mixture condensation are severely lacking. Thus, we compare some of the data obtained by Ge et al. [27] with the simulation results in Figure 3. In Ge et al.’s study [27], an experiment was carried out on a 5 × 5 cm vertical plate surface, and the saturated mixture gas had a CO2 percentage of 34.3%, with a total pressure of 101.325 kPa. As shown in Figure 3, the maximum deviation between experimental data and calculation results is 17.75%.

3.2. Data Trends of Steam–CO2 Mixture Condensation

At a total pressure of 1 atm, Figure 4a shows the variation trend in the hx of the steam–CO2 mixture condensation along the horizontal plate under different flow rates and inlet CO2 percentages, at a subcooling rate of 7 K. Figure 4b shows the accumulation of liquid film on the surface of a horizontal plate under different flow rate conditions and α = 5%. The data results show that the liquid film builds and the hx drops sharply at the leading edge of the horizontal plate. However, for the downstream function, both the film thickness and hx curves become flattened along the plate. The CO2 percentage is a significant factor influencing the heat transfer coefficient. While there is a distinguishable effect of mixed gas velocity (in the calculation range) on δ, an increase in the velocity can reduce the accumulation of condensation liquid, thus reducing the thickness of the liquid film. With a decrease in the CO2 percentage and an increase in the gas velocity, the heat transfer coefficient on the surface of the horizontal plate is enhanced. The simulation data trend of condensation on the horizontal plate can be expected.
Figure 4c shows the variation trend of the liquid film thickness and flow velocity of the gas–liquid interface along the vertical plate with different inlet and infinity CO2 percentages. The curve trends in Figure 4 are consistent with the existing study of the mixed gas condensation, further indicating that it is better to use the boundary layer model to analyze the features of the steam–CO2 mixture condensation on the horizontal and vertical plate surface.

4. Results and Discussion

4.1. The Effect of CO2 Solubility on Steam–CO2 Mixture Condensation

The effect of CO2 solubility on the condensation of the steam–CO2 mixture is analyzed by comparing the simulation results of boundary layer models with and without considering CO2 solubility, which have been marked as the CG model and the NG model, respectively. By taking the horizontal plate as an example, Figure 5 gives the average HTC and q/q0 curves versus the subcooling. The inlet CO2 mass percentages shown in Figure 5 are taken at 2%, 30%, and 80%, and the simulation results of CG and NG models are both given and compared under total pressure conditions of 1 atm, 1 MPa, and 10 MPa. Steam–CO2 mixture condensation data curves calculated by CG (solid curve) and NG (dash curve) boundary layer models are different. Although the solid and dash curves are closer to each other at total pressure conditions of 1 atm and 1 MPa, they can be distinguished under 10 MPa conditions. All the distinguishable curves presented in Figure 5 show that the average HTC and q/q0 curves calculated by the CG model are larger than those calculated by the NG model under the same conditions. Moreover, under 10 MPa, increases in the concentration and differences in the average HTC curve between the CG and NG models are found. This indicates that the solubility of CO2 is an enhancement or reduction alleviation effect of the steam–CO2 condensation heat transfer.
The effect of CO2 solubility on steam–CO2 mixture condensation is different under different conditions, partly because of the different CO2 solubility itself, which increases with increasing CO2 partial pressure and decreases with temperature. For the small calculation range of subcooling in this study (≤20 K), the solubility of CO2 is mainly controlled by the total pressure and the CO2 percentage. When the total pressure is 1 atm and 1 MPa, the solid and dash curves are almost coincident in Figure 5, especially for Figure 5a,b; this is because the solubility of CO2 is small under this condition. Thus, it is reasonable to regard CO2 as an NCG under low-pressure conditions. As the pressure increases, the partial pressure of CO2 on the surface of the condensate also gradually increases. At a total pressure of 10 MPa, the CG and NG models show significantly different simulation results. By comparing the curves in Figure 5a,c,e, with an increase in the CO2 concentration, the differences in the average HTC values of the CG and NG models become more significant under the same total pressure levels. When the concentration reaches 80%, a reduction in the condensate heat transfer is significant, and q/q0 is not discussed here.
The effect of CO2 solubility on mixture condensation will influence the average HTC data. Figure 6a shows the average HTC deviations between CG and NG models and Figure 6b shows CO2 solubility values under a total pressure of 10 MPa, in which ΔHTC is:
Δ H T C = h ¯ C G h ¯ N G h ¯ C G
For the inlet condition of 2% CO2, the CO2 solubility is around 0.6508–1.6401 g/kg water, and ΔHTC is the range of 0.218–4.4%; meanwhile, for the inlet condition of 30% CO2, the CO2 solubility is around 10.8774–25.3173 g/kg water and ΔHTC is around 1.4–3.5%; and for the inlet condition of 80% CO2, the CO2 solubility is around 36.1272–39.6 g/kg water and ΔHTC is around 2.117–3.005%. It can be found that with the increase in subcooling, the growth rate of low α is faster. This is because the mixed gas with lower α has better condensation and heat transfer conditions, the condensation mass flux is larger, and ΔHTC may show a faster growth trend. In conclusion, when subcooling increases, the solubility of CO2 reduces the reduction in condensation heat transfer more obviously. The contrast of ΔHTC and S shows an increase in the total pressure and CO2 concentration, and the low solubility of CO2 is not negligible.
Figure 7a,b show the average temperature (Ti) and average mass percentage (αi) curves at the gas–liquid interface under total pressure conditions of 1 atm and 10 MPa. Figure 7c shows a variation in the condensate film thickness along the horizontal plate at a subcooling rate of 7 K. The solubility of CO2 is larger under a total pressure condition of 10 MPa, and its influence on Ti and αi is more obvious. Similarly, CO2 solubility is larger for a 30% CO2 inlet condition than a 2% inlet CO2 at the inlet because of the larger partial pressure, so it also has an influence on Ti and αi. Through Figure 7b, the αi calculated by the CG model is smaller than that calculated by the NG model, which indicates that the mechanism of CO2 can promote condensation in order to reduce the accumulation of CO2 at the gas–liquid interface. In Figure 7c, the effect of CO2 solubility on the accumulation of condensate film on the plate can also be seen, but it is not very remarkable. The thickness of the liquid film in the CG model is greater than that in the NG model.
Therefore, when the total pressure is higher than 1 MPa, it is wrong to regard CO2 as an NCG. It is necessary to consider the solubility of CO2 in order to estimate the heat transfer of steam–CO2 mixture condensation under high-pressure conditions. CO2 solubility has an enhanced effect on condensation.

4.2. Comparison of the Steam–CO2 Mixture Condensation Features with Steam–Air

Figure 8 shows the variation in condensation parameters of steam–CO2 (CG model) and steam–air condensation mixtures at 1 atm and 10 MPa on the surface of a 0.2 m long horizontal plate. The velocity of the mixture is 1 m/s, and the CO2 and air mass percentages at the inlet mixture gas are taken as 2% and 30%. As shown in Figure 8a, for the same mass percentage and other conditions of steam–CO2 and steam–air mixtures, there is a remarkable difference in their condensation heat transfer parameters, which is only caused by the different components of CO2 and air. The gas–liquid interface temperature of steam–CO2 mixture condensation is higher than that of a steam–air mixture. Figure 8c shows that the percentage of CO2 at the gas–liquid interface is higher than that of air. Since the solubility of CO2 is small, its effect on condensation can be negligible in low-pressure conditions. Moreover, 10 MPa is selected as the total pressure of this part for the following data comparison. Figure 8b uses ΔNDTi (the deviation in non-dimensional Ti) to compare the ratio of thermal resistance in the gas phase of NG or CG with air at 10 MPa.
Δ N D T i = ( T air T NG / CG ) ( T NG / CG T air ) T T w = ( T T i T T w ) air ( T T i T T w ) NG / CG
The physical significance of the first term represents the proportion of the gas-phase thermal resistance of steam–air total thermal resistance, and the second term represents the proportion of the gas-phase thermal resistance of steam–CO2 (CG model or NG model) in the total resistance. Both terms are dimensionless parameters. Due to the differences in physical properties, under the same condition in α, the gas-phase thermal resistance of steam–air mixture condensation is greater than that of the steam–CO2 mixture. Given the micro-soluble nature of CO2, the gas-phase thermal resistance calculated by the CG model is lower than that calculated by the NG model, as shown in Figure 8a,b, due to the entry of CO2 into the liquid phase. Figure 8d shows ΔNDαi (the deviation in non-dimensional αi), as given below:
Δ N D α i = α i , air α i , NG / CG ( α i α ) air = 1 ( α i α ) air ( α i α ) NG / CG
ΔNDαi represents the CO2 fraction difference at the steam–CO2 and steam–air gas–liquid interface. In Figure 8d, the difference in ΔNDαi increases with an increase in the CO2 concentration and subcooling rate. Figure 8a–d indicate that the impact of CO2 solubility at low pressures or concentrations is small, and the difference between mixture condensation is mainly reflected by the difference in the physical properties of air and CO2. The physical parameters of the mixture gas considered in the boundary layer model are viscosity, thermal conductivity, specific heat capacity, density, and the diffusion coefficient. Moreover, the low solubility of CO2 at high pressures and concentrations is one of the most influential physical properties.
The results in Figure 8 show that different gas components also have a significant effect on the condensation features, which can explain why some existing correlations cannot be used to predict the condensation of steam–CO2 [27]. Figure 8e,f show the h ¯ and Δ h ¯ of the mixtures of steam–CO2 (both the NG model and the CG model) and steam–air on the horizontal plate surface. The heat transfer coefficient of the steam–CO2 mixture is higher than that of the steam–air mixture. Moreover, the difference in the average HTC becomes more obvious with an increase in the subcooling and concentration of CO2 or air; when the subcooling is 20 K and α is 30%, the Δ h ¯ reaches 16%.

4.3. Comparison of Steam–CO2 Mixture Condensation Features on Horizontal and Vertical Plates

Due to the different inlet velocity settings and driving forces of the liquid film, the boundary layer models of the steam–CO2 mixture condensation on the surfaces of the vertical and horizontal plates are not the same.
Figure 9 shows a comparison of the condensation ΔHTC values at CO2 inlet fractions of 2% and 10% on horizontal and vertical plates at 10 MPa. The result shows that the ΔHTC caused by the solubility of CO2 on a vertical plate is larger than that on a horizontal plate under the same α. Thus, the effect of CO2 solubility on improvements in average HTC values on the vertical plate is larger than those on the horizontal plate. On the vertical plate, the gas phase accounts for the main thermal resistance. However, it becomes smaller when CO2 dissolves in the liquid phase.
Figure 10 shows a comparison of some other calculation results on the condensation on horizontal and vertical plate surfaces. Figure 10a shows the calculation results of the average HTC values at the inlet CO2 mass fractions of 2% and 10% under total pressure conditions of 1 atm and 10 MPa. The average HTC condensation value on the vertical plate is higher than that on the horizontal plate under a very small subcooling regime. This shows the advantages of using gravity to drain the condensate compared with using gas flow.
Figure 10b shows the heat transfer reduction (q/q0) curves of the steam–CO2 mixture condensation, and the label is same as Figure 10a. Condensation on the horizontal plate surface is significantly less affected by CO2 than that on the vertical surface at the same CO2 mass fraction. Due to the influence of gravity, the thickness of the liquid film accumulated in the horizontal plate is thicker than the vertical plate, which can be confirmed by the calculation of the thickness of the local liquid film in Figure 10c. With an increase in the CO2 mass fraction, the condensation heat transfer is reduced and the accumulation of liquid film decreases, so the thickness of the liquid film becomes thinner for larger α conditions in Figure 10c. In addition, with the increase in pressure, an improvement is shown in the reduction in condensation heat transfer, and the condensate quality also increases.
Through the above analysis, at medium and low concentrations, the concentration of condensate is the main factor determining the thermal resistance to condensation on the horizontal plate. As shown in Figure 10c, compared with the vertical plate, the heat transfer reduction in the horizontal plate is not sensitive to changes in CO2 concentration. As the thickness of the liquid film covered by the vertical plate is thinner, the thermal resistance in the liquid phase is smaller than that in the gas phase. Figure 10d shows the ratio of the temperature difference in the gas boundary layer and the liquid film, which represents the relative thermal resistance of both the gas phase and the liquid phase. Figure 10d shows that the thermal resistance in the gas phase is higher than that in the liquid phase on the vertical plate. Thus, a small mass fraction of CO2 will have a significant impact.

5. Conclusions

We carried out the theoretical modeling and simulation of the boundary layer for the condensation of the steam–CO2 mixture on the surface of horizontal and vertical plates. (i) Based on the boundary layer theory, CO2 is regarded as an NCG, and the solubility of CO2 also matters. The boundary layer model is modified to obtain a new CG model. At the same time, a new iterative calculation method is used in the iterative calculation to extend the calculation range of the NG model, including extending the NG model to calculate the steam–air mixed condensation to calculate the steam–CO2 mixed condensation and expanding the calculation range of α and pressure. (ii) The influence of CO2 solubility on the heat transfer of mixed gas condensation is clarified, which adds a consideration factor for the next research study on steam–CO2 condensation heat transfer. Comparing the differences between steam–CO2 and steam–air, besides the difference in gas physical properties, solubility was also the factor that caused these differences. (iii) The effect of CO2 solubility on the condensation heat transfer of mixed gas on the horizontal plate and vertical plate was compared. The main conclusions are as follows:
Due to the slight solubility of CO2, part of CO2 can pass through the gas–liquid interface and enter the liquid film, and the CO2 accumulation at the surface of the liquid film decreases, which reduces the condensation heat transfer. However, this effect is only manifested under high-pressure conditions, where the solubility of CO2 is non-negligible. In this paper, the data show that when the total pressure of the steam–CO2 mixture is higher than 1 MPa, the influence of the slightly soluble CO2 on the condensation of the mixed gas should be considered.
The same mass fraction condition was selected to compare the condensation features with different mixture gases. The differences in condensation heat transfer of steam–CO2 and steam–air mixtures under the same conditions are caused by different physical properties of different components at low pressures and low concentrations. CO2 reduces condensation less compared to air. When the total pressure or concentration is high, the slight solubility of CO2 can also be regarded as a physical property, which has an effect on condensation heat transfer.
Within the calculation range, the effect of the CO2 concentration on the condensation of the horizontal plate surface is smaller than that of the vertical plate surface. However, the effect of CO2 solubility in terms of improvements in the heat transfer coefficient of the vertical plate is better than that of the horizontal plate.
For the NG and CG models of the developed boundary layer theory, the inclination angle β of the plate was introduced to unify one set of models of the horizontal plate and vertical plate to determine the influence of geometric position on condensation heat transfer. Based on the limitations of the model, CFD was used in the following calculation to simulate a new scheme containing a high mass fraction of CO2 in steam–CO2 condensation, and the influence of solubility was further compared. The rationality of simulation calculation is to show the difference between the two calculation schemes.

Author Contributions

Conceptualization, Y.C. and J.W.; Methodology, J.W.; Investigation, Y.J. and H.G.; Writing—original draft, B.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number: 51976033).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Some or all of the data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Nomenclature

cconstant
CPspecific heat (J/(kg·K))
Ddiffusion coefficient (m2/s)
Esolution enthalpy (KJ/mol)
F, fdimensionless function of η and ξ
ggravity (m/s2)
h ¯ average heat transfer coefficient (W/(m2·K))
hfglatent heat of condensation (J/kg)
hsosolution enthalpy of CO2 (kJ/kg)
hxlocal heat transfer coefficient (W/(m2·K))
jdiffusive mass flux (kg/(m2·s))
Llength of the plate (m)
mcondensation mass flux (kg/(m2·s))
Mmolar mass (g/mol or kg/mol)
Ppressure (Pa)
Parparameter depending on T and P
PrPrandtl number
Psisteam partial pressure (Pa)
qheat flux (W/m2)
Rfugacity coefficient of CO2
Ssolubility of CO2 (kg/kg water or g/kg water)
ScSchmidt number
Ttemperature (K)
u, vvelocity (m/s)
Uinlet velocity (m/s)
x, ydirection or coordinate (m)
Greek symbols
αmass fraction of CO2
δfilm thickness (mm or m)
ηdimensionless liquid coordinate
θdimensionless temperature
μdynamic viscosity (Pa·s)
ξdimensionless gas coordinate
νkinematic viscosity (m/s2)
ρdensity (kg/m3)
φconcentration function
ψstream function
Acronyms
CGcondensable gas model
ΔHTCaverage HTC deviation
ΔNDTideviation in non-dimensional Ti
ΔNDWideviation in non-dimensional Wi
NCGnon-condensable gas
NGnon-condensable gas model
Subscripts
caliterative calculation value
iliquid and gas interface
lliquid phase
wwall
vgas phase
xlocal parameter
infinity
0pure steam

References

  1. Mualim, A.; Huda, H.; Altway, A.; Sutikno, J.P.; Handogo, P. Evaluation of multiple time carbon capture and storage network with capital-carbon trade-off. J. Clean. Prod. 2021, 291, 125710. [Google Scholar] [CrossRef]
  2. Freire, A.L.; Jose, H.J.; Moreira, R.F.P.M. Potential applications for geopolymers in carbon capture and storage. Int. J. Greenh. Gas Control. 2022, 118, 103687. [Google Scholar] [CrossRef]
  3. Intergovernmental Panel on Climate Change (IPCC). Contribution of Working Group I. Contribution to the Sixth Assessment Report, the Physical Science Basis; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  4. Chen, Y.P.; Zhu, Z.L.; Wu, J.F.; Yang, S.F.; Zhang, B.H. A novel LNG/CO2 combustion gas and steam mixture cycle with energy storage and CO2 capture. Energy 2017, 120, 128–137. [Google Scholar] [CrossRef]
  5. Zhao, W.J.; Chen, Y.P.; Wu, J.F.; Zhu, Z.L.; Wang, J. Study on CO2 capture upgrading of existing coal fired power plants with gas-H2O mixture cycle and supercritical water coal gasification. Int. J. Greenh. Gas Control. 2021, 112, 103482. [Google Scholar] [CrossRef]
  6. High, M.; Patzschke, C.F.; Zheng, L.Y.; Zeng, D.W.; Gavalda-Diaz, O.; Ding, N.; Chien, K.H.H.; Zhang, Z.L.; Wilson, G.E.; Berenov, A.V.; et al. Precursor engineering of hydrotalcite-derived redox sorbents for reversible and stable thermochemical oxygen storage. Nat. Commun. 2022, 13, 5109. [Google Scholar] [CrossRef]
  7. Yi, Q.J.; Tian, M.C.; Yan, W.J.; Qu, X.H.; Chen, X.B. Visualization study of the influence of non-condensable gas on steam condensation heat transfer. Appl. Therm. Eng. 2016, 106, 13–21. [Google Scholar] [CrossRef]
  8. Kuhn, S.Z.; Schrock, V.E.; Peterson, P.F. An investigation of condensation from steam-gas mixtures flowing downward inside a vertical tube. Nucl. Eng. Des. 1997, 177, 53–69. [Google Scholar] [CrossRef]
  9. Rose, J.W. Approximate equations for forced-convection condensation in the presence of a non-condensing gas on a flat-plate and horizontal tube. Int. J. Heat Mass Transf. 1980, 23, 539–546. [Google Scholar] [CrossRef]
  10. Huhtiniemi, I.K.; Corradini, M.L. Condensation in the presence of noncondensable gases. Nucl. Eng. Des. 1993, 141, 429–446. [Google Scholar] [CrossRef]
  11. Kang, H.C.; Kim, M.H. Characteristics of film condensation of supersaturated steam-air mixture on a flat plate. Int. J. Multiph. Flow 1999, 25, 1601–1618. [Google Scholar] [CrossRef]
  12. Aldiwany, H.K.; Rose, J.W. Free convection film condensation of steam in the presence of non-condensing gases. Int. J. Heat Mass Transf. 1973, 16, 1359–1369. [Google Scholar] [CrossRef]
  13. Morrison, J.N.A.; Deans, J. Augmentation of steam condensation heat transfer by addition of ammonia. Int. J. Heat Mass Transf. 1997, 40, 765–772. [Google Scholar] [CrossRef]
  14. Sparrow, E.M.; Yang, J.W.; Scott, C.J. Free convection in an air-water vapor boundary layer. Int. J. Heat Mass Transf. 1966, 9, 53. [Google Scholar] [CrossRef]
  15. Sparrow, E.M.; Minkowycz, W.J.; Saddy, M. Forced convection condensation in the presence of noncondensables and interfacial resistance. Int. J. Heat Mass Transf. 1967, 10, 1829–1845. [Google Scholar] [CrossRef]
  16. Minkowycz, W.J.; Sparrow, E.M. Condensation heat transfer in the presence of noncondensables, interfacial resistance, superheating, variable properties, and diffusion. Int. J. Heat Mass Transf. 1966, 9, 1125–1144. [Google Scholar] [CrossRef]
  17. Millsaps, K.; Pohlhausen, K. The Laminar-free-convective-heat transfer from the outer surface of a vertical circular cylinder. J. Aeronaut. Sci. 1958, 25, 357–360. [Google Scholar] [CrossRef]
  18. Popiel, C.O. Free convection heat transfer from vertical slender cylinders: A review. Heat Transf. Eng. 2008, 29, 521–536. [Google Scholar] [CrossRef]
  19. Angelino, M.; Boghi, A.; Gori, F. Numerical solution of three-dimensional rectangular submerged jets with the evidence of the undisturbed region of flow. Numer. Heat Transf. Part A 2016, 70, 815–830. [Google Scholar] [CrossRef]
  20. Boghi, A.; Venuta, I.D.; Gori, F. Passive scalar diffusion in the near field region of turbulent rectangular submerged free jets. Int. J. Heat Mass Transf. 2017, 112, 1017–1031. [Google Scholar] [CrossRef]
  21. Boghi, A.; Angelino, M.; Gori, F. Numerical Evidence of an Undisturbed Region of Flow in Turbulent Rectangular Submerged Free Jet. Numer. Heat Transf. 2016, 70, 14–29. [Google Scholar] [CrossRef]
  22. Hossein, S.; Nima, A.; Ali, J. Effects of Applying Brand-New Designs on the Performance of PEM Fuel Cell and Water Flooding Phenomena. Iran. J. Chem. Chem. Eng.-Int. Engl. Ed. 2022, 41, 618–635. [Google Scholar]
  23. Dehbi, A.; Janasza, F.; Bell, B. Prediction of steam condensation in the presence of noncondensable gases using a CFD-based approach. Nucl. Eng. Des. 2013, 258, 199–210. [Google Scholar] [CrossRef]
  24. Li, J.-D. CFD simulation of water vapor condensation in the presence of non-condensable gas in vertical cylindrical condensers. Int. J. Heat Mass Transf. 2013, 57, 708–721. [Google Scholar] [CrossRef] [PubMed]
  25. Choudhury, R.; Das, U.J. Viscoelastic effects on the three- Dimensional hydrodynamic flow past a vertical porous plate. Int. J. Heat Technol. Assoc. 2013, 31, 1–8. [Google Scholar]
  26. Ge, M.H.; Zhao, J.; Wang, S.X. Experimental investigation of steam condensation with high concentration CO2 on a horizontal tube. Appl. Therm. Eng. 2013, 61, 334–343. [Google Scholar] [CrossRef]
  27. Ge, M.H.; Wang, S.X.; Zhao, J.; Zhao, Y.L.; Liu, L.S. Condensation of steam with high CO2 concentration on a vertical plate. Exp. Therm. Fluid Sci. 2016, 75, 147–155. [Google Scholar] [CrossRef]
  28. Lu, J.H.; Cao, H.S.; Li, J.M. Experimental study of condensation heat transfer of steam in the presence of non-condensable gas CO2 on a horizontal tube at sub-atmospheric pressure. Exp. Therm. Fluid Sci. 2019, 105, 278–288. [Google Scholar] [CrossRef]
  29. Takami, K.M.; Mahmoudi, J.; Time, R.W. A simulated H2O/CO2 condenser design for oxy-fuel CO2 capture process. Energy Procedia 2009, 1, 1443–1450. [Google Scholar] [CrossRef]
  30. Lu, J.H.; Ren, K.X.; Tang, J.J. Numerical simulations of laminar film condensation of H2O/Air or H2O/CO2 on a vertical plate. Heat Transf. Eng. 2022, 44, 720–733. [Google Scholar] [CrossRef]
  31. Han, J.M.; Shin, H.Y.; Min, B.M. Measurement and correlation of high pressure phase behavior of carbon dioxide+water system. Ind. Eng. Chem. Eng. Data 2009, 15, 212–216. [Google Scholar] [CrossRef]
  32. Lucile, F.; Cezac, P.; Contamine, F.; Serin, J.P.; Houssin, D.; Arpentinier, P. Solubility of Carbon Dioxide in Water and Aqueous Solution Containing Sodium Hydroxide at Temperatures from (293.15 to 393.15) K and Pressure up to 5 MPa: Experimental Measurements. J. Chem. Eng. Data 2012, 57, 784–789. [Google Scholar] [CrossRef]
  33. Hou, S.X.; Maitland, G.C.; Trusler, J.P.M. Measurement and modeling of the phase behavior of the (carbon dioxide+water) mixture at temperature from 298.15 K to 448.15 K. J. Supercrit. Fluids 2013, 73, 87–96. [Google Scholar] [CrossRef]
  34. Meyer, C.W.; Harvey, A.H. Dew-point measurements for water in compressed carbon dioxide. AICHE J. 2015, 61, 2913–2925. [Google Scholar] [CrossRef]
  35. Mao, S.D.; Zhang, D.H.; Li, Y.Q.; Liu, N.Q. An improved model for calculating CO2 solubility in aqueous NaCl solutions and the application to CO2-H2O-NaCl fluid inclusions. Chem. Geol. 2013, 347, 43–58. [Google Scholar] [CrossRef]
  36. Liu, G.Q.; Ma, L.X. Chemical Physics Data; Chemical Industry Press: Beijing, China, 2002; pp. 84–85. [Google Scholar]
  37. Duan, Z.H.; Sun, R.; Zhu, C.; Chou, I.M. An improved model for the calculation of CO2 solubility in aqueous solutions containing Na+, K+, Ca2+, Mg2+, Cl, SO42−. Mar. Chem. 2006, 98, 131–139. [Google Scholar] [CrossRef]
  38. Valtz, A.; Chapoy, A.; Coquelet, C.; Paricaud, P.; Richon, D. Vapor-liquid equilibria in the carbon dioxide–water system, measurement and modelling from 278.2 to 318.2 K. Fluid Phase Equilibria 2004, 226, 333–344. [Google Scholar] [CrossRef]
  39. Fronk, B.M.; Garimella, S. Condensation of ammonia and high-temperature-glide zeotropic ammonia/water mixtures in minichannels-Part II: Heat transfer models. Int. J. Heat Mass Transf. 2016, 101, 1357–1373. [Google Scholar] [CrossRef]
  40. Wang, X.Y.; Ming, Z.Y. Measurements of gas enthalpy of dissolution. China Acad. J. Electron. Publ. House 1987, 1, 62–64. [Google Scholar]
  41. Koschel, D.; Coxam, J.Y.; Rodier, L.; Majer, V. Enthalpy and solubility data of CO2 in water and NaCl(aq) at conditions of interest for geological sequestration. Fluid Phase Equilibria 2006, 247, 107–120. [Google Scholar] [CrossRef]
  42. NIST. NIST Reference Fluid Thermodynamic and Transport Properties-REFPROP (Version 9.1); NIST Stand Ref Database 2007:23; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2013.
  43. Othmer, D.F. The condensation of steam. Ind. Eng. Chem. 1929, 21, 576–583. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of a boundary layer model of steam–CO2 mixture condensation.
Figure 1. Schematic diagram of a boundary layer model of steam–CO2 mixture condensation.
Applsci 13 05747 g001
Figure 2. Iterative for Tix and αix.
Figure 2. Iterative for Tix and αix.
Applsci 13 05747 g002
Figure 3. Comparison between the boundary layer model data and data from Ge et al. [27] study (P = 101.325 kPa, α = 34.3%, T = 367.74 K).
Figure 3. Comparison between the boundary layer model data and data from Ge et al. [27] study (P = 101.325 kPa, α = 34.3%, T = 367.74 K).
Applsci 13 05747 g003
Figure 4. Data trends of steam–CO2 condensation on the flat plate. (a) Variation in hx along the horizontal plate; (b) Variation in δ along the horizontal plate; (c) Condensation film thickness and gas–liquid interface velocity along the vertical plate.
Figure 4. Data trends of steam–CO2 condensation on the flat plate. (a) Variation in hx along the horizontal plate; (b) Variation in δ along the horizontal plate; (c) Condensation film thickness and gas–liquid interface velocity along the vertical plate.
Applsci 13 05747 g004
Figure 5. Comparison of CG and NG models for steam–CO2 mixture condensation (U = 1 m/s; L = 0.2 m). (a) α=2% average HTC; (b) α=2% q/q0; (c) α=30% average HTC; (d) α=30% q/q0; (e) α=80% average HTC.
Figure 5. Comparison of CG and NG models for steam–CO2 mixture condensation (U = 1 m/s; L = 0.2 m). (a) α=2% average HTC; (b) α=2% q/q0; (c) α=30% average HTC; (d) α=30% q/q0; (e) α=80% average HTC.
Applsci 13 05747 g005
Figure 6. Comparison of ΔHTC and CO2 solubility values at CO2 inlet fractions of 2%, 30%, and 80% (U = 1 m/s; L = 0.2 m). (a) Comparison of ΔHTC at different α; (b) Comparison of S at different α.
Figure 6. Comparison of ΔHTC and CO2 solubility values at CO2 inlet fractions of 2%, 30%, and 80% (U = 1 m/s; L = 0.2 m). (a) Comparison of ΔHTC at different α; (b) Comparison of S at different α.
Applsci 13 05747 g006
Figure 7. Comparison of Ti, αi, and δ curves at total pressures of 10 MPa. (a) Ti versus subcooling; (b) αi versus subcooling; (c) Film thickness along the horizontal plate at 7 K subcooling (U = 1 m/s; L = 0.2 m).
Figure 7. Comparison of Ti, αi, and δ curves at total pressures of 10 MPa. (a) Ti versus subcooling; (b) αi versus subcooling; (c) Film thickness along the horizontal plate at 7 K subcooling (U = 1 m/s; L = 0.2 m).
Applsci 13 05747 g007
Figure 8. Comparison of steam–CO2 and steam–air mixture condensation on horizontal plate at the same inlet CO2 and air mass percentage (U = 1 m/s; L = 0.2 m). (a) Variation in Ti with subcooling; (b) Variation in ΔNDTi with subcooling; (c) Variation in αi with subcooling; (d) Variation in ΔNDαi with subcooling; (e) Variation in h ¯ with subcooling; (f) Variation in Δ h ¯ with subcooling.
Figure 8. Comparison of steam–CO2 and steam–air mixture condensation on horizontal plate at the same inlet CO2 and air mass percentage (U = 1 m/s; L = 0.2 m). (a) Variation in Ti with subcooling; (b) Variation in ΔNDTi with subcooling; (c) Variation in αi with subcooling; (d) Variation in ΔNDαi with subcooling; (e) Variation in h ¯ with subcooling; (f) Variation in Δ h ¯ with subcooling.
Applsci 13 05747 g008
Figure 9. Comparison of condensation ΔHTC at CO2 inlet fractions of 2% and 10% on horizontal and vertical plates (U = 1 m/s; L = 0.2 m).
Figure 9. Comparison of condensation ΔHTC at CO2 inlet fractions of 2% and 10% on horizontal and vertical plates (U = 1 m/s; L = 0.2 m).
Applsci 13 05747 g009
Figure 10. Comparison of steam–CO2 mixture condensation on the horizontal and vertical surfaces. (a) Average HTC; (b) q/q0; (c) δ; (d) Ratio of temperature differences in gas and liquid-phases.
Figure 10. Comparison of steam–CO2 mixture condensation on the horizontal and vertical surfaces. (a) Average HTC; (b) q/q0; (c) δ; (d) Ratio of temperature differences in gas and liquid-phases.
Applsci 13 05747 g010
Table 1. Governing equations and boundary conditions of the NG model.
Table 1. Governing equations and boundary conditions of the NG model.
PhaseLiquidGas
PlateHorizontalVerticalHorizontalVertical
Continuity u l x + v l y = 0 (1a) u v x + v v y = 0 (1b)
Momentum conservation 2 u l y 2 = 0   (2a) g + ν l 2 u l y 2 = 0   (2b) u v u v x + v v u v y = ν v 2 u v y 2   (2c) u v u v x + v v u v y = g ( 1 ρ ρ v ) + ν v 2 u v y 2   (2d)
Energy conservation or mass conservation 2 T l y 2 = 0 (3a) u v α v x + v v α v y = D 2 α v y 2 (3b)
Note: u and v are the velocities in the x and y directions, respectively. g is gravity. ρ is the density. ν is the kinematic viscosity. D is diffusion coefficients. The subscript l represents the liquid phase, and v represents the gas phase.
Table 2. Derivation and solution of the NG model.
Table 2. Derivation and solution of the NG model.
PlatePhaseDerivation and Solution
HorizontalLiquid Introduce   intermediate   variables , (4a) η = y U υ l x ; ψ l = U ν l x f ( η ) ; θ = T l T w T i x T w u l = ψ l y = U f , v l = ψ l x = 1 2 U ν l x ( η f f ) Ordinary   differential   equation   equations , (4b) f ( η ) = 0 ; θ ( η ) = 0 Solution , (4c) f = 1 2 f ( 0 ) η 2 ;   θ = η η i ; η i = δ U ν l x
Gas Introduce   intermediate   variables , (5a) ξ = ( y δ ) U ν v x , φ v ( ξ ) = α v α α i x α , ψ v = U ν v x F ( ξ ) u v = ψ v y = U F , v v = ψ v x = 1 2 U ν v x ( ξ F F ) + U d δ d x F Ordinary   differential   equation   equations , (5b) F + 1 2 F F = 0 ; ϕ + 1 2 S c F ϕ = 0 Solution , (5c) F ( 0 ) = 0 ; ρ l μ l ρ v μ v C p ( T i x T w ) h f g P r = F ( 0 ) 3 2 F ( 0 ) ; η i = 2 F ( 0 ) F ( 0 ) F ( 0 ) = 2 u i U x U ν ; α i x = α 1 + 1 / 2 S c F ( 0 ) φ v ( 0 )
VerticalLiquid Solve   directly , (6) u i x = [ C p ( T i x T w ) h f g P r ] 1 / 2 ( g x ) 1 / 2 ; m = q h f g
Gas Introduce   intermediate   variables , (7a) ξ = c ( y δ ) x 1 / 4 ; c = [ g ( M CO 2 M H 2 O ) / 4 ν v 2 M CO 2 ( M CO 2 M H 2 O ) α ] 1 / 4 ψ v = 4 ν v c x 3 / 4 F ( ξ ) ; φ v = α v α Ordinary   differential   equation   equations , (7b) F + 3 F F 2 ( F ) 2 + ϕ = 0 ϕ + 3 S c F ϕ = 0 Solution , (7c) F ( 0 ) = g 4 ν v c 2 [ C p ( T i x T w ) h f g P r ] 1 / 2 ; F ( 0 ) = g 4 3 2 ν v c ρ l μ l ρ v μ v [ C p ( T i x T w ) h f g P r ] 3 / 4 (7d) φ v ( 0 ) = 3 S c F ( 0 ) [ φ v ( 0 ) + α ]
Note: f (for liquid) and F (for gas) are the functions of η and ξ, and η and ξ are the dimensionless coordinates. θ is the dimensionless temperature, φ is the concentration function, ψ is the stream function, CP is the specific heat, hfg is the latent heat, Pr is the Prandtl number, Sc is the Schmidt number, q is the heat flux, m is the condensation mass flux, M is the molar mass, and μ is the dynamic viscosity.
Table 3. The solution enthalpy of CO2.
Table 3. The solution enthalpy of CO2.
Pressure
(MPa)
Temperature
(K)
Solution Enthalpy
(kJ/mol)
Solution Enthalpy
(kJ/kg), hso
hso/hfghso∙S/hfg
0.101325298.1519.75 [40]448.8640.18380.0002698
2.06323.114.8 [41]336.3640.14120.002223
5.10323.113.1 [41]297.7270.12500.004401
10.53323.17.5 [41]170.4550.071550.003569
5.05373.16.3 [41]143.1820.063450.001235
10.08373.14.6 [41]104.5450.046330.001561
Table 4. Data comparison with Sparrow.
Table 4. Data comparison with Sparrow.
Numerical Results F ( 0 ) F ( 0 ) ρ l μ l ρ v μ v C p ( T i x T w ) h f g P r 1 η i α α i
Sparrow [15]0.10.368670.0368271.35770.90510
0.20.406120.0992441.00760.82332
0.40.483250.257330.77720.69002
0.60.563000.437990.68490.58654
10.728870.828250.60370.43803
21.169431.84950.54070.24075
This work0.10.368650.0368281.35770.90509
0.20.406090.0992471.00760.82330
0.40.483210.257340.77720.68998
0.60.56290.438000.68490.58648
10.728800.828290.60380.43794
21.169291.84960.54080.24064
Table 5. Data comparison with Othmer (TTw = 5 K).
Table 5. Data comparison with Othmer (TTw = 5 K).
T (K)αP (kPa)q/q0 [43]q/q0 of This WorkDeviation
383.150.0171450.44560.385713.44%
383.150.0311146.30.32550.30576.08%
383.150.0457147.750.26420.24796.17%
373.150.0546105.130.23770.202814.68%
373.150.0226102.910.36230.33756.85%
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

Jiang, B.; Jiang, Y.; Gu, H.; Chen, Y.; Wu, J. Numerical Study of Steam–CO2 Mixture Condensation over a Flat Plate Based on the Solubility of CO2. Appl. Sci. 2023, 13, 5747. https://doi.org/10.3390/app13095747

AMA Style

Jiang B, Jiang Y, Gu H, Chen Y, Wu J. Numerical Study of Steam–CO2 Mixture Condensation over a Flat Plate Based on the Solubility of CO2. Applied Sciences. 2023; 13(9):5747. https://doi.org/10.3390/app13095747

Chicago/Turabian Style

Jiang, Bingran, Yi’ao Jiang, Huaduo Gu, Yaping Chen, and Jiafeng Wu. 2023. "Numerical Study of Steam–CO2 Mixture Condensation over a Flat Plate Based on the Solubility of CO2" Applied Sciences 13, no. 9: 5747. https://doi.org/10.3390/app13095747

APA Style

Jiang, B., Jiang, Y., Gu, H., Chen, Y., & Wu, J. (2023). Numerical Study of Steam–CO2 Mixture Condensation over a Flat Plate Based on the Solubility of CO2. Applied Sciences, 13(9), 5747. https://doi.org/10.3390/app13095747

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