Local modification of subdiffusion by initial Fickian diffusion: Multiscale modeling, analysis and computation

Xiangcheng Zheng,Yiqun Li,Wenlin Qiu

2024-01-30

Abstract:We propose a local modification of the standard subdiffusion model by introducing the initial Fickian diffusion, which results in a multiscale diffusion model. The developed model resolves the incompatibility between the nonlocal operators in subdiffusion and the local initial conditions and thus eliminates the initial singularity of the solutions of the subdiffusion, while retaining its heavy tail behavior away from the initial time. The well-posedness of the model and high-order regularity estimates of its solutions are analyzed by resolvent estimates, based on which the numerical discretization and analysis are performed. Numerical experiments are carried out to substantiate the theoretical findings.

Numerical Analysis

What problem does this paper attempt to address?

The paper aims to address the issue of incompatibility between nonlocal operators in subdiffusion models and local initial conditions. This incompatibility leads to initial singularities in the solutions of such models. To resolve this issue, the authors propose a modification to the standard subdiffusion model by incorporating an initial Fickian diffusion term. This results in a multiscale diffusion model that combines Fickian diffusion near the initial time and subdiffusion behavior away from the initial time.

The main contributions of the paper are summarized as follows:

1. **Model Proposal**: The authors propose a local modification of the standard subdiffusion model by introducing an initial Fickian diffusion term. This leads to a multiscale diffusion model that resolves the incompatibility between nonlocal operators in subdiffusion and local initial conditions, thereby eliminating the initial singularity of the solutions of the subdiffusion model.

2. **Mathematical Analysis**: The paper analyzes the well-posedness and regularity of the proposed multiscale diffusion model. It proves that the model is well-posed and shows that the solutions have higher-order regularity estimates. This analysis is based on resolvent estimates.

3. **Numerical Discretization and Computation**: The paper develops a fully discrete finite element scheme for the proposed multiscale diffusion model. It also provides error estimates for the numerical scheme, showing that the numerical solution