Deutscher Wetterdienst
|
|
- Alberta Watkins
- 6 years ago
- Views:
Transcription
1 S-LC-RK3 (3-stage, 2nd order) + up5 (Wicker, Skamarock, 2002) TVD-RK3 (Shu, Osher, 1988) + up5 different small timesteps ΔT/Δt = 6 possible better efficiency for S-LC-RK3 method used in COSMO-DE, planned for COSMO-EU 34
2 Artificial horizontal diffusion 4th order dim.-less diffusion coeff.: D = K 4 Δt / Δx 4 stable + 'non-oszillating' sinus-waves: für 0 D 1/64 (LM-Doku: 0 D 1/128, Leapfrog (2 Δt)) Amplification factor for D=0.2*1/195 used in COSMO-model up to now: D = 1/ (2π 4 ) 1/195 (COSMO-EU: whole domain, COSMO-DE: only in boundary zone) k * dx planned for COSMO-EU (dx=7 km): D = 0.2 * 1/195 reduction of amplitude for (2 Δx, 2 Δy)-waves: 7% / time step Diffusion 4th order is not monotone! flux limitation (Doms, 2001) additional orography-limitation: diffus. flux=0, if slope of a coordinate plane > 250 m / Δx 35
3 Advantages of p't'-dynamics over p't-dynamics A reference state improves the calculation of horizontal pressure gradient terms in terrain following coordinates p -dynamics 1. Improved representation of T-advection in terrain-following coordinates Terms (a) and (b) cancel analytically, but not numerically 2. Better representation of buoyancy term in fast waves solver Buoyancy term alone generates an oscillation equation: using T: ω = g/c s using T': ω = ω a = acoustic cut-off frequency 36
4 point 1.): 'improved T-advection'... Idealised test case: Steady atmosphere with mountain base state: T 0, p 0 deviations from base state: T', p' 0 introduces spurious circulations! 37
5 Leapfrog Runge-Kutta old p'-t-dynamics contours: vertical velocity w isolines: potential temperature Θ 38
6 Runge-Kutta new p'-t'-dynamik Runge-Kutta old p'-t-dynamik contours: vertical velocity w isolines: potential temperature Θ 39
7 Deep / shallow atmosphere Momentum equations for deep atmosphere (spherical coordinates): shallow atmosphere approximation: r ~ a neglect terms in advection and Coriolis force additionally: introduce a hydrostatic, steady base state transformation to terrain following coordinates deep atmosphere terms are implemented since COSMO 3.21 (but only for RK) lcori_deep, ladv_deep diploma thesis R. Petrik, 2006, Univ. Leipzig White, Bromley (1995), QJRMS Davies et al. (2005), QJRMS 40
8 Time integration of Coriolis terms Old: Coriolis terms are handled like physical tendencies = Euler forward numerically unstable, nevertheless the instability is very slow, the amplification factor is 1 + i Δt f not important in weather forecast applications, but perhaps in climate runs? l_coriolis_every_rk_substep =.FALSE. in src_runge_kutta.f90 (since COSMO 4.8) New: Coriolis terms are handled like horizontal advection, i.e. tendencies are calculated and added in every RK substep stable l_coriolis_every_rk_substep =.TRUE. Not expensive (Subr. coriolis needs ~0.2% comp. time) 41
9 Test case Weisman, Klemp (1982): warm bubble in a base flow with vertical velocity shear + Coriolis force w max dx= 2 km RR precipitation distribution deep (shaded), shallow (isolines) RR deep -RR shallow (shaded) 42
10 Case study (Diploma thesis R. Petrik, 2006) summary for precipitation forecast in deep, convection resolving models: additional advection terms: not relevant additional Coriolis terms: have a certain influence, but don't seem to be important for COSMO-DE application could be important for simulations near the equator 43
11 fields ϕ n =(u, v, w, pp, T,...) n Physics-Dynamics-coupling Descr. of Advanced Research WRF, v. 2 (2005) Physics (I) Radiation Convection Coriolis force ('old' scheme) Turbulence Physics (I) -tendencies Δϕ n (phys I) + Δϕ n-1 (phys II) NEW! Dynamics Runge-Kutta [Δϕ (phys) + Δϕ (adv) --> fast_waves ] fields ϕ * =(u, v, w, pp, T,...)* - Δϕ n-1 (phys II) NEW! Physics (II) Cloud Physics Physics (II) -tendencies Δϕ n (phys II) fields ϕ n+1 =(u, v, w, pp, T,...) n+1 44
12 Tracer transport General conservation equation for moisture variables: ρ total density of air [kg/m 3 ] q x = ρ x /ρ mixing ratio [kg/kg] F x turbulent flux of x [kg/m 2 /s] P x sedimentation flux of x (x=r,s,g) [kg/m 2 /s] S x Sources/sinks of x [kg/m 3 /s] x=v x=c x=i x=r x=s x=g water vapour cloud water cloud ice rain, v sedi ~ 5 m/s snow, v sedi ~ 1 m/s graupel 45
13 Tracer transport Semi-Lagrange advection scheme lsl_adv_qx =.true. Bott advection scheme lsl_adv_qx =.false. COSMO-EU: qv, qc: centered diff. 2nd order qi: 2nd ord. flux-form advection scheme qr, qs: Semi-Lagrange (tri-linear interpol.) planned: all qx: Semi-Lagrange (tri-cubic interpol.) COSMO-DE: Courant-number-independent (CNI)-advection for all variables - Bott (1989) (2., 4. order), in conservation form - Semi-Lagrange (tri-cubic interpol.) 46
14 Semi-Lagrange-Advection advection eq. (1-dim.) rewritten as step 1: calculation of backward trajectory in principle any ODE-solver can be used (here: 2nd order) Staniforth, Côté (1991) MWR Baldauf, Schulz (2004) COSMO-Newsl. Förstner, Baldauf, Seifert (2006) COSMO-Newsl. 47
15 Semi-Lagrange Advection 2nd step: Interpolation from neighbouring points i,j,k = -1,0 for tri-linear interpol. 8 grid points i,j,k = -2,...,1 for tri-cubic interpol. 64 grid points x,y,z = position in the grid cell (from backtrajectory calculation) q i,j,k = grid point value of q linear weighting polynomials: cubic weighting polynomials: 48
16 properties of Semi-Lagrange advection + unconditionally stable (at least for u=const.) + fully multi-dimensional scheme (no directional splitting necessary) + increased efficiency if used for many tracers (calculation of backtrajectory only once) + can be implemented also in unstructured grids + no non-linear instability if used for velocity advection - non-conserving scheme (but conservation properties are not bad without clipping) - multi-cubic interpolation not positive definite clipping; multiplicative filling (multi-linear interpolation monotone, but highly diffusive) 49
17 Multiplicative Filling (Rood, 1987) clipped values are globally summed and distributed over the whole field easy fast but only global conservation Problem of reproducibility: a sum of 'real' (=floating point) numbers is not associative: (a + b) + c a + ( b + c ) solution: a sum of integer numbers is associative map the Real number space to the Integer number space 50
18 SL-Advection with v = (100, 50, 0) m/s, tri-cubic interpolation Deutscher Clipping Wetterdienst Mult. Filling 51
19 SL-Advection with solid body rotation (200 time steps), tri-cubic-interpolation Deutscher Clipping Wetterdienst Mult. Filling 52
20 Bott advection scheme 1.) finite volume discretization: with flux (case u j+1/2 > 0 ): to evaluate the integral: interpolation of ϕ(x) from ϕ j n 2.) limitation of fluxes to get positive definiteness 3.) use for Courant numbers C=u dt/dx > 1: apply this idea only to the fractional part of the Courant number 'fractional flux' (the 'integer flux' is a simple summing up of the intermediate cell values) 4.) directional splitting ('x-y-z') using the mass-conservation idea of Easter (1993) 5.) Strang-splitting ('x-y-z' + 'z-y-x') to achieve 2nd order in time Bott (1989) MWR Easter (1993) MWR Skamarock (2004, 2006) MWR Förstner, Baldauf, Seifert (2006) COSMO-Newsl. 53
21 Properties of the Bott advection scheme + in conservative form + positive definite ('flux limitation') + useable for Courant-numbers C > 1 if using a splitting of fluxes into integer and non-integer part (Skamarock, 2004, 2006) (but loss of conservation) - one-dimensional scheme directional splitting necessary use splitting idea of Easter (1993), Skamarock (2004, 2006) mass consistency: parallel computation of the density equation necessary - Strang splitting necessary to get a 2nd order temporal scheme 54
22 Idealised tests of the tracer advection schemes Tracer initialisation: '3D cone' (figure) or '3D sphere' terrain following coordinate system over a series of hills 55
23 Test: constant velocity v=(10, 10, 0) m/s, Cx=Cy=0.107 tracer: '3D cone'; difference to analytic solution Semi-Lagrange Bott 56
24 Test: constant velocity v=(70, 70, 0) m/s, Cx+Cy=1.5 tracer: '3D sphere'; difference to analytic solution Semi-Lagrange Bott 57
25 Test LeVeque (1996); tracer: '3D sphere'; difference to analytic solution Semi-Lagrange Bott 58
26 Transport of Tracer in a Real Case Flow Field init Bott (2 nd ) Flux Form - DIV + Clipping semi- Lagrange (tri-cubic) + Clipping Bott (2 nd ) Conserv. Form 59
27 Problems found with Bott (1989)-scheme in the meanwhile: 1.) Directional splitting of the scheme: Parallel Marchuk-splitting of conservation equation for density can lead to a complete evacuation of cells Solution: Easter (1993), Skamarock (2004, 2006), mass-consistent splitting 2.) Strang-splitting ( 'x-y-z' and 'z-y-x' in 2 time steps) produces 2*dt oscillations Solution: proper Strang-Splitting ('x-y-2z-y-x') in every time step solves the problem, but nearly doubles the computation time 3.) metric terms prevent the scheme to be properly mass conserving <-- Schär test case of an unconfined jet and tracer=1 initialisation (remark: exact mass conservation is already violated by the 'flux-shifting' to make the Bott-scheme Courant-number independent) 60
28 Diagnostic tools Possibility to calculate potential vorticity (PV) (on model levels) and curl v in COSMO 4.12 Possibility to perform volume integrations and integrations of surface fluxes for budget calculations (src_integrals.f90) Baldauf (2008) COSMO-Newsl. Other numerical stuff not mentioned here: New reference state boundary conditions bottom (physical) BC's lateral BC's upper BC's Carry out idealised test cases (exercises in the afternoon) 61
29 Vorhersage von Leewellen am Sat.-bild Vis., DLR, , 10 UTC COSMO-DE-Vorhersage 62
30 Vorhersage von Leewellen am v y yz-schnitt (j=293) 63
31 Vorhersage von Leewellen am v y yz-schnitt (j=293) 64
32 Vorhersage von Leewellen am stabile Schichtung und mit der Höhe anwachsendes Windprofil führen zur Bildung einer Leewelle ('Resonanzwelle') v y yz-schnitt (j=293) 65
33 Vorhersage von Leewellen am v y yz-schnitt (j=293) 66
Mass Conserving Courant Number Independent Eulerian Advection of the Moisture Quantities for the LMK
Mass Conserving Courant Number Independent Eulerian Advection of the Moisture Quantities for the LMK Jochen Förstner, Michael Baldauf, Axel Seifert Deutscher Wetterdienst, Kaiserleistraße 29/35, 63067
More informationThe Runge-Kutta method
The Runge-Kutta method Advanced Numerics Seminar Langen, 08.03.2006 Michael Baldauf Deutscher Wetterdienst, Offenbach, Germany Outline View of LMK-Dynamics, -Numerics Time-splitting schemes General remarks
More informationDeutscher Wetterdienst
Stability Analysis of the Runge-Kutta Time Integration Schemes for the Convection resolving Model COSMO-DE (LMK) COSMO User Seminar, Langen 03.+04. March 2008 Michael Baldauf Deutscher Wetterdienst, Offenbach,
More informationStabilitätsanalyse des Runge-Kutta- Zeitintegrationsschemas für das konvektionserlaubende Modell COSMO-DE
Stabilitätsanalyse des Runge-Kutta- Zeitintegrationsschemas für das konvektionserlaubende Modell COSMO-DE DACH-Tagung 20.-24.09.200, Bonn Michael Baldauf GB FE, Deutscher Wetterdienst, Offenbach M. Baldauf
More informationPrognostic Precipitation in the Lokal-Modell (LM)
Prognostic Precipitation in the Lokal-Modell (LM) Michael Baldauf, Günther Doms German Weather Service (DWD), Offenbach, Germany COSMO-Meeting, Mailand, 22.09.2004 COSMO, Work Package 2.1.1: Task: replacement
More informationSplit explicit methods
Split explicit methods Almut Gassmann Meteorological Institute of the University of Bonn Germany St.Petersburg Summer School 2006 on nonhydrostatic dynamics and fine scale data assimilation Two common
More informationFlux Diagnosis in the COSMO model and their problems
Flux Diagnosis in the COSMO model and their problems Ronny Petrik 1 Michael Baldauf 2 Almut Gassmann 1 1 Max-Planck-institute for meteorology 2 German Weather Service 2nd March 2008 Outline Motivation
More informationThe 'Conservative Dynamical Core' priority. developments in the COSMO model
1 The 'Conservative Dynamical Core' priority project and other current dynamics developments in the COSMO model 33 rd EWGLAM and 18 th SRNWP meeting 10-13 Oct. 2011, Tallinn, Estonia M. Baldauf, U. Blahak
More informationWeather Research and Forecasting Model. Melissa Goering Glen Sampson ATMO 595E November 18, 2004
Weather Research and Forecasting Model Melissa Goering Glen Sampson ATMO 595E November 18, 2004 Outline What does WRF model do? WRF Standard Initialization WRF Dynamics Conservation Equations Grid staggering
More informationThe Advanced Research WRF (ARW) Dynamics Solver
Dynamics: Introduction The Advanced Research WRF (ARW) Dynamics Solver 1. What is a dynamics solver? 2. Variables and coordinates 3. Equations 4. Time integration scheme 5. Grid staggering 6. Advection
More informationA semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate
A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate Juan Simarro and Mariano Hortal AEMET Agencia Estatal de
More informationDevelopment of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA.
Development of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA. M. Sakamoto, J. Ishida, K. Kawano, K. Matsubayashi, K. Aranami, T. Hara, H. Kusabiraki, C. Muroi, Y. Kitamura Japan Meteorological
More informationIdealized Nonhydrostatic Supercell Simulations in the Global MPAS
Idealized Nonhydrostatic Supercell Simulations in the Global Joe Klemp, Bill Skamarock, and Sang-Hun Park National Center for Atmospheric Research Boulder, Colorado Typical characteristics: Supercell Thunderstorms
More informationFormulation and performance of the Variable-Cubic Atmospheric Model
Formulation and performance of the Variable-Cubic Atmospheric Model John McGregor CSIRO Marine and Atmospheric Research Aspendale, Melbourne Southern Hemisphere PDEs on the Sphere NCAR 11 April 2014 CSIRO
More informationConvection-Resolving Model Simulations: Process-Based Comparison of LM Results with Observations
Convection-Resolving Model Simulations: Process-Based Comparison of LM Results with Observations Jörg Trentmann, Britta Wecker, Marcus Paulat, Heini Wernli, Ulrich Corsmeier, Jan Handwerker Goal Improve
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationMOX EXPONENTIAL INTEGRATORS FOR MULTIPLE TIME SCALE PROBLEMS OF ENVIRONMENTAL FLUID DYNAMICS. Innsbruck Workshop October
Innsbruck Workshop October 29 21 EXPONENTIAL INTEGRATORS FOR MULTIPLE TIME SCALE PROBLEMS OF ENVIRONMENTAL FLUID DYNAMICS Luca Bonaventura - Modellistica e Calcolo Scientifico Dipartimento di Matematica
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More information3D-Transport of Precipitation
3D-Transport of Precipitation Almut Gassmann DWD/FE 13 October 17, 2000 1 Current diagnostic scheme for gridscale precipitation The current scheme for gridscale precipitation assumes column equilibrium
More informationMODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction
MODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction Grid point and spectral models are based on the same set of primitive equations. However, each type formulates and solves the equations
More informationHorizontal nonlinear Smagorinsky diffusion
1 Working Group on Numerical Aspects 3 Horizontal nonlinear Smagorinsky diffusion M. Baldauf and G. Zängl Deutscher Wetterdienst, Offenbach, Germany 1 Introduction The COSMO model uses several diffusion
More informationJP1.15 PHYSICAL INITIALIZATION TO INCORPORATE RADAR PRECIPITATION DATA INTO A NUMERICAL WEATHER PREDICTION MODEL (Lokal Model)
JP1.15 PHYSICAL INITIALIZATION TO INCORPORATE RADAR PRECIPITATION DATA INTO A NUMERICAL WEATHER PREDICTION MODEL (Lokal Model) M. Milan 1, F. Ament 1, V. Venema 1, A. Battaglia 1, C. Simmer 1 1 Meteorological
More informationFlorian Meier, Peter Knippertz Johannes Gutenberg-University Mainz
Sensitivity experiments with the COSMO-Model - the influence of upstream latent heating and wave amplification on a heavy precipitation event over West Africa Florian Meier, Peter Knippertz Johannes Gutenberg-University
More informationExperiments with Idealized Supercell Simulations
Latent Heat Nudging in almo: Experiments with Idealized Supercell Simulations Daniel Leuenberger and Andrea Rossa MeteoSwiss Talk Outline Introduction Supercell Storms Characteristics Reference simulation
More informationFor the operational forecaster one important precondition for the diagnosis and prediction of
Initiation of Deep Moist Convection at WV-Boundaries Vienna, Austria For the operational forecaster one important precondition for the diagnosis and prediction of convective activity is the availability
More informationUsing hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling
Using hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling Michael D. Toy NCAR Earth System Laboratory/ Advanced Study Program Boulder, Colorado, USA 9th International SRNWP
More informationMechanisms and prediction of precipitation over complex terrain
Mechanisms and prediction of precipitation over complex terrain S. Wassermann, M. Kunz, Ch. Kottmeier Institut für Meteorologie und Klimaforschung, Universität Karlsruhe (TH) / International PhD-students
More informationImpact of Turbulence on the Intensity of Hurricanes in Numerical Models* Richard Rotunno NCAR
Impact of Turbulence on the Intensity of Hurricanes in Numerical Models* Richard Rotunno NCAR *Based on: Bryan, G. H., and R. Rotunno, 2009: The maximum intensity of tropical cyclones in axisymmetric numerical
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationMS-GWaves / GWING: Towards UA-ICON A non-hydrostatic global model for studying gravity waves from troposphere to thermosphere
MS-GWaves / GWING: Towards UA-ICON A non-hydrostatic global model for studying gravity waves from troposphere to thermosphere Sebastian Borchert, Günther Zängl, Michael Baldauf (1), Guidi Zhou, Hauke Schmidt,
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationThe prognostic deep convection parametrization for operational forecast in horizontal resolutions of 8, 4 and 2 km
The prognostic deep convection parametrization for operational forecast in horizontal resolutions of 8, 4 and 2 km Martina Tudor, Stjepan Ivatek-Šahdan and Antonio Stanešić tudor@cirus.dhz.hr Croatian
More information1 INTRODUCTION 2 MODEL SETUP
8.1 INFLUENCE OF AMBIENT ENVIRONMENTAL CONDITIONS AND OROGRAPHY ON THE CHARACTERISTICS OF DEEP CONVECTIVE CELLS AS SIMULATED WITH A SOPHISTICATED TWO-MOMENT (BULK) MICROPHYSICAL SCHEME Ulrich Blahak, H.
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More information5. General Circulation Models
5. General Circulation Models I. 3-D Climate Models (General Circulation Models) To include the full three-dimensional aspect of climate, including the calculation of the dynamical transports, requires
More informationTheories on the Optimal Conditions of Long-Lived Squall Lines
Theories on the Optimal Conditions of Long-Lived Squall Lines References: Thorpe, A. J., M. J. Miller, and M. W. Moncrieff, 1982: Two -dimensional convection in nonconstant shear: A model of midlatitude
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationway and atmospheric models
Scale-consistent consistent two-way way coupling of land-surface and atmospheric models COSMO-User-Seminar 9-11 March 2009 Annika Schomburg, Victor Venema, Felix Ament, Clemens Simmer TR / SFB 32 Objective
More informationConvection Permitting Simulation of Extreme Precipitation Events Lessons Learned
Convection Permitting Simulation of Extreme Precipitation Events Lessons Learned A. F. Prein, M. Suklitsch, H. Truhetz, and A. Gobiet Wegener Center for Climate and Global Change (WEGC) and Institute for
More informationDescription of the ET of Super Typhoon Choi-Wan (2009) based on the YOTC-dataset
High Impact Weather PANDOWAE Description of the ET of Super Typhoon Choi-Wan (2009) based on the YOTC-dataset ¹, D. Anwender¹, S. C. Jones2, J. Keller2, L. Scheck¹ 2 ¹Karlsruhe Institute of Technology,
More information9D.3 THE INFLUENCE OF VERTICAL WIND SHEAR ON DEEP CONVECTION IN THE TROPICS
9D.3 THE INFLUENCE OF VERTICAL WIND SHEAR ON DEEP CONVECTION IN THE TROPICS Ulrike Wissmeier, Robert Goler University of Munich, Germany 1 Introduction One does not associate severe storms with the tropics
More informationDescription of. Jimy Dudhia Dave Gill. Bill Skamarock. WPS d1 output. WPS d2 output Real and WRF. real.exe General Functions.
WPS d1 output WPS d2 output Real and WRF Description of real.exe General Functions Jimy Dudhia Dave Gill wrf d01 input wrf d01 bdy wrf d02 input Bill Skamarock wrf.exe Real program in a nutshell Function
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationPV Generation in the Boundary Layer
1 PV Generation in the Boundary Layer Robert Plant 18th February 2003 (With thanks to S. Belcher) 2 Introduction How does the boundary layer modify the behaviour of weather systems? Often regarded as a
More information3. Midlatitude Storm Tracks and the North Atlantic Oscillation
3. Midlatitude Storm Tracks and the North Atlantic Oscillation Copyright 2006 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without permission. EFS 3/1 Review of key results
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationRecent developments of the global semi-lagrangian atmospheric model
Recent developments of the global semi-lagrangian atmospheric model Mikhail Tolstykh, Vladimir Shashkin (Inst. of Numerical Math. RAS, Hydrometcentre of Russia), Alla Yurova (Hydrometcentre of Russia)
More informationρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More informationICON. The Icosahedral Nonhydrostatic model: Formulation of the dynamical core and physics-dynamics coupling
ICON The Icosahedral Nonhydrostatic model: Formulation of the dynamical core and physics-dynamics coupling Günther Zängl and the ICON deelopment team PDEs on the sphere 2012 Outline Introduction: Main
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationSimulating orographic precipitation: Sensitivity to physics parameterizations and model numerics
Simulating orographic precipitation: Sensitivity to physics parameterizations and model numerics 2nd COPS-Meeting, 27 June 2005 Günther Zängl Overview A highly idealized test of numerical model errors
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationPreliminary results with very high resolution COSMO model for the forecast of convective events. Antonella Morgillo. Arpa-Simc.
2 Working Group on Physical Aspects 52 Preliminary results with very high resolution COSMO model for the forecast of convective events. Antonella Morgillo Arpa-Simc amorgillo@arpa.emr.it 1 Introduction
More informationM.Sc. in Meteorology. Physical Meteorology Prof Peter Lynch. Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield.
M.Sc. in Meteorology Physical Meteorology Prof Peter Lynch Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield. Climate Change???????????????? Tourists run through a swarm of pink
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationSome remarks on climate modeling
Some remarks on climate modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Selected overheads by Doug Nychka Outline Hierarchy of atmospheric modeling strategies
More informationALARO-0 at the RMI. Alex, Daan, Geert, Josette, Luc, Olivier, Rafiq, Piet( s first operational talk)
ALARO-0 at the RMI Alex, Daan, Geert, Josette, Luc, Olivier, Rafiq, Piet( s first operational talk) Operational use of Alaro-0 0 at RMIB The model (cy35t1) is run 4 times a day, at: 1. Resolution 4km (square
More informationTorben Königk Rossby Centre/ SMHI
Fundamentals of Climate Modelling Torben Königk Rossby Centre/ SMHI Outline Introduction Why do we need models? Basic processes Radiation Atmospheric/Oceanic circulation Model basics Resolution Parameterizations
More informationP1.1 THE QUALITY OF HORIZONTAL ADVECTIVE TENDENCIES IN ATMOSPHERIC MODELS FOR THE 3 RD GABLS SCM INTERCOMPARISON CASE
P1.1 THE QUALITY OF HORIZONTAL ADVECTIVE TENDENCIES IN ATMOSPHERIC MODELS FOR THE 3 RD GABLS SCM INTERCOMPARISON CASE Fred C. Bosveld 1*, Erik van Meijgaard 1, Evert I. F. de Bruijn 1 and Gert-Jan Steeneveld
More informationDeutscher Wetterdienst
Deutscher Wetterdienst Modelling the Volcanic Ash Episode: Experiences with COSMO-ART Detlev Majewski (FE1) Bernhard Vogel, Heike Vogel (KIT) Thomas Hanisch, Jochen Förstner (FE13), Ulrich Pflüger (FE15)
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationThe perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh
z = The perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh, that is, p = p h + p nh. (.1) The former arises
More informationWaVaCS summerschool Autumn 2009 Cargese, Corsica
Introduction Part I WaVaCS summerschool Autumn 2009 Cargese, Corsica Holger Tost Max Planck Institute for Chemistry, Mainz, Germany Introduction Overview What is a parameterisation and why using it? Fundamentals
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationProf. Simon Tett, Chair of Earth System Dynamics & Modelling: The University of Edinburgh
SAGES Scottish Alliance for Geoscience, Environment & Society Modelling Climate Change Prof. Simon Tett, Chair of Earth System Dynamics & Modelling: The University of Edinburgh Climate Modelling Climate
More informationThe Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science
The Fifth-Generation NCAR / Penn State Mesoscale Model (MM5) Mark Decker Feiqin Xie ATMO 595E November 23, 2004 Department of Atmospheric Science Outline Basic Dynamical Equations Numerical Methods Initialization
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationCAM-SE: Lecture I. Peter Hjort Lauritzen
CAM-SE: Lecture I Peter Hjort Lauritzen Atmospheric Modeling and Predictability Section Climate and Global Dynamics Laboratory National Center for Atmospheric Research 2nd WCRP Summer School on Climate
More informationThunderstorm-Scale EnKF Analyses Verified with Dual-Polarization, Dual-Doppler Radar Data
Thunderstorm-Scale EnKF Analyses Verified with Dual-Polarization, Dual-Doppler Radar Data David Dowell and Wiebke Deierling National Center for Atmospheric Research, Boulder, CO Ensemble Data Assimilation
More informationA Large-Eddy Simulation Study of Moist Convection Initiation over Heterogeneous Surface Fluxes
A Large-Eddy Simulation Study of Moist Convection Initiation over Heterogeneous Surface Fluxes Song-Lak Kang Atmospheric Science Group, Texas Tech Univ. & George H. Bryan MMM, NCAR 20 th Symposium on Boundary
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationEvaluation of three spatial discretization schemes with the Galewsky et al. test
Evaluation of three spatial discretization schemes with the Galewsky et al. test Seoleun Shin Matthias Sommer Sebastian Reich Peter Névir February 22, 2 Abstract We evaluate the Hamiltonian Particle Methods
More informationComparing the formulations of CCAM and VCAM and aspects of their performance
Comparing the formulations of CCAM and VCAM and aspects of their performance John McGregor CSIRO Marine and Atmospheric Research Aspendale, Melbourne PDEs on the Sphere Cambridge 28 September 2012 CSIRO
More informationA Time-Split Nonhydrostatic Atmospheric Model for Weather Research and Forecasting Applications. William C. Skamarock and Joseph B.
A Time-Split Nonhydrostatic Atmospheric Model for Weather Research and Forecasting Applications William C. Skamarock and Joseph B. Klemp National Center for Atmospheric Research, Boulder, Colorado Submitted
More informationBMeteorologische Zeitschrift, Vol.23, No. 4, (published online August 27, 2014)
BMeteorologische Zeitschrift, Vol.23, No. 4, 449 464 (published online August 27, 214) 214 The authors Open Access Article On discontinuous Galerkin approach for atmospheric flow in the mesoscale with
More informationNon-local Momentum Transport Parameterizations. Joe Tribbia NCAR.ESSL.CGD.AMP
Non-local Momentum Transport Parameterizations Joe Tribbia NCAR.ESSL.CGD.AMP Outline Historical view: gravity wave drag (GWD) and convective momentum transport (CMT) GWD development -semi-linear theory
More informationEvaluating parameterisations of subgridscale variability with satellite data
Evaluating parameterisations of subgridscale variability with satellite data Johannes Quaas Institute for Meteorology University of Leipzig johannes.quaas@uni-leipzig.de www.uni-leipzig.de/~quaas Acknowledgements
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationNWP Equations (Adapted from UCAR/COMET Online Modules)
NWP Equations (Adapted from UCAR/COMET Online Modules) Certain physical laws of motion and conservation of energy (for example, Newton's Second Law of Motion and the First Law of Thermodynamics) govern
More informationAtmospheric dynamics and meteorology
Atmospheric dynamics and meteorology B. Legras, http://www.lmd.ens.fr/legras III Frontogenesis (pre requisite: quasi-geostrophic equation, baroclinic instability in the Eady and Phillips models ) Recommended
More informationChapter 1. Introduction
Chapter 1. Introduction In this class, we will examine atmospheric phenomena that occurs at the mesoscale, including some boundary layer processes, convective storms, and hurricanes. We will emphasize
More informationImplementation and Evaluation of WSR-88D Reflectivity Data Assimilation for WRF-ARW via GSI and Cloud Analysis. Ming Hu University of Oklahoma
Implementation and Evaluation of WSR-88D Reflectivity Data Assimilation for WRF-ARW via GSI and Cloud Analysis Ming Hu University of Oklahoma 1. Previous work and Goal of visiting Previous work: Radar
More informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationCONVECTIVE CLOUD MICROPHYSICS IN A HIGH-RESOLUTION NWP MODEL
CONVECTIVE CLOUD MICROPHYSICS IN A HIGH-RESOLUTION NWP MODEL J. Trentmann 1, A. Seifert 2, H. Wernli 1 1 Institute for Atmospheric Physics, Johannes Gutenberg University Mainz, Germany 2 German Weather
More informationConvection-Resolving NWP with WRF. Section coordinator Ming Xue University of Oklahoma
Convection-Resolving NWP with WRF Section coordinator Ming Xue University of Oklahoma Convection-resolving NWP Is NWP that explicitly treats moist convective systems ranging from organized MCSs to individual
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationDaniel J. Jacob, Models of Atmospheric Transport and Chemistry, 2007.
1 0. CHEMICAL TRACER MODELS: AN INTRODUCTION Concentrations of chemicals in the atmosphere are affected by four general types of processes: transport, chemistry, emissions, and deposition. 3-D numerical
More informationModel equations for planetary and synoptic scale atmospheric motions associated with different background stratification
Model equations for planetary and synoptic scale atmospheric motions associated with different background stratification Stamen Dolaptchiev & Rupert Klein Potsdam Institute for Climate Impact Research
More informationAtmospheric Dynamics with Polyharmonic Spline RBFs
Photos placed in horizontal position with even amount of white space between photos and header Atmospheric Dynamics with Polyharmonic Spline RBFs Greg Barnett Sandia National Laboratories is a multimission
More informationWRF Model Simulated Proxy Datasets Used for GOES-R Research Activities
WRF Model Simulated Proxy Datasets Used for GOES-R Research Activities Jason Otkin Cooperative Institute for Meteorological Satellite Studies Space Science and Engineering Center University of Wisconsin
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationQuasi-Geostrophic Implications
Chapter 10 Quasi-Geostrophic Implications When you look at a weather chart with all its isolines and plotted data, you need a framework upon which to interpret what you see. Quasi-geostrophic theory provides
More informationt tendency advection convergence twisting baroclinicity
RELATIVE VORTICITY EQUATION Newton s law in a rotating frame in z-coordinate (frictionless): U + U U = 2Ω U Φ α p U + U U 2 + ( U) U = 2Ω U Φ α p Applying to both sides, and noting ω U and using identities
More informationMeteorology 6150 Cloud System Modeling
Meteorology 6150 Cloud System Modeling Steve Krueger Spring 2009 1 Fundamental Equations 1.1 The Basic Equations 1.1.1 Equation of motion The movement of air in the atmosphere is governed by Newton s Second
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationz g + F w (2.56) p(x, y, z, t) = p(z) + p (x, y, z, t) (2.120) ρ(x, y, z, t) = ρ(z) + ρ (x, y, z, t), (2.121)
= + dw dt = 1 ρ p z g + F w (.56) Let us describe the total pressure p and density ρ as the sum of a horizontally homogeneous base state pressure and density, and a deviation from this base state, that
More informationNon-hydrostatic sound-proof equations of motion for gravity-dominated compressible fl ows. Thomas Dubos
Non-hydrostatic sound-proof of for gravity-dominated compressible fl ows Thomas Dubos Laboratoire de Météorologie Dynamique/IPSL, École Polytechnique, Palaiseau, France Fabrice Voitus Centre National de
More information