522 research outputs found
An hp-Local Discontinuous Galerkin method for Parabolic\ud Integro-Differential Equations
In this article, a priori error analysis is discussed for an hp-local discontinuous Galerkin (LDG) approximation to a parabolic integro-differential equation. It is shown that the L2 -norm of the gradient and the L2 -norm of the potential are optimal in the discretizing parameter h and suboptimal in the degree of polynomial p. Due to the presence of the integral term, an introduction of an expanded mixed type Ritz-Volterra projection helps to achieve optimal estimates. Further, it is observed that a negative norm estimate of the gradient plays a crucial role in our convergence analysis. As in the elliptic case, similar results on order of convergence are established for the semidiscrete method after suitably modifying the numerical fluxes. The optimality of these theoretical results is tested in a series of numerical experiments on two dimensional domains
An hp-version discontinuous Galerkin method for integro-differential equations of parabolic type
We study the numerical solution of a class of parabolic integro-differential equations with weakly singular kernels. We use an -version discontinuous Galerkin (DG) method for the discretization in time. We derive optimal -version error estimates and show that exponential rates of convergence can be achieved for solutions with singular (temporal) behavior near caused by the weakly singular kernel. Moreover, we prove that by using nonuniformly refined time steps, optimal algebraic convergence rates can be achieved for the -version DG method. We then combine the DG time-stepping method with a standard finite element discretization in space, and present an optimal error analysis of the resulting fully discrete scheme. Our theoretical results are numerically validated in a series of test problems
Discontinuous Galerkin method for an integro-differential equation modeling dynamic fractional order viscoelasticity
An integro-differential equation, modeling dynamic fractional order
viscoelasticity, with a Mittag-Leffler type convolution kernel is considered. A
discontinuous Galerkin method, based on piecewise constant polynomials is
formulated for temporal semidiscretization of the problem. Stability estimates
of the discrete problem are proved, that are used to prove optimal order a
priori error estimates. The theory is illustrated by a numerical example.Comment: 16 pages, 2 figure
Optimal L2 estimates for semidiscrete Galerkin methods for\ud parabolic integro-differential equations with nonsmooth data
In this article, we discuss an alternate approach to a priori error estimates for the semidiscrete Galerkin approximation to a time dependent parabolic integro-differential equation with nonsmooth initial data. It is based on energy arguments and on a repeated use of time integration, but without using parabolic type duality technique. Optimal L2-error estimate is derived for the semidiscrete approximation, when the initial data is in L2
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
Artificial boundary conditions for parabolic Volterra integro-differential equations on unbounded two-dimensional domains
AbstractIn this paper we study the numerical solution of parabolic Volterra integro-differential equations on certain unbounded two-dimensional spatial domains. The method is based on the introduction of a feasible artificial boundary and the derivation of corresponding artificial (fully transparent) boundary conditions. Two examples illustrate the application and numerical performance of the method
- ā¦