Turbulent transport of the energy in the entrainment interface layer Marta K. Kopec 1, S. P. Malinowski1, Z. P. Piotrowski2 1 nstitute 2 nstitute of Geophysics, Faculty of Physics, University of Warsaw, Poland of Meteorology and Water Management, National Research nstitute, Warsaw, Poland 5/1/216, Reading Workshop on numerical and computational methods for simulation of all-scale geophysical flows 1/25
Outline Physics of Stratocumulus Top (POST) Setup Results Summary Workshop on numerical and computational methods for simulation of all-scale geophysical flows 2/25
Outline Physics of Stratocumulus Top (POST) Setup Results Summary Workshop on numerical and computational methods for simulation of all-scale geophysical flows 3/25
Physics of Stratocumulus Top (POST) POST (Physics of Stratocumulus Top) aircraft field study off the California Coast in 28. Single aircraft investigating unbroken stratocumulus clouds (Sc) near cloud top Aircraft instrumented with a full suite of probes for measuring state parameters of the atmosphere, drop spectra, CCN, irradiances, wind velocity and turbulence Workshop on numerical and computational methods for simulation of all-scale geophysical flows 4/25
Flight TO13 date: 9-8-28 evening flight weak inversion small surface fluxes: 5 W/m2 (sensible heat) and 1 W/m2 (latent heat) no bouayancy reversal due to entrainment/mixing Workshop on numerical and computational methods for simulation of all-scale geophysical flows 5/25
Outline Physics of Stratocumulus Top (POST) Setup Results Summary Workshop on numerical and computational methods for simulation of all-scale geophysical flows 6/25
LES setup EULAG model - anelastic equations domain size: 4 4 1.2 km grid: 2 2 2.5 m time step:.3 s simulated time: 4 h 4 simulations Workshop on numerical and computational methods for simulation of all-scale geophysical flows 7/25
LES setup EULAG model - anelastic equations domain size: 4 4 1.2 km grid: 2 2 2.5 m time step:.3 s simulated time: 4 h 4 simulations ref - referential simulation (without wind shear and radiative cooling) Workshop on numerical and computational methods for simulation of all-scale geophysical flows 7/25
LES setup EULAG model - anelastic equations domain size: 4 4 1.2 km grid: 2 2 2.5 m time step:.3 s simulated time: 4 h 4 simulations ref - referential simulation (without wind shear and radiative cooling) sh - with wind shear Workshop on numerical and computational methods for simulation of all-scale geophysical flows 7/25
LES setup EULAG model - anelastic equations domain size: 4 4 1.2 km grid: 2 2 2.5 m time step:.3 s simulated time: 4 h 4 simulations ref - referential simulation (without wind shear and radiative cooling) sh - with wind shear rc - with radiative cooling Workshop on numerical and computational methods for simulation of all-scale geophysical flows 7/25
LES setup EULAG model - anelastic equations domain size: 4 4 1.2 km grid: 2 2 2.5 m time step:.3 s simulated time: 4 h 4 simulations ref - referential simulation (without wind shear and radiative cooling) sh - with wind shear rc - with radiative cooling shrc - with wind shear and radiative cooling. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 7/25
LES initial profiles potential temperature, velocity components and water vapour mixing ratio Workshop on numerical and computational methods for simulation of all-scale geophysical flows 8/25
LES radiative cooling ( Frad = R F e Rz kqc dz z < zi z kqc dz 1/2 F e +C (z zi ) z zi Workshop on numerical and computational methods for simulation of all-scale geophysical flows 9/25
Outline Physics of Stratocumulus Top (POST) Setup Results Summary Workshop on numerical and computational methods for simulation of all-scale geophysical flows 1/25
Turbulent kinetic energy (TKE) e =.5(u 2 + v 2 + w 2 ) Workshop on numerical and computational methods for simulation of all-scale geophysical flows 11/25
Turbulent kinetic energy (TKE) e =.5(u 2 + v 2 + w 2 ) TKE is larger in cases with radiative cooling. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 11/25
Turbulent kinetic energy (TKE) e =.5(u 2 + v 2 + w 2 ) TKE is larger in cases with radiative cooling. TKE profile has two local maximum near cloud top in wind shear cases. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 11/25
TKE budget e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv buoyant production/ consumption Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv buoyant production/ consumption shear production/ loss Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv turbulent transport buoyant production/ consumption shear production/ loss Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv turbulent transport buoyant production/ consumption shear production/ loss pressure correlation Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget tendency advection e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv turbulent transport buoyant production/ consumption dissipation shear production/ loss pressure correlation Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget e g u v w e 1 w p e +w = w θv u w +v w +. t ρ θv buoyant production/ consumption shear production/ loss Workshop on numerical and computational methods for simulation of all-scale geophysical flows 12/25
TKE budget - buoyant and shear production Workshop on numerical and computational methods for simulation of all-scale geophysical flows 13/25
TKE budget - buoyant and shear production Radiative cooling prevents decoupling. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 13/25
TKE budget - buoyant and shear production Workshop on numerical and computational methods for simulation of all-scale geophysical flows 14/25
Richardson number Flux Richardson number (dots) Rf = g w θv θv v u w u + v w Gradient Richardson number (grey solid lines) Ri = g θv θv 2 u 2 + v Workshop on numerical and computational methods for simulation of all-scale geophysical flows 15/25
TKE budget e e g u v w e 1 w p +w = w θv u w +v w +. t ρ θv turbulent transport pressure correlation Workshop on numerical and computational methods for simulation of all-scale geophysical flows 16/25
TKE budget - turbulent transport and pressure correlation Workshop on numerical and computational methods for simulation of all-scale geophysical flows 17/25
TKE budget - turbulent transport and pressure correlation Workshop on numerical and computational methods for simulation of all-scale geophysical flows 18/25
TKE - vertical transport Workshop on numerical and computational methods for simulation of all-scale geophysical flows 19/25
TKE - vertical transport Stronger vertical transport of TKE in radiative cooling cases. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 19/25
TKE - vertical transport Stronger vertical transport of TKE in radiative cooling cases. Upward transport of TKE ends 5 m below cloud top. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 19/25
Comparison with measurements - layers Workshop on numerical and computational methods for simulation of all-scale geophysical flows 2/25
Comparison with measurements TSL CTMSL m SHRC m SHRC Thickness [m] 14.3 ± 14.3 19.7 ± 27.7 74.2 ± 35.5 54.3 ± 47. TKE dissipation rate [m2 s 3 1 3 ].32 ±.92 5.3 ±.33.85 ±.45 1.22 ± 1.22 Shear [s 1 ].11 ±.6.7 ±.3.5 ±.2.3 ±.1 Workshop on numerical and computational methods for simulation of all-scale geophysical flows 21/25 Corrsin scale [m].59 ±.45 3.62 ±.6 3.6 ± 1.72 4.5 ±.1
Comparison with measurements TSL CTMSL m SHRC m SHRC Thickness [m] 14.3 ± 14.3 19.7 ± 27.7 74.2 ± 35.5 54.3 ± 47. TKE dissipation rate [m2 s 3 1 3 ].32 ±.92 5.3 ±.33.85 ±.45 1.22 ± 1.22 Shear [s 1 ].11 ±.6.7 ±.3.5 ±.2.3 ±.1 n simulation TSL was much more wavier. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 21/25 Corrsin scale [m].59 ±.45 3.62 ±.6 3.6 ± 1.72 4.5 ±.1
Comparison with measurements TSL CTMSL m SHRC m SHRC Thickness [m] 14.3 ± 14.3 19.7 ± 27.7 74.2 ± 35.5 54.3 ± 47. TKE dissipation rate [m2 s 3 1 3 ].32 ±.92 5.3 ±.33.85 ±.45 1.22 ± 1.22 Shear [s 1 ].11 ±.6.7 ±.3.5 ±.2.3 ±.1 n simulation TSL was much more wavier. Turbulence in both layers is anisotropic. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 21/25 Corrsin scale [m].59 ±.45 3.62 ±.6 3.6 ± 1.72 4.5 ±.1
Comparison with measurements TSL CTMSL m SHRC m SHRC Thickness [m] 14.3 ± 14.3 19.7 ± 27.7 74.2 ± 35.5 54.3 ± 47. TKE dissipation rate [m2 s 3 1 3 ].32 ±.92 5.3 ±.33.85 ±.45 1.22 ± 1.22 Shear [s 1 ].11 ±.6.7 ±.3.5 ±.2.3 ±.1 Corrsin scale [m].59 ±.45 3.62 ±.6 3.6 ± 1.72 4.5 ±.1 n simulation TSL was much more wavier. Turbulence in both layers is anisotropic. Corsin scale in TSL in simulation was much larger than in the measurements, in CTMSL was comparable. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 21/25
Comaprison with measurements TKE production by shear near cloud top has similar shape in measurements and in SHRC simulation. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 22/25
Outline Physics of Stratocumulus Top (POST) Setup Results Summary Workshop on numerical and computational methods for simulation of all-scale geophysical flows 23/25
Summary Two adjacent to cloud top layers were found (TSL and CTMSL). Both layers are turbulent. Turbulence in both layers is anisotropic. TKE transport in cloud top region is downward. TKE production by shear near cloud top has similar shape in measurements and in SHRC simulation. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 24/25
Thank you for your attention. Workshop on numerical and computational methods for simulation of all-scale geophysical flows 25/25