Semi-Implicit Scheme for the DWD Lokal-Modell

Summary The fourth generation of numerical weather prediction (NWP) models is currently under development at the Deutscher Wetterdienst (DWD) consisting of a global grid point model (GME) and limited-area Lokal-Modell (LM). The nonhydrostatic fully compressible LM has been designed to meet high-resol ution regional fore-cast requirements at meso-β and meso-γ scales. The initial LM implementation is based on the NCAR/Penn State MM5 with the addition of a novel generalized terrain-following coordinate and rotated lat-lon grid. A fully 3D semi-implicit time-stepping scheme has been implemented by retaining the full buoyancy term instead of the approximate form found in MM5. In contrast with earlier schemes, mass-lumping is not applied to simplify the elliptic operator on an Arakawa-C/Lorenz grid. The resulting variable-coefficient elliptic problem is solved using a minimal residual Krylov iterative method with line relaxation preconditione rs. The new semi-implicit scheme is compared with a variant of the Klemp–Wilhelmson split-explicit scheme (horizontal explicit, vertical implicit) on the basis of computational efficiency and accuracy at resolutions ranging from 7 km to 400 m. Both idealized 3D mountain wave flows and naturally occuring flows are analyzed. Below the tropopause, the 3D semi-implicit scheme can be more efficient for low Mach number M ≪ 1 flows when the number of small time steps Δts of the split-explicit approach increases with the sound-speed Courant number.

[1]  Louis J. Wicker,et al.  A Time-Splitting Scheme for the Elastic Equations Incorporating Second-Order Runge–Kutta Time Differencing , 1998 .

[2]  S. Eisenstat,et al.  Variational Iterative Methods for Nonsymmetric Systems of Linear Equations , 1983 .

[3]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[4]  J. Dudhia.  A Nonhydrostatic Version of the Penn State–NCAR Mesoscale Model: Validation Tests and Simulation of an Atlantic Cyclone and Cold Front , 1993 .

[5]  G. Doms,et al.  The Nonhydrostatic Limited-Area Model LM (Lokal-Modell) of DWD: Part I: Scientific Documentation (Ve , 1999 .

[6]  Stephen J. Thomas,et al.  A new adiabatic kernel for the MC2 model , 1998 .

[7]  Kazuo Saito,et al.  3-D Mountain Waves by the Lokal-Modell of DWD and the MRI Mesoscale Nonhydrostatic Model. , 1998 .

[8]  Piotr K. Smolarkiewicz,et al.  Preconditioned Conjugate-Residual Solvers for Helmholtz Equations in Nonhydrostatic Models , 1997 .

[9]  Motohki Ikawa,et al.  Comparison of Some Schemes for Nonhydrostatic Models with Orography , 1988 .

[10]  D. Durran,et al.  A Compressible Model for the Simulation of Moist Mountain Waves , 1983 .

[11]  T. Clark.  A small-scale dynamic model using a terrain-following coordinate transformation , 1977 .

[12]  H. Davies,et al.  A lateral boundary formulation for multi-level prediction models. [numerical weather forecasting , 1976 .

[13]  Richard C. J. Somerville,et al.  On the use of a coordinate transformation for the solution of the Navier-Stokes equations , 1975 .

[14]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[15]  M. Cullen.  A test of a semi‐implicit integration technique for a fully compressible non‐hydrostatic model , 1990 .

[16]  J. Klemp,et al.  The Simulation of Three-Dimensional Convective Storm Dynamics , 1978 .

[17]  Ronald B. Smith.  Linear theory of stratified hydrostatic flow past an isolated mountain , 1980 .

[18]  Len G. Margolin,et al.  On Forward-in-Time Differencing for Fluids: Stopping Criteria for Iterative Solutions of Anelastic Pressure Equations , 1997 .