An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

• Published in: Advances in Atmospheric Sciences Vol. 39; no. 3; pp. 487 - 501
• Main Authors: Huang, Pei, Chen, Chungang, Li, Xingliang, Shen, Xueshun, Xiao, Feng
• Format: Journal Article
• Language: English
• Published: Heidelberg Springer Science and Business Media LLC 01.03.2022; Science Press; Springer Nature B.V; State Key Laboratory for Strength and Vibration of Mechanical Structures,Xi'an Jiaotong University,Xi'an 710049,China%Center of Numerical Weather Prediction,China Meteorological Administration,Beijing 100081,China%Department of Mechanical Engineering,Tokyo Institute of Technology,Tokyo 226-8502,Japan
• Subjects: Algorithms; Atmospheric dynamics; Atmospheric Sciences; Computer applications; Conservation; Constraints; Earth and Environmental Science; Earth Sciences; Finite element method; Finite volume method; Geophysics/Geodesy; Grid refinement (mathematics); Interpolation; Mathematical models; Meteorology; Modelling; Numerical experiments; Original Paper; 多矩限制有限体积法; nonhydrostatic model; 自适应网格加密; 动力框架; dynamical core; 非静力模式; multi-moment constrained finite-volume method; adaptive mesh refinement; high-order methods; 高精度格式
• ISSN: 0256-1530; 1861-9533
• DOI: 10.1007/s00376-021-1185-9

An adaptive 2D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume (MCV) scheme and the Berger-Oliger adaptive mesh refinement (AMR) algorithm. The MCV scheme takes several point-wise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average (VIA) and point value (PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrain-following grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarse-fine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.