A Reduced Basis Method For Fractional Diffusion Operators II

05/07/2020
by   Tobias Danczul, et al.
0

We present a novel numerical scheme to approximate the solution map s u(s) := L^-sf to partial differential equations involving fractional elliptic operators. Reinterpreting L^-s as interpolation operator allows us to derive an integral representation of u(s) which includes solutions to parametrized reaction-diffusion problems. We propose a reduced basis strategy on top of a finite element method to approximate its integrand. Unlike prior works, we deduce the choice of snapshots for the reduced basis procedure analytically. Avoiding further discretization, the integral is interpreted in a spectral setting to evaluate the surrogate directly. Its computation boils down to a matrix approximation L of the operator whose inverse is projected to a low-dimensional space, where explicit diagonalization is feasible. The universal character of the underlying s-independent reduced space allows the approximation of (u(s))_s∈(0,1) in its entirety. We prove exponential convergence rates and confirm the analysis with a variety of numerical examples. Further improvements are proposed in the second part of this investigation to avoid inversion of L. Instead, we directly project the matrix to the reduced space, where its negative fractional power is evaluated. A numerical comparison with the predecessor highlights its competitive performance.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/26/2021

On Rational Krylov and Reduced Basis Methods for Fractional Diffusion

We establish an equivalence between two classes of methods for solving f...
research
11/11/2020

Exponential Convergence of hp FEM for Spectral Fractional Diffusion in Polygons

For the spectral fractional diffusion operator of order 2s∈ (0,2) in bou...
research
05/26/2019

Approximation of a fractional power of an elliptic operator

Some mathematical models of applied problems lead to the need of solving...
research
06/28/2023

Randomized Reduced Basis Methods for Parameterized Fractional Elliptic PDEs

This paper is interested in developing reduced order models (ROMs) for r...
research
10/13/2020

Approximation of Integral Fractional Laplacian and Fractional PDEs via sinc-Basis

Fueled by many applications in random processes, imaging science, geophy...
research
03/24/2021

A Unified Rational Krylov Method for Elliptic and Parabolic Fractional Diffusion Problems

We present a unified framework to efficiently approximate solutions to f...
research
09/02/2021

Matrix-oriented FEM formulation for stationary and time-dependent PDEs on x-normal domains

When numerical solution of elliptic and parabolic partial differential e...

Please sign up or login with your details

Forgot password? Click here to reset