A kernel compression method for distributed-order fractional partial differential equations
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , |
|---|---|
| Format: | Preprint |
| Publié: |
2025
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866911111060652032 |
|---|---|
| author | Beddrich, Jonas Wohlmuth, Barbara |
| author_facet | Beddrich, Jonas Wohlmuth, Barbara |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_13631 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A kernel compression method for distributed-order fractional partial differential equations Beddrich, Jonas Wohlmuth, Barbara Numerical Analysis 34A08, 35R11, 65R10, 74S40 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. |
| title | A kernel compression method for distributed-order fractional partial differential equations |
| topic | Numerical Analysis 34A08, 35R11, 65R10, 74S40 |
| url | https://arxiv.org/abs/2508.13631 |