Structure-preserving scheme for fractional nonlinear diffusion equations

Hélène Hivert,Florian Salin
2024-09-27
Abstract:In this paper, we introduce and analyze a numerical scheme for solving the Cauchy-Dirichlet problem associated with fractional nonlinear diffusion equations. These equations generalize the porous medium equation and the fast diffusion equation by incorporating a fractional diffusion term. We provide a rigorous analysis showing that the discretization preserves main properties of the continuous equations, including algebraic decay in the fractional porous medium case and the extinction phenomenon in the fractional fast diffusion case. The study is supported by extensive numerical simulations. In addition, we propose a novel method for accurately computing the extinction time for the fractional fast diffusion equation and illustrate numerically the convergence of rescaled solutions towards asymptotic profiles near the extinction time.
Numerical Analysis
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to develop and analyze numerical solution schemes for the Cauchy - Dirichlet problem of fractional nonlinear diffusion equations (FNDEs). Specifically, the paper addresses the following problems: 1. **Generalize classical models**: - The paper generalizes the classical porous medium equation (PME) and fast diffusion equation (FDE) to forms that include fractional - order diffusion terms. These equations can describe more complex diffusion phenomena, especially in systems with long - range interactions or heterogeneity. 2. **Maintain the main properties of continuous equations**: - The researchers propose a discretization scheme to ensure that the numerical method can maintain the main characteristics of the continuous equations. For example, in the fractional porous medium equation (FPME), the solution decays at an algebraic rate; while in the fractional fast diffusion equation (FFDE), the solution disappears within a finite time (extinction phenomenon). 3. **Accurately calculate the extinction time**: - For FFDE, the paper proposes a novel method to accurately calculate the time when the solution disappears (extinction time), which is difficult to achieve with traditional methods. This not only helps to understand the behavior of the solution but also can be used to study the long - time asymptotic behavior. 4. **Verify the results of numerical simulations**: - Through extensive numerical simulations, the researchers verify the effectiveness and accuracy of the proposed discretization scheme. In particular, they show that the rescaled solution converges to the asymptotic profile as it approaches the extinction time and determine the convergence rate. 5. **Combine theory with practice**: - The research does not stop at the theoretical level but also verifies the theoretical results through numerical experiments, making the research results more practical. This is very important for practical application fields that need to handle complex diffusion processes, such as materials science and environmental science. In summary, the core problem of this paper is to develop a numerical method that can effectively handle fractional nonlinear diffusion equations and ensure that this method can maintain the key characteristics of the original equations, thereby providing a powerful tool for scientific research in related fields.