Numerical Investigation on Turbulence and Bubbles Distribution in Bubbly Flow Under Normal Gravity and Microgravity Conditions

Two-phase flows of gas and liquid are increasingly paid much attention to space application due to excellent properties of heat and mass transfer, so it is very meaningful to develop studies on them in microgravity. In this paper, gas-phase distribution and turbulence characteristics of bubbly flow in normal gravity and microgravity were investigated in detail by using Euler–Lagrange two-way model. The liquid-phase velocity field was solved by using direct numerical simulations (DNS) in Euler frame of reference, and the bubble motion was tracked by using Newtonian motion equations that took into account interphase interaction forces including drag force, shear lift force, wall lift force, virtual mass force and inertia force, etc. in Lagrange frame of reference. The coupling between gas–liquid phases was made with regarding interphase forces as a momentum source term in the momentum equation of the liquid phase. Under the normal gravity condition, a great number of bubbles accumulate near the walls under the influence of the shear lift force, and addition of bubbles reduces turbulence of the liquid phase. Different from the normal gravity condition, in microgravity, an overwhelming majority of bubbles migrate towards the centre of the channel driven by the pressure gradient force, and bubbles have little effect on the turbulence of the liquid phase.     

Two-phase bubbly flows in microgravity: Some open questions

Several studies on gas-liquid pipe flows in micro gravity have been performed. They were motivated by the technical problems arising in the design of the thermohydraulic loops for the space applications. Most of the studies were focused on the determination of the flow pattern, wall shear stress, heat transfer and phase fraction and provided many empirical correlations. Unfortunately some basic mechanism are not yet well understood in micro gravity. For example the transition from bubbly to slug flow is well predicted by a critical value of the void fraction depending on an Ohnesorge number, but the criteria of transition cannot take into account the pipe length and the bubble size at the pipe inlet. To improve this criteria, a physical model of bubble coalescence in turbulent flow is used to predict the bubble size evolution along the pipe in micro gravity, but it is still limited to bubble smaller than the pipe diameter and should be extended to larger bubbles to predict the transition to slug flow. Another example concerns the radial distribution of the bubbles in pipe flow, which control the wall heat and momentum transfers. This distribution is very sensitive to gravity. On earth it is mainly controlled by the action of the lift force due to the bubble drift velocity. In micro gravity in absence of bubble drift, the bubbles are dispersed by the turbulence of the liquid and the classical model fails in the prediction of the bubble distribution. The first results of experiments and numerical simulations on isolated bubbles in normal and micro gravity conditions are presented. They should allow in the future improving the modelling of the turbulent bubbly flow in micro gravity but also on earth.     

Specificity of laminar-turbulent transition in upward monodispersed microbubbly flow

Results of experimental investigation of the wall shear stress in the upward monodispersed bubbly flow in a vertical tube are presented. The bubble generator based on the flow focusing technique has been developed for monodispersed submillimeter bubbles production. The results of investigation prove that submillimeter bubbles significantly increase the flow mass transfer with the wall. Some peculiarities of the inherent liquid turbulence interaction with pseudoturbulence induced by submillimeter bubbles in transitional flow regime have been detected.     

Effect of the size of air bubbles on enhancement of heat transfer in an impinging liquid jet

The numerical modeling of heat transfer in a bubbly impinging jet is carried out. The axisymmetric system of RANS equations that take into account the two-phase nature of the flow is resolved based on the Euler approach. The turbulence of the liquid phase is described by the Reynolds stress transport model with taking into account the effect of bubbles on modification of the turbulence. The effect of the gas volumetric flow rate ratio and the bubble size on the flow structure and the heat transfer in a gas–liquid impact stream is studied. It is shown that the addition of the gas phase in a turbulent fluid causes an increase up to 1.5-fold in heat transfer. The comparison of the simulation results with experimental data showed that the developed model enables the simulation of turbulent bubbly impinging jet with heat transfer with the pipe wall in a wide range of gas fraction.     

Two-phase Flow in Fuel Cells in Short-term Microgravity Condition

A visual observation of liquid–gas two-phase flow in anode channels of a direct methanol proton exchange membrane fuel cells in microgravity has been carried out in a drop tower. The anode flow bed consisted of 2 manifolds and 11 parallel straight channels. The length, width and depth of single channel with rectangular cross section was 48.0 mm, 2.5 mm and 2.0 mm, respectively. The experimental results indicated that the size of bubbles in microgravity condition is bigger than that in normal gravity. The longer the time, the bigger the bubbles. The velocity of bubbles rising is slower than that in normal gravity because buoyancy lift is very weak in microgravity. The flow pattern in anode channels could change from bubbly flow in normal gravity to slug flow in microgravity. The gas slugs blocked supply of reactants from channels to anode catalyst layer through gas diffusion layer. When the weakened mass transfer causes concentration polarization, the output performance of fuel cells declines.     

Numerical Investigation of Phase Distribution and Liquid Turbulence Modulation in Dilute Particle-Laden Flow

Studies on the motion of particles in turbulence and interactions between particles and turbulence are extremely significant, which can help us to improve the efficiency of industrial processes. In this article, we investigated the particle distribution and particle-turbulence interaction in a solid-liquid channel flow with the Euler-Lagrange two-way model. The liquid phase was solved using direct numerical simulation (DNS), and the particle motion was tracked by Newtonian equations of motion considering effects of drag force, pressure gradient force, and gravity. Two-way coupling was used to explain the effect of particles on the turbulence structure. The results show that the local void fraction of particles indicates the wall-peaked profile, particles scatter uniformly in the spanwise direction, and the injection of particles suppresses the turbulence activities in the near wall region. Suppression of the liquid turbulence is mainly caused by vortexes decay of different sizes.     

The emission of sound by statistically homogeneous bubble layers

This paper is concerned with the flow of a bubbly fluid along a wavy wall, which is one Fourier component of a linearized hydrofoil. The bubbles are dispersed, not throughout the whole of the liquid, but only over a certain distance from the wall, as occurs in practice with cavitation bubbles. Outside the bubbly regime there is pure liquid.The interface between the bubbly fluid and pure liquid fluctuates for various reasons. One of these is the relative motion between bubbles and liquid. This is considered here in detail. A calculation is made of the sound emitted by the bubbly layer into pure liquid as a result of this stochastic motion of the interface.     

Two-phase bubbly flows in microgravity: Some open questions

Several studies on gas-liquid pipe flows in micro gravity have been performed. They were motivated by the technical problems arising in the design of the thermohydraulic loops for the space applications. Most of the studies were focused on the determination of the flow pattern, wall shear stress, heat transfer and phase fraction and provided many empirical correlations. Unfortunately some basic mechanism are not yet well understood in micro gravity. For example the transition from bubbly to slug flow is well predicted by a critical value of the void fraction depending on an Ohnesorge number, but the criteria of transition cannot take into account the pipe length and the bubble size at the pipe inlet. To improve this criteria, a physical model of bubble coalescence in turbulent flow is used to predict the bubble size evolution along the pipe in micro gravity, but it is still limited to bubble smaller than the pipe diameter and should be extended to larger bubbles to predict the transition to slug flow.     

液氢空间贮存过程膜态沸腾数值模拟

为实现液氢在空间中安全高效应用,针对微重力条件下液氢膜态沸腾现象,建立了加热细丝浸没在过冷液氢池中的数值计算模型.采用VOF方法捕捉相界面,相变模型选取Lee模型,利用文献中的实验数据验证了模型的准确性.从气泡运动行为和换热特性两方面开展研究,结果发现液体过冷度和重力水平是影响换热机理的两个重要因素.在高重力水平、低液...     

Modeling the Effect of Bubbles on a Pattern and Heat Transfer in a Turbulent Polydisperse Upward Two-Phase Flow after Sudden Enlargement in a Tube

In this work, the numerical modeling of the flow pattern and heat transfer in a polydisperse bubbly turbulent flow after sudden enlargement in a tube is performed. The pattern of average and fluctuation twophase flows at small volumetric gas flow rate ratios (β ≤ 10%) is qualitatively similar to the one-phase liquid flow pattern. It is shown that small bubbles are present almost throughout the entire cross section of a tube, while great bubbles generally pass through the flow core and the shear mixing layer. The addition of air bubbles to a one-phase liquid flow appreciably intensifies heat transfer (up to two times), and these effects become stronger with an increase in the diameter of bubbles and the volumetric gas flow rate ratios gasratios.     

Numerical Simulation of Bubble Dynamics in Microgravity

A numerical method for the simulation of two-phase flows under microgravity conditions is presented in this paper. The level set method is combined with the moving mesh method in a collocated grid to capture the moving interfaces of the two-phase flow, and a SIMPLER-based method is employed to numerically solve the complete incompressible Navier-Stokes equations, and the surface tension force is modeled by a continuum surface force approximation. Based on the numerical results, the coalescence process of two bubbles under microgravity conditions (10^???2×g) is compared to that under normal gravity, and the effect of gravities on the bubbles coalescence dynamics is analyzed. It is showed that the velocity fields inside and around the bubbles under different gravity conditions are quite similar, but the strength of vortices behind the bubbles in the normal gravity is much stronger than that under microgravity conditions. It is also found that under microgravity conditions, the time for two bubbles coalescence is much longer, and the deformation of bubbles is much less, than that under the normal gravity.     

应用双头电导探针技术测量 气液两相泡状流局部参数

本文研究应用双头电导探针技术测量气泡局部参数,从而揭示了气液两相泡状流的内部流动规律。成功地设计了一种能够快速可靠测量气泡局部统计参数,包括空隙率、气泡速度、气泡尺寸、界面浓度等的电导探针系统。发现探针尖部的导通距离、沿流动方向两探针间的距离和两个探针针尖的间隙是设计电导探针的关键尺寸。     

Semianalytical solution of unsteady quasi-one-dimensional cavitating nozzle flows

Unsteady quasi-one-dimensional bubbly cavitating nozzle flows are considered by employing a homogeneous bubbly liquid flow model, where the nonlinear dynamics of cavitating bubbles is described by a modified Rayleigh–Plesset equation. The model equations are uncoupled by scale separation leading to two evolution equations, one for the flow speed and the other for the bubble radius. The initial-boundary value problem of the evolution equations is then formulated and a semianalytical solution is constructed. The solution for the mixture pressure, the mixture density, and the void fraction are then explicitly related to the solution of the evolution equations. In particular, a relation independent of flow dimensionality is established between the mixture pressure, the void fraction, and the flow dilation for unsteady bubbly cavitating flows in the model considered. The steady-state compressible and incompressible limits of the solution are also discussed. The solution algorithm is first validated against the numerical solution of Preston et al. [Phys Fluids 14:300–311, 2002] for an essentially quasi-one-dimensional nozzle. Results obtained for a two-dimensional nozzle seem to be in good agreement with the mean pressure measurements at the nozzle wall for attached cavitation sheets despite the observed two-dimensional cavitation structures.     

应用双头电导探针技术测量气液两相泡状流局部参数

本研究应用双头电导探针技术测量气泡局部参数,从而揭示了气液两相泡状流的内部流动规律。     

Predictions of Phase Distribution in Liquid–Liquid Two-Component Flow

Ground-based liquid–liquid two-component flow can be used to study reduced-gravity gas-liquid two-phase flows provided that the two liquids are immiscible with similar densities. In this paper, we present a numerical study of phase distribution in liquid–liquid two-component flows using the Eulerian two-fluid model in FLUENT, together with a one-group interfacial area transport equation (IATE) that takes into account fluid particle interactions, such as coalescence and disintegration. This modeling approach is expected to dynamically capture changes in the interfacial structure. We apply the FLUENT-IATE model to a water-Therminol 59^® two-component vertical flow in a 25-mm inner diameter pipe, where the two liquids are immiscible with similar densities (3% difference at 20°C). This study covers bubbly (drop) flow and bubbly-to-slug flow transition regimes with area-averaged void (drop) fractions from 3 to 30%. Comparisons of the numerical results with the experimental data indicate that for bubbly flows, the predictions of the lateral phase distributions using the FLUENT-IATE model are generally more accurate than those using the model without the IATE. In addition, we demonstrate that the coalescence of fluid particles is dominated by wake entrainment and enhanced by increasing either the continuous or dispersed phase velocity. However, the predictions show disagreement with experimental data in some flow conditions for larger void fraction conditions, which fall into the bubbly-to-slug flow transition regime. We conjecture that additional fluid particle interaction mechanisms due to the change of flow regimes are possibly involved.     

Effects of Reduced Gravity Conditions on Bubble Dispersion Characteristics in the Bubble Column

Bubble-liquid turbulent flow has an excellent heat and mass transfer behaviors than single gas or liquid flow. In order to analyze the effects of normal and reduced gravity on cold bubble-liquid two-phase turbulent flow in bubble column a second-order moment cold bubble-liquid two-phase turbulent model was developed to disclose the bubble dispersion characteristics. Under the reduced gravity condition, volume fraction caused by the decrease of buoyance force is larger than normal gravity level due to bigger bubble solid volume. In addition, bubble frequency is also decreased by in decrease of buoyance force. Normal and shear stresses have strongly anisotropic characteristics at every directions and have larger values under normal gravity than reduced gravity. The liquid turbulent kinetic energy has the two-peak bimodal distribution and weaker than bubble turbulent kinetic energy with one peak unimodal, which is caused by vigorous wake fluctuations. The correlation of fluctuation velocities between bubble and liquid has clearly anisotropic behaviors Under reduced gravity, the bubble motion has a little impact on liquid turbulent flow caused by slight buoyancy force, however, it will greatly reduce the liquid turbulent intensity due to energy cascade transport, which was transformed into bubbles or dissipated by interface friction. Bubble formation and detachment mechanisms affected by gravity conditions lead to the different levels of bubble dispersion distributions.     

Effects of higher wall roughness on dense particle–laden dispersion behaviors

To analyze the effects of higher wall roughness on dense particle–laden dispersion behaviors under reduced gravity environments, a dense gas–particle two-phase second-order-moment turbulent model are developed. In this model, the wall roughness function and the kinetic theory of granular flows are coupled and closed. Anisotropy of gas–solid two-phase stresses and the interaction between gas–particle are fully considered using two-phase Reynolds stress model and the two-phase correlation transport equation. Numerical simulation test is validated by Sommerfeld and Kussin (2003) experiments data with higher wall roughness 8.32 μm. Predicted results showed that the particle concentration distribution, particle fluctuation velocity, particle temperature and particle collision frequency are greatly affected by higher wall roughness, as well as particle Reynolds stress and interactions between gas and particle turbulent flows are redistributed. Under microgravity conditions, particle temperature and collision frequency are greatly less than those of earth and lunar gravity.     

The Influence of the Gravity Force on Longwave Marangoni Patterns in Two-Layer Films

An Euler–Euler two-fluid model based on the second-order-moment closure approach and the granular kinetic theory of dense gas-particle flows was presented. Anisotropy of gas-solid two-phase stress and the interaction between two-phase stresses are fully considered by two-phase Reynolds stress model and the transport equation of two-phase stress correlation. Under the microgravity space environments, hydrodynamic characters and particle dispersion behaviors of dense gas-particle turbulence flows are numerically simulated. Simulation results of particle concentration and particle velocity are in good agreement with measurement data under earth gravity environment. Decreased gravity can decrease the particle dispersion and can weaken the particle–particle collision as well as it is in favor of producing isotropic flow structures. Moreover, axial–axial fluctuation velocity correlation of gas and particle in earth gravity is approximately 3.0 times greater than those of microgravity and it is smaller than axial particle velocity fluctuation due to larger particle inertia and the larger particle turbulence diffusions.     

Performance analysis of no-vent fill process for liquid hydrogen tank in terrestrial and on-orbit environments

Two finite difference computer models, aiming at the process predictions of no-vent fill in normal gravity and microgravity environments respectively, are developed to investigate the filling performance in a liquid hydrogen (LH₂) tank. In the normal gravity case model, the tank/fluid system is divided into five control volume including ullage, bulk liquid, gas–liquid interface, ullage-adjacent wall, and liquid-adjacent wall. In the microgravity case model, vapor–liquid thermal equilibrium state is maintained throughout the process, and only two nodes representing fluid and wall regions are applied. To capture the liquid–wall heat transfer accurately, a series of heat transfer mechanisms are considered and modeled successively, including film boiling, transition boiling, nucleate boiling and liquid natural convection. The two models are validated by comparing their prediction with experimental data, which shows good agreement. Then the two models are used to investigate the performance of no-vent fill in different conditions and several conclusions are obtained. It shows that in the normal gravity environment the no-vent fill experiences a continuous pressure rise during the whole process and the maximum pressure occurs at the end of the operation, while the maximum pressure of the microgravity case occurs at the beginning stage of the process. Moreover, it seems that increasing inlet mass flux has an apparent influence on the pressure evolution of no-vent fill process in normal gravity but a little influence in microgravity. The larger initial wall temperature brings about more significant liquid evaporation during the filling operation, and then causes higher pressure evolution, no matter the filling process occurs under normal gravity or microgravity conditions. Reducing inlet liquid temperature can improve the filling performance in normal gravity, but cannot significantly reduce the maximum pressure in microgravity. The presented work benefits the understanding of the no-vent fill performance and may guide the design of on-orbit no-vent fill system.