Abstract:An adaptive finite difference scheme for variable-order fractional-time subdiffusion equations in the Caputo form is studied. The fractional time derivative is discretized by the L1 procedure but using nonhomogeneous timesteps. The size of these timesteps is chosen by an adaptive algorithm in order to keep the local error bounded around a preset value, a value that can be chosen at will. For some types of problems this adaptive method is much faster than the corresponding usual method with fixed timesteps while keeping the local error of the numerical solution around the preset values. These findings turns out to be similar to those found for constant-order fractional diffusion equations.
What problem does this paper attempt to address?
This paper attempts to solve the numerical solution problem of variable - order fractional - time subdiffusion equations. Specifically, the author has studied an adaptive finite - difference method for solving variable - order fractional - time subdiffusion equations in Caputo form. The following are the main problems that the paper attempts to solve:
1. **Improve the efficiency of numerical solution**: Traditional fixed - step - size methods have low computational efficiency when dealing with fractional - order differential equations, especially in the case of anomalous diffusion. To improve computational efficiency, the author proposes an adaptive algorithm that can dynamically adjust the time step according to the behavior of the solution.
2. **Keep the local error within the preset range**: By using non - uniform time steps and an adaptive algorithm, ensure that the local error of the numerical solution is controlled around a preset value. This method can significantly reduce the computation time while ensuring accuracy.
3. **Handle different types of variable - order fractional - order diffusion equations**: The paper not only discusses the general form of variable - order fractional - order diffusion equations, but also pays special attention to the case where the diffusion order γ depends on time and space variables. This makes the method applicable to a wider range of physical and engineering problems.
4. **Verify the effectiveness of the adaptive method**: Through several specific numerical experiments, the effectiveness and superiority of the adaptive method under different conditions are verified, especially in long - time simulations and the rapidly changing initial stage.
### Formula Summary
The key formulas involved in the paper include:
- **Discretization of Caputo derivative**:
\[
\frac{\partial^\gamma u(x_j, t_n)}{\partial t^\gamma} = \frac{1}{\Gamma(1-\gamma_j^n)} \sum_{m = 0}^{n - 1} \frac{u_j^{m + 1}-u_j^m}{\tau_{m + 1}} \int_{t_m}^{t_{m + 1}} \frac{dt'}{(t - t')^{\gamma_j^n}}+R_j(t_n)
\]
where,
\[
T_{j,m,n}=\frac{(t_n - t_m)^{1-\gamma_j^n}-(t_n - t_{m + 1})^{1-\gamma_j^n}}{\tau_{m + 1}}, \quad m\leq n - 1
\]
- **Spatial discretization**:
\[
\frac{\partial^2 u(x, t)}{\partial x^2}\bigg|_{(x_j, t_n)}=\frac{u_j^{n + 1}-2u_j^n + u_j^{n - 1}}{(\Delta x)^2}+R_n(x_j)
\]
- **Final finite - difference scheme**:
\[
- S_j^n U_j^{n + 1}+(1 + 2S_j^n)U_j^n - S_j^n U_j^{n - 1}=M[U_j^{(n)}]+\tilde{F}_j^n
\]
where,
\[
S_j^n=\frac{\Gamma(2-\gamma_j^n)K(t_n - t_{n - 1})^{\gamma_j^n}}{(\Delta x)^2}
\]
\[
M[U_j^{(n)}]=U_j^{n - 1}-\sum_{m = 0}^{n - 2}\tilde{T}_{j,m,n}[U_j^{m + 1}-U_j^m]
\]
\[
\tilde{T}_{j,m,n}=(t_n - t_{n - 1})^{\gamma_j^n}T_{j,m,n}, \quad \tilde{F}_j^n=(t_n - t_{n - 1})^{\gamma_j^n}F(x_j, t_n)
\]
Through these