Index of content:
Volume 17, Issue 3, March 2005
- Interfacial Flows
17(2005); http://dx.doi.org/10.1063/1.1850691View Description Hide Description
Coaxial jet electrospray is a technique to generate microencapsules, which uses electric forces to create a coaxial microjet from two immiscible liquids. Compound droplets with narrow size distribution are produced after the jet breaks up. In this paper, the spraying modes are investigated experimentally with proper flow rates of the inner and outer liquids. Ethanol/glycerol/tween mixture (outer liquid) and cooking oil (inner liquid) are fed into the gap between outer and inner capillaries and the inner capillary, respectively. The spraying modes presented in our experiments are “dripping mode,” “dripping mode in spindle,” “cone-jet mode,” “pulse mode in cone,” and “multijets mode” sequentially, as the applied voltage increases. The region of stable cone-jet mode extends with decrease of the outer liquidflow rate and increase of the inner one. It is found that the spray phenomena are mainly determined by properties of the outer liquid, which is viscous and electric conductive enough. A rudimentary physical model is developed, in which both the viscosity and liquid interface tension are taken into account.
17(2005); http://dx.doi.org/10.1063/1.1846791View Description Hide Description
Liquid bridges are usually encapsulated with another liquid in float-zone crystal growth processes to ensure the containment of volatile components. This paper is concerned with the stability of an inviscid liquid bridge that is off-centered with respect to its encapsulant. Perturbation theory is used to study the stability of such a bridge subject to inertial disturbances. It is concluded that while the off-centered nature does not change the neutral point it does affect the rate of growth and decay of the disturbances causing the unstable regions to become less unstable and stable regions to become less stable. Limiting conditions are considered in order to provide a better understanding of the physics of off-centering.
17(2005); http://dx.doi.org/10.1063/1.1856532View Description Hide Description
A curvatureless three-dimensional boundary-integral algorithm has been developed for tangential Marangoni stresses and used to study breakup and capture of two deformable drops for arbitrary drop-to-medium viscosity and thermal conductivity ratios in parallel and antiparallel arrangements of gravity and an applied vertical temperature gradient. When the driving forces are opposed, the previously observed inhibition of breakup by a weak thermocapillary effect for drops with equal viscosity and thermal conductivity ratios is shown to be almost exclusively the result of changing interfacial tension, with Marangoni stresses having virtually no influence. Alignment of gravity and the temperature gradient in the same direction enhances breakup more than opposing driving forces reduce it, with the limitation that the drops are moving toward a region of zero interfacial tension. The thermal conductivity ratio has negligible impact on these interactions. For bubbles, the effect of a temperature gradient on gravitational results is much less pronounced than for drops. Under certain conditions, parallel orientation of the driving forces weakly inhibits the capture interaction for bubbles, while antiparallel orientation enhances the phenomenon.
17(2005); http://dx.doi.org/10.1063/1.1852459View Description Hide Description
The stability of the interface between two thin leaky dielectric liquid layers bounded between two flat electrodes is considered. A coupled system of evolution equations is derived for the interfacial location and charge density using lubrication theory. This system is parametrized by the dielectric constants of the two fluids in addition to ratios of their conductivities, viscosities, and thicknesses. A linear stability analysis is conducted and the behavior of the system in the nonlinear regime is also examined. The system is destabilized by electrical stresses that are resisted by capillarity and modified by viscous dissipation. Our results suggest that decreasing the thickness ratio is destabilizing, giving rise to periodic structures of decreasing wavelength. Decreasing the viscosity ratio was also found to lead to the formation of sharp-edged structures whose vertical extent is virtually equal to the gap width between the electrodes. Similar structures were also determined upon increasing the ratio of the dielectric constants and electric conductivities.
17(2005); http://dx.doi.org/10.1063/1.1859191View Description Hide Description
In reduced gravity, the stability of cylindrical liquid bridges and other systems having free surfaces is affected by ambient vibrations of the spacecraft. Such vibrations are expected to excite capillary modes. The lowest-order unstable mode of a liquid bridge is particularly susceptible to vibration as the length of the bridge approaches the stability limit. This mode is known as the (2,0) mode and is an axisymmetric varicose mode of one wavelength in the axial direction. In this work, an optical system is used to detect the (2,0)-mode amplitude. The derivative of the error signal produced by this detector is used to produce the appropriate voltages on a pair of annular disk electrodes which are concentric with the bridge. A mode-coupled Maxwell stress profile is thus generated in proportion to the modal velocity. Depending on the sign of the gain, the damping of the capillary oscillation can be either increased or decreased. This effect has been demonstrated in Plateau-tank experiments. Increasing the damping of the capillary modes on free liquid surfaces in space could be beneficial for containerless processing and other technologies.
17(2005); http://dx.doi.org/10.1063/1.1862234View Description Hide Description
We consider the temporal instability of parallel two-phase mixing layers. The viscous case is examined using a composite error-function velocity profile. The inviscid case is considered for the broken-line velocity profile, where the thickness of the boundary layer in each fluid next to the interface is chosen to match the viscous error-function profile at the interface and far away from it. Viscosity modifies the inviscid stability properties quantitatively, but we can also discern an additional unstable mode exclusively related to viscous shear. In the absence of interfacial tension, this mode dominates at large wavenumbers when the Reynolds number is sufficiently high. The various viscous modes cannot generally be attributed to either one of the phases due to mode mixing or exchange. For parameters resembling those of atomization experiments and applications, the most unstable wavelength and growthrate in the viscous case can exceed the inviscid values significantly. The viscous stability analysis also provides better agreement with recent experimental results for air and water than inviscid stability calculations.
17(2005); http://dx.doi.org/10.1063/1.1864132View Description Hide Description
The problem of nonperturbative description of stationary flames with arbitrary gas expansion is considered. On the basis of the Thomson circulation theorem an implicit integral of the flow equations is constructed. With the help of this integral, a simple explicit expression for the vortex mode of the burnt gas flow near the flame front is obtained. Furthermore, a dispersion relation for the potential mode at the flame front is written down, thus reducing the initial system of bulk equations and jump conditions for the flow variables to a set of integrodifferential equations for the flame front position and the flow velocity at the front. The developed approach is applied to the case of thin flames. Finally, an asymptotic expansion of the derived equations is carried out in the case where is the gas expansion coefficient, and a single equation for the front position is obtained in the second post-Sivashinsky approximation. It is demonstrated, in particular, how the well-known problem of correct normalization of the front velocity is resolved in our approach. It is verified also that in the first post-Sivashinsky approximation, the equation reduces to the Sivashinsky–Clavin equation corrected according to Cambray and Joulin. Analytical solutions of the derived equations are found, and compared with the results of numerical simulations.
- Viscous and Non-Newtonian Flows
17(2005); http://dx.doi.org/10.1063/1.1844911View Description Hide Description
Free surface flows of thixotropic fluids such as paints, self-compacting concrete, or natural mudflows are of noticeable practical interest. Here we study the basic characteristics of the uniform flow of a layer of thixotropic fluid under gravity. A theoretical approach relying on a simple thixotropy constitutive equation shows that after some time at rest over a small slope angle the fluid layer should start to flow rather abruptly beyond a new, larger, critical slope angle. The theory also predicts that the critical time at which the layer velocity should significantly increase is proportional to the duration of the preliminary rest and tends to infinity when the new slope approaches the critical slope. Experiments carried out with different suspensions show that the qualitative trends of the flows are in very good agreement with the theoretical predictions, except that the critical time for flow start appears to be proportional to a power 0.6 of the time of rest whereas the theory predicts a linear dependence. We show that this indicates a restructuration process at rest differing from the restructuration process under flow.
17(2005); http://dx.doi.org/10.1063/1.1850148View Description Hide Description
Experimental studies are made on the collapse of macroscopic cavity regions placed in an otherwise homogeneous granular material due to viscousflow. First, an initially circular two-dimensional hole is exposed to a uniform flow at infinity and the process of boundary shape deformation due to locally enhanced stress components is clarified. Above a certain critical velocity, particles on the upstream-side boundary lose contact with neighboring ones, and are carried to the other side of the boundary. At the same time the fluidized region develops towards upstream direction in the granular material. A simple scaling law on the reduction of the void area is obtained. Second, attention is paid to the interaction of two circular cavities of equal radii, whose center-to-center distance and angle of attack are varied. In some configuration of cavities, lowering of the critical velocity of collapse is recognized. Successive processes on the collapse of upstream-side and downstream-side cavities are classified depending on the magnitude of velocity and the configuration of cavities. A possible scenario on the growth of global scale patterns in granular material such as water vein formation and landslides is suggested.
17(2005); http://dx.doi.org/10.1063/1.1850151View Description Hide Description
In the present study, we apply a distributed (i.e., spatially varying) forcing to flow over a circular cylinder for drag reduction. The distributed forcing is realized by a blowing and suction from the slots located at upper and lower surfaces of the cylinder. The forcing profile from each slot is sinusoidal in the spanwise direction but is steady in time. We consider two different phase differences between the upper and lower blowing/suction profiles: zero (in-phase forcing) and (out-of-phase forcing). The Reynolds numbers considered are from 40 to 3900 covering various regimes of flow over a circular cylinder. For all the Reynolds numbers larger than 47, the present in-phase distributed forcing attenuates or annihilates the Kármán vortex shedding and thus significantly reduces the mean drag and the drag and lift fluctuations. The optimal wavelength and amplitude of the in-phase forcing for maximum drag reduction are also obtained for the Reynolds number of 100. It is shown that the in-phase forcing produces the phase mismatch along the spanwise direction in the vortex shedding, weakens the strength of vortical structures in the wake, and thus reduces the drag. Unlike the in-phase forcing, the out-of-phase distributed forcing does not reduce the drag at low Reynolds numbers, but it reduces the mean drag and the drag and lift fluctuations at a high Reynolds number of 3900 by affecting the evolution of the separating shear layer, although the amount of drag reduction is smaller than that by the in-phase forcing.
17(2005); http://dx.doi.org/10.1063/1.1850431View Description Hide Description
The presence of small amounts of polymer in Newtonian solvents can have a significant impact on the flow behavior of these fluids in extension-dominated flows. This study investigates the effect of elasticity on the on-demand drop formation through the use of low viscosityelastic liquids. A high-speed camera is employed to observe drops ejected from a nozzle. The drops are created using a piezoelectric sleeve that contracts around a nozzle forcing fluid out. We observe that the satellite drops commonly produced with Newtonian fluids of identical shear viscosity in the same geometry can be suppressed when polymers with sufficient molecular weight are added. For our nozzle, the minimum required molecular weight is 300k PEO at a concentration of 25 ppm. However, the resultant increased elasticity in the solution requires a greater pulse strength to eject the drop. In addition, the fluids containing polymers have a longer thread, a longer time to separation, and a lower velocity than the Newtonian fluids with similar shear viscosity.
17(2005); http://dx.doi.org/10.1063/1.1863320View Description Hide Description
Ferrofluid pipe flow in an oscillating magnetic field along the pipe axis is studied theoretically in a wide range of the flow rate. The field-dependent part of viscosity (it can be positive or negative) reveals significant dependence on the flowvorticity, i.e., ferrofluids exhibit non-Newtonian behavior. This is manifested in an alteration of the velocity profile—it ceases to be parabolic—and deviation of the flow rate from the value prescribed by Poiseuille’s formula. The presented model based on the conventional ferrohydrodynamic equations and an assumption of the ferrofluid structure fits well experimental data recently obtained by Schumacher, Sellien, Konke, Cader, and Finlayson [“Experiment and simulation of laminar and turbulent ferrofluid pipe flow in an oscillating magnetic field,” Phys. Rev. E67, 026308 (2003)].
- Laminar Flows
17(2005); http://dx.doi.org/10.1063/1.1850374View Description Hide Description
Liquid microdroplets represent a convenient system for studies of mixing by chaotic advection in discrete microscopic volumes. The mixing properties of the flows in microdroplets are governed by their symmetries, which give rise to invariant surfaces serving as barriers to transport. Thorough mixing via chaotic advection requires destruction of all such invariant surfaces. To illustrate this idea, we demonstrate that quick and thorough mixing inside a spherical microdroplet suspended in a layer of substrate fluid can be obtained by moving the droplet along a two-dimensional path using temperature-induced surface tension gradients. The use of flow invariants also provides a convenient way to analyze the mixing properties of flows in many other experimental implementations.
Sharp scalar and tensor bounds on the hydrodynamic friction and mobility of arbitrarily shaped bodies in Stokes flow17(2005); http://dx.doi.org/10.1063/1.1852315View Description Hide Description
We prove rigorous inequalities for the hydrodynamic translational friction and mobility matrices and of an arbitrarily shaped rigid particle in terms of the electrostaticcapacitance of a conducting particle of identical shape. Specifically, we derive the scalar and matrix inequalities and , where all quantities are normalized by the corresponding values for a sphere, and the mobility matrix is evaluated in the center-of-mobility reference frame. These bounds are obtained using a variational approach with the energy dissipation functional expressed in terms of the induced force distribution on the surface of the particle. To relate the hydrodynamic problem to the solution of the corresponding electrostatic problem, the trial force field is expressed in terms of the charge distribution on the equipotential particle surface. This procedure yields the first rigorous bounds on hydrodynamicfriction that apply to bodies with translation-rotation coupling. We demonstrate that the error of the Hubbard–Douglas approximation , corresponding to our scalar bound, is quadratic in the deviation of the trial induced-force field from the exact form—which explains why this relation is highly accurate for many particle shapes. Our numerical results confirm that the Hubbard–Douglas approximation is accurate for a variety of objects, including helices with translational–rotational coupling. In addition, we establish a rigorous, sharp bound on the effective (scalar) Brownian diffusion coefficient of an arbitrarily shaped particle.
- Instability and Transition
17(2005); http://dx.doi.org/10.1063/1.1848547View Description Hide Description
We consider the effect of a shock passing through an arbitrarily shaped interface between two fluids. The evolution of the interface into a new shape, written formally as , is found by applying the linear, classical Richtmyer–Meshkov instability result to each mode in the Fourier expansion of the original interface. We provide several examples where the new shape can be found analytically. For any interface we define an associated dual interface and show that . Representing a shock by a new mathematical operator we find how , and transform under the effect of a shock. Kink-singularities are found in when and where has a discontinuous change in its first derivative. These are the locations where jetting occurs. We briefly discuss the effects of nonlinearity, compressibility, viscosity, etc., all of which suppress kink-singularities, and present hydrocode simulations of shock tube and high-explosive-driven experiments to highlight the influence of compressibility, nonlinearity, and material strength.
17(2005); http://dx.doi.org/10.1063/1.1851451View Description Hide Description
The normal mode linear analysis is applied to investigate the stability of a circular isolated compressible vortex with the emphasis on studying the effect of entropy stratification of the basic flow on stability properties. We study a family of vortices that have zero total circulation, with the swirl velocity being presented by a Taylor-type distribution in the radial direction. The stratification of entropy is modeled by a Gaussian two-parametric profile whose parameters control the maximal deviation from the homentropic distribution and the extent of the entropic zone. Results presented concern the effect of these parameters and the vortex intensity on instability characteristics. In particular, in the case of homentropic flow,vortices considered are found to be stable only for sufficiently high intensities and cease to be stable as the intensity weakens. As an example of the situation where unstable normal modes can be excited, the scattering of sound waves by the vortexflow is considered. By simulating this process numerically, we show that the scattered field becomes unstable in the course of time and acquires a typical periodical pattern in polar angle. The characteristics of this instability (viz., the growth rate and the azimuthal phase frequency) coincide with those of the linear analysis very closely. This fact shows conclusively that the instability is caused by the induction of unstable normal modes in the vortex by means of sound irradiation.
17(2005); http://dx.doi.org/10.1063/1.1852575View Description Hide Description
Characteristics and evolution processes of the traveling coherent flow structure in the shear layer of an elevated round jet in crossflow are studied experimentally in an open-loop wind tunnel. Streak pictures of the smoke flow patterns illuminated by the laser-light sheet in the median and horizontal planes are recorded with a high speed digital camera. Time histories of the instantaneous velocity of the vortical flows in the shear layer are digitized by a hot-wire anemometer through a high-speed data acquisition system. By analyzing the streak pictures of the smoke flow visualization, five characteristic flow structures, mixing-layer type vortices, backward-rolling vortices, forward-rolling vortices, swing-induced mushroom vortices, and jet-type vortices, are identified in the shear layer evolving from the up-wind edge of the jet exit. The behaviors and mechanisms of the vortical flow structure in the bent shear layer are prominently distinct in different flow regimes. The frequency characteristics, Strouhal number, power-spectrum density functions, autocorrelation coefficient, as well as the time and length scales of the coherent structure and the Lagrangian integral scales are obtained by processing the measured instantaneous velocity data. The Strouhal number is found to decay exponentially with the increase of the jet-to-crossflow momentum flux ratio. The autocorrelation coefficients provide the information for calculating the statistical time scales of the coherent structure and the integral time scales of turbulencefluctuations. The corresponding length scales of the vortical structure and the integral length scales of turbulence in the shear layer are therefore obtained and discussed.
17(2005); http://dx.doi.org/10.1063/1.1852576View Description Hide Description
The dynamo effect is demonstrated numerically in precession driven flow in a spherical container. Both laminar as well as unstable flows act as dynamos. At low Ekman numbers in unstable flows, the dynamo mechanism relies predominantly on the components of the flow excited by instabilities. All calculations are performed in a frame of reference attached to the boundaries. In this frame, the rotation axis of the fluid executes a periodic motion with a period equal to the rotation period of the boundaries, whereas the magnetic dipole moment undergoes slower variations interrupted by reversals.
17(2005); http://dx.doi.org/10.1063/1.1852574View Description Hide Description
The late-time development of Richtmyer–Meshkov instability is studied in shock tubeexperiments. This investigation makes use of the experimental apparatus and visualization methods utilized in the earlier study of Collins and Jacobs [J. Fluid Mech.464, 113 (2002)] but employs stronger shocks and initial perturbations with shorter wavelengths to obtain much later-time (in the dimensionless sense) images of the single-mode instability. These modifications produce a very detailed look at the evolution of the late-time single-mode instability, revealing the transition and development of turbulence in the vortex cores that eventually results in the disintegration of the laminar vortex structures into small scale features. Amplitude measurements taken from these images are shown to be effectively collapsed when plotted in dimensionless variables defined using the wave number and the initial growth rate. The amplitude measurements are compared with several late-time nonlinear models and solutions. The best agreement is obtained with the model of Sadot et al. [Phys. Rev. Lett.80, 1654 (1998)] which can be slightly improved by modifying the expression for the late-time asymptotic growth rate.
17(2005); http://dx.doi.org/10.1063/1.1863285View Description Hide Description
A discussion is presented on the role of limited conductivity and permittivity on the behavior of electrified jets. Under certain conditions, significant departures with respect to the perfect-conductor limit are to be expected. In addition, an exploration is undertaken concerning the validity of one-dimensional average models in the description of charged jets. To that end, a temporal linear modal stability analysis is carried out of poor-conductor viscousliquid jets flowing relatively to a steady radial electric field. Only axisymmetric perturbations, leading to highest quality aerosols, are considered. A grounded coaxial electrode is located at variable distance. Most available studies in the literature are restricted to the perfect-conductor limit, while the present contribution is an extension to moderate and low electrical conductivity and permittivityjets, in an effort to describe a situation increasingly prevalent in the sector of small-scale free-surface flows. The influence of the electrode distance , a parameter defined as the ratio of the electric relaxation time scale to the capillary time scale, and the relative permittivity on the growth rate has been explored yielding results on the stability spectrum. In addition, arbitrary viscosity and electrification parameters are contemplated. In a wide variety of situations, the perfect-conductor limit provides a good approximation; however, the influence of and on the growth rate and most unstable wavelength cannot be neglected in the general case. An interfacial boundary layer in the axial velocity profile occurs in the low-viscosity limit, but this boundary layer tends to disappear when or are large enough. The use of a one-dimensional (1D) averaged model as an alternative to the 3D approach provides a helpful shortcut and a complementary insight on the nature of the jet’s perturbative behavior. Lowest-order 1D approximations (average model), of widespread application in the literature of electrified jets, are shown to be inaccurate in low-viscosity imperfect-conductor jets.