AbstractThis paper investigates the use of the weighted essentially nonoscillatory (WENO) space discretization methods of third and fifth order for momentum transport in the Meso-NH meteorological model, and their association with explicit Runge–Kutta (ERK) methods, with the specific purpose of finding an optimal combination in terms of wall-clock time to solution. A linear stability analysis using von Neumann theory is first conducted that considers six different ERK time integration methods. A new graphical representation of linear stability is proposed, which allows a first discrimination between the ERK methods. The theoretical analysis is then completed by tests on numerical problems of increasing complexity (linear advection of high wind gradient, orographic waves, density current, large eddy simulation of fog, and windstorm simulation), using a fourth-order-centered scheme as a reference basis. The five-stage third-order and fourth-order ERK combinations appear as the time integration methods of choice for coupling with WENO schemes in terms of stability. An explicit time-splitting method added to the ERK temporal scheme for WENO improves the stability properties slightly more. When the spatial discretizations are compared, WENO schemes present the main advantage of maintaining stable, nonoscillatory transitions with sharp discontinuities, but WENO third order is excessively damping, while WENO fifth order provides better accuracy. Finally, WENO fifth order combined with the ERK method makes the whole physics of the model 3 times faster compared to the classical fourth-order centered scheme associated with the leapfrog temporal scheme.
Read full abstract