An efficient implicit FEM scheme for fractional-in-space reaction-diffusion equations

SIAM Journal on Scientific Computing

Volumn 34, Issue 4, 2012, Pages

  Burrage, Kevin ^a,b   Hale, Nicholas ^a   Kay, David ^a  

^a UNIVERSITY OF OXFORD   (United Kingdom)

^b QUEENSLAND UNIVERSITY OF TECHNOLOGY   (Australia)

Author keywords

Finite elements; Fractional diffusion; Numerical solvers

Indexed keywords

DIFFUSION IN LIQUIDS; LINEAR EQUATIONS; MATHEMATICAL OPERATORS; PARTIAL DIFFERENTIAL EQUATIONS; VECTOR SPACES;

ANOMALOUS DIFFUSION; FINITE ELEMENT; FRACTIONAL DIFFERENTIAL EQUATIONS; FRACTIONAL DIFFERENTIAL OPERATORS; FRACTIONAL DIFFUSION; FRACTIONAL POWER; MODELLING TOOLS; NUMERICAL SOLVERS; REACTION DIFFUSION EQUATIONS; SPATIAL HETEROGENEITY;

FINITE ELEMENT METHOD;

EID: 84866419689     PISSN: 10648275     EISSN: 10957200     Source Type: Journal    
DOI: 10.1137/110847007     Document Type: Article

References (71)

1. E. E. ADAMS AND L. GELHAR, Field study of dispersion in a heterogeneous aquifer: 2. Spatial moment analysis, Water Resources Res., 28 (1992), pp. 3293-3307.
2. N. I. AKHIEZER, Elements of the Theory of Elliptic Functions, Transl. Math. Monogr. 79, AMS, Providence, RI, 1990.
3. F. ANDREU, J. M. MAZÓN, J. D. ROSSI, AND J. TOLEDO, A nonlocal p-Laplacian evolution equation with Neumann boundary conditions, J. Math. Pures Appl., 90 (2008), pp. 201-227.
4. F. ANDREU, J. M. MAZÓN, J. D. ROSSI, AND J. TOLEDO, A nonlocal p-Laplacian evolution equation with nonhomogeneous Dirichlet boundary conditions, SIAM J. Math. Anal., 40 (2009), pp. 1815-1851.
5. A. BAYLISS, C. I. GOLDSTEIN, AND E. TURKEL, An iterative method for the Helmholtz equation, J. Comput. Phys., 49 (1983), pp. 443-457.
6. P. BECKER-KERN, M. M. MEERSCHAERT, AND H. P. SCHEFFLER, Limit theorem for continuous time random walks with two time scales, J. Appl. Probab., 41 (2004), pp. 455-466.
7. B. BECKERMANN AND S. GÜTTEL, Superlinear convergence of the rational Arnoldi method, Numer. Math., 121 (2012), pp. 205-236.
8. B. BECKERMANN AND L. REICHEL, Error estimates and evaluation of matrix functions via the Fa bertrans form, SIAM J. Numer. Anal., 47 (2009), pp. 3849-3883.
9. D. A. BENSON, S. W. WHEATCRAFT, AND M. M. MEERSCHAERT, Application of a fractional advection-dispersion equation, Water Resources Res., 36 (2000), pp. 1403-1412.
10. J. P. BOUCHAUD AND A. GEORGES, Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications, Phys. Rep., 195 (1990), pp. 127-293.
11. N. BURCH AND R. LEHOUCQ, Classical, nonlocal, and fractional diffusion equations on bounded domains, Internat. J. Multiscale Comp. Eng., to appear.
12. C. CHEN, F. LIU, AND K. BURRAGE, Finite difference methods and a Fourier analysis for the fractional reaction-subdiffusion equation, Appl. Math. Comput., 198 (2008), pp. 754-769.
13. C. CHEN, F. LIU, I. TURNER, AND V. ANH, Fourier method for the fractional diffusion equation describing sub-diffusion, J. Comput. Phys., 227 (2007), pp. 886-897.
14. T. A. DRISCOLL, The Schwarz-Christoffel toolbox, available online at http://www.math.udel.edu/~driscoll/software/SC/.
15. V. DRUSKIN AND L. KNIZHNERMAN, Extended Krylov subspaces: Approximation of the matrix square root and related functions, SIAM J. Matrix Anal. Appl., 19 (1998), pp. 755-771.
16. Q. DU AND K. ZHOU, Mathematical and numerical analysis of linear peridynamic models with nonlocal boundary conditions, SIAM J. Numer. Anal., 48 (2010), pp. 1759-1780.
17. H. ENGLER, On the speed of spread for fractional reaction-diffusion equations, Int. J. Difference Equ. (2010), 31542.
18. Y. A. ERLANGGA, C. W. OOSTERLEE, AND C. VUIK, Comparison of multigrid and incomplete LU shifted-Laplace preconditioners for the inhomogeneous Helmholtz equation, Appl. Numer. Math., 56 (2006), pp. 648-666.
19. Y. A. ERLANGGA, C. VUIK, AND C. W. OOSTERLEE, On a class of preconditioners for the Helmholtz equation, Appl. Numer. Math., 20 (2005), pp. 409-425.
20. M. B. VAN GIJZEN, Y. A. ERLANGGA, AND C. VUIK, Spectral analysis of the discrete Helmholtz operator preconditioned with a shifted Laplacian, SIAM J. Sci. Comput., 29 (2007), pp. 1942-1958.
21. Q. Y. GUAN AND Z. M. MA, Boundary problems for fractional Laplacians, Stoch. Dyn., 5 (2005), pp. 385-424.
22. M. GUNZBURGER AND R. B. LEHOUCQ, A nonlocal vector calculus with application to nonlocal boundary value problems, Multiscale Model. Simul., 8 (2010), pp. 1581-1598.
23. S. GÜTTEL, Rational Krylov Methods for Operator Functions, Ph.D. thesis, TU Bergakademie Freiberg, Freiberg, Germany, 2010.
24. α,log(A), and related matrix functions by contour integrals, SIAM J. Numer. Anal., 46 (2008), pp. 2505-2523.
25. N. J. HIGHAM, Functions of Matrices: Theory and Computation, SIAM, Philadelphia, 2008.
26. U. HORNUNG, Homogenization and Porous Media, Springer-Verlag, New York, 1997.
27. M. ILIĆ, F. LIU, I. TURNER, AND V. ANH, Numerical approximation of a fractional-in-space diffusion equation, I, Frac. Calc. and App. Anal., 8 (2005), pp. 323-341.
28. M. ILIĆ, F. LIU, I. TURNER, AND V. ANH, Numerical approximation of a fractional-in-space diffusion equation (II) - with nonhomogeneous boundary conditions, Fract. Calc. Appl. Anal., 9 (2006), pp. 333-349.
29. M. ILIĆ, I. TURNER, F. LIU, AND V. ANH, Analytical and numerical solutions of a one-dimensional fractional-in-space diffusion equation in a composite medium, Appl. Math. Comput., 216 (2010), pp. 2248-2262.
30. M. ILIĆ AND I. W. TURNER, Approximating functions of a large sparse positive definite matrix using a spectral splitting method, in Proceedings of the 12th Computational Techniques and Applications Conference CTAC-2004, 2005, pp. C472-C487.
31. M. ILIĆ, I. W. TURNER, AND V. ANH, A numerical solution using an adaptively preconditioned Lanczos method for a class of linear systems related with the fractional Poisson equation, J. Appl. Math. Stoch. Anal. (2008), 104525.
32. D. KAY AND R. WELFORD, A multigrid finite element solver for the Cahn-Hilliard equation, J. Comput. Phys., 212 (2006), pp. 288-304.
33. L. KNIZHNERMAN AND V. SIMONCINI, A new investigation of the extended Krylov subspace method for matrix function evaluations, Linear Algebra Appl., 17 (2010), pp. 615-638.
34. A. L. LAIRD AND M. B. GILES, Preconditioning Harmonic Unsteady Potential Flow Calculations, Technical report 02/12, Mathematical Institute, University of Oxford, Oxford, UK, 2002.
35. T. A. M. LANGLANDS AND B. HENRY, The accuracy and stability of an implicit solution method for the fractional diffusion equation, J. Comput. Phys., 205 (2005), pp. 719-736.
36. A. LEVIN AND E. SAFF, Optimal ray sequences of rational functions connected with the Zolotarev problem, Constr. Approx., 10 (1994), pp. 235-273.
37. X. LI AND C. XU, Existence and uniqueness of the weak solution of the space-time fractional diffusion equation and a spectral method approximation, Commun. Comput. Phys., 8 (2010), pp. 1016-1051.
38. Y. LIN AND C. XUE, Finite difference/spectral approximations for the time-fractional diffusion equation, J. Comput. Phys., 225 (2007), pp. 15333-1552.
39. F. LIU, V. ANH, AND I. TURNER, Numerical solution of the space fractional Fokker-Planck equation, J. Comput. Appl. Math., 166 (2004), pp. 209-219.
40. F. LIU AND K. BURRAGE, Parameter estimation for fractional dynamical models in biological systems, in Proceedings of the Fourth IFAC Workshop on Fractional Differentiation and Its Applications, Badajoz, 2010.
41. F. LIU, I. TURNER, AND M. COX, A finite volume method for saltwater intrusion with heavy-tailed motions, in review.
42. F. LIU, P. ZHUANG, V. ANH, I. TURNER, AND K. BURRAGE, Stability and convergence of the difference methods for the space-time fractional advection-diffusion equation, Appl. Math. Comput., 191 (2007), pp. 12-20.
43. R. L. MAGIN, O. ABDULLAH, D. BALEANU, AND X. J. ZHOU, Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation, J. Mag. Res., 190 (2008), pp. 255-270.
44. M. M. MEERSCHAERT, D. A. BENSON, AND S. W. WHEATCRAFT, Subordinated advection-dispersion equation for contaminant transport, Water Resource Res., 37 (2001), pp. 1543-1550.
45. M. M. MEERSCHAERT, J. MORTENSEN, AND S. W. WHEATCRAFT, Fractional vector calculus for fractional advection-dispersion, Phys. A, 367 (2006), pp. 181-190.
46. M. M. MEERSCHAERT AND C. TADJERAN, Finite difference approximations for two-sided space-fractional partial differential equations, Appl. Numer. Math, 56 (2006), pp. 80-90.
47. R. METZLER AND J. KLAFTER, The random walk's guide to anomalous diffusion: A fractional dynamics approach, Phys. Rep., 339 (2000), pp. 1-77.
48. L. REICHEL, Newton interpolation at Leja points, BIT, 30 (1990), pp. 332-346.
49. M. RIESZ, L'intégrale de Rieman-Liouville et le problème de Cauchy, Acta Math., 81 (1949), pp. 1-223.
50. 2, J. Comput. Appl. Math., 193 (2005), pp. 243-268.
51. A. RUHE, Rational Krylov sequences for eigenvalue computations, Linear Algebra Appl, 58 (1984), pp. 391-405.
52. A. RUHE, Rational Krylov: A practical algorithm for large sparse nonsymmetric matrix pencils, SIAM J. Sci. Comput., 19 (1998), pp. 1535-1551.
53. Y. SAAD, Preconditioning techniques for non symmetric and indefinite linear system, J. Comput. Appl. Math., 24 (1988), pp. 89-105.
54. M. SAXTON, Anomalous subdiffusion in fluorescence photobleaching recovery: A Monte Carlo study, J. Biophys., 81 (2001), pp. 2226-2240.
55. E. SCALAS, R. GORENFLO, AND F. MAINARDI, Fractional calculus and continuous time finance, Phys. A, 284 (2000), pp. 376-384.
56. R. SCHUMER, D. A. BENSON, M. M. MEERSCHAERT, AND B. BAEUMER, Fractal mobile/immobile solute transport, Water Res., 39 (2003), pp. 1296-1307.
57. V. SIMONCINI, A new iterative method for solving large-scale Lyapunov matrix equations, SIAM J. Sci. Comput., 29 (2007), pp. 1268-1288.
58. D. P. SIMPSON, Krylov Subspace Methods for Approximating Functions of Symmetric Positive Definite Matrices with Applications to Applied Statistics and Anomalous Diffusion, Ph.D. thesis, Queensland University of Technology, Brisbane, 2008.
59. Z. SUN AND X. WU, A fully discrete difference scheme for a diffusion-wave system, Appl. Numer. Math., 56 (2006), pp. 193-209.
60. I. TURNER, M. ILIĆ, AND P. PERR, The use of fractional-in-space diffusion equations for describing microscale diffusion in porous media, in Proceedings of the 11th International Drying Conference, Magdeburg, Germany, 2010.
61. H. A. VAN DER VORST, An iterative solution method for solving f(A)x = b, using Krylov subspace information obtained for the symmetric positive definite matrix A, J. Comput. Appl. Math., 18 (1987), pp. 249-263.
62. H. A. VAN DER VORST, Bi-CGSTAB: A fast and smoothly converging variant of Bi-CG for solution of nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 13 (1992), pp. 631-644.
63. A. J. WATHEN, Realistic eigenvalue bounds for the Galerkin mass matrix, IMA J. Numer. Anal., 7 (1987), pp. 449-457.
64. S. WHITAKER, Flow in porous media I: A theoretical derivation of Darcy's law, Transp. Porous Media, 1 (1986), pp. 3-25.
65. Q. YANG, Novel Analytical and Numerical Methods for Solving Fractional Dynamical Systems, Ph.D. thesis, Queensland University of Technology, Brisbane, 2005.
66. Q. YANG, I. TURNER, F. LIU, AND M. ILIĆ, Novel numerical methods for solving the timespace fractional diffusion equation in two dimensions, SIAM J. Sci. Comput., 33 (2011), pp. 1159-1180.
67. Q. YANG, F. LIU, AND I. TURNER, Numerical methods for fractional partial differential equations with Riesz space fractional derivatives, Appl. Math. Model., 34 (2010), pp. 200-218.
68. S. B. YUSTE AND L. ACEDO, An explicit finite difference method and a new von Neumanntype stability analysis for fractional diffusion equations, SIAM J. Numer. Anal., 42 (2005), pp. 1862-1874.
69. S. B. YUSTE, L. ACEDO, AND K. LINDENBERG, Reaction front in an A + B - C reaction-subdiffusion process, Phys. Rev. E, 69 (2004), 036126.
70. Y. ZHANG, D. A. BENSON, AND D. M. REEVES, Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications, Adv. Water Res., 32 (2009), pp. 561-581.
71. P. ZHUANG, F. LIU, V. ANH, AND I. TURNER, Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term, SIAM J. Numer. Anal., 47 (2009), pp. 1760-1781.