The influence of film structure on the interfacial friction in annular two-phase flow under microgravity and normal gravity conditions

Results for the interfacial friction factor and relative interfacial roughness on the gas-liquid interface are reported for an air-water annular flow in a small inner diameter tube (9.53 mm i.d.). The film structure was obtained through processing the time trace signal of film thickness measurements using conductance probes. The interfacial friction factor and the wave height were altered through changing the gravity level and gas Reynolds number. It was found that the wave height decreased with increasing the gas Reynolds number. The wave height in microgravity is less than half of that in normal gravity, while the friction factor was about 10% smaller in microgravity than that in normal gravity. It was shown that the annular two-phase flow friction factor decreased less dramatically as the relative interfacial roughness decreased compared to the single-phase case. It is interesting to note that the interfacial shear stress values at microgravity were very close (or even larger than) those at normal gravity. This was attributed to the thicker substrate at microgravity.     

Particle accumulation structures in time-dependent thermocapillary flow in a liquid bridge under microgravity

The formation of a dynamic spiral string of particles with larger density than the fluid was investigated for time-dependent thermocapillary flow in liquid bridges under various gravity conditions including microgravity. The dynamic spiral string forms after approximately 20–60 oscillation periods from the homogeneous dilute particle suspension. It was found that the action of gravity is not decisive in the process of the particle accumulation structure (PAS) but gravity influences the flow field for PAS-formation. We could realize and observe PAS with modal structure m=3 under μ-g but modal structure m=2 occurred only during a transient of the operating parameters in an aspect ratio-range different from that under normal gravity. The correlation of the optically observed PAS structure with the temperature structure of the azimuthally rotating hydrothermal wave on the free surface is the same under microgravity as under normal gravity indicating that PAS is a pure Marangoni effect.     

Flow physics of upward cocurrent gas-liquid annular flow in a vertical small diameter tube

In this paper, a steady RNG k-ε model, in conjunction with enhanced wall treatment method, was applied to the gas core in order to simulate the flow physics of annular two-phase flow. The model incorporated a physical model of wave characteristics and included the liquid entrainment influence on the flow. Based on the simulation results, flow features in the gas core were quantitatively presented and a model of the liquid entrainment mechanism was proposed. In addition, a parametric study was conducted to determine the impact of changing wave velocity, pressure, and gravitational force on the liquid film flow. The results were validated using a large set of experimental data at normal and microgravity conditions. Also, the law of the wall was applied to previously-collected experimental data. Analysis yielded different flow features of the liquid film at microgravity and normal gravity conditions.     

A Computational Study of Soot Formation in Methane Air Co-Flow Diffusion Flame Under Microgravity Conditions

An in-house developed code has been used to predict soot formation in a methane air co flow diffusion flame at normal gravity and at lower gravity levels of 0.5 G, and 0.0001 G (microgravity). There is an augmentation of soot formation at lower gravity levels because of lower buoyancy induced acceleration leading to an increased residence time. The peak temperature at microgravity is reduced by about 50 K than that at normal gravity level. The axial velocity under normal gravity and reduced gravity show negative values (relatively small in magnitude) near the wall at axial height beyond 15 cm; but axial velocity is never negative in microgravity condition. Peak value of soot volume fraction at 0.5 G and microgravity multiplies by a factor of ～3 and ～7, respectively of that at normal gravity. The zone of peak soot volume fraction shifts away from the axis towards the wings, as gravity level is lowered. In comparison to soot volume fraction, the factors of amplification of soot number density at reduced gravity and at microgravity are comparatively lower at 1.2 and 1.5 of that at normal gravity respectively. On the other hand, mean soot particle sizes at reduced gravity and microgravity increase to 1.5 and 2 times of that at normal gravity respectively.     

Experimental Study of Subcooled Flow Boiling Heat Transfer on a Smooth Surface in Short-Term Microgravity

The flow boiling heat transfer characteristics of subcooled air-dissolved FC-72 on a smooth surface (chip S) were studied in microgravity by utilizing the drop tower facility in Beijing. The heater, with dimensions of 40 × 10 × 0.5 mm³ (length × width × thickness), was combined with two silicon chips with the dimensions of 20 × 10 × 0.5 mm³. High-speed visualization was used to supplement observation in the heat transfer and vapor-liquid two-phase flow characteristics. In the low and moderate heat fluxes region, the flow boiling of chip S at inlet velocity V =?0.5 m/s shows almost the same regulations as that in pool boiling. All the wall temperatures at different positions along the heater in microgravity are slightly lower than that in normal gravity, which indicates slight heat transfer enhancement. However, in the high heat flux region, the pool boiling of chip S shows much evident deterioration of heat transfer compared with that of flow boiling in microgravity. Moreover, the bubbles of flow boiling in microgravity become larger than that in normal gravity due to the lack of buoyancy Although the difference of the void fraction in x-y plain becomes larger with increasing heat flux under different gravity levels, it shows nearly no effect on heat transfer performance except for critical heat flux (CHF). Once the void fraction in y-z plain at the end of the heater equals 1, the vapor blanket will be formed quickly and transmit from downstream to upstream along the heater, and CHF occurs. Thus, the height of channel is an important parameter to determine CHF in microgravity at a fixed velocity. The flow boiling of chip S at inlet velocity V =?0.5 m/s shows higher CHF than that of pool boiling because of the inertia force, and the CHF under microgravity is about 78–92% of that in normal gravity.     

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.     

A model of interfacial waves in annular two-phase flow at different gravity levels

Microgravity provides an ideal environment to study interfacial waves (small scale phenomena) without the dominating effects of gravity. In this paper, conductance probes were used to measure film thickness and liquid velocities. Physical models of interfacial waves are presented for vertical upward flow (+1g), vertical downward flow (?1g), and microgravity (μg), at different gas and liquid flow rates. Based on these studies, the predominant parameters and the features of the wave structure are discussed. A preliminary mathematical model was proposed for the interfacial waves.     

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.     

Simulation of Jet-Induced Geysers in Reduced Gravity

Control of cryogenic propellant tank pressure during tank refueling and expulsion in low gravity is an important technical challenge to overcome for future long duration missions in space. One method proposed to control tank pressurization involves the use of jet-induced geysers. Two-dimensional computational models have been developed and used with limited success in previous efforts to predict geyser heights in microgravity. A three-dimensional flow simulation is used to model jet-induced geysers in reduced gravity. Geyser flows are commonly characterized by the presence of turbulent jets, transient flow, deforming free surfaces, and surface tension effects. As is the case for many turbulent flow applications, accuracy in simulating complex turbulent flows is critically dependent on the selection of a suitable turbulence model. The sensitivity of the simulation geyser predictions to a suite of popular turbulence models is assessed. Simulation results are compared to available experiment results. By expanding upon the work already completed, the model is used to simulate a broad range of cases within the experiment test matrix. Simulation results suggest the two dimensional simulation using the k-ε turbulence model provides the most accurate results for jet-induced geysers in reduced gravity when compared to available experiment data.     

Convective heat transfer coefficients for near-horizontal two-phase flow of nitrogen and hydrogen at low mass and heat flux

Correlations for convective heat transfer coefficients are reported for two-phase flow of nitrogen and hydrogen under low mass and heat flux conditions. The range of flowrates, heat flux and tube diameter are representative of thermodynamic vent systems (TVSs) planned for propellant tank pressure control in spacecraft operating over long durations in microgravity environments. Experiments were conducted in normal gravity with a 1.5° upflow configuration. The Nusselt number exhibits peak values near transition from laminar to turbulent flow based on the vapor Reynolds number. This transition closely coincides with a flow pattern transition from plug to slug flow. The Nusselt number was correlated using components of the Martinelli parameter and a liquid-only Froude number. Separate correlating equations were fitted to the laminar liquid/laminar vapor and laminar liquid/turbulent vapor flow data. The correlations give root-mean-squared (rms) prediction errors within 15%.     

Sloshing in microgravity

The International Space Station provides a low gravity environment for experiments that require very low acceleration. The steady component of acceleration due to the gravity gradient is in the microgravity range. It is possible to achieve microgravity levels for the variable component by using isolation racks. For experiments cooled by liquid cryogens sloshing may increase the variable acceleration at the experiment beyond acceptable levels. Sloshing of cryogens in microgravity can be predicted using a surface wave model. The model should include: a calculation of the shape of the unperturbed liquid–gas interface; a listing of the normal modes and resonant frequencies for the container; a prediction of the amplitude of the modes in response to the motion of the container; and a test to detect the breakdown of linear theory. A model is presented that contains these components. The shape of the interface is calculated and it is found that for most anticipated applications the interface is nearly cylindrical or spherical. Since gravity is not aligned with the symmetry of the container, the depth of the liquid is variable. Examples are presented to show how to estimate the extent of variable depth and curved interface on the normal modes and resonant frequencies. Equations are derived for the dynamic interaction of the isolation rack, the dewar and the sloshing motion. Damping is introduced by using boundary layer theory. Random vibration theory is applied to the incoherent component of the driving spectrum while standard resonance formalism is used for the coherent component. The model cannot be used if the wave amplitude becomes so large that linear theory does not apply. A procedure is developed to check for nonlinear difficulties.     

Numerical Simulation of Transient Development of Flame,Temperature and Velocity under Reduced Gravity in a Methane Air Diffusion Flame

A methane air co flow diffusion flame has been numerically simulated with the help of an in-house developed code at normal gravity, 0.5 G, and 0.0001 G (microgravity) for the study of transient behavior of the flame in terms of flame shape, temperature profile and velocity (streamlines). The study indicates that lower is the gravity level, the higher is the time of early transience. The flame developments during transience are marked by the formation of a secondary flamelet at different heights above the primary flame at all gravity levels. The development of temperature profile at microgravity takes a much longer time to stabilize than the flame development. At normal gravity and 0.5 G gravity level, streamlines, during transience, show intermediate vortices which are finally replaced by recirculation of ambient air from the exit plane. At microgravity, neither any vortex nor any recirculation at any stage is observed. Centerline temperature plots, at all gravity levels during transience, demonstrate a secondary peak at some instants as a consequence of the secondary flamelet formation. The centerline velocity at microgravity decreases gradually during transience, unlike at other two gravity levels where the fall is very sharp and is indicative of negligible buoyancy at microgravity.     

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.     

Hydrothermal wave in a shallow liquid layer

The oscillatory thermocapillary convection and hydrothermal wave in a shallow liquid layer, where a temperature difference is applied between two parallel sidewalls, have been numerically investigated in a two-dimensional model. The oscillatory thermocapillary convection and hydrothermal wave appear if the Marangoni number is larger than a critical value. The critical phase speed and critical wave number of the hydrothermal wave agree with the ones given analytically by Smith and Davis in the microgravity environment, and it travels in the direction opposed to the surface flow. Another wave traveled downstream in addition to the hydrothermal wave traveled upstream was observed in the case of earth gravity condition.     

Nonlinear Wave Simulation on a Surface of Liquid Film Entrained by Turbulent Gas Flow at Weightlessness

The article considers the joint flow of a liquid film entrained by turbulent gas. Since liquid velocity is small in comparison with gas velocity, the problem is reduced to the calculation of pressure and shear stresses produced by the gas flowing over a wavy wall with small amplitude. Further, these data are used at the boundary conditions, when the flow of a liquid film is considered separately. As a result, we obtain a new system of equations for modeling the dynamics of long-wave perturbations on the surface of a viscous liquid film entrained by turbulent gas at microgravity. At small Reynolds numbers typical to the condition of microgravity, it was proved that this system is reduced to one evolution equation for the film thickness. Some numerical solutions of this equation have been received in this work.     

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.     

Granular Flow Under Microgravity: A Preliminary Review

The complex macroscopic rheological behavior of granular flow contains elements of both solid and liquid flow. Furthermore, under microgravity, granular flow exhibits novel flow features. To overcome a lack of comprehensive analyses of granular flow under microgravity, this study reviews the microgravity platforms and devices under which granular flow can be observed, the experimental findings made in such settings, and the range of numerical simulations that can be used to examine granular flow under microgravity. Differences in experimental research between normal gravity and microgravity are highlighted. These differences are found in the modifications made to conventional granular flow experimental devices, in new or unique granular flow behaviors, and in the numerical simulation methods needed for microgravity modeling. Additionally, the benefits of numerical simulation methods for examining rapid and dense flows under microgravity are also discussed. This study may have wide-ranging implications in such fields as investigations of the surface geology of asteroids or the efficient design and development of anchoring mechanisms or space vehicles.     

Mathematical and Numerical Modelling of Microconvection in the Long Rectangular Domains

The thermal gravitational convection of liquids under conditions of microgravity is studied in a 2D long rectangular domain elongated in the direction of the gravity force. The liquid is located between two solid regions of equal thickness. The solid parts are heat conducting. The mathematical modelling of the coupled problem is presented. Two mathematical models of convection are used to describe a motion of liquid: the classical Oberbeck-Boussinesq model and the microconvection model of isothermally incompressible liquid. The numerical experiments of convection are performed and demonstrate the qualitative and quantitative differences in the flow characteristics.     

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.     

Effects of Gravity on Ignition and Combustion Characteristics of Externally Heated Polyethylene Film

The objective of this research is to investigate the effects of gravity on the ignition and the combustion characteristics of the Polyethylene (PE) film by outer heating. Combustion experiments of PE film were carried out in a normal gravity field and the microgravity field. In the microgravity experiments, it was carried out in 50 m-class drop facility. Here it can be realized 10^??4G microgravity field in about 2.5-3.0 second. The PE film is heated by the inserted high-temperature chamber. In the experiments, the PE was used film type. The chamber temperature was fixed at 900 K and 1000 K. In the case of microgravity field, the ignition delay period has become about 50 percent shorter than that in the case of the normal gravitational field. In the normal gravity field, since the PE surface layer is cooled by natural convection, the ignition delay period is considered to be longer than that in the microgravity field. The combustion time in the normal gravity was about 0.8 sec. In the microgravity field, the combustion time was more than 2 sec, and it could not be measured during the free fall period.