Abstract:We establish a uniform-in-scaling error estimate for the asymptotic preserving scheme proposed in \cite{XW21} for the Lévy-Fokker-Planck (LFP) equation. The main difficulties stem from not only the interplay between the scaling and numerical parameters but also the slow decay of the tail of the equilibrium state. We tackle these problems by separating the parameter domain according to the relative size of the scaling $\epsilon$: in the regime where $\epsilon$ is large, we design a weighted norm to mitigate the issue caused by the fat tail, while in the regime where $\epsilon$ is small, we prove a strong convergence of LFP towards its fractional diffusion limit with an explicit convergence rate. This method extends the traditional AP estimates to cases where uniform bounds are unavailable. Our result applies to any dimension and to the whole span of the fractional power.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to conduct a unified error estimate on the Asymptotic Preserving (AP) numerical scheme of the Lévy - Fokker - Planck (LFP) equation. Specifically, the authors hope to establish an error estimate for the AP scheme under different scale parameters \(\varepsilon\), ensuring that this estimate is independent of \(\varepsilon\) and the time step \(\Delta t\).
### Main Challenges
1. **Interaction between scale parameters and numerical parameters**: When the scale parameter \(\varepsilon\) changes, how to ensure that the error estimate of the numerical solution remains valid is a key issue.
2. **Slow decay of the equilibrium tail**: The equilibrium state of the LFP equation has an algebraically decaying tail, which makes traditional error estimation methods no longer applicable.
### Solutions
To address the above challenges, the authors adopted the following strategies:
1. **Application of weighted norms**: In the case of a large \(\varepsilon\), a weighted norm was designed to alleviate the problems caused by the slow decay of the equilibrium tail.
2. **Proof of strong convergence**: In the case of a small \(\varepsilon\), the strong convergence of the LFP equation to its fractional - order diffusion limit was proven, and the explicit convergence rate was given.
3. **Partition of the parameter domain**: The parameter domain was divided into two regions (Regime I and Regime II) according to the size of \(\varepsilon\), and they were processed separately in each region.
### Specific Contributions
- **Unified error estimate**: A unified error estimate applicable to all scale parameters \(\varepsilon\) was established, expanding the application range of traditional AP estimates.
- **Applicability to multiple dimensions and fractional powers**: The results are applicable to any dimension and the entire range of fractional powers.
### Conclusions
Through these methods, the authors successfully solved the problem of unified error estimation for the AP scheme in the LFP equation, providing a solid theoretical foundation for numerical simulation. This achievement not only represents a theoretical breakthrough but also provides reliable guidance for numerical calculations in practical applications.
### Formula Summary
The main formulas involved in the paper include:
- Lévy - Fokker - Planck equation:
\[
\partial_t f + v \cdot \nabla_x f=\nabla_v \cdot (vf)-(-\Delta_v)^s f := L_s(f), \quad s \in (0,1)
\]
- Definition of the fractional - order Laplacian operator:
\[
(-\Delta_v)^s \phi(v)=\mathcal{F}^{-1}(|k|^{2s} \mathcal{F}[\phi](k))
\]
- Properties of the equilibrium state \(M(v)\):
\[
L_s(M) = 0, \quad \int_{\mathbb{R}^d} M(v) dv = 1, \quad M(v)\sim\frac{C}{|v|^{d + 2s}}\text{ as }|v|\to\infty
\]
These formulas help us understand the structure of the LFP equation and the characteristics of its solutions, thereby providing a theoretical basis for error estimation.