VORTEX BRAJWOWN: A TWO-STAGE TRANSITION M.P. Escudier and J.J. Keller Brown Boveri Research Centre CH-5405 Baden, Switzerland r 1

Similar documents
Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction

Article published by EDP Sciences and available at or

In which of the following scenarios is applying the following form of Bernoulli s equation: steady, inviscid, uniform stream of water. Ma = 0.

Evolution and transition mechanisms of internal swirling flows with tangential entry

D. C. Lewellen and W. S. Lewellen West Virginia University, Morgantown, WV.

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015

Primary oscillatory instability in low-aspect-ratio rotating disk cylinder system (rotor stator cavity)

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation

Detailed Outline, M E 521: Foundations of Fluid Mechanics I

Central recirculation zones and instability waves in internal swirling flows with an annular entry

Axisymmetric vortex breakdown Part 1. Confined swirling flow

Basic Fluid Mechanics

SPECTRAL ELEMENT STABILITY ANALYSIS OF VORTICAL FLOWS

6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s

SIMULATION OF PRECESSION IN AXISYMMETRIC SUDDEN EXPANSION FLOWS

F11AE1 1. C = ρν r r. r u z r

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

Intermittency in spiral Poiseuille flow

Iran University of Science & Technology School of Mechanical Engineering Advance Fluid Mechanics

5. Secondary Current and Spiral Flow

Friction Factors and Drag Coefficients

Figure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m

THREE-DIMENSIONAL VERSUS TWO-DIMENSIONAL AXISYMMETRIC ANALYSIS FOR DECELERATED SWIRLING FLOWS

Influence of outlet geometry on strongly swirling turbulent flow through a circular tube

Fundamentals of Fluid Mechanics

The onset of absolute instability of rotating Hagen Poiseuille flow: A spatial stability analysis

Analysis of the Swirling Flow Downstream a Francis Turbine Runner

New features of swirling jets

Aerothermodynamics of high speed flows

SHORT WAVE INSTABILITIES OF COUNTER-ROTATING BATCHELOR VORTEX PAIRS

UNIT II Real fluids. FMM / KRG / MECH / NPRCET Page 78. Laminar and turbulent flow

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation

Part A: 1 pts each, 10 pts total, no partial credit.

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

Chapter IX BREAKDOWN OF SLENDER VORTICES. W. ALTHAUS, CH. BRUCKER, M. WEIMER Aerodynamisches Institut RWTH Aachen GERMANY

Fig. 1. The coordinate system and the directions of the velocities.

A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER

VARIED FLOW IN OPEN CHANNELS

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition

Capillary Collapse and Rupture. Abstract

L. Pratt and J. Whitehead 6/25/ Parabolic Bottom

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes

Effects of Free-Stream Vorticity on the Blasius Boundary Layer

Fluid Dynamics Problems M.Sc Mathematics-Second Semester Dr. Dinesh Khattar-K.M.College

Shock and Expansion Waves

Strongly nonlinear long gravity waves in uniform shear flows

Structure and mixing of a transverse jet in incompressible flow

Notes #4a MAE 533, Fluid Mechanics

Mathematical Modeling of Peristaltic Flow of Chyme in Small Intestine

Some Basic Plane Potential Flows

2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;

FLOW MEASUREMENT. INC 102 Fundamental of Instrumentation and Process Control 2/2560

Self-Excited Vibration in Hydraulic Ball Check Valve

Lecture-2. One-dimensional Compressible Fluid Flow in Variable Area

Helical Coil Flow: a Case Study

Thin airfoil theory. Chapter Compressible potential flow The full potential equation

Evolution of the pdf of a high Schmidt number passive scalar in a plane wake

Turbulence Modeling Applied to Flow over a Hydraulic Ball Check Valve

Differential criterion of a bubble collapse in viscous liquids

To study the motion of a perfect gas, the conservation equations of continuity

OUTLINE FOR Chapter 3

Acoustics flow analysis in circular duct using sound intensity and dynamic mode decomposition

DEPOSITION FROM A CURVED SHALLOW FLOW TREATED AS A MOVING BOUNDARY OF CONSTANT FORM. R.A. Wooding

On a variational principle for Beltrami flows

Experimental and Numerical Investigation of Flow over a Cylinder at Reynolds Number 10 5

THE COMPUTATION OF FLUID-INDUCED FORCES ON CENTRIFUGAL IMPELLERS ROTATING AND WHIRLING IN A VOLUTE CASING. N s Specific speed (Q 1 2 =(gh) 3 4 )

General introduction to Hydrodynamic Instabilities

Velocity and temperature measurements in a large-scale Rayleigh-Bénard experiment using LDA and micro thermistors

Contents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1

CE Final Exam. December 12, Name. Student I.D.

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

The Rossby-inertia-buoyancy instability in baroclinic vortices

Given the water behaves as shown above, which direction will the cylinder rotate?

Offshore Hydromechanics Module 1

Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus

CONVECTIVE HEAT TRANSFER

Experiments on capillary-gravity waves of solitary type on deep water

Time dependent singularities in incompressible potential flow

Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows

STUDY OF FLOW & HEAT TRANSFER IN PLATE FIN HEAT EXCHANGER AT VARYING REYNOLD S NO

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS ON CHOKING PHENOMENA OF AXISYMMETRIC CONVERGENT NOZZLE FLOW

HIGH SPEED GAS DYNAMICS HINCHEY

Numerical study of 2D heat transfer in a scraped surface heat exchanger

Temperature, Pressure and Velocity measurements on the Ranque-Hilsch Vortex Tube

The effect of disturbances on the flows under a sluice gate and past an inclined plate

Rotational viscometers

Experimental and numerical investigation of the competition between axisymmetric time-periodic modes in an enclosed swirling flow

Visualization of flow pattern over or around immersed objects in open channel flow.

3 Generation and diffusion of vorticity

Transactions on Engineering Sciences vol 5, 1994 WIT Press, ISSN

Turbulence Laboratory

Breakdown in a boundary layer exposed to free-stream turbulence

Dynamics of Transient Liquid Injection:

In this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor

APPH 4200 Physics of Fluids

6. Basic basic equations I ( )

Steady waves in compressible flow

Application of COMSOL Multiphysics in Transport Phenomena Educational Processes

Transcription:

. 5 VORTEX BRAJWOWN: A TWO-STAGE TRANSITION M.P. Escudier and J.J. Keller Brown Boveri Research Centre CH-5405 Baden, Switzerland r 1 0 I SU R 0t is demonstrated that a large-scale isentropic transition between conjugate ' swirling flow states can occur with no change in the flow force and that both flow states are supercritical. It is argued that such a transition represents the first stage of vortex breakdown in a tube, the second stage being a non-isentropic transition in the nature of a hydraulic jump to the downstream subcritical state. The intermediate (supercritical) state consists of a zone of stagnant fluid surrounded by a region of potential flow. These two zones are separated by a layer of rotational fluid originating in the upstream vortex core. An outline is given of the analysis for an upstream flow modelled as a Rankine vortex. It is found that for any ratio of core-to-tube raii, breakdown (i.e. the first transition) occurs for a unique value of the-wirl number/,/r-u In the limiting case of an infitesimally small core, the value is ex compared witfi he critical value 2.405. It is argued that this limit cannot represent free breakdown, which in consequence must have a different character from the tube-flow breakdowns generally observed. I NTRODUCTlON Previous efforts to explain the phenomenon of vortex breakdown have been centred upon the concept of a single transition, either in the sense of Benjamin's [1] conjugate-state analysis or stability theories such as those of Ludwieg 12] and, most recently, Leibovich 131. For the most part, these efforts have been limited to weak transitions. The experimental evidence reported by Harvey [4], Sarpkaya (5] and others (6,7], in contrast, shows breakdown to be invariably a strong perturbation of the flow, as has been emphasized by, for example, Hall [8] and Leibovich [9]. We present a new and simple approach to the problem of vortex breakdown in tube flow which yields a simple breakdown criterion and is consistent with the observed characteristics of the phenomenon. complete vortex breakdown. ~tube wall I!i: frrotatione -- srbrcrrtical I f stlowaeisat I I I- Svortexý _ I. -, I - S~~stagnation zone _-- I first supercriticol second supercriticai I Jsubcriticot f low state / I f flow state!_._1flow state first (isentropic) second transition transition (hydraulic jump) Figure 1. Schematic diagram illustrating the idealized process of vortex breakdown in a cylindrical tube. The key features of the process of vortex breakdown in a tube are illustrated in figure 1. We propose that the transition may involve two stages, the first isentropic from the initial supercritical state to an intermediate state which is also supercritical. The second stage of the transition, to the downstream subcritical state, is nonisentropic, much like a hydraulic jump or shock wave. The crucial new idea here is that the first transition is both isentropic and also involves no change in the flow force even for a large-scale transition. The latter possibility has been overlooked until recently [10], most significantly by Benjamin (1]. A striking consequence of the analy-

25-2 sis is that for a given upstream vortex structure, breakdown occurs only for a unique value of the swirl number r/nr U (an inverse Rossby number): e.g. 42 for a Rankine vortex with an infinitesimally *mill core radius r. Another prediction is that breakdown is also possible for a purely irrotational swirlfng flow: the example considered is that of flow in an annulus. complete details of the analysis are given in two recent papers [10,11] and only the barest essentials are excerpted here to emphasize the principal assumptions and results. we also present the results of a number of experiments in support of the analysis and speculate on the nature of free vortex breakdown and on the real and apparent differences between the axisymmetric and spiral forms of breakdown. VORTEX BREAKDOWN IN A CYLINDRICAL TUBE - ANALYSIS A convenient model for the upstream flow state (1) (see figure 1) is the Rankine vortex: v = rr/2nr 2 if 0 < r c r/2nr if rc < r < rt (1)1 w = constant = U where v and w are the swirl and axial components of velocity, r is the core radius. r the tube radius and r the circulation. The intermediate state (5) is "-so assumed to b cylindrical with an inner stagnation region of radius r= separated from the outer potential flow by a layer r < r < rb of rotational fluid o iginating in the upstream vortex core. The flow between +,hese two states is assumed to be steady, incompressible, axisymmetric and inviscid. The equation governing the intermediate flow state is then [12] where ro is d2 94-1 d ff - 1Y= - k2* if ra < r < rb (2) r 0 if rb < r < rt the departure of the stream function from its upstream form and k r/nr 2 U. (3) Specification of the intermediate state is completed by introducing appropriate matching conditions at r, r and rt. i.e. constancy of the stream function on stream surfaces, continuity of tre v~locity across the interface r - rh and zero velocity at r = r. The solution of equation (2) and calculation of the distributions of v, w and the Itatic pressure p (from dp/dr = pv 2 /r) is then straightforward (see 110] - [12]). To answer the question which pair of conjugate flow states corresponding to the preceding analysis is physically realistic, we must consider the momentum, equation. We define first the flow force S by S = 2nfr (pw2 + p)rdr (4) where p is the fluid density. Using the results already obtained, the flow force difference between the first and second flow states may be shown to be A S ru2k r 2[ + 1 1r 4-ra4 4 _ r1r 3r i l l b = [3-r ) + + n(-)]. _r (5) Since we are considering flow in a cylindrical tube, with no neans of applying ar -- ternal force to the flow, the difference AS =_ 0. The values of krc (= r/nr U) for..iich this identity is satisfied are plotted in figure 2 as a function of rc/rt (ogether with the corresponding critical values of krc obtained from [1] 1 2 ;.rc Jl(krc) Jo(krc) = (rt/rc)=_l -1 (6) It is found that th, upstream flow (for AS = 0) is always supercritical (i.e. kr less than the critical value), as is the intermediate flow state which can be demonstraied on the basis of a variational principle [1,10]. As mentioned in the INTRODUCTION, this possibility of a non-trivial supercritical-supercritical transition was only recently recognized (10] and is suggested to be the key to understanding vortex breakdown. For the Rankine vortex, it turns out that for any value of r /r the first transition occurs for a unique value of krc. More generally it may be shown hl0] tl.at for given distributions of w(r) and v(r) the speed at which the vortex breakdown structure moves is unique. The subsequent transition to the downstream flow state muot be dissipative, and can be treated essentially as a hydraulic jump.

25-3 krc 3.852 5 2.405-2 rc / rt 0 R I I I I O 0 0.5 1 Figure 2. Swirl parameter kr for breakdown (solid curve) and critical-flow (broken curve) conditions is a function of normalized core radius. A remarkable result of the analysis is that a non-trivial transition is predicted even for the asymptotic situation rc/rt * 0. in this case krc M "2 krc = 2.405 for breakdown for the critical state At first sight it might appear that the asymptote r_/r. - 0 is relevant to the case of a free vortex (r - P). However, an asymptotic analysib, which was the subject of (10], shows that thik is not the case. In particular, it is found that for r /r * 0, r_/r. 4 0 whereas for a free vortex it must be the case that rj/r, assumes aqailae different from zero. It has to be concluded that the free vortex reptestnta a second type of transition compared to that which we have considered, and that this is ". direct supercritical - subcritical transition necessarily involving dissipation. A corollary is that investiga ions of vortex breakdown in tubes are not directly representative of free-vortex breakdown. Within the small-core approximation, the ideas contained here have been extended to the analysis of general vortex flow in a tube, including the possibility of area changes. In many circumstances it is possible and convenient to relate the actual flow to a fictitious reference Rankine vortex. An interesting result for flows with a strong overshoot in the axial velocity, as is produced by a strong area reduction, is that the critical and breakdown values of kr are proportional. This is suggested to be the main reason why Escudier and Zehnder were able to make quite accurate predictions of breakdown on the basis of the criticality condition. i i VORTEX BREAKDOWN IN A CYLINDRICAL TUBE - OBSERVATIONS The two examples of axisymmetric vortex breakdown shown in figure 3 are typical of those reported by Harvey [4), Sarpkaya (5] and others (6,7]. Here advantage has been taken of the laser-induced fluorescence visualisation technique to reveal the inner structure of a breakdown bubble. Fluorescein dye has been injected on the tube axis into the upstream flow. The flow was illuminated by a rapidly oscillating Argon-ion laser beam sweeping through a diametral plane of the tube (further details of the experiment are given in [10] and [13]). Two features may be observed from figure 3 which support the interpretation of vortex breakdown proposed here. First, the fluid entering the bubble in each case clearly emanates from a reqion much smaller in radius than that of A the bubble itself (and also much smaller in radius than the core), and the latter must J therefore be a zone of essentially stagnant fluid. secondly, the smooth appearance of -Ii

"t e " bul~b e c I 1Wen, tot' :: " 'hen t ' :nt t 'r;! 1: t1ll 2:t 'td. -;t:]]. :; tici l ),. ; I : - tin' K i :-1 1 :t.:; :1t:v1:}.eI- dd tt l t i tutab J t I. -:, auld 11d Io jul I 1 up +at, :.gi't. 1L)II - lo II. Id d;:~~l:' 11L : i: st: ' -;L, I' ti.t: ti. ' I k Redk t- it dh. iýxe : :.. a) Ur. 575 (b) U 1 t,.= 1650 Figure 3. Exa.sples of vortex breakdown in a cylirdrical tube. (Flow from left to right) The internal structure which may be inferred from flow visual;sation is confirmed by laser Doppler anemometry measurem~ents: axial velocity profiles w(r:x) together wilth the corresponding streamline map for a low Reynolds number breakdown bubble are show i, figure 4. It may be noted that a consequence of the subcritical nature of the flow downstream of such a bubble is that the detailed structurc of the bubble itself and the downstream flow are strongly influenced by the downstream geometry arid conditions. The problem of masking of the intermediate flow state, tintioned above, may be overcome by injecting air into the bubble as in the situation shown in figure 5. The smooth character of the first transition is again evident as is the more disturbed nature of the second (the hydraulic junj,). The flow conditions ale those under whlich a normal breakdown would occuir in the absence of air inlectron.

25-5 I 10 l'i/ I'0 ~~ 20 2 _30 I0 ~iizi~i~~ 50 10050 50I W.UJJ~t 4JJ -'-70 S 80 980 90I 100 100 1000 ~110 / l! } I120 r(mm) 255_[52 ml Figure 4. -4 LA - Measured axill velocity profiles w(r;n) and streamline map *(r x) for vortex breakdown in a cylindrical tube with~ Ur t/v =265. Figure S. An aie-filled vortex-breakdown bubble. (Flow from left to right. Arrows indicate transition locations)

25-6 BREAKDOWN OF POTENTIAL FLOW IN AN ANNULUS ANALYSIS We consider here a swirling potential flow confined between concentric cylinders. The upstream flow etate is defined by v = r/2nr WU rc < r rt (8) i.e. we car. imagine the viscous core in figure 1 to be replaced by a solid rod. The second or intermediate flow state is assumed to consist of an annulus of stagnant fluid of outer radius rb surrounded by a potential vortex i.e. v = r/2nr and w = W if rb < r <r (9) v = 0 and w = 0 if rc r < rb rh/r kr 12t 0.5 10 8_ 0 0 0 0.5 r /rt 6k--6 4 2 Figure 6. Normalized bubble radius and breakdown swirl parameter krc versus normalized inner radius for annular breakdown. The assumption of an isentropic transition between these two flow states then leads to (101 ( 1 b( ] (10) ~ rt-rjt

whilst for the flow-force difference we have AS = Lcu - k,) 1- l..()? 42 (Sc)ý,ýf(ES) I b -t _1b b b b Setting AS = 0 as before leads to krc and r ar as functions ofa,r'r as plotted in figure 6. In this case the criticality condat ontfur the downstreamtfi8w 1141 is 'k t~ Ci(i2) and it is again found that thils flow state is always supetctitical. BREAKDOWN OF POTENTIAL FLOW IN AN ANNULUS - OBSERVATIONS Fio Vii~u~lrur".Vn conflrms the occurence of large-scale transitions for swirlling flow in an annulus. In the first example (figule 7) multiple breakdowns are evident. Althouah the viscous cote has been "replaced" by a,olld rod, viscous influences are undoubtedly responsible for producing the conditions which allow successive breakdowns- The introduction of ail into the breakdown legion again produces a breakdown bubble similar to that predicted (figure 8). It. is found that unless conditions closely match those corresponding to the analysis, air introduced into the flow either penetrates far upstream (k:c too ldrgc or is swept away (kre too small). Figure 7. Successive annular breakdowns. (Flow from left to right) "--.,,.. -. $... -.. c i - ~i Figure 8. An air-filled annular breakdown. (Flow from left to right) CONCLUDING REMARKS he concept proposed here for the explanation of vortex breakdown has been worked out in dctail using the Rankine vortex as a model for the upstream flow. The advantage of this formulation is that the equation for the stream function (2) is then linear both within the rotational layer (r, r < r ) and also within the surrounding potential flow (r r < r ). As has been in~icated, More complex flow situations can be analysed with reverence to a fictitious reference Rankine vortex. Such an analysis has been carried out in [101 within the small-core approximation, a restriction which can also be dis- SI! "!' "!

"2S4 pensed with as shown in [111. Any more general analysis would evidently require numerical integration of the equation governing the stream function * (the Long [151-Squire [16] equation) for a cylindrical flow 2. 1 =r2 dh _ dk (3 wherein H is the dynamic head and K = rv. Such an exercise would undoubtedly provide more accurate results for particular flows, but would be unlikely to contribute significantly to further understanding of the breakdown phenomenon. We have already speculated that free-vortex breakdown represents a different type of transition to that for tube flow. By way of further speculation we suggest that there is no fundamental difference between the so-called spiral type of breakdown and the axisymmetric bubble type we consider here. The sniral character is suggested to be a consequence of azimuthal instability leading to roll up and detachment of th: shear layer surrounding a bubble. This detachment process would then lead to a precessive motion of the near-stagnant interior fluid. Strong evidence in support of this point of view comes from the experiments of Escudier and Zehnder [7] who found the correlation of experimental breakdown data in terms of the simple criterion Fz/UvL = constant to be independent of the breakdown's appearance. It is also well known [8] that the spiral and axisymmetric breakdown forms are found in the same continuous range and are even interchangable under appropriate conditions. A further point is that were the spiral form inherently ron-axisymnettic, the spiral would have to remain stationary, as is observed for the non-axisymmetric standing-wave patterns for hollow-core vortices [14]. REFERENCES I. T.B.Benjamin. Theory of the vortex breakdown phenomenon. J. Fluid Mech. 14, 593-629 (1962). 2. H. Ludwieg. Zur Erklirung der Instabilitht der diber angestelliten D6ltaflugeln auftretenden freien Wirbelkerne. Z. Flugwiss. 10, 242-249 (1962). 3. S. Leibovich. Wave propagation, instability, and breakdown of vortices. In "Vortex Flows", 50-67, Ed. Hornung and MUller, Vieweg (1982). 4. J.K. Harvey. Some observations of the vortex breakdown phenomenon. J. Fluid Mach. 14, 585-592 (1962). S. T. Sarpkayd. On stationary and travelling vortex breakdowns. J. Fluid Mech. 45, 545-559 (1971). 6. J.H. Faler and S. Leibovich. Disrupted states of vortex flows and vortex breakdown. Phys. Fluids 20, 1385-1400 (1977). 7. M.P. Escudier and N. Zehnder. Vortex flow regimes. J. Fluid Mech. 115, 105-121 (1982). 8. M.G. Hall. Vortex breakdown. Ann. Rev. Fluid Mech. 4, 195-218 (1972). 9. S. Leibovich. The structure of vortex breakdown. Ann. Rev. Fluid Mech. 10, 221-246 (1978). 10. J.J. Keller, M.P. Escudier and J. Bornstein. The physics of vortex breakdown. (1982). Submitted for publication. 11. J.J. Keller. Breakdown of a Rankine vortex. (1983). Submitted for publication. 12. G.K. Batchelor. An introduction to fluid dynamics. Cambiidge University Press. (1967). 13. J. Bornstein and M.P. Escudier. LDA measurements within a vortex.-breakdown bubble. Proc. Int. Symp. on Applications of Laser Doppler Anemometry. Lisbon, (1982). 14. J.J. Keller and M.P. Escudier. Theory and observations of wave on hollow-core vortices. J. Fluid Mech. 99, 495-511 (1980). 15. R.R. Long. Steady motion around a symmetrical obstacle moving along the axis of a rotating fluid. J. Meteor 10, 197-203 (1953). 16. H.B. Squire. Rotating fluids. In "Surveys in mechanics". 139-161,. Ed. Batchelor and Davies, Cambridge University Press (1956).