Generalized moving least squares and moving kriging least squares approximations for solving the transport equation on the sphere
This work provides meshless numerical methods for solving transport equations on the sphere, which is relevant for geophysical and atmospheric modeling.
The authors apply two meshless methods (generalized moving least squares and moving kriging least squares) to solve the time-dependent transport equation on the sphere, achieving accurate solutions for solid body rotation, vortex roll-up, and deformational flow test problems.
In this work, we apply two meshless methods for the numerical solution of the time-dependent transport equation defined on the sphere in spherical coordinates. The first technique, which was introduced by Mirzaei (BIT Numerical Mathematics, 54 (4) 1041-1063, 2017) in Cartesian coordinates is a generalized moving least squares approximation, and the second one, which is developed here, is moving kriging least squares interpolation on the sphere. These methods do not depend on the background mesh or triangulation, and they can be implemented on the transport equation in spherical coordinates easily using different distribution points. Furthermore, the time variable is approximated by a second-order backward differential formula. The obtained fully discrete scheme is solved via the biconjugate gradient stabilized algorithm with zero-fill incomplete lower-upper (ILU) preconditioner at each time step. Three well-known test problems namely solid body rotation, vortex roll-up, and deformational flow are solved to demonstrate our developments.