ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compre...ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.展开更多
The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux...The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.展开更多
The paper is devised to combine the approximated semi-Lagrange weighted essentially non-oscillatory scheme and flux vector splitting. The approximated finite volume semi-Lagrange that is weighted essentially non-oscil...The paper is devised to combine the approximated semi-Lagrange weighted essentially non-oscillatory scheme and flux vector splitting. The approximated finite volume semi-Lagrange that is weighted essentially non-oscillatory scheme with Roe flux had been proposed. The methods using Roe speed to construct the flux probably generates entropy-violating solutions. More seriously, the methods maybe perform numerical instability in two-dimensional cases. A robust and simply remedy is to use a global flux splitting to substitute Roe flux. The combination is tested by several numerical examples. In addition, the comparisons of computing time and resolution between the classical weighted essentially non-oscillatory scheme (WENOJS-LF) and the semi-Lagrange weighted essentially non-oscillatory scheme (WENOEL-LF) which is presented (both combining with the flux vector splitting).展开更多
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 presents a new Lagrangian type scheme for solving the Euler equations of compressible gas dynamics.In this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method...This paper presents a new Lagrangian type scheme for solving the Euler equations of compressible gas dynamics.In this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method,and the mesh moves with the fluid flow.The scheme is conservative for the mass,momentum and total energy and maintains second-order accuracy.The scheme avoids solving the geometrical part and has free parameters.Results of some numerical tests are presented to demonstrate the accuracy and the non-oscillatory property of the scheme.展开更多
The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is base...The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.展开更多
A new approach for treating the mesh with Lagrangian scheme of finite volume method is presented. It has been proved that classical Lagrangian method is difficult to cope with large deformation in tracking material pa...A new approach for treating the mesh with Lagrangian scheme of finite volume method is presented. It has been proved that classical Lagrangian method is difficult to cope with large deformation in tracking material particles due to severe distortion of cells, and the changing connectivity of the mesh seems especially attractive for solving such issues. The mesh with large deformation based on computational geometry is optimized by using new method. This paper develops a processing system for arbitrary polygonal unstructured grid,the intelligent variable grid neighborhood technologies is utilized to improve the quality of mesh in calculation process, and arbitrary polygonal mesh is used in the Lagrangian finite volume scheme. The performance of the new method is demonstrated through series of numerical examples, and the simulation capability is efficiently presented in coping with the systems with large deformations.展开更多
基金jointly sponsored by the Key Project of the Chinese National Programs for Fundamental Research and Development ("973 Program" Grant No.2013CB430106)+1 种基金the Key Project of the Chinese National Science & Technology Pillar Program during the Twelfth Five-year Plan Period (Grant No.2012BAC22B01)the National Natural Science Foundation of China ( Grant No.41375108)
文摘ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.
文摘The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.
文摘The paper is devised to combine the approximated semi-Lagrange weighted essentially non-oscillatory scheme and flux vector splitting. The approximated finite volume semi-Lagrange that is weighted essentially non-oscillatory scheme with Roe flux had been proposed. The methods using Roe speed to construct the flux probably generates entropy-violating solutions. More seriously, the methods maybe perform numerical instability in two-dimensional cases. A robust and simply remedy is to use a global flux splitting to substitute Roe flux. The combination is tested by several numerical examples. In addition, the comparisons of computing time and resolution between the classical weighted essentially non-oscillatory scheme (WENOJS-LF) and the semi-Lagrange weighted essentially non-oscillatory scheme (WENOEL-LF) which is presented (both combining with the flux vector splitting).
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.11171038)the Science Foundation of China Academy of Engineering Physics,China(Grant No.2013A0202011)the National Council for Scientific and Technological Development of Brazil(CNPq).
文摘This paper presents a new Lagrangian type scheme for solving the Euler equations of compressible gas dynamics.In this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method,and the mesh moves with the fluid flow.The scheme is conservative for the mass,momentum and total energy and maintains second-order accuracy.The scheme avoids solving the geometrical part and has free parameters.Results of some numerical tests are presented to demonstrate the accuracy and the non-oscillatory property of the scheme.
文摘The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.
基金supported in part by the National Natural Science Foundation of China under Grant 11372051,Grant 11475029part by the Fund of the China Academy of Engineering Physics under Grant 20150202045
文摘A new approach for treating the mesh with Lagrangian scheme of finite volume method is presented. It has been proved that classical Lagrangian method is difficult to cope with large deformation in tracking material particles due to severe distortion of cells, and the changing connectivity of the mesh seems especially attractive for solving such issues. The mesh with large deformation based on computational geometry is optimized by using new method. This paper develops a processing system for arbitrary polygonal unstructured grid,the intelligent variable grid neighborhood technologies is utilized to improve the quality of mesh in calculation process, and arbitrary polygonal mesh is used in the Lagrangian finite volume scheme. The performance of the new method is demonstrated through series of numerical examples, and the simulation capability is efficiently presented in coping with the systems with large deformations.
