Asymptotically Compatible Variational Energy Dissipation Laws and Structure-Preserving Variable-Step Linear Relaxation Schemes for a Subdiffusive Tumor-Nutrient Model
Abstract
In this work, we develop a mass-conservative subdiffusive tumor-nutrient model coupled through chemotaxis and memory-dependent reaction exchange, together with structure-preserving variable-step linear relaxation schemes. A weighted Riemann-Liouville identity is proved and used to construct a variational energy and the associated energy dissipation law for the proposed model. Both the model and the associated variational energy dissipation law are asymptotically compatible with their integer-order counterparts as the fractional order approaches one. For numerical approximation, we discretize the Riemann-Liouville fractional derivatives using the L1R formula on nonuniform time meshes and derive the associated discrete weighted inequality. We first construct a variable-step linear relaxation scheme for the subdiffusive model. Then, to mitigate alternating perturbations induced by the auxiliary-variable error and amplified by the L1R convolution at small fractional orders, we construct a temporally damped linear relaxation scheme. Since temporal damping leaves the auxiliary-variable recurrence unchanged, we further introduce an a posteriori relaxation correction to reduce alternating errors in the auxiliary variable. All three schemes preserve total mass, are proven to satisfy discrete variational energy dissipation laws using the discrete weighted inequality, and remain consistent with the integer-order limit. Several numerical examples are presented to validate the accuracy and stability of the proposed schemes.
Related
- A nonlocal nonlinear Schrödinger model: well-posedness, local limit, and structure-preserving asymptotically compatible Fourier approximations
- Dissipation Bounds and Hardness in Globally Asymptotically Stable Quadratic Systems
- Fast and scalable multi-GPU sliced FFT solvers for asymptotically compatible discretizations of 3D nonlocal diffusion problems
- A Structure-Preserving and Thermodynamically Compatible Cell-Centered Lagrangian Finite Volume Scheme for Continuum Mechanics
- Asymptotically Valid Laplace Variational Inference for Bayesian Nonparametric Clustering of Marked Point Patterns
BibTeX
@misc{yu2026asymptotically,
title = {{Asymptotically Compatible Variational Energy Dissipation Laws and Structure-Preserving Variable-Step Linear Relaxation Schemes for a Subdiffusive Tumor-Nutrient Model}},
author = {Hui Yu and Zhaoyang Wang and Ping Lin},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.40120},
eprint = {2610.40120},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.40120},
}
RIS
TY - UNPB
TI - Asymptotically Compatible Variational Energy Dissipation Laws and Structure-Preserving Variable-Step Linear Relaxation Schemes for a Subdiffusive Tumor-Nutrient Model
AU - Hui Yu
AU - Zhaoyang Wang
AU - Ping Lin
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.40120
UR - https://arxiv.org/abs/2610.40120
AB - In this work, we develop a mass-conservative subdiffusive tumor-nutrient model coupled through chemotaxis and memory-dependent reaction exchange, together with structure-preserving variable-step linear relaxation schemes. A weighted Riemann-Liouville identity is proved and used to construct a variational energy and the associated energy dissipation law for the proposed model. Both the model and the associated variational energy dissipation law are asymptotically compatible with their integer-order counterparts as the fractional order approaches one. For numerical approximation, we discretize the Riemann-Liouville fractional derivatives using the L1R formula on nonuniform time meshes and derive the associated discrete weighted inequality. We first construct a variable-step linear relaxation scheme for the subdiffusive model. Then, to mitigate alternating perturbations induced by the auxiliary-variable error and amplified by the L1R convolution at small fractional orders, we construct a temporally damped linear relaxation scheme. Since temporal damping leaves the auxiliary-variable recurrence unchanged, we further introduce an a posteriori relaxation correction to reduce alternating errors in the auxiliary variable. All three schemes preserve total mass, are proven to satisfy discrete variational energy dissipation laws using the discrete weighted inequality, and remain consistent with the integer-order limit. Several numerical examples are presented to validate the accuracy and stability of the proposed schemes.
ER -