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.展开更多
A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with...A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.展开更多
A partial Runge-Kutta Discontinuous Galerkin(RKDG)method which preserves the exactly divergence-free property of the magnetic field is proposed in this paper to solve the two-dimensional ideal compressible magnetohydr...A partial Runge-Kutta Discontinuous Galerkin(RKDG)method which preserves the exactly divergence-free property of the magnetic field is proposed in this paper to solve the two-dimensional ideal compressible magnetohydrodynamics(MHD)equations written in semi-Lagrangian formulation on moving quadrilateral meshes.In this method,the fluid part of the ideal MHD equations along with zcomponent of the magnetic induction equation is discretized by the RKDG method as our previous paper[47].The numerical magnetic field in the remaining two directions(i.e.,x and y)are constructed by using the magnetic flux-freezing principle which is the integral form of the magnetic induction equation of the ideal MHD.Since the divergence of the magnetic field in 2D is independent of its z-direction component,an exactly divergence-free numerical magnetic field can be obtained by this treatment.We propose a new nodal solver to improve the calculation accuracy of velocities of the moving meshes.A limiter is presented for the numerical solution of the fluid part of the MHD equations and it can avoid calculating the complex eigen-system of the MHD equations.Some numerical examples are presented to demonstrate the accuracy,non-oscillatory property and preservation of the exactly divergence-free property of our method.展开更多
Shape-from-shading(SFS) is one of the important approaches of 3-D surface reconstruction in computer vision. Since reflectance map equation in SFS is a nonlinear partial differential equation(PDE) with two unknown var...Shape-from-shading(SFS) is one of the important approaches of 3-D surface reconstruction in computer vision. Since reflectance map equation in SFS is a nonlinear partial differential equation(PDE) with two unknown variables, SFS with one image is ill-posed in mathematical sense. A linear perspective SFS method with two images is proposed to deal with the problem. We assume that two images with different light source directions are captured firstly. Orthogonal projection is not as accurate as perspective one to simulate imaging processes. Two reflectance map equations are established based on the Lambertian model under perspective projection, and the equations are further transformed into one linear PDE. Then the iterative semi-Lagrangian algorithm is used to approximate the solution. Finally, 3-D height values of pixel points in imaging planes are solved by the numerical interpolation method. Experimental results of both hemisphere and complex surfaces show that the proposed method can reconstruct surfaces accurately.展开更多
A 1D1 V hybrid Vlasov-fluid model was developed for this study to elucidate discharge current oscillations of Hall thrusters(HTs).The Vlasov equation for ions velocity distribution function with ionization source term...A 1D1 V hybrid Vlasov-fluid model was developed for this study to elucidate discharge current oscillations of Hall thrusters(HTs).The Vlasov equation for ions velocity distribution function with ionization source term is solved using a constrained interpolation profile conservative semiLagrangian method.The fourth-order weighted essentially non-oscillatory(4 th WENO)limiter is applied to the first derivative value to minimize numerical oscillation in the discharge oscillation analyses.The fourth-order accuracy is verified through a 1 D scalar test case.Nonoscillatory and high-resolution features of the Vlasov model are confirmed by simulating the test cases of the Vlasov–Poisson system and by comparing the results with a particle-in-cell(PIC)method.A1 D1 V HT simulation is performed through the hybrid Vlasov model.The ionization oscillation is analyzed.The oscillation amplitude and plasma density are compared with those obtained from a hybrid PIC method.The comparison indicates that the hybrid Vlasov-fluid model yields noiseless results and that the steady-state waveform is calculable in a short time period.展开更多
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.展开更多
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).展开更多
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.展开更多
In order to improve the practicality of spectral method and the efficiency of computation, the multi-spectrum method is proposed on the basis of multi-grid method. Coarse spectra are used to compute the slow nonlinear...In order to improve the practicality of spectral method and the efficiency of computation, the multi-spectrum method is proposed on the basis of multi-grid method. Coarse spectra are used to compute the slow nonlinear part (including physical process), while fine spectra are used to compute the fast linear part. This method not only can reduce computation time, but also can obtain computational efficiency similar to that from only fine spectra. Thus, it is an economical numerical method. Both explicit complete-square-conservation scheme and multispectrum scheme are used to improve IAP L 9 T 42 spectral climate models, with and without physical forcings respectively, and the advantage of reducing computation time is obtained satisfactorily. In order to overcome the difficulty that vapor equation is very sensitive to the change of time step, the square-conservation semi-Lagrangian scheme is used to solve vapor equation. Because the semi-Lagrangian scheme has the property of square-conservation, computational instability can be avoided. When time step becomes longer with the semi-Lagrangian Scheme, through numerical examples, the vapor transportation can be depicted objectively and the effect of precipitation simulation can be modified.展开更多
基金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.
基金supported by the Chinese Academy of Science Strategic Priority Research Program (Grant No. XDA05110303)"973" Program (Grant Nos. 2010CB950403,2012CB417203,and 2013CB955803)+1 种基金"863" Program(Grant No. 2010AA012305)the National Natural Science Foundation of China (Grant Nos. 40925015,40875034,and 41023002)
文摘A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.
基金supported by National Natural Science Foundation of China(Nos.12071046,11671049,91330107,11571002 and 11702028)China Postdoctoral Science Foundation(No.2020TQ0013).
文摘A partial Runge-Kutta Discontinuous Galerkin(RKDG)method which preserves the exactly divergence-free property of the magnetic field is proposed in this paper to solve the two-dimensional ideal compressible magnetohydrodynamics(MHD)equations written in semi-Lagrangian formulation on moving quadrilateral meshes.In this method,the fluid part of the ideal MHD equations along with zcomponent of the magnetic induction equation is discretized by the RKDG method as our previous paper[47].The numerical magnetic field in the remaining two directions(i.e.,x and y)are constructed by using the magnetic flux-freezing principle which is the integral form of the magnetic induction equation of the ideal MHD.Since the divergence of the magnetic field in 2D is independent of its z-direction component,an exactly divergence-free numerical magnetic field can be obtained by this treatment.We propose a new nodal solver to improve the calculation accuracy of velocities of the moving meshes.A limiter is presented for the numerical solution of the fluid part of the MHD equations and it can avoid calculating the complex eigen-system of the MHD equations.Some numerical examples are presented to demonstrate the accuracy,non-oscillatory property and preservation of the exactly divergence-free property of our method.
基金supported by the National Natural Science Foundation of China(61005015)the Third National Post-Doctoral Special Foundation of China(201003280)
文摘Shape-from-shading(SFS) is one of the important approaches of 3-D surface reconstruction in computer vision. Since reflectance map equation in SFS is a nonlinear partial differential equation(PDE) with two unknown variables, SFS with one image is ill-posed in mathematical sense. A linear perspective SFS method with two images is proposed to deal with the problem. We assume that two images with different light source directions are captured firstly. Orthogonal projection is not as accurate as perspective one to simulate imaging processes. Two reflectance map equations are established based on the Lambertian model under perspective projection, and the equations are further transformed into one linear PDE. Then the iterative semi-Lagrangian algorithm is used to approximate the solution. Finally, 3-D height values of pixel points in imaging planes are solved by the numerical interpolation method. Experimental results of both hemisphere and complex surfaces show that the proposed method can reconstruct surfaces accurately.
基金supported by the China Scholarship Council(No.201708050185)。
文摘A 1D1 V hybrid Vlasov-fluid model was developed for this study to elucidate discharge current oscillations of Hall thrusters(HTs).The Vlasov equation for ions velocity distribution function with ionization source term is solved using a constrained interpolation profile conservative semiLagrangian method.The fourth-order weighted essentially non-oscillatory(4 th WENO)limiter is applied to the first derivative value to minimize numerical oscillation in the discharge oscillation analyses.The fourth-order accuracy is verified through a 1 D scalar test case.Nonoscillatory and high-resolution features of the Vlasov model are confirmed by simulating the test cases of the Vlasov–Poisson system and by comparing the results with a particle-in-cell(PIC)method.A1 D1 V HT simulation is performed through the hybrid Vlasov model.The ionization oscillation is analyzed.The oscillation amplitude and plasma density are compared with those obtained from a hybrid PIC method.The comparison indicates that the hybrid Vlasov-fluid model yields noiseless results and that the steady-state waveform is calculable in a short time period.
基金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.
文摘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).
文摘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.
文摘In order to improve the practicality of spectral method and the efficiency of computation, the multi-spectrum method is proposed on the basis of multi-grid method. Coarse spectra are used to compute the slow nonlinear part (including physical process), while fine spectra are used to compute the fast linear part. This method not only can reduce computation time, but also can obtain computational efficiency similar to that from only fine spectra. Thus, it is an economical numerical method. Both explicit complete-square-conservation scheme and multispectrum scheme are used to improve IAP L 9 T 42 spectral climate models, with and without physical forcings respectively, and the advantage of reducing computation time is obtained satisfactorily. In order to overcome the difficulty that vapor equation is very sensitive to the change of time step, the square-conservation semi-Lagrangian scheme is used to solve vapor equation. Because the semi-Lagrangian scheme has the property of square-conservation, computational instability can be avoided. When time step becomes longer with the semi-Lagrangian Scheme, through numerical examples, the vapor transportation can be depicted objectively and the effect of precipitation simulation can be modified.
