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

Flickr Twitter iResearch App Facebook

Search Issue | RSS Feeds RSS
Previous Issue Next Issue

Nov 2007

Volume 19, Issue 11, Articles (11xxxx)

Issue Cover Spotlight Figure

Phys. Fluids 19, 114108 (2007); http://dx.doi.org/10.1063/1.2800371 (11 pages)

J. P. Kubitschek and P. D. Weidman
back to top
RSS Feeds
back to top Instability and Transition

The development of asymmetry for oscillatory flow within a tube containing sharp edge periodic baffles

Mingzhi Zheng, Jie Li, M. R. Mackley, and Jianjun Tao

Phys. Fluids 19, 114101 (2007); http://dx.doi.org/10.1063/1.2799553 (15 pages)

Online Publication Date: 5 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
This paper investigates the evolution of asymmetric patterns for oscillatory flow in a baffled tube. A numerical simulation for three-dimensional flows in an axisymmetric geometry was developed and compared with experimental results obtained using particle image velocimetry (PIV). Sharp edged baffles were used for both numerical simulations and experiments. From the numerical simulation, a stability map of the flow symmetry was obtained as a function of Reynold-Strouhal numbers. The simulations show that for all Strouhal numbers, the flow was axisymmetric at Reynolds numbers less than 100 and asymmetric at Reynolds numbers larger than 225. The flow was less stable to asymmetric disturbances at small or large Strouhal numbers when compared to St = 1.0. In particular, the flow in the region St<0.5 and Re>100 was asymmetric. Two mechanisms for vortex instability transition into three dimensions has been identified. At small Strouhal numbers, the primary mechanism is a shear (Kelvin-Helmholtz) instability. At larger Strouhal numbers, the axisymmetry of the flow is broken because of the collision of travelling eddies that have been shed from opposite baffles. The numerical results are in general in agreement quantitatively with the experimental observations and both experiment and simulation assist in understanding the development of unsteadiness in periodic reversing flows.
Show PACS
47.60.-i Flow phenomena in quasi-one-dimensional systems
47.11.-j Computational methods in fluid dynamics
47.32.-y Vortex dynamics; rotating fluids
47.20.-k Flow instabilities
47.32.cd Vortex stability and breakdown

Stability study of the floating zone with respect to the Prandtl number value

Othman Bouizi, Claudine Delcarte, and Guillaume Kasperski

Phys. Fluids 19, 114102 (2007); http://dx.doi.org/10.1063/1.2798810 (14 pages) | Cited 5 times

Online Publication Date: 6 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Thermocapillary convection in a laterally heated liquid bridge is studied numerically, using a Chebyshev spectral method. The stability of the axisymmetric basic state, with respect to 3D perturbations, is characterized over a large range of Prandtl number values (Pr ∊ [10−3,102]), thanks to the choice of a sufficiently sharp regularizing function of the vorticity singularities. Criteria are established to ensure the compatibility of this function with mass conservation, and to choose a correct model of the lateral heat flux in order to reach a spectral convergence of the results. First 3D nonlinear spectral computations are presented.
Show PACS
81.10.Fq Growth from melts; zone melting and refining

Three-dimensional localized coherent structures of surface turbulence. I. Scenarios of two-dimensional–three-dimensional transition

E. A. Demekhin, E. N. Kalaidin, S. Kalliadasis, and S. Yu. Vlaskin

Phys. Fluids 19, 114103 (2007); http://dx.doi.org/10.1063/1.2793148 (15 pages) | Cited 10 times

Online Publication Date: 9 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The evolution of naturally excited disturbances on a thin liquid film falling down an inclined planar substrate undergoes several transitions between different wave regimes starting from two-dimensional (2D, one spatial dimension and height) solitary pulses at small Reynolds numbers to the “surface turbulence” stage for sufficiently large Reynolds numbers where the surface is randomly covered by localized three-dimensional (3D, two spatial dimensions and height) coherent structures. In this study, we analyze for the first time the instability of 2D pulses to 3D disturbances and the transitions of 2D pulses to fully developed 3D waves. The main instability mechanism responsible for the transition to 3D localized patterns is the Rayleigh instability of well-separated (isolated) 2D solitary waves. On the other hand, the physical mechanism for the 3D instability of 2D periodic waves and wave trains of 2D solitary waves sufficiently close to each other is not related to the Rayleigh instability but is due to wave-wave interaction and mass exchange between neighboring waves. These instabilities are characteristic of small inclination angles but under special conditions can be realized also for the vertical flow.
Show PACS
68.15.+e Liquid thin films
47.27.De Coherent structures
47.27.Cn Transition to turbulence
47.35.Fg Solitary waves
47.20.Ib Instability of boundary layers; separation
47.27.Jv High-Reynolds-number turbulence

Three-dimensional localized coherent structures of surface turbulence. II. Λ solitons

E. A. Demekhin, E. N. Kalaidin, S. Kalliadasis, and S. Yu. Vlaskin

Phys. Fluids 19, 114104 (2007); http://dx.doi.org/10.1063/1.2793149 (15 pages) | Cited 8 times

Online Publication Date: 9 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We numerically construct Λ solitons as a function of the generalized Reynolds number δ. The numerical scheme is based on an impulse response analysis in which the nonlinear hump region is replaced with a Dirac delta function. We also examine the linear stability of Λ solitons with respect to three-dimensional disturbances. It is shown that the operator of the linearized system has both a discrete and a continuous spectrum. The discrete spectrum is always stable, while the continuous spectrum can be destabilized leading to a convective instability of Λ solitons. We demonstrate that the region of existence and stability of Λ solitons is 0.054<δ<0.51.
Show PACS
47.35.Fg Solitary waves
68.15.+e Liquid thin films
05.45.Yv Solitons
47.27.N- Wall-bounded shear flow turbulence
47.20.Ib Instability of boundary layers; separation

Spatially convective global modes in a boundary layer

Frédéric Alizard and Jean-Christophe Robinet

Phys. Fluids 19, 114105 (2007); http://dx.doi.org/10.1063/1.2804958 (12 pages) | Cited 9 times

Online Publication Date: 12 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The linear stability of a weakly nonparallel flow, the case of a flat plate boundary layer, is revisited by a linear global stability approach where the two spatial directions are taken as inhomogeneous, leading to a fully nonparallel stability method. The resulting discrete eigenvalues obtained by the fully nonparallel approach seem to be in agreement with classical Tollmien–Schlichting waves. Then the different modes are compared with classical linear stability approach and weakly nonparallel method based on linear parabolized stability equations (PSEs). It is illustrated that the nonparallel correction provided by the linear global stability approach is well matched by linear PSE. Furthermore, physical interpretation of these spatio-temporal global modes is given where a real pulsation, which has more physical interest, is considered. In particular the use of a Gaster transformation and the pseudospectrum illustrate the local and global properties of these Tollmien–Schlichting modes. Finally, the contribution of different components of global modes (normal and streamwise) in the transient amplifying behavior associated with the convectively unstable boundary layer is analyzed and compared with a classical steepest descent method. Then, a discussion of an equivalent of the continuous branch is given.
Show PACS
47.20.Ib Instability of boundary layers; separation
47.35.-i Hydrodynamic waves
47.27.nb Boundary layer turbulence

Oblique axisymmetric stagnation flows in magnetohydrodynamics

M. Amaouche, F. Naït-Bouda, and H. Sadat

Phys. Fluids 19, 114106 (2007); http://dx.doi.org/10.1063/1.2804957 (7 pages) | Cited 3 times

Online Publication Date: 13 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The present paper deals with the influence of the Lorentz force associated with an applied radial magnetic field on the axisymmetric stagnation flow impinging obliquely onto a uniformly rotating circular cylinder. It is found that the boundary layer flow is described by an exact solution of the Navier–Stokes equations when Hall effects are ignored. The stability of this basic solution is then considered in the framework of Görtler–Hammerlin assumption according to which linear disturbances inherit the underlying symmetry of the basic flow. The resulting eigenvalue problem is solved numerically by means of a pseudospectral method using Laguerre’s polynomials. The scheme is specifically designed to solve boundary layer equations. The numerical experiments indicate that the effect of cylinder rotation is to reduce the stability of the basic flow and most importantly the magnetic field acts to either increase or decrease it, depending on whether the rotation rate is smaller or greater than some critical value that changes with the Hartmann number. At criticality, the basic flow undergoes Hopf bifurcations leading to branching off solutions in the form of azimuthally travelling waves. In the case of axisymmetric disturbances the bifurcation remains of Hopf type provided that the Hartmann number is small enough, a saddle-node bifurcation is encountered in the opposite case.
Show PACS
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking
47.65.-d Magnetohydrodynamics and electrohydrodynamics
47.10.ad Navier-Stokes equations

Bifurcation and stability of near-critical compressible swirling flows

Z. Rusak, J. J. Choi, and J.-H. Lee

Phys. Fluids 19, 114107 (2007); http://dx.doi.org/10.1063/1.2801508 (13 pages) | Cited 3 times

Online Publication Date: 14 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
The bifurcation and global nonlinear stability of near-critical states of a compressible and axisymmetric swirling flow of a perfect gas in a finite-length straight, circular pipe is studied. This work extends the bifurcation and stability analyses of Wang and Rusak [Phys. Fluids 8, 1007 (1996) ; Wang and Rusak Phys. Fluids8, 1017 (1996) ] to include the influence of Mach number on the flow dynamics. The first- and second-order equations of motion for the evolution of small axially symmetric perturbations on a base columnar state are developed. These equations are reduced to an eigenvalue problem for the perturbation shape function and critical swirl ratio and a model ordinary differential equation for the nonlinear evolution of the perturbations’ amplitude as function of swirl level and Mach number. It is found that noncolumnar equilibrium states bifurcate from the branch of the base columnar equilibrium states at the critical swirl ratio of a compressible vortex flow in the form of a transcritical bifurcation, where both the critical swirl and the bifurcation slope ratio are functions of Mach number. It is also shown that this critical swirl ratio is a point of exchange of global stability for both the columnar and noncolumnar states as the swirl ratio increases across this critical level. When the swirl ratio of the incoming flow is below the critical level the columnar states have an asymptotically decaying mode of perturbation whereas the noncolumnar states are unstable. On the other hand, when the swirl ratio of the incoming flow is greater than the critical level, the columnar states are unstable whereas the noncolumnar states have an asymptotically decaying mode of disturbance. The effect of Mach number on the bifurcation behavior and on the stability characteristics of the various states is presented. The relationship between the present results and the breakdown of a compressible vortex flow is discussed.
Show PACS
47.32.cd Vortex stability and breakdown
47.40.-x Compressible flows; shock waves
47.60.-i Flow phenomena in quasi-one-dimensional systems
02.30.Hq Ordinary differential equations
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking

Helical instability of a rotating viscous liquid jet

J. P. Kubitschek and P. D. Weidman

Phys. Fluids 19, 114108 (2007); http://dx.doi.org/10.1063/1.2800371 (11 pages) | Cited 4 times

Online Publication Date: 15 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Vertical rotating viscous liquid jet experiments show a clear preference for helical instabilities that evolve from initially planar disturbances at large rotation rates for fixed fluid properties. The laboratory setup for the experiments described herein was chosen as the nearest earth-based equivalent to a uniformly rotating viscous liquid column in the absence of gravity. In the ideal situation with stress-free boundaries, the preferred modes of linear temporal instability are theoretically known over the entire physical domain spanned by the Hocking parameter L = γ/ρa3Ω2 and the rotational Reynolds number Re = a2Ω/ν, where a is the column radius, Ω is its uniform angular velocity, and ρ, ν, and γ are, respectively, the fluid density, kinematic viscosity, and surface tension. The theoretical results show that instability in L-Re parameter space is dominated by three mode types: The axisymmetric mode, the n ≥ 2 planar modes, and the first n = 1 spiral mode. Experiments reveal that, in the L-Re region for which the uniformly rotating liquid column is dominated by planar modes of instability, the rotating liquid jet spontaneously gives rise to planar disturbances of mode n ≥ 2 that rapidly evolve into helical instabilities. However, these observed instabilities are not the spiral normal modes that exist for n ≥ 1 as posited in linear stability theory. In spite of obvious fundamental differences between the rotating liquid jet and the uniformly rotating liquid column, some remarkable similarities associated with initial growth rates, angular frequencies, and mode transitions between the two systems are found.
Show PACS
47.20.-k Flow instabilities
47.32.Ef Rotating and swirling flows

Rotating disk flow stability in electrochemical cells: Effect of the transport of a chemical species

N. Mangiavacchi, J. Pontes, O. E. Barcia, O. R. Mattos, and B. Tribollet

Phys. Fluids 19, 114109 (2007); http://dx.doi.org/10.1063/1.2805844 (15 pages) | Cited 2 times

Online Publication Date: 20 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
We consider the stability of rotating disk flow coupled, through the fluid viscosity, to the mass concentration field of a chemical species. This configuration refers to an electrochemical cell where the working electrode consists of an iron rotating rod, which is dissolved in the 1 M H2SO4 solution of the electrolyte. Polarization curves obtained in these cells present a current instability at the beginning of the region where the current is controlled by the mass transport. The instability appears at a certain value of the applied potential and is suppressed beyond another value. Dissolution of the electrode gives rise to a thin concentration boundary layer, due to a Schmidt number Sc = 2000 of the setup. This boundary layer, together with the potential applied to the electrode, leads to an increase in the fluid viscosity and to a decrease in the diffusion coefficient, both affecting the chemical species field. Since the current is proportional to the normal derivative of the species concentration at the interface, an instability of the coupled fields at sufficiently low Reynolds numbers would result in a current instability. This work deals with the question of whether the coupling reduces the critical Reynolds number to values comparable to those attained in experimental setups, and if a possible field instability would be large enough to drive a detectable current instability. A phenomenological law is assumed, relating the fluid viscosity to the concentration of the chemical species. Parameters appearing in this law are evaluated on the basis of experimental electrochemical data. The steady-state solution is obtained by solving the coupled hydrodynamic and mass concentration equations. A temporal stability analysis is made, showing that small variations in the fluid viscosity significantly affect the stability of the flow. The analysis reveals the existence of a new unstable region, not found in the case of constant viscosity fluids. We call modes in this new unstable region chemical modes, in contrast with the hydrodynamic modes, which are amplified in the case of constant viscosity fluids. The chemical modes are destabilized at much lower Reynolds numbers than the hydrodynamic ones and close to values attained in electrochemical setups, in the most relevant cases. Hydrodynamic modes are strongly affected by the coupling in three aspects: the critical Reynolds number of this region is of order of 50% smaller than in the case of constant viscosity fluids, the unstable region is enlarged to a wider range of wave numbers, and the rate of growth of unstable modes is 30% larger for comparable Reynolds numbers. Concentration eigenmodes in the new unstable region show a combination of properties including rate of growth, amplitude higher than the amplitude of the hydrodynamic variables, and high normal derivative at the interface, sufficiently strong to drive detectable current oscillations. The destabilizing effect of a higher interface viscosity attains a maximum when the ratio between the interface and the bulk viscosities, ν0/ν, takes a value close to 1.5. Below this value, as the viscosity stratification diminishes and the condition of uniform viscosity is restored, the unstable region of chemical modes moves to higher Reynolds and eventually disappears. Conversely, as the ratio ν0/ν increases beyond the value 1.5, the new unstable region also collapses and the neutral curve of hydrodynamic modes tends to the one of constant viscosity fluids. A sustained increase of the interfacial viscosity and the high Schmidt number results in an in facto field discontinuity with a thin high viscosity layer at the interface and restored constant viscous hydrodynamic boundary layer. A link between the current instability and the stability of the coupled fields may be inferred from the present analysis.
Show PACS
47.20.Ib Instability of boundary layers; separation
47.32.Ef Rotating and swirling flows
47.85.Dh Hydrodynamics, hydraulics, hydrostatics

Stability analysis of a class of unsteady nonparallel incompressible flows via separation of variables

Georgy I. Burde, Ildar Sh. Nasibullayev, and Alexander Zhalij

Phys. Fluids 19, 114110 (2007); http://dx.doi.org/10.1063/1.2814296 (14 pages) | Cited 1 time

Online Publication Date: 26 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Stability of some unsteady three-dimensional flows (exact solutions of the viscous incompressible Navier–Stokes equations in cylindrical coordinates) is studied via separation of variables in the linearized equations for the flow perturbations. The flows in an expanding rotating porous cylinder and in a gap between two coaxial rotating cylinders are considered. Converting the stability equations to the new variables allows perturbation forms (counterparts of normal modes of the steady state parallel flow stability problem) such that the linear stability problems are exactly reduced to eigenvalue problems of ordinary differential equations. The eigenvalue problems are solved numerically with the help of the spectral collocation method based on Chebyshev polynomials. The results showing dependence of the stability threshold on the parameters of the problems and a spatial structure of the unstable perturbation modes are presented. For some classes of perturbations, exact analytical solutions of the eigenvalue problems are available. A combination of analytical and numerical solutions can provide useful testing for numerical methods used in the hydrodynamic stability studies. It may also provide a basis for a well-grounded discussion of some problematic points of (numerical) stability analysis. In particular, in the present paper, a problem of formulation of the boundary conditions for perturbations at the axis r = 0 is discussed on the basis of the solutions obtained.
Show PACS
47.32.Ef Rotating and swirling flows
47.20.-k Flow instabilities

Equilibrium shapes and stability of a liquid film subjected to a nonuniform electric field

Hak Koon Yeoh, Qi Xu, and Osman A. Basaran

Phys. Fluids 19, 114111 (2007); http://dx.doi.org/10.1063/1.2798806 (22 pages) | Cited 7 times

Online Publication Date: 30 November 2007

Full Text: Read Online (HTML) | Download PDF

Show Abstract
Stresses induced by a spatially nonuniform electric field acting on an initially flat fluid-fluid interface can (i) be exploited beneficially to pattern polymer microstructures without the use of resists, exposure, development, and etching, but (ii) cause undesirable nonuniformity in film thickness in precision coating processes. The equilibrium shape of an interface separating a liquid film from an ambient fluid subjected to a uniform electric field is flat so long as the field strength is below a critical value. A nonuniform electric field, however, results in the deformation of the interface no matter how small its strength, an important difference which previous theoretical studies have not addressed satisfactorily. Hence, whereas under a uniform field loss of stability occurs via a bifurcation from the flat film solution, under a nonuniform field destabilization may occur at a turning point at which the film profile already exhibits a finite-amplitude deformation. This deficiency in understanding is remedied here by analyzing a model problem in which a gas overlying a perfect dielectric liquid film is sandwiched between two electrodes wherein the top electrode is grounded and the electric potential of the bottom electrode varies sinusoidally with distance measured along it. The equilibrium shapes and stability of the liquid-gas interface are determined directly in the present work by simultaneously solving the augmented Young-Laplace equation governing the shape of the free surface and the Laplace equation governing electric potentials theoretically by domain perturbation analysis and numerically by finite element analysis. For small nonuniformities in the electric field, analytical solutions are reported for the profile of the free surface. The computational predictions are shown to be in excellent accord with these small-deformation results. Moreover, computations are used to extend the investigations into the nonlinear regime where nonuniformities in the electric field and deformations of the free surface are large, and loss of stability may occur. The variation of the equilibrium shapes and the limits of stability with the governing dimensionless groups are investigated thoroughly. It is shown that the rich response exhibited by the system can be rationalized by interrogating the computed solutions and scrutinizing the balance of stresses due to the normal component of the electric field, which are destabilizing, and those due to its tangential component, which are stabilizing.
Show PACS
68.15.+e Liquid thin films
77.55.-g Dielectric thin films
68.03.Hj Liquid surface structure: measurements and simulations
62.10.+s Mechanical properties of liquids
Close
Google Calendar
ADVERTISEMENT

close