Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under sat...Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.展开更多
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.展开更多
It is known that the Allen-Chan equations satisfy the maximum principle. Is this true for numerical schemes？ To the best of our knowledge, the state-of-art stability framework is the nonlinear energy stability which ...It is known that the Allen-Chan equations satisfy the maximum principle. Is this true for numerical schemes？ To the best of our knowledge, the state-of-art stability framework is the nonlinear energy stability which has been studied extensively for the phase field type equations. In this work, we will show that a stronger stability under the infinity norm can be established for the implicit-explicit discretization in time and central finite difference in space. In other words, this commonly used numerical method for the Allen-Cahn equation preserves the maximum principle.展开更多
In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with...In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with adjustable time intervals, which is proved to be effective and remarkably time-saving in numerical tests and applications.展开更多
In this article,a Crank-Nicolson/Explicit scheme is designed and analyzed for the time-dependent natural convection problem with nonsmooth initial data.The Galerkin finite element method(FEM)with stable MINI element i...In this article,a Crank-Nicolson/Explicit scheme is designed and analyzed for the time-dependent natural convection problem with nonsmooth initial data.The Galerkin finite element method(FEM)with stable MINI element is used for the velocity and pressure and linear polynomial for the temperature.The time discretization is based on the Crank-Nicolson scheme.In order to simplify the computations,the nonlinear terms are treated by the explicit scheme.The advantages of our numerical scheme can be list as follows:(1)The original problem is split into two linear subproblems,these subproblems can be solved in each time level in parallel and the computational sizes are smaller than the origin one.(2)A constant coefficient linear discrete algebraic system is obtained in each subproblem and the computation becomes easy.The main contributions of this work are the stability and convergence results of numerical solutions with nonsmooth initial data.Finally,some numerical results are presented to verify the established theoretical results and show the performances of the developed numerical scheme.展开更多
We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the...We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.展开更多
A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncatio...A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.展开更多
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.展开更多
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.展开更多
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.展开更多
A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segme...A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.展开更多
In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local trunc...In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.展开更多
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)).展开更多
Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations. For time discretization, a three-stage explicit Runge-Kutta method with TVD property is used at pr...Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations. For time discretization, a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage, a cubic spline function is adopted at correcting stage, which made the time discretization accuracy up to fourth order; For spatial discretization, a three-point explicit compact difference scheme with arbitrary order accuracy is employed. The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme. The numerical results agree well with the experimental data. At the same time, the comparisons of the two numerical results between the present scheme and low accuracy difference method are made, which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations. As a valid sample, the wave propagation on the rectangular step is formulated by the present scheme, the modelled results are in better agreement with the experimental data than those of Kittitanasuan.展开更多
The time distances of the completely square-conservative difference schemes in an explicit way are still limited by the stability criteria, though the schemes keep the advantages of traditional square-conservative sch...The time distances of the completely square-conservative difference schemes in an explicit way are still limited by the stability criteria, though the schemes keep the advantages of traditional square-conservative schemes systhetically. How to lengthen the time distances is the central problem in this note. To solve this problem, the aforesaid explicit schemes are improved to be unconditionally stable in theory, according to the thought of upstream differ-展开更多
A kind of explicit square-conserving scheme is proposed for the Landau-Lifshitz equation with Gilbert component. The basic idea was to semidiscrete the Landau-Lifshitz equation into the ordinary differential equation...A kind of explicit square-conserving scheme is proposed for the Landau-Lifshitz equation with Gilbert component. The basic idea was to semidiscrete the Landau-Lifshitz equation into the ordinary differential equations. Then the Lie group method and the Runge-Kutta (RK) method were applied to the ordinary differential equations. The square conserving property and the accuracy of the two methods were compared. Numerical experiment results show the Lie group method has the good accuracy and the square conserving property than the RK method.展开更多
In this paper. a three explicit difference shcemes with high order accuracy for solving the equations of two-dimensional parabolic type is proposed. The stability condition is r=△t/△x￣ 2=△t/△y￣2≤1/4 and the...In this paper. a three explicit difference shcemes with high order accuracy for solving the equations of two-dimensional parabolic type is proposed. The stability condition is r=△t/△x￣ 2=△t/△y￣2≤1/4 and the truncation error is O (△t￣2 + △x￣4 ).展开更多
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 class of two-level high-order accuracy explicit difference scheme for solving 3-D parabolic P.D.E is constructed. Its truncation error is (Δt2+Δx4) and the stability condition is r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/6.
基金the support from the University Transportation Center for Underground Transportation Infrastructure at the Colorado School of Mines for partially funding this research under Grant No. 69A3551747118 of the Fixing America's Surface Transportation Act (FAST Act) of U.S. DoT FY2016
文摘Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.
文摘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.
文摘It is known that the Allen-Chan equations satisfy the maximum principle. Is this true for numerical schemes？ To the best of our knowledge, the state-of-art stability framework is the nonlinear energy stability which has been studied extensively for the phase field type equations. In this work, we will show that a stronger stability under the infinity norm can be established for the implicit-explicit discretization in time and central finite difference in space. In other words, this commonly used numerical method for the Allen-Cahn equation preserves the maximum principle.
基金Partly supported by the State Major Key Project for Basic Researches of China
文摘In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with adjustable time intervals, which is proved to be effective and remarkably time-saving in numerical tests and applications.
文摘In this article,a Crank-Nicolson/Explicit scheme is designed and analyzed for the time-dependent natural convection problem with nonsmooth initial data.The Galerkin finite element method(FEM)with stable MINI element is used for the velocity and pressure and linear polynomial for the temperature.The time discretization is based on the Crank-Nicolson scheme.In order to simplify the computations,the nonlinear terms are treated by the explicit scheme.The advantages of our numerical scheme can be list as follows:(1)The original problem is split into two linear subproblems,these subproblems can be solved in each time level in parallel and the computational sizes are smaller than the origin one.(2)A constant coefficient linear discrete algebraic system is obtained in each subproblem and the computation becomes easy.The main contributions of this work are the stability and convergence results of numerical solutions with nonsmooth initial data.Finally,some numerical results are presented to verify the established theoretical results and show the performances of the developed numerical scheme.
基金funded by the SNF project 200020_204917 entitled"Structure preserving and fast methods for hyperbolic systems of conservation laws".
文摘We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.
基金NSF of the Education Department of Henan Province(20031100010)
文摘A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.
基金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 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 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.
基金Project supported by the National Natural Science Foundation of China(No.10671113)the Natural Science Foundation of Shandong Province of China(No.Y2003A04)
文摘A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.
文摘In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.
文摘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)).
基金The project was financially supported by the National Natural Science Foundation of China (Grant No50479053)
文摘Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations. For time discretization, a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage, a cubic spline function is adopted at correcting stage, which made the time discretization accuracy up to fourth order; For spatial discretization, a three-point explicit compact difference scheme with arbitrary order accuracy is employed. The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme. The numerical results agree well with the experimental data. At the same time, the comparisons of the two numerical results between the present scheme and low accuracy difference method are made, which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations. As a valid sample, the wave propagation on the rectangular step is formulated by the present scheme, the modelled results are in better agreement with the experimental data than those of Kittitanasuan.
基金Project partly supported by the State Major Key Project for Basic Researches.
文摘The time distances of the completely square-conservative difference schemes in an explicit way are still limited by the stability criteria, though the schemes keep the advantages of traditional square-conservative schemes systhetically. How to lengthen the time distances is the central problem in this note. To solve this problem, the aforesaid explicit schemes are improved to be unconditionally stable in theory, according to the thought of upstream differ-
文摘A kind of explicit square-conserving scheme is proposed for the Landau-Lifshitz equation with Gilbert component. The basic idea was to semidiscrete the Landau-Lifshitz equation into the ordinary differential equations. Then the Lie group method and the Runge-Kutta (RK) method were applied to the ordinary differential equations. The square conserving property and the accuracy of the two methods were compared. Numerical experiment results show the Lie group method has the good accuracy and the square conserving property than the RK method.
文摘In this paper. a three explicit difference shcemes with high order accuracy for solving the equations of two-dimensional parabolic type is proposed. The stability condition is r=△t/△x￣ 2=△t/△y￣2≤1/4 and the truncation error is O (△t￣2 + △x￣4 ).
基金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 class of two-level high-order accuracy explicit difference scheme for solving 3-D parabolic P.D.E is constructed. Its truncation error is (Δt2+Δx4) and the stability condition is r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/6.
