Breakup of a droplet at high velocity impacting a solid surface

We have studied the collision between a droplet of different liquids with high impact energy and a solid plate with varied surface roughness, which is characterized by a dimensionless Weber number (We, defined as the impact inertia of the droplet normalized by its surface force) extending up to 12,000 for water. To make such collision, a technique was developed to generate a single droplet with speed up to 42 m/s, which was initially driven by upstream air flow through a nozzle and accelerated to nearly the same velocity of the high-speed flow downstream. Via a high-speed photographing system, the various splashing mechanisms were investigated and a specific prompt splash on a smooth plate was found at sufficiently high We, which was different somehow from the conventionally defined one that was generally believed to occur only on a rough surface. The radius when multiple secondary droplets were shed out of the rim of the expanding lamella was found to scale almost invariantly with We at large values, whereas the coupled effect of liquid viscosity might affect the ultimate value.     

Droplet impact dynamics for two liquids impinging on anisotropic superhydrophobic surfaces

Droplet impingement experiments were performed on grooved hydrophobic surfaces with cavity fractions of 0, 80, and 93?% using droplets of water and a 50?%/50?% water/glycerol mixture. The influence of liquid viscosity, cavity fraction, and spreading direction, relative to the surface grooves, is explored qualitatively and quantitatively. The maximum droplet spread diameter, velocity of the rebounding jet, and the time delay between droplet impact and jet emission were characterized for Weber numbers, We, based on droplet impact speed and diameter, up to 500. The unequal shear stresses and contact angles influence the maximum spread diameters in the two primary spread directions. At We?>?100, the ratio of the spread diameter along the direction of the grooves to the spread diameter perpendicular to the grooves increases above unity with increasing We. The maximum droplet spread diameter is compared to recent predictive models, and the data reveal differing behavior for the two fluids considered. The results also reveal the existence of very high relative jet velocities in the range 5????We????15 for water droplets, while such jets were not observed for the more viscous mixture. Further, in the range 115????We????265, the water/glycerol jet formation dynamics are radically different from the water behavior. Most evident is the existence of two-pronged jets, which arise from the anisotropy of the surface and the unequal shear stresses and contact angles that prevail on the surfaces. It is these influences that give rise to differences in the maximum spread diameters in the two primary spread directions. Similar two-pronged jet emission was observed for water over the very narrow range of We from 91 to 96. The issuing jet velocities were also observed to increase with increasing cavity fraction for both fluids and over the entire range of We explored. Lastly, the elapsed time between droplet impact and jet emission decreased with increasing cavity fraction.     

Flow behaviour of negatively buoyant jets in immiscible ambient fluid

In this paper we investigate experimentally the injection of a negatively buoyant jet into a homogenous immiscible ambient fluid. Experiments are carried out by injecting a jet of dyed fresh water through a nozzle in the base of a cylindrical tank containing rapeseed oil. The fountain inlet flow rate and nozzle diameter were varied to cover a wide range of Richardson Ri (8 × 10⁻⁴ < Ri < 1.98), Reynolds Re (467 < Re < 5,928) and Weber We (2.40 < We < 308.56) numbers. Based on the Re, Ri and We values for the experiments, we have determined a regime map to define how these values may control the occurrence of the observed flow types. Whereas Ri plays a stronger role when determining the maximum penetration height, the effect of the Reynolds number is stronger predicting the flow behaviour for a specific nozzle diameter and injection velocity.     

Production of secondary drops during the single water drop impact onto a plane water surface

The normal impact of single water drops onto a plane water surface was studied experimentally to reveal the amount of secondary drops produced from the rim of crown-like interfacial structure. Within the experimental ranges tested, the ratio of the total mass of secondary drops to the mass of primary drop was approximately within 0–1 and correlated well as the function of dimensionless parameter K that consisted of the impact Weber number We and the Ohnesorge number Oh (K = We Oh ^−0.4). The dependences of the number and the mean diameter of secondary drops on K and dimensionless film thickness were also investigated.     

Investigation of droplets impinging on a deep pool: transition from coalescence to jetting

An experimental investigation of droplets impinging vertically on a deep liquid pool of the same fluid was conducted. Coalescence and jetting as two of the main regimes were identified and studied. Five fluids, distilled water, technical ethanol, n-pentane, methanol and 1-propanol were used for providing different liquid-phase physical properties with density from 600 to 1,000 kg/m³, viscosity from 0.20 to 2.00 mPa s, and surface tension from 13.7 to 72.0 mN/m. Except for the experimental run of n-pentane, which was carried out in n-pentane saturated vapor, the ambient gas for the other experiments was air. The impact processes of micro-level (diameter below 1 mm) droplets were captured using a high-speed camera with a backlight. The observations, velocity and diameter ranges of the experimental runs were described, and based on them, the effects of the liquid-phase properties were studied. It was found that both low viscosity and low surface tension can increase the instability during impact processes. By curve-fitting, the transition from coalescence to jetting was characterized by using two models, one employing the Weber number (We) and the Ohnesorge number (Oh), and one employing the Froude number (Fr) and the Capillary number (Ca). Both models characterize the coalescence-jetting threshold well. The We-Oh model was based on a commonly used model from Cossali et al. (in Exp Fluids 22:463–472, 1997) for characterizing coalescence-splashing. For the small droplet diameters (below 1 mm) considered in this study, it was required to modify the We-Oh model with a diameter-dependent term to fit the sharp change in thresholds for fluids with relatively high viscosity. The Fr-Ca model has not previously been presented in the literature. A comparison of the two models with literature data (Rodriguez and Mesler, J Colloid Interface Sci 106(2):347–352, 1985) indicates that they are also valid for impacts of droplets with diameters above 1mm. Calculation methods to generalize the two models were proposed.     

Measurement of liquid film thickness in a micro parallel channel with interferometer and laser focus displacement meter

In the present study, liquid film thicknesses in parallel channels with heights of H = 0.1, 0.3 and 0.5 mm are measured with two different optical methods, i.e., interferometer and laser focus displacement meter. Ethanol is used as a working fluid. Liquid film thicknesses obtained from two optical methods agree very well. At low capillary numbers, dimensionless liquid film thickness is in accordance with Taylor’s law. However, as capillary number increases, dimensionless liquid film thickness becomes larger than Taylor’s law for larger channel heights. It is attributed to the dominant inertial effect at high capillary numbers. Using channel height H for dimensionless liquid film thickness δ₀/H and hydraulic diameter D_h = 2H as the characteristic length for Reynolds and Weber numbers, liquid film thickness in a parallel channel can be predicted well by the circular tube correlation previously proposed by the authors. This is because curvature differences between bubble nose and flat film region are identical in circular tubes and parallel channels.     

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.     

Transient conjugated mixed-convective heat transfer in a vertical plate channel with one wall heated discretely

 This study presents a numerical solution of the unsteady conjugated mixed-convection heat transfer in a vertical plate channel with one wall suddenly subjected to either isoflux or isothermal discrete heat sources. The effects of the dimensionless heat source length H ₁, the dimensionless spacing between heat sources H ₂, the dimensionless channel length L, the dimensionless heated-plate thickness B _l, the wall-to-fluid conductivity ratio K and the ratio of Grashof number to Reynolds number Gr/Re on the interface heat flux, Nusselt number and bulk fluid temperature are discussed in detail. Results show that the discrete heating can cause the heat transfer direction conversely from the fluid to the heated plate during the transient period, which is more significant for the cases with larger L and H ₂. For the system with isoflux discrete heat sources, the time required to reach the steady-state is shorter for larger H ₂. While the trend is reverse for system with isothermal discrete heat sources. Additionally, a higher ratio of the input energy is axially conducted through the plate wall from heated sections to unheated regions for a larger H ₂ and B _l or smaller L. Received on 9 November 1998     

Investigation of binary drop rebound and coalescence in liquids using dual-field PIV technique

Experiments on binary drop collisions within an index-matched liquid were conducted for Weber numbers (We) in the range of 1–50. Drop pairs of water/glycerin mixture were injected horizontally into silicone oil and, due to gravitational effects, travelled on downward trajectories before colliding. A dual-field high-speed PIV measurement system was employed to quantify drop trajectories and overall collision conditions while simultaneously examining detailed velocity fields at the collision interface. Sequences of velocity and vorticity fields were computed for both larger and smaller fields of view. In the We range examined, both rebounding and coalescing behavior occurred. Coalescence was found to result from a combination of vortical flow within drops and strong drop deformation characteristic of higher We. Flow through the centers of opposing ring vortices, strengthened by drop deformation, enhanced drainage of the thin film in the impact region, leading to film rupture and coalescence. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Impact dynamics of drops on thin films of viscoelastic wormlike micelle solutions

The impact dynamics of water drops on thin films of viscoelastic wormlike micelle solutions is experimentally studied using a high-speed digital video camera at frame rates up to 4000 frame/s. The composition and thickness of the thin film is modified to investigate the effect of fluid rheology on the evolution of crown growth, the formation of satellite droplets and the formation of the Worthington jet. The experiments are performed using a series of wormlike micelle solutions composed of a surfactant, cetyltrimethylammonium bromide (CTAB), and a salt, sodium salicylate (NaSal), in deionized water. The linear viscoelastic shear rheology of the wormlike micelle solutions is well described by a Maxwell model with a single relaxation time while the steady shear rheology is found to shear thin quite heavily. In transient homogeneous uniaxial extension, the wormlike micelle solutions demonstrate significant strain hardening. The size and velocity of the impacting drop is varied to study the relative importance of Weber, Ohnesorge, and Deborah numbers on the impact dynamics. The addition of elasticity to the thin film fluid is found to suppress the crown growth and the formation of satellite drops with the largest effects observed at small film thicknesses. A new form of the splashing threshold is postulated which accounts for the effects of viscoelasticity and collapses the satellite droplet data onto a single master curve dependent only on dimensionless film thickness and the underlying surface roughness. Additionally, a plateau is observed in the growth of the maximum height of the Worthington jet height with increasing impact velocity. It is postulated that the complex behavior of the Worthington jet growth is the result of a dissipative mechanism stemming from the scission of wormlike micelles.     

A Simplified Model for Calculating Hydrodynamic Lubrication Film Thickness in Elastoplastic Line Contacts

The present paper proposes a simplified model for calculating hydrodynamic lubrication film thickness in elastoplastic line contacts. According to the Saint-Venant’s principle, the pressure in the contact is taken as uniformly distributed, this gives the contact surface elastic deformations in the inlet zone far away from the contact center close to real ones while gives those close to the contact center greater than real ones. This treatment is validated for hydrodynamic lubricated elastic contacts for relatively light loads and high rolling speeds. It gives the film thickness at the contact center a little higher than that calculated based on the real elastic model. The treatment is extended to a hydrodynamic lubricated elastoplastic line contact. The contact surfaces in the inlet zone are assumed as elastic and their deformations are calculated based on the uniform pressure distribution in the elastoplastic contact area. An inlet zone analysis is taken for obtaining the calculating equation of the hydrodynamic film thickness at the contact center. The equation overestimates the central film thickness but gives a satisfactory film thickness prediction for the heavy load which gives significant plastic deformations in the elastoplastic contact. It is found that when the load is lighter than 0.6 w _pc , the contact can be taken as elastic when calculating the central film thickness, while when the load is heavier than 0.6 w _pc , the contact can be taken as fully plastic; Here w _pc is the critical load for the contact fully plastic deformation. The plastic deformation in an elastoplastic line contact is found to reduce the hydrodynamic lubrication film thickness in the contact. This reduction is greater for higher rolling speeds and heavier loads. However, it is significantly dropped with increasing surface hardness.     

Effect of Droplet Size and Atomization on Spray Formation: A Priori Study Using Large-Eddy Simulation

The paper is mainly focused to the vast number of researchers who work within direct injection (DI) engine fuel spray simulations. The most common simulation framework today within the community is the Reynolds Averaged Navier Stokes (RANS) approach together with the Lagrangian Particle Tracking (LPT) method. In fact, this study is one of the first studies where high resolution LES/LPT diesel spray modeling is considered. The potential of LES to deepen the present day multidimensional LPT fuel spray simulations is discussed. Spray evolution is studied far from an injector by modeling a spray as a particle laden jet (PLJ). The effect of d on mixing in non-atomizing and atomizing sprays is thoroughly investigated at jet inlet Reynolds number Re?=?10⁴ and Mach number Ma?=?0.3. Based on and justified by rather recent and also quite old ideas, novel and compact views on droplet breakup in turbulent flows are pointed out from the literature. We use LES/LPT to illustrate that even in a low Weber number flow (We?<?13) the droplet breakup modeling may need considerable attention in contrast to what is typically assumed in the present-day breakup models. LES and LPT techniques are first applied to essentially confirm certain expected droplet size effects on spray shape in non-atomizing monodisperse sprays. In the simulations LES e.g. produces an expected turbulent dispersion pattern that depends on droplet diameter (d) without a droplet dispersion model in contrast to RANS. A new compact droplet breakup model is formulated and tested for droplets that break with a natural resonance time rate according to the Poisson process. As a result of the study: 1) the analysis gives a rigorous and enriching proof of currently existing views on droplet size effects on mixing, and 2) the presented a priori analysis points out the importance of modeling the resonance breakup even at a low We.     

Characterisation of gas-liquid two-phase flow in minichannels with co-flowing fluid injection inside the channel,part I: unified mapping of flow regimes

Present research highlights the potential of apparatuses with integrated minichannel packings to intensify gas-liquid-solid contacting. Especially an operation of these devices within the Taylor flow regime gained extraordinary attention due to its excellent heat and mass transfer and the segmented flow characteristics. However, criteria for flow regime transitions are mainly developed from water-similar fluids and are contradictory which hinders uniform flow regime prediction.This work presents a systematic analysis of adiabatic gas-liquid downflow in a square minichannel of 1.0 mm hydraulic diameter. In the mixing zone located within the flow channel, gas was injected into the co-flowing liquid by so-called capillary injectors with variable inner diameter (0.184, 0.317, 0.490 mm). Experiments were conducted using water, water-glycerol, and water-ethanol mixtures to cover a broad range of material properties. The gas and liquid superficial velocities were varied between 9.81·10^-4…2.72 m/s and 1.7·10⁻⁴…0.80 m/s, respectively. Taylor flow, Taylor-annular flow, annular flow, churn flow, and bubbly flow were observed. Using the Pi-theorem, 8 significant dimensionless groups dictating the flow transition were identified, namely u_{G, s}/u_{L, s}, Re_G, Re_L, We_G, We_L, Θ*, d_{In, CI}/d_h, and d_{Ou, CI}/d_h. Based on more than 1500 experimental data, criteria for the regime transitions of Taylor flow are provided. The derived flow regime map shows good agreement for all applied liquids and for the two larger injector geometries.     

High-resolution imaging of turbulence structures in jet flames and non-reacting jets with laser Rayleigh scattering

A comparative study of the length scales and morphology of dissipation fields in turbulent jet flames and non-reacting jets provides a quantitative analysis of the effects of heat release on the fine-scale structure of turbulent mixing. Planar laser Rayleigh scattering is used for highly resolved measurements of the thermal and scalar dissipation in the near fields of CH₄/H₂/N₂ jet flames (Re _d  = 15,200 and 22,800) and non-reacting propane jets (Re _d  = 7,200–21,700), respectively. Heat release increases the dissipation cutoff length scales in the reaction zone of the flames such that they are significantly larger than the cutoff scales of non-reacting jets with comparable jet exit Reynolds numbers. Fine-scale anisotropy is enhanced in the reaction zone. At x/d = 10, the peaks of the dissipation angle PDFs in the Re _d  = 15,200 and 22,800 jet flames exceed those of non-reacting jets with corresponding jet exit Reynolds numbers by factors of 2.3 and 1.8, respectively. Heat release significantly reduces the dissipation layer curvature in the reaction zone and in the low-temperature periphery of the jet flames. These results suggest that the reaction zone shields the outer regions of the jet flame from the highly turbulent flow closer to the jet axis.     

Study of airborne particles produced by normal impact of millimetric droplets onto a liquid film

The aim of this work is to carry out an experimental investigation into the generation of airborne microparticles when millimetric droplets of aqueous solutions impact onto a liquid film. Impact experiments using 3.9 mm diameter droplets were carried out for Weber numbers between 159 and 808, with a fixed Ohnesorge number of 2 × 10⁻³ and film parameters S _f (the ratio between the thickness of the liquid film h _film and the diameter of the impacting droplet d _i) between 0.3 and 1. Observed results show that the deposition/splashing threshold is independent of the parameter S _f in agreement with the data in the literature. The aerosol measurement results demonstrate the production of solid particles from the evaporation of secondary microdroplets with diameters less than 30 μm formed when splash occurs. The median diameter of these microdroplets is around 20 μm, corresponding to a value of d ₅₀/d _i = 5 × 10⁻³. Taken together, the results show that the mass and the number of particles emitted increase as the Weber number increases. Moreover, at a Weber number of 808, the results show that the mass and number of particles emitted increases as the parameter S _f decreases. In this case, the mean number of microdroplets emitted per impact is equal to 14 for S _f = 1 and equal to 76 for S _f = 0.3.     

Experimental study on surface wave and film breakdown of falling liquid film flow

In order to evaluate characteristics of the liquid film flow and their influences on heat and mass transfer, measurements of the instantaneous film thickness using a capacitance method and observation of film breakdown are performed. Experimental results are reported in the paper. Experiments are carried out at Re = 250–10000, T _in = 20–50°C and three axial positions of vertically falling liquid films for film thickness measurements. Instantaneous surface waveshapes are given by the interpretation of the test data using the cubic spline method. The correlation of the mean film thickness versus the film Reynolds number is also given by fitting the test data. It is revealed that the surface wave has nonlinear behavior. Observation of film breakdown is performed at Re = 1.40 × 10³–1.75 × 10⁴ and T _in = 85–95°C. From experimental results, the correlation of the film breakdown criterion can be obtained as follows: Bd = 1.567 × 10⁻⁶ Re ^1.183