Introduction to EULAG (and cloud modeling in general) Wojciech W. - - PowerPoint PPT Presentation

introduction to eulag and cloud modeling in general
SMART_READER_LITE
LIVE PREVIEW

Introduction to EULAG (and cloud modeling in general) Wojciech W. - - PowerPoint PPT Presentation

Introduction to EULAG (and cloud modeling in general) Wojciech W. Grabowski NCAR, Boulder, USA SMALL-SCALE DYNAMICS: In the spirit of the Boussinesq approximation, moisture and condensate affect gas dynamics equations only though the


slide-1
SLIDE 1

Introduction to EULAG (and cloud modeling in general) Wojciech W. Grabowski NCAR, Boulder, USA

slide-2
SLIDE 2

SMALL-SCALE DYNAMICS: In the spirit of the Boussinesq approximation, moisture and condensate affect gas dynamics equations only though the buoyancy term

p=ρRT p’/p = ρ’/ρ + T’/T p~1000 hPa p’ ~ ρu2 ~ 1 hPa

slide-3
SLIDE 3

T, q and Q – thermodynamics (and much more!)

slide-4
SLIDE 4

Water vapor is a minor constituent:

mass loading is typically smaller than 1%; thermodynamic properties (e.g., specific heats etc) only slightly modified;

Suspended small particles (cloud droplets, cloud ice):

mass loading is typically smaller than a few tenths of 1%, particles are much smaller than the smallest scale of the flow; multiphase approach is not required, but sometimes used (e.g., DNS with suspended droplets, Lagrangian Cloud Model)

Precipitation (raindrops, snowflakes, graupel, hail):

mass loading can reach a few %, particles are larger than the smallest scale the flow; multiphase approach needed only for very-small-scale modeling

slide-5
SLIDE 5

Continuous medium approach: apply density as the main field variable (density of water vapor, density of cloud water, density of rainwater, etc…) In practice, mixing ratios are typically used. Mixing ratio is the ratio between the density (of water vapor, cloud water…) and the air density.

slide-6
SLIDE 6

Mixing ratios versus specific humidities…

slide-7
SLIDE 7
slide-8
SLIDE 8

And we also need equation for the temperature:

slide-9
SLIDE 9

And we also need equation for the temperature:

slide-10
SLIDE 10

Modeling of cloud microphysics: solving a system of PDEs (advection/diffusion type) coupled through the source terms:

slide-11
SLIDE 11
slide-12
SLIDE 12

A very simple (but useful) model: rising adiabatic parcel… Take a parcel from the surface and move it up… … by solving these equations.

slide-13
SLIDE 13

θk , qv

k , qc k

θk+1 , qv

k+1 , qc k+1

The nonlinear equation for Δq can be solved using the Newton-Raphson method…

slide-14
SLIDE 14

F(x) = 0 x = ?

slide-15
SLIDE 15

qv qc

Look not only on the patterns (i.e., processes), but also on specific numbers (e.g., temperature change, mixing ratios, etc).

slide-16
SLIDE 16

Invariant variables: total water, liquid water potential temperature, equivalent potential temperature.

slide-17
SLIDE 17

Adding rain or drizzle:

slide-18
SLIDE 18

We need something more complicated than a rising parcel as rain has to fall

  • ut. One possibility is to use the kinematic (prescribed flow) framework…
slide-19
SLIDE 19

Cloud water and rain (drizzle) fields after 2 hrs (almost quasi- equilibrium…)

slide-20
SLIDE 20

BabyEULAG documentation and code

slide-21
SLIDE 21
slide-22
SLIDE 22
slide-23
SLIDE 23
slide-24
SLIDE 24
slide-25
SLIDE 25
slide-26
SLIDE 26
slide-27
SLIDE 27
slide-28
SLIDE 28
slide-29
SLIDE 29
slide-30
SLIDE 30
slide-31
SLIDE 31
slide-32
SLIDE 32
slide-33
SLIDE 33
slide-34
SLIDE 34
slide-35
SLIDE 35

Time to look at and run the code…