摘要
A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. High- precision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.
A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. High- precision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.
基金
supported by the Shandong Post-Doctoral Innovation Fund(Grant No.201303064)
the Qingdao Post-Doctoral Application Research Project
the National Basic Research(973) Program of China(Grant No.2012CB417402 and 2010CB950402)
the National Natural Science Foundation of China(Grant No.41176017)