HOPE: An Automatically Differentiable High-Order Non-Oscillatory Finite-Volume Shallow-Water Dynamic Core

Lilong Zhou
2024-12-03
Abstract:An automatically differentiable, high-order non-oscillatory finite volume shallow water dynamic core has been constructed on a cubed sphere grid. This dynamic core has four advantageous properties: high order accuracy, essential non-oscillation, mass conservation, and scalability. Besides, the code development is based on PyTorch, enabling the model to run seamlessly on both CPU and GPU, while naturally possessing the capability of automatic differentiation. We named the new dynamic core as High Order Prediction Environment (HOPE). The spatial reconstruction is based on the two-dimensional tensor product polynomial and the genuine two-dimensional Weighted Essentially Non-Oscillatory (WENO) scheme. A novel panel boundary approach ensures that the accuracy can reach arbitrary order. The algorithm has very high degree of compatibility with GPU architecture, allowing the computational overhead to be mitigated through the utilization of GPU. The LMARS (Low Mach number Approximate Riemann Solver) scheme is adopted as Riemann solvers to determine fluxes on the Gaussian points on edges. Flux across the interface between each cell edge is computed using linear Gaussian quadrature, and the tendency of prognostic values is obtained by summing all fluxes around the cell. This shallow water dynamic core demonstrates outstanding performance in ideal shallow water test cases. In the steady-state geostrophic flow, the 11th order scheme reduces errors to nearly round-off error even on coarse grids (2 degrees). Furthermore, the Rossby-Haurwitz wave preserves its shape for over 90 days. To test the non-oscillation property, we designed a cylinder dam break case. The WENO approach effectively suppresses non-physical oscillation, and the genuine two-dimensional reconstruction exhibits significantly better isotropy than the dimension-by-dimension scheme.
Fluid Dynamics,Geophysics
What problem does this paper attempt to address?