A novel numerical scheme is proposed to solve the shallow water equations (SWEs) on arbitrary rotating curved surfaces. Based on the method of moving frames (MMF) in which the geometry is represented by orthonormal vectors, the proposed scheme not only has the fewest dimensionality both in space and time, but also does not require either of metric tensors, composite meshes, or the ambient space. The MMF–SWE formulation is numerically discretized using the discontinuous Galerkin method of arbitrary polynomial order p in space and an explicit Runge–Kutta scheme in time. The numerical model is validated against six standard tests on the sphere and the optimal order of convergence of p+1 is numerically demonstrated. The MMF–SWE scheme is also demonstrated for its efficiency and stability on the general rotating surfaces such as ellipsoid, irregular, and non-convex surfaces.
Bibliographical noteFunding Information:
The authors would like to thank Professor Janusz Pudykiewicz (Environment Canada) for inspiring discussions and valuable comments. The research of the first author was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (No. 2016R1D1A1A02937255).
© 2016 Elsevier Inc.
All Science Journal Classification (ASJC) codes
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy(all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics