1,083 research outputs found
Numerical Methods for the Fractional Laplacian: a Finite Difference-quadrature Approach
The fractional Laplacian is a non-local operator which
depends on the parameter and recovers the usual Laplacian as . A numerical method for the fractional Laplacian is proposed, based on
the singular integral representation for the operator. The method combines
finite difference with numerical quadrature, to obtain a discrete convolution
operator with positive weights. The accuracy of the method is shown to be
. Convergence of the method is proven. The treatment of far
field boundary conditions using an asymptotic approximation to the integral is
used to obtain an accurate method. Numerical experiments on known exact
solutions validate the predicted convergence rates. Computational examples
include exponentially and algebraically decaying solution with varying
regularity. The generalization to nonlinear equations involving the operator is
discussed: the obstacle problem for the fractional Laplacian is computed.Comment: 29 pages, 9 figure
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
Chebyshev Collocation Method for Parabolic Partial Integrodifferential Equations
An efficient technique for solving parabolic partial integrodifferential equation is presented. This technique is based on Chebyshev polynomials and finite difference method. A priori error estimate for the proposed technique is deduced. Some examples are presented to illustrate the validity and efficiency of the presented method
- ā¦