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.     

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.     

Collision dynamics of high-speed droplets upon layers of variable thickness

The collision dynamics between a droplet and a film has been studied with high-impact energy that can be grouped in a dimensionless Weber number, We, as normalized by surface energy. To accomplish this, we have developed a technique based on cutting of a high-speed jet, which can generate a single droplet with speed up to 23 m/s and We on the order of thousands. It was found that the boundaries indicating the occurrence of a central jet and that of a secondary droplet disintegrated from the jet decreased monotonically with increased dimensionless film thickness, H, and remained constant when the film thickness was larger than the crater depth. However, the transition designating multiple droplets that are originated from a central jet shows a non-monotonic trend with the variation of H, with a minimum We being at H ≈ 3, which is about the maximum crater depth, owing to a tuning behavior. The critical We for splashing that occurs at an early phase immediately after the impact is relatively sensitive to the film thickness only when H is between 1 and 2, which increases with reduced H. At large We (≳2,570 for high H), the ejected crown is closed to form a bubble and the transition boundary reveals a similar dependence on H as that for creation of a central jet.     

Influence of surfactant upon air entrainment hysteresis in curtain coating

The onset of air entrainment for curtain coating onto a pre-wetted substrate was studied experimentally in similar parameter regimes to commercial coating (Re = ρQ/μ = O(1), We = ρQ u _c /σ = O(10), Ca = μU/σ = O(1)). Impingement speed and viscosity were previously shown to be critical parameters in correlating air entrainment data with three qualitatively different regimes of hydrodynamic assist identified (Marston et al. in Exp Fluids 42(3):483–488, 2007a). The interaction of the impinging curtain with the pre-existing film also led to a significant hysteretic effect throughout the flow rate-substrate speed parameter space. For the first time, results considering the influence of surfactants are presented in attempt to elucidate the relative importance of surface tension in this inertia-dominated system. The results show quantitative and qualitative differences to previous results with much more complex hysteretic behaviour which has only been reported previously for rough surfaces.     

Determination of the aerodynamic droplet breakup boundaries based on a total force approach

The determination of the critical We_g number separating the different breakup regimes has been extensively studied in several experimental and numerical works, while empirical and semi-analytical approaches have been proposed to relate the critical We_g number with the Oh_l number. Nevertheless, under certain conditions, the Re_g number and the density ratio ε may become important. The present work provides a simple but reliable enough methodology to determine the critical We_g number as a function of the aforementioned parameters in an effort to fill this gap in knowledge. It considers the main forces acting on the droplet (aerodynamic, surface tension and viscous) and provides a general criterion for breakup to occur but also for the transition among the different breakup regimes. In this light, the present work proposes the introduction of a new set of parameters named as We_g,eff and Ca_l monitored in a new breakup plane. This plane provides a direct relation between gas inertia and liquid viscosity forces, while the secondary effects of Re_g number and density ratio have been embedded inside the effective We_g number (We_g,eff)     

Deformation and breakup of liquid drops in a gas stream at nearly critical Weber numbers

Experimentally determined values of the critical Weber number available from the literature are scattered over a very wide range ofWe _c from 2.2 to 99.6. To study one possible source of these discrepancies an experimental investigation was made of the deformation and breakup of water droplets at nearly critical Weber numbers. Experiments were conducted in a small horizontal wind tunnel. A continuous stream of uniform water droplets was allowed to fall perpendicularly to the continuous stream of air. The time histories of water droplets were recorded by using a high speed camera. Five different basic behaviours of water droplets were recorded in the range ofWe = 11 to 14. It was found that an increase in the Weber number in this region resulted in an increased percentage of droplets with regular bag type breakup.     

Characterization of a system generating a homogeneous isotropic turbulence field by free synthetic jets

A facility inspired by Hwang and Eaton (2004a, b) for generating a homogeneous isotropic turbulence was built, the objective being to study evaporating droplets in the presence of turbulence. Turbulence was produced by the mixing of six synthetic jets, in ambient atmosphere. Combined PIV and LDA techniques were used to measure the statistical turbulence properties. The turbulence produced was found to be homogeneous isotropic with a small mean flow within a domain having an average size of 50 mm × 50 mm × 50 mm. The rms fluctuations were of the order of 0.9 m/s, corresponding to a Taylor Reynolds number of 240 and an integral length scale of about 40 mm. This apparatus proved to be well suited to the study of the evaporation of droplets in a controlled turbulence field.     

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.     

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.     

A blood-mimicking fluid for particle image velocimetry with silicone vascular models

For accurate particle image velocimetry measurements in hemodynamics studies, it is important to use a fluid with a refractive index (n) matching that of the vascular models (phantoms) and ideally a dynamic viscosity matching human blood. In this work, a blood-mimicking fluid (BMF) composed of water, glycerol, and sodium iodide was formulated for a range of refractive indices to match most common silicone elastomers (n = 1.40–1.43) and with corresponding dynamic viscosity within the average cited range of healthy human blood (4.4 ± 0.5 cP). Both refractive index and viscosity were attained at room temperature (22.2 ± 0.2°C), which eliminates the need for a temperature-control system. An optimally matched BMF, suitable for use in a vascular phantom (n = 1.4140 ± 0.0008, Sylgard 184), was demonstrated with composition (by weight) of 47.38% water, 36.94% glycerol (44:56 glycerol–water ratio), and 15.68% sodium iodide salt, resulting in a dynamic viscosity of 4.31 ± 0.03 cP.     

A numerical study of non-isothermal turbulent coaxial jets

In this work, we propose to study non isothermal air–air coaxial jets with two different approaches: parabolic and elliptic approaches. The standard k−ε model and the RSM model were applied in this study. The numerical resolution of the equations governing this flow type was carried out for: the parabolic approach, by a “home-made” CFD code based on a finite difference method, and the elliptic approach by an industrial code (FLUENT) based on a finite volume method. In forced convection mode (Fr = ∞), the two turbulence models are valid for the prediction of the mean flow. But for turbulent sizes, k−ε model gives results closer to those achieved in experiments compared to RSM Model. Concerning the limit of validity of the parabolic and elliptic approaches, we showed that for velocities ratio r lower than 1, the results of the two approaches were satisfactory. On the other hand, for r > 1, the difference between the results became increasingly significant. In mixed convection mode (Fr ≅ 20), the results obtained by the two turbulence models for the mean axial velocity were very different even in the plume region. For the temperature and the turbulent sizes the two models give satisfactory results which agree well with the correlations suggested by the experimenters for X ≥ 20. Thus, the second order model with σ _t = 0.85 is more effective for a coaxial jet study in a mixed convection mode.     

Exact and numerical solutions of a fully developed generalized second-grade incompressible fluid with power-law temperature-dependent viscosity

We compute exact and numerical solutions of a fully developed flow of a generalized second-grade fluid, with power-law temperature-dependent viscosity (μ=θ^-M), down an inclined plane. Analytical solutions are found for the case when M=m+1, m≠1, m being a constant that models shear thinning (m<0) or shear thickening (m>0). The exact solutions are given in terms of Bessel functions. The numerical solutions indicate that both the velocity and the temperature increase with decreasing Froude number and that there is a critical value of Fr below which temperature “overshoots” its free surface value of unity. This phenomena is not reported in the work of Massoudi and Phuoc [Fully developed flow of a modified second grade fluid with temperature dependent viscosity, Acta Mech. 150 (2001) 23-37.] for viscosity that depends exponentially on temperature.     

The interchain pressure effect in shear rheology

Recently, the tube diameter relaxation time in the evolution equation of the molecular stress function (MSF) model (Wagner et al., J Rheol 49: 1317–1327, 2005) with the interchain pressure effect (Marrucci and Ianniruberto, Macromolecules 37:3934–3942, 2004) included was shown to be equal to three times the Rouse time in the limit of small chain stretch. From this result, an advanced version of the MSF model was proposed, allowing modeling of the transient and steady-state elongational viscosity data of monodisperse polystyrene melts without using any nonlinear parameter, i.e., solely based on the linear viscoelastic characterization of the melts (Wagner and Rolón-Garrido 2009a, b). In this work, the same approach is extended to model experimental data in shear flow. The shear viscosity of two polybutadiene solutions (Ravindranath and Wang, J Rheol 52(3):681–695, 2008), of four styrene-butadiene random copolymer melts (Boukany et al., J Rheol 53(3):617–629, 2009), and of four polyisoprene melts (Auhl et al., J Rheol 52(3):801–835, 2008) as well as the shear viscosity and the first and second normal stress differences of a polystyrene melt (Schweizer et al., J Rheol 48(6):1345–1363, 2004), are analyzed. The capability of the MSF model with the interchain pressure effect included in the evolution equation of the chain stretch to model shear rheology on the basis of linear viscoelastic data alone is confirmed.     

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.     

In situ observations of unusual drop deformation and wobbling in simple shear flow

Deformation and wobbling of a liquid drop immersed in a liquid matrix were studied under mild shear conditions for various viscosity ratios. In situ visualization experiments were conducted on a homemade transparent Couette cell incorporated to the Paar Physica MCR500 shear rheometer. The effect of drop or matrix elasticity was examined and was found to play a major role in both deformation and wobbling processes. Experimental results were compared to Jackson and Tucker (J Rheol 47:659–682, 2003), Maffettone and Minale (J Non-Newton Fluid Mech 78:227–241, 1998) and Yu and Bousmina (J Rheol 47:1011–1039, 2003) ellipsoidal models. It was found that the agreement between the Newtonian models and the experimental results required an increase in the drop viscosity. Such increment in viscosity was found to scale with the first normal stress difference.     

Experimental assessment of scale effects affecting two-phase flow properties in hydraulic jumps

A hydraulic jump is the rapid transition from a supercritical to subcritical free-surface flow. It is characterised by strong turbulence and air bubble entrainment. New air–water flow properties were measured in hydraulic jumps with partially developed inflow conditions. The data set together with the earlier data of Chanson (Air bubble entrainment in hydraulic jumps. Similitude and scale effects, 119 p, 2006) yielded similar experiments conducted with identical inflow Froude numbers Fr ₁ = 5 and 8.5, but Reynolds numbers between 24,000 and 98,000. The comparative results showed some drastic scale effects in the smaller hydraulic jumps in terms of void fraction, bubble count rate and bubble chord time distributions. The present comparative analysis demonstrated quantitatively that dynamic similarity of two-phase flows in hydraulic jumps cannot be achieved with a Froude similitude. In experimental facilities with Reynolds numbers up to 10⁵, some viscous scale effects were observed in terms of the rate of entrained air and air–water interfacial area.     

Group classification of one-dimensional nonisentropic equations of fluids with internal inertia

A systematic application of the group analysis method for modeling fluids with internal inertia is presented. The equations studied include models such as the nonlinear one-velocity model of a bubbly fluid (with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski Zhurnal Prikladnoj Mekhaniki i Tekhnitheskoj Fiziki 3, 102–111, 1960; Kogarko Dokl. AS USSR 137, 1331–1333, 1961; Wijngaarden J. Fluid Mech. 33, 465–474, 1968), and the dispersive shallow water model (Green and Naghdi J. Fluid Mech. 78, 237–246, 1976; Salmon 1988). These models are obtained for special types of the potential function W(r,[(r)\dot],S){W(\rho,\dot \rho,S)} (Gavrilyuk and Teshukov Continuum Mech. Thermodyn. 13, 365–382, 2001). The main feature of the present paper is the study of the potential functions with W _S  ≠ 0. The group classification separates these models into 73 different classes.     

Computations of breakup modes in laminar compound liquid jets in a coflowing fluid

We present a numerical investigation of breakup modes of an axisymmetric, laminar compound jet of immiscible fluids, which flows in a coflowing immiscible outer fluid. We use a front-tracking/finite difference method to track the unsteady evolution and breakup of the compound jet, which is governed by the Navier–Stokes equations for incompressible Newtonian fluids. Numerical results show that depending on parameters such as the Reynolds number Re (in the range of 5–30) and Weber Number We (in the range of 0.1–0.7), based on the inner jet radius and inner fluid properties, the compound jet can break up into drops in various modes: inner dripping–outer dripping (dripping), inner jetting–outer jetting (jetting), and mixed dripping–jetting. Decreasing Re or increasing We promotes the jetting mode. The transition from dripping to jetting is also strongly affected by the velocity ratios, U₂₁ (intermediate to inner velocities) and U₃₁ (outer to inner velocities). Increasing U₂₁ makes the inner jet thinner and stretches the outer jet and thus promotes jetting. In contrast, increasing U₃₁ thins the outer jet, and thus, when the inner jet is dripping, the outer jet can break up into drops in the mixed dripping–jetting mode. Continuously increasing U₃₁ results in thinning both inner and outer jets and thus produces small drops in the jetting mode. In addition, starting from dripping, a decrease in the interfacial tension ratio of the outer to inner interfaces results in the mixed dripping–jetting and jetting modes. These modes produce various types of drops: simple drops, and compound drops with a single inner drop (single-core compound drops) or a few inner drops (multi-core compound drops).