Insight into instabilities in burning droplets

The complex multiscale physics of nano-particle laden functional droplets in a reacting environment is of fundamental and applied significance for a wide variety of applications ranging from thermal sprays to pharmaceutics to modern day combustors using new brands of bio-fuels. Formation of homogenous nucleated bubbles at the superheat limit inside vaporizing droplets (with or without nanoparticles) represents an unstable system. Here we show that self-induced boiling in burning functional pendant droplets can produce severe volumetric shape oscillations. Internal pressure build-up due to ebullition activity ejects bubbles from the droplet domain causing undulations on the droplet surface and oscillations in bulk. Through experiments, we establish that the degree of droplet deformation depends on the frequency and intensity of these bubble expulsion events. In a distinct regime of single isolated bubble residing in the droplet, however, pre-ejection transient time is identified by Darrieus-Landau evaporative instability, where bubble-droplet system behaves as a synchronized driver-driven system with bulk bubble-shape oscillations being imposed on the droplet. The agglomeration of nanophase additives modulates the flow structures within the droplet and also influences the bubble inception and growth leading to different levels of instabilities. (C) 2014 AIP Publishing LLC.

Flame surface statistics of constant-pressure turbulent expanding premixed flames

In this paper we investigate the local flame surface statistics of constant-pressure turbulent expanding flames. First the statistics of local length ratio is experimentally determined from high-speed planar Mie scattering images of spherically expanding flames, with the length ratio on the measurement plane, at predefined equiangular sectors, defined as the ratio of the actual flame length to the length of a circular-arc of radius equal to the average radius of the flame. Assuming isotropic distribution of such flame segments we then convolute suitable forms of the length-ratio probability distribution functions (pdfs) to arrive at the corresponding area-ratio pdfs. It is found that both the length ratio and area ratio pdfs are near log-normally distributed and shows self-similar behavior with increasing radius. Near log-normality and rather intermittent behavior of the flame-length ratio suggests similarity with dissipation rate quantities which stimulates multifractal analysis. (C) 2014 AIP Publishing LLC.

Transition in vortex breakdown modes in a coaxial isothermal unconfined swirling jet

This paper reports first observations of transition in recirculation pattern from an open-bubble type axisymmetric vortex breakdown to partially open bubble mode through an intermediate, critical regime of conical sheet formation in an unconfined, co-axial isothermal swirling flow. This time-mean transition is studied for two distinct flow modes which are characterized based on the modified Rossby number (Ro(m)), i.e., Ro(m) <= 1 and Ro(m) > 1. Flow modes with Ro(m) <= 1 are observed to first undergo cone-type breakdown and then to partially open bubble state as the geometric swirl number (S-G) is increased by similar to 20% and similar to 40%, respectively, from the baseline open-bubble state. However, the flow modes with Ro(m) > 1 fail to undergo such sequential transition. This distinct behavior is explained based on the physical significance associated with Ro(m) and the swirl momentum factor (xi). In essence, xi represents the ratio of angular momentum distributed across the flow structure to that distributed from central axis to the edge of the vortex core. It is observed that xi increases by similar to 100% in the critical swirl number band where conical breakdown occurs as compared to its magnitude in the S-G regime where open bubble state is seen. This results from the fact that flow modes with Ro(m) <= 1 are dominated by radial pressure gradient due to swirl/rotational effect when compared to radial pressure deficit arising from entrainment (due to the presence of co-stream). Consequently, the imparted swirl tends to penetrate easily towards the central axis causing it to spread laterally and finally undergo conical sheet breakdown. However, the flow modes with Ro(m) > 1 are dominated by pressure deficit due to entrainment effect. This blocks the radial inward penetration of imparted angular momentum thus preventing the lateral spread of these flow modes. As such these structures fail to undergo cone mode of vortex breakdown which is substantiated by a mere 30%-40% rise in xi in the critical swirl number range. (C) 2014 AIP Publishing LLC.

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.

Kinetics of phase separation in polymer mixtures: A molecular dynamics study

We present detailed results from a molecular dynamics (MD) simulation of phase-separation kinetics in polymer mixtures. Our MD simulations naturally incorporate hydrodynamic effects. We find that polymeric phase separation (with dynamically symmetric components) is in the same universality class as segregation of simple fluids: the degree of polymerization only slows down the segregation kinetics. For d = 2 polymeric fluids, the domain growth law is L(t) similar to t(phi) with phi showing a crossover from 1/3 -> 1/2 -> 2/3. For d = 3 polymeric fluids, we see the crossover phi = 1/3 -> 1. Our MD simulations do not yet access the inertial hydrodynamic regime (with L similar to t(2/3)) of phase separation in 3-d fluids. (C) 2014 AIP Publishing LLC.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Shear-imposed falling film

The study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469-485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Effect of surfactants on the instability of a two-layer film flow down an inclined plane

The effect of insoluble surfactants on the instability of a two-layer film flow down an inclined plane is investigated based on the Orr-Sommerfeld boundary value problem. The study, focusing on Stokes flow P. Gao and X.-Y. Lu, ``Effect of surfactants on the inertialess instability of a two-layer film flow,'' J. Fluid Mech. 591, 495-507 (2007)], is further extended by including the inertial effect. The surface mode is recognized along with the interface mode. The initial growth rate corresponding to the interface mode accelerates at sufficiently long-wave regime in the presence of surface surfactant. However, the maximum growth rate corresponding to both interface and surface modes decelerates in the presence of surface surfactant when the upper layer is more viscous than the lower layer. On the other hand, when the upper layer is less viscous than the lower layer, a new interfacial instability develops due to the inertial effect and becomes weaker in the presence of interfacial surfactant. In the limit of negligible surface and interfacial tensions, respectively, two successive peaks of temporal growth rate appear in the long-wave and short-wave regimes when the interface mode is analyzed. However, in the case of the surface mode, only the long-wave peak appears. (C) 2014 AIP Publishing LLC.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

Investigation of the separated region ahead of three-dimensional protuberances on plates and cones in hypersonic flows with laminar boundary layers

Heat transfer rate and pressure measurements were made upstream of surface pro-tuberances on a flat plate and a sharp cone subjected to hypersonic flow in a conventional shock tunnel. Heat flux was measured using platinum thin-film sensors deposited on macor substrate and the pressure measurements were made using fast acting piezoelectric sensors. A distinctive hot spot with highest heat flux was obtained near the foot of the protuberance due to heavy vortex activity in the recirculating region. Schlieren flow visualization was used to capture the shock structures and the separation distance ahead of the protrusions was quantitatively measured for varying protuberance heights. A computational analysis was conducted on the flat plate model using commercial computational fluid dynamics software and the obtained trends of heat flux and pressure were compared with the experimental observation. Experiments were also conducted by physically disturbing the laminar boundary layer to check its effect on the magnitude of the hot spot heat flux. In addition to air, argon was also used as test gas so that the Reynolds number can be varied. (C) 2014 AIP Publishing LLC.

Simple three parameter equations for correlating liquid phase compositions in subcritical and supercritical systems

Two Chrastil type expressions have been developed to model the solubility of supercritical fluids/gases in liquids. The three parameter expressions proposed correlates the solubility as a function of temperature, pressure and density. The equation can also be used to check the self-consistency of the experimental data of liquid phase compositions for supercritical fluid-liquid equilibria. Fifty three different binary systems (carbon-dioxide + liquid) with around 2700 data points encompassing a wide range of compounds like esters, alcohols, carboxylic acids and ionic liquids were successfully modeled for a wide range of temperatures and pressures. Besides the test for self-consistency, based on the data at one temperature, the model can be used to predict the solubility of supercritical fluids in liquids at different temperatures. (C) 2014 Elsevier B.V. All rights reserved.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

Wall-mode instability in plane shear flow of viscoelastic fluid over a deformable solid

The linear stability analysis of a plane Couette flow of an Oldroyd-B viscoelastic fluid past a flexible solid medium is carried out to investigate the role of polymer addition in the stability behavior. The system consists of a viscoelastic fluid layer of thickness R, density rho, viscosity eta, relaxation time lambda, and retardation time beta lambda flowing past a linear elastic solid medium of thickness HR, density rho, and shear modulus G. The emphasis is on the high-Reynolds-number wall-mode instability, which has recently been shown in experiments to destabilize the laminar flow of Newtonian fluids in soft-walled tubes and channels at a significantly lower Reynolds number than that for flows in rigid conduits. For Newtonian fluids, the linear stability studies have shown that the wall modes become unstable when flow Reynolds number exceeds a certain critical value Re c which scales as Sigma(3/4), where Reynolds number Re = rho VR/eta, V is the top-plate velocity, and dimensionless parameter Sigma = rho GR(2)/eta(2) characterizes the fluid-solid system. For high-Reynolds-number flow, the addition of polymer tends to decrease the critical Reynolds number in comparison to that for the Newtonian fluid, indicating a destabilizing role for fluid viscoelasticity. Numerical calculations show that the critical Reynolds number could be decreased by up to a factor of 10 by the addition of small amount of polymer. The critical Reynolds number follows the same scaling Re-c similar to Sigma(3/4) as the wall modes for a Newtonian fluid for very high Reynolds number. However, for moderate Reynolds number, there exists a narrow region in beta-H parametric space, corresponding to very dilute polymer solution (0.9 less than or similar to beta < 1) and thin solids (H less than or similar to 1.1), in which the addition of polymer tends to increase the critical Reynolds number in comparison to the Newtonian fluid. Thus, Reynolds number and polymer properties can be tailored to either increase or decrease the critical Reynolds number for unstable modes, thus providing an additional degree of control over the laminar-turbulent transition.