Churn turbulent flow is a two-phase gas/liquid flow regime characterized by a highly-agitated flow where gas bubbles are sufficient in numbers to both interact with each other and, while interacting, coalesce to form larger distorted bubbles with unique shapes and behaviors in the system. This flow regime is created when there is a large gas fraction in a system with a high gas and low liquid velocity. It is an important flow regime to understand and model because of its predictive value in nuclear reactor vessel boiling flow.

Occurrence

A flow in which the number of bubbles is low is called ideally-separated bubble flow. The bubbles don’t interact with each other. As the number of bubbles increase they start colliding each other. A situation then arises where they tend to coalesce to form cap bubbles, and the new flow pattern formed is called churn turbulent flow. The bubbles occurring in such a flow can be classified in small, large, and distorted bubbles. The small bubbles are generally spherical or elliptical and are encountered in a major concentration in the wake of large and distorted bubbles and close to the walls. Large, ellipsoidal or cap bubbles can be found in the core region of the flow as well as the distorted bubbles with a highly deformed interface.

Churn turbulent flow is commonly encountered in industrial applications. A typical example is boiling flow in nuclear reactors.

Numerical simulation of bubble column flows in churn turbulent regime

Numerical simulations of cylindrical bubble columns operating in the churn-turbulent regime have been carried out using an Euler–Euler approach incorporated with the RNG k–ε model for liquid turbulence. Several approaches have been carried out, including single-sized bubble modeling, double-sized bubble modeling, and the multiple sizes group modeling (MUSIG).

Breakup mass conserved formulations and coalescence rates mass conserved formulation was used in the computation of bubble size distributions. For single size modelling the Schiller–Naumann drag force was used, and for the modelling of MUSIG the Ishii–Zuber drag force was used. An empirical drag formulation was used for the double size bubble model. The simulation results of time-averaged axial velocity and gas holdup obtained with the three models were compared with reported experimental data in the resulting literature. After the comparison of all the three results it gets very clear that only MUSIG models with some lift force can replicate the measured radial distribution of gas holdup in the fully developed flow regime. The inhomogeneous MUSIG model gives a little better result than other models in the prediction of axial liquid velocity. For all the simulations the RNG k–ε model was used, and the results showed that this version of k–ε model did yield comparatively high rate of turbulence dissipation and high bubble breakup and, hence, a rational bubble size distribution formed. Here the ad hoc manipulation of the breakup rates was ignored. Mutual effects of drag force, mean bubble sizes, and turbulence characteristics profound from the simulation results. A decrease in the relative velocity between two phases is encounters due to an increase in the drag force, and this could result in decrease in k and ε. Low breakup rates results a large Sauter diameter which was directly connected to the dissipation rates of turbulence. Drag force is directly influenced by the change of Sauter diameter.

References

• Montoya, G.; Liao, Y.; Lucas, D.; Krepper, E. "Analysis and Applications of a Two-Fluid Multi-Field Hydrodynamic Model for Churn-Turbulent Flows", 21st International Conference on Nuclear Engineering – ICONE 21. China (2013)
• Montoya, G.; Baglietto, E.; Lucas, D.; Krepper, E. "A Generalized Multi-Field Two-Fluid Approach for Treatment of Multi-Scale Interfacial Structures in High Void-Fraction Regimes", MIT Energy Night 2013. Cambridge, Massachusetts, USA (2013)

•Montoya, G.; Lucas, D.; Krepper, E.; Hänsch, S.; Baglietto, E. Analysis and Applications of a Generalized Multi-Field Two-Fluid Approach for Treatment of Multi-Scale Interfacial Structures in High Void-Fraction Regimes 2014 International Congress on Advances in Nuclear Power Plants – ICAPP 2014. USA (2014) •Montoya, G.; Baglietto, E.; Lucas, D.; Krepper, E.; Hoehne, T. Comparative Analysis of High Void Fraction Regimes using an Averaging Euler-Euler Multi-Fluid Approach and a Generalized Two-Phase Flow (GENTOP) Concept 22nd International Conference on Nuclear Engineering – ICONE 22. Czech Republic (2014)

• Montoya, G.; Baglietto, E.; Lucas, D.; Krepper, E.

Development and Analysis of a CMFD Generalized Multi-Field Model for Treatment of Different Interfacial Scales in Churn-Turbulent and Transitional Flows CFD4NRS-5 – Application of CFD/CMFD Codes to Nuclear Reactor Safety Design and their Experimental Validation. Switzerland (2014)

• https://www.hzdr.de/db/!Publications?pSelTitle=18077&pSelMenu=-1&pNid=3016
• Shuiqing Zhan, Mao Li, Jiemin Zhou, Jianhong YangYiwen Zhou Applied Thermal Engineering 2014, 73, 803–816 [CrossRef]
• T. T. DeviB. Kumar Thermophysics and Aeromechanics 2014, 21, 365–382 [CrossRef]
• R.M.A. MasoodA. Delgado Chemical Engineering Science 2014, 108, 154–168 [CrossRef]

^[1]

• ×M. Pourtousi, J.N. SahuP. Ganesan Chemical Engineering and Processing: Process Intensification 2014, 75, 38–47 [CrossRef]