A semi‐implicit non‐hydrostatic dynamical kernel using finite elements in the vertical discretization

This work is a first step in the direction of implementing a high-order finite-element discretization in the vertical in the non-hydrostatic version of the HARMONIE model. The present dynamical core of the HARMONIE model is shared with the ECMWF and the ALADIN models and uses a horizontal spectral discretization and a semi-implicit semi-Lagrangian time stepping scheme, all of which are maintained in this work. Trying to implement a finite-element discretization in the non-hydrostatic version of the HARMONIE model has been found to be very difficult due to the set of prognostic variables used and the mass-based vertical coordinate. A different set of prognostic variables and a hybrid vertical coordinate based on height are tested here on a vertical slice non-hydrostatic kernel. A stability analysis of the linear model has been done. To evaluate the model stability and accuracy, a set of test cases from the literature are presented in the linear and nonlinear regimes, with and without orography. An iterative centred-implicit scheme can be applied to avoid instability related to steep orography, although this reduces the efficiency of the model. The novel aspects with respect to existing non-hydrostatic model kernels are the use of cubic finite elements in the vertical discretization, the use of a height-based vertical coordinate in conjunction with a spectral discretization in the horizontal, and the coordinate-independent formulation of each element of the model including the semi-Lagrangian advection. Copyright © 2011 Royal Meteorological Society

[1]  Louis J. Wicker,et al.  Numerical solutions of a non‐linear density current: A benchmark solution and comparisons , 1993 .

[2]  Pierre Bénard,et al.  Dynamical kernel of the Aladin–NH spectral limited‐area model: Revised formulation and sensitivity experiments , 2010 .

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

[4]  Xiaohua Yang.  A nonhydrostatic model for simulation of airflow over mesoscale bell-shaped ridges , 1993 .

[5]  Thomas Bewley,et al.  On the contravariant form of the Navier-Stokes equations in time-dependent curvilinear coordinate systems , 2004 .

[6]  René Laprise,et al.  The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable , 1992 .

[7]  Oliver Fuhrer,et al.  Numerical consistency of metric terms in terrain-following coordinates , 2003 .

[8]  A. Robert.  A Semi-Lagrangian and Semi-Implicit Numerical Integration Scheme for the Primitive Meteorological Equations , 1982 .

[9]  Jean Côté,et al.  Stability of Vertical Discretization Schemes for Semi-Implicit Primitive Equation Models: Theory and Application , 1983 .

[10]  Roger A. Pielke,et al.  The Derivation of a Terrain-Following Coordinate System for Use in a Hydrostatic Model , 1981 .

[11]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[12]  Robert Sharman,et al.  Incompressible and anelastic flow simulations on numerically generated grids , 1988 .

[13]  Monique Tanguay,et al.  A Semi-implicit Send-Lagrangian Fully Compressible Regional Forecast Model , 1990 .

[14]  Pierre Bénard,et al.  Integration of the fully elastic equations cast in the hydrostatic pressure terrain-following coordinate in the framework of the ARPEGE/Aladin NWP system , 1995 .

[15]  Richard Asselin,et al.  Frequency Filter for Time Integrations , 1972 .

[16]  John Thuburn,et al.  Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes , 2005 .

[17]  William C. Skamarock,et al.  Efficiency and Accuracy of the Klemp-Wilhelmson Time-Splitting Technique , 1994 .

[18]  René Laprise,et al.  A Semi-Implicit Semi-Lagrangian Regional Climate Model: The Canadian RCM , 1999 .

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

[20]  A. Simmons,et al.  Implementation of the Semi-Lagrangian Method in a High-Resolution Version of the ECMWF Forecast Model , 1995 .

[21]  Slobodan Nickovic,et al.  An Alternative Approach to Nonhydrostatic Modeling , 2001 .

[22]  J. Soulis,et al.  Computer‐aided design of channel expansions and contractions , 1992 .

[23]  Pierre Bénard,et al.  Stability of Leapfrog Constant-Coefficients Semi-Implicit Schemes for the Fully Elastic System of Euler Equations: Case with Orography , 2005 .

[24]  Pierre Benard,et al.  Stability of Semi-Implicit and Iterative Centered-Implicit Time Discretizations for Various Equation Systems Used in NWP , 2003, physics/0304114.

[25]  Jozef Vivoda,et al.  Stability of Leapfrog Constant-Coefficients Semi-Implicit Schemes for the Fully Elastic System of Euler Equations: Flat-Terrain Case , 2004 .

[26]  Mariano Hortal,et al.  A finite‐element scheme for the vertical discretization of the semi‐Lagrangian version of the ECMWF forecast model , 2004 .

[27]  D. Lüthi,et al.  A new terrain-following vertical coordinate formulation for atmospheric prediction models , 2002 .