Drag coefficients of spherical liquid droplets Part 2: Turbulent gaseous fields

The drag of non-evaporating, spherical, liquid droplets was measured in turbulent flow fields at parametric ranges relevant to spray combustion, characterized by the droplet Reynolds number, and the intensity and spatial scales of turbulence. The experimental apparatus comprised a wind-tunnel and a piezo-electric droplet generator. The procedure was to inject water droplets of uniform size co-currently and continuously with vertical turbulent air flows while droplet velocity was measured at different elevations using laser-Doppler velocimetry. Turbulence was characterized using hot-wire anemometry prior to droplet injection. Drag coefficients were calculated using these main measurements and the law of conservation of mechanical energy. Reynolds numbers were investigated in the range 10–100, in terms of the equivalent spherical diameter of a droplet, and the mean relative speed between the ambient gaseous field and the droplets. Weber numbers were much less than unity so droplets were effectively spherical. Relative intensities of turbulence were investigated in the range 20–65 percent, in terms of the mean relative speed. Spatial scales of turbulence were large in comparison to the droplets; the ratio between the spatial integral scale and the droplet diameter was in the range 11–38, and the Kolmogorov scale was comparable in size or smaller than the droplet diameter. Experimental data showed that the drag in turbulent fields under these conditions is not significantly different than that of solid spheres in a quiescent field at the same Reynolds number.The financial support of the Natural Sciences and Engineering Research Council of Canada and the Manufacturing Research Corporation of Ontario is gratefully acknowledged.     

Evaporation of bi-component droplets in a heated, highly turbulent flow

This work aims to understand the phenomena that occur in a combustion chamber where multi-component fuel droplets are injected. Many evaporation models exist but the influence of turbulence on spray vaporization is not yet well understood. This study gives a useful database to improve these models. The objective of the work is to measure the dispersion and the evaporation of bi-component (octane/3-pentanone) droplets and the resulting vapor mixing in a well-known, heated, highly turbulent channel flow. The carrier flow shows high turbulence levels, flat profiles for the mean velocity and the velocity fluctuations. The injected droplets have a large variety of behaviors due to the large polydispersion and to the turbulence. The evolution of 3-pentanone liquid concentration, mass flux, and droplet clusters are described. Mean concentration, fluctuations of concentration, and mixing of the vapor phase are characterized.     

Modeling and Simulation of Effects of Turbulence on Vaporization, Mixing and Combustion of Liquid-Fuel Sprays

The objective of this work is twofold. Firstly, the effects of turbulence intensity variations on the turbulent droplet dispersion, vaporization and mixing for non-reacting sprays (with and without swirl) are pointed out. Secondly, the effects of the coupling of the turbulence modulation with external parameters, such as swirl intensity, on turbulent spray combustion are analyzed in configurations of engineering importance. This is achieved by using advanced models for turbulence, evaporation and turbulence modulation implemented into FASTEST-LAG3D-codes: (1) To highlight the influence of turbulence modulation on some spray properties, a thermodynamically consistent modulation model has been considered besides the standard assumption and the well known Crowe's model. For turbulent droplet dispersion, we rely on the Markov-sequence formulation. (2) In order to characterize phase transition processes ongoing on droplets surfaces, a non-equilibrium evaporation model shows better agreement with experiments in comparison with the quasi-equilibrium-based evaporation models often used. (3) The results of turbulence intensity variations reveal the existence of a limited range out of which the increase or decrease of the turbulence intensity affects no more the efficiency of the heat and mass transfer. A derived characteristic number, a vaporization Damkhöler number, possesses a critical value which separates two different behavior regimes with respect to the turbulence/droplet vaporization interactions. (4) Under reacting conditions, it is shown how the evaporation characteristics, mixing rate and combustion process are strongly influenced by swirl intensity and turbulence modulation. In particular, the turbulence modulation modifies the evaporation rate, which in turn influences the mixing and the species concentration distribution. In the case under investigation, it is demonstrated that this effect cannot be neglected for low swirl intensities (Sw.Nu. ≤ 1) in the region far from the nozzle, and close to the nozzle for high swirl number intensities. In providing these particular characteristics, a reliable control of the mixing of gaseous fuel and air in evaporating and reacting sprays, and a possible optimization of the mixing process can tentatively be achieved.     

An experimental investigation of droplet evaporation and coalescence in a simple jet flow

An experimental study of a simple jet flow, which contains a dispersion of fine droplets, has been carried out in order to investigate the effect of turbulence, evaporation and coalescence on the droplet size distributions within the jet. Very little evaporation occurs in the potential core of the jet, while in the far-field, where the potential core has vanished and the droplets disperse more readily, evaporation occurs predominantly in the outer portions of the spray. Evidently, turbulence enhances the evaporation rate of droplets at the edges of the spray, and fresh air entrained from the outer regions increases the evaporative driving force. Coalescence has also been observed within the spray, although this effect is rather subtle compared to the evaporation effect in the dilute jets investigated here. Nevertheless, sufficient measurements have been taken to validate, at least partially, any coalescence models, in addition to any turbulence and evaporation models for dilute poly-disperse sprays.     

Numerical modeling of electrostatic spray painting transfer processes in rotary bell cup for automotive painting

Here we implemented a 3D comprehensive Eulerian-Lagrangian model in order to investigate the electrostatic spray transfer processes in the high-speed rotary bell sprayer. This efficient algorithm contains spray dynamics, airflow, paint droplets tracking and an electrostatic effect to simulate atomization. The algorithm is implemented using the OpenFOAM package. A solver for the particle trajectory was used to illustrate the process of spray transport and also the interaction of the airflow and the particle that is solved by momentum coupling. Creating an initial condition of the particle approach has been proposed that is matched with practical applications. The fluid-dynamics is simulated by solving the unsteady 3D compressible Navier-Stokes equations. Unsteady flow is computed by using a Large eddy simulation (LES) turbulence approach, while the motion of the particles is simulated by tracking the droplet size distribution approach. The model correctly predicts that the bell cup spin forces the paint particles to fall off from the bell surface towards the high-velocity airflow. The present work illustrates a tentative benchmark and contains a systematic analysis of the recirculation zone length, the toroidal vortex, the overspray phenomena and the flowfield characteristics like mean velocity, pressure, turbulent kinetic energy and velocity fluctuation. The results indicate as dominant operating parameter the air-paint flow rate with voltage level deeply affecting the spray shape. A more uniform distribution of the coating is obtained by growing this high-velocity shaping airflow, although the values of the transfer efficiency (TE) are reduced. The distribution of the particle size is very sensitive to changes in the rotational speed. Experimental results obtained in this study put forward a clear link between the shaping air flow rate and the rotation frequency on the aerodynamics and also provide valuable insights to design modern ERBS. The paint spray distribution obtained in the present work is validated against coating experimental results with suitable accuracy.     

A two-equation turbulence model for jet flows laden with vaporizing droplets

A two-equation turbulence model for steady incompressible two-phase flows including phase change has been recently developed by Mostafa & Elghobashi (1984). This model is tested for the flow of a turbulent axisymmetric gaseous jet laden with evaporating liquid droplets. To avoid the problem of density fluctuations of the carrier phase at this stage, only isothermal flow is considered and vaporization is assumed to be due to the vapor concentration gradient. The continuous size distribution of the droplets is approximated by finite size groups. Each group is considered as a continuous phase interpenetrating and interacting with the carrier phase. Two test cases have been predicted by the model. The first is for a Freon-11 spray issuing from a round nozzle, where experimental data are available at distances equal to or greater than 170 nozzle diameters. Good agreement between the data and the predictions was achieved. The second is for a methanol spray where no experiments are available yet and the predictions consider the flow region close to the nozzle ($\text{z/D < 40}$). The results of the methanol spray include distributions of the mean velocity, volume fractions of the different phases, concentration of the evaporated material in the carrier phase, turbulence intensity and shear stress of the carrier phase, droplet diameter distribution, and the jet spreading rate. In this case the results are analyzed based on a qualitative comparison with the corresponding single phase jet flow.     

Numerical analysis of swirling turbulent droplet-laden flow and heat transfer in a sudden pipe expansion

The effect of swirling intensity on the structure and heat transfer of a turbulent gas–droplet flow after a sudden pipe expansion has been numerically simulated. Air is used as the carrier phase, and water, ethanol, and acetone are used as the dispersed phase. The Eulerian approach is applied to simulate the dynamics and heat transfer in the dispersed phase. The gas phase is described by a system of Reynolds-averaged Navier-Stokes (RANS) equations, taking into account the effect of droplets on mean transport and turbulent characteristics in the carrier phase. Gas phase turbulence is predicted using the second-moment closure. A swirling droplet-laden flow is characterized by an increase in the number of small particles on the pipe axis due to their accumulation in the zone of flow recirculation and the action of the turbulent migration (turbophoresis) force. A rapid dispersion of fine droplets over the pipe cross-section is observed without swirling. With an increase in swirling intensity, a significant reduction in the length of the separation region occurs. The swirling of a two-phase flow with liquid droplets leads to an increase in the level of turbulence for all three types of liquid droplets investigated in this work due to their intensive evaporation. It is shown that the addition of droplets leads to a significant increase in heat transfer in comparison with a single-phase swirling flow. The greatest effect of flow swirling on heat transfer intensification in a two-phase gas-droplet flow is obtained for the droplets of ethanol and water and smallest effect is for the acetone droplets.     

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.     

Evaporation of a droplet larger than the Kolmogorov length scale immersed in a relative mean flow

An experimental effort to understand the contribution of turbulence to the evaporation rate of fuel droplets has been performed with particular attention to conditions when the turbulence scale is smaller than the droplet diameter. N-heptane has been chosen as working fluid to give measurable evaporation rates from droplet images over relatively short experiment times. An active turbulence grid wind tunnel is built for the requirements of this experiment. A camera triggered by a pulse generator takes images of the droplets pinned on wires across the tunnel. The results show a small increase in evaporation rate with increasing turbulence intensity, and that mean flow around the droplets has more impact on evaporation than does the turbulence state.     

A detailed experimental investigation of well-defined,turbulent evaporating spray jets of acetone

This paper aims at investigating the detailed structure of turbulent non-reacting dilute spray flows using advanced laser diagnostics. A simple spray jet nozzle is designed to produce a two-phase slender shear flow in a co-flowing air stream with well-defined boundary conditions. The carrier flow is made intentionally simple and easy to model so that the focus can be placed on the important aspects of droplet dispersion and evaporation, as well as turbulence–droplet interactions. Phase Doppler interferometry is employed to record droplet quantities, while planar laser-induced fluorescence imaging is applied separately to obtain acetone vapour data. Measurements are conducted for four acetone spray jets in air at several axial stations starting from the nozzle exit. The combined liquid and vapour mass fluxes of acetone integrated across the jet at downstream locations agree satisfactorily with the total mass flow rate of acetone injected.     

Droplet dynamics and size characterization of high-velocity airblast atomization

Airblast atomizers are especially useful and commonplace in liquid fuel combustion applications. However, the spray formation processes, the droplet dynamics and the final drop size distributions are still not sufficiently understood due to the coupled gas-liquid interactions and turbulence generation. Therefore, empirical and semi-empirical approaches are typically used to estimate the global spray parameters. To develop a physical understanding of the spray evolution, a plain-jet airblast atomizer was investigated in an atmospheric spray rig using the phase-Doppler technique. The simultaneous drop size and axial and radial velocity components were measured on radial traverses across the spray at various axial distances from the nozzle for a range of atomizing pressures. The droplet turbulent and mean kinetic energies were found to be proportional to the atomizing pressure. Hence, the scatter of the radial motion of the droplets increased with the atomizing pressure. A droplet stability analysis was performed to locate the regions characterized by ongoing secondary atomization. The volume-to-surface diameter, D₃₂, of the fully developed spray was compared with estimates provided by five published formulae. The role of liquid viscosity, hence the Ohnesorge number, was found to be negligible in the investigated regime. Three commonly used size distribution functions were fitted to the measured data to analyze their dependence on the atomizing pressure. The Gamma distribution function was found to give the best approximation to the atomization process.     

Fuel injection model for Euler–Euler and Euler–Lagrange large-eddy simulations of an evaporating spray inside an aeronautical combustor

Large-Eddy Simulations (LES) of an evaporating two-phase flow in an experimental burner are investigated. Two different numerical approaches for the simulation of the dispersed phase are coupled to the same gaseous solver: a mesoscopic Eulerian method and a Lagrangian particle tracking technique. The spray is represented by a single droplet size owing to the locally monodisperse formulation of the employed mesoscopic Eulerian approach. Both approaches use the same drag and evaporation models. They do not take into account the atomization process and a simplified injection model is applied instead. The presented methodology, referred as FIM-UR (Fuel Injection Method by Upstream Reconstruction) defines injection profiles for the monodisperse spray produced by a pressure-swirl atomizer. It is designed so as to ensure similar spray characteristics for both approaches and allows for a direct comparison between them. After a validation of the purely gaseous flow in the burner, liquid-phase dynamics and droplet dispersion are qualitatively and quantitatively evaluated for the Eulerian and Lagrangian simulations. Results obtained for both approaches are in very good agreement and compare reasonably with experiments, indicating that simplified injection methods are appropriate for the simulation of realistic combustor geometries.     

Modeling the evolution of droplet size distribution in two-phase flows

A theoretical model is developed in the present study to simulate droplet motion and the evolution of droplet size distribution (DSD) in two-phase air/dispersed water spray flows. The model takes into account several processes which influence DSD and droplet trajectory: droplet collision and coalescence, evaporation and cooling, gravitational settling, and turbulent dispersion of dispersed phase. The DSDs determined by the model at different locations in a two-phase flow are evaluated by comparing them to experimental observations obtained in an icing wind tunnel. The satisfactory coincidence between simulation and experimental results proves that the model is reliable when modeling two-phase flows under icing conditions. The model is applied for two particular examples in which the modification of DSD is calculated in two-phase flows under conditions describing in-cloud icing and freezing drizzle.     

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.     

Flow characteristics of a large-size pressure-atomized spray using DTV

The purpose of this study is to characterize the atomization of a jet of water sprayed into the air at high velocity through a commercial nozzle widely used for sprinkler irrigation. The typical diameter of the droplets present in the spray is in the range of several tens of micrometers to several millimeters. They are visualized by ombroscopy. A specific Droplet Tracking Velocimetry (DTV) technique is developed to estimate the size and velocity of these highly polydispersed droplets that are distinctly non spherical. This analysis is performed from the rupture of the liquid core region (about a distance of 550 nozzle diameters) to the dispersed zone (about a distance of 900 nozzle diameters). With this technique, we obtain joint size-velocity measurements that are rarely produced. Especially two velocity components and also a large diameter range are characterized at the same time; while with other techniques, such as Particle Doppler Anemometry (PDA), the diameter range is quite reduced and requires specific settings. Additional measurements of the liquid volume fraction are performed using a single mode fiber-optic probe. In the light of our experimental data, it appears that the turbulent droplet motion in the spray is strongly anisotropic. This anisotropy is quite unexpected because other studies on sprays (generally concerned with engine applications) show a relatively low anisotropy. We attribute this increase of anisotropy to the fact that, for this type of spray, the droplet relaxation time is long in comparison to the characteristic time of the turbulence and that biggest droplets are still submitted to atomization process. This strong anisotropy is responsible for the poor radial dispersion of the spray.     

The Effect of Droplets Evaporation on Turbulence Modification and Heat Transfer Enhancement in a Two-Phase Mist Flow Downstream of a Pipe Sudden Expansion

Turbulent droplet-laden flow downstream of a sudden pipe expansion is numerically studied using an Eulerian two-fluid model. The model is used to investigate the effect of droplet evaporation on the particle dispersion and on the gas phase turbulence modification. Turbulence suppression in the case of evaporating droplets is hardly observed near the wall, and the level of turbulence tends to the corresponding value for the single-phase flow regime. In the flow core, where evaporation is insignificant, a decrease in the level of gas turbulence (to 20 % as compared to a single-phase flow) can be observed. The maximal effect of droplet evaporation is obtained in the wall region of the tube. A considerable increase in the maximal value of heat exchange on adding the evaporating droplets to the separated flow is shown (more than 1.5-fold as compared to the single-phase flow at a small value of droplet mass concentration of M _L1≤ 0.05). The addition of the solid non-evaporating particles causes a slight increase in the maximum value of heat transfer in the case of small particles and a decrease in heat transfer in the case of large particles.     

Acetone Droplet Behavior in Reacting and Non Reacting Turbulent Flow

Acetone droplet characteristics in reacting and non-reacting turbulent flow are predicted and compared to experimental data. Investigations are conducted to study the effects of surrounding environment properties on the velocities, dispersion, and evaporation of a relatively volatile spray fuel that featured a wide range of Stokes numbers. The simulations are performed in the framework of Reynolds Averaged Navier Stokes equations along with the Eulerian-Lagrangian approach in which 12 different classes of the dispersed phase. The phase transition is modeled by the Langmuir-Knudsen law that accounts for non equilibrium effects based on a consistent determination of the molar mass fraction on the droplet surfaces. For the droplet dispersion, the Markov sequence model is improved by adding a correction drift term to the fluid fluctuation velocity at the parcel position along the droplet trajectory. This correction term aimed at accounting for the non-homogeneity effects in the turbulent flow. The combustion is captured using the Bray-Moss-Libby model that is extended to account for the partially premixed spray combustion. The chemistry is described with the flamelet model using a recent detailed reaction mechanism that involves 84 species and 409 reactions for which the Lewis number is not set to the unity. Mean droplet velocities for reacting and non-reacting test cases are compared with experimental data. Good agreement is observed. The spray is interacting with the nozzle edge developing new classes and relatively dense region. Hence the RMS-velocities close to the nozzle exit plan demonstrate discrepancies. The droplets group combustion effect is found to be important in the modeling of the burning velocity which influences the flame propagation. Reasonable agreements between the numerical and the experimental results are also observed in the spray flux and temperature profiles.     

Turbulent time and length scale measurements in high-velocity open channel flows

In high-velocity open channel flows, the measurements of air–water flow properties are complicated by the strong interactions between the flow turbulence and the entrained air. In the present study, an advanced signal processing of traditional single- and dual-tip conductivity probe signals is developed to provide further details on the air–water turbulent level, time and length scales. The technique is applied to turbulent open channel flows on a stepped chute conducted in a large-size facility with flow Reynolds numbers ranging from 3.8E+5 to 7.1E+5. The air water flow properties presented some basic characteristics that were qualitatively and quantitatively similar to previous skimming flow studies. Some self-similar relationships were observed systematically at both macroscopic and microscopic levels. These included the distributions of void fraction, bubble count rate, interfacial velocity and turbulence level at a macroscopic scale, and the auto- and cross-correlation functions at the microscopic level. New correlation analyses yielded a characterisation of the large eddies advecting the bubbles. Basic results included the integral turbulent length and time scales. The turbulent length scales characterised some measure of the size of large vortical structures advecting air bubbles in the skimming flows, and the data were closely related to the characteristic air–water depth Y ₉₀. In the spray region, present results highlighted the existence of an upper spray region for C > 0.95–0.97 in which the distributions of droplet chord sizes and integral advection scales presented some marked differences with the rest of the flow.     

Rotation effects on a fully-developed turbulent pipe flow

A fully-developed turbulent pipe flow is allowed to pass through a rotating pipe section, whose axis of rotation coincides with the pipe axis. At the exit end of the rotating section, the flow passes into a stationary pipe. As a result of the relaxation of surface rotation, the turbulent flow near the pipe wall is affected by extra turbulence production created by the large circumferential shear strain set up by the rapid decrease of the rotational velocity to zero at the wall. However, the flow in the most part of the pipe is absent of this extra turbulence production because the circumferential strain is zero as a result of the solid-body rotation imparted to the flow by the rotating pipe section. The combined effect of these two phenomena on the flow is investigated in detail using hot-wire anemometry techniques. Both mean and turbulence fields are measured, together with the wall shear and the turbulent burst behavior at the wall. A number of experiments at different rotational speeds are carried out. Therefore, the effects of rotation on the behavior of wall shear, turbulent burst at the wall, turbulence production and the near-wall flow can be documented and analysed in detail.