Search

In this paper, we study the two-dimensional Volterra integro-differential equations for viscoelastic rods and membranes in a bounded smooth domain. The memory kernel of this equation makes it not easy to construct an efficient numerical scheme and perform theoretical analysis. The numerical method is considered by the finite difference approach for spatial discretization and Crank–Nicolson (CN) alternating direction implicit (ADI) scheme in the time direction. The integral terms are approximated by the fractional convolution quadrature. We prove that the proposed method is unconditional stable and derive error estimates in L 2 norm. The literature reported on the ADI finite difference method of this model is extremely sparse. Numerical results support the theoretical analysis. [ABSTRACT FROM AUTHOR]

In this paper, a high order finite difference scheme for a two dimensional nonlinear fractional Klein–Gordon equation in polar coordinates is proposed. The difficulty caused by the variable coefficient and the nonlinear term in the equation is carefully handled. Stability and convergence are analyzed using the matrix form of the scheme. Numerical examples are provided to demonstrate the theoretical results. [ABSTRACT FROM AUTHOR]

In this paper, a high order compact (HOC) scheme for time fractional Fokker–Planck equations with variable convection is constructed. The scheme is studied using its matrix form by the energy method. We find that the difficulty arising from the variable coefficient can be overcome by simple modifications of the coefficient matrices. The scheme is shown to be stable and convergent with order τ 2 − α + h 4 which is higher than some recently studied schemes. Numerical examples are given to justify the theoretical analysis. [ABSTRACT FROM AUTHOR]

In this paper, we study the numerical method for time-fractional KdV–Burgers' equation with initial singularity. The famous L 2 - 1 σ formula on graded meshes is adopted to approximate the Caputo derivative. Meanwhile, a nonlinear finite difference method on uniform grids is deduced for spatial discretization. The proposed method is second order in time and first order in space. With the help of the fractional Grönwall inequality, the unconditional stability and convergence of the current scheme are analyzed based on some skills. To raise the accuracy in spatial direction, a second order method is then carefully deduced. At last, theoretical results are verified by numerical experiments. [ABSTRACT FROM AUTHOR]