Yue-Kin Tsang School of Mathematics and Statistics University of St - - PowerPoint PPT Presentation

yue kin tsang
SMART_READER_LITE
LIVE PREVIEW

Yue-Kin Tsang School of Mathematics and Statistics University of St - - PowerPoint PPT Presentation

Ellipsoidal vortices beyond the quasi-geostrophic approximation Yue-Kin Tsang School of Mathematics and Statistics University of St Andrews David G .Dritschel and Jean N. Reinaud Rotating, Continuously Stratified Flows D u u = 1 (


slide-1
SLIDE 1

Ellipsoidal vortices beyond the quasi-geostrophic approximation

Yue-Kin Tsang

School of Mathematics and Statistics University of St Andrews David G .Dritschel and Jean N. Reinaud

slide-2
SLIDE 2

Rotating, Continuously Stratified Flows

D u Dt + f0 ˆ k × u = − 1 ρ0 ∇Φ + bˆ k Db Dt + N2w = 0 ∇ · u = 0

ρ( x) = ρ0+ ¯ ρ(z)+ρ′( x) b = −ρ−1

0 g ρ′

  • ω = ∇ ×

u

rotating: f-plane approximation, Ω = 1

2 f0 ˆ

k stratified: constant buoyancy frequency, N2 = −ρ−1

0 g d ¯ ρ dz

σ = f0 N = 0.1 potential vorticity anomaly: Q ≡ ω + f0ˆ k f0 · ∇b + N2ˆ k N2 = 1 + q ∼ ( ω + 2 Ω) · ∇ρ

slide-3
SLIDE 3

Computational Domain and Initial Conditions

periodic three-dimensional domain: L × L × H small aspect ratio: H = σL initial conditions: an ellipsoidal volume of constant PV anomaly q0 in a near balanced state aspect ratios λ and µ: λ = a b (a < b) µ = c σ √ ab

slide-4
SLIDE 4

Numerical methods

change of variables: {b, uh} → {q, Ah}

  • Ah =

ωh/f0 + ∇hb/f 2 equations of motion: Dq Dt = 0 D Ah Dt + f0ˆ k × Ah = N( Ah, q) q is materially conserved ⇒ equations can be solved efficiently by the Contour-Advective Semi-Lagrangian (CASL) algorithm (Dritschel & Viúdez 2003 , JFM 488, 123-150)

slide-5
SLIDE 5

Beyond the QG approximation

In the limit of asymptotically strong rotation and stratification, geostrophic-hydrostatic balance leads to the quasi-geostrophic (QG) approximation DqQG Dt = 0

  • Ah : a measure of the leading order imbalance

in contrast to the QG system, the ellipsoid is not an exact solution to the full non-hydrostatic system extend previous works on ellipsoidal vortices in the QG approximation (much less computational intensive) problem parameters: q0, λ and µ

slide-6
SLIDE 6

Effects of horizontal aspect ratio λ: q0 = 0.5

slide-7
SLIDE 7

Effects of vertical aspect ratio µ: q0 = 0.5

slide-8
SLIDE 8

Various phases for q0 = 0.5

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

unstable shape oscillation tumbling quasi-stable

slide-9
SLIDE 9

Stability and vortex geometry: q0 = 0.5

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

  • blate (small µ) vortices are more stable

vortices with close to circular cross-section (large λ) are more stable

slide-10
SLIDE 10

Stability and the Rossby number q0

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = 0.5

slide-11
SLIDE 11

Stability and the Rossby number q0

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = 0.25

slide-12
SLIDE 12

Stability and the Rossby number q0

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = -

  • 0.25
slide-13
SLIDE 13

Stability and the Rossby number q0

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = -

  • 0.5
slide-14
SLIDE 14

Stability and the Rossby number q0

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = 0.5

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = 0.25

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = -

  • 0.25

0.2 0.3 0.4 0.5 0.6 0.7 0.8

λ0

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4

µ0

q0 = -

  • 0.5

For a given (λ0, µ0), stability generally decreases with q0 cyclones are more stable than anti-cyclones

slide-15
SLIDE 15

Where do the unstable vortices go?

0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4

µ λ q0 = 0.5 q0 = 0.25 q0 = −0.25 q0 = −0.5

slide-16
SLIDE 16

Vortex rotation rate Ω

0.5 1 1.5 2 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.5 1 1.5 2 −0.04 −0.03 −0.02 −0.01 0.5 1 1.5 2 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.5 1 1.5 2 2.5 −0.1 −0.08 −0.06 −0.04 −0.02

Ω µ q0 = 0.5 q0 = 0.25 q0 = −0.25 q0 = −0.5

slide-17
SLIDE 17

Summary

study the stability and evolution of an ellipsoidal vortex in a rotating stratified fluid using the full Boussinesq’s equations

  • blate vortices with almost circular cross-section are the

most stable increase in Rossby number q0 enhances stability cyclones are more stable than anti-cyclones the vortex rotation rate scales with the vertical aspect ratio