期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 449, 期 -, 页码 -出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2021.110792
关键词
Geophysical fluid dynamics; Space-time tensors; Direct flux reconstruction; Jacobian-free methods; Exponential integrators
资金
- National Science Foundation, Computational Mathematics Program [1115978]
- Direct For Mathematical & Physical Scien [1115978] Funding Source: National Science Foundation
- Division Of Mathematical Sciences [1115978] Funding Source: National Science Foundation
A novel numerical approach for solving shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. The method uses a space-time tensor formalism and direct flux reconstruction for spatial discretization. Exponential integration is employed for larger time step sizes and improved efficiency in time integration. New multistep-type exponential propagation iterative methods allow high-order accuracy in time integration without significant increases in computational time.
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion covariantly and to describe the geometry of the rotated cubed-sphere grid. The spatial discretization is done with the direct flux reconstruction method, which is an alternative formulation to the discontinuous Galerkin approach. The equations of motion are solved in differential form and the resulting discretization is free from quadrature rules. It is well known that the time step of traditional explicit methods is limited by the phase velocity of the fastest waves. Exponential integration is employed to enable integrations with significantly larger time step sizes and improve the efficiency of the overall time integration. New multistep-type exponential propagation iterative methods of orders 4, 5 and 6 are constructed and applied to integrate the shallow-water equations in time. These new schemes enable time integration with high-order accuracy but without significant increases in computational time compared to low-order methods. The exponential matrix functions-vector products used in the exponential schemes are approximated using the complex-step approximation of the Jacobian in the Krylov-based KIOPS (Krylov with incomplete orthogonalization procedure solver) algorithm. Performance of the new numerical methods is evaluated using a set of standard benchmark tests. Crown Copyright (C) 2021 Published by Elsevier Inc.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据