modeling orographic flows on unstructured meshes

15
Modeling orographic flows on unstructured meshes Piotr Smolarkiewicz, National Center for Atmospheric Research*, USA; Joanna Szmelter, Loughborough University, United Kingdom * The National Center for Atmospheric Research is supported by the National The National Center for Atmospheric Research is supported by the National Science Foundation Science Foundation

Upload: joanna

Post on 06-Jan-2016

27 views

Category:

Documents


1 download

DESCRIPTION

Modeling orographic flows on unstructured meshes Piotr Smolarkiewicz , National Center for Atmospheric Research*, USA; Joanna Szmelter , Loughborough University, United Kingdom. * The National Center for Atmospheric Research is supported by the National Science Foundation. - PowerPoint PPT Presentation

TRANSCRIPT

Page 1: Modeling  orographic  flows on unstructured meshes

Modeling orographic flows on unstructured meshes

Piotr Smolarkiewicz, National Center for Atmospheric Research*, USA; Joanna Szmelter, Loughborough University, United Kingdom

**The National Center for Atmospheric Research is supported by the National Science The National Center for Atmospheric Research is supported by the National Science Foundation Foundation

Page 2: Modeling  orographic  flows on unstructured meshes

NCAR

The twofold aim of this talk:

highlight the progress with the development of a nonhydrostatic, soundproof, unstructured-mesh model for atmospheric flows;

assess the accuracy of unstructured-mesh discretization relative to established structured-grid methods for orographic flows;

Page 3: Modeling  orographic  flows on unstructured meshes

NCAR

EULAG’s key features(www.mmm.ucar.edu/eulag/)

• A suit of governing PDEs;

numerical laboratory

• Conservative, nonoscillatory,

forward in time (NFT) semi-

implicit numerics

• Robust elliptic solver; “exact

projection”

• Static /dynamic grid stretching

with 2nd order accuracy

Cloud cover and precipitation over the Alps

Gray and blue volumes denote cloud water and rainfall, respectively.

Page 4: Modeling  orographic  flows on unstructured meshes

NCARUnstructured-mesh framework for

atmospheric flows Smolarkiewicz Szmelter, pubs in JCP, IJNMF,Comp. Fluids, 2005-2011

• Differential manifolds formulation

• Finite-volume NFT numerics with a fully unstructured spatial discretization, heritage of EULAG and its predecessors (A = MPDATA)

• Focus (so far) on wave phenomena across a range of scales and Mach, Froude & Rossby numbers

Page 5: Modeling  orographic  flows on unstructured meshes

NCARThe edge-based discretisation

Edges Dual mesh, finite volumes

Page 6: Modeling  orographic  flows on unstructured meshes

NCAR

Nonhydrostatic Boussinesq mountain waveSzmelter & Smolarkiewicz , Comp. Fluids, 2011

Comparison with the EULAG’s results and the linear theories (Smith 1979, Durran 2003): 3% in wavelength; 8% in propagation angle; wave amplitude loss 7% over 7 wavelenghts

NL/Uo = 2.4

Page 7: Modeling  orographic  flows on unstructured meshes

NCAR

For [anelastic, pseudo-incompressible]:

Soundproof generalizations Smolarkiewicz & Szmelter, Acta Geophysica, 2011

Page 8: Modeling  orographic  flows on unstructured meshes

NCAR

Numerics:

Page 9: Modeling  orographic  flows on unstructured meshes

NCAR

Non-Boussinesq amplification and breakingof deep stratospheric gravity wave

Prusa et al., JAS 1996; Smolarkiewicz & Margolin, Atmos. Ocean,1997; Klein, Ann. Rev. Fluid Dyn., 2010

NL/Uo ≈ 1 , Fr ≈ 1.6;o = 2 km Hρ A(H/2)=10ho = o

isothermal reference profiles ; Hθ = 3.5 Hρ

t=1h

t=1.5h

t=2h

EULAG “reference” solution using terrain-following coordinates

Page 10: Modeling  orographic  flows on unstructured meshes

NCARunstructured-mesh “EULAG” solution, using mesh (left) mimicking terrain-following coordinates, and (right) a fully unstructured mesh

t=1.5h

Pseudo-incmpressible equations (Durran 1989)

Anelastic equations (Lipps & Hemler 1982, Lipps 1990)

Page 11: Modeling  orographic  flows on unstructured meshes

NCAR

Hθ Hρ (RE: Achatz et al., JFM, 2010)

t=0.75h

ln θ

Page 12: Modeling  orographic  flows on unstructured meshes

NCARA global hydrostatic sound-proof model

Isentropic model:

Isosteric/isopycninc model: = ρ-1 , =p

,

M ≡ gh +

(Szmelter & Smolarkiewicz, J. Comput. Phys. 2010)

Page 13: Modeling  orographic  flows on unstructured meshes

NCAR

Fr=2

Fr=1

Fr=0.5

Stratified (mesoscale) flow past an isolated hillon a reduced planet

4 hours

Hunt & Snyder J. Fluid Mech. 1980; Smolar. & Rotunno, J. Atmos. Sci. 1989 ; Wedi & Smolar., QJR 2009

Page 14: Modeling  orographic  flows on unstructured meshes

NCAR

Smith, Advances in Geophys 1979; Hunt, Olafsson & Bougeault, QJR 2001

h=8-1o

Page 15: Modeling  orographic  flows on unstructured meshes

NCAR Conclusions:

Unstructured-mesh discretization sustains the accuracy of

structured-grid discretization and offers full flexibility in spatial resolution.

Future atmospheric models may blend structured/unstructured discretization techniques.