Abstract

A new numerical method that solves a time-space-fractional reaction-diffusion equation effectively is presented. The space-fractional Riesz operator is first discretized using fourth order compact finite differences, resulting in a system with a linear stiff term. The system of differential equations is then solved through time integration using an exponential time difference method, which explicitly handles the nonlinear non-stiff term. Restricted Padé rational approximation of a single real pole is used to approximate the matrix exponential e-A. The problems concerning the stability and computational efficiency of this new approach are tackled by means of a splitting technique. Based on compact finite differences with third-order accuracy in time and restricted Padé approximation, this novel efficient technique reduces computing cost and effectively tackles the problems with non-smooth initial data. The superiority of new method in terms of accuracy, reliability, and computational efficiency is finally shown by numerical examples.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.