Stability and convergence of multistep schemes for 1D and 2D fractional model with nonlinear source term

Stable multistep schemes based on Caputo fractional derivative approximation are pre-sented for solving 1D and 2D nonlinear fractional model arising from dielectric media. We approximate Caputo fractional derivatives in time with a multistep scheme of order O(tau(3-alpha)) & O(tau(3-beta)), 1 < beta < 2 , spatial Laplacian operator with a central difference scheme, and nonlinear source term g(B) by using Taylor series. The discretization of the problem results in a linear system of equations that is tridiagonal and penta-diagonal for 1D and 2D case, respectively. The unique solvability and unconditional stability are derived for both cases. The convergence of schemes is established with the help of optimal error bounds. Further, we establish that the order of convergence for 1D case is O(tau(3-alpha) + tau(3-beta) + h(2)) and for 2D case is O(tau(3-alpha) + tau(3-beta) + h(x)(2) + h(y)(2)). Moreover, the stability of our schemes are verified numerically by adding some linear and nonlinear noisy inputs. Finally, four test functions are investigated to show the effectiveness and stability of our schemes. The method is simple, easy to implement, and yields very accurate results. (C) 2020 Elsevier Inc. All rights reserved.

Automated summary

Stable multistep schemes based on Caputo fractional derivative approximation are proposed for solving 1D and 2D nonlinear fractional model arising from dielectric media. The schemes are numerically verified to be effective and stable through test functions.