19 research outputs found
Correction of high-order BDF convolution quadrature for fractional evolution equations
We develop proper correction formulas at the starting steps to restore
the desired -order convergence rate of the -step BDF convolution
quadrature for discretizing evolution equations involving a fractional-order
derivative in time. The desired -order convergence rate can be
achieved even if the source term is not compatible with the initial data, which
is allowed to be nonsmooth. We provide complete error estimates for the
subdiffusion case , and sketch the proof for the
diffusion-wave case . Extensive numerical examples are provided
to illustrate the effectiveness of the proposed scheme.Comment: 22 pages, 3 figure
Stability analysis of a second-order difference scheme for the time-fractional mixed sub-diffusion and diffusion-wave equation
This study investigates a class of initial-boundary value problems pertaining
to the time-fractional mixed sub-diffusion and diffusion-wave equation (SDDWE).
To facilitate the development of a numerical method and analysis, the original
problem is transformed into a new integro-differential model which includes the
Caputo derivatives and the Riemann-Liouville fractional integrals with orders
belonging to (0,1). By providing an a priori estimate of the solution, we have
established the existence and uniqueness of a numerical solution for the
problem. We propose a second-order method to approximate the fractional
Riemann-Liouville integral and employ an L2 type formula to approximate the
Caputo derivative. This results in a method with a temporal accuracy of
second-order for approximating the considered model. The proof of the
unconditional stability of the proposed difference scheme is established.
Moreover, we demonstrate the proposed method's potential to construct and
analyze a second-order L2-type numerical scheme for a broader class of the
time-fractional mixed SDDWEs with multi-term time-fractional derivatives.
Numerical results are presented to assess the accuracy of the method and
validate the theoretical findings
Numerical methods for time-fractional evolution equations with nonsmooth data: a concise overview
Over the past few decades, there has been substantial interest in evolution
equations that involving a fractional-order derivative of order
in time, due to their many successful applications in
engineering, physics, biology and finance. Thus, it is of paramount importance
to develop and to analyze efficient and accurate numerical methods for reliably
simulating such models, and the literature on the topic is vast and fast
growing. The present paper gives a concise overview on numerical schemes for
the subdiffusion model with nonsmooth problem data, which are important for the
numerical analysis of many problems arising in optimal control, inverse
problems and stochastic analysis. We focus on the following aspects of the
subdiffusion model: regularity theory, Galerkin finite element discretization
in space, time-stepping schemes (including convolution quadrature and L1 type
schemes), and space-time variational formulations, and compare the results with
that for standard parabolic problems. Further, these aspects are showcased with
illustrative numerical experiments and complemented with perspectives and
pointers to relevant literature.Comment: 24 pages, 3 figure
Optimal and a posteriori error estimates for the fully discrete approximations of time fractional parabolic differential equations
We derive optimal order a posteriori error estimates in the
and -norms for the fully discrete approximations of time fractional
parabolic differential equations. For the discretization in time, we use the
methods, while for the spatial discretization, we use standard conforming
finite element methods. The linear and quadratic space-time reconstructions are
introduced, which are generalizations of the elliptic space reconstruction.
Then the related a posteriori error estimates for the linear and quadratic
space-time reconstructions play key roles in deriving global and pointwise
final error estimates. Numerical experiments verify and complement our
theoretical results.Comment: 22 page