A technique for measuring Lagrangian and Eulerian particle statistics in a turbulent flow

An experimental technique is described which has been developed to study particle dispersion in a round turbulent jet. Droplets are injected on the jet axis, and a laser sheet and position sensitive photomultiplier tube are used to track their radial displacement. Data processing is greatly simplified compared to video or photo imaging techniques which provide similar measurements. Statistically large samples are used to calculate dispersion and axial velocity as a function of axial downstream distance or particle time-of-flight. Dispersion and velocity statistics can be computed which are Lagrangian or Eulerian in nature. The technique has been demonstrated with 69 m droplets of hexadecane in a jet of air with a Reynolds number of 15,000; in principle it could be used to study the motion of very small, quasi-fluid particles.     

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.     

Flame Structure and Propagation in Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Analysis

Three-dimensional Direct Numerical Simulations (DNS) in canonical configuration have been employed to study the combustion of mono-disperse droplet-mist under turbulent flow conditions. A parametric study has been performed for a range of values of droplet equivalence ratio ?_d, droplet diameter a_d and root-mean-square value of turbulent velocity u^′. The fuel is supplied entirely in liquid phase such that the evaporation of the droplets gives rise to gaseous fuel which then facilitates flame propagation into the droplet-mist. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d>1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion: the premixed mode ocurring mainly under fuel-lean conditions and the non-premixed mode under stoichiometric or fuel-rich conditions. The prevalence of premixed combustion was seen to decrease with increasing droplet size. Furthermore, droplet-fuelled turbulent flames have been found to be thicker than the corresponding turbulent stoichiometric premixed flames and this thickening increases with increasing droplet diameter. The flame thickening in droplet cases has been explained in terms of normal strain rate induced by fluid motion and due to flame normal propagation arising from different components of displacement speed. The statistical behaviours of the effective normal strain rate and flame stretching have been analysed in detail and detailed physical explanations have been provided for the observed behaviour. It has been found that the droplet cases show higher probability of finding positive effective normal strain rate (i.e. combined contribution of fluid motion and flame propagation), and negative values of stretch rate than in the stoichiometric premixed flame under similar flow conditions, which are responsible for higher flame thickness and smaller flame area generation in droplet cases.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.     

Hydrodynamics of a drying multicomponent liquid droplet

Using microscopy methods on light and dark fields, the flow patterns developing in drying droplets of pure transparent liquids, solutions, and suspensions of micro- and nanoparticles are investigated. The flow patterns inside drying droplets of real, colloid, and mixed solutions containing nanoand microparticles-markers are studied by means of video and photo registration of microscopic images. The analysis of particle displacements indicates the existence of a global convective flow which forms a toroidal circulation with an ascending jet at the droplet center. The typical types of the structures depending on the droplet composition are distinguished. It is shown that the intensity of the flow inside the droplet affects the surface convection. The effect of the hydrodynamic flow on the transport of a substance, forming the dry-deposit texture, is studied.     

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.     

Heat transfer in droplet-laden turbulent channel flow with phase transition in the presence of a thin film of water

We present results of a numerical study of turbulent droplet-laden channel flow with phase transition. Previous studies of the same system did not take into account the presence of gravity. Here, we do so introducing a thin film of water at the bottom wall and permitting droplets to fall into and merge with it. We treat the carrier phase with the Eulerian approach. Each droplet is considered separately in the Lagrangian formulation, adopting the point–particle approximation. We maintain the film thickness constant by draining water from the bottom wall to compensate for (a) the droplets that fall onto the film and (b) evaporation/condensation. We also maintain on average the total mass of water in the channel by inserting new droplets at the top wall to compensate for the water that has been drained from the bottom wall. We analyze the behavior of the statistically averaged gas and droplet quantities focusing on the heat exchange between the two phases. We increase (a) the initial droplet diameter keeping the same initial droplet volume fraction and (b) the initial number of droplets in the channel keeping their diameter the same. In both parameter studies we find that droplets grow less than in the reference case. In case (a) this is explained by the larger velocity with which they travel to the bottom wall and in case (b) by the lower rate of condensation of vapor due to the presence of neighboring droplets.     

Drop trajectory predictions and their importance in the design of spray dryers

A knowledge of the trajectories of atomized droplets in both the nozzle zone (where the droplets are rapidly decelerating from their initial high velocity) and in the free-entrainment zone (where the droplets are conveyed by the drying gas) is required for the design of spray dryers, since it governs the evaporative capacity and thermal efficiency of the chamber, while affecting the moisture content and general quality of the product through the control of the drying time.

The trajectories of droplets in three-dimensional motion were determined theoretically in both zones. In the case of two-fluid pneumatic atomizers, the characteristics of the jet of atomizing fluid were found to be important in both the zones.

Predictions of droplet trajectories were tested in an experimental circular cocurrent spray-drying chamber with a conical bottom, in which the drying air was introduced tangentially near the top. Water was used as the feed material. A study was made of the effects of liquid feed rate and temperature, drying air flowrate and temperature, and of nozzle position on the thermal efficiency and evaporative capacity of the chamber. The results were interpreted in the light of the droplet trajectories predicted.     

Breakup of Newtonian and non-Newtonian fluids in air jets

The breakup of droplets of non-Newtonian fluids has been investigated by high speed photography and impaction following preliminary results of Newtonian fluids, which confirmed the suitability of the measurement techniques. Single droplets with diameters from 2.4 to 3.3 mm, were arranged to fall under gravity into a jet of air with velocities up to 36o m/s. The droplets of Newtonian fluids, water and Diesel oil, were atomised in the expected manner within three main regimes characterised by the Weber number of the droplet and air jet conditions, while similar droplets of non-Newtonian fluids were found not to atomise but to develop under shear and stretching into ligaments of fluid separated from a local region of their surface; these ligaments were elongated until breakup occurred, though not into small droplets as with the Newtonian fluids. Some of the non-Newtonian fluids (TEP with 7.5% and lo% K125, with and without water) were found not to break up at the maximum speed of the tests and they will be re-examined at higher jet velocities. Increase in the concentration of K125 in TEP resulted in higher critical speed for a given droplet diameter.The authors would like to thank Prof. J. H. Whitelaw of Imperial College for many discussions and useful suggestions during the course of this work, and Dr. G. Cambray of CBDE for his valuable administrative support     

Interaction of a polydisperse spray with vortices

The objective of the present work is to provide, through the association of optical diagnostics on a well-chosen experimental configuration, new insights into the coupling of a vortical gaseous flow with a polydisperse evaporating spray representative of practical injections. A cloud of droplets is injected in an inert laminar round jet, axisymmetric and pulsated, enabling the study of the interaction of strong-vorticity structures with a polydisperse spray. The experiment is a laboratory-scale representation of realistic injection configurations such as in engine combustion chambers or industrial burners. The chosen set-up leads to a well-controlled configuration and allows the coupling of two optical diagnostics, particle imaging velocimetry (PIV) and interferometric particle imaging (IPI), which leads to the study of both the flow dynamic and the droplet size distribution. The behaviour of droplets is analysed regarding their relaxing and evaporating properties. Size-conditioned preferential concentration of both weakly evaporating and strongly evaporating sprays is investigated. Droplet trajectories are also analysed by means of high-rate tomographic visualizations. The time history between their ejection from the nozzle and their interaction with the vortex is strongly related to the droplet preferential concentration and the observed heterogeneous repartition in the gas flow.     

Large eddy simulation and experiment of shear breakup in liquid-liquid jet: Formation of ligaments and droplets

Understanding the shear breakup in jet flows and the formation of droplets from ligaments is important to determine the final droplet size distribution (DSD). The initial droplet size, which affects the final DSD, is considered to be generated by the shear breakup. Large eddy simulation (LES) was performed to investigate the shear breakup in liquid-liquid jet flows. The explicit Volume of Fluid (VOF) model with the geometric reconstruction scheme was used to capture the oil-water interface. The estimated oil distribution including wave peaks, ligaments, droplets and water streaks were compared to the experiments with a good agreement. The estimated DSD matched with the measurements favorably well. In the simulation, the formation of droplets with a smooth and curved surface from ligaments or sheet-like structures was obtained. Different mechanisms were observed along with the shear layer including the formation of droplets from ligament through the capillary forces, breakage of a droplet into smaller ones and attachment of a droplet to a ligament. The destructive shear forces and resisting surface tension forces were quantified on stretching and retracting ligaments. The influence of internal viscous force was found to be negligible due to low oil viscosity. The critical capillary number was found to be larger than 5.0 for ligaments breaking with the shear breakup. The capillary number was below unity for retracting ligaments. The coalescence of two equal-sized droplets was obtained in the shear breakup region. The shear stress magnitude at the contact region increased more than two folds. The total surface area decreased nearly 20% after the coalescence.     

Microgravity experiments and numerical studies on ethanol/air spray flames

Spray flames are known to exhibit amazing features in comparison with single-phase flames. The weightless situation offers the conditions in which the spray characteristics can be well controlled before and during combustion. The article reports on a joint experimental/numerical work that concerns ethanol/air spray flames observed in a spherical chamber using the condensation technique of expansion cooling (based on the Wilson cloud chamber principle), under microgravity.We describe the experimental set-up and give details on the creation of a homogeneous and nearly monosized aerosol. Different optical diagnostics are employed successfully to measure the relevant parameters of two-phase combustion. A classical shadowgraphy system is used to track the flame speed propagation and allow us to observe the flame front instability. The complete characterization of the aerosol is performed with a laser diffraction particle size analyser by measuring the droplet diameter and the droplet density number, just before ignition. A laser tomography device allows us to measure the temporal evolution of the droplet displacement during flame propagation, as well as to identify the presence of droplets in the burnt gases. The numerical modelling is briefly recalled. In particular, spray-flame propagation is schematized by the combustion spread in a 2-D lattice of fuel droplets surrounded by an initial gaseous mixture of fuel vapour and air.In its spherical expansion, the spray flame presents a corrugated front pattern, while the equivalent single-phase flame does not. From a numerical point of view, the same phenomena of wrinkles are also observed in the simulations. The front pattern pointed out by the numerical approach is identified as of Darrieus–Landau (DL) type. The droplets are found to trigger the instability. Then, we quantitatively compare experimental data with numerical predictions on spray-flame speed. The experimental results show that the spray-flame speed is of the same order of magnitude as that of the single-phase premixed flame. On the other hand, the numerical results exhibit the role played by the droplet radius in spray-flame propagation, and retrieve the experiments only when the droplets are small enough and when the Darrieus–Landau instability is triggered. A final discussion is developed to interpret the various patterns experimentally observed for the spray-flame front.     

DROPLET COLLISION AND COALESCENCE MODEL

A new droplet collision and coalescence model was presented,a quick-sort method for locating collision partners was also devised and based on theoretical and experimental results,further advancement was made to the droplet collision outcome. The advantages of the two implementations of smoothed particle hydrodynamics (SPH) method were used to limit the collision of droplets to a given number of nearest droplets and define the probability of coalescence,numerical simulations were carried out for model validation.Results show that the model presented is mesh-independent and less time consuming,it can not only maintains the system momentum conservation perfectly,but not susceptible to initial droplet size distribution as well.     

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.     

Consistency issues of Lagrangian particle tracking applied to a spray jet in crossflow

Numerical simulations are performed for multiphase jets in crossflow. The flow solver uses an Eulerian/Lagrangian approach. Turbulence in the gas phase is modeled in the framework of large eddy simulation. The dispersed phase is handled using Lagrangian particle tracking. The model assumptions of solvers for Lagrangian particle tracking are critically assessed for typical flow conditions of spray jets in crossflow. The droplets are assumed to be spherical and isolated. It is shown that several model assumptions are apparently inconsistent in larger portions of the flow field. Firstly, average Weber numbers can be so large that the model assumption to regard droplets as spherical is questionable, not only near the nozzle, but also in the far-field. Secondly, the average droplet spacing can be so low that droplets directly interact with each other, again also in the far-field. Thirdly, the average Stokes numbers in the jet region can be so large that the phase coupling between the dispersed and continuous phase is weak. Some remedies to these deficiencies are proposed.     

An experimental and numerical study into turbulent condensing steam jets in air

 Temperatures, velocities, and droplet sizes are measured in turbulent condensing steam jets produced by a facial sauna, for varying nozzle diameters and varying initial velocities (Re=3,600–9,200). The release of latent heat due to droplet condensation causes the temperature in the two-phase jet to be significantly higher than in a single-phase jet. At some distance from the nozzle, droplets reach a maximum size and start to evaporate again, which results in a change in sign of latent heat release. The distance of maximum size is determined from droplet size measurements. The experimental results are compared with semi-analytical expressions and with a fully coupled numerical model of the turbulent condensing steam jet. The increase in centreline temperature due to droplet condensation is successfully predicted. Received: 5 April 2000 / Accepted: 15 November 2000     

超声速气流中雾化燃料喷射的三维数值研究

首次用双流体模型对雾化燃料在扩张形超燃室中沿九喷嘴顺流喷射的混合问题进行了数值研究。气相用迎风 TVD格式求解三维全 Navier- Stokes方程 ,液相用预估、校正 NND格式求解三维 Euler方程。相间相互作用的常微分方程用预估、校正Runge- Kutta法求解。用三维 Poisson方程生成网格。结果表明 :气相较液相的扩散效果好 ,小直径液滴的扩散效果好。相间速度滑移、改变气相喷射压力和喷射速度对液相扩散的贡献不大 ,但调整喷射角度会明显地增强液相的扩散、混合 ,本文结果未出现阻塞。     

Experiments examining drag in linear droplet packets

This paper presents an experimental study of vertically traveling droplet packets, where the droplets in each packet are aligned linearly, one behind another. The paper describes in detail, an experimental apparatus that produces repeatable, linearly aligned, and isolated droplet packets containing 1–6 droplets per packet. The apparatus is suitable for examining aerodynamic interactions between droplets within each packet. This paper demonstrates the performance of the apparatus by examining the drag reduction and collision of droplets traveling in the wake of a lead droplet. Comparison of a calculated single droplet trajectory with the detailed droplet position versus time data for a droplet packet provides the average drag reduction experienced by the trailing droplets due to the aerodynamic wake of the lead droplet. For the conditions of our experiment (4 droplet packet, 145 m methanol droplets, 10 m/s initial velocity, initial droplet spacing of 5.2 droplet diameters, Reynolds number approx. 80) the average drag on the first trailing droplet was found to be 75% of the drag on the lead droplet.     

The use of single-point measurements to characterize dynamic behavior in sprays

In an effort to derive dynamic information from a single-particle counting device such as the phase-Doppler interferometer (PDI), a one-dimensional motion model and a sequential sorting algorithm have been developed to calculate the distance and velocity difference between consecutive droplets in sprays, as well as the number of droplets in close proximity (grouplets) and the frequency of grouplet formation. By focusing on sprays that are primarily uni-directional, this model can be used to highlight dynamic behavior and provide information that may indicate tendencies for clustering, collision, coalescence or evaporation among the droplets – information that is not available when considering the histograms of droplet size and velocity normally produced by PDIs. The motion model and sorting algorithm are validated for a mono-sized droplet stream generated by piezo-electric excitation, and are then applied to two other sprays to determine dynamic behavior. For a pressurized-liquid nozzle emitting water, the expected trends of droplet deceleration and spatial dispersion are found moving away from the exit of the nozzle, and significant in velocity differences between closely-spaced droplet pairs are observed. No dominant frequencies are present in the grouping of the droplets. For a droplet-laden air jet excited at a known frequency, similar results are calculated for droplet spacing, grouping and velocity differences, and the driving frequency is also obtained from analysis of the sequentially-sorted PDI data. Although the models developed assume an idealized PDI system, real operating characteristics of PDI systems place limits on the accuracy of these analyses. Received: 10 October 1998/Accepted: 5 October 1999