5 research outputs found
A second-order and nonuniform time-stepping maximum-principle preserving scheme for time-fractional Allen-Cahn equations
In this work, we present a second-order nonuniform time-stepping scheme for
the time-fractional Allen-Cahn equation. We show that the proposed scheme
preserves the discrete maximum principle, and by using the convolution
structure of consistency error, we present sharp maximum-norm error estimates
which reflect the temporal regularity. As our analysis is built on nonuniform
time steps, we may resolve the intrinsic initial singularity by using the
graded meshes. Moreover, we propose an adaptive time-stepping strategy for
large time simulations. Numerical experiments are presented to show the
effectiveness of the proposed scheme. This seems to be the first second-order
maximum principle preserving scheme for the time-fractional Allen-Cahn
equation.Comment: 22pages, 22 figures, 2 table
A discrete Gr\"{o}nwall inequality with application to numerical schemes for subdiffusion problems
We consider a class of numerical approximations to the Caputo fractional
derivative. Our assumptions permit the use of nonuniform time steps, such as is
appropriate for accurately resolving the behavior of a solution whose
derivatives are singular at~. The main result is a type of fractional
Gr\"{o}nwall inequality and we illustrate its use by outlining some stability
and convergence estimates of schemes for fractional reaction-subdiffusion
problems. This approach extends earlier work that used the familiar L1
approximation to the Caputo fractional derivative, and will facilitate the
analysis of higher order and linearized fast schemes.Comment: 15 pages, 2 figure
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
Sharp -norm error estimates of two time-stepping schemes for reaction-subdiffusion problems
Due to the intrinsically initial singularity of solution and the discrete
convolution form in numerical Caputo derivatives, the traditional -norm
analysis (corresponding to the case for a classical diffusion equation) to the
time approximations of a fractional subdiffusion problem always leads to
suboptimal error estimates (a loss of time accuracy). To recover the
theoretical accuracy in time, we propose an improved discrete Gr\"{o}nwall
inequality and apply it to the well-known L1 formula and a fractional
Crank-Nicolson scheme. With the help of a time-space error-splitting technique
and the global consistency analysis, sharp -norm error estimates of the
two nonuniform approaches are established for a reaction-subdiffusion problems.
Numerical experiments are included to confirm the sharpness of our analysis.Comment: 22 pages, 8 table
A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem
Stability and convergence of a time-weighted discrete scheme with nonuniform
time steps are established for linear reaction-subdiffusion equations. The
Caupto derivative is approximated at an offset point by using linear and
quadratic polynomial interpolation. Our analysis relies on two tools: a
discrete fractional Gr\"{o}nwall inequality and the global consistency
analysis. The new consistency analysis makes use of an interpolation error
formula for quadratic polynomials, which leads to a convolution-type bound for
the local truncation error. To exploit these two tools, some theoretical
properties of the discrete kernels in the numerical Caputo formula are crucial
and we investigate them intensively in the nonuniform setting. Taking the
initial singularity of the solution into account, we obtain a sharp error
estimate on nonuniform time meshes. The fully discrete scheme generates a
second-order accurate solution on the graded mesh provided a proper grading
parameter is employed. An example is presented to show the sharpness of our
analysis.Comment: 23 pages, 4 table