Numerical prediction of breakup mode of contracting gas filament in liquid

The breaking up of gas filament in liquid is important in many industrial and scientific applications. In this study, a transient axisymmetric model with the level set method is built up to examine the dynamics of a contracting gas filament, and to determine the effects of the aspect ratio, Ohnesorge (Oh) number, and viscosity ratio on its breakup mode. The filament undergoes no break, middle break, or end-pinching modes with increasing aspect ratio at either a low or a high Oh number, and one critical initial aspect ratio is observed for each Oh number. The fate of the filament is determined by the interaction of capillary waves on its surface, and can be predicted accurately by using the one-dimensional wave superposition method. The decreasing viscosity ratio of liquid over gas reduces the critical initial aspect ratio for the fate transition between the no break and breakup modes, and this effect is reduced at a low viscosity ratio. These findings may be helpful in fabricating gas bubbles and their breakup suppression.     

Filament-breaking length—a measure of elasticity in extension

A technique has been developed to assess elasticity in extension of fluids which are weakly elastic. The technique is based on stretching a fluid sample held between two small closely spaced co-axial disks until the fluid bridge or filament breaks. The distance between the disks on breakup, ‘the breaking length’, is measured. When the fluid is elastic, this length is greater than that of an equivalent inelastic fluid. An inelastic baseline was established by measuring the breaking lengths of a wide variety of Newtonian fluids. Measurements with weakly elastic fluids reveal that the extra breaking length increases with polymer type and concentration and with the rate of stretching, the expected behavior for elastic liquids. The breaking lengths of model paper coatings were measured and found to correlate with the degree of misting in a roll coating machine, indicating that droplet formation is related to extensional elasticity.     

Droplets splashing upon films of the same fluid of various depths

We explore the effects of fluid films of variable depths on droplets impacting into them. Corresponding to a range of fluid “film” depths, a non-dimensional parameter—H*, defined as the ratio of the film thickness to the droplet diameter—is varied in the range 0.1≤H*≤10. In general, the effect of the fluid film imposes a dramatic difference on the dynamics of the droplet–surface interaction when compared to a similar impact on a dry surface. This is illustrated by the size distribution and number of the splash products. While thin fluid films (H*≈0.1) promote splashing, thicker films (1≤H*≤10) act to inhibit it. The relative roles of surface tension and viscosity are investigated by comparison of a matrix of fluids with low and high values of these properties. Impingement conditions, as characterized by Reynolds and Weber numbers, are varied by velocity over a range from 1.34 to 4.22 m/s, maintaining a constant droplet diameter of 2.0 mm. The dependence of splashing dynamics, characterized by splash product size and number, on the fluid surface tension and viscosity and film thickness are discussed.     

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.     

Low-Reynolds-number droplet motion in a square microfluidic channel

In this study, we investigate computationally the low-Reynolds-number droplet motion in a square micro-channel, a problem frequently encountered in microfluidic devices, enhanced oil recovery and coating processes. The droplet deformation and motion are determined via a three-dimensional spectral boundary element method for wall-bounded flows. The effects of the flow rate, viscosity ratio and droplet size on the interfacial dynamics are identified for droplets smaller and larger than the channel size and for a wide range of viscosity ratio. Owing to the stronger hydrodynamic forces in the thin lubrication film between the droplet interface and the solid walls, large droplets exhibit larger deformation and smaller velocity. Under the same average velocity, a droplet in a channel shows a significantly smaller deformation and higher velocity than in a cylindrical tube with the same size, owing to the existence of the corners’ area in the channel which permits flow of the surrounding fluid. A suitable periodic boundary implementation for our spectral element method is developed to study the dynamics of an array of identical droplets moving in the channel. In this case, the droplet deformation and velocity are reduced as their separation decreases; the reduction is influenced by the flow rate, viscosity ratio and more significantly the droplet size.     

Producing droplets smaller than the nozzle diameter by using a pneumatic drop-on-demand droplet generator

A pneumatic droplet generator to produce water/glycerin droplets smaller than the nozzle diameter is described. The generator consists of a T-junction with a nozzle fit into one opening, the second opening connected to a gas cylinder through a solenoid valve and the third connected to a length of steel tubing. The droplet generator is filled with liquid. Opening the valve for a preset time creates a pulse of alternating negative and positive pressure in the gas above the surface of the liquid, ejecting a single droplet through the nozzle. Droplet formation was photographed and the pressure variation in the droplet generator recorded. The effect of various experimental parameters, such as nozzle size, pressure pulse width and liquid properties on droplet formation was investigated. Small droplets could not be generated when liquid viscosity was too low or too high. For pure water, droplet diameters were several times that of the nozzle. Using more viscous glycerin mixtures, droplets with diameters as small as 65% of the nozzle diameter could be produced.     

Droplet evaporation on porous and non-porous ceramic solids heated from top

 This paper presents the results of an experimental investigation, of the effect of radiation heat, on the evaporation of five droplet sizes of pure water, softly deposited on porous and non-porous ceramic solids, at temperature ranging from 75 to 250 °C. Both solids were instrumented with several surface and in-depth thermocouples, and had the same thermal properties. Results show that, the droplet evaporation time, and the surface recovery time for the porous solid were shorter than that of non-porous solid for the same droplet size under identical conditions. Also, smaller droplets were more efficient for cooling both solids. The results were compared with data for the evaporation of water droplets on similar ceramic solids heated from bottom (Abu-Zaid M; Atreya A (1994) J Heat Transfer 116: 694–701). The comparison shows that, the heat radiation has a significant effect of reducing evaporation time, recovery time, and droplet volume of influence for both solids, at the same initial surface temperature. Received on 6 December 1999 / Published online: 29 November 2001     

Spreading of inkjet droplet of non-Newtonian fluid on solid surface with controlled contact angle at low Weber and Reynolds numbers

The dynamics of inkjet droplet of non-Newtonian fluid on glass substrates was investigated experimentally and compared with that of Newtonian fluid. The non-Newtonian fluids used here were 100 ppm solutions of polyethylene oxide (300k, 600k and 900k) dissolved in the 1:1 mixture of water and glycerin. Weber number (We) was 2–35 and Ohnesorge number was fixed at 0.057 ± 0.003. The wettability of solid substrate was also varied. The diameter of inkjet droplets in the present study was about 50 μm and was much smaller than the size of the previous studies on drop impact. Due to the development of a thin and long thread at the rear of the main drop the jetting window of polymer solution was much narrower than that of Newtonian fluid, and hence the experimental range of Weber number was restricted. The impact scenarios of non-Newtonian inkjet droplets were found to be qualitatively different from those of Newtonian droplets during the receding phase while they were almost the same as the Newtonian fluid case during the kinematic phase. The spreading diameter at the equilibrium was well correlated with the modified Weber number (We′ = We/(1 − cos θ_eq)) as in the case of Newtonian fluid, where θ_eq is the equilibrium contact angle. The similarity or disparity between the Newtonian and non-Newtonian cases was discussed considering the conformation of polymer chains during each stage of drop deformation.     

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.     

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     

Attenuation of a shock wave passing through a cloud of water droplets

The mitigation of a planar shock wave caused by a cloud of calibrated water droplets was studied both experimentally and numerically. Experiments were carried out, with different shock wave Mach numbers ranging from 1.1 to 1.8, in a vertical shock tube coupled with a droplet generator which produced a well-characterized cloud of droplets of 120, 250 and 500 μm in diameter. By exploiting such an experimental set-up, we successfully measured the attenuation of a normal shock wave when passing through the water droplet cloud. This series of experiments allowed to identify the main parameters of this investigation and a clear dependence between the attenuation of the shock wave and terms governing the regimes of droplet breakup has been found. On the other hand, to support this experimental approach, 1D unsteady calculations were performed in similar configurations. Although the mathematical model based on an Eulerian/Eulerian approach was actually incomplete, the first comparisons between the experiments and the simulations were rather interesting and pointed out the need to improve the physical model, by taking into account the fragmentation and the vaporization of the droplets submitted to the shock wave as well as the size distribution of the water spray.     

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.     

Spray characteristics of a gasoline direct injector injector with short durations of injection

The sprays into atmosphere from a GDI injector were visualised and the velocity and droplet characteristics measured at an injection pressure of 50 bar and at different injection durations, with emphasis on short injection periods. The images show that the initial and closing delay times were 0.225 and 0.2 ms, respectively, and that the cone angle increased with injection duration to a constant value of 62° at 0.5 ms. They also revealed large droplets ahead of the main spray with its smaller droplets. An injection duration of 0.15 ms led to fuel leaving the injector with little atomisation, but a 30% increase led to the formation of the cone, which was present for times greater than 0.5 ms. The poor atomisation associated with short injection durations and the initial phase of longer injections, was due to low swirl velocities. The droplet velocities were higher in the initial phase of injection than in the main phase, with values up to 50 m/s. The Sauter mean diameters of the initial and main-spray droplets were approximately 55 and 35 μm respectively and with a tendency to decrease with time from the start of injection. Received: 20 October 2000/Accepted: 30 March 2001     

液滴对爆炸冲击波的衰减作用

为分析液滴对舰船舱内爆炸冲击波的耗散与衰减作用,通过有限元分析方法,建立冲击波作用于不同尺寸单个液滴和多排液滴的模型,分析冲击波与单个及多个液滴的作用过程及液滴形态变化,对冲击波衰减规律进行分析总结。得到结论如下:单个液滴模型中,小液滴破碎更迅速,破碎的规律性强;大液滴抛撒现象发生较早,抛撒出的小液滴数目多,但整体变化规律性偏差;不同尺寸单个液滴对冲击波有一定的衰减作用,衰减率随液滴尺寸增大而增大,线性规律较明显;成排液滴对冲击波有明显的衰减作用,相同液滴密度下衰减率随着液滴数量的增多而增大,呈现明显的线性特性。     

Droplet streams for serial crystallography of proteins

Serial diffraction of proteins requires an injection method to deliver analyte molecules—preferably uncharged, fully hydrated, spatially oriented, and with high flux—into a focused probe beam of electrons or X-rays that is only a few tens of microns in diameter. This work examines conventional Rayleigh sources and electrospray-assisted Rayleigh sources as to their suitability for this task. A comparison is made and conclusions drawn on the basis of time-resolved optical images of the droplet streams produced by these sources. Straight-line periodic streams of monodisperse droplets were generated with both sources, achieving droplet diameters of 4 and 1 micrometer, respectively, for the conventional and electrospray-assisted versions. Shrinkage of droplets by evaporation is discussed and quantified. It is shown experimentally that proteins pass undamaged through a conventional Rayleigh droplet source.     

Investigation of droplet collisions for solutions with different solids content

The collision behaviour of droplets and the collision outcome are investigated for high viscous polymer solutions. For that purpose, two droplet chains produced by piezoelectric droplet generators are directed towards each other at a certain angle so that individual droplet pairs collide. For recording the collision event, one double-image and one high-speed CCD camera were used. One camera is positioned perpendicular to the collision plane recording the outcome of the collision, and the second camera is aligned parallel to the collision plane to assure that the droplet chains are exactly in one plane. A new approach for tracking droplets in combination with an extended particle tracking velocimetry algorithm has been developed. Time-resolved series of pictures were used to analyse the dynamics of droplet collisions. The three different water soluble substances were saccharose and 1-Ethenyl-2-pyrrolidone (PVP) with different molecular weights (K17, K30). The solvent was demineralised water. The solids contents ranged from 20 to 60 %, 5 to 25 % and 5 to 35 %, yielding dynamic viscosities in the range of 2–60 mPa s. Results were collected for different pairs of impact angles and Weber numbers in order to establish common collision maps for characterising the outcomes. Here, relative velocities between 0.5 and 4 m/s and impact parameters in the interval from 0 to 1 for equal-sized droplets (Δ = 1) have been investigated. Additionally, satellite formation will be discussed exemplarily for K30. A comparison with common models of different authors (Ashgriz and Poo in J Fluid Mech 221:183–204, 1990; Estrade et al. in Int J Heat Fluid Flow 20:486–491, 1999) mainly derived for low viscous droplets revealed that the upper limit of their validity is given by an Ohnesorge number of Oh = 0.115 and a capillary number of Ca = 0.577. For higher values of these non-dimensional parameters and hence higher dynamic viscosities, these models are unable to predict correctly the boundaries between collision scenarios. The model proposed by Jiang et al. (J Fluid Mech 234:171–190, 1992), which includes viscous dissipation, is able to predict the boundary between coalescence and stretching separation for higher viscosities (i.e. Oh > 0.115 and Ca > 0.577). However, the model constants are not identical for different solution properties. As a conclusion, an alteration of the collision appearance takes place because of the relative importance between surface tension and viscosity.     

Water in water emulsions: phase separation and rheology of biopolymer solutions

Partially miscible polymers in solution do not separate into two macroscopic phases; in general they behave as viscoelastic fluids containing droplets of the minority phase dispersed into a continuous majority phase (emulsion type systems). Both phases contain two types of polymers and solvent in variable amounts. With time, the smaller droplets tend to merge into larger ones and eventually sediment. Provided the time stability of the emulsion is long enough and the size of the droplets does not exceed a few tens of microns, the emulsion can be characterized by conventional rheological methods as an effective medium, both in the linear regime (viscoelasticity) and under flow. We investigated a ternary system composed by two biopolymers, a protein (caseinate) and a polysaccharide (alginate) in aqueous solution and established an analogy between these phase separated solutions and immiscible blends of polymers. We first characterized the biopolymers and determined the phase diagram at room temperature that we interpreted within the framework of the Edmond and Ogston model. For the rheological investigations, starting with an initial composition of the system, we separated the two phases by centrifugation. The individual phases were then characterized through their viscoelastic and flow behaviors. By recombining variable amounts of these phases, thereby varying only their volume fractions, we were able to prepare stable emulsions with constituents having constant compositions. The effect of shear on these emulsions was investigated. After different shearing protocols, the size of the droplets was derived from the Palierne model and the flow curves were analyzed. The droplet sizes were compared to the critical capillary numbers and coalescence predictions. The flow curves and the dynamic viscosities of the emulsions were interpreted with a model recently proposed by Kroy et al. that extends earlier work of Oldroyd (1953), Schowalter et al. (1968), and Frankel and Acrivos (1970). Received: 11 September 2000/Accepted: 21 December 2000     

Simultaneous planar measurement of droplet velocity and size with gas phase velocities in a spray by combined ILIDS and PIV techniques

A new approach for simultaneous planar measurement of droplet velocity and size with gas phase velocities is reported, which combines the out-of-focus imaging technique ‘Interferometric Laser Imaging Droplet Sizing’ (ILIDS) for planar simultaneous droplet size and velocity measurements with the in-focus technique ‘Particle Image Velocimetry’ (PIV) for gas velocity measurements in the vicinity of individual droplets. Discrimination between the gas phase seeding and the droplets is achieved in the PIV images by removing the glare points of focused droplet images, using the droplet position obtained through ILIDS processing. Combination of the two optical arrangements can result in a discrepancy in the location of the centre of a droplet, when imaging through ILIDS and PIV techniques, of up to about 1 mm, which may lead to erroneous identification of the glare points from droplets on the PIV images. The magnitude of the discrepancy is a function of position of the droplet’s image on the CCD array and the degree of defocus, but almost independent of droplet size. Specifically, it varies approximately linearly across the image along the direction corresponding to the direction of propagation of the laser sheet for a given defocus setting in ILIDS. The experimental finding is supported by a theoretical analysis, which was based on geometrical optics for a simple optical configuration that replicates the essential features of the optical system. The discrepancy in the location was measured using a monodisperse droplet generator, and this was subtracted from the droplet centres identified in the ILIDS images of a polydisperse spray without ‘seeding’ particles. This reduced the discrepancy between PIV and ILIDS droplet centres from about 1 mm to about 0.1 mm and hence increased the probability of finding the corresponding fringe patterns on the ILIDS image and glare points on the PIV image. In conclusion, it is shown that the proposed combined method can discriminate between droplets and ‘seeding’ particles and is capable of two-phase measurements in polydisperse sprays.     

The motion of long bubbles through viscoelastic fluids in capillary tubes

The penetration of long gas bubble through a viscoelastic fluid in a capillary tube has been studied in order to investigate the influence of viscoelastic material properties on the hydrodynamic coating thickness and local flow kinematics. Experiments are conducted for three tailored ideal elastic (Boger) fluids, designed to exhibit similar steady shear properties but substantially different elastic material functions. This allows for the isolation of elastic and extensional material effects on the bubble penetration process. The shear and extensional rheology of the fluid is characterized using rotational and filament stretching rheometers (FSR). The fluids are designed such that the steady-state extensional viscosity measured by the FSR at a Deborah number (De) greater than 1 differs over three orders of magnitude (Trouton ratio = 10³–10⁶). The experiment set up to measure the hydrodynamic coating thickness is designed to provide accurate data over a wide range of capillary numbers (0.01 < Ca < 100). The results indicate that the coating thickness in this process increases with an increase in the extensionally thickening nature of the fluid. Experiments are also conducted using several different capillary tube diameters (0.1 < D < 1 cm), in order to compare responses at similar Ca but different flow De. Suitable scaling methods and nonlinear viscoelastic constitutive equations are explored to characterize the displacement process for polymeric fluids. Bubble tip shapes at different De are recorded using a CCD camera, and measured using an edge detection algorithm. The influence of the mixed flow field on the bubble tip shape is examined. Particle tracking velocimetry experiments are conducted to compare the influence of viscoelastic properties on the velocity field in the vicinity of the bubble tip. Local shear and extension rates are calculated in the vicinity of the bubble tip from the velocity data. The results provide quantitative information on the influence of elastic and extensional properties on the bubble penetration process in gas-assisted injection molding. The bubble shape and velocity field information provides a basis for evaluating the performance of constitutive equations in mixed flow. Received: 19 January 1999 Accepted: 30 June 1999     

Deformation and orientation of single droplets during shear flow: combined effects of confinement and compatibilization

Combined effects of geometrical confinement and compatibilization on the deformation and orientation of single droplets during steady-state shear flow are investigated in a counter-rotating cell by means of microscopic observations. The model system consists of polydimethylsiloxane droplets of varying sizes and viscosities in a polyisobutylene matrix. To this system, a premade polyisobutylene–polydimethylsiloxane block copolymer is added as compatibilizer in different concentrations. For each droplet, the equilibrium interfacial tension is determined by comparing droplet axes with the predictions of the confined Minale model for uncompatibilized droplets at the appropriate degree of confinement. Although large reductions in interfacial tension are seen for all compatibilized droplets, it is shown that the effect of compatibilization on droplet deformation and orientation can efficiently be taken into account in the equilibrium capillary number. This way, for all viscosity ratios and confinement ratios, steady-state data for compatibilized and uncompatibilized droplets coincide, and agree well with the predictions of the confined Minale model at sub-critical conditions. For near-critical capillary numbers, compatibilization slightly reduces droplet deformation and postpones breakup, irrespective of the degree of confinement.