Partially implicit peer methods for the compressible Euler equations

When cut cells are used for the representation of orography in numerical weather prediction models this leads to very small cells. On one hand this results in very harsh time step restrictions for explicit methods due to the CFL criterion. On the other hand cut cells only appear in a small region of the domain. Therefore we consider a partially implicit method: In cut cells the Jacobian incorporates advection, diffusion and acoustics while in the full cells of the free atmosphere only the acoustic part is used, i.e. the method is linearly implicit in the cut cell regions and semi-explicit in the free regions and computes with time step sizes restricted only by the CFL condition in the free atmosphere. Furthermore we use a simplified Jacobian in the cut cell regions in order to save storage and gain computational efficiency. While the method retains the order independently of the Jacobian we present a linear stability theory which takes the effects of the simplifications of the Jacobian on stability into account. The presented method is as stable and accurate as the underlying split-explicit method but furthermore it can compute with cut cells with nearly no additional effort.

[1]  Jordan G. Powers,et al.  A Description of the Advanced Research WRF Version 2 , 2005 .

[2]  Rüdiger Weiner,et al.  Explicit Two-Step Peer Methods for the Compressible Euler Equations , 2009 .

[3]  Roni Avissar,et al.  The Ocean-Land-Atmosphere Model (OLAM). Part II: Formulation and Tests of the Nonhydrostatic Dynamic Core , 2008 .

[4]  George H. Bryan,et al.  A Benchmark Simulation for Moist Nonhydrostatic Numerical Models , 2002 .

[5]  K. Burrage,et al.  Non-linear stability of a general class of differential equation methods , 1980 .

[6]  Helmut Podhaisky,et al.  Explicit two-step peer methods , 2008, Comput. Math. Appl..

[7]  Louis J. Wicker,et al.  Time-Splitting Methods for Elastic Models Using Forward Time Schemes , 2002 .

[8]  H. Munthe-Kaas,et al.  Mathematics and Computation, a Contemporary View , 2008 .

[9]  G. Powers,et al.  A Description of the Advanced Research WRF Version 3 , 2008 .

[10]  Rüdiger Weiner,et al.  Parameter optimization for explicit parallel peer two-step methods , 2009 .

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

[12]  Joseph B. Klemp,et al.  Behavior of flow over step orography , 2000 .

[13]  Helmut Podhaisky,et al.  Two-step W-methods and their application to MOL-systems , 2003 .

[14]  Eitan Tadmor,et al.  Energy-Preserving and Stable Approximations for the Two-Dimensional Shallow Water Equations , 2008 .

[15]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[16]  R. Pielke,et al.  A comprehensive meteorological modeling system—RAMS , 1992 .

[17]  John C. Butcher,et al.  General linear methods for ordinary differential equations , 2009, Math. Comput. Simul..

[18]  Oswald Knoth,et al.  Non-dissipative cloud transport in Eulerian grid models by the volume-of-fluid (VOF) method , 2005 .

[19]  R. Rannacher,et al.  Benchmark Computations of Laminar Flow Around a Cylinder , 1996 .

[20]  Helmut Podhaisky,et al.  High-order linearly implicit two-step peer -- finite element methods for time-dependent PDEs , 2009 .

[21]  Helmut Podhaisky,et al.  Rosenbrock-type 'Peer' two-step methods , 2005 .

[22]  Rüdiger Weiner,et al.  Linearly implicit peer methods for the compressible Euler equations , 2012 .

[23]  Gerd Steinebach,et al.  Order-reduction of ROW-methods for DAEs and method of lines applications , 1995 .

[24]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[25]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

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

[27]  Jürgen Steppeler,et al.  Nonhydrostatic Atmospheric Modeling using az-Coordinate Representation , 2002 .

[28]  William C. Skamarock,et al.  A time-split nonhydrostatic atmospheric model for weather research and forecasting applications , 2008, J. Comput. Phys..

[29]  Jens Lang,et al.  ROS3P—An Accurate Third-Order Rosenbrock Solver Designed for Parabolic Problems , 2000 .

[30]  Helmut Podhaisky,et al.  Superconvergent explicit two-step peer methods , 2009 .

[31]  Joachim Rang,et al.  New Rosenbrock W-Methods of Order 3 for Partial Differential Algebraic Equations of Index 1 , 2005 .

[32]  N. Nikiforakis,et al.  Well-balanced compressible cut-cell simulation of atmospheric flow , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[33]  Lutz Angermann,et al.  New Rosenbrock methods of order 3 for PDAEs of index 2 , 2007 .

[34]  Hilding Sundqvist,et al.  On vertical interpolation and truncation in connexion with use of sigma system models , 1976 .