The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-ori...The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components （called Yin and Yang respectively） that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid.展开更多
A new non-interpolating semi-Lagrangian scheme has been proposed, which can eliminate any interpolation,and consequently numerical smoothing of forecast fields. Here the new scheme is applied to KdV equation and its p...A new non-interpolating semi-Lagrangian scheme has been proposed, which can eliminate any interpolation,and consequently numerical smoothing of forecast fields. Here the new scheme is applied to KdV equation and its performance is assessed by comparing the numerical results with those produced by Ritchie＇s scheme (1986).The comparison shows that the non-interpolating semi-Lagrangian scheme appears to have efficiency advantages.展开更多
Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directl...Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile （CIP） method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.展开更多
基金This paper is sponsored by the National Natural Science Foundation of China （No. 40575050） the National Key Program for Developing Basic Reseach （＂973＂） （No. 2004CB418306）.
文摘The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components （called Yin and Yang respectively） that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid.
文摘A new non-interpolating semi-Lagrangian scheme has been proposed, which can eliminate any interpolation,and consequently numerical smoothing of forecast fields. Here the new scheme is applied to KdV equation and its performance is assessed by comparing the numerical results with those produced by Ritchie＇s scheme (1986).The comparison shows that the non-interpolating semi-Lagrangian scheme appears to have efficiency advantages.
基金supported by National Natural Science Foundation of China (NSFC) projects (Grant Nos. 40875065 and 40805045)the research projects 2008R001 at Chinese Academy of Meteorological Sciences (CAMS) and 2008 LASWZI05 at the State Key Laboratory of Severe Weather, CAMS
文摘Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile （CIP） method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.
