ORIGINAL RESEARCH article

Front. Energy Res., 13 April 2022

Sec. Nuclear Energy

Volume 10 - 2022 | https://doi.org/10.3389/fenrg.2022.890795

Comparison of Intergroup Mass Transfer Coefficient Correlations in Two-Group IATE for Subcooled Boiling Flow

  • 1. Key Laboratory of Low-Grade Energy Utilization Technologies and Systems (Chongqing University), Ministry of Education, Chongqing, China

  • 2. Postdoctoral Station of Power Engineering and Engineering Thermophysics, Chongqing University, Chongqing, China

  • 3. Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, United States

Abstract

The interfacial area concentration (IAC) is of vital importance in determining the interfacial transfer terms of mass, momentum, and energy between phases in the two-fluid model. The two-group (2-G) interfacial area transport equation (IATE) dynamically models the IAC evolution of large and small bubbles, and the intergroup transfer bookkeeps the mass transferring between two groups. Different intergroup mass transfer coefficient correlations are employed in the 2-G IATE model, which include one old correlation and three new correlations previously developed by the authors. The two-group IATE is benchmarked with subcooled boiling experimental dataset. The group-1 interfacial area concentration results show that the three modified correlations improve the physically incorrect prediction by the old correlation. With the modified correlations, the 2-G IATE is foundationally capable of predicting the group-1 interfacial area concentration in subcooled boiling flow.

Introduction

The two-fluid model is now widely used to simulate two-phase flow systems like in RELAP5 () and (). The model is formulated with two sets of conservation equations to separately consider the vapor phase and the liquid phase () and the interface exists between the two phases. The interfacial area concentration, ai, plays a key role in the interfacial drag closure and the interfacial heat transfer closures in the two-fluid model, as it determines the interfacial transfer terms of mass, momentum, and energy between phases () (). Figure 1 shows the diagram of the interface between the liquid and the vapor phase. The interfacial transfer terms can be written in the following form,

FIGURE 1

Many researchers () developed empirical interfacial area concentration correlations, including , , , and . The correlations can be divided into two kinds, namely geometry-based semi-empirical correlations and empirical correlations (). Different correlations are used in bubbly, slug, and the following regimes because they are dependent on the flow regime. As a result, numerical instabilities can occur at the transitions between the flow regimes. However, the interfacial area transport equation (IATE) is developed to dynamically model the interfacial area concentration across the flow regimes. The successful closure of interfacial area concentration in the two-fluid model relies on the interfacial area transport equation.

developed the one-group IATE based on the fluid particle number density transport equations analogous to Boltzmann’s transport equation. The one-group IATE disregards the type of bubbles and treats them equally as spherical bubbles. provided the formulation of one-group IATE for subcooled boiling flow with the bubble layer thickness model. proposed a complete framework of the one-dimensional one-group (1-G) IATE model for phase change flow, then further validate the 1-G IATE using extensive boiling flow datasets and improved the condensation closure in the model. At low void fractions, the one-group two-fluid model is sufficient for describing cases with uniform bubble size (). However, for high void fraction cases, large bubbles in cap shape or slug shape come into existence as the consequence of bubble coalescence, bubble expansion, or heating from the wall. Under the high-void-fraction circumstance, the one-group IATE does not accurately predict ai due to the spherical-shape and uniform-size assumption being applied to large bubbles ().

As the void fraction increases, the differences in shape, size, and number density between larger bubbles and small bubbles necessitate the separate consideration of bubbles in the two-fluid model (). In order to account for the different mechanisms for heat and mass transfer of bubbles with different shapes and sizes, the two-group (2-G) IATE model is introduced in the work of . Spherical and distorted bubbles are categorized as group-1 bubbles while cap and slug bubbles are categorized as group-2 bubbles in the two-group model (). The separation of bubble groups thereby requires an additional set of conservation equations for the gas phase. proposed the one-dimensional two-group IATE for subcooled boiling flow considering the uniform phase distribution pattern across the channel. Group-1 bubbles are averaged within a bubble layer while the group-2 bubbles span all the way across the channel.

Intergroup transfer happens between group-1 bubbles and group-2 bubbles: small bubbles can grow into large ones contributed from the coalescence, expansion, and nucleation, while the large bubbles can shrink into small ones as a result of the disintegration, condensation, and pressurization. Figure 2 shows the transfer between small and big bubbles. Intergroup transfer of the interfacial area concentration between groups is important for two-group IATE as it tracks and partitions of the gas phase into the bubble groups. This bookkeeping of bubble groups is facilitated by the intergroup mass transfer coefficient. The intergroup mass transfer correlation was first developed by . provided a thorough review of the intergroup mass transfer coefficient and proposed three modified correlations for the term. Two analytical intergroup mass transfer coefficient correlations are also developed based on the bubble size distribution functions, including the Nukiyama-Tanasawa distribution and the Rosin-Rammler distribution . Also, one empirical intergroup mass transfer coefficient correlation is obtained from the experimental data.

FIGURE 2

This paper focuses on the benchmark of the 2-G IATE with subcooled boiling experimental data using different intergroup mass transfer coefficient correlations. The modeling of two-group IATE and the various correlations of intergroup mass transfer coefficient are presented in section 2. The prediction results of interfacial area concentration and the sink/source terms are presented in section 3, accompanied by the discussion about the performance of the old and three modified correlations.

Modeling

proposed the one-dimensional two-group IATE for subcooled boiling flow with the consideration of bubble layer thickness model; however, only intergroup mass transfer from group-1 bubbles to group-2 bubbles was modeled. developed the model to include the intergroup mass transfer from group-2 to group-1. Therefore, incorporating the intergroup mass transfer from both group-1 and group-2 bubble sides, the change in group ai in the one-dimensional, steady-state two-group IATE model for subcooled boiling flow is,andwhere subscripts g, 1, 2, and B denote vapor, group-1, group-2, and averaging over the bubble layer, respectively. The brackets <> and <<>> represent cross-sectional area average and void-weighted area average. The bracket <>B represent the bubble-layer average, which is calculated by the improved flat bubble layer thickness model in . The terms ai, z, v, and ϕ are interfacial area concentration, axial coordinate, velocity, and interfacial area concentration source/sink rate, respectively. There are nine source/sink terms in Eqs 2, 3: the wall nucleation term, the two expansion terms, the two bubble interaction terms, the two convection terms, the condensation term, and the intergroup mass transfer term. The wall nucleation term, <ϕWN,1>, and the expansion terms, <ϕEXP>, are modeled and described in . The bubble interaction terms, <ϕj>, and convection term, <ϕConv>, are described in . The condensation term for subcooled boiling flow, <ϕCO,1>, is modified by . The intergroup transfer term, <ϕΔṁ12>, is modeled by and is written as,where,where,where the terms α, ρ, η, χ, Dsm, Dc are void fraction, density, volume change rate per unit mixture volume, intergroup mass transfer coefficient, Sauter mean diameter, and critical bubble diameter at the group boundary, respectively. The subscripts V, Γ, DP, 12, and 21 denote volume change, interphase mass transfer, pressure change, group-1 to group-2, and group-2 to group-1, respectively.

In Eqs 6, 8, the interphase mass transfer rate per unit volume for group-n vapor, Γg,n, is a key term which cannot be measured directly and is hard to estimate in subcooled boiling flow with the theoretical models. However, this term can be determined based on the one-dimensional two-group mass balance equations (),

The intergroup mass transfer term, Δ12, is difficult to model in subcooled boiling flow, so it is assumed to be small relative to Γg,1 in Eq. 9. The relative magnitude of Δ12 to Γg,2 in Eq. 10 is inconsequential to the calculation of <ϕΔṁ12> due to the suppression using the min function in Eq. 4. Therefore, Δ12 is assumed to be negligible in Eqs 9, 10. Even though in some instances (e.g., when χ1 is nonzero) Δ12 will not be negligible relative to Γg,1, the analysis uses the same Γg,1 to isolate the effect of χ1 on <ϕΔṁ12>. More rigorous analysis to the two-group IATE can be applied once Γg,n is modeled correctly.

In Eqs. 6, 8, the intergroup mass transfer coefficient for group-n vapor, χn, is another important parameter in the multigroup formulation and needs to be modeled accurately. The group-1 intergroup mass transfer coefficient, χ1, and the group-2 intergroup mass transfer coefficient, χ2, can be related analytically () as χ2=(χ1·nb,1)/nb,2 where nb,n is the number density of group-n bubbles. The group-1 intergroup mass transfer coefficient was first developed by , then proposed three new correlations: an analytical correlation with Nukiyama-Tanasawa distribution, an analytical correlation with Rosin-Rammler distribution, and an empirical correlation. The correlations for the group-1 intergroup mass transfer coefficient are summarized in Table 1. In the equations, Dc* is the dimensionless critical bubble diameter which is calculated Dc* = 1.5Dc/Dsm,1. G (z,x) and G(z) are the upper incomplete Gamma function and the complete Gamma function. In addition, a constant value of 0, 1, or 0.01 has been assigned to χ1 in some works ().

TABLE 1

ModelCorrelation
Analytical correlation with Nukiyama-Tanasawa distribution ()
Analytical correlation with Rosin-Rammler distribution ()
Empirical Correlation ()
where, , , ,

Correlations for the group-1 intergroup mass transfer coefficient.

The two analytical correlations are based on the distributions proposed by , and by . Both are single-peaked, right-skewed distributions; but the height and the location of the peak are different with the same characteristic parameter, i.e., the bubble nondimensional diameter. The factors in the power-law and exponential function of the nondimensional diameter are different in the two analytical correlations. Consequentially, in the range of Dsm,1/Dc∈[0.2,0.8], the intergroup mass transfer coefficient calculated by the analytical correlation with Nukiyama-Tanasawa distribution is larger than the correlation with Rosin-Rammler distribution. The empirical correlation is developed using an annulus condensing dataset (), an annulus flashing dataset (), an annulus boiling dataset (), and a pipe flashing dataset (). The empirical correlation predicts larger intergroup mass transfer coefficient values in the range of Dsm,1/Dc∈[0.2,0.5]; thus it better captures the data-points trend of the four datasets compared with the analytical correlations.

Results and Discussion

The subcooled boiling database by is used to validate the two-group IATE model with four different group-1 intergroup mass transfer coefficient correlations: correlation , analytical-NT correlation, analytical-RR correlation, and empirical correlation (). The experiment was conducted in an internally heated annulus channel with inner and outer diameters of 19.1 and 38.1 mm, respectively. The test section has a heated length of 2.85 m and an unheated length of 1.63 m. Four-sensor conductivity probes at three ports are used to measure interfacial area concentration and other two-phase parameters for group-1 and group-2 bubbles in the heated section. The conditions by considered for the benchmark of two-group IATE model are listed in Table 2.

TABLE 2

Run#Inlet pressure (kPa)Heat flux (kW/m2)Inlet liquid velocity (m/s)Inlet subcooling (°C)Max void fraction
1232.6680.955.350.61
2273.0901.255.950.52
3248.7560.686.470.55
4489.9901.325.120.31
5900.8901.334.040.24

Summary of conditions considered in the benchmark.from the database of .

Figures 36 present the calculation results of Run 1 in the database of with two-group IATE model using the correlation , the analytical-NT correlation, the analytical-RR correlation, and the empirical correlation (). The subfigure A) in each figure is the calculated group-1 intergroup mass transfer coefficients. The χ1 calculated by the correlation is observed to exceed the physical limit, 2, stated by , in the lower region of the heated section where void fraction is small; thus, it is manually set to be two in the calculation. The cause of the nonphysical value is the power-law relation between χ1 and α1 in the correlation , shown in the equation in Table 1. The χ1 calculated by the correlation is also observed to decrease when the void fraction increases along the heated length. On the contrary, χ1 calculated by analytical-NT correlation, analytical-RR correlation, and empirical correlation () increases along the heated length in the lower test section. This trend is more physically appropriate: the group-1 intergroup mass transfer coefficient should increase when more group-1 bubbles are present as χ1 evaluates the amount of group-1 bubbles transferring to group-2 bubbles. Compared to the χ1 calculated by the correlation , χ1 calculated by the other three correlations also has a smaller range, the maximum of which is ∼0.18 for the three cases.

FIGURE 3

.

FIGURE 4

).

FIGURE 5

).

FIGURE 6

) (fourth row).

In the figures, the subfigure B) shows the source terms of group-1 interfacial area concentration: the intergroup mass transfer, wall nucleation, condensation, expansion, convection, and bubble interaction. The sum of the terms is shown with a solid blue line in the subfigures, and the intergroup mass transfer term is shown as a dashed black line. The subfigure C) is the comparison of simulated interfacial area concentration with experimental data for group-1 bubbles presented as red dots and red line. In the source term subfigures with the correlation , the sum of the terms is negative or even infinite, dominated by a strong intergroup mass transfer term caused by the large χ1 mentioned above. The intergroup term is so large on account of the unphysically large χ1 term that any ai,1 generated by other means is immediately transferred to group-2, resulting in unphysical growth in group-2 ai while keeping ai,1 near-zero. This highlights the necessity of a new group-1 intergroup mass transfer coefficient correlation. With the new analytical-NT correlation, and analytical-RR correlation, the sums of terms calculated by 2-G IATE are positive values in the whole heated section. The calculated ai,1 is positive and increases along the heated length. The results show an improvement of the new analytical and empirical group-1 intergroup mass transfer coefficient correlations compared to the correlation . In addition, compared with the χ1 scale by analytical correlations, the constant value of 0 or 1 should not be assumed to χ1, which give under- and over-prediction of intergroup transfer term. The constant value of 0.01 should not be assigned to the whole heated section. It overpredicts in the low void fraction section and affects the section later. Thus, the correlations which can dynamically predict χ1 are needed.

Figure 7 shows the group-1 interfacial area concentration benchmark of 2-G IATE with four group-1 intergroup mass transfer coefficient correlations, i.e., the correlation , the analytical-NT correlation, the analytical-RR correlation, and the empirical correlation (). The results calculated by the correlation are not physically correct: the group-1 values of ai are zero as explained above, while the group-2 part is the contributor to the ai values here. However, when using the analytical-NT correlation, the analytical-RR correlation, and the empirical correlation, the 2-G IATE is foundationally capable to handle the transition to group-2 bubbles in subcooled boiling flow. A lot more works are still required to improve the two-group IATE in subcooled boiling flow, such as an accurate two-group interphase mass transfer model for Γg,n.

FIGURE 7

Conclusion

At high void fractions, the two-group (2-G) IATE model is developed to account for the different mechanisms for heat and mass transfer of bubbles with different shapes and sizes. Spherical and distorted bubbles are categorized as group-1 bubbles while cap and slug bubbles are categorized as group-2 bubbles. Intergroup transfer of the IAC between groups tracks and partitions the gas phase into the bubble groups, and the accurate modeling of the intergroup mass transfer coefficient is important for the two-group IATE.

In this paper, the 2-G IATE is benchmarked with subcooled boiling experimental data from . Different intergroup mass transfer coefficient correlations are employed in the 2-G IATE model, which include the correlation , the analytical correlation with Nukiyama-Tanasawa distribution (), the analytical correlation with Rosin-Rammler distribution (), and the empirical correlation by . The group-1 interfacial area concentration results show that the three modified correlations improve the physically incorrect prediction by the correlation . With the modified correlations, the 2-G IATE is foundationally capable of predicting the group-1 interfacial area concentration in subcooled boiling flow. More rigorous analysis to the two-group IATE can be applied once two-group Γg is modeled correctly. The IATE model will be coupled with system codes in the next phase to present the numerical improvement across flow regimes.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

LoZ: conceptualization, investigation, methodology, validation, formal analysis, visualization, writing—original draft, project administration; JB: writing—review and editing; CB: methodology, writing—review and editing, supervision; LuZ: writing—review and editing.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Nomenclature

  • interfacial area concentration [m−1]

  • diameter [m]

  • gap width [m]

  • intergroup mass transfer rate per unit mixture volume [kg/m3/s]

  • bubble number density [m−3]

  • axial position [m]

Greek

  • void fraction [-]

  • vapor generation rate to the k-phase per unit volume [kg/m3/s]

  • difference

  • volume source rate [s−1]

  • density [kg/m3]

  • Interfacial area concentration source or sink term [1/m/s]

  • intergroup mass transfer coefficient

Superscripts

  • nondimensional

  • boundary value between cap bubble and slug bubble

Subscripts

  • group-1

  • 12

    from group-1 to group-2

  • group-2

  • 21

    from group2 to group-1

  • critical value at the boundary of group-1 and group-2 bubbles

  • condensation

  • convection

  • pressure change

  • expansion

  • g

    vapor phase

  • interaction

  • Sauter mean

  • volume change

  • wall nucleation

Brackets

  • area-averaged

  • bubble layer averaged

  • void-weighted area-averaged

References

Summary

Keywords

two-phase flow, boiling flow, IATE, two-group, intergroup transfer

Citation

Zhu L, Bottini JL, Brooks CS and Zhang L (2022) Comparison of Intergroup Mass Transfer Coefficient Correlations in Two-Group IATE for Subcooled Boiling Flow. Front. Energy Res. 10:890795. doi: 10.3389/fenrg.2022.890795

Received

06 March 2022

Accepted

28 March 2022

Published

13 April 2022

Volume

10 - 2022

Edited by

Mingjun Wang, Xi’an Jiaotong University, China

Reviewed by

Zhaoming Meng, Harbin Engineering University, China

Yacine Addad, Khalifa University, United Arab Emirates

Updates

Copyright

*Correspondence: Longxiang Zhu,

This article was submitted to Nuclear Energy, a section of the journal Frontiers in Energy Research

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics