229
Views
11
CrossRef citations to date
0
Altmetric
Original Articles

A fast high-order compact difference method for the fractal mobile/immobile transport equation

, &
Pages 1860-1883 | Received 16 Oct 2017, Accepted 12 Jul 2019, Published online: 20 Sep 2019
 

ABSTRACT

In this paper, we present a high-order compact finite difference scheme for fractal mobile/immobile transport model with a Caputo fractional time derivative. The compact finite difference scheme is stable and convergent with convergence order O(τ2α+h4) in L2-norm. Furthermore, because of the non-local property of fractional differential operators, it is necessary to find a fast technique to reduce the computational cost. We develop a fast solution technique which is based on a fast Fourier transform. The computational work will reduce to O(MNlog2N) while the direct method requires the computational cost of O(MN2), where N=τ1 and τ is the size of time step, M=h1 and h is the size of space step. Moreover, numerical results are consistent with the theoretical analysis.

2010 AMS Subject Classifications:

Acknowledgments

We would like to acknowledge the assistance of volunteers in putting together this example manuscript and supplement.

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

This work is supported by the Postdoctoral Science Foundation of China [grant numbers BX20190187 and 2019M650152], by National Natural Science Foundation of China [grant numbers 11931003, 41974133, 11901489, 11971276]; National Science Foundation DMS-1216923; OSD/ARO MURI W911NF-15-1-0562; Shandong Provincial Natural Science Foundation, China ZR2011AM015; National Science and Technology Major Project of China 2011ZX05052, 2011ZX05011-004.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.