Multiscale simulations and experiments on water jet atomization

Numerical simulation of primary atomization at high Reynolds number is still a challenging problem. In this work a multiscale approach for the numerical simulation of liquid jet primary atomization is applied, using an Eulerian-Lagrangian coupling. In this approach, an Eulerian volume of fluid (VOF) method, where the Reynolds stresses are closed by a Reynolds stress model is applied to model the global spreading of the liquid jet. The formation of the micro-scale droplets, which are usually smaller than the grid spacing in the computational domain, is modelled by an energy-based sub-grid model. Where the disruptive forces (turbulence and surface pressure) of turbulent eddies near the surface of the jet overcome the capillary forces, droplets are released with the local properties of the corresponding eddies. The dynamics of the generated droplets are modelled using Lagrangian particle tracking (LPT). A numerical coupling between the Eulerian and Lagrangian frames is then established via source terms in conservation equations. As a follow-up study to our investigation in Saeedipour et al. (2016a), the present paper aims at modelling drop formation from liquid jets at high Reynolds numbers in the atomization regime and validating the simulation results against in-house experiments. For this purpose, phase-Doppler anemometry (PDA) was used to measure the droplet size and velocity distributions in sprays produced by water jet breakup at different Reynolds numbers in the atomization regime. The spray properties, such as droplet size spectra, local and global Sauter-mean drop sizes and velocity distributions obtained from the simulations are compared with experiment at various locations with very good agreement.     

The atomization of water–oil emulsions

The paper presents the results of experimental studies on atomization of the emulsions flowing through twin-fluid atomizers obtained by the use of the digital microphotography method. The main elements of the test installation were: nozzle, reservoir, pump and measurement units of liquid flow. The photographs were taken by a digital camera with automatic flash at exposure time of 1/8000 s and subsequently analyzed using Image Pro-Plus. The oils used were mineral oils 20–90, 20–70, 20–50 and 20–30. The studies were performed at flow rates of liquid phase changed from 0.0014 to 0.011 (dm³/s) and gas phase changed from 0.28 to 1.4 (dm³/s), respectively. The analysis of photos shows that the droplets being formed during the liquid atomization have very different sizes. The smallest droplets have diameters of the order of 10 μm. The experimental results showed that the changes in physical properties of a liquid phase lead to the significant changes in the spray characteristics. The analysis of the photos of water and emulsions atomization process showed that the droplet sizes are dependent on gas and liquid flow rates, construction of nozzle and properties of liquid. The differences between characteristics of atomization for water and emulsions have been observed. Analysis of photos on forming the droplets in air–water and air-emulsions systems showed that droplets are bigger in air-emulsion system (at the same value of gas to liquid mass ratio). The values of Sauter mean diameter (SMD) increased with increase of volume fraction of oil in emulsion. The droplet size increased with emulsion viscosity.     

An Eulerian–Lagrangian hybrid model for the coarse-grid simulation of turbulent liquid jet breakup

In this paper we present a numerical model for the coarse-grid simulation of turbulent liquid jet breakup using an Eulerian–Lagrangian coupling. To picture the unresolved droplet formation near the liquid jet interface in the case of coarse grids we considered a theoretical model to describe the unresolved flow instabilities leading to turbulent breakup. These entrained droplets are then represented by an Eulerian–Lagrangian hybrid concept. On the one hand, we used a volume of fluid method (VOF) to characterize the global spreading and the initiation of droplet formation; one the other hand, Lagrangian droplets are released at the liquid–gas interface according to the theoretical model balancing consolidating and disruptive energies. Here, a numerical coupling was required between Eulerian liquid core and Lagrangian droplets using mass and momentum source terms. The presented methodology was tested for different liquid jets in Rayleigh, wind-induced and atomization regimes and validated against literature data. This comparison reveals fairly good qualitative agreement in the cases of jet spreading, jet instability and jet breakup as well as relatively accurate size distribution and Sauter mean diameter (SMD) of the droplets. Furthermore, the model was able to capture the regime transitions from Rayleigh instability to atomization appropriately. Finally, the presented sub-grid model predicts the effect of the gas-phase pressure on the droplet sizes very well.     

Transient high-frequency ultrasonic water atomization

An experimental study was performed to improve the understanding of the characteristics of ultrasonic water atomization when excited with waves in the MHz range. In the present experiments, small volumes of water were atomized, observing the temporal evolution of the process. Typical diameters of the resulting droplets are of the order of a few microns. To visualize them, images were acquired with very high magnification. Appropriate lenses were used to enable high resolution at a distance from the flow. Droplet size distributions were also calculated with a Malvern diffractometer. Droplet exit velocity was measured using particle image velocimetry. It was noticeable that, as the remaining liquid mass deposited over the ultrasonic transducer decreased, the atomization characteristics changed, and a second peak of larger droplets appeared in the size distribution function. This phenomenon is related to the change in the curvature of the liquid surface. Although results are not conclusive, it appears that, under the conditions in this study, some observations about droplet formation are better described by cavitation phenomena rather than by the simplified surface wave theory usually invoked to explain these processes.     

Instantaneous and simultaneous planar velocity field measurements of two phases for turbulent mixing of high pressure sprays

This paper describes the tests of accuracy and the first application of a combined planar visualization technique. Its goal is two-phase flow discrimination, i.e. simultaneous measurements of velocity of droplets and ambient gas in the case of two-phase flow mixing, at the same location and with possible conditioning by “apparent diameter” (AD) of the droplets. It combines the mature techniques of particle image velocimetry (PIV), planar Mie scattering diffusion (PMSD), planar laser-induced fluorescence (PLIF), and it necessitates two synchronized cross-correlation systems, digital image treatment and analysis. This technique was developed with the objective of better describing the mixing between liquid and gaseous phases as in the case of high-pressure spray atomization in quiescent ambient gas. The basic principle of separation is to seed the ambient gas with micrometer particles and to tag the liquid with fluorescent dye. We use digital image treatment and analysis to discriminate between the phases. We use two cross-correlation PIV systems in order to obtain the velocity field of the droplets and gas simultaneously and separately at the same location. The digital image processing for separating the phases involves geometric measurement of droplet shapes. This leads to measurement of droplet parameters close to their real diameter, which could be used for analysis of actual mixing. A synchronized system composed of two CCD cameras is used for image recording, and two Nd:YAG lasers are used for generating pulsed light sheets at times t and t + δt. Tests were performed to check for different sources of errors. The combined technique was applied to measurements in high-pressure spray flow atomizing in a quiescent ambient gas, and first results are presented.     

Numerical investigation of unsteady shock waves in vapor-gas-droplet mixtures

The results of mathematical modeling of the evolution of unsteady shock waves in two-phase mixtures of inert gas, vapor and suspended liquid droplets with allowance for dynamic, thermal and mass phase interaction processes are presented. The influence of interphase mass transfer effects (droplet breakdown and evaporation, vapor condensation) on the structure of unsteady shock waves in vapor-gas-droplet mixtures is analyzed. The important influence of phase mass transfer and, in particular, droplet breakdown as a result of surface layer stripping by the gas flow on the distribution of the parameters of the carrier and dispersed components of the mixture behind the shock front is demonstrated. The effect of the principal governing parameters of the two-phase mixture on the unsteady shock wave propagation process is analyzed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 67–75, July–August, 1992.     

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.     

Scaling the effects of surface topography in the secondary atomization resulting from droplet/wall interactions

The impact of droplets onto micro-structured surfaces has been the focus of numerous recent studies, under the perspective of many different applications. However, much is still to be known about the effects of surface patterning in order to devise realistic physical models to accurately predict interfacial transfer rates. In this context, the present paper addresses the question of how to scale the effects of the surface topography to find adequate parameters, which can be easily obtained a priori. The approach is based on the characterization of the hydrodynamic and thermal behaviors of individual droplets impacting onto smooth and micro-structured heated surfaces, with the objective of quantifying the effects of the modified wettability associated with the topography of the surface. The focus is put on the thermal-induced mechanisms of secondary atomization as these are of particular interest for spray-cooling applications. The analysis suggests that different wetting properties lead to particular characteristics of the thermal-induced atomization, which can be related with the ratio between the roughness amplitude and the fundamental wavelength of the surface topography R _a/λ_R. This hypothesis is consistent with the theoretical prediction of the wetting behavior of the surfaces. The results also show a good correlation between the mean sizes of the secondary droplets generated by thermal-induced atomization and the ratio R _a/λ_R.     

Spray features in the near field of a flow-blurring injector investigated by high-speed visualization and time-resolved PIV

In a flow-blurring (FB) injector, atomizing air stagnates and bifurcates at the gap upstream of the injector orifice. A small portion of the air penetrates into the liquid supply line to create a turbulent two-phase flow. Pressure drop across the injector orifice causes air bubbles to expand and burst thereby disintegrating the surrounding liquid into a fine spray. In previous studies, we have demonstrated clean and stable combustion of alternative liquid fuels, such as biodiesel, straight vegetable oil and glycerol by using the FB injector without requiring fuel pre-processing or combustor hardware modification. In this study, high-speed visualization and time-resolved particle image velocimetry (PIV) techniques are employed to investigate the FB spray in the near field of the injector to delineate the underlying mechanisms of atomization. Experiments are performed using water as the liquid and air as the atomizing gas for air to liquid mass ratio of 2.0. Flow visualization at the injector exit focused on a field of view with physical dimensions of 2.3 mm × 1.4 mm at spatial resolution of 7.16 µm per pixel, exposure time of 1 µs, and image acquisition rate of 100 k frames per second. Image sequences illustrate mostly fine droplets indicating that the primary breakup by FB atomization likely occurs within the injector itself. A few larger droplets appearing mainly at the injector periphery undergo secondary breakup by Rayleigh–Taylor instabilities. Time-resolved PIV is applied to quantify the droplet dynamics in the injector near field. Plots of instantaneous, mean, and root-mean-square droplet velocities are presented to reveal the secondary breakup process. Results show that the secondary atomization to produce fine and stable spray is complete within a few diameters from the injector exit. These superior characteristics of the FB injector are attractive to achieve clean combustion of different fuels in practical systems.     

PTV investigation of phase interaction in dispersed liquid–liquid two-phase turbulent swirling flow

An investigation of dispersed liquid–liquid two-phase turbulent swirling flow in a horizontal pipe is conducted using a particle tracking velocimetry (PTV) technique and a shadow image technique (SIT). Silicone oil with a low specific gravity is used as immiscible droplets. A swirling motion is given to the main flow by an impeller installed in the pipe. Fluorescent tracer particles are applied to flow visualization. Red/green/blue components extracted from color images taken with a digital color CCD camera are used to simultaneously estimate the liquid and droplet velocity vectors. Under a relatively low swirl motion, a large number of droplets with low specific gravity tend to accumulate in the central region of the pipe. With increasing droplet volume fraction, the liquid turbulence intensity in the axial direction increases while that in the wall-normal direction decreases in the central region of the pipe. In addition, the turbulence modification in the present flow is strongly dependent on the droplet Reynolds number; however, the interaction of droplet-induced turbulences is significant due to vortex shedding, particularly at high droplet Reynolds numbers and higher droplet volume fraction.     

用格子Boltzmann方法模拟液滴撞击固壁动力学行为

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟.详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响.通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象.另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合.     

Large-Eddy Simulation of Interactions Between a Reacting Jet and Evaporating Droplets

Large-eddy simulation of a turbulent reactive jet with and without evaporating droplets is performed to investigate the interactions among turbulence, combustion, heat transfer and evaporation. A hybrid Eulerian–Lagrangian approach is used for the gas–liquid flow system. Arrhenius-type finite-rate chemistry is employed for the chemical reaction. To capture the highly local interactions, dynamic procedures are used for all the subgrid-scale models, except that the filtered reaction rate is modelled by a scale similarity model. Various representative cases with different initial droplet sizes (St ₀) and mass loading ratios (MLR) have been simulated, along with a case without droplets. It is found that compared with the bigger, slow responding droplets (St ₀ = 16), smaller droplets (St ₀ = 1) are more efficient in suppressing combustion due to their preferential concentration in the reaction zones. The peak temperature and intensity of temperature fluctuations are found to be reduced in all the droplet cases, to a varying extent depending on the droplet properties. Detailed analysis on the contributions of respective terms in a transport equation for grid-scale kinetic energy (GSKE) shows that the droplet evaporation effect on GSKE is small, while the droplet momentum effect depends on St ₀. When the MLR is sufficiently high, the bigger (St ₀ = 16) droplets can have profound influence on GSKE, and consequently on the formation and evolution of large-scale flow structures. On the other hand, the turbulence level is found to be lower in the droplet cases than in the pure flame case, due to the dissipative droplet dynamic effect.     

Secondary atomization of water and isooctane drops impinging on tilted heated surfaces

The present paper reports an experimental study aimed at characterizing the effects of heat transfer on the secondary atomization, which occurs during droplet impact on hot surfaces at conditions reproducing those occurring at fuel injection in internal combustion engines. The experiments consider single isooctane and water droplets impacting at different angles on a stainless steel surface with known roughness and encompass a range of Weber numbers from 240 to 600 and heat transfer regimes from the film-vaporization up to the Leidenfrost regime. The mechanisms of secondary breakup are inferred from the temporal evolution of the morphology of the impact imaged with a CCD camera, together with instantaneous measurements of droplet size and velocity. The combination of a technique for image processing with a phase Doppler instrument allows evaluating extended size distributions from 5.5 μm up to a few millimetres and to cover the full range of secondary droplet sizes observed at all heat transfer regimes and impaction angles. Temporal evolution of the size and velocity distributions are then determined. The experiments are reported at impact conditions at which disintegration does not occur at ambient temperature. So, any alteration observed in droplet impact behavior is thermally induced. The analysis is relevant for port fuel injection systems, where droplets injected to impact on the back surface of the valves, behave differently depending on fuel properties, particularly when the use of alcohols is considered, even as an additive to gasoline.     

Atomization of liquids by disintegrating thin liquid films using gas jets

Liquid atomization is useful in many applications, such as engineering, science, pharmaceutics, medicine, forensics and others. In the present research, an innovative methodology and a new device for atomization of liquids into mists of micron and submicron droplets have been developed. The new liquid-atomization method exploits the physical phenomenon of fragmentation of thin liquid films into fine micron and submicron droplets by gas jets. For several tested prototypes, the direct observations using a high-speed visualization technique have demonstrated that bubbles were generated within a liquid and their shells have been subsequently destroyed by applying a mechanical impulse (pressure of a compressed air) once the bubbles came over the liquid surface. The main characteristics of the generated tap water mists have been experimentally measured by means of the laser diffraction technique under various conditions for each prototype. One of the prototype devices allowed obtaining mists containing 90–99% of droplets smaller than 1 µm, with the minimum arithmetic and Sauter mean droplet diameters of 1.48 µm and 2.66 µm, and the 2.64 ml/min of droplet flow rate for 3.5 bar manometer pressure of atomizing air. The gas to liquid mass ratios (GLR) in the new device are depending on the atomizing tube length and the number of perforated orifices in the tube: more the tube length, hence more the number of perforated orifices, and therefore more liquid droplets will form for the same gas flow rate. The measured GLR values related to 1 m length of the utilized atomizing tube were in the range of 0.65–1.06, and for the specifically utilized atomizing tube of 72 mm length were among 9.07–14.67. The results of this study demonstrate that the developed method of generation of very fine droplet mists has many advantages over the existing techniques and can be perspective for many practical applications.     

竖直振动无黏液滴的法拉第不稳定性分析

由于外部周期性的振动而在液滴表面产生的Faraday不稳定效应广泛存在于超声雾化、喷涂加工等应用中,对Faraday不稳定性进行分析对研究振动液滴的表面动力学有着重要意义.本文将Faraday不稳定性问题从径向振动拓展到竖直振动,研究了竖直振动无黏液滴表面波的不稳定性.竖直方向的振动使得液滴动量方程为含有空间相关项和时...     

The analysis of silica suspensions atomization

The paper contains the results of experimental investigation of air–water and air–silica suspension atomization process in effervescent nozzles with internal mixing obtained by the use of the digital microphotography method. In experiments the different aqueous solutions of silica Aerosil 300 of different concentration have been used. The suspensions containing up to 0.04 (kg solid particles/kg solution) have Newtonian rheological properties. The observations were carried out at liquid flow rates changed from 0.0014 to 0.011 (kg/s) and gas flow rates from 0.00015 to 0.0065 (kg/s). It corresponded to gas to liquid mass ratios (GLR) values from 0.014 to 0.46. The analysis of photos shows that the droplets which have been formed during the liquid atomization have very different sizes. The differences between characteristics of effervescent atomization for water and suspensions used have not been observed. The present study confirmed the previous reports which suggested that the small particles added to solution do not change spray characteristics. The experimental results show that C_D and SMD are non-linear functions of GLR. Their values are decreasing rapidly as GLR is increased from zero to around 0.07 and thereafter decreasing at a slower rate with further increase in GLR. In the same point (GLR = 0.07) the value of α is maximal. The first regime is characteristic for bubbly flow. The second is typical of annular flow regime. Boundary between bubbly and annular flow regime is observed at GLR = 0.07 for investigated systems. The correlations for C_D and Sauter mean diameter were proposed. The results may be used for example to verify numerical models or comparisons with respect to similar atomization processes.