An accurate parallel finite difference method for the variable-order fractional Laplacian in heterogeneous media
Abstract
Anomalous diffusion in heterogeneous media is ubiquitous across physical and biological systems, where spatially varying microstructures give rise to distinct transport regimes. Variable-order nonlocal operators provide a natural modeling framework for describing heterogeneity and transitions between diffusion regimes. However, the position-dependent long-range interactions inherent in these operators impose prohibitive memory costs and make numerical methods developed for constant-order operators computationally impractical. Consequently, accurate and computationally efficient numerical methods for variable-order nonlocal operators remain scarce. In this paper, we propose accurate finite difference discretizations for the variable-order fractional Laplacian (−Δ)α(x)/2 and develop a fast Fourier transform (FFT)-accelerated scheme for their efficient implementation. The computational and memory costs are O (MNlog(2N)) and O(MN), respectively, where N is the total number of spatial points and M ≪ N. The method is also naturally parallelizable. We establish its consistency and accuracy through a rigorous error analysis that quantifies the local interplay between solution regularity and the spatially varying order α(x). To the best of our knowledge, this is the first work to quantify how these two factors jointly determine convergence rates, with numerical experiments demonstrating this dependence. Extensive numerical experiments confirm the accuracy and efficiency of the proposed method. Applications to anomalous diffusion in composite materials further demonstrate its ability to capture the effects of material heterogeneity.
Related
- On Faber-Krahn inequality for Dirichlet Fractional-Logarithmic Laplacian
- A Note on the Smoluchowski Coagulation-Fragmentation Equation with Fractional Laplacian
- Optimal design problems for the fractional g-Laplacian at prescribed energy level
- Almost L∞-stability of the Ritz projection for the integral fractional Laplacian with 0<s≤ 1/2
- Finite difference methods for a continuous-time heterogeneous agent model with recursive utility
BibTeX
@misc{zhou2026accurate,
title = {{An accurate parallel finite difference method for the variable-order fractional Laplacian in heterogeneous media}},
author = {Shiping Zhou and Yanzhi Zhang},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.17063},
eprint = {2610.17063},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.17063},
}
RIS
TY - UNPB
TI - An accurate parallel finite difference method for the variable-order fractional Laplacian in heterogeneous media
AU - Shiping Zhou
AU - Yanzhi Zhang
PY - 2026
DA - 2026/10/12
PB - arXiv
DO - 10.48550/arxiv.2610.17063
UR - https://arxiv.org/abs/2610.17063
AB - Anomalous diffusion in heterogeneous media is ubiquitous across physical and biological systems, where spatially varying microstructures give rise to distinct transport regimes. Variable-order nonlocal operators provide a natural modeling framework for describing heterogeneity and transitions between diffusion regimes. However, the position-dependent long-range interactions inherent in these operators impose prohibitive memory costs and make numerical methods developed for constant-order operators computationally impractical. Consequently, accurate and computationally efficient numerical methods for variable-order nonlocal operators remain scarce. In this paper, we propose accurate finite difference discretizations for the variable-order fractional Laplacian (−Δ)^α(x)/2 and develop a fast Fourier transform (FFT)-accelerated scheme for their efficient implementation. The computational and memory costs are O (MNlog(2N)) and O(MN), respectively, where N is the total number of spatial points and M ≪ N. The method is also naturally parallelizable. We establish its consistency and accuracy through a rigorous error analysis that quantifies the local interplay between solution regularity and the spatially varying order α(x). To the best of our knowledge, this is the first work to quantify how these two factors jointly determine convergence rates, with numerical experiments demonstrating this dependence. Extensive numerical experiments confirm the accuracy and efficiency of the proposed method. Applications to anomalous diffusion in composite materials further demonstrate its ability to capture the effects of material heterogeneity.
ER -