Dynamic modeling of interfacial structures via interfacial area transport equation

Full text of publication follows:In the current thermal-hydraulic system analysis codes using the two-fluid model, the empirical correlations that are based on the two-phase flow regimes and regime transition criteria are being employed as closure relations for the interfacial transfer terms. Due to its inherent shortcomings, however, such static correlations are inaccurate and present serious problems in the numerical analysis. In view of this, a new dynamic approach employing the interfacial area transport equation has been studied. The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Therefore, the interfacial area transport equation can make a leapfrog improvement in the current capability of the two-fluid model from both scientific and practical point of view. Accounting for the substantial differences in the transport phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right hand-side of the transport equations have been established by mechanistically modeling the creation and destruction of bubbles due to major bubble interaction mechanisms. The coalescence mechanisms include the random collision driven by turbulence, and the entrainment of trailing bubbles in the wake region of the preceding bubble. The disintegration mechanisms include the break-up by turbulence impact, shearing-off at the rim of large cap bubbles and the break-up of large cap bubbles due to surface instability. In the present paper, the interfacial area transport equations currently available are reviewed to address the feasibility and reliability of the model. Results from extensive benchmark experiments for the model evaluation are also present. These include the data from adiabatic upward air-water two-phase flow in round tubes of various sizes, from a rectangular duct, and from adiabatic co-current downward air-water two-phase flow in round pipes of two different sizes. Furthermore, some guidelines for the future study on interfacial area transport equation are discussed. (authors)