• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Search Issue | RSS Feeds RSS
Previous Issue Next Issue

Jun 2010

Volume 22, Issue 6, Articles (06xxxx)

Issue Cover Spotlight Figure

Phys. Fluids 22, 062003 (2010); http://dx.doi.org/10.1063/1.3450319 (10 pages)

Indira Sriram, Alexander Meyer, and Eric M. Furst
back to top
RSS Feeds
back to top Biofluid Mechanics

Flexible sheets falling in an inviscid fluid

Silas Alben

Phys. Fluids 22, 061901 (2010); http://dx.doi.org/10.1063/1.3432128 (13 pages)

Online Publication Date: 16 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We use inviscid simulations to study falling flexible sheets in the two-parameter space of sheet density and bending rigidity. The basic behavior is a repeated series of accelerations to a critical speed at which the sheet flexes and rapidly decelerates, shedding large vortices. The maximum and average speeds of the sheet are closely related to the critical flutter speed. The sheet trajectories also show persistent circling, quasiperiodic flapping, and more complex repeated patterns. For small bending rigidity, the motion becomes less regular. At intermediate bending rigidity, trajectories show a well-defined falling angle relative to the vertical. For larger sheet density and bending rigidity, the overall motion is more horizontal.
Show PACS
47.11.-j Computational methods in fluid dynamics
47.15.ki Inviscid flows with vorticity
47.32.-y Vortex dynamics; rotating fluids

A boundary element model of the transport of a semi-infinite bubble through a microvessel bifurcation

Andres J. Calderon, Brijesh Eshpuniyani, J. Brian Fowlkes, and Joseph L. Bull

Phys. Fluids 22, 061902 (2010); http://dx.doi.org/10.1063/1.3442829 (11 pages) | Cited 4 times

Online Publication Date: 29 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Motivated by a developmental gas embolotherapy technique for selective occlusion of blood flow to tumors, we examined the transport of a pressure-driven semi-infinite bubble through a liquid-filled bifurcating channel. Homogeneity of bubble splitting as the bubble passes through a vessel bifurcation affects the degree to which the vascular network near the tumor can be uniformly occluded. The homogeneity of bubble splitting was found to increase with bubble driving pressure and to decrease with increased bifurcation angle. Viscous losses at the bifurcation were observed to affect the bubble speed significantly. The potential for oscillating bubble interfaces to induce flow recirculation and impart high stresses on the vessel endothelium was also observed.
Show PACS
87.19.rh Fluid transport and rheology
47.60.Dx Flows in ducts and channels
47.55.dd Bubble dynamics
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking
47.11.Hj Boundary element methods
47.63.Cb Blood flow in cardiovascular system
back to top Micro- and Nanofluid Mechanics

Gas-flow animation by unsteady heating in a microchannel

A. Manela and N. G. Hadjiconstantinou

Phys. Fluids 22, 062001 (2010); http://dx.doi.org/10.1063/1.3437602 (12 pages) | Cited 5 times

Online Publication Date: 18 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We study the flow-field generated in a one-dimensional wall-bounded gas layer due to an arbitrary small-amplitude time variation in the temperature of its boundaries. Using the Fourier transform technique, analytical results are obtained for the slip-flow/Navier–Stokes limit. These results are complemented by low-variance simulations of the Boltzmann equation, which are useful for establishing the limits of the slip-flow description, as well as for bridging the gap between the slip-flow analysis and previously developed free-molecular analytical predictions. Results are presented for both periodic (sinusoidal) and nonperiodic (step-jump) heating profiles. Our slip-flow solution is used to elucidate a singular limit reported in the literature for oscillatory heating of a dynamically incompressible fluid.
Show PACS
47.45.Gx Slip flows and accommodation
47.60.Dx Flows in ducts and channels
47.11.-j Computational methods in fluid dynamics
47.10.ad Navier-Stokes equations
02.30.Nw Fourier analysis

Absolute to convective instability transition in charged liquid jets

José M. López-Herrera, Alfonso M. Gañán-Calvo, and Miguel A. Herrada

Phys. Fluids 22, 062002 (2010); http://dx.doi.org/10.1063/1.3446972 (9 pages) | Cited 2 times

Online Publication Date: 21 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We show that the presence of electric charge at the interface of a capillary liquid jet plays a secondary role concerning the onset of an absolute or a convective instability in common operational conditions for cone-jet electrosprays, compared to other factors such as the convective velocity, jet diameter, surface tension γ, density ρ, or viscosity μ. Thus, in most situations, the critical convective velocity (or its related dimensionless number, the critical Weber number Wecr) at the threshold between the dripping and the jetting regimes depends mainly on the viscosity of the fluid, scaled as a Reynolds number Re, and not so importantly on the electric forces at the interface of the jet. Accordingly, for any liquid, the classical curve of Leib and Goldstein [Phys. Fluids 29, 952 (1986)] for Wecr versus Re is accurate enough to explore the parametrical conditions where a steady cone-jet mode is to be expected, linked to the convectively unstable nature of the issued jet. However, at the limit of low Reynolds numbers, the stability behavior becomes strongly sensitive to the electrical conductivity of the liquid. Thus, a parametrical region where a charged capillary jet becomes strongly stabilized by the viscous damping against the destabilizing surface electrical forces is described in detail in this work. The “unconditional jetting” limit previously described for a capillary jet surrounded by a coflowing liquid [ A. M. Gañán-Calvo, Phys. Rev. E 78, 026304 (2008) ] is here recovered in the absence of a coflowing fluid when “frozen” surface charges are present.
Show PACS
47.20.Dr Surface-tension-driven instability
73.40.-c Electronic transport in interface structures
68.03.Cd Surface tension and related phenomena
66.20.-d Viscosity of liquids; diffusive momentum transport
47.60.Kz Flows and jets through nozzles
47.27.te Turbulent convective heat transfer

Active microrheology of a colloidal suspension in the direct collision limit

Indira Sriram, Alexander Meyer, and Eric M. Furst

Phys. Fluids 22, 062003 (2010); http://dx.doi.org/10.1063/1.3450319 (10 pages) | Cited 12 times

Online Publication Date: 22 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The single-point active nonlinear microrheology of a colloidal suspension is measured using laser tweezers in the limit that the diameter of the probe particle approaches the diameter of the bath suspension particles. The microviscosity thins as the probe velocity (and corresponding microrheological Péclet number) increases. This thinning behavior correlates with the development of a nonequilibrium suspension microstructure surrounding the probe particle, in which a boundary layer forms on the upstream face of the probe and a wake depleted of bath particles trails the probe. The magnitude of the microviscosities and the thinning behavior are in good agreement with Brownian dynamics simulations reported by Carpen and Brady [J. Rheol. 49, 1483 (2005)] . The microviscosity increment collapses onto a single curve for all volume fractions when scaled by the contact distribution of bath particles around the probe. Scaling the microviscosity increment yields values lower than the dilute theory; furthermore, it plateaus at significantly higher Péclet numbers. The latter effect is corrected by rescaling the Péclet number with the suspension collective diffusion coefficient in place of the bath particle self-diffusivity. The magnitude of the microviscosity increment suggests the theory overestimates the frequency of bath-probe collisions. The presence and role of hydrodynamic interactions and the effect of the soft repulsive potential are discussed.
Show PACS
47.57.Qk Rheological aspects
47.57.eb Diffusion and aggregation
66.20.Ej Studies of viscosity and rheological properties of specific liquids
47.55.-t Multiphase and stratified flows
82.70.Dd Colloids
82.70.Kj Emulsions and suspensions

Variational derivation of second-order slip coefficients on the basis of the Boltzmann equation for hard-sphere molecules

Carlo Cercignani and Silvia Lorenzani

Phys. Fluids 22, 062004 (2010); http://dx.doi.org/10.1063/1.3435343 (8 pages) | Cited 11 times

Online Publication Date: 24 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The objective of the present paper is to provide an analytic expression for the first- and second-order velocity slip coefficients. Therefore, gas flow rates in microchannels have been rigorously evaluated in the near-continuum limit by means of a variational technique which applies to the integrodifferential form of the Boltzmann equation based on the true linearized collision operator. The diffuse-specular reflection condition of Maxwell’s type has been considered in order to take into account the influence of the accommodation coefficient on the slip parameters. The polynomial form of the Knudsen number obtained for the Poiseuille mass flow rate and the values of the velocity slip coefficients, found on the basis of our variational solution of the linearized Boltzmann equation for hard-sphere molecules, are analyzed in the frame of potential applications of classical continuum numerical tools in simulations of microscale flows.
Show PACS
47.85.Np Fluidics
47.11.-j Computational methods in fluid dynamics
47.45.-n Rarefied gas dynamics
47.60.Dx Flows in ducts and channels
back to top Interfacial Flows

Drop dynamics after impact on a solid wall: Theory and simulations

Jens Eggers, Marco A. Fontelos, Christophe Josserand, and Stéphane Zaleski

Phys. Fluids 22, 062101 (2010); http://dx.doi.org/10.1063/1.3432498 (13 pages) | Cited 14 times

Online Publication Date: 11 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We study the impact of a fluid drop onto a planar solid surface at high speed so that at impact, kinetic energy dominates over surface energy and inertia dominates over viscous effects. As the drop spreads, it deforms into a thin film, whose thickness is limited by the growth of a viscous boundary layer near the solid wall. Owing to surface tension, the edge of the film retracts relative to the flow in the film and fluid collects into a toroidal rim bounding the film. Using mass and momentum conservation, we construct a model for the radius of the deposit as a function of time. At each stage, we perform detailed comparisons between theory and numerical simulations of the Navier–Stokes equation.
Show PACS
47.55.D- Drops and bubbles
47.11.-j Computational methods in fluid dynamics
47.27.nb Boundary layer turbulence

A numerical study of thermocapillary migration of a small liquid droplet on a horizontal solid surface

Huy-Bich Nguyen and Jyh-Chen Chen

Phys. Fluids 22, 062102 (2010); http://dx.doi.org/10.1063/1.3432848 (12 pages) | Cited 6 times

Online Publication Date: 15 June 2010

Full Text: Read Online (HTML) | Download PDF


See Also: Erratum

Show Abstract
In the present study, the transient thermocapillary migration of a small liquid droplet on a horizontal solid surface is numerically investigated. The droplet has a large static contact angle and a high aspect ratio of the maximum height of the droplet to its footprint. The Navier–Stokes and energy equations for both the droplet and surrounding air are solved through the finite element method. The evolution of the isotherms, the flow fields and the contact angle hysteresis are presented. Two asymmetric thermocapillary vortices appear inside the droplet. The variation of the size of the thermocapillary vortex during the migration process causes the speed of the droplet to first increase significantly, and then decrease gradually to approach a constant value. The higher imposed temperature gradient causes the droplet velocity to reach its maximal value earlier and have a higher final speed. If the static contact angle of the droplet is less than (or higher) than 90°, the droplet speed is lower (or higher) since the net thermocapillary momentum in the horizontal direction is diminished (or enhanced) by the presence of capillary force. The present results for the migration velocity and the contact angle hysteresis for a squalane droplet are also in good agreement with the previous experimental results.
Show PACS
47.55.nb Capillary and thermocapillary flows
47.32.-y Vortex dynamics; rotating fluids
47.11.Fg Finite element methods
47.10.ad Navier-Stokes equations
68.03.Cd Surface tension and related phenomena
47.55.dm Thermocapillary effects
back to top Particulate, Multiphase, and Granular Flows

Dispersion and temperature statistics of inertial particles in isotropic turbulence

Saensuk Wetchagarun and James J. Riley

Phys. Fluids 22, 063301 (2010); http://dx.doi.org/10.1063/1.3392772 (15 pages)

Online Publication Date: 4 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The dispersion and temperature distribution of inertial particles are important in many turbulent, multiphase flow problems. In order to understand these better, direct numerical simulations (DNSs) are performed for inertial particles in a fluid with a constant temperature gradient and whose motion is either statistically stationary or decaying, isotropic turbulence. It is found that, for long times, the dispersion of inertial particles is the greatest when the Stokes number, Stη = τp/τη, is of order 1, where τp and τη are, respectively, the particle response time and the flow Kolmogorov time scale. A similar result is found for the long time behavior of the time rate of change of the mean-square particle temperature fluctuations, dTp′2〉/dt. To understand the DNS results better, an evolution equation for Tp′2, along with the short and long time limits, is derived analytically from the thermal energy equation for inertial particles.
Show PACS
47.55.Kf Particle-laden flows
47.27.Gs Isotropic turbulence; homogeneous turbulence
47.27.eb Statistical theories and models

Investigation and modeling of bubble-bubble interaction effect in homogeneous bubbly flows

Jung Hee Seo, Sanjiva K. Lele, and Gretar Tryggvason

Phys. Fluids 22, 063302 (2010); http://dx.doi.org/10.1063/1.3432503 (18 pages) | Cited 4 times

Online Publication Date: 15 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The effect of bubble-bubble interaction in homogeneous bubbly flow is investigated by direct numerical simulation and a bubbly mixture model for bubbly shock flows at void fraction 0.4%–13%. It is found that the bubble-bubble interaction effect is significant at void fraction higher than O(1)% and decreases the amplitude and wavelength of the macroscale oscillations in the dispersive shock structure. For the modeling of bubble-bubble interaction effect, the locally volume averaged Rayleigh–Plesset (LVARP) equation, which is an extended version of the original Rayleigh–Plesset equation, is proposed in the present study. The results of bubbly mixture model using LVARP agree well with the direct simulation results for bubbly shock flows at void fraction up to 13%. The bubble-bubble interaction in nonuniform bubbly flows is also investigated in bubbly flows with randomized initial bubble positions. It is found that the LVARP model predicts the ensemble averaged behavior with reasonable accuracy.
Show PACS
47.55.D- Drops and bubbles
47.11.-j Computational methods in fluid dynamics
47.40.Nm Shock wave interactions and shock effects
back to top Laminar Flows

Counterflow driven by swirl decay

Vladimir N. Shtern and Anatoli A. Borissov

Phys. Fluids 22, 063601 (2010); http://dx.doi.org/10.1063/1.3407646 (8 pages) | Cited 2 times

Online Publication Date: 4 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The global meridional circulation of a viscous fluid, caused by swirl decay in a cylindrical container, is studied. To this end, a new solution to the Navier–Stokes equations is obtained, and simple experiments are performed to verify the predictions of the theory. The swirl decay mechanism explains elongated counterflows in hydrocyclones and vortex tubes sometimes extending over a hundred diameters.
Show PACS
47.10.ad Navier-Stokes equations
47.32.-y Vortex dynamics; rotating fluids
47.27.N- Wall-bounded shear flow turbulence

Acoustic streaming in simplified liquid rocket engines with transverse mode oscillations

Sean R. Fischbach, Gary A. Flandro, and Joseph Majdalani

Phys. Fluids 22, 063602 (2010); http://dx.doi.org/10.1063/1.3407663 (21 pages) | Cited 2 times

Online Publication Date: 7 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
This study considers a simplified model of a liquid rocket engine in which uniform injection is imposed at the faceplate. The corresponding cylindrical chamber has a small length-to-diameter ratio with respect to solid and hybrid rockets. Given their low chamber aspect ratios, liquid thrust engines are known to experience severe tangential and radial oscillation modes more often than longitudinal ones. In order to model this behavior, tangential and radial waves are superimposed onto a basic mean-flow model that consists of a steady, uniform axial velocity throughout the chamber. Using perturbation tools, both potential and viscous flow equations are then linearized in the pressure wave amplitude and solved to the second order. The effects of the headwall Mach number are leveraged as well. While the potential flow analysis does not predict any acoustic streaming effects, the viscous solution carried out to the second order gives rise to steady secondary flow patterns near the headwall. These axisymmetric, steady contributions to the tangential and radial traveling waves are induced by the convective flow motion through interactions with inertial and viscous forces. We find that suppressing either the convective terms or viscosity at the headwall leads to spurious solutions that are free from streaming. In our problem, streaming is initiated at the headwall, within the boundary layer, and then extends throughout the chamber. We find that nonlinear streaming effects of tangential and radial waves act to alter the outer solution inside a cylinder with headwall injection. As a result of streaming, the radial wave velocities are intensified in one-half of the domain and reduced in the opposite half at any instant of time. Similarly, the tangential waves are either enhanced or weakened in two opposing sectors that are at 90° angle to the radial velocity counterparts. The second-order viscous solution that we obtain clearly displays both an oscillating and a steady flow component. The steady part can be an important contributor to wave steepening, a mechanism that is often observed during the onset of acoustic instability.
Show PACS
47.32.Ef Rotating and swirling flows
47.27.nb Boundary layer turbulence

Particle image velocimetry measurements of the interaction of synthetic jets with a zero-pressure gradient laminar boundary layer

Mark Jabbal and Shan Zhong

Phys. Fluids 22, 063603 (2010); http://dx.doi.org/10.1063/1.3432133 (17 pages) | Cited 3 times

Online Publication Date: 21 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
An experimental investigation of the interaction between a synthetic jet actuator and a zero-pressure gradient laminar boundary layer is reported. The aim of this study is to quantify the impact of synthetic jet vortical structures; namely, hairpin vortices, stretched vortex rings and tilted vortex rings on a boundary layer, and to assess their relative potential for flow separation control. Streamwise particle image velocimetry was employed in a water flume (free stream boundary layer thickness Reynolds number of 500 and boundary layer thickness-to-jet orifice diameter ratio of 4) to obtain phase- and time-averaged boundary layer profile information of the impact of synthetic jets near the wall. The potential for flow control was assessed by analyzing near wall fluid mixing, realized by the measure of increase in wall shear stress produced by a passing vortex. Hairpin vortices (produced at a jet-to-free stream velocity ratio, VR = 0.32 and dimensionless stroke length, L = 1.6) and stretched vortex rings (VR = 0.27; L = 2.7) exhibit characteristics akin to a streamwise vortex pair with a common upwash. Conversely, tilted vortex rings (VR = 0.54; L = 2.7) induce a streamwise vortex pair in the near wall region with a common downwash. Wall shear stress measurements show that synthetic jets composed of stretched vortex rings offer the best combination of near wall fluid mixing, persistency, and low rms fluctuations for potential applications of flow separation control.
Show PACS
47.15.Uv Laminar jets
47.80.Jk Flow visualization and imaging
47.85.L- Flow control
47.32.Ff Separated flows
47.15.Cb Laminar boundary layers
back to top Instability and Transition

Effect of thermally induced perturbation in supersonic boundary layers

Hong Yan and Datta Gaitonde

Phys. Fluids 22, 064101 (2010); http://dx.doi.org/10.1063/1.3432513 (16 pages) | Cited 3 times

Online Publication Date: 15 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
This paper investigates the mechanism of steady and unsteady thermal perturbation (also denoted as thermal bump) in a Mach 1.5 flat plate boundary layer. A high-fidelity upwind-biased third-order Roe scheme is used with the compressive van Leer harmonic limiter on a suitably refined mesh. The study consists of two parts. In the first part, the effects of the steady and pulsed thermal bumps are explored. It is shown that the finite-span thermal bumps generate streamwise vortices. With steady heating, the disturbance decays downstream. However, when the thermal bump is pulsed, vortex shedding is observed and the streamwise vortical disturbance grows with downstream distance, consistent with linear stability analysis. The integrated disturbance energy indicates that streamwise kinetic disturbance energy growth dominates over those associated with other two velocity and thermodynamic components. The second part of this paper explores the physical consequences of the nonlinear dynamics between the vortices produced by the pulsed bump and the compressible boundary layer. The resulting three-dimensional flow distortion generates hairpin structures which are aligned in the streamwise direction, suggesting that the transition process bears some similarity to K-type breakdown. The arrangement of these vortices is connected to the low-speed streaks observed in the evolving boundary layer. The shape factor, velocity, and Reynolds stress profiles suggest that the perturbed flow shows initiation of transition to turbulence, but remains transitional at the end of the plate.
Show PACS
47.40.Ki Supersonic and hypersonic flows
47.27.nb Boundary layer turbulence
47.32.-y Vortex dynamics; rotating fluids

Rayleigh–Taylor instability in dielectric fluids

Amey Joshi, M. C. Radhakrishna, and N. Rudraiah

Phys. Fluids 22, 064102 (2010); http://dx.doi.org/10.1063/1.3435342 (10 pages) | Cited 1 time

Online Publication Date: 16 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Force on dielectric fluids in the presence of a nonuniform electric field is shown to reduce their specific weights. An appropriately chosen field gradient makes the specific weights of superposed fluids equal and prevents Rayleigh–Taylor instability. We derive the dispersion relation for perturbation at the interface of superposed dielectric fluids, within limits of linear theory, successively for ideal, Newtonian, and those with stratified viscosity. A dimensionless dielectric number is shown to determine the stability of the arrangement.
Show PACS
47.20.Ma Interfacial instabilities (e.g., Rayleigh-Taylor)
47.65.-d Magnetohydrodynamics and electrohydrodynamics

Electrohydrodynamic instabilities at interfaces subjected to alternating electric field

P. Gambhire and R. M. Thaokar

Phys. Fluids 22, 064103 (2010); http://dx.doi.org/10.1063/1.3431043 (16 pages) | Cited 7 times

Online Publication Date: 23 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Instabilities at the interface of two immiscible fluids, either perfect or leaky dielectrics, subjected to alternating electric fields, is studied using a linear stability analysis in the limit of the electrode spacing being large compared to the wavelength of the perturbation. The Floquet analysis of the stability of this system indicates a significant effect of the frequency on the value of smax, the growth rate of the fastest growing instabilities and ETaylor, the minimum field required to excite an instability. It is seen that alternating fields act to damp the system instabilities compared to the direct current (dc) case. Moreover, the growth rate of the instabilities can be tuned from that of leaky dielectric fluids subjected to dc fields, in the low frequency limit, to that of perfect dielectrics in the high frequency limit. It is also observed that for a leaky dielectric-leaky dielectric interface, the alternating current (ac) fields can induce instabilities in a system which is stable at zero frequency, by increasing the frequency of the applied voltage.
Show PACS
47.65.-d Magnetohydrodynamics and electrohydrodynamics
47.55.N- Interfacial flows
47.20.Ma Interfacial instabilities (e.g., Rayleigh-Taylor)
47.11.-j Computational methods in fluid dynamics
02.60.-x Numerical approximation and analysis
77.84.Nh Liquids, emulsions, and suspensions; liquid crystals

Absolute lateral instability in capillary coflowing jets

Miguel A. Herrada, Conrado Ferrera, José M. Montanero, and Alfonso M. Gañán-Calvo

Phys. Fluids 22, 064104 (2010); http://dx.doi.org/10.1063/1.3447800 (10 pages) | Cited 6 times

Online Publication Date: 24 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We examine the stability of coflowing capillary jets under lateral (m = 1) perturbations of small amplitude. Two models are considered for the perturbed basic flow: the Kelvin–Helmholtz (KH) and the outer boundary layer (OBL) models. We revisit the temporal analysis of the KH model and show that the flow is unstable if and only if the (conveniently defined) Weber number is greater than unity. On the contrary, the OBL flow becomes unstable for Weber numbers much smaller than unity, although the growth rate of the perturbations is very small in that case. The spatiotemporal analysis of the dispersion relations shows that both the KH and OBL flows become absolutely unstable (absolute whipping) for sufficiently large values of the Weber number and the ratio between the outer and inner stream velocities. Absolute whipping dominates the behavior of high-viscosity jets for large velocity ratios and prevents the jetting regime from being reached even when varicose perturbations are convected downstream. For sufficiently large values of the Reynolds number, the flow becomes absolutely unstable if the velocity ratio exceeds a critical value, which is almost independent of the Weber number. For small values of the velocity ratio, the flow is stable or at most convectively unstable independently of the Reynolds and Weber numbers. For sufficiently large values of the velocity ratio, there is a critical Reynolds number above which jetting can not be reached because the flow becomes absolutely unstable due to the modes m = 0 and/or m = 1. That critical Reynolds number decreases as the velocity ratio increases. These results have important implications in technological applications such as steady high-viscosity liquid microjet production and fiber spinning using coflowing gas conformation.
Show PACS
47.55.nb Capillary and thermocapillary flows
47.60.Kz Flows and jets through nozzles
47.20.Ib Instability of boundary layers; separation
47.20.Ft Instability of shear flows (e.g., Kelvin-Helmholtz)

Global and local instability of flow focusing: The influence of the geometry

Emilio J. Vega, José M. Montanero, Miguel A. Herrada, and Alfonso M. Gañán-Calvo

Phys. Fluids 22, 064105 (2010); http://dx.doi.org/10.1063/1.3450321 (10 pages) | Cited 9 times

Online Publication Date: 30 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
In the flow focusing technique, a liquid flow rate Q is injected through a microcapillary to form a meniscus attached to its edge. The meniscus is stretched until a thin jet tapers from its tip due to the action of a gas stream driven by a pressure drop Δp. Both the liquid jet and the gas stream cross the orifice of a plate located in front of the capillary at a distance H. In the present work, the stability of both the tapering liquid meniscus and the emitted jet is analyzed experimentally. Three regimes are identified: (i) the steady jetting regime, where the liquid meniscus is stable and the jet is convectively unstable; (ii) the local instability regime, where the liquid meniscus is stable and the jet is absolutely unstable; and (iii) the global instability regime, where the liquid meniscus is unstable. The mechanisms responsible for the transitions between those regimes are described. The experiments show the existence of a minimum value Qmin of the flow rate Q below which flow focusing is globally unstable independent of the pressure drop Δp applied to the gas stream. The dependence of the stability threshold Qmin with respect to the capillary-to-orifice distance H is analyzed considering different liquids. If the rest of the geometrical parameters are fixed, there is an optimum value Hopt of the capillary-to-orifice distance H for which the stability threshold Qmin is minimum. We also determine the dependence of Hopt and the corresponding minimum flow rate Qopt with respect to the capillary diameter. In addition, we find that Qmin diverges as the capillary-to-orifice distance H decreases and approaches a certain critical value, at which the transition from flow focusing to “flow blurring” takes place. We confirm our interpretation of the experimental results by conducting numerical simulations for the aforementioned three regimes.
Show PACS
47.20.Ma Interfacial instabilities (e.g., Rayleigh-Taylor)
47.60.Kz Flows and jets through nozzles
47.11.-j Computational methods in fluid dynamics
back to top Turbulent Flows

On the measurement of vortex filament lifetime statistics in turbulence

Luca Biferale, Andrea Scagliarini, and Federico Toschi

Phys. Fluids 22, 065101 (2010); http://dx.doi.org/10.1063/1.3431660 (5 pages) | Cited 2 times

Online Publication Date: 4 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
A numerical study of turbulence seeded with light particles is presented. We analyze the statistical properties of coherent, small-scale structures by looking at the trapping events of light particles inside vortex filaments. We study the properties of particles attracting set, measuring its fractal dimension and the probability that the separation between two particles remains within the dissipative scale, even for time lapses as long as the large-scale correlation time, TL. We show how to estimate the vortex lifetime by studying the moment of inertia of bunches of particles, showing the presence of an exponential lifetime distribution, with events up to TL.
Show PACS
47.32.-y Vortex dynamics; rotating fluids
47.27.eb Statistical theories and models
47.55.-t Multiphase and stratified flows
47.53.+n Fractals in fluid dynamics
02.60.Cb Numerical simulation; solution of equations

An analysis of superhydrophobic turbulent drag reduction mechanisms using direct numerical simulation

Michael B. Martell, Jonathan P. Rothstein, and J. Blair Perot

Phys. Fluids 22, 065102 (2010); http://dx.doi.org/10.1063/1.3432514 (13 pages) | Cited 9 times

Online Publication Date: 11 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Superhydrophobic surfaces combine hydrophobic surface chemistry with topological microfeatures. These surfaces have been shown to provide drag reduction in laminar and turbulent flows. In this work, direct numerical simulation is used to investigate the drag reducing performance of superhydrophobic surfaces in turbulent channel flow. Slip velocities, wall shear stresses, and Reynolds stresses are determined for a variety of superhydrophobic surface microfeature geometry configurations at friction Reynolds numbers of Reτ ≈ 180, Reτ ≈ 395, and Reτ ≈ 590. This work provides evidence that superhydrophobic surfaces are capable of reducing drag in turbulent flow situations by manipulating the laminar sublayer. For the largest microfeature spacing, an average slip velocity over 80% of the bulk velocity is obtained, and the wall shear stress reduction is found to be greater than 50%. The simulation results suggest that the mean velocity profile near the superhydrophobic wall continues to scale with the wall shear stress and the log layer is still present, but both are offset by a slip velocity that is primarily dependent on the microfeature spacing.
Show PACS
47.11.-j Computational methods in fluid dynamics
47.27.nd Channel flow
47.60.Dx Flows in ducts and channels
47.85.lb Drag reduction
68.08.Bc Wetting
Author Select

Wall-bounded turbulent flows at high Reynolds numbers: Recent advances and key issues

I. Marusic, B. J. McKeon, P. A. Monkewitz, H. M. Nagib, A. J. Smits, and K. R. Sreenivasan

Phys. Fluids 22, 065103 (2010); http://dx.doi.org/10.1063/1.3453711 (24 pages) | Cited 63 times

Online Publication Date: 29 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Wall-bounded turbulent flows at high Reynolds numbers have become an increasingly active area of research in recent years. Many challenges remain in theory, scaling, physical understanding, experimental techniques, and numerical simulations. In this paper we distill the salient advances of recent origin, particularly those that challenge textbook orthodoxy. Some of the outstanding questions, such as the extent of the logarithmic overlap layer, the universality or otherwise of the principal model parameters such as the von Kármán “constant,” the parametrization of roughness effects, and the scaling of mean flow and Reynolds stresses, are highlighted. Research avenues that may provide answers to these questions, notably the improvement of measuring techniques and the construction of new facilities, are identified. We also highlight aspects where differences of opinion persist, with the expectation that this discussion might mark the beginning of their resolution.
Show PACS
47.27.nb Boundary layer turbulence
47.27.nd Channel flow
47.15.Cb Laminar boundary layers
47.15.Rq Laminar flows in cavities, channels, ducts, and conduits
47.60.Dx Flows in ducts and channels

Direct numerical simulations of turbulence subjected to a straining and destraining cycle

P. Gualtieri and C. Meneveau

Phys. Fluids 22, 065104 (2010); http://dx.doi.org/10.1063/1.3453709 (15 pages) | Cited 1 time

Online Publication Date: 30 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
In many turbulent flows, significant interactions between fluctuations and mean velocity gradients occur in nonequilibrium conditions, i.e., the turbulence does not have sufficient time to adjust to changes in the velocity gradients applied by the large scales. The simplest flow that retains such physics is the time dependent homogeneous strain flow. A detailed experimental study of initially isotropic turbulence subjected to a straining and destraining cycle was reported by Chen et al. [“Scale interactions of turbulence subjected to a straining-relaxation-destraining cycle,” J. Fluid Mech. 562, 123 (2006)] . Direct numerical simulation (DNS) of the experiment of Chen et al. [“Scale interactions of turbulence subjected to a straining-relaxation-destraining cycle,” J. Fluid Mech. 562, 123 (2006)] is undertaken, applying the measured straining and destraining cycle in the DNS. By necessity, the Reynolds number in the DNS is lower. The DNS study provides a complement to the experimental one including time evolution of small-scale gradients and pressure terms that could not be measured in the experiments. The turbulence response is characterized in terms of velocity variances, and similarities and differences between the experimental data and the DNS results are discussed. Most of the differences can be attributed to the response of the largest eddies, which, even if are subjected to the same straining cycle, evolve under different conditions in the simulations and experiment. To explore this issue, the time evolution of different initial conditions parametrized in terms of the integral scale is analyzed in computational domains with different aspect ratios. This systematic analysis is necessary to minimize artifacts due to unphysical confinement effects of the flow. The evolution of turbulent kinetic energy production predicted by DNS, in agreement with experimental data, provides a significant backscatter of kinetic energy during the destraining phase. This behavior is explained in terms of Reynolds stress anisotropy and nonequilibrium conditions. From the DNS, a substantial persistency of anisotropy is observed up to small scales, i.e., at the level of velocity gradients. Due to the time dependent deformation, we find that the major contribution in the Reynolds stresses budget is provided by the production term and by the pressure/strain correlation, resulting in large time variation of velocity intensities. The DNS data are compared with predictions from the classical Launder–Reece–Rodi isoptropic production [ B. E. Launder et al., “Progress in the development of a Reynolds stress turbulence closure,” J. Fluid Mech. 68, 537 (1975) ] Reynolds stress model, showing good agreement with some differences for the redistribution term.
Show PACS
47.27.Gs Isotropic turbulence; homogeneous turbulence
47.27.ek Direct numerical simulations
47.11.-j Computational methods in fluid dynamics
47.10.ad Navier-Stokes equations
back to top Compressible Flows

Critical ignition in rapidly expanding self-similar flows

Matei I. Radulescu and Brian M. Maxwell

Phys. Fluids 22, 066101 (2010); http://dx.doi.org/10.1063/1.3432132 (9 pages) | Cited 1 time

Online Publication Date: 9 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The generic problem of ignition of a particle undergoing an expansion given by a power law rate of decay behind a decaying shock is addressed in the present study. It is demonstrated, using a one-step Arrhenius irreversible reaction, that a sufficiently strong expansion wave can quench the reaction. The critical conditions for extinction are obtained in closed form in terms of the time scale for the expansion process and the thermochemical properties of the gas, yielding a critical Damkohler number, i.e., the ratio of the expansion time scale to the homogeneous ignition time scale, given by (γ−1)(Ea/RT)−1/n, where n is the power law exponent of the self-similar expansion. The critical ignition criteria, which are valid in the asymptotic limit n(γ−1)(Ea/RT) = O(1), were found in excellent agreement with numerical results. The applicability of the results obtained are discussed for ignition in rapidly expanding flows which occur behind decaying shock waves, as encountered in problems of detonation initiation by a Taylor–Sedov blast wave, and reacting jet startup, and for reactions in steady hypersonic flows around projectiles.
Show PACS
82.33.Vx Reactions in flames, combustion, and explosions
47.70.Fw Chemically reactive flows
47.40.Nm Shock wave interactions and shock effects
47.40.Rs Detonation waves
47.20.-k Flow instabilities
47.40.Ki Supersonic and hypersonic flows
back to top Geophysical Flows

Spherical single-roll dynamos at large magnetic Reynolds numbers

Henrik Latter and David Ivers

Phys. Fluids 22, 066601 (2010); http://dx.doi.org/10.1063/1.3453712 (12 pages)

Online Publication Date: 29 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
This paper concerns kinematic helical dynamos in a spherical fluid body surrounded by an insulator. In particular, we examine their behavior in the regime of large magnetic Reynolds number Rm, for which dynamo action is usually concentrated on a simple resonant stream surface. The dynamo eigensolutions are computed numerically for two representative single-roll flows using a compact spherical harmonic decomposition and fourth-order finite differences in radius. These solutions are then compared with the growth rates and eigenfunctions of the Gilbert and Ponty large Rm asymptotic theory [ Geophys. Astrophys. Fluid Dyn. 93, 55 (2000) ]. We find good agreement between the growth rates when Rm>104 and between the eigenfunctions when Rm>105.
Show PACS
47.65.-d Magnetohydrodynamics and electrohydrodynamics
47.32.Ef Rotating and swirling flows
47.27.Jv High-Reynolds-number turbulence
47.11.Bc Finite difference methods
back to top Others

Slow evaporation and condensation on a spherical droplet in the presence of a noncondensable gas

Shingo Kosuge, Kazuo Aoki, and Masatake Hatano

Phys. Fluids 22, 067101 (2010); http://dx.doi.org/10.1063/1.3432130 (14 pages) | Cited 1 time

Online Publication Date: 14 June 2010

Full Text: Read Online (HTML) | Download PDF

Show Abstract
A spherical droplet is placed in a binary mixture composed of the vapor of the droplet and another gas which neither evaporates nor condenses (a noncondensable gas). The mixture is in an equilibrium state at rest at infinity. A slow steady flow of the vapor caused by weak evaporation or condensation, under the influence of the noncondensable gas, is investigated on the basis of a linearized model Boltzmann equation. Numerical analyses by means of a finite-difference method are carried out for a wide range of the Knudsen number (i.e., from a large to small droplet compared to the molecular mean free path). The numerical results, together with analytical solutions for small and large Knudsen numbers, clarify the behavior the mixture, i.e., the mass- and heat-flow rates from or onto the droplet as well as spatial distributions of the macroscopic quantities, in the entire range of gas rarefaction. The solution for the steady heat transfer problem between a solid sphere and a binary gas mixture is also obtained as a byproduct.
Show PACS
47.45.-n Rarefied gas dynamics
47.55.D- Drops and bubbles
47.11.Bc Finite difference methods
Close
Google Calendar
ADVERTISEMENT

close