Numerical investigation of two models of nonlinear fractional reaction subdiffusion equations

Abstract We consider new numerical schemes to solve two different systems of nonlinear fractional reaction subdiffusion equations. These equations model the reversible $$A+B \rightleftharpoons C$$ A + B ⇌ C in presence anomalous subdiffusion. The first is based on Henry & Wearne [1] where term added equation. second by Angstmann, Donnelly [2] which involves a modified differential operator. For both models Keller Box method [3] along with L1 scheme (ML1), adapted from Oldham and Spanier [4], are used approximate spatial derivatives respectively. Numerical prediction were compared for number examples given same initial boundary conditions exponents. From results, we see similar short time behaviour predicted. However long times solution remains positive whilst based–model predictions may become negative.