A Rotational View of the Atmoshere Chapter 4
Total Page:16
File Type:pdf, Size:1020Kb
A Rotational View of the Atmoshere Chapter 4 Paul A. Ullrich [email protected] Part 2: The Vorticity Equation Vorticity: Key Questions Question: How is positive / negative vorticity generated? Question: How do we describe the time rate of change of vorticity? Question: How do we describe conservation of vorticity? Is vorticity actually conserved following the fow? Question: What is the role of the Earth’s rotation? Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation Start from the scaled horizontal momentum equation in z coordinates, no viscosity: Du 1 @p = + fv Dt −⇢ @x Dv 1 @p = fu Dt −⇢ @y − Take derivatives on both sides and compute: @ Dv @ @v = + u v @ Dv @ Du @x Dt @x @t · r ✓ ◆ ✓ ◆ @x Dt − @y Dt @ Du @ @u ✓ ◆ ✓ ◆ = + u u @y Dt @y @t · r ✓ ◆ ✓ ◆ Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation @ @u @ @ 1 @p @ (u u)=+ (fv) −@y @t − @y · r @y ⇢ @x − @y ✓ ◆ @ @v @ @ 1 @p @ + + (u v)= (fu) @x @t @x · r −@x ⇢ @y − @x ✓ ◆ A rather long derivaon @⇣ @u @v + u ⇣ +(⇣ + f) + @t · r @x @y ✓ ◆ @w @v @w @u @f @ 1 @p @ 1 @p + + v = + @x @z − @y @z @y −@x ⇢ @y @y ⇢ @x ✓ ◆ ✓ ◆ ✓ ◆ Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation @⇣ @u @v + u ⇣ +(⇣ + f) + @t · r @x @y ✓ ◆ @w @v @w @u @f @ 1 @p @ 1 @p + + v = + @x @z − @y @z @y −@x ⇢ @y @y ⇢ @x ✓ ◆ ✓ ◆ ✓ ◆ Vorticity Equation D @w @v @w @u (⇣ + f)= (⇣ + f)( u) ( ↵ p) k Dt − rh · − @x @z − @y @z − r ⇥r · ✓ ◆ 1 With specifc volume ↵ = ⇢ Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation What are these terms? D @w @v @w @u (⇣ + f)= (⇣ + f)( u) ( ↵ p) k Dt − rh · − @x @z − @y @z − r ⇥r · ✓ ◆ Tilting or twisting Divergence term Solenoidal term term Rate of change of the absolute vorticity following the motion Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation Take a closer look at the divergence term… (⇣ + f)( u) − rh · Horizontal Absolute vorticity divergence • Remember that divergence is related to vertical motion (column integrated divergence gives the local vertical velocity) • We now see that it is also related to changes in vorticity… Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vertical Velocity Obtained from horizontal divergence Continuity Equation in @u @v @! + + =0 Pressure Coordinates @x @y @p ✓ ◆p ⇥u ⇥v dω = + dp − ⇥x ⇥y ✓ ◆p ω @ p=0 p=0 ⇥u ⇥v dω = + dp − ⇥x ⇥y Zω @ ps Zps ✓ ◆p p=0 ⇥u ⇥v ω(p = 0) ω(p )= + dp − s − ⇥x ⇥y Zps ✓ ◆p ps ω(p )= ( u)dp s − ⇥p · Zp=0 Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vorticity and Divergence Vertical wind is related to the divergence of the horizontal wind. (which requires an ageostrophic component to the wind) Changes in vorticity are partially driven by divergence of the horizontal wind. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vorticity and Divergence This fow feld is purely vortical (zero divergence) Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vorticity and Divergence This fow feld is purely divergent (zero vorticity) Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vorticity and Divergence Question: What is the effect of divergence on vorticity? We can think about this like angular momentum Paul Ullrich Introduction to Atmospheric Dynamics March 2014 A Spinning Skater Motion is in the (x,y) plane Axis of rotation is in the vertical plane Paul Ullrich Introduction to Atmospheric Dynamics March 2014 A Spinning Skater Motion is in the (x,y) plane Question: A skater that is spinning with a given radius then brings in her arms. What happens to her angular velocity? Paul Ullrich Introduction to Atmospheric Dynamics March 2014 A Spinning Skater Recall: Angular momentum Motion is in the is a conserved quantity. Its (x,y) plane magnitude is given by L = rm! | | Answer: Since her radius decreases, conservation of angular momentum says that her angular velocity must increase. Question: How is this like divergence? Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Vorticity and Divergence Question: What is the effect of divergence on vorticity? Answer: A divergent fow must lead to a decrease in angular velocity. A convergent fow must lead to an increase in angular velocity. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation Divergence term only D (⇣ + f)= (⇣ + f)( u) Dt − rh · Horizontal Absolute vorticity divergence Question: What happens to absolute vorticity when the fow is (a) Divergent? ( h u) > 0 r · (b) Convergent? ( u) < 0 rh · Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation What are these terms? D @w @v @w @u (⇣ + f)= (⇣ + f)( u) ( ↵ p) k Dt − rh · − @x @z − @y @z − r ⇥r · ✓ ◆ Tilting or twisting Divergence term Solenoidal term term Rate of change of the absolute vorticity following the motion Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting and Twisting Recall how relative vorticity is obtained: Relative velocity: u =(u, v, w) @w @v @u @w @v @u 3D vorticity vector: u = , , r⇥ @y − @z @z − @x @x − @y ✓ ◆ @v @u Relative vorticity: ⇣ = k ( u)= · r⇥ @x − @y ✓ ◆ Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting and Twisting Tilting or twisting term @w @v @w @u @x @z − @y @z ✓ ◆ Tilting represents the change in the vertical component of vorticity by tilting horizontal vorticity into the vertical. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting or twisting term @w @v @w @u @x @z − @y @z ✓ ◆ z Consider a long, thin cylinder of fuid aligned with the y-axis. y x Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting or twisting term @w @v @w @u @x @z − @y @z ✓ ◆ =0 > 0 z Initial axis of rotation Consider a long, thin cylinder of fuid aligned with the y-axis. y The fow satisfes v=0 and is sheared in the vertical so that the vorticity vector is along the y axis. x Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting or twisting term @w @v @w @u @x @z − @y @z ✓ ◆ =0 > 0 > 0 z Consider a long, thin The fow also satisfes cylinder of fuid dw/dy>0 leading to an aligned with the y-axis. upward tilting of the y cylindrical tube The fow satisfes v=0 and is sheared in the vertical so that the vorticity vector is along the y axis. x Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting and Twisting z Tilting or twisting term y @w @v @w @u @x @z − @y @z ✓ ◆ =0 > 0 > 0 x Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Tilting and Twisting Rotation in the (y,z) plane. Vorticity vector points along the x axis. As the wheel is turned there is a component of vorticity in the z plane. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation What are these terms? D @w @v @w @u (⇣ + f)= (⇣ + f)( u) ( ↵ p) k Dt − rh · − @x @z − @y @z − r ⇥r · ✓ ◆ Tilting or twisting Divergence term Solenoidal term term Rate of change of the absolute vorticity following the motion Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Solenoidal / Baroclinic Terms @ 1 @p @ 1 @p @↵ @p @↵ @p + = + −@x ⇢ @y @y ⇢ @x − @x @y @y @x ✓ ◆ ✓ ◆ The solenoidal / baroclinic terms capture changes in vorticity due to the alignment of surfaces of constant density and pressure. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Solenoidal / Baroclinic Terms Defnition: In a barotropic fuid density depends only on pressure. By the ideal gas law, this implies that surfaces of constant density are surfaces of constant pressure are surfaces of constant temperature. Defnition: In a baroclinic fuid density depends on pressure and temperature. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Solenoidal / Baroclinic Terms Barotropic: surfaces of constant density parallel to surfaces of constant pressure y ↵ p =0 p-2Δp r ⇥r ρ-2Δρ p-Δp ρ-Δρ ρ p A barotropic fuid does not x induce changes in vorticity Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Solenoidal / Baroclinic Terms Baroclinic: surfaces of constant density intersect surfaces of constant pressure y ↵ p =0 p-2Δp r ⇥r 6 ρ-2Δρ p-Δp ρ-Δρ ρ p A baroclinic fuid induces x a change in vorticity. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Solenoidal / Baroclinic Terms In the real atmosphere, the solenoidal terms are a primary driver of the development of low-level low pressure systems in the middle latitudes. This term is also important in understanding circulation around fronts. Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Aside: Baroclinic Fronts Note that this case has to do with the development of horizontal vorticity (perpendicular to x-z plane) z p-2Δp ρ-2Δρ p-Δp ρ-Δρ ρ p Cold air sinking Warm air rising x behind front ahead of front Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation What are these terms? D @w @v @w @u (⇣ + f)= (⇣ + f)( u) ( ↵ p) k Dt − rh · − @x @z − @y @z − r ⇥r · ✓ ◆ Tilting or twisting Divergence term Solenoidal term term Rate of change of the absolute vorticity following the motion Paul Ullrich Introduction to Atmospheric Dynamics March 2014 Advection of Vorticity D ⇥ ⇥ ⇥ ⇥ (ζ + f)= (ζ + f)+u (ζ + f)+v (ζ + f)+w (ζ + f) Dt ⇥t ⇥x ⇥y ⇥z @f @f @f =0 =0 =0 @t @x @z D @⇣ @f (⇣ + f)= + u ⇣ + v Dt @t · r @y Advection of Advection of relative vorticity planetary vorticity Paul Ullrich Introduction to Atmospheric Dynamics March 2014 The Vorticity Equation @⇣ @u @v + u ⇣ +(⇣ + f) + @t · r @x @y ✓ ◆ @w @v @w @u @f @ 1 @p @ 1 @p + + v = + @x @z − @y @z @y −@x ⇢ @y @y ⇢ @x