Abstract:In this paper,a new numerical scheme was proposed to solve the advection equation in a multi-mo-ment nonhydrostatic dynamical core.To guarantee the shape-preserving property,the limiting operations are de-vised for a hybrid discretization framework adopted by the multi-moment dynamical core,consisting of the multi-moment finite-volume and the conservative finite-difference schemes for the horizontal and vertical discret-izations respectively.In the horizontal direction,a nonoscillaory scheme is accomplished by adjusting the slope of the multi-moment reconstruction polynomial at the cell center with the application of a WENO(weighted es-sentially non-oscillatory)algorithm.The resulting multi-moment scheme can achieve the fourth-order accuracy in the convergence test.In the vertical direction,a TVD(total variation diminishing)slope limiter is applied in the finite-difference discretization to remove the non-physical oscillations around the discontinuities.To accom-plish the time marching in the proposed advection model,a second-order paired explicit Runge-Kutta scheme is adopted,which is expected to be an efficient and practical method for the advection solvers in the atmospheric models with very high spatial resolutions.The explicit time marching,without the dimension splitting,is useful to avoid the divergence errors in the advection transport calculations.Two Runge-Kutta schemes,requiring dif-ferent times of conducting the spatial discretization within a time step,are combined,and used for the time marching in the different directions.The finite-difference discretization is called for six times within a time step in order to increase the maximum available CFL(Courant-Friedrichs-Lewy)number in the vertical direction,while the horizontal multi-moment spatial discretization is conducted for two times as the regular second-order schemes.As a result,the difference between the maximum time steps determined by the horizontal and vertical discretizations,due to the very large aspect ratio of the computational cells in atmospheric modeling,can be di-minished.The non-negativity property of the proposed advection scheme is assured by devising a new flux-cor-rection algorithm.It improves the existing positivity-preserving algorithm through further considering the mass flowing into the computational cell in an iterative procedure during the flux-correction operations.The proposed flux-correction algorithm can approach the necessary and sufficient condition for assuring the non-negative solu-tions and is more accurate for the advection calculations with CFL numbers larger than one.The widely used two-dimensional benchmark tests were checked in this study and the numerical results verified the performance the proposed advection scheme,which has the practical potential to build an accurate and efficient advection equa-tion solver for the scalable high-resolution nonhydrostatic atmospheric models.

A Performance Study of Horizontally Explicit Vertically Implicit (HEVI) Time-Integrators for Non-Hydrostatic Atmospheric Models

An energetically balanced, quasi-Newton integrator for non-hydrostatic vertical atmospheric dynamics

A Comparison of Rosenbrock-Wanner and Crank-Nicolson Time Integrators for Atmospheric Modelling

Applicability Study of Euler–Lagrange Integration Scheme in Constructing Small-Scale Atmospheric Dynamics Models

A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics

Parallel-in-time integration of the shallow water equations on the rotating sphere using Parareal and MGRIT

A comparative computational study of different formulations of the compressible Euler equations for mesoscale atmospheric flows in a finite volume framework

A Global Non-Hydrostatic Atmospheric Model with a Mass and Energy Conserving Vertically-Implicit-Correction (VIC) Scheme

High-Order Dynamic Integration Method (HODIM) for Modeling Turbulent Fluid Dynamics

Comparison of the performance of high-order schemes based on the gas-kinetic and HLLC fluxes

Finite-element schemes for extended integrations of atmospheric models

Exponential time integration for 3D compressible atmospheric models

Filter stabilization for the mildly compressible Euler equations with application to atmosphere dynamics simulations

An Advection Scheme Using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling

Current and Emerging Time-Integration Strategies in Global Numerical Weather and Climate Prediction

An Atmospheric Model with Semi-Lagrangian and Semi-Inplicit Eulerian Equations Preliminary Numeical Tests on Deep Convection

Adaptively Implicit Advection for Atmospheric Flows

On the performances of standard and kinetic energy preserving time-integration methods for incompressible-flow simulations

Implicit high order discontinuous Galerkin method for compressible laminar and turbulent flow simulation

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport