Direct numerical simulation of a liquid sheet in a compressible gas stream in axisymmetric and planar configurations

A thin liquid sheet present in the shear layer of a compressible gas jet is investigated using an Eulerian approach with mixed-fluid treatment for the governing equations describing the gas–liquid two-phase flow system, where the gas is treated as fully compressible and the liquid as incompressible. The effects of different topological configurations, surface tension, gas pressure and liquid sheet thickness on the flow development of the gas–liquid two-phase flow system have been examined by direct solution of the compressible Navier–Stokes equations using highly accurate numerical schemes. The interface dynamics are captured using volume of fluid and continuum surface force models. The simulations show that the dispersion of the liquid sheet is dominated by vortical structures formed at the jet shear layer due to the Kelvin–Helmholtz instability. The axisymmetric case is less vortical than its planar counterpart that exhibits formation of larger vortical structures and larger liquid dispersion. It has been identified that the vorticity development and the liquid dispersion in a planar configuration are increased at the absence of surface tension, which when present, tends to oppose the development of the Kelvin–Helmholtz instability. An opposite trend was observed for an axisymmetric configuration where surface tension tends to promote the development of vorticity. An increase in vorticity development and liquid dispersion was observed for increased liquid sheet thickness, while a decreasing trend was observed for higher gas pressure. Therefore surface tension, liquid sheet thickness and gas pressure factors all affect the flow vorticity which consequently affects the dispersion of the liquid.      

Application of maximum entropy method for droplet size distribution prediction using instability analysis of liquid sheet

This paper describes the implementation of the instability analysis of wave growth on liquid jet surface, and maximum entropy principle (MEP) for prediction of droplet diameter distribution in primary breakup region. The early stage of the primary breakup, which contains the growth of wave on liquid–gas interface, is deterministic; whereas the droplet formation stage at the end of primary breakup is random and stochastic. The stage of droplet formation after the liquid bulk breakup can be modeled by statistical means based on the maximum entropy principle. The MEP provides a formulation that predicts the atomization process while satisfying constraint equations based on conservations of mass, momentum and energy. The deterministic aspect considers the instability of wave motion on jet surface before the liquid bulk breakup using the linear instability analysis, which provides information of the maximum growth rate and corresponding wavelength of instabilities in breakup zone. The two sub-models are coupled together using momentum source term and mean diameter of droplets. This model is also capable of considering drag force on droplets through gas–liquid interaction. The predicted results compared favorably with the experimentally measured droplet size distributions for hollow-cone sprays.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.     

Experimental study of the gas flow in an air-blasted liquid sheet

An experimental study has been performed to improve the understanding of the initial air–liquid interaction in the near field of an air-blasted breaking water sheet. For the first time, planar laser-induced fluorescence (PLIF) has been used to visualize the air-flow field, seeding the air streams with acetone vapor. Mie scattering from the liquid sheet, together with the acetone fluorescence signal has enabled simultaneous determination of the instantaneous water sheet location and the air-flow structures. The two-phase flow visualization has revealed detachment of the air boundary layer over the air–water interface behind the zones of strong curvature. The pressure field induced by these vortices has been identified as a cause of the enhanced sheet flapping and the instability growth. Received: 30 October 2000/Accepted: 29 March 2001     

Flow visualization of shock/water column interactions

The breakup of a liquid droplet induced by a high speed gas stream is a typical multiphase flow problem. The shock/droplet interaction is the beginning stage of the droplet breakup. Therefore, investigation of the shock/droplet interactions would be a milestone for interpreting the mechanism of the droplet breakup. In this study, a compressible multiphase solver with a five-equation model is successfully developed to study shock/water column interactions. For code validation, interface-only, gas–gas shock tube, and gas–liquid shock tube problems are first computed. Subsequently, a planar shock wave interacting with a water column is simulated. The transmitted wave and the alternative appearances of local high- and low-pressure regions inside the water column are observed clearly. Finally, a planar shock wave interacting with two water columns is investigated. In this work, both horizontal and vertical arrangements of two water columns are studied. It is found that different arrangements can result in the diversity of the interacting process. The complex flow structures generated by shock/water column interactions are presented by flow-visualization techniques.      

Turbulent flow simulation of liquid jet emanating from pressure-swirl atomizer

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) for the simulation of turbulent flow inside and outside pressure-swirl atomizer are evaluated by comparing the interface position with available experimental data and by comparing the turbulence intensity profiles at the atomizer exit. It is found that the turbulence models investigated exhibit zonal behaviors, i.e. none of the models investigated performs well throughout the entire flow field. The turbulence intensity has a significant influence on the global characteristics of the flow field. The turbulence models with better predictions of the turbulence intensity, such as Gatski-Speziale’s ARSM model, can yield better predictions of the global characteristics of the flow field, e.g. the reattachment lengths for the backward-facing step flow and the sudden expansion pipe flow, or the discharge coefficient, film thickness and the liquid sheet outer surface position for the atomizer flows. The standard k–ε model predicts stronger turbulence intensity as compared to the other models and therefore yields smaller film thickness and larger liquid sheet outer surface position. In average, the ARSM model gives both quantitatively and qualitatively better results as compared to the standard k–ε model and the low Reynolds number models.     

Molecular dynamics simulation of annular flow boiling with the modified Lennard-Jones potential function

Molecular dynamics simulation of annular flow boiling in a nanochannel is numerically investigated. In this research, an annular flow model is developed to predict the superheated flow boiling heat transfer characteristics in a nanochannel. To characterize the forced annular boiling flow in a nanochannel, an external driving force F^?_\textext \overrightarrow {F}_{\text{ext}} ranging from 1 to 12 PN (PN = pico newton) is applied along the flow direction to inlet fluid particles during the simulation. Based on an annular flow model analysis, it is found that saturation condition and superheat degree have great influences on the liquid–vapor interface. Also, the results show that due to the relatively strong influence of the surface tension in small channels, the interface between the liquid film and the vapor core is fairly smooth, and the mean velocity along the stream-wise direction does not change anymore. Also, it is found that the heat flux values depend on the boundary conditions. Finally, the Green–Kubo formula is used to calculate the thermal conductivity of liquid Argon. The simulations predict thermal conductivity of liquid Argon quite well.     

Foam Self-Clarification Phenomenon: An Experimental Investigation

This paper is focused on the capabilities of gas–liquid foams to attenuate acoustic waves. It is postulated that the sound attenuation phenomenon in foams is largely governed by the hydrodynamic resistance of the Plateau-Gibbs channels (PGC) to the flow of liquid through them. It is shown that the addition of solid particles to gas–liquid foams has opposite effects depending on the concentration of the added solid particles. As long as the concentration of the added solid particles is smaller than a certain critical value the sound attenuation coefficient increases and as a result in the sound velocity decreases. However, if the concentration of the added solid particles becomes larger than this critical value the attenuation coefficient decreases and the sound velocity increases. When the concentration of the solid particles reaches some critical value, the particles block the Plateau-Gibbs channels and stop the filtration. As a result the attenuation coefficient of the sound wave decreases while the sound velocity, in such three-phase foams, increases. The point at which the sound wave stops attenuating and its velocity starts to increase is known as the point of self-clarification. Based on this postulate and on the results of our preliminary tests the present study provides a plausible explanation to the above-mentioned contradicting effect, and the self- clarification phenomenon.     

Buoyancy-driven motion of a two-dimensional bubble or drop through a viscous liquid in the presence of a vertical electric field

The effect of an electric field on the buoyancy-driven motion of a two-dimensional gas bubble rising through a quiescent liquid is studied computationally. The dynamics of the bubble is simulated numerically by tracking the gas–liquid interface when an electrostatic field is generated in the vertical gap of the rectangular enclosure. The two phases of the system are assumed to be perfect dielectrics with constant but different permittivities, and in the absence of impressed charges, there is no free charge in the fluid bulk regions or at the interface. Electric stresses are supported at the bubble interface but absent in the bulk and one of the objectives of our computations is to quantify the effect of these Maxwell stresses on the overall bubble dynamics. The numerical algorithm to solve the free-boundary problem relies on the level-set technique coupled with a finite-volume discretization of the Navier–Stokes equations. The sharp interface is numerically approximated by a finite-thickness transition zone over which the material properties vary smoothly, and surface tension and electric field effects are accounted for by employing a continuous surface force approach. A multi-grid solver is applied to the Poisson equation describing the pressure field and the Laplace equation governing the electric field potential. Computational results are presented that address the combined effects of viscosity, surface tension, and electric fields on the dynamics of the bubble motion as a function of the Reynolds number, gravitational Bond number, electric Bond number, density ratio, and viscosity ratio. It is established through extensive computations that the presence of the electric field can have an important effect on the dynamics. We present results that show a substantial increase in the bubble’s rise velocity in the electrified system as compared with the corresponding non-electrified one. In addition, for the electrified system, the bubble shape deformations and oscillations are smaller, and there is a reduction in the propensity of the bubble to break up through increasingly larger oscillations.     

Characteristics of liquid sheets formed by splash plate nozzles

An experimental study was conducted to identify the effect of viscosity on the characteristics of liquid sheets formed by a splash plate nozzle. Various mixtures of corn syrup and water are used to obtain viscosities in the range 1–170 mPa.s. Four different splash plates with nozzle diameters of 0.5, 0.75, 1, and 2 mm, with a constant plate angle of 55° were tested. Liquid sheets formed under various operating conditions were directly visualized. The sheet atomization process for the range of parameters studied here is governed by two different mechanisms: Rayleigh–Plateau (R–P) and Rayleigh–Taylor (R–T) instabilities. R–P occurs at the rim and R–T occurs on the thin sheet. The rim instability can be laminar or turbulent, depending on the jet Reynolds number. The R–T instability of the sheet is observed at the outer edges of the radially spreading sheet, where the sheet is the thinnest. It can also occur inside the sheet, due to formation of holes and ruptures.     

Linear Stability Analysis of a Viscous Liquid Sheet in a High-Speed Viscous Gas

Air-assisted atomizers in which a thin liquid sheet is deformed under the action of a high-speed air flow are extensively used in industrial applications, e.g., in aircraft turbojet injectors. Primary atomization in these devices is a consequence of the onset and growth of instabilities on the air/liquid interfaces. To better understand this process, a temporal linear instability analysis is applied to a thin planar liquid sheet flowing between two semi-infinite streams of a high-speed viscous gas. This study includes the full viscous effects both in the liquid and gas basic states and perturbations. The relevant dimensionless groups entering the non-dimensional Orr–Sommerfeld equations and boundary conditions are the liquid and gas stream Reynolds numbers, the gas to liquid momentum flux ratio, the gas/liquid velocity ratio, the Weber number and the equivalent gas boundary layer to liquid sheet thickness ratio. Growth rates and temporal frequencies as a function of the wave number, varying the different dimensionless parameters are presented, together with neutral stability curves. From the results of this parametric study it is concluded that when the physical properties of gas and liquid are fixed, the momentum flux ratio is especially relevant to determine the instability conditions. It is also observed that the gas boundary layer thickness strongly influences the wave propagation, and acts by damping sheet oscillation frequency and growth. This is especially important because viscosity in the basic gas velocity profile has always been ignored in instability analysis applied to the geometry under study. This revised version was published online in July 2006 with corrections to the Cover Date.     

Startup analysis of oil-fired furnace – the smoothing Monte Carlo model approach

 A multi-dimensional mathematical thermal model of the oil-fired furnace has been developed. Radiation calculation of this model is based on the Hottel's zone method incorporating with the statistical Monte Carlo method to determine the total exchange areas. The whole furnace is divided into 4 gaseous zones and 18 surface zones. The Monte Carlo method integrated with the least square smoothing technique considering both conservation and reciprocity is used to evaluate the total exchange areas, i.e. surface–surface, surface–gas, gas–gas and gas–surface, directly for both the absorbing and emitting media. A better accuracy is achieved in the determination of the total exchange areas by using the proposed smoothing Monte Carlo method. Received on 10 April 2000     

Axisymmetric thermocapillary convection in a slowly rotating annular cylindrical container

In a slowly rotating annular cylindrical container the free liquid surface (liquid-gas interface) is subjected to a temperature gradient in radial direction. The temperature dependent surface tension creates a shear stress on the interface which is transmitting a thermocapillary convection in the bulk of the liquid. For constant temperature T ₁ of the inner and T ₂ of the outer wall a steady Marangoni convection takes place, exhibiting a double vortex ring of equal directional flow. For time-oscillatory temperatures of the walls a time-dependent thermocapillary convection appears, which will create on the free liquid surface various wave patterns. They shall, depending on the forcing frequency of the temperature, exhibit resonance peaks. The velocity distribution and the response magnitude inside the container has been determined. Received on 3 September 1997     

Thermocapillary convection in an annular cylindrical container

The liquid–gas interface of a liquid in an annular container is subjected to a temperature gradient. Since shear stress on the free liquid surface depends on the temperature it shall create a variable shear stress on the surface which in turn yields by viscous traction a thermocapillary convection in the bulk of the liquid. For constant temperature T ₂ > T ₁ on the outer cylindrical wall and T ₁ on the inner wall a steady convection shall emerge. The morphology of the ensuing flow is investigated as a function of the diameter ratio and the liquid height ratio and exhibits the growth of bottom rim vortex rings, which unite to single vortex rings of opposite flow direction to the original single strong vortex ring above it. The results also show how such a system depends on the magnitude of the diameter ratio and the liquid height ratio. The evaluation of the results exhibits the evolution of vortex rings.     

Liquid film characterization in horizontal,annular, two-phase,gas–liquid flow using time-resolved laser-induced fluorescence

A non-intrusive optical technique was developed to provide time-resolved longitudinal and cross-sectional images of the liquid film in horizontal annular pipe flow of air and water, revealing the interfacial wave behavior. Quantitative information on the liquid film dynamics was extracted from the time-resolved images. The planar laser-induced fluorescence technique was utilized to allow for optical separation of the light emitted by the film from that scattered by the air–water interface. The visualization test section was fabricated from a tube presenting nearly the same refractive index as water, which allowed the visualization of the liquid film at regions very close to the pipe wall. Longitudinal images of the liquid film were captured using a high-frame-rate digital video camera synchronized with a high-repetition-rate laser. An image processing algorithm was developed to automatically detect the position of the air–water interface in each image frame. The thickness of the liquid film was measured at two axial stations in each processed image frame, providing time history records of the film thickness at two different positions. Wave frequency information was obtained by analyzing the time-dependent signals of film thickness for each of the two axial positions recorded. Wave velocities were measured by cross-correlating the amplitude signals from the two axial positions. For the film cross-section observations, two high-speed digital video cameras were used in a stereoscopic arrangement. Comparisons with results from different techniques available in literature indicate that the technique developed presents equivalent accuracy in measuring the liquid film properties. Time-resolved images of longitudinal and cross-section views of the film were recorded, which constitute valuable information provided by the technique implemented.     

Application of a borescope to studies of gas–liquid flow in downward inclined pipes

The aim of the present study is to investigate stratified downward gas–liquid pipe flow with a non-intrusive measurement technique that is based on a borescope connected to a digital video camera. The borescope-based technique enables to determine the instantaneous cross-sectional distribution of both phases within the pipe. Water and air were used as working fluids. Quantitative data was extracted from sequences of recorded video images by applying a developed data processing technique for instantaneous gas–liquid interface boundaries determination. Experiments were performed for a wide range of downward pipe inclinations and gas and liquid flow rates. The instantaneous and time-average cross-sectional holdup for each set of flow parameters was calculated. Particular attention was given to the study of the interface shape that in many occasions was not flat and was characterized by the penetration of the liquid along the pipe periphery. Temporal variation of the surface elevation was also studied and various regimes characterizing interfacial waves were defined using both the recorded time series of the instantaneous depth of the water layer and the Fourier analysis of those records.     

On the Mutual Penetrations of Two Fluids Whose Interface is Accelerated by a Shock Wave

A shock tube experimental investigation and numerical simulations are undertaken to study the evolution of a perturbed interface of two different gases accelerated by a shock wave. The experimental method is based on a high-speed camera laser sheet diagnostic technique, and simulations are provided by our code CARBUR based on a finite volume discretization of Navier–Stokes’s equations. Two gas pairs are used to illustrate both the heavy/light (air/He) and the light/ heavy (air/SF₆) cases. Two simultaneous large initial perturbations, one positive and one negative, are tested for an incident shock wave Mach number in air of about 1.3. The thin membrane (less than 1 μ) which materializes the initial interface between the two test gases presents 2D perturbations whose wave number is close to 1 in order to rapidly reach the non-linear regime. The development of the perturbations is captured at a frequency of 10 kHz after the interface acceleration, and the experiments are complemented with a numerical simulation to validate the interface deformations. Results show an asymmetric mutual gas penetration increasing with the absolute value of the Atwood’s number. Furthermore, they confirm that the heavier gas penetrates the lighter as thin spikes and the lighter gas penetrates the heavier as large bubbles. Moreover, we show that the spike moves faster than the bubble in the heavy/light case and slightly faster in the light/heavy one. Finally, numerical and experimental results are in agreement.     

圆环旋转粘性液膜射流三维不稳定性研究

利用线性稳定性理论进行了射流液体粘性对圆环旋转液膜射流稳定性影响的研究,推导出了三维扰动下具有固体旋涡型速度分布的圆环旋转粘性液膜射流的色散方程;在此基础上进行了类反对称模式与类对称模式下的圆环旋转粘性液膜射流的三维不稳定性分析。研究结果表明,在类反对称模式下,液体粘性超过一定值后,射流最大扰动增长率随液体粘性的增加而迅速减小;轴对称模态的射流特征频率产生一个突降变化;随液体粘性增加,轴对称模态不稳定波数范围减小,非轴对称模态不稳定波数范围呈现出先减小后增大趋势。在类对称模式下,液体粘性对射流最大扰动增长率的影响主要体现在对非轴对称模态的影响上;液体粘性只在粘性较大时才会对非轴对称模态射流特征频率产生一定影响;液体粘性超过一定值后,轴对称模态与非轴对称模态的不稳定波数范围都会快速下降。     

Marangoni convection in an annular cylindrical tank

If the free liquid-gas interface of a liquid in an annular container is subjected to a temperature gradient the temperature dependent shear stress on the free liquid surface creates by viscous traction a thermocapillary convection in the bulk of the liquid. For constant temperature T ₁ of the inner and T ₂ of the outer wall a steady Marangoni convection takes place, while for time-oscillatory temperatures of the walls a time-dependent thermocapillary convection appears, which will create on the free liquid surface wave patterns. They shall, depending on the forcing frequency of the temperature, exhibit resonance peaks. The velocity distribution and the response magnitude inside the container has been determined. Received on 11 July 1997