In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two ...In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.展开更多
The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part o...The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.展开更多
The improved discrete velocity method(IDVM),a multiscale kinetic approach,solves both the Boltzmann model equation and the macroscopic governing equations.It maintains the simplicity of the conventional discrete veloc...The improved discrete velocity method(IDVM),a multiscale kinetic approach,solves both the Boltzmann model equation and the macroscopic governing equations.It maintains the simplicity of the conventional discrete velocity method while significantly enhancing accuracy and efficiency in the continuum flow regime.However,previous implementations of IDVM were limited to monatomic gases,rendering it inapplicable under real atmospheric conditions dominated by diatomic gases.In this paper,the fully implicit IDVM is extended to the Boltzmann-Rykov model equation to simulate multiscale gas flows involving the rotational non-equilibrium effects of diatomic molecules.This approach incorporates macroscopic governing equations to predict the equilibrium distribution function,enabling the fully implicit discretization of the BoltzmannRykov model equation and ensuring rapid convergence.Numerical simulations are performed for several cases,including the one-dimensional nitrogen shock tube,two-dimensional flows around a flat plate and a blunt circular cylinder,and threedimensional supersonic flow over a sphere.The results of the present algorithm show remarkable alignment with those from similar kinetic algorithms and experimental data,demonstrating its ability to simulate non-equilibrium diatomic gas flows.Relative to the conventional semi-implicit discrete velocity method and implicit kinetic methods with a macroscopic forecasting technique,the proposed approach attains varying degrees of acceleration.展开更多
A monotone compact implicit finite difference scheme with fourth-order accuracy in space and second-order in time is proposed for solving nonlinear reaction-diffusion equations.An accelerated monotone iterative method...A monotone compact implicit finite difference scheme with fourth-order accuracy in space and second-order in time is proposed for solving nonlinear reaction-diffusion equations.An accelerated monotone iterative method for the resulting discrete problem is presented.The sequence of iteration converges monotonically to the unique solution of the discrete problem,and the convergence rate is either quadratic or nearly quadratic,depending on the property of the nonlinear reaction.The numerical results illustrate the high accuracy of the proposed scheme and the rapid convergence rate of.the iteration.展开更多
A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed fro...A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed from given macroscopic variables. Based on this fact, we propose a technique to reconstruct distribution function at discrete level, and employ it to develop an implicit numerical method for kinetic equations. The new implicit method only stores the macroscopic quantities which appear in the collision term, and does not store the distribution functions. As a result, enormous memory requirement for solving kinetic equations is totally relieved. Several boundary conditions, such as, inlet, outlet and isothermal boundaries, are discussed. Some numerical tests demonstrate the validity and efficiency of the technique.The new implicit solver provides nearly identical solution as the explicit kinetic solver, while the memory requirement is on the same order as the Navier–Stokes solver.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
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...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 deals with an efficient two-step time split explicit/implicit scheme applied to a two-dimensional nonlinear unsteady convection-diffusionreaction equation.The computational cost of the new algorithm at each...This paper deals with an efficient two-step time split explicit/implicit scheme applied to a two-dimensional nonlinear unsteady convection-diffusionreaction equation.The computational cost of the new algorithm at each time level is equivalent to solving a pentadiagonalmatrix equation with strictly dominant diagonal elements.Such a bandwidth matrix can be easily inverted using the Gaussian Decomposition and the corresponding linear system should be solved by the back substitutionmethod.The proposed approach is unconditionally stable,temporal second-order accuracy and fourth-order convergence in space.These results suggest that the developed technique is faster and more efficient than a large class of numerical methods studied in the literature for the considered initial-boundary value problem.Numerical experiments are carried out to confirm the theoretical analysis and to demonstrate the performance of the constructed numerical scheme.展开更多
Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme...Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.展开更多
An implicit scheme discretized by the finite element method is developed for the numerical solution of the time-dependent head/disk dynamics. The scheme is second-order timeaccurate, which offers the flexibility of em...An implicit scheme discretized by the finite element method is developed for the numerical solution of the time-dependent head/disk dynamics. The scheme is second-order timeaccurate, which offers the flexibility of employing time steps compatible with the system dynamics while minimizing the numerical instability presented in most explicit methods. Good agreement is orreryed between the present calculations and D. K. Miu's work.展开更多
Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly im...Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.展开更多
Many interesting applications of hyperbolic systems of equations are stiff,and require the time step to satisfy restrictive stability conditions.One way to avoid small time steps is to use implicit time integration.Im...Many interesting applications of hyperbolic systems of equations are stiff,and require the time step to satisfy restrictive stability conditions.One way to avoid small time steps is to use implicit time integration.Implicit integration is quite straightforward for first-order schemes.High order schemes instead also need to control spurious oscillations,which requires limiting in space and time also in the linear case.We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor,which is used both to set up the nonlinear weights of a standard high order space reconstruction,and to achieve limiting in time.In this preliminary work,we concentrate on the case of a third-order scheme,based on diagonally implicit Runge Kutta(DIRK)integration in time and central weighted essentially non-oscillatory(CWENO)reconstruction in space.The numerical tests involve linear and nonlinear scalar conservation laws.展开更多
Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid.A drawback of such an explicit app...Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid.A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability.An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed.This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.展开更多
The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity. To solve this problem, implicit schemes should be ...The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity. To solve this problem, implicit schemes should be adopted since explicit schemes may bring about a great increase in computation quantity according to the Courant-FrledrichsLewy condition. Whereas the adoption of implicit schemes is difficult to be realized because of the existence of two partial derivatives of surface elevations with respect to variables of alternative direction coordinates in each momentum equation in non-rectangular coordinates. With an aim to design an implicit scheme in non-reetangular ccordinates in the present paper, new momentum equations with the contravariant components of velocity vector are derived based on the shallow water dynamic equations in generalized curvilinear coordinates. In each equation, the coefficients before the two detivatives of surface elevations have different orders of magnitude, i. e., the derivative with the larger ceefficient rnay play a more important role than that with the smaller one. With this advantage, the ADI scheme can then be easily employed to improve the numerical stability and decrease the calculating quantity. The calculation in a harbour and a channel in Macau nearshore area shows that the implicit model is effective in calculating current fields in small size areas.展开更多
High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discont...High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discontinuous Galerkin(DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry representation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomerated high-order mesh works well with DG method in dealing with flow problems of curved geometries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant efficiency improvements compared to its explicit counterparts.展开更多
The tempered fractional calculus has been successfully applied for depicting the time evolution of a system describing non-Markovian diffusion particles.The related governing equations are a series of partial differen...The tempered fractional calculus has been successfully applied for depicting the time evolution of a system describing non-Markovian diffusion particles.The related governing equations are a series of partial differential equations with tempered fractional derivatives.Using the polynomial interpolation technique,in this paper,we present three efficient numerical formulas,namely the tempered L1 formula,the tempered L1-2 formula,and the tempered L2-1σformula,to approximate the Caputo-tempered fractional derivative of orderα∈(0,1).The truncation error of the tempered L1 formula is of order 2-α,and the tempered L1-2 formula and L2-1σformula are of order 3-α.As an application,we construct implicit schemes and implicit ADI schemes for one-dimensional and two-dimensional time-tempered fractional diffusion equations,respectively.Furthermore,the unconditional stability and convergence of two developed difference schemes with tempered L1 and L2-1σformulas are proved by the Fourier analysis method.Finally,we provide several numerical examples to demonstrate the correctness and effectiveness of the theoretical analysis.展开更多
This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ2+hx4+hy4),where τ,h_(x )and hy are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
摘要In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.
摘要The aim of this paper is to give an appropriate numerical method to solve Allen-Cahn equation, with Dirichlet or Neumann boundary condition. The time discretization involves an explicit scheme for the nonlinear part of the operator and an implicit Euler discretization of the linear part. Finite difference schemes are used for the spatial part. This finally leads to the numerical solution of a sparse linear system that can be solved efficiently.
基金support by the National Key Basic Research Projects(Grant No.2022JCJQZD20600).
摘要The improved discrete velocity method(IDVM),a multiscale kinetic approach,solves both the Boltzmann model equation and the macroscopic governing equations.It maintains the simplicity of the conventional discrete velocity method while significantly enhancing accuracy and efficiency in the continuum flow regime.However,previous implementations of IDVM were limited to monatomic gases,rendering it inapplicable under real atmospheric conditions dominated by diatomic gases.In this paper,the fully implicit IDVM is extended to the Boltzmann-Rykov model equation to simulate multiscale gas flows involving the rotational non-equilibrium effects of diatomic molecules.This approach incorporates macroscopic governing equations to predict the equilibrium distribution function,enabling the fully implicit discretization of the BoltzmannRykov model equation and ensuring rapid convergence.Numerical simulations are performed for several cases,including the one-dimensional nitrogen shock tube,two-dimensional flows around a flat plate and a blunt circular cylinder,and threedimensional supersonic flow over a sphere.The results of the present algorithm show remarkable alignment with those from similar kinetic algorithms and experimental data,demonstrating its ability to simulate non-equilibrium diatomic gas flows.Relative to the conventional semi-implicit discrete velocity method and implicit kinetic methods with a macroscopic forecasting technique,the proposed approach attains varying degrees of acceleration.
基金supported in part by NSF of China No.10571059E-Institutes of Shanghai Municipal Education Commission No.E03004+4 种基金Shanghai Priority Academic Discipline,and the Scientific Research Foundation for the Returned Overseas Chinese Scholars of the State Education MinistrySF of Shanghai No.04JC14062the fund of Chinese Education Ministry No.20040270002the Shanghai Leading Academic Discipline Project No.T0401the fund for E-Institutes of Shanghai Municipal Education Commission No.E03004 and the fund No.04DB15 of Shanghai Municipal Education Commission
摘要A monotone compact implicit finite difference scheme with fourth-order accuracy in space and second-order in time is proposed for solving nonlinear reaction-diffusion equations.An accelerated monotone iterative method for the resulting discrete problem is presented.The sequence of iteration converges monotonically to the unique solution of the discrete problem,and the convergence rate is either quadratic or nearly quadratic,depending on the property of the nonlinear reaction.The numerical results illustrate the high accuracy of the proposed scheme and the rapid convergence rate of.the iteration.
基金supported by the National Natural Science Foundation of China(11602091 and 91530319)the National Key Research and Development Plan(2016YFB0600805)
摘要A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed from given macroscopic variables. Based on this fact, we propose a technique to reconstruct distribution function at discrete level, and employ it to develop an implicit numerical method for kinetic equations. The new implicit method only stores the macroscopic quantities which appear in the collision term, and does not store the distribution functions. As a result, enormous memory requirement for solving kinetic equations is totally relieved. Several boundary conditions, such as, inlet, outlet and isothermal boundaries, are discussed. Some numerical tests demonstrate the validity and efficiency of the technique.The new implicit solver provides nearly identical solution as the explicit kinetic solver, while the memory requirement is on the same order as the Navier–Stokes solver.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
摘要Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
摘要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 deals with an efficient two-step time split explicit/implicit scheme applied to a two-dimensional nonlinear unsteady convection-diffusionreaction equation.The computational cost of the new algorithm at each time level is equivalent to solving a pentadiagonalmatrix equation with strictly dominant diagonal elements.Such a bandwidth matrix can be easily inverted using the Gaussian Decomposition and the corresponding linear system should be solved by the back substitutionmethod.The proposed approach is unconditionally stable,temporal second-order accuracy and fourth-order convergence in space.These results suggest that the developed technique is faster and more efficient than a large class of numerical methods studied in the literature for the considered initial-boundary value problem.Numerical experiments are carried out to confirm the theoretical analysis and to demonstrate the performance of the constructed numerical scheme.
摘要Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.
摘要An implicit scheme discretized by the finite element method is developed for the numerical solution of the time-dependent head/disk dynamics. The scheme is second-order timeaccurate, which offers the flexibility of employing time steps compatible with the system dynamics while minimizing the numerical instability presented in most explicit methods. Good agreement is orreryed between the present calculations and D. K. Miu's work.
摘要Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.
基金MIUR(Ministry of University and Research)PRIN2017 project number 2017KKJP4XProgetto di Ateneo Sapienza,number RM120172B41DBF3A.
摘要Many interesting applications of hyperbolic systems of equations are stiff,and require the time step to satisfy restrictive stability conditions.One way to avoid small time steps is to use implicit time integration.Implicit integration is quite straightforward for first-order schemes.High order schemes instead also need to control spurious oscillations,which requires limiting in space and time also in the linear case.We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor,which is used both to set up the nonlinear weights of a standard high order space reconstruction,and to achieve limiting in time.In this preliminary work,we concentrate on the case of a third-order scheme,based on diagonally implicit Runge Kutta(DIRK)integration in time and central weighted essentially non-oscillatory(CWENO)reconstruction in space.The numerical tests involve linear and nonlinear scalar conservation laws.
基金supported by the US National Science Foundation(DMS-0713718)
摘要Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid.A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability.An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed.This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.
摘要The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity. To solve this problem, implicit schemes should be adopted since explicit schemes may bring about a great increase in computation quantity according to the Courant-FrledrichsLewy condition. Whereas the adoption of implicit schemes is difficult to be realized because of the existence of two partial derivatives of surface elevations with respect to variables of alternative direction coordinates in each momentum equation in non-rectangular coordinates. With an aim to design an implicit scheme in non-reetangular ccordinates in the present paper, new momentum equations with the contravariant components of velocity vector are derived based on the shallow water dynamic equations in generalized curvilinear coordinates. In each equation, the coefficients before the two detivatives of surface elevations have different orders of magnitude, i. e., the derivative with the larger ceefficient rnay play a more important role than that with the smaller one. With this advantage, the ADI scheme can then be easily employed to improve the numerical stability and decrease the calculating quantity. The calculation in a harbour and a channel in Macau nearshore area shows that the implicit model is effective in calculating current fields in small size areas.
基金co-supported by the Aeronautical Science Foundation of China (No. 20152752033)the National Natural Science Foundation of China (No. 11272152)the Open Project of Key Laboratory of Aerodynamic Noise Control
摘要High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discontinuous Galerkin(DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry representation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomerated high-order mesh works well with DG method in dealing with flow problems of curved geometries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant efficiency improvements compared to its explicit counterparts.
摘要The tempered fractional calculus has been successfully applied for depicting the time evolution of a system describing non-Markovian diffusion particles.The related governing equations are a series of partial differential equations with tempered fractional derivatives.Using the polynomial interpolation technique,in this paper,we present three efficient numerical formulas,namely the tempered L1 formula,the tempered L1-2 formula,and the tempered L2-1σformula,to approximate the Caputo-tempered fractional derivative of orderα∈(0,1).The truncation error of the tempered L1 formula is of order 2-α,and the tempered L1-2 formula and L2-1σformula are of order 3-α.As an application,we construct implicit schemes and implicit ADI schemes for one-dimensional and two-dimensional time-tempered fractional diffusion equations,respectively.Furthermore,the unconditional stability and convergence of two developed difference schemes with tempered L1 and L2-1σformulas are proved by the Fourier analysis method.Finally,we provide several numerical examples to demonstrate the correctness and effectiveness of the theoretical analysis.
摘要This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ2+hx4+hy4),where τ,h_(x )and hy are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.