Prototype Instabilities
|
|
- Laureen Lawrence
- 5 years ago
- Views:
Transcription
1 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 of the mean flow. The first type of instability is sometimes called mechanical, and is associated with the shear of the mean flow. The second is called convective, and is associated with the vertical distribution of temperature or density of the mean state. Both types of instability are important for the boundary layer and for small-scale dynamics in general. The two varieties of mechanical instability that we will consider are the inflectional or Rayleigh instability, which is always associated with maximum in the vorticity distribution, and inertial instability, which can only occur in a rotating fluid, and which results from an unstable distribution of angular momentum. Inflectional instability is certainly responsible for the maintenance of turbulence in the stable PBL, and may also play an important role in some other boundary layer phenomena, such as cloud streets and billows. Inertial instability can also occur within both the boundary layer and the cumulus layer. Ekman layers occur in viscous, rotating fluids, and are observed to be subject to Ekman layer instability. The relationships of Ekman layer instability to inflectional instability and inertial instability have been the subject of some discussion in the literature. Convective instabilities play a prominent role in this course, because they maintain the convective turbulence of the deep, unstable PBL, and also lead to cumulus convection. The growth rates and characteristic shapes of convective disturbances are strongly influenced by shear in the atmosphere, and may occasionally be significantly modified by rotation. Both types of instabilities tend to relieve the conditions which lead to their growth, by stabilizing the basic state. Shearing instability transports momentum down gradient and so tends to reduce the shear. Convective instability transports heat upward and so tends to reduce the lapse rate. The two types of instability often interact with each other. For example, a convective thermal creates shear zones with inflectional instability near the its edge. As a second example, finite-amplitude shearing instability can create local patches of convective instability in the form of breaking waves. In this chapter, we use linear theory to investigate the essential mechanisms of each instability. 1
2 Shearing Instability The simplest model of inflectional instability is based on consideration of a vortex sheet, as pictured in Figure 1. Part (a) of the figure depicts a stable configuration, in which the vorticity increases monotonically from the lower to the upper part of the figure. A vortex displaced, as indicated, from its equilibrium position will tend to be carried back by the circulations induced by the neighboring vortices. But if the vorticity has a local maximum, as in part (b) of the figure, then a displaced vortex tends to be carried further away. Such a vorticity distribution is therefore unstable. Infinite vortex sheet (balanced) Perturb one vortex Instability Figure 1: Diagram illustrating the instability of a vortex sheet. The argument above is based on a qualitative appeal to intuition. The following more quantitative analysis is based on a portion of Chapter XI of a book by S. Chandrasekhar (1961). 2
3 Consider stratified motion in the x - direction, and sheared in the z -direction, so that the mean flow is given by U z ( ). We adopt the linearized Boussinesq system of equations, which have the following form: u u +U x + w du dz = p x, v (1) v +U x = p y, (2) w +U w x = p z gρ, (3) ρ ρ +U x + w z = 0, (4) u x + v y + w z = 0. (5) In (4), the stratification enters in the form of the vertical gradient of the mean density. We assume that decreases upward, so that the stratification gives a gravitational restoring force that permits gravity waves. We seek separable solutions of the form ( ) = ( ~ )exp i( kx + ly + σt). With this assumed form, the amplitudes depend only on z and instability can occur only if σ is complex. Substituting (6) into (1-5), we obtain i( σ + ku )u + w du dz = ikp, (7) (6) i( σ + ku )v = ilp, (8) 3
4 i( σ + ku )w = p z gρ, (9) i( σ + ku )ρ + w d dz = 0, (10) i( ku + lv) + w z = 0. (11) Using (7-9) to form the divergence equation, and then substituting from (5), we obtain i σ + ku ( ) w ( ) p. z k du dz w = k 2 + l 2 This is the anelastic pressure equation, which can be used to determine the pressure diagnostically given appropriate boundary conditions, although we will not actually solve (12) here. Elimination of the perturbation density between Eqs. (9) and (10) gives (12) i( σ + ku) 2 w = ( σ + ku) p z igw d dz, (13) which describes gravity waves. Elimination of the pressure between (12) and (13) gives an equation in which the only unknown is the vertical velocity: z ρ 0 σ + ku ( ) w z kw du dz ( k 2 + l 2 ) ( σ + ku )w = g k 2 + l 2 ( ) d dz w σ + ku. Before we can solve this equation, we need some far-field boundary conditions on the vertical velocity. The simplest case is an unbounded domain in which the motions die out far from the shear zone: w 0 as z ±. We also want to allow the possibility of discontinuous changes of U and at an internal boundary defined by z = z s ( x, y,t), which is assumed to be a material surface. To derive the form of the kinematic boundary condition at z = z s ( x, y,t), we revert to the full continuity equation, which is (14) (15) 4
5 ρ + i ρv ( ) = 0. Integrating (16) across a thin layer defined by z s ε z z s + ε, we obtain (16) z s +ε ρ + i( ρv) dz = 0. z s ε Here ε is a small constant that we will shrink to zero just below. Taking the derivatives inside the integrals, and using Leibniz rule, we find that z s +ε z s ε ρ s+ z s + ε +ρ s z s ε ρ dz + i z s +ε z s ε ρvdz ( ) + V i ( z s + ε ) w s+ ( ) + V i ( z s ε ) w s (17) = 0. (18) As ε 0, the integrals on the first line of (18) go to zero, and their derivatives also go to zero, so the terms drop out. We re left with where Δ( ) ( ) s+ ( ) s. We can write Δ ρ z s + V z w s = 0, (19) and z ρ s+ s+ + V s+ z s+ w s+ = F ρ s z s + V s z s w s = F, (20) (21) which automatically satisfy (19). Here F is the mass flux across z s ( x, y,t). If z s is a material surface, then F = 0, so that 5
6 and z s+ + V z s w s+ = 0, (22) z s + V z s w s = 0. If we require that z s+ and z s stay together, so that z s+ = z s = z s, then follows from (22) and (23) that Δ( V z s w) = 0. Eq. (19) says that the mass flux must be continuous, while (24) says that the normal component of the velocity must be continuous. We now return to our earlier analysis. Substitution of (6) into (22) and (23) gives (23) (24) and i( σ + ku + )z s w s+ = 0, (25) respectively. These can be combined to obtain i( σ + ku )z s w s = 0, w Δ σ + ku = 0. Similarly, by vertical integration of (14) we find that Δ σ + ku ( ) w z kw du dz = g ( k 2 + l 2 w )Δ σ + ku s. (26) (27) (28) Here we have used (27) to justify writing w σ + ku two interfacial conditions have to be satisfied, namely (27) and (28). s on the right-hand side of (28). Note that 6
7 We now consider two cases. First, suppose that both U and are discontinuous at the material surface z s, but also that U and are uniform throughout the rest of the domain. Then (14) simplifies to 2 z k ( 2 + l 2 ) 2 w = 0, which applies everywhere except at the interface. Let subscripts 1 and 2 denote the regions below and above the interface, respectively. In view of the far-field boundary conditions given by (17) and the first interfacial condition (27), the solution for w can be written as w 1 = A( σ + ku 1 )exp k 2 + l 2 z w 2 = A( σ + ku 2 )exp k 2 + l 2 z (29) ( ) for z < 0, and ( ) for z > 0, (30) where A is an arbitrary constant. Substitution of (30) into the second interfacial condition, (28), gives, after a few lines of algebra, where we define σ 2 + 2k( α 1 U 1 + α 2 U 2 )σ + k 2 α 1 U 2 2 ( 1 + α 2 U 2 ) = g k 2 + l 2 ( α 1 α 2 ), (31) α i = ρ i ρ 1 + ρ 2, i = 1,2. We can solve (31) as a quadratic equation for σ. After some algebra, the result can be written as (32) ( ) ± g k 2 + l 2 ( α 1 α 2 ) k 2 α 1 α 2 ( U 1 U 2 ) 2 σ = k α 1 U 1 + α 2 U 2 To obtain (33), we have used α 1 + α 2 = 1, which follows from (32). 1 2 (33) Recall that, in order to have instability, σ must have an imaginary part. Eq. (33) shows that instability will definitely occur if the density increases upward, i.e., if α 1 α 2 < 0. In the following discussion we restrict ourselves to the statically stable case α 1 α 2 > 0. For k = 0, U 1, and U 2 drop out of (33), and we just get a gravity wave. This implies that when the interface 7
8 is statically stable all growing disturbances have k > 0. The condition for instability with α 1 α 2 > 0 can be written as g α k 2 + l 2 > 1 α 2 α 1 α 2 U 1 U 2 ( ) ( ) 2 cos 2 θ, (34) where θ is the angle between the wave number vector and the x -direction. The minimum wave number occurs for cos 2 θ = 1, which means that variations occur only in the x -direction. No matter how small U 1 U 2 is made, instability will occur for sufficiently large k 2 + l 2. If the density difference is set to zero, then all modes are unstable. We conclude that the stable stratification stabilizes the longer waves. Now consider the case of continuous U and, but still with a shear layer centered at z = 0. First, consider the following simple analysis. Consider two parcels at heights z and z + dz. We mix the two parcels together, conserving the mean density and mean wind. The increase in the potential energy of the system, due to the work done against the buoyancy force, is g dρ ( )dz, per unit volume. The decrease in the kinetic energy of the mean state is ρ 1 2 U ( 2 U + du ) U +U + du 2 ( )2 = 1 4 ρ( du )2. Recall that the mechanism of shearing instability involves conversion of the kinetic energy of the mean state into the kinetic energy of the growing disturbance. In order to have instability in the present example, the decrease in mean-state kinetic energy must be larger than the work done against buoyancy. This condition can be written as (35) R i < 1 4, (36) where R i g ρ z ρ ( U z) 2 (37) is the Richardson number, or more precisely the gradient Richardson number. Note that the sign of the Richardson number is determined by the sign of d dz. Since we are assuming that the density decreases upward, the Richardson number will be positive. The simple argument 8
9 given above does not show that (36) is a sufficient condition for instability - only that it is a necessary condition. We now consider the problem more carefully. The following analysis is due to Drazin (1958), but is summarized by Chandrasekhar (1961). Let the mean-state density decrease exponentially upward, so that and assume that the wind profile takes the form ρ = 0 e βz, U = U 00 tanh( z d), where β and d are constants, so that the shear is concentrated near z = 0. Then (38) (39) ( R i ) 0 = gβd 2 U 2 00 is a representative value of R i in the vicinity of z = 0. Taking l = 0 for simplicity, we expand (14) as ( σ + ku ) 2 z k 2 2 w k d 2 U dz w gk 2 d w 2 ρ o dz σ + ku + 1 d σ + ku dz Two of the terms in (41) involve d ( ) w (40) z k du dz w = 0. (41) dz. In geophysical applications, the second-last term, which is proportional to g, is much larger than the last term. The physical interpretation is that d dz matters for its effects on the potential energy, but not so much for its effects on the momentum. Neglecting the last term, (41) can be simplified to σ k +U 2 z k 2 w d 2 U 2 dz w g d w 2 dz σ +U 0. k We now nondimensionalize by scaling velocities with U 00 and lengths with d, and define a nondimensional phase speed (42) c σ ku 00. (43) 9
10 Using the definition (43), we can rewrite (42) as ( U c) 2 z k 2 2 w d 2 U dz w + R 2 i w ( ) 0 U c = 0, where all variables are now nondimensional and the symbol U now stands for tanh z. We use the boundary conditions w 0 as z ±. Recall from (6) that instability occurs when the imaginary part of σ is negative. Since c is a negative real number times σ, we can also say that instability occurs when the imaginary part of c is positive. Instability may or may not occur, depending on the values of the parameters. Imagine that we start with parameter values for which instability does not occur, and then progressively change the parameters until we reach a point beyond which instability begins. In the stable regime, disturbances are damped through the action of a restoring force. In the unstable regime, the restoring force no longer restores; instead, it makes the disturbance grow. In the marginal state that separates stability from instability, the disturbance is neutrally stable, and the restoring force is equal to zero. We want to find the parameter values that are characteristic of the marginal state. It is going to turn out that the marginal state occurs for ( R i ) 0 =1/ 4. In the marginal state, the imaginary parts of σ and c are equal to zero. What about their real parts? If the real parts of of σ and c were also equal to zero in the marginal state, then we would simply have σ = c = 0, which would allow us to simplify (44). If the real part of σ is different from zero in the marginal state, then the unstable disturbance appears as a growing oscillation, which means that it grows through the action of an overshooting restoring force. This is called overstability. In contrast, if the real part of σ is zero in the marginal state, then the unstable disturbance appears as a stationary (i.e., nonoscillatory) mode. When the growing disturbance is non-oscillatory, the principle of the exchange of stabilities is said to apply, and the marginal state has σ = c = 0. The terms introduced above are due to Eddington (1926, p. 201). We now assume that the principle of the exchange of stabilities applies to the problem at hand. We therefore set c = 0 in (44), giving U 2 z k 2 2 w d 2 U dz 2 (44) (45) w + R w ( i ) 0 U = 0. (46) 10
11 Because (46) has been obtained from (44) by setting c = 0, solutions of (46) will describe the marginal state. To find the solution of (46), we first make a change of variables. Since U = tanh z is a single-valued function of z, we can use U instead of z as our vertical coordinate. After making this change of variables, (46) becomes d du which is subject to the boundary conditions ( ) 0 ( ) dw 1 U2 du + 2 k 2 1 U 2 + R i U 2 1 U 2 w = 0 at U = ±1. ( ) w = 0, You can confirm by substitution that a solution of (47) that satisfies (48) is (47) (48) where χ is a constant, and we define w = χu µ ( 1 U 2 ) υ, (49) ( ) 0 µ R i, (50) and υ 1 2 k 2 ( R i ) 0. (51) Upon substituting (49) into (47), you will find that ( 2υ + µ + 2) ( 2υ + µ 1) = 0 (52) is the dispersion relation that relates k to ( R i ) 0. Substitution for υ and µ shows that the second factor on the left-hand side of (52) will be equal to zero if we set ( R i ) 0 = k 2 ( 1 k 2 ). (53) 11
12 From (53), we see that ( R i ) 0 has a maximum of 1/ 4 for k 2 =1/ 2. If ( R i ) 0 is larger than 1/ 4, then (53) does not have a solution, which means that instability cannot occur. Suppose that we run a series of experiments, starting with ( R i ) 0 >1/ 4, and then testing progressively smaller values of ( R i ) 0. Instability will not occur until we reach ( R i ) 0 =1/ 4. We say that ( R i ) 0 =1/ 4 is characteristic of the marginal state of the system, and that ( R i ) 0 =1/ 4 is the critical value of the Richardson number. The first sign of instability, at ( R i ) 0 =1/ 4, will be a ( stationary, i.e., non-oscillatory) growing disturbance with k 2 = 1 2. Howard (1961; see also Miles, 1959) proved that if R i exceeds a critical value (essentially 1/4) throughout the fluid, then the motion must be stable in the inflectional sense, and that if the motion is unstable, then R i must fall below the critical value somewhere in the fluid. Notice that this does not prove that instability must occur for R i less than critical. Dutton (1976) discusses Howard s theorem in some detail. An analogous theorem exists for the instability of quasigeostrophic flow. References and Bibliography Chandrasekhar, S., 1961: Hydrodynamic and hydromagnetic stability. Clarendon Press, Oxford, 652 pp. Drazin, P. G., 1958: The stability of a shear layer in an unbounded heterogeneous inviscid fluid. J. Fluid Mech., 4, Dutton, J. A., 1976: The ceaseless wind. McGraw-Hill, New York, 579 pp. Eddington, A. S., 1926: The internal constitution of the stars. Cambridge University Press, 407 pp. Miles, J. W., 1959: On the generation of surface waves by shear flows. III. Kelvin-Helmholtz instability. J. Fluid Mech., 3,
General 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 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 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 informationShear instabilities. Chapter Energetics of shear instabilities
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
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 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 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 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 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 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 informationTheory of linear gravity waves April 1987
April 1987 By Tim Palmer European Centre for Medium-Range Weather Forecasts Table of contents 1. Simple properties of internal gravity waves 2. Gravity-wave drag REFERENCES 1. SIMPLE PROPERTIES OF INTERNAL
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 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 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 informationBulk Boundary-Layer Model
Bulk Boundary-Layer Model David Randall Ball (1960) was the first to propose a model in which the interior of the planetary boundary layer (PBL) is well-mixed in the conservative variables, while the PBL
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 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 informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
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 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 informationConvection and buoyancy oscillation
Convection and buoyancy oscillation Recap: We analyzed the static stability of a vertical profile by the "parcel method"; For a given environmental profile (of T 0, p 0, θ 0, etc.), if the density of an
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 informationReferences: Parcel Theory. Vertical Force Balance. ESCI Cloud Physics and Precipitation Processes Lesson 3 - Stability and Buoyancy Dr.
References: ESCI 340 - Cloud Physics and Precipitation Processes Lesson 3 - Stability and Buoyancy Dr. DeCaria Glossary of Meteorology, 2nd ed., American Meteorological Society A Short Course in Cloud
More informationFigure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ
6. Internal waves Consider a continuously stratified fluid with ρ o (z) the vertical density profile. z p' ξ p ρ ρ ο (z) Figure 1. Figure by MIT OpenCourseWare. At a point P raise a parcel of water by
More informationMagnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium
Mechanics and Mechanical Engineering Vol. 12, No. 3 (2008) 267 274 c Technical University of Lodz Magnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium Emad E. Elmahdy Mathematics
More informationThe Atmospheric Boundary Layer. The Surface Energy Balance (9.2)
The Atmospheric Boundary Layer Turbulence (9.1) The Surface Energy Balance (9.2) Vertical Structure (9.3) Evolution (9.4) Special Effects (9.5) The Boundary Layer in Context (9.6) Fair Weather over Land
More informationThe atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:
The atmospheric boundary layer: Utrecht Summer School on Physics of the Climate System Carleen Tijm-Reijmer IMAU The atmospheric boundary layer: Where the atmosphere meets the surface Photo: Mark Wolvenne:
More informationNewton's Laws You should be able to state these laws using both words and equations.
Review before first test Physical Mechanics Fall 000 Newton's Laws You should be able to state these laws using both words and equations. The nd law most important for meteorology. Second law: net force
More informationSAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany
WAVES IN THE OCEANS Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany Keywords: Wind waves, dispersion, internal waves, inertial oscillations, inertial waves,
More informationInstabilities due a vortex at a density interface: gravitational and centrifugal effects
Instabilities due a vortex at a density interface: gravitational and centrifugal effects Harish N Dixit and Rama Govindarajan Abstract A vortex placed at an initially straight density interface winds it
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationAtmospheric Dynamics: lecture 2
Atmospheric Dynamics: lecture 2 Topics Some aspects of advection and the Coriolis-effect (1.7) Composition of the atmosphere (figure 1.6) Equation of state (1.8&1.9) Water vapour in the atmosphere (1.10)
More informationGravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017
Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature
More informationFloquet Theory for Internal Gravity Waves in a Density-Stratified Fluid. Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 2012
Floquet Theory for Internal Gravity Waves in a Density-Stratified Fluid Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 212 Density-Stratified Fluid Dynamics Density-Stratified
More informationρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
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 informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationGeneralizing the Boussinesq Approximation to Stratified Compressible Flow
Generalizing the Boussinesq Approximation to Stratified Compressible Flow Dale R. Durran a Akio Arakawa b a University of Washington, Seattle, USA b University of California, Los Angeles, USA Abstract
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More informationModule 3: "Thin Film Hydrodynamics" Lecture 12: "" The Lecture Contains: Linear Stability Analysis. Some well known instabilities. Objectives_template
The Lecture Contains: Linear Stability Analysis Some well known instabilities file:///e /courses/colloid_interface_science/lecture12/12_1.htm[6/16/2012 1:39:16 PM] Linear Stability Analysis This analysis
More information( u,v). For simplicity, the density is considered to be a constant, denoted by ρ 0
! Revised Friday, April 19, 2013! 1 Inertial Stability and Instability David Randall Introduction Inertial stability and instability are relevant to the atmosphere and ocean, and also in other contexts
More informationInstabilities and Basic Convection
Instabilities and Basic Convection Buoyant Instability (gravity is restoring force) Assume a stationary incompressible fluid (like water), so that ρ = ρ 0 + ρ/ z z and also it is in hydrostatic equilibrium
More informationReynolds Averaging. Let u and v be two flow variables (which might or might not be velocity components), and suppose that. u t + x uv ( ) = S u,
! Revised January 23, 208 7:7 PM! Reynolds Averaging David Randall Introduction It is neither feasible nor desirable to consider in detail all of the small-scale fluctuations that occur in the atmosphere.
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 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 informationUniversity of Bristol - Explore Bristol Research. Link to publication record in Explore Bristol Research PDF-document.
Dobra, T., Lawrie, A., & Dalziel, S. B. (2016). Nonlinear Interactions of Two Incident Internal Waves. 1-8. Paper presented at VIIIth International Symposium on Stratified Flows, San Diego, United States.
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 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 informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
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 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 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 informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationATMOSPHERIC AND OCEANIC FLUID DYNAMICS
ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-scale Circulation G E O F F R E Y K. V A L L I S Princeton University, New Jersey CAMBRIDGE UNIVERSITY PRESS An asterisk indicates more advanced
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 informationQuasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States
Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States See CalTech 2005 paper on course web site Free troposphere assumed to have moist adiabatic lapse rate (s* does
More information( ) = 1005 J kg 1 K 1 ;
Problem Set 3 1. A parcel of water is added to the ocean surface that is denser (heavier) than any of the waters in the ocean. Suppose the parcel sinks to the ocean bottom; estimate the change in temperature
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationDynamics of Strongly Nonlinear Internal Solitary Waves in Shallow Water
Dynamics of Strongly Nonlinear Internal Solitary Waves in Shallow Water By Tae-Chang Jo and Wooyoung Choi We study the dynamics of large amplitude internal solitary waves in shallow water by using a strongly
More informationSymmetric Instability and Rossby Waves
Chapter 8 Symmetric Instability and Rossby Waves We are now in a position to begin investigating more complex disturbances in a rotating, stratified environment. Most of the disturbances of interest are
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More informationDynamic Meteorology: lecture 2
Dynamic Meteorology: lecture 2 Sections 1.3-1.5 and Box 1.5 Potential temperature Radiatively determined temperature (boxes 1.1-1.4) Buoyancy (-oscillations) and static instability, Brunt-Vaisala frequency
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationThe Stable Boundary layer
The Stable Boundary layer the statistically stable or stratified regime occurs when surface is cooler than the air The stable BL forms at night over land (Nocturnal Boundary Layer) or when warm air travels
More information6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) instability
1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) instability in a porous layer 6-6-Lapwood.tex May 11, 3 Nield & Bejan, Chapter 6 Convection
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationStability of meridionally-flowing grounded abyssal currents in the ocean
Advances in Fluid Mechanics VII 93 Stability of meridionally-flowing grounded abyssal currents in the ocean G. E. Swaters Applied Mathematics Institute, Department of Mathematical & Statistical Sciences
More informationThermodynamics Review [?] Entropy & thermodynamic potentials Hydrostatic equilibrium & buoyancy Stability [dry & moist adiabatic]
Thermodynamics Review [?] Entropy & thermodynamic potentials Hydrostatic equilibrium & buoyancy Stability [dry & moist adiabatic] Entropy 1. (Thermodynamics) a thermodynamic quantity that changes in a
More informationLINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION
LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION É. FONTANA 1, E. MANCUSI 1,2, A. A. ULSON DE SOUZA 1, S. M. A. GUELLI U. SOUZA 1 1 Universidade Federal
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 informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationPAPER 345 ENVIRONMENTAL FLUID DYNAMICS
MATHEMATICAL TRIPOS Part III Monday, 11 June, 2018 9:00 am to 12:00 pm PAPER 345 ENVIRONMENTAL FLUID DYNAMICS Attempt no more than THREE questions. There are FOUR questions in total. The questions carry
More informationLinear stability analysis of a dual viscous layer
1 Appendix to: A hydrodynamic instability is used to create aesthetically appealing patterns in painting Sandra Zetina 1, Francisco A. Godínez 2, Roberto Zenit 3, 1 Instituto de Investigaciones Estéticas,
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationThe Equations of Motion in a Rotating Coordinate System. Chapter 3
The Equations of Motion in a Rotating Coordinate System Chapter 3 Since the earth is rotating about its axis and since it is convenient to adopt a frame of reference fixed in the earth, we need to study
More informationNon-local Momentum Transport Parameterizations. Joe Tribbia NCAR.ESSL.CGD.AMP
Non-local Momentum Transport Parameterizations Joe Tribbia NCAR.ESSL.CGD.AMP Outline Historical view: gravity wave drag (GWD) and convective momentum transport (CMT) GWD development -semi-linear theory
More informationAST242 LECTURE NOTES PART 5
AST242 LECTURE NOTES PART 5 Contents 1. Waves and instabilities 1 1.1. Sound waves compressive waves in 1D 1 2. Jeans Instability 5 3. Stratified Fluid Flows Waves or Instabilities on a Fluid Boundary
More informationLogarithmic velocity profile in the atmospheric (rough wall) boundary layer
Logarithmic velocity profile in the atmospheric (rough wall) boundary layer P =< u w > U z = u 2 U z ~ ε = u 3 /kz Mean velocity profile in the Atmospheric Boundary layer Experimentally it was found 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 information14. Instability of Superposed Fluids
4. Instability of Superposed Fluids Figure 4.: Wind over water: A layer of fluid of density ρ + moving with relative velocity V over a layer of fluid of density ρ. Define interface: h(x,y,z) = z η(x,y)
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
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 information1 The Richardson Number 1 1a Flux Richardson Number b Gradient Richardson Number c Bulk Richardson Number The Obukhov Length 3
Contents 1 The Richardson Number 1 1a Flux Richardson Number...................... 1 1b Gradient Richardson Number.................... 2 1c Bulk Richardson Number...................... 3 2 The Obukhov
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 informationNonlinear Balance on an Equatorial Beta Plane
Nonlinear Balance on an Equatorial Beta Plane David J. Raymond Physics Department and Geophysical Research Center New Mexico Tech Socorro, NM 87801 April 26, 2009 Summary Extension of the nonlinear balance
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
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 informationStructure Formation and Particle Mixing in a Shear Flow Boundary Layer
Structure Formation and Particle Mixing in a Shear Flow Boundary Layer Matthew Palotti palotti@astro.wisc.edu University of Wisconsin Center for Magnetic Self Organization Ellen Zweibel University of Wisconsin
More informationConvective instability
Chapter 2 Convective instability 2.1 Introduction The simplest problem we will explore is that of convective instability. The mechanism is intuitively straightforward: dense fluid sinks under the action
More informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in 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 information( ) where the phase! is given by! = kx + mz!"t. We also know
GFD I, Final Exam Solutions 3/7/1 Parker MacCready 1.(a) The expression for the pressure perturbation is found from the vertical momentum equation: Z-MOM w t! 1! p' z b which may be rearranged to give:
More informationCorrection of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves
PHYSICS OF FLUIDS 19, 082105 2007 Correction of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves J. C. Padrino and D. D. Joseph Aerospace Engineering and Mechanics
More informationInstability of an Axisymmetric Vortex in a Stably Stratified, Rotating Environment
Theoret. Comput. Fluid Dynamics (1998) 11: 35 322 Theoretical and Computational Fluid Dynamics Springer-Verlag 1998 Instability of an Axisymmetric Vortex in a Stably Stratified, Rotating Environment W.D.
More information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More informationInternal Wave Generation and Scattering from Rough Topography
Internal Wave Generation and Scattering from Rough Topography Kurt L. Polzin Corresponding author address: Kurt L. Polzin, MS#21 WHOI Woods Hole MA, 02543. E-mail: kpolzin@whoi.edu Abstract Several claims
More information