A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx...A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).展开更多
A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By usi...A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.展开更多
A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and t...A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).展开更多
In this paper, the first boundary problem of quasilinear parabolic system of second order is studied by the finite difference method with intrinsic parallelism. for the problem, the stability of the difference schemes...In this paper, the first boundary problem of quasilinear parabolic system of second order is studied by the finite difference method with intrinsic parallelism. for the problem, the stability of the difference schemes with intrinsic parallelism are justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete data of the original problem, without assuming the existence of the smooth solutions for the origillal problem.展开更多
A new method was proposed for constructing total variation diminishing (TVD) upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can e...A new method was proposed for constructing total variation diminishing (TVD) upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD (CTVD) schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.展开更多
The computational stability of the explicit difference schemes of the forced dissipative nonlinear evolution equations is analyzed and the computational quasi-stability criterion of explicit difference schemes of the ...The computational stability of the explicit difference schemes of the forced dissipative nonlinear evolution equations is analyzed and the computational quasi-stability criterion of explicit difference schemes of the forced dissipative nonlinear atmospheric equations is obtained on account of the concept of computational quasi-stability, Therefore, it provides the new train of thought and theoretical basis for designing computational stable difference scheme of the forced dissipative nonlinear atmospheric equations. Key words Computational quasi-stability - Computational stability - Forced dissipative nonlinear evolution equation - Explicit difference scheme This work was supported by the National Outstanding Youth Scientist Foundation of China (Grant No. 49825109), the Key Innovation Project of Chinese Academy of Sciences (KZCX1-10-07), the National Natural Science Foundation of China (Grant Nos, 49905007 and 49975020) and the Outstanding State Key Laboratory Project (Grant No. 40023001).展开更多
A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh rat...A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.展开更多
Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmos...Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmospheric and oceanic equations is analyzed and the computationally stable explicit complete square conservative difference schemes are constructed. The theoretical analysis and numerical experiment prove that the explicit complete square conservative difference schemes are computationally stable and deserve to be disseminated.展开更多
In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O(τ^2+ h%4), and the stability c...In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O(τ^2+ h%4), and the stability condition is given. Finally, the numerical examples and numerical results are presented to show the advantage of the schemes and the correctness of theoretical analysis.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, th...Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.展开更多
The nature of infinite-dimensional Hamiltonian systems are studied for the purpose of further study on some generalized Hamiltonian systems equipped with a given Poisson bracket. From both theoretical and practical vi...The nature of infinite-dimensional Hamiltonian systems are studied for the purpose of further study on some generalized Hamiltonian systems equipped with a given Poisson bracket. From both theoretical and practical viewpoints, we summarize a general method of constructing symplectic-like difference schemes of these kinds of systems. This study provides a new algorithm for the application of the symplectic geometry method in numerical solutions of general evolution equations.展开更多
This paper proposes a new method to improve the stability condition of difference scheme of a parabolic equation. Necessary and sufficient conditions of the stability of this new method are given and proved. Some nume...This paper proposes a new method to improve the stability condition of difference scheme of a parabolic equation. Necessary and sufficient conditions of the stability of this new method are given and proved. Some numerical examples show that this method has some calculation advantages.展开更多
This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matri...This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matrix and because of the stability of the M-matrix, it is convergent by the asynchronous iterative method on multiprocessors. Then this paper gives a class of differeifce schemes for linear elliptic PDEs so that their difference matrixes are all M-matrixes and their asynchronous parallel computation are convergent.展开更多
In this paper,two fourth-order compact finite difference schemes are derived to solve the nonlinear fourth-order wave equation which can be viewed as a generalized model from the nonlinear beam equation.Differing from...In this paper,two fourth-order compact finite difference schemes are derived to solve the nonlinear fourth-order wave equation which can be viewed as a generalized model from the nonlinear beam equation.Differing from the existing compact finite difference schemes which preserve the total energy in a recursive sense,the new schemes are proved to per-fectly preserve the total energy in the discrete sense.By using the standard energy method and the cut-off function technique,the optimal error estimates of the numerical solutions are established,and the convergence rates are of O(h^(4)+τ^(2))with mesh-size h and time-step τ.In order to improve the computational efficiency,an iterative algorithm is proposed as the outer solver and the double sweep method for pentadiagonal linear algebraic equations is introduced as the inner solver to solve the nonlinear difference schemes at each time step.The convergence of the iterative algorithm is also rigorously analyzed.Several numerical results are carried out to test the error estimates and conservative properties.展开更多
In the present paper two contents are enclosed .First ,the Fourier analysis approach of the dispersion relation and group velocity effect of finite difference schemes is discussed.the defects of the approach is pointe...In the present paper two contents are enclosed .First ,the Fourier analysis approach of the dispersion relation and group velocity effect of finite difference schemes is discussed.the defects of the approach is pointed out and the correction is made;Second,a new systematic analysis method -remaider -effect analysis (abbr.REAM)is proposed by means of the modified partial differential equations (abbr MPDE)of finite difference schemes.The analysis is based on the synthetical study of the rational dispersion-and dissipation relations of finite difference schemes.And the method clearly possesses constructivity展开更多
In this paper, a singularly perturbed boundary value problem for second order self-adjoint ordinary differential equation is discussed. A class of variational difference schemes is constructed by the finite element me...In this paper, a singularly perturbed boundary value problem for second order self-adjoint ordinary differential equation is discussed. A class of variational difference schemes is constructed by the finite element method. Uniform convergence about small parameter is proved under a weaker smooth condition with respect to the coefficients of the equation. The schemes studied in refs. [1], [3], [4] and [51 belong to the cllass.展开更多
In this paper, we construct a class of difference schemes with fitted factors for a singular perturbation problem of a self-adjoint ordinary differential equation. Using a different method from [1], by analyzing the t...In this paper, we construct a class of difference schemes with fitted factors for a singular perturbation problem of a self-adjoint ordinary differential equation. Using a different method from [1], by analyzing the truncation errors of schemes, we give the sufficient conditions under which the solution of lite difference scheme converges uniformly to the solution of the differential equation. From this we propose several specific schemes under weaker conditions, and give much higher order of uniform convergence, and applying them to example, obtain the numerical results.展开更多
In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic fi...In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.展开更多
A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and im...A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.展开更多
文摘A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).
文摘A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.
文摘A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).
文摘In this paper, the first boundary problem of quasilinear parabolic system of second order is studied by the finite difference method with intrinsic parallelism. for the problem, the stability of the difference schemes with intrinsic parallelism are justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete data of the original problem, without assuming the existence of the smooth solutions for the origillal problem.
基金Project supported by the National Natural Science Foundation of China (Nos. 10172015 and 90205010)
文摘A new method was proposed for constructing total variation diminishing (TVD) upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD (CTVD) schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.
基金the National Outstanding Youth Scientist Foundation of China (GrantNo. 49825109), the Key Innovation Project of Chinese Academ
文摘The computational stability of the explicit difference schemes of the forced dissipative nonlinear evolution equations is analyzed and the computational quasi-stability criterion of explicit difference schemes of the forced dissipative nonlinear atmospheric equations is obtained on account of the concept of computational quasi-stability, Therefore, it provides the new train of thought and theoretical basis for designing computational stable difference scheme of the forced dissipative nonlinear atmospheric equations. Key words Computational quasi-stability - Computational stability - Forced dissipative nonlinear evolution equation - Explicit difference scheme This work was supported by the National Outstanding Youth Scientist Foundation of China (Grant No. 49825109), the Key Innovation Project of Chinese Academy of Sciences (KZCX1-10-07), the National Natural Science Foundation of China (Grant Nos, 49905007 and 49975020) and the Outstanding State Key Laboratory Project (Grant No. 40023001).
文摘A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.
基金the Outstanding State Key Laboratory Project of National Science Foundation of China (Grant No. 40023001 )the Key Innovatio
文摘Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmospheric and oceanic equations is analyzed and the computationally stable explicit complete square conservative difference schemes are constructed. The theoretical analysis and numerical experiment prove that the explicit complete square conservative difference schemes are computationally stable and deserve to be disseminated.
基金Supported by NSF of the Education Department of Henan Province(20031100010)
文摘In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O(τ^2+ h%4), and the stability condition is given. Finally, the numerical examples and numerical results are presented to show the advantage of the schemes and the correctness of theoretical analysis.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
基金supported by the National Natural Science Foundation of China (Grant No. U1811464)。
文摘Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.
基金Acknowledgments. This work was supported by the China National Key Development Planning Project for Ba-sic Research (Abbreviation: 973 Project Grant No. G1999032801), the Chinese Academy of Sciences Key Innovation Direction Project (Grant No. KZCX2208)
文摘The nature of infinite-dimensional Hamiltonian systems are studied for the purpose of further study on some generalized Hamiltonian systems equipped with a given Poisson bracket. From both theoretical and practical viewpoints, we summarize a general method of constructing symplectic-like difference schemes of these kinds of systems. This study provides a new algorithm for the application of the symplectic geometry method in numerical solutions of general evolution equations.
文摘This paper proposes a new method to improve the stability condition of difference scheme of a parabolic equation. Necessary and sufficient conditions of the stability of this new method are given and proved. Some numerical examples show that this method has some calculation advantages.
文摘This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matrix and because of the stability of the M-matrix, it is convergent by the asynchronous iterative method on multiprocessors. Then this paper gives a class of differeifce schemes for linear elliptic PDEs so that their difference matrixes are all M-matrixes and their asynchronous parallel computation are convergent.
基金supported by the National Natural Science Foundation of China under Grant No.11571181the Natural Science Foundation of Jiangsu Province of China under Grant No.BK20171454.
文摘In this paper,two fourth-order compact finite difference schemes are derived to solve the nonlinear fourth-order wave equation which can be viewed as a generalized model from the nonlinear beam equation.Differing from the existing compact finite difference schemes which preserve the total energy in a recursive sense,the new schemes are proved to per-fectly preserve the total energy in the discrete sense.By using the standard energy method and the cut-off function technique,the optimal error estimates of the numerical solutions are established,and the convergence rates are of O(h^(4)+τ^(2))with mesh-size h and time-step τ.In order to improve the computational efficiency,an iterative algorithm is proposed as the outer solver and the double sweep method for pentadiagonal linear algebraic equations is introduced as the inner solver to solve the nonlinear difference schemes at each time step.The convergence of the iterative algorithm is also rigorously analyzed.Several numerical results are carried out to test the error estimates and conservative properties.
文摘In the present paper two contents are enclosed .First ,the Fourier analysis approach of the dispersion relation and group velocity effect of finite difference schemes is discussed.the defects of the approach is pointed out and the correction is made;Second,a new systematic analysis method -remaider -effect analysis (abbr.REAM)is proposed by means of the modified partial differential equations (abbr MPDE)of finite difference schemes.The analysis is based on the synthetical study of the rational dispersion-and dissipation relations of finite difference schemes.And the method clearly possesses constructivity
文摘In this paper, a singularly perturbed boundary value problem for second order self-adjoint ordinary differential equation is discussed. A class of variational difference schemes is constructed by the finite element method. Uniform convergence about small parameter is proved under a weaker smooth condition with respect to the coefficients of the equation. The schemes studied in refs. [1], [3], [4] and [51 belong to the cllass.
文摘In this paper, we construct a class of difference schemes with fitted factors for a singular perturbation problem of a self-adjoint ordinary differential equation. Using a different method from [1], by analyzing the truncation errors of schemes, we give the sufficient conditions under which the solution of lite difference scheme converges uniformly to the solution of the differential equation. From this we propose several specific schemes under weaker conditions, and give much higher order of uniform convergence, and applying them to example, obtain the numerical results.
文摘In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.
文摘A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.