A Riesz basis Galerkin method for the tempered fractional Laplacian
arXiv:1709.10415 · doi:10.1137/17M1151791
Abstract
The fractional Laplacian is the generator of -stable Lévy process, which is the scaling limit of the Lévy fight. Due to the divergence of the second moment of the jump length of the Lévy fight it is not appropriate as a physical model in many practical applications. However, using a parameter to exponentially temper the isotropic power law measure of the jump length leads to the tempered Lévy fight, which has finite second moment. For short time the tempered Lévy fight exhibits the dynamics of Lévy fight while after sufficiently long time it turns to normal diffusion. The generator of tempered -stable Lévy process is the tempered fractional Laplacian [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., in press, 2017]. In the current work, we present new computational methods for the tempered fractional Laplacian equation, including the cases with the homogeneous and nonhomogeneous generalized Dirichlet type boundary conditions. We prove the well-posedness of the Galerkin weak formulation and provide convergence analysis of the single scaling B-spline and multiscale Riesz bases finite element methods. We propose a technique for efficiently generating the entries of the dense stiffness matrix and for solving the resulting algebraic equation by preconditioning. We also present several numerical experiments to verify the theoretical results.
28 pages, 2 figures
References in corpus (4)
- A short FE implementation for a 2d homogeneous Dirichlet problem of a Fractional Laplacian
- Boundary problems for the fractional and tempered fractional operators
- Variational formulation and efficient implementation for solving the tempered fractional problems
- Numerical schemes of the time tempered fractional Feynman-Kac equation
Cited by in corpus (8)
- Analytic approaches of the anomalous diffusion: a review
- Numerical approximations for the tempered fractional Laplacian: Error analysis and applications
- Tempered Fractional LES Modeling
- Nonlocal Diffusion Operators for Normal and Anomalous Dynamics
- Numerical approximation for fractional diffusion equation forced by a tempered fractional Gaussian noise
- A Multigrid method for nonlocal problems: non-diagonally dominant Toeplitz-plus-tridiagonal systems
- Algorithm implementation and numerical analysis for the two-dimensional tempered fractional Laplacian
- Numerical algorithms of the two-dimensional Feynman-Kac equation for reaction and diffusion processes