77 research outputs found
Local discontinuous Galerkin methods for fractional ordinary differential equations
This paper discusses the upwinded local discontinuous Galerkin methods for
the one-term/multi-term fractional ordinary differential equations (FODEs). The
natural upwind choice of the numerical fluxes for the initial value problem for
FODEs ensures stability of the methods. The solution can be computed element by
element with optimal order of convergence in the norm and
superconvergence of order at the downwind point of each
element. Here is the degree of the approximation polynomial used in an
element and () represents the order of the one-term
FODEs. A generalization of this includes problems with classic 'th-term
FODEs, yielding superconvergence order at downwind point as
. The underlying mechanism of the
superconvergence is discussed and the analysis confirmed through examples,
including a discussion of how to use the scheme as an efficient way to evaluate
the generalized Mittag-Leffler function and solutions to more generalized
FODE's.Comment: 17 pages, 7 figure
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
- …