Structure-Preserving Numerical Schemes for a Time-Fractional Porous Medium Equation with Nonlocal Pressure
Abstract
In this work, we study numerical approximations of the equation ∂tαu−∇·(um−1∇(−Δ)−σu)=0, where 0<α<1, 0<σ<1 and m>1. The Caputo derivative is discretized by the L1 formula, while the nonlinear nonlocal flux is approximated by a conservative upwind finite-difference scheme. We establish several structural properties of the implicit scheme in arbitrary space dimension, including a discrete maximum principle, non-negativity, mass conservation, ℓ∞ and ℓ2 stability, existence of non-negative discrete solutions, and a discrete energy-dissipation estimate. Our assumptions include natural symmetry, a sign condition on the sum of the directional coefficients, and decay on the discrete weights. We also derive a consistency estimate for the complete nonlinear spatial operator. For rigorous convergence analysis, we restrict ourselves to one space dimension and introduce the cumulative variable associated with the density. This transforms the discrete problem into a monotone approximation of an equation involving the fractional Laplacian. Interestingly, for semigroup-based pressure-gradient weights, the resulting discrete operator is identified with the fractional discrete Laplacian. Next, using comparison and contraction properties of the integrated scheme, for m≥2 and under the stated spatial regularity and weighted temporal consistency assumptions, we prove uniform convergence of the cumulative numerical solution and convergence of the corresponding density against compactly supported smooth test functions. The temporal estimates allow a weak singularity at the initial time. We also analyze an explicit variant of the scheme and derive a sufficient CFL condition covering the full range 0<σ<1.
Related
- Structure-Preserving Numerical Schemes for Two-Stage Local and Nonlocal Dispersal Systems
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
- A structure-preserving finite volume scheme for the time-fractional Allen-Cahn equation with nondegenerate state-dependent mobility
- On a Nonlinear Fractional Delay Differential Equation with Nonlocal Boundary Conditions
- Numerical Simulation of the Time-Fractional Nonlinear Schrödinger Equation and the Time-Fractional Schrödinger–Poisson System
BibTeX
@misc{plociniczak2026structurepreserving,
title = {{Structure-Preserving Numerical Schemes for a Time-Fractional Porous Medium Equation with Nonlocal Pressure}},
author = {Łukasz Płociniczak and Hubert Woszczek},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.27172},
eprint = {2610.27172},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.27172},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Structure-Preserving Numerical Schemes for a Time-Fractional Porous Medium Equation with Nonlocal Pressure
AU - Łukasz Płociniczak
AU - Hubert Woszczek
PY - 2026
DA - 2026/10/24
PB - arXiv
DO - 10.48550/arxiv.2610.27172
UR - https://arxiv.org/abs/2610.27172
AB - In this work, we study numerical approximations of the equation ∂_t^αu−∇·(u^m−1∇(−Δ)^−σu)=0, where 0<α<1, 0<σ<1 and m>1. The Caputo derivative is discretized by the L1 formula, while the nonlinear nonlocal flux is approximated by a conservative upwind finite-difference scheme. We establish several structural properties of the implicit scheme in arbitrary space dimension, including a discrete maximum principle, non-negativity, mass conservation, ℓ^∞ and ℓ^2 stability, existence of non-negative discrete solutions, and a discrete energy-dissipation estimate. Our assumptions include natural symmetry, a sign condition on the sum of the directional coefficients, and decay on the discrete weights. We also derive a consistency estimate for the complete nonlinear spatial operator. For rigorous convergence analysis, we restrict ourselves to one space dimension and introduce the cumulative variable associated with the density. This transforms the discrete problem into a monotone approximation of an equation involving the fractional Laplacian. Interestingly, for semigroup-based pressure-gradient weights, the resulting discrete operator is identified with the fractional discrete Laplacian. Next, using comparison and contraction properties of the integrated scheme, for m≥2 and under the stated spatial regularity and weighted temporal consistency assumptions, we prove uniform convergence of the cumulative numerical solution and convergence of the corresponding density against compactly supported smooth test functions. The temporal estimates allow a weak singularity at the initial time. We also analyze an explicit variant of the scheme and derive a sufficient CFL condition covering the full range 0<σ<1.
ER -