Criteria for thermally induced atomization and catastrophic breakup of acoustically levitated droplet

We report thermally induced instability leading to catastrophic breakup in acoustically levitated vaporizing fuel droplets. Change in surface tension and viscosity with increase in droplet temperature causes wide fluctuations in droplet aspect ratio. If the viscous damping of aspect ratio oscillation is not strong enough, the droplet goes through unbounded stretching. If the droplet exceeds a critical Weber number locally, a bag type and capillary wave induced atomization can occur, which leads to catastrophic breakup. A stability criterion has been established based on the inhomogeneity of Bernoulli (acoustic) pressure and surface tension of the droplet in terms of a local Weber number and Ohnesorge number. This instability is thermally induced in a droplet which does not experience instabilities without heating. (C) 2012 Elsevier Ltd. All rights reserved.

Experimental studies on air-assisted impinging jet atomization

In the present study, a novel air-assisted impinging jet atomization is demonstrated. A configuration in which a gas jet is directed on to the impinging point of two liquid jets is used to improve the atomization. The effect of liquid properties such as viscosity and surface tension, angle between liquid jets and gas injection orifice diameter on spray characteristics has been experimentally studied. Backlit imaging and particle/droplet imaging and analysis techniques are utilized to characterize the sprays. The experimental results indicate that the effect of liquid viscosity is significant on the liquid sheet break up formed by the impinging jets. However, surface tension does not affect the spray structure significantly in this mode of atomization. At low liquid jet velocity, the prompt mode of atomization is observed where as atomization occurs in classical mode at higher liquid jet velocity. Results showed that variation in the angle between liquid jets do not affect the breakup phenomenon significantly. The spray angle is computed by finding the angle between the lines joining the impinging point and spray edge at an axial distance of 15 mm downstream of the impinging point from the ensemble-averaged data over 100 spray images. It was observed that effect of liquid jets impinging angle on the spray angle is higher at higher liquid velocity. Higher viscosity liquids exhibit lower spray angles. Droplet size measurements indicate a radial variation in the spray. An overall Sauter Mean Diameter (SMD) value is obtained by combining the droplet statistics at all radial locations at a fixed axial location. A very interesting trend is that the SMD is constant beyond a critical Gas to Liquid Ratio (GLR) and momentum ratio for a large variation in liquid viscosity and surface tension. This observation has important ramifications for fuel flexible systems. (C) 2013 Elsevier Ltd. All rights reserved.

Nucleation of a Stable Solid from Melt in the Presence of Multiple Metastable Intermediate Phases: Wetting, Ostwald's Step Rule, and Vanishing Polymorphs

In many systems, nucleation of a stable solid may occur in the presence of other (often more than one) metastable phases. These may be polymorphic solids or even liquid phases. Sometimes, the metastable phase might have a lower free energy minimum than the liquid but higher than the stable-solid-phase minimum and have characteristics in between the parent liquid and the globally stable solid phase. In such cases, nucleation of the solid phase from the melt may be facilitated by the metastable phase because the latter can ``wet'' the interface between the parent and the daughter phases, even though there may be no signature of the existence of metastable phase in the thermodynamic properties of the parent liquid and the stable solid phase. Straightforward application of classical nucleation theory (CNT) is flawed here as it overestimates the nucleation barrier because surface tension is overestimated (by neglecting the metastable phases of intermediate order) while the thermodynamic free energy gap between daughter and parent phases remains unchanged. In this work, we discuss a density functional theory (DFT)-based statistical mechanical approach to explore and quantify such facilitation. We construct a simple order-parameter-dependent free energy surface that we then use in DFT to calculate (i) the order parameter profile, (ii) the overall nucleation free energy barrier, and (iii) the surface tension between the parent liquid and the metastable solid and also parent liquid and stable solid phases. The theory indeed finds that the nucleation free energy barrier can decrease significantly in the presence of wetting. This approach can provide a microscopic explanation of the Ostwald step rule and the well-known phenomenon of ``disappearing polymorphs'' that depends on temperature and other thermodynamic conditions. Theory reveals a diverse scenario for phase transformation kinetics, some of which may be explored via modem nanoscopic synthetic methods.

Experimental analysis of shape deformation of evaporating droplet using Legendre polynomials

Experiments involving heating of liquid droplets which are acoustically levitated, reveal specific modes of oscillations. For a given radiation flux, certain fluid droplets undergo distortion leading to catastrophic bag type breakup. The voltage of the acoustic levitator has been kept constant to operate at a nominal acoustic pressure intensity, throughout the experiments. Thus the droplet shape instabilities are primarily a consequence of droplet heating through vapor pressure, surface tension and viscosity. A novel approach is used by employing Legendre polynomials for the mode shape approximation to describe the thermally induced instabilities. The two dominant Legendre modes essentially reflect (a) the droplet size reduction due to evaporation, and (b) the deformation around the equilibrium shape. Dissipation and inter-coupling of modal energy lead to stable droplet shape while accumulation of the same ultimately results in droplet breakup. (C) 2013 Elsevier B.V. All rights reserved.

Reduced graphene oxide induced phase miscibility in polystyrene-poly(vinyl methyl ether) blends

Graphene oxide and reduced graphene oxide (r-GO) were synthesized by wet chemistry and the effect of r-GO in PS-PVME blends was investigated here with respect to phase miscibility, intermolecular cooperativity in the glass transition region and concentration fluctuation variance by shear rheology and dielectric spectroscopy. The spinodal decomposition temperature (T-s) and correlation length were evaluated from isochronal temperature scans in shear rheology. The r-GO is shown to induce miscibility in the blends, which may lead to increased local heterogeneity in the blends, though the length of cooperatively re-arranged regions (xi) at T-g is more or less unaltered. The evolution of the phase morphology as a function of temperature was assessed using polarized optical microscopy (POM). In the case of the 60/40 PS-PVME blends with 0.25 wt% r-GO, apart from significant refinement in the morphology, retention of the interconnected ligaments of PVME was observed, even in the late stages of phase separation suggesting that the coarsening of the phase morphology has been slowed down in the presence of r-GO. This phenomenon was also supported by AFM. Surface enrichment of PVME, owing to its lower surface tension, in the demixed samples was supported by XPS scans. The interconnected network of PVME has resulted in significantly higher permittivity in the bi-phasic blends, although the concentration of r-GO is below the percolation threshold.

Multimodal shape oscillations of droplets excited by an air stream

The shape dynamics of droplets exposed to an air jet at intermediate droplet Reynolds numbers is investigated. High speed imaging and hot-wire anemometry are employed to examine the mechanism of droplet oscillation. The theory that the vortex shedding behind the droplet induces oscillation is examined. In these experiments, no particular dominant frequency is found in the wake region of the droplet. Hence the inherent free-stream disturbances prove to be driving the droplet oscillations. The modes of droplet oscillation show a band of dominant frequencies near the corresponding natural frequency, further proving that there is no particular forcing frequency involved. In the frequency spectrum of the lowest mode of oscillation for glycerol at the highest Reynolds number, no response is observed below the threshold frequency corresponding to the viscous dissipation time scale. This selective suppression of lower frequencies in the case of glycerol is corroborated by scaling arguments. The influence of surface tension on the droplet oscillation is studied using ethanol as a test fluid. Since a lower surface tension reduces the natural frequency, ethanol shows lower excited frequencies. The oscillation levels of different fluids are quantified using the droplet aspect ratio and correlated in terms of Weber number and Ohnesorge number. (C) 2014 Elsevier Ltd. All rights reserved.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

STM verification of the reduction of the Young's modulus of CdS nanoparticles at smaller sizes

We demonstrate the first STM evaluation of the Young's modulus (E) of nanoparticles (NPs) of different sizes. The sample deformation induced by tip-sample interaction has been determined using current-distance (I-Z) spectroscopy. As a result of tip-sample interaction, and the induced surface deformations, the I-z curves deviates from pure exponential dependence. Normally, in order to analyze the deformation quantitatively, the tip radius must be known. We show, that this necessity is eliminated by measuring the deformation on a substrate with a known Young's modulus (Au(111)) and estimating the tip radius, and afterwards, using the same tip (with a known radius) to measure the (unknown) Young's modulus of another sample (nanoparticles of CdS). The Young's modulus values found for 3 NP's samples of average diameters of 3.7, 6 and 7.5 nm, were E similar to 73%, 78% and 88% of the bulk value, respectively. These results are in a good agreement with the theoretically predicted reduction of the Young's modulus due to the changes in hydrostatic stresses which resulted from surface tension in nanoparticles with different sizes. Our calculation using third order elastic constants gives a reduction of E which scales linearly with 1/r (r is the NP's radius). This demonstrates the applicability of scanning tunneling spectroscopy for local mechanical characterization of nanoobjects. The method does not include a direct measurement of the tip-sample force but is rather based on the study of the relative elastic response. (C) 2014 Elsevier B.V. All rights reserved.

Water-responsive carbon nanotubes for selective detection of toxic gases

Ammonia plays an important role in our daily lives and hence its quantitative and qualitative sensing has become necessary. Bulk structure of carbon nanotubes (CNTs) has been employed to detect the gas concentration of 10 ppm. Hydrophobic CNTs were turned to hydrophilic via the application of a ramp electric field that allowed confinement of a controlled amount of water inside CNT microstructure. These samples were then also used to detect different gases. A comparative study has been performed for sensing three reducing gases, namely, ammonia, sulphur-di-oxide, and hydrogen sulphide to elaborate the selectivity of the sensor. A considerable structural bending in the bulk CNT was observed on evaporation of the confined water, which can be accounted to the zipping of individual nanotubes. However, the rate of the stress induced on these bulk microstructures increased on the exposure of ammonia due to the change in the surface tension of the confined solvent. A prototype of an alarm system has been developed to illustrate sensing concept, wherein the generated stress in the bulk CNT induces a reversible loss in electrical contact that changes the equivalent resistance of the electrical circuit upon exposure to the gas. (C) 2015 AIP Publishing LLC.

THREE-LAYER FLUID FLOW OVER A SMALL OBSTRUCTION ON THE BOTTOM OF A CHANNEL

Many boundary value problems occur in a natural way while studying fluid flow problems in a channel. The solutions of two such boundary value problems are obtained and analysed in the context of flow problems involving three layers of fluids of different constant densities in a channel, associated with an impermeable bottom that has a small undulation. The top surface of the channel is either bounded by a rigid lid or free to the atmosphere. The fluid in each layer is assumed to be inviscid and incompressible, and the flow is irrotational and two-dimensional. Only waves that are stationary with respect to the bottom profile are considered in this paper. The effect of surface tension is neglected. In the process of obtaining solutions for both the problems, regular perturbation analysis along with a Fourier transform technique is employed to derive the first-order corrections of some important physical quantities. Two types of bottom topography, such as concave and convex, are considered to derive the profiles of the interfaces. We observe that the profiles are oscillatory in nature, representing waves of variable amplitude with distinct wave numbers propagating downstream and with no wave upstream. The observations are presented in tabular and graphical forms.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.

Numerical study of laminar, standing hydraulic jumps in a planar geometry

We solve the two-dimensional, planar Navier-Stokes equations to simulate a laminar, standing hydraulic jump using a Volume-of-Fluid method. The geometry downstream of the jump has been designed to be similar to experimental conditions by including a pit at the edge of the platform over which liquid film flows. We obtain jumps with and without separation. Increasing the inlet Froude number pushes the jump downstream and makes the slope of the jump weaker, consistent with experimental observations of circular jumps, and decreasing the Reynolds number brings the jump upstream while making it steeper. We study the effect of the length of the domain and that of a downstream obstacle on the structure and location of the jump. The transient flow which leads to a final steady jump is described for the first time to our knowledge. In the moderate Reynolds number regime, we obtain steady undular jumps with a separated bubble underneath the first few undulations. Interestingly, surface tension leads to shortening of wavelength of these undulations. We show that the undulations can be explained using the inviscid theory of Benjamin and Lighthill (Proc. R. Soc. London, Ser. A, 1954). We hope this new finding will motivate experimental verification.

THE EFFECTS OF SURFACTANT ON SIMPLEX NOZZLE SPRAY BEHAVIOR AND ITS COMPARISON TO LIQUID FUELS

Pressure-swirl nozzles (simplex nozzles) are used in various field applications such as aero-engines, power generation, spray painting and agricultural irrigation. For this particular nozzle, research in the past decade has dealt with the development of numerical models for predicting droplet distribution profiles. Although these results have been valuable, the experimental results have been contradictory, therefore fundamental understanding of the influence of properties in nozzle is important. This paper experimentally investigates the effect of surfactants on breakup and coalescence. Since most of the fuels and biofuels have low surface tension compared to water, a comparative analysis between a surfactant solution and a liquid fuel is imperative. For this experimental study, a simplex nozzle characterized as flow number 0.4 will be utilized. The injection pressures will range from 0.3 - 4Mpa while altering the surface tension from 72 to 28mN/m. By applying Phase Doppler Particle Anemometry (PDPA) which is a non-intrusive laser diagnostic technique, the differences in spray characteristics due to spray surface tension can be highlighted. The average droplet diameter decreases for a low surface tension fluid in the axial direction in comparison to pure water. The average velocity of droplets is surprisingly lower in the same spray zone. Measurements made in the radial direction show no significant changes, but at the locations close to the nozzle, water droplets have larger diameter and velocity. The results indicate the breakup and coalescence regimes have been altered when surface tension is lowered. A decrease in surface tension alters the breakup length while increasing the spray angle. Moreover, higher injection pressure shortens the breakup length and decrease in overall diameter of the droplets. By performing this experimental study the fundamentals of spray dynamics, such as spray formation, liquid breakup length, and droplet breakup regimes can be observed as a function of surface tension and how a surrogate fuel compares with a real fuel for experimental purposes. This knowledge potentially will lead to designing a better atomizer or new biofuels.