Abstract
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.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.