A new class of semi-implicit methods with linear complexity for nonlinear fractional differential equations
This work provides efficient and stable numerical methods for solving nonlinear fractional differential equations, which are important in modeling anomalous diffusion and other phenomena.
The authors propose a new class of semi-implicit methods for nonlinear fractional differential equations that are unconditionally stable with suitable parameters, and develop a fast convolution strategy reducing memory and computational cost to O(N log n_T) and O(N n_T log n_T), respectively. Numerical simulations verify effectiveness.
We propose a new class of semi-implicit methods for solving nonlinear fractional differential equations and study their stability. Several versions of our new schemes are proved to be unconditionally stable by choosing suitable parameters. Subsequently, we develop an efficient strategy to calculate the discrete convolution for the approximation of the fractional operator in the semi-implicit method and we derive an error bound of the fast convolution. The memory requirement and computational cost of the present semi-implicit methods with a fast convolution are about $O(N\log n_T)$ and $O(Nn_T\log n_T)$, respectively, where $N$ is a suitable positive integer and $n_T$ is the final number of time steps. Numerical simulations, including the solution of a system of two nonlinear fractional diffusion equations with different fractional orders in two-dimensions, are presented to verify the effectiveness of the semi-implicit methods.