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.展开更多
The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechan...The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics (English Edition), 2007, 28(7), 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.展开更多
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.展开更多
In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(...In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.展开更多
This paper presents an explicit difference scheme with accuracy and branching stability for solving onedimensional parabolic type equation by the method of undetermined parameters and its truncation error is O(△t4+△...This paper presents an explicit difference scheme with accuracy and branching stability for solving onedimensional parabolic type equation by the method of undetermined parameters and its truncation error is O(△t4+△x4). The stability condition is r=a△t/△x2<1/2.展开更多
The heat equation is a second-order parabolic partial differential equation, which can be solved in many ways using numerical methods. This paper provides a numerical solution that uses the finite difference method li...The heat equation is a second-order parabolic partial differential equation, which can be solved in many ways using numerical methods. This paper provides a numerical solution that uses the finite difference method like the explicit center difference method. The forward time and centered space (FTCS) is used to a problem containing the one-dimensional heat equation and the stability condition of the scheme is reported with different thermal conductivity of different materials. In this study, results obtained for different thermal conductivity of distinct materials are compared. Also, the results reveal the well-behavior properties of the materials in good agreement.展开更多
In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite differ...In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.展开更多
A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordin...A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.展开更多
The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the ...The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.展开更多
This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation (2.1) is set up by means of introducing incr...This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation (2.1) is set up by means of introducing incremental unknowns methods. Through the stability analyzing for the scheme, it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.展开更多
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.展开更多
A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorth...A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided.The Delaunay triangulation method is improved to generate the unstructured grid.The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow.The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively.In the former case,a good agreement is achieved in the main velocity and drag coefficient curve.In the latter case,the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment.It is indicated that the factor difference scheme is of higher accuracy,and feasible to be applied to Navier-Stokes equation.展开更多
In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical...In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical observation data.The schemes have three advantages:high-order accuracy in time,generalized square conservation,and smart use of historical observations.Numerical tests based on the one-dimensional linear advection equations suggest that reasonable consideration of accuracy,square conservation,and inclusion of historical observations is critical for good performance of a finite difference scheme.展开更多
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.展开更多
Much effort has been devoted to researching the common Rosenau equation, but the numerical method of it has not been studied. In this paper, a conservative Crank-Nicolson difference scheme for an initial-boundary valu...Much effort has been devoted to researching the common Rosenau equation, but the numerical method of it has not been studied. In this paper, a conservative Crank-Nicolson difference scheme for an initial-boundary value problem of the generalized Rosenau equation is proposed. Existence and uniqueness of numerical solutions are derived. By method of discrete energy, the second order convergence and stability are discussed. Numerical examples demonstrate the theoretical results.展开更多
The study by H.X.Zhang shows that in order to suppress the spurious oscillation at both upstream and downstream of the shock, the coefficient of the third_order derivative on the right hand side of the modified equati...The study by H.X.Zhang shows that in order to suppress the spurious oscillation at both upstream and downstream of the shock, the coefficient of the third_order derivative on the right hand side of the modified equation of the difference scheme must be positive upstream and negative downstream of the shock. According to this principle, a new non_oscillatory, containing no free parameters and dissipative difference scheme of second_order both in time and space is proposed. It is proved that this scheme possesses TVD property and is generalized Gudunov scheme of second_order. In the presence of the shock wave in the flow field, this scheme is the generalization and improvement of the Lax_Wendroff scheme. Several numerical examples are given which demonstrate that the proposed scheme is non_oscillatory of high order accuracy and high resolution. It also has the advantages of compact form, greater maximum allowable Courant number and convenient to use.展开更多
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.展开更多
文摘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.
基金Supported by the National Natural Science Foundation of China (Nos.50876114 and 10602043)the Program for New Century Excellent Talents in University,and the Scientific Research Key Project Fund of Ministry of Education (No.106142)
文摘The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics (English Edition), 2007, 28(7), 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.
文摘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.
文摘In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.
文摘This paper presents an explicit difference scheme with accuracy and branching stability for solving onedimensional parabolic type equation by the method of undetermined parameters and its truncation error is O(△t4+△x4). The stability condition is r=a△t/△x2<1/2.
文摘The heat equation is a second-order parabolic partial differential equation, which can be solved in many ways using numerical methods. This paper provides a numerical solution that uses the finite difference method like the explicit center difference method. The forward time and centered space (FTCS) is used to a problem containing the one-dimensional heat equation and the stability condition of the scheme is reported with different thermal conductivity of different materials. In this study, results obtained for different thermal conductivity of distinct materials are compared. Also, the results reveal the well-behavior properties of the materials in good agreement.
文摘In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.
基金The National Natural Science Foundation of China (No10471023)
文摘A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.
文摘The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.
文摘This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation (2.1) is set up by means of introducing incremental unknowns methods. Through the stability analyzing for the scheme, it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.
基金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.
基金Supported by National Natural Science Foundation of China (No. 10632050)
文摘A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided.The Delaunay triangulation method is improved to generate the unstructured grid.The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow.The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively.In the former case,a good agreement is achieved in the main velocity and drag coefficient curve.In the latter case,the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment.It is indicated that the factor difference scheme is of higher accuracy,and feasible to be applied to Navier-Stokes equation.
基金the Ministry of Science and Technology of China for funding the National Basic Research Program of China (973 Program,Grant No.2011CB309704)
文摘In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical observation data.The schemes have three advantages:high-order accuracy in time,generalized square conservation,and smart use of historical observations.Numerical tests based on the one-dimensional linear advection equations suggest that reasonable consideration of accuracy,square conservation,and inclusion of historical observations is critical for good performance of a finite difference scheme.
基金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.
基金The National Natural Science Foundation of China (No.40701014)the Scientific Research Fund of Sichuan Provincial Education Department (No.09ZB081)the Research Fund of key Discipline of Xihua University:Applied Mathe-matics (No.XZD0910-09-1)
文摘Much effort has been devoted to researching the common Rosenau equation, but the numerical method of it has not been studied. In this paper, a conservative Crank-Nicolson difference scheme for an initial-boundary value problem of the generalized Rosenau equation is proposed. Existence and uniqueness of numerical solutions are derived. By method of discrete energy, the second order convergence and stability are discussed. Numerical examples demonstrate the theoretical results.
文摘The study by H.X.Zhang shows that in order to suppress the spurious oscillation at both upstream and downstream of the shock, the coefficient of the third_order derivative on the right hand side of the modified equation of the difference scheme must be positive upstream and negative downstream of the shock. According to this principle, a new non_oscillatory, containing no free parameters and dissipative difference scheme of second_order both in time and space is proposed. It is proved that this scheme possesses TVD property and is generalized Gudunov scheme of second_order. In the presence of the shock wave in the flow field, this scheme is the generalization and improvement of the Lax_Wendroff scheme. Several numerical examples are given which demonstrate that the proposed scheme is non_oscillatory of high order accuracy and high resolution. It also has the advantages of compact form, greater maximum allowable Courant number and convenient to use.
基金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.