A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

JOURNAL OF COMPUTATIONAL PHYSICS

卷:229

Article

Ii, Satoshi²  Xiao, Feng^1,3 

[1] Tokyo Inst Technol, Dept Energy Sci, Midori Ku, Yokohama, Kanagawa 2268502, Japan

[2] Univ Tokyo, Dept Mech Engn, Bunkyo Ku, Tokyo 1138656, Japan

[3] Chinese Acad Sci, Inst Mech, LHD, Beijing 100080, Peoples R China

关键词: Global shallow water model;    Icosahedral grid;    High order accuracy;    Finite volume method;    Multi-moment;    Compact-stencil;   

DOI  :  10.1016/j.jcp.2009.11.008

来源: Elsevier

【 摘 要 】

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

【 授权许可】

Free