Computational study of bouncing and non-bouncing droplets impacting on superhydrophobic surfaces

We numerically investigate bouncing and non-bouncing of droplets during isothermal impact on superhydrophobic surfaces. An in-house, experimentally validated, finite element method-based computational model is employed to simulate the droplet impact dynamics and transient fluid flow within the droplet. The liquid–gas interface is tracked accurately in Lagrangian framework with dynamic wetting boundary condition at three-phase contact line. The interplay of kinetic, surface and gravitational energies is investigated via systematic variation of impact velocity and equilibrium contact angle. The numerical simulations demonstrate that the droplet bounces off the surface if the total droplet energy at the instance of maximum recoiling exceeds the initial surface and gravitational energy, otherwise not. The non-bouncing droplet is characterized by the oscillations on the free surface due to competition between the kinetic and surface energy. The droplet dimensions and shapes obtained at different times by the simulations are compared with the respective measurements available in the literature. Comparisons show good agreement of numerical data with measurements, and the computational model is able to reconstruct the bouncing and non-bouncing of the droplet as seen in the measurements. The simulated internal flow helps to understand the impact dynamics as well as the interplay of the associated energies during the bouncing and non-bouncing. A regime map is proposed to predict the bouncing and non-bouncing on a superhydrophobic surface with an equilibrium contact angle of 155°, using data of 86 simulations and the measurements available in the literature. We discuss the validity of the computational model for the wetting transition from Cassie to Wenzel state on micro- and nanostructured superhydrophobic surfaces. We demonstrate that the numerical simulation can serve as an important tool to quantify the internal flow, if the simulated droplet shapes match the respective measurements utilizing high-speed photography.     

PDA measurements of droplet size and mass flux in the three-dimensional atomisation region of water jet in air cross-flow

The main objective of this paper is to develop a technique to measure the global droplet properties in the atomisation region of a water jet issuing in a high-speed air cross-flow. Knowledge of these global properties allows comparison of the break-up outcome of geometrically different water nozzles. This is achieved by extending a PDA system to enable measurements in three-dimensional droplet flows. First, the droplet size and the spatial droplet distribution are measured by the PDA method. The global droplet properties are then obtained by using the measured local mass flux as a weighting factor in integrating the local droplet size. To facilitate the measurement of mass flux in three-dimensional flows, the PDA method is extended so that the reference area for the mass flux is derived as a function of both the geometry of the measurement volume and the flow direction. In the present application of three-dimensional droplet flow (a water jet in air cross-flow), a simple method is developed to measure the three velocity components of droplets by means of a two-component PDA system. The paper outlines the measurement technique and the procedure of estimating the global droplet size and the global droplet size spectra from local droplet properties and local mass flux. Received: 26 July 1998/Accepted: 23 February 1999     

Accurate analysis of multicomponent fuel spray evaporation in turbulent flow

The aim of this paper is to perform an accurate analysis of the evaporation of single component and binary mixture fuels sprays in a hot weakly turbulent pipe flow by means of experimental measurement and numerical simulation. This gives a deeper insight into the relationship between fuel composition and spray evaporation. The turbulence intensity in the test section is equal to 10%, and the integral length scale is three orders of magnitude larger than the droplet size while the turbulence microscale (Kolmogorov scales) is of same order as the droplet diameter. The spray produced by means of a calibrated droplet generator was injected in a gas flow electrically preheated. N-nonane, isopropanol, and their mixtures were used in the tests. The generalized scattering imaging technique was applied to simultaneously determine size, velocity, and spatial location of the droplets carried by the turbulent flow in the quartz tube. The spray evaporation was computed using a Lagrangian particle solver coupled to a gas-phase solver. Computations of spray mean diameter and droplet size distributions at different locations along the pipe compare very favorably with the measurement results. This combined research tool enabled further investigation concerning the influencing parameters upon the evaporation process such as the turbulence, droplet internal mixing, and liquid-phase thermophysical properties.     

Measurement of drop size in horizontal annular flow with the immersion method

A modified version of the oil immersion sampling system developed by Okada was exploited to obtain measurements of drop diameter in a horizontal annular flow of air and water in a 9.53-cm pipe. Drop samples were captured in a high-viscosity oil and photographic images of the samples were used to measure the distribution of diameters. Advantages of the immersion method are that the drops are spherical when photographed and that the measurement does not require removal of the film from the wall. Drop size distributions are found to be similar to distributions measured using a photographic technique for vertical annular flow in a 4.2-cm pipe. Drop size distributions show some differences from measurements that used a laser diffraction system in the same horizontal pipeline. In particular, the diffraction method indicates the existence of larger maximal drop sizes than are determined with the immersion technique.     

Abnormal rotation of a deformed liquid crystal droplet immersed in an isotropic fluid after transient flow

The morphology evolution of liquid crystal droplets immersed in an isotropic fluid in flow field is found to be different from flexible polymer droplets. In this paper, we investigated the retraction of a liquid crystal droplet after transient flow. It is found that the liquid crystal droplet will rotate during the shape recovery, which has never been observed for an isotropic droplet. The factors that influence the rotational angle of a single liquid crystal droplet during retraction progress were studied, including the temperature, the dimension of the droplets, the time of shear flow, the shear rate, the flow type, and the properties of liquid crystal molecules. The rotation of liquid crystal droplet during shape recovery is ascribed to both the bulk elasticity of liquid crystal droplets and the anisotropic properties of the interface between liquid crystal and isotropic fluid.     

Particle velocimetry inside Newtonian and non-Newtonian droplets impacting a hydrophobic surface

A particle velocimetry technique is described which enables the measurement of the fluid velocity inside impacting drops. Using high speed photography of 2 μm fluorescent tracer particles suspended in the fluid, the velocity field was measured as a function of time and radial position. The potential of the technique is illustrated using velocimetry measurements of drops of pure water and aqueous solutions of 200 ppm poly-(ethylene oxide) (PEO). Dilute solutions of PEO have been known for some time to suppress the rebound of water from hydrophobic surfaces. The dissipation has traditionally been attributed to an increased extensional viscosity as the polymers stretch in the extensional flow of the droplet. Our results enable us to infer that the extensional viscosity of PEO drops, during both the spreading and retraction phase, is similar to that of pure water. The data suggest that the true source of dissipation lies at the droplet edge. We also show, by analysing the spreading of water drops, that the Roisman-Yarin theory for a droplet spreading on a surface is valid in the bulk of the droplet prior to the final stages of spreading.     

Experimental determination of droplet size and density field in condensing flows

 We report a detailed experimental characterization of the process of homogeneous condensation in supersonic expanding flow. In our experiments, the supersaturated mixture expands in a Laval nozzle, where, depending on the initial conditions, a steady or periodically oscillating flow may evolve due to the non-linear interaction of nucleation and droplet growth rate with the flow field. Two experimental techniques are utilized: holographic interferometry for the determination of the density field and a time-resolved white-light extinction method. The latter is employed to derive the evolution in time of the droplet cloud (i.e., modal radius, number density, and relative width) and to measure the frequency of oscillations. In combination with the wide-field density data, droplet size measurements provide additional physical insights in the mechanism of interaction in condensing flows and serve as an excellent test case for the critical assessment of nucleation and droplet growth theories. To this purpose, the accuracy of the measurements is carefully reviewed due to the difficulties of characterizing dense sub-micron droplet clouds by means of light-scattering techniques. An important byproduct of this analysis is an evaluation of the applicability of single-scattering approximations, i.e., Lambert-Beer law, for a variety of experimental configurations. Received: 24 April 2001 / Accepted: 29 August 2001     

Numerical modeling and experimental measurements of a high speed solid-cone water spray for use in fire suppression applications

Experimental measurements and numerical simulations of a high-speed water spray are presented. The numerical model is based on a stochastic separated flow technique that includes submodels for droplet dynamics, heat and mass transfer, and droplet–droplet collisions. Because the spray characteristics near the nozzle are difficult to ascertain, a new method for initialization of particle diameter size is developed that assumes a Rosin–Rammler distribution for droplet size, which correctly reproduces experimentally measured Sauter and arithmetic mean diameters. By relating the particle initialization to lower moments of the droplet statistics, it is possible to take advantage of measurements without substantial penalties associated with the greater experimental uncertainty of individual droplet measurements. Overall, very good agreement is observed in the comparisons of experimental measurements to computational predictions for the streamwise development of mean drop size and velocity. In addition, the importance of modeling droplet–droplet collisions is highlighted with comparison of selected droplet–droplet collision models.     

Shape and structure recovery of an LCP droplet under a large step strain: observation and stress calculation

Shape recovery of a droplet of liquid crystalline polymer (LCP) hydroxypropylcellulose in a matrix of poly(dimethyl siloxane) subjected to a step shear strain has been studied via optical microscopy. Just after application of a large strain, the LCP droplet shape is flat ellipsoid, and then the droplet takes cylindrical shape and band texture perpendicular to the flow direction appears. The band texture fades away before emergence of poly-domain structure. In the final process with the shape of spheroid, poly-domain structure recovers very slowly. Except for the final process, the shape change is identical with that of isotropic droplet at strains smaller than 3, when the LCP viscosity in Region II is taken as an equivalent viscosity for normalization. For a 20:80 blend, the excess relaxation modulus is calculated based on the Doi-Ohta theory, taking account of the distribution of droplet size and compared with experimental modulus data.     

Instantaneous and simultaneous planar velocity field measurements of two phases for turbulent mixing of high pressure sprays

This paper describes the tests of accuracy and the first application of a combined planar visualization technique. Its goal is two-phase flow discrimination, i.e. simultaneous measurements of velocity of droplets and ambient gas in the case of two-phase flow mixing, at the same location and with possible conditioning by “apparent diameter” (AD) of the droplets. It combines the mature techniques of particle image velocimetry (PIV), planar Mie scattering diffusion (PMSD), planar laser-induced fluorescence (PLIF), and it necessitates two synchronized cross-correlation systems, digital image treatment and analysis. This technique was developed with the objective of better describing the mixing between liquid and gaseous phases as in the case of high-pressure spray atomization in quiescent ambient gas. The basic principle of separation is to seed the ambient gas with micrometer particles and to tag the liquid with fluorescent dye. We use digital image treatment and analysis to discriminate between the phases. We use two cross-correlation PIV systems in order to obtain the velocity field of the droplets and gas simultaneously and separately at the same location. The digital image processing for separating the phases involves geometric measurement of droplet shapes. This leads to measurement of droplet parameters close to their real diameter, which could be used for analysis of actual mixing. A synchronized system composed of two CCD cameras is used for image recording, and two Nd:YAG lasers are used for generating pulsed light sheets at times t and t + δt. Tests were performed to check for different sources of errors. The combined technique was applied to measurements in high-pressure spray flow atomizing in a quiescent ambient gas, and first results are presented.     

Effect of fumed silica nanoparticles on the morphology and rheology of immiscible polymer blends

We examined the effect of interfacially active particles on the morphology and rheology of droplet/matrix blends of two immiscible homopolymers. Experiments were conducted on polybutadiene/polydimethylsiloxane (10/90) blend and the inverse system. The effects of fumed silica nanoparticles, at low particle loadings (0.1–2.0 wt%), were examined by direct flow visualization and by rheology. Fumed silica nanoparticles were found to significantly affect the morphology of polymer blends, inducing droplet cluster structure and decreasing the droplet size, regardless of which phase wets the particles preferentially. This is surprising in light of much past research that shows that particles are capable of bridging and thus induce droplet cluster structure in droplet/matrix systems only when they are preferentially wetted by the continuous phase. Therefore, there should exist other possible mechanisms responsible for these droplet cluster structures except for the bridging mechanism. We proposed a particle-flocculating mechanism based on the fact that fumed silica particles readily flocculate due to their high aspect ratio, fractal-like shape, or interparticle attractions. Optical microscopy also reveals that the clustering structure becomes more extensive, and the droplet sizes in the clusters become smaller when the particle loading is increased. Rheologically, the chief effect of particles is to change the flow behavior from a liquid-like rheology to gel-like behavior. This gel-like behavior can be attributed to droplet clustering. Moreover, it should be emphasized that such gel-like behavior can be seen in the blends regardless of which phase wets the particles preferentially, suggesting that, once again, bridging is not the only cause of droplet clustering.     

3D phase field modeling of electrohydrodynamic multiphase flows

A 3D phase field model is developed to investigate the electrohydrodynamic (EHD) two phase flows. The explicit finite difference method, enhanced by parallel computing, is employed to solve the coupled nonlinear governing equations for the electric field, the fluid flow field and free surface deformation. Numerical tests indicate that an appropriate interpolation of densities within the interface is critical in ensuring numerical stability for highly stratified flows. The 3D phase field model compares well with the Taylor theory for the deformation of a single dielectric droplet in an electric field. Computed results show that the deformation of a leaky dielectric droplet in an electric field undergoes various stages before it reaches the final oblate shape. This is caused by the free charge relaxation near the fluid–fluid interface. The coalescence of four droplets in an electric field illustrates a truly 3D deformation behavior and a complex evolving fluid flow field associated with the participating droplets. The coalescence is a result of combined actions produced by the global electric force, the circulatory flows generated by the local electrohydrodynamic stress and the electrically-induced deformation. The 3D phase field model is also applied in modeling of an electrohydrodynamic patterning process for manufacturing nanoscaled structures, in which complex 3D flow structures develop as the electrically-induced deformation evolves.     

Droplet heat transfer on micro-post arrays: Effect of droplet size on droplet thermal characteristics

Heat transfer towards a water droplet from hydrophobic micro-post array surface is considered while mimicking the environmental temperatures. Micro-post arrays are created on a silicon wafer surface via lithography technique. The textured surfaces are replicated by polydimethylsiloxane (PDMS) to achieve an optical transmittance. The droplet adhesion on micro-post array surface is presented and the influence of droplet size on the heat transfer and droplet internal flow characteristics is examined. The flow predictions are validated via the particle image velocimetry data. It is found that adhesion force between the water droplet and the micro-post arrays surface depends on the geometric size and the orientation of the micro-post arrays on the surface. Temperature and flow fields are influenced by the droplet size. The Nusselt and the Bond numbers increase with the droplet volume; however, the Bond number remains less than unity indicating that the Marangoni current dominates over the buoyancy current in the droplet. The Nusselt number attains larger values for micro-post array surface than that of the plain surface. This is because of temperature and velocity oscillations along the contact lines at the droplet bottom due to the pitches of the micro-post arrays.     

A STUDY ON DROPLET VELOCITY OF ETHANOL DURING ELECTROSPRAYING PROCESS AT CONE-JET MODE

研究液滴在静电喷雾下的速度特性是理解喷雾形态的形成及演化的关键.结合锥射流模式下乙醇静电喷雾实验数据,建立了静电喷雾二维轴对称模型.基于离散相液滴运动方程、连续相空气运动方程、电场方程以及用户自定义函数,进行了数值求解,获得了锥射流模式下的乙醇静电喷雾形态、空间电场分布以及液滴速度场分布.考虑了不同空气入口流速的影响,得到了乙醇/空气同轴射流静电喷雾形态的变化规律.结果表明,喷雾外围液滴与空气流场有较强的相互作用,导致喷雾中轴线附近的液滴速度分布变化较小,而在喷雾外围处的液滴速度分布沿径向剧烈变化;随着空气入口速度的增大,乙醇/空气同轴射流静电喷雾形态先趋于发散,当空气入口速度大于喷雾外围液滴轴向速度时,喷雾形态则趋于聚拢.因此,除改变施加电压、液体流量和电极结构外,通过控制空气入口速度来影响喷雾液滴速度场,也可获得不同的静电喷雾效果.     

Measurements of the vaporization dynamics in the development zone of a burning spray by planar laser induced fluorescence and Raman scattering

A linear measurement technique based on simultaneous planar imaging of laser induced dye fluorescence and Raman scattering in the liquid phase is reported. Calibrations in a stream of monosized droplets doped with weak concentrations of rhodamin show that the intensities on the droplet images are proportional to the actual droplet volume for Raman scattering and to the initial volume of the droplet for fluorescence, as the mass of dissolved dye does not vaporize. Thus, the mass fraction of liquid fuel that has vaporized before the probing event can be derived from these simultaneous measurements. Experiments are performed in the early development of a burning spray to derive cumulative information on the vaporization dynamics in terms of mass fraction or evaporation constant. Size distribution from conjoined phase-Doppler measurements are also used to derive the rate of droplet consumption along the axis of the burning spray.     

Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution

A Hamiltonian framework of non-equilibrium thermodynamics is adopted to construct a set of dynamical continuum equations for a polymer blend with matrix viscoelasticity and a narrow droplet size distribution that is assumed to obey a Weibull distribution function. The microstructure of the matrix is described in terms of a conformation tensor. The variable droplet distribution is described in terms of two thermodynamic variables: the droplet shape tensor and the number density of representative droplets. A Hamiltonian functional in terms of the thermodynamic variables is introduced and a set of time evolution equations for the system variables is derived. Sample calculations for homogenous flows and constant droplet distribution are compared with data of a PIB/PDMS blend and a HPC/PDMS blend with high viscoelastic contrast. For the PIB/PDMS blend, satisfactory predictions of the flow curves are obtained. Sample calculations for a blend with variable droplet distribution are performed and the effect of flow on the rheology, droplet morphology, and on the droplet distribution are discussed. It is found that deformation can increase or decrease the dispersity of the droplet morphology for the flows investigated herein.     

Experimental study on effect of relative humidity on heat transfer of an evaporating water droplet in air flow

Dispersed water droplets are often seen in environmental air flows in rain, cloud, mist, sea spray and so on. It is therefore of great importance to precisely estimate heat transfer between water droplets and atmospheric air in developing a reliable climate model. The purpose of this study is to fabricate the measurement system for the temperature of a small water droplet in air flow under the controlled relative humidity condition and to investigate the effect of relative humidity on heat transfer across the surface of an evaporating water droplet in air flow. The results show that the droplet temperature decreases in the low-relative-humidity condition, whereas it increases in the high-relative-humidity condition. Nusselt number on the droplet surface is not affected by the relative humidity.     

Droplet relaxation in blends with one viscoelastic component: bulk and confined conditions

Using a counter rotating parallel plate shear flow cell, the shape relaxation of deformed droplets in a quiescent matrix is studied microscopically. Both the effects of geometrical confinement and component viscoelasticity are systematically explored at viscosity ratios of 0.45 and 1.5. The flow conditions are varied from a rather low to a nearly critical Ca number. Under all conditions investigated, viscoelasticity of the droplet phase has no influence on shape relaxation, whereas matrix viscoelasticity and geometrical confinement result in a slower droplet retraction. Up to high confinement ratios, the relaxation curves for ellipsoidal droplets can be superposed onto a master curve. Confined droplets with a sigmoidal shape relax in two stages: the first consists of a shape change to an ellipsoid with a limited amount of retraction, and the second is the retraction of this ellipsoid. The latter stage can be described by means of one single relaxation time that can be obtained from the relaxation of initially ellipsoidal droplets. The experimental results are compared to the predictions of a recently published phenomenological model for droplet dynamics in confined systems with viscoelastic components (Minale et al., Langmuir 26:126–132, 2010). However, whereas the model predicts additive effects of geometrical confinement and component viscoelasticity, the experimental data reveal more complex interactions.     

Low-pressure twin-fluid atomization: Effect of mixing process on spray formation

The present work comparably examines four different twin-fluid atomizers operated under the same operating conditions. Spray formation was examined by several approaches. The internal flow pattern was estimated using a simplified analytical approach, and the results were supported by the observation of the liquid discharge in the near-nozzle region. A high-speed back illumination was used for visualisation of the primary breakup. In the region of fully developed spray, the dynamics of droplets was studied using a phase-Doppler analyser (PDA). The information obtained from all methods was then correlated. Results show that the spray formation process depends mainly on the internal design of twin-fluid atomizer at low gas to liquid ratios (GLR). The amount of gas influences the character of the internal two-phase flow, a mechanism of the liquid breakup, droplet dynamics and a resulting drop size distribution. Differences among the atomizers are reduced with the increase in GLR. Moreover, it was shown that a certain mixing process can inherently create the annular internal flow which generates a stable spray characterized by relatively low mean droplet size.