Shear instabilities. Chapter Energetics of shear instabilities
|
|
- Cameron Dean
- 6 years ago
- Views:
Transcription
1 Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever two adjacent fluid parcels move in parallel directions, but at different velocities. The shear is defined as the amplitude of the local velocity gradient perpendicular to the motion. Shear flows are ubiquitous in nature and can occur on any scale. Flows pumped though pipes by some pressure gradient along the pipe called Poiseuille flows) are present everywhere in natural or engineered systems: blood flow through the body, from small capillaries to arteries, fluid flow through underground river systems, magma flows and pyroclastic flows through volcano chimneys, water flowing through a hose, a kitchen faucet, oil in a car engine, in a pipeline, etc.. These are often subject to strong shear if the wall boundaries are no-slip so fluid is moving in the center of the pipe, but not on the sides. Shear flows can also be driven by differential pressure gradients or any other forces) in open systems, and are found in the ocean, in the atmospheric wind patterns, in the surface and subsurface flows of the Sun, giant planets, other stars, in the orbital motion of gas in accretion disks, etc.. The reason why shear drives instabilities is most easily understood by considering the energetics of a sheared fluid, as shown in the next Section. The technique introduced there is actually widely applicable to any fluid instabilities, but it is worth bearing in mind that it is not very rigorous nor is it meant to be), by contrast with all the other techniques introduced in the previous chapter. 7. Energetics of shear instabilities Energetic considerations for the development of instabilities are based on the following idea: given a background fluid flow, is that background the lowest possible energy state of the system, or can energy be extracted from the background by mixing material around? If it can, then instabilities are indeed energetically 39
2 40 CHAPTER 7. SHEAR INSTABILITIES favorable because a perturbation can use the energy extracted from the background to amplify itself, in a positive feedback loop. These ideas are illustrated below, first on unstratified shear flows, then on convection, then on stratified shear flows. 7.. Energetics of unstratified inviscid shear flows Consider a simple background fluid flow, with constant density ρ, subject to some shear. Without loss of generality, the shear can be modeled locally as ū = ūz)e x = Sze x 7.) where S is the constant) shearing rate, by moving into a suitable frame of reference. This expression can, for instance, be viewed as a Taylor-expansion of the actual shear flow near a certain point. Consider two parcels of fluid, one at z = 0, and one at z = ɛ, where ɛ is small enough for 7.) to be a good representation of the local flow. The initial total kinetic energy contained in the two parcels is E i = ρ ū0) + ρ ūɛ) = ρ S ɛ 7.) Meanwhile, the momentum of the lower parcel is zero, while that of the upper parcel is ρūɛ) = ρsɛ. Next, suppose we switch the two parcels around, but in the process, mix and equalize their momenta. If the total momentum remains conserved in the process, then the momenta of each of the two switched parcels becomes ρsɛ/, and they are both moving at the average velocity Sɛ/. The total kinetic energy of the switched parcels is then E s = 4 ρs ɛ 7.3) We see that the new background flow consisting of the switched parcels has a lower energy than the original one. The difference E = E i E s = ρs ɛ /4 > 0 is available to amplify initial perturbations and thus drive instabilities. This is the main reason why shear instabilities are so common the shear itself is a reservoir of kinetic energy that can, under some circumstances, be tapped into to drive instabilities. 7.. Energetics of stratified fluids A very similar argument can be applied to determine whether switching two elements in a stratified fluid at rest is energetically favorable or unfavorable. Indeed, consider a fluid whose density is stratified, with ρz) = ρ m + ρ 0z z 7.4)
3 7.. ENERGETICS OF SHEAR INSTABILITIES 4 where ρ 0z is the constant density gradient. As before, we consider two parcels, one at z = 0, and one at z = ɛ. The initial total potential energy of the two parcels is E i = 0 + ρɛ)gɛ = gɛρ m + ρ 0z ɛ) 7.5) We then switch the two parcels after equalizing their density assuming mass is conserved). The density of each parcel, after the switch, is ρ m + ρ 0z ɛ/, so the potential energy is E s = gɛ ρ m + ρ ) 0zɛ 7.6) The difference in potential energy before and after the switch is E = E i E s = gρ 0zɛ 7.7) In this case, we see that if ρ 0z < 0, then the switch is energetically defavorable: E s is higher than E i, which means that we have to pay energy to make it happen. Meanwhile, if ρ 0z > 0, then the switch is energy-favorable as in the previous section, and the remaining E can be used to drive fluid instabilities. This result is not surprising at all: in the absence of dissipation, ρ 0z < 0 corresponds to stratified fluids which are unstable to convection, as studied in the last Chapter Energetics of stratified shear flows Finally, we can combine these two ideas to study the energetics of stratified shear flows. Starting with the background shear 7.) and background density profile 7.4), we consider again two parcels, and exchange them after equalizing their momenta and densities respectively. The total energy before the switch is now the sum of the kinetic and potential energies of both parcels: E i = ρ m + ρ 0z ɛ)s ɛ + gɛρ m + ρ 0z ɛ) 7.8) After the switch, both parcels have density ρ m + ρ 0z ɛ/, and their momentum is the average of the momenta of the two parcels before the switch which must now be calculated more carefully since the density is no longer constant): ρ m + ρ 0zɛ ) v = ρ m + ρ 0z ɛ)sɛ 7.9) which then yields v, the velocity of both parcels after the switch: v = ρ m + ρ 0z ɛ)sɛ ρ m + ρ0zɛ The total energy after the switch is then E s = ρ m + ρ ) 0zɛ v + gɛ ρ m + ρ ) 0zɛ 4 = ρ m + ρ 0z ɛ) S ɛ ρ m + ρ0zɛ + gɛ ρ m + ρ ) 0zɛ 7.0) 7.)
4 4 CHAPTER 7. SHEAR INSTABILITIES The energy gained or lost is E = E i E s = g ρ 0zɛ + ρ m + ρ 0z ɛ)s ɛ 4 ρ m + ρ 0z ɛ) S ɛ ρ m + ρ0zɛ 7.) Reducing this to the same denominator, and keeping only the lowest-order terms in ɛ, we get successively E = g ρ0zɛ ρm + ρ0zɛ ) + ρ m + ρ 0z ɛ)s ɛ ρ m + ρ0zɛ ) 4 ρ m + ρ 0z ɛ) S ɛ ɛ ) ρ ms + gρ 0z ρ m + ρ0zɛ In comparison with the previous sections, this expression can be interpreted in two equivalent ways. We saw that in the absence of stratification, shear instabilities are always energetically favorable. We now see that they remain energetically favorable in the presence of an unstable density gradient ρ 0z > 0 which is really not surprising), but can be stabilized by a sufficiently strong negative background density gradient. Alternatively, 7.3) can also be interpreted to mean that instabilities on a stably stratified flow can be energetically favorable provided the shear S is sufficiently strong The Richardson number for stratified shear flows The results of the previous section suggest that whether stratified shear flows are stable or unstable depends on the non-dimensional ratio J = gρ 0z ρ m S = N S 7.4) recalling that N can also be written as N = gρ 0z /ρ m. This is called the Richardson number after Richardson). Because J is based on local gradients, it is often called the gradient Richardson number. Global Richardson numbers can be defined as an alternative, see later on. The Richardson number, as we saw, is related to the ratio of the potential energy lost while moving density around, to the kinetic energy gained from the shear. For instability to occur, J should be smaller than a factor of order unity, while stability is likely when J is larger than unity. It is important to realize, however, that these are just approximate arguments based on energetic considerations. They do not replace a thorough linear stability analysis, nor do they provide any rigorous results on energy stability. However, they give us a good idea of what to expect next! 7. Local linear stability analysis As for the case of convection, we now attempt to perform a local stability analysis of shear flows. We do not need to specify the origin of the background 7.3)
5 7.3. GLOBAL ANALYSIS OF UNSTRATIFIED PLANE PARALLEL FLOWS43 shear ū see 7.), or of the background density profile ρ see 7.4), but merely assume they exist. We also assume for simplicity that we are in the limit where the Boussinesq approximation is valid. If we let then the perturbations evolve according to ρ = ρ + ρ and u = ū + ũ 7.5) ũ = 0 + ū ũ + ũ ū = p ρge z + ν ũ t ρ m ρ = α ρ T m T t + ū T + ũ T = κ T T 7.6) where T = ρ/αρ m = T 0z e z which defines T 0z ). Once expanded into components and assuming D flow, for instance), this becomes x + w z = 0 t w w + Sz t x T t p + Sz + S w = x x + ν ũ = p z + αg T + ν w + Sz T x + wt 0z = κ T T 7.7) We immediately detect a crucial problem: this equation does not have constant coefficients, so the standard local analysis which would have us use the ansatz q = ˆq expik x x + ik z z + λt) will not work. For this reason, we press on to a global mode analysis. 7.3 Global analysis of unstratified plane parallel flows 7.3. Linear stability analysis in the inviscid limit We now consider a global model of fluid flow between two horizontal parallel plates. For simplicity, we assume in this section that there is no viscosity. As we saw in the very first lecture, this implies that any fluid flow of the form ū = ūz)e x is a solution of the momentum equation subject to the boundary
6 44 CHAPTER 7. SHEAR INSTABILITIES conditions w = 0 at both plates. Perturbations to this background flow satisfy x + w z = 0 t + ūz) x w w + ūz) t x + w dū dz = p x = p z 7.8) The coefficients of this system of PDEs are independent of time and of x, so we seek solutions of the form q = ˆqz) expik x x + λt). We get Eliminating ˆp, we get λ + ik x ūz)) dû dz + ik x û + dŵ dz = 0 λû + ik x ūz)û + ŵ dū dz = ik x ˆp λŵ + ik x ūz)ŵ = dˆp dz 7.9) ik x û + dŵ ) dū dz dz + ŵ d ū dz = ik x λ + ik x ūz)) ŵ 7.0) We then use the continuity equation to eliminate û, to get ) dŵ λ + ik x ūz)) dz k xŵ ik x ŵ d ū dz = 0 7.) Finally, if we define a new variable c such that λ = ik x c 7.) then the vertical velocity perturbation must satisfy ) dŵ ūz) c) dz k xŵ ŵ d ū dz = 0 7.3) This equation is called the Rayleigh stability equation. It is an ODE which needs to be solved for the function ŵz) and the eigenvalue c, subject to the boundary conditions ŵ = 0 at the top and bottom of the domain. It usually needs to be solved numerically, and is not a standard Sturm-Liouville eigenvalue problem, so that most of the theorems we have learned about the latter do not apply here. Nevertheless, there are still an interesting number of strict results that can be derived from 7.3). First, note that for any solution with k x > 0, there is another solution with k x < 0 and c c. Without loss of generality, we can therefore take k x > 0. Next, note that unstable modes for k x > 0 are characterized by c values that have a positive imaginary part. Also it s easy to show that since ŵ can be
7 7.3. GLOBAL ANALYSIS OF UNSTRATIFIED PLANE PARALLEL FLOWS45 complex, if ŵ is a solution then its complex conjugate ŵ is also a solution with eigenvalue c. To see this, simply take the complex conjugate of the entire equation 7.3). This means that there are only two possibilities: either the solution ŵ has c real the mode is called neutrally stable since it neither grows nor decays, but merely oscillates), or there is a complex-conjugate pair of solutions ŵ and ŵ with one of the two being an unstable mode. If c is real, then there is a possibility that the quantity ūz) c) could be zero somewhere in the domain. If this is the case, then the ODE is singular at this point. We shall return to that case later. For growing modes, however, c has a non-zero imaginary part, and the problem is regular everywhere. Based on this consideration, we can then derive one of the most important results on plane parallel shear flows: Rayleigh s inflection point theorem. This theorem states that a necessary condition for the existence of an unstable mode is that the flow profile must have an inflection point in the domain, that is, some point where ū z) = d ū/dz = 0. To show this, first assume that c has a non-zero imaginary part, and rewrite 7.3) as d ŵ dz k xŵ ŵ ūz) c = 0 7.4) ū z) Then multiply this equation by ŵ, and integrate the result over z, from the bottom boundary at z = z ) to the top boundary at z = z ). This yields z z [ ŵ d ŵ dz k x ŵ ŵ ū z) ūz) c ] dz = 0 7.5) The first term in the integral can then be integrated by parts, to give [ z ] dŵ dz + kx ŵ + ŵ ū z) dz = 0 7.6) ūz) c z The imaginary part of this equation is then Ic) z ŵ ū z) dz = 0 7.7) z ūz) c Since all of the terms within the integral are strictly positive except ū z), there are two possibilities: either ū z) changes sign somewhere in the domain, or Ic) = 0 which we had assumed is not the case). This shows that Ic) 0 requires the presence of an inflection point! Note that the presence of an inflection point, by this analysis, is only a necessary condition for instability, but not a sufficient one, which means that there could be flows with an inflection point that are still stable. However, it does imply that flows without infection point are always stable in this inviscid limit). Rayleigh s inflection point theorem has a really important, and rather peculiar consequence: it implies that linear shear flows that is, flows of the kind
8 46 CHAPTER 7. SHEAR INSTABILITIES described by 7.)) are linearly stable to shear instabilities! This result is quite remarkable in its stark contrast with our energetic argument of Section 7., which suggested that linear shear flows in the inviscid limit should always be unstable. As it turns out, shear flows are quite different from the previously studied case of convection in that they are subject to finite amplitude instabilities, that is, they are often linearly stable, but can be destabilized by a strong enough perturbation. More on this later. In contrast, if the flow has an inflection point, then this is usually the position around which the shear instability will develop, and where the amplitude of the perturbation is the largest. A number of additional interesting theorems can be derived from this Rayleigh s instability equation. These include: Fjortoft s theorem: This is a stronger form of Rayleigh s inflection point criterion, which states that a necessary condition for instability is that ū z)ū ū i ) < 0 somewhere in the fluid, where ū i = ūz i ) and where z i is the point at which ū z i ) = 0. Howard s semicircle theorem: This is a very useful theorem that bounds the possible values of the real and imaginary parts of c, and therefore places upper limits on the growth rate of the unstable modes. The theorem states that all unstable modes that is with Ic) > 0) have an eigenvalue c that satisfies Rc) ) max ū + min ū + Ic) max ū min ū) < 7.8) or in other words, the complex number c lies on the complex plane within the semi-circle of radius R = max ū min ū)/4, centered on the point max ū + min ū)/ on the real axis. The proof of these theorems can be found in the textbook Hydrodynamic Stability by Drazin and Reid, for instance, or in the original papers. Finally, the combination of these theorems leads to the following strong statements whose proofs are beyond the scope of this class). If ūz) is a smooth flow, then eigenmodes of Rayleigh s instability equation are of two kinds: Neutrally stable modes, with Ic) = 0. There is a continuum of them, one for each possible real) value of c in the interval [min ū, max ū]. The modes are singular at the point z where ūz) = c, and their first derivative is discontinuous there. The position of this singularity is called a critical layer. Pairs of complex-conjugate modes, with Ic) 0. There are only very few of them, at most one pair for each inflection point in the flow, and c for these complex-conjugate pairs satisfies Howard s semicircle theorem.
Prototype Instabilities
Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies
More informationChapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow
Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already
More information4 Shear flow instabilities
4 Shear flow instabilities Here we consider the linear stability of a uni-directional base flow in a channel u B (y) u B = 0 in the region y [y,y 2. (35) 0 y = y 2 y z x u (y) B y = y In Sec. 4. we derive
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationPAPER 331 HYDRODYNAMIC STABILITY
MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal
More informationFor example, for values of A x = 0 m /s, f 0 s, and L = 0 km, then E h = 0. and the motion may be influenced by horizontal friction if Corioli
Lecture. Equations of Motion Scaling, Non-dimensional Numbers, Stability and Mixing We have learned how to express the forces per unit mass that cause acceleration in the ocean, except for the tidal forces
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationKelvin Helmholtz Instability
Kelvin Helmholtz Instability A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram November 00 One of the most well known instabilities in fluid
More information196 7 atmospheric oscillations:
196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationHydrodynamic Instability. Agenda. Problem definition. Heuristic approach: Intuitive solution. Mathematical approach: Linear stability theory
Environmental Fluid Mechanics II Stratified Flow and Buoyant Mixing Hydrodynamic Instability Dong-Guan Seol INSTITUTE FOR HYDROMECHANICS National Research Center of the Helmholtz Association www.kit.edu
More informationCoherent structures in stably stratified plane Couette flow
Coherent structures in stably stratified plane Couette flow D. Olvera * & R. R. Kerswell School of Mathematics, University of Bristol, Bristol, UK. * do2542@bristol.ac.uk Abstract A large body of recent
More information2.5 Sound waves in an inhomogeneous, timedependent
.5. SOUND WAVES IN AN INHOMOGENEOUS, TIME-DEPENDENT MEDIUM49.5 Sound waves in an inhomogeneous, timedependent medium So far, we have only dealt with cases where c was constant. This, however, is usually
More informationτ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization
τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization Note that w *3 /z i is used to normalized the TKE equation in case of free
More informationMAE 210C FLUID MECHANICS III SPRING 2015 Rayleigh s equation for axisymmetric flows Forthestabilityanalysisofparallelaxisymmetricflowswithbasicvelocity profileū = (U x,u r,u θ ) = (U(r),0,0) weintroducemodalperturbationsoftheform(v
More informationLinear Hydrodynamic Stability Analysis Summary and Review. Example of Stability Problems. Start with Base Velocity Profile
Linear Hydrodynamic Stability Analysis Summary and Review OCEN 678 Fall 2007 Scott A. Socolofsky Example of Stability Problems Shear Flows: flows with strong velocity gradients Classic examples Wakes Jets
More information5.3 Inviscid instability mechanism of parallel flows
1 Lecture Notes on Flui Dynamics 1.63J/2.21J) by Chiang C. Mei, 27 5.3 Invisci instability mechanism of parallel flows We now turn to an oler problem of the instability of parallel flow without stratification
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More information5.3 Inviscid instability mechanism of parallel flows
1 Lecture Notes on Flui Dynamics 1.63J/2.21J) by Chiang C. Mei, 22 5.3 Invisci instability mechanism of parallel flows We now turn to an oler problem of the instability of parallel flow without stratification
More informationCHAPTER 19. Fluid Instabilities. In this Chapter we discuss the following instabilities:
CHAPTER 19 Fluid Instabilities In this Chapter we discuss the following instabilities: convective instability (Schwarzschild criterion) interface instabilities (Rayleight Taylor & Kelvin-Helmholtz) gravitational
More informationCHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY
1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 00 CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY References: Drazin: Introduction to Hydrodynamic Stability Chandrasekar: Hydrodynamic
More informationz g + F w (2.56) p(x, y, z, t) = p(z) + p (x, y, z, t) (2.120) ρ(x, y, z, t) = ρ(z) + ρ (x, y, z, t), (2.121)
= + dw dt = 1 ρ p z g + F w (.56) Let us describe the total pressure p and density ρ as the sum of a horizontally homogeneous base state pressure and density, and a deviation from this base state, that
More informationCoupled Reduced Equations for Strongly Stratified Flows
Coupled Reduced Equations for Strongly Stratified Flows Cesar B. Rocha October 15, 2015 Abstract We present a set of reduced equations in the limit of strong stratification. The asymptotics lead to the
More informationLecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev
Lecture Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Introduction In many cases in nature, like in the Earth s atmosphere, in the interior of stars and planets, one sees the appearance
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More information3 Hydrostatic Equilibrium
3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose
More informationGoals of this Chapter
Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence
More information9 Fluid Instabilities
9. Stability of a shear flow In many situations, gaseous flows can be subject to fluid instabilities in which small perturbations can rapidly flow, thereby tapping a source of free energy. An important
More informationThe perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh
z = The perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh, that is, p = p h + p nh. (.1) The former arises
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationP. Garaud & L. Acevedo-Arreguin
On the Penetration of Meridional Circulation below the Solar Convection Zone II: Models with Convection Zone, the Taylor-Proudman constraint and Applications to Other Stars. P. Garaud & L. Acevedo-Arreguin
More informationINTERNAL GRAVITY WAVES
INTERNAL GRAVITY WAVES B. R. Sutherland Departments of Physics and of Earth&Atmospheric Sciences University of Alberta Contents Preface List of Tables vii xi 1 Stratified Fluids and Waves 1 1.1 Introduction
More informationHorizontal buoyancy-driven flow along a differentially cooled underlying surface
Horizontal buoyancy-driven flow along a differentially cooled underlying surface By Alan Shapiro and Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, OK, USA 6th Baltic Heat Transfer
More informationThe Reynolds experiment
Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationAnalysis of Least Stable Mode of Buoyancy Assisted Mixed Convective Flow in Vertical Pipe Filled with Porous Medium
Proceedings of the World Congress on Engineering 2 Vol I WCE 2, July 6-8, 2, London, U.K. Analysis of Least Stable Mode of Buoyancy Assisted Mixed Convective Flow in Vertical Pipe Filled with Porous Medium
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationIn two-dimensional barotropic flow, there is an exact relationship between mass
19. Baroclinic Instability In two-dimensional barotropic flow, there is an exact relationship between mass streamfunction ψ and the conserved quantity, vorticity (η) given by η = 2 ψ.the evolution of the
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationDynamics Rotating Tank
Institute for Atmospheric and Climate Science - IACETH Atmospheric Physics Lab Work Dynamics Rotating Tank Large scale flows on different latitudes of the rotating Earth Abstract The large scale atmospheric
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More informationBoundary Layers. Lecture 2 by Basile Gallet. 2i(1+i) The solution to this equation with the boundary conditions A(0) = U and B(0) = 0 is
Boundary Layers Lecture 2 by Basile Gallet continued from lecture 1 This leads to a differential equation in Z Z (A ib) + (A ib)[ C + Uy Uy 2i(1+i) ] = 0 with C = 2. The solution to this equation with
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationOnset of convection of a reacting fluid layer in a porous medium with temperature-dependent heat source
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2011, Science Huβ, http://www.scihub.org/ajsir ISSN: 2153-649X doi:10.5251/ajsir.2011.2.6.860.864 Onset of convection of a reacting fluid layer in
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationNonmodal Growth and the Unstratified MRI Dynamo
MPPC general meeting, Berlin, June 24 Nonmodal Growth and the Unstratified MRI Dynamo Jonathan Squire!! and!! Amitava Bhattacharjee MRI turbulence Observed accretion rate in disks in astrophysical disks
More informationBoundary Layers: Stratified Fluids
Boundary Layers: Stratified Fluids Lecture 3 by Jeroen Hazewinkel continued from lecture 2 Using w I = /(2σS) 2 T I, the interior of the cylinder is be described by 2σS 2 T I + 1/2 4 2 This result can
More informationPhysics 212: Statistical mechanics II Lecture IV
Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most
More informationLecture 6.1 Work and Energy During previous lectures we have considered many examples, which can be solved using Newtonian approach, in particular,
Lecture 6. Work and Energy During previous lectures we have considered many examples, which can be solved using Newtonian approach, in particular, Newton's second law. However, this is not always the most
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationChapter 7. Linear Internal Waves. 7.1 Linearized Equations. 7.2 Linear Vertical Velocity Equation
Chapter 7 Linear Internal Waves Here we now derive linear internal waves in a Boussinesq fluid with a steady (not time varying) stratification represented by N 2 (z) andnorotationf =0. 7.1 Linearized Equations
More informationStratified Shear Flow: Instability and Wave Radiation
Stratified Shear Flow: Instability and Wave Radiation B. R. Sutherland Dept. Mathematical and Statistical Sciences, University of Alberta Edmonton, Canada Abstract The stability and evolution of a two-dimensional
More informationStability of Vertically and Radially Stratified Protoplanetary Disks
Stability of Vertically and Radially Stratified Protoplanetary Disks Figure from Stoll & Kley (2014) Glen Stewart Laboratory for Atmospheric and Space Physics University of Colorado Vertical Shear instability
More informationChapter 13 Instability on non-parallel flow Introduction and formulation
Chapter 13 Instability on non-parallel flow. 13.1 Introduction and formulation We have concentrated our discussion on the instabilities of parallel, zonal flows. There is the largest amount of literature
More informationES265 Order of Magnitude Phys & Chem Convection
ES265 Order of Magnitude Phys & Chem Convection Convection deals with moving fluids in which there are spatial variations in temperature or chemical concentration. In forced convection, these variations
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More information3.8 The First Law of Thermodynamics and the Energy Equation
CEE 3310 Control Volume Analysis, Sep 30, 2011 65 Review Conservation of angular momentum 1-D form ( r F )ext = [ˆ ] ( r v)d + ( r v) out ṁ out ( r v) in ṁ in t CV 3.8 The First Law of Thermodynamics and
More informationDonald Slinn, Murray D. Levine
2 Donald Slinn, Murray D. Levine 2 Department of Civil and Coastal Engineering, University of Florida, Gainesville, Florida College of Oceanic and Atmospheric Sciences, Oregon State University, Corvallis,
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationVII. Hydrodynamic theory of stellar winds
VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
More informationBounds on Surface Stress Driven Flows
Bounds on Surface Stress Driven Flows George Hagstrom Charlie Doering January 1, 9 1 Introduction In the past fifteen years the background method of Constantin, Doering, and Hopf has been used to derive
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More information1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics
1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics Theorem. Consider the equation u t = Du xx + f(u) with f(0) = f(1) = 0, f(u) > 0 on 0 < u < 1, f (0)
More informationSG2221 Wave Motion and Hydrodinamic Stability. MATLAB Project on 2D Poiseuille Flow Alessandro Ceci
SG2221 Wave Motion and Hydrodinamic Stability MATLAB Project on 2D Poiseuille Flow Alessandro Ceci Base Flow 2D steady Incompressible Flow Flow driven by a constant pressure gradient Fully developed flow
More informationNonlinear evolution of linear optimal perturbations of strongly stratified shear layers
Under consideration for publication in J. Fluid Mech. Nonlinear evolution of linear optimal perturbations of strongly stratified shear layers A. K. Kaminski,, C. P. Caulfield 3, and J. R. Taylor Department
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 5 NATURAL CONVECTION HEAT TRANSFER BASIC CONCEPTS MECHANISM OF NATURAL
More informationKinematics of fluid motion
Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 27 May, 2004 1.30 to 3.30 PAPER 64 ASTROPHYSICAL FLUID DYNAMICS Attempt THREE questions. There are four questions in total. The questions carry equal weight. Candidates
More informationDirect numerical simulation of salt sheets and turbulence in a double-diffusive shear layer
Direct numerical simulation of salt sheets and turbulence in a double-diffusive shear layer Satoshi Kimura, William Smyth July 3, 00 Abstract We describe three-dimensional direct numerical simulations
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More information1. The vertical structure of the atmosphere. Temperature profile.
Lecture 4. The structure of the atmosphere. Air in motion. Objectives: 1. The vertical structure of the atmosphere. Temperature profile. 2. Temperature in the lower atmosphere: dry adiabatic lapse rate.
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationA note concerning the Lighthill sandwich model of tropical cyclones
Applied Mathematics A note concerning the Lighthill sandwich model of tropical cyclones G.I. Barenblatt, A.J. Chorin, V.M. Prostokishin The authors dedicate this work to the glowing memory of the great
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More information2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven Ocean Fronts
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. 2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven
More informationSurface tension driven oscillatory instability in a rotating fluid layer
J. Fluid Mech. (1969), wol. 39, part 1, pp. 49-55 Printed in Great Britain 49 Surface tension driven oscillatory instability in a rotating fluid layer By G. A. McCONAGHYT AND B. A. FINLAYSON University
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More informationUNIT II Real fluids. FMM / KRG / MECH / NPRCET Page 78. Laminar and turbulent flow
UNIT II Real fluids The flow of real fluids exhibits viscous effect that is they tend to "stick" to solid surfaces and have stresses within their body. You might remember from earlier in the course Newtons
More informationTheory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 7 Instability in Rotor Systems Lecture - 2 Fluid-Film Bearings
More informationHurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.
Hurricanes: Observations and Dynamics Houze Section 10.1. Holton Section 9.7. Emanuel, K. A., 1988: Toward a general theory of hurricanes. American Scientist, 76, 371-379 (web link). http://ww2010.atmos.uiuc.edu/(gh)/guides/mtr/hurr/home.rxml
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationCalculus Math 21B, Winter 2009 Final Exam: Solutions
Calculus Math B, Winter 9 Final Exam: Solutions. (a) Express the area of the region enclosed between the x-axis and the curve y = x 4 x for x as a definite integral. (b) Find the area by evaluating the
More informationCan a Simple Two-Layer Model Capture the Structure of Easterly Waves?
Can a Simple Two-Layer Model Capture the Structure of Easterly Waves? Cheryl L. Lacotta 1 Introduction Most tropical storms in the Atlantic, and even many in the eastern Pacific, are due to disturbances
More information4.7 Dispersion in an oscillatory shear flow
Lecture notes in Fluid ynamics.63j/.0j) by Chiang C. Mei, MIT, Spring, 007 4-6dispersion.tex March, 007 [Refs]:. Aris:. Fung, Y. C. Biomechanics 4.7 ispersion in an oscillatory shear flow Relevant to the
More informationVortices in accretion discs: formation process and dynamical evolution
Vortices in accretion discs: formation process and dynamical evolution Geoffroy Lesur DAMTP (Cambridge UK) LAOG (Grenoble) John Papaloizou Sijme-Jan Paardekooper Giant vortex in Naruto straight (Japan)
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationOn the turbulent Prandtl number in homogeneous stably stratified turbulence
J. Fluid Mech. (2010), vol. 644, pp. 359 369. c Cambridge University Press 2010 doi:10.1017/s002211200999293x 359 On the turbulent Prandtl number in homogeneous stably stratified turbulence SUBHAS K. VENAYAGAMOORTHY
More informationNonlinear Vertical Scale Selection in Equatorial Inertial Instability
1APRIL 003 GRIFFITHS 977 Nonlinear Vertical Scale Selection in Equatorial Inertial Instability STEPHEN D. GRIFFITHS Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge,
More informationGlobal Magnetorotational Instability with Inflow
Global Magnetorotational Instability with Inflow Evy Kersalé PPARC Postdoctoral Research Associate Dept. of Appl. Maths University of Leeds Collaboration: D. Hughes & S. Tobias (Appl. Maths, Leeds) N.
More information