Hamiltonian Particle-Mesh Method for Two-Layer Shallow-Water Equations Subject to the Rigid-Lid Approximation
Hamiltonian Particle-Mesh Method for Two-Layer Shallow-Water Equations Subject to the Rigid-Lid Approximation
复制标题
DOI：
10.1137/030600076
复制
发表时间：
2004
期刊：
影响因子：
--
通讯作者：
C. Cotter;J. Frank;S. Reich
中科院分区：
文献类型：
--
作者：
C. Cotter;J. Frank;S. Reich
本文发展了一种求解两层浅水方程的颗粒网格方法。该方法是基于最近提出的哈密顿粒子网格（HPM）方法和解释的刚性盖近似作为一组完整的约束。建议的空间离散化导致一个受约束的Hamilton系统的常微分方程，这是集成在时间上使用的一个变体的辛SHAKE/RATTLE算法。结果表明，用刚盖近似消除外部重力波是一种计算稳定和有效的方法。
We develop a particle-mesh method for two-layer shallow-water equations subject to the rigid-lid approximation. The method is based on the recently proposed Hamiltonian particle-mesh (HPM) method and the interpretation of the rigid-lid approximation as a set of holonomic constraints. The suggested spatial discretization leads to a constrained Hamiltonian system of ODEs which is integrated in time using a variant of the symplectic SHAKE/RATTLE algorithm. It is demonstrated that the elimination of external gravity waves by the rigid-lid approximation can be achieved in a computationally stable and efficient way.