On the numerical solutions of coupled nonlinear time-fractional reaction-diffusion equations

In this paper, we investigate the numerical solutions of coupled nonlinear time-fractional reaction-diffusion equations obtained by applying a procedure that combines the Lie symmetry analysis with the numerical methods. By Lie symmetries, the model, governed by two fractional differential equations defined in terms of the Riemann-Liouville fractional derivatives, is reduced into nonlinear fractional ordinary differential equations that, by introducing the Caputo derivative, are numerically solved by the implicit trapezoidal method. The solutions of the original model are computed by the inverse transformations. Numerical examples are performed in order to show the efficiency and the reliability of the proposed approach applied for solving a wide class of fractional models.