A fast linearized numerical method for nonlinear time-fractional diffusion equations

Numerical Algorithms - In this paper, we study a fast linearized numerical method for solving nonlinear time-fractional diffusion equations. A new weighted method is proposed to construct...     

Stability of fully discrete schemes with interpolation-type fractional formulas for distributed-order subdiffusion equations

Two fully discrete methods are investigated for simulating the distributed-order sub-diffusion equation in Caputo’s form. The fractional Caputo derivative is approximated by the Caputo’s BDF1 (called L1 early) and BDF2 (or L1-2 when it was first introduced) approximations, which are constructed by piecewise linear and quadratic interpolating polynomials, respectively. It is shown that the first scheme, using the BDF1 formula, possesses the discrete minimum-maximum principle and nonnegativity preservation property such that it is stable and convergent in the maximum norm. The method using the BDF2 formula is shown to be stable and convergent in the discrete H ¹ norm by using the discrete energy method. For problems of distributed order within a certain region, the method is also proven to preserve the discrete maximum principle and nonnegativity property. Extensive numerical experiments are provided to show the effectiveness of numerical schemes, and to examine the initial singularity of the solution. The applicability of our numerical algorithms to a problem with solution which lacks the smoothness near the initial time is examined by employing a class of power-type nonuniform meshes.     

High‐order compact schemes for fractional differential equations with mixed derivatives

In this article, we consider two‐dimensional fractional subdiffusion equations with mixed derivatives. A high‐order compact scheme is proposed to solve the problem. We establish a sufficient condition and show that the scheme converges with fourth order in space and second order in time under this condition.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2141–2158, 2017     

A high‐order compact scheme for the nonlinear fractional Klein–Gordon equation

In this article, a high‐order finite difference scheme for a kind of nonlinear fractional Klein–Gordon equation is derived. The time fractional derivative is described in the Caputo sense. The solvability of the difference system is discussed by the Leray–Schauder fixed point theorem, while the stability and L_∞ convergence of the finite difference scheme are proved by the energy method. Numerical examples are provided to demonstrate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 706–722, 2015     

A compact ADI scheme for the two dimensional time fractional diffusion‐wave equation in polar coordinates

In this article, we consider finite difference schemes for two dimensional time fractional diffusion‐wave equations on an annular domain. The problem is formulated in polar coordinates and, therefore, has variable coefficients. A compact alternating direction implicit scheme with accuracy order is derived, where τ, h₁, h₂ are the temporal and spatial step sizes, respectively. The stability and convergence of the proposed scheme are studied using its matrix form by the energy method. Numerical experiments are presented to support the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1692–1712, 2015     

A nonuniform L2 formula of Caputo derivative and its application to a fractional Benjamin–Bona–Mahony-type equation with nonsmooth solutions

To recover the full accuracy of discretized fractional derivatives, nonuniform mesh technique is a natural and simple approach to efficiently resolve the initial singularities that always appear in the solutions of time-fractional linear and nonlinear differential equations. We first construct a nonuniform L2 approximation for the fractional Caputo's derivative of order 1 < α < 2 and present a global consistency analysis under some reasonable regularity assumptions. The temporal nonuniform L2 formula is then utilized to develop a linearized difference scheme for a time-fractional Benjamin–Bona–Mahony-type equation. The unconditional convergence of our scheme on both uniform and nonuniform (graded) time meshes are proven with respect to the discrete H¹-norm. Numerical examples are provided to justify the accuracy.     

High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

In this work, an effective and fast finite element numerical method with high-order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two-level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H¹-norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.     

A linearized second-order scheme for nonlinear time fractional Klein-Gordon type equations

We consider difference schemes for nonlinear time fractional Klein-Gordon type equations in this paper. A linearized scheme is proposed to solve the problem. As a result, iterative method need not be employed. One of the main difficulties for the analysis is that certain weight averages of the approximated solutions are considered in the discretization and standard energy estimates cannot be applied directly. By introducing a new grid function, which further approximates the solution, and using ideas in some recent studies, we show that the method converges with second-order accuracy in time.     

A two-step modulus-based matrix splitting iteration method for horizontal linear complementarity problems

Numerical Algorithms - In this paper, for solving horizontal linear complementarity problems, a two-step modulus-based matrix splitting iteration method is established. The convergence analysis of...     

A modified modulus-based matrix splitting iteration method for solving implicit complementarity problems

In this paper, a modified modulus-based matrix splitting iteration method is established for solving a class of implicit complementarity problems. The global convergence conditions are given when the system matrix is a positive definite matrix or an H₊-matrix, respectively. In addition, some numerical examples show that the proposed method is efficient.