A kernel compression method for distributed-order fractional partial differential equations
By: Jonas Beddrich, Barbara Wohlmuth
Potential Business Impact:
Solves hard math problems faster using simpler steps.
We propose a kernel compression method for solving Distributed-Order (DO) Fractional Partial Differential Equations (DOFPDEs) at the cost of solving corresponding local-in-time PDEs. The key concepts are (1) discretization of the integral over the order of the fractional derivative and (2) approximation of linear combinations of integral kernels with exponential sums, expressing the non-local history term as a sum of auxiliary variables that solve a weakly coupled, local in time system of PDEs. For the second step, we introduce an improved algorithm that approximates the occurring integral kernels with double precision accuracy using only a moderate number (<100) of exponential terms. After temporal discretization using implicit Runge--Kutta methods, we exploit the inherent structure of the PDE system to obtain the solution at each time step by solving a single PDE. At the same time, the auxiliary variables are computed by a linear update, not even requiring a matrix-vector multiplication. Choosing temporal meshes with a grading factor corresponding to the convergence order of the Runge--Kutta schemes, we achieve the optimal decay of the temporal discretization error. The flexibility and robustness of our numerical scheme are illustrated by recreating well-studied test cases and solving linear and nonlinear DOFPDEs in 2D and 3D with up to 40 million spatial degrees of freedom.
Similar Papers
A novel class of arbitrary high-order numerical schemes for fractional differential equations
Numerical Analysis
Solves hard math problems faster and with less computer memory.
A fast and memoryless numerical method for solving fractional differential equations
Numerical Analysis
Solves hard math problems for science and engineering.
Vectorised Parallel in Time methods for low-order discretizations with application to Porous Media problems
Numerical Analysis
Speeds up computer simulations of tricky science problems.