1,387 research outputs found
An efficient implementation of an implicit FEM scheme for fractional-in-space reaction-diffusion equations
Fractional differential equations are becoming increasingly used as a modelling tool for processes with anomalous diffusion or spatial heterogeneity. However, the presence of a fractional differential operator causes memory (time fractional) or nonlocality (space fractional) issues, which impose a number of computational constraints. In this paper we develop efficient, scalable techniques for solving fractional-in-space reaction diffusion equations using the finite element method on both structured and unstructured grids, and robust techniques for computing the fractional power of a matrix times a vector. Our approach is show-cased by solving the fractional Fisher and fractional Allen-Cahn reaction-diffusion equations in two and three spatial dimensions, and analysing the speed of the travelling wave and size of the interface in terms of the fractional power of the underlying Laplacian operator
Fourier spectral methods for fractional-in-space reaction-diffusion equations
Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is computationally demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reactiondiffusion equations. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is show-cased by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models,together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator
Numerical solving unsteady space-fractional problems with the square root of an elliptic operator
An unsteady problem is considered for a space-fractional equation in a
bounded domain. A first-order evolutionary equation involves the square root of
an elliptic operator of second order. Finite element approximation in space is
employed. To construct approximation in time, regularized two-level schemes are
used. The numerical implementation is based on solving the equation with the
square root of the elliptic operator using an auxiliary Cauchy problem for a
pseudo-parabolic equation. The scheme of the second-order accuracy in time is
based on a regularization of the three-level explicit Adams scheme. More
general problems for the equation with convective terms are considered, too.
The results of numerical experiments are presented for a model two-dimensional
problem.Comment: 21 pages, 7 figures. arXiv admin note: substantial text overlap with
arXiv:1412.570
An unstructured mesh control volume method for two-dimensional space fractional diffusion equations with variable coefficients on convex domains
In this paper, we propose a novel unstructured mesh control volume method to
deal with the space fractional derivative on arbitrarily shaped convex domains,
which to the best of our knowledge is a new contribution to the literature.
Firstly, we present the finite volume scheme for the two-dimensional space
fractional diffusion equation with variable coefficients and provide the full
implementation details for the case where the background interpolation mesh is
based on triangular elements. Secondly, we explore the property of the
stiffness matrix generated by the integral of space fractional derivative. We
find that the stiffness matrix is sparse and not regular. Therefore, we choose
a suitable sparse storage format for the stiffness matrix and develop a fast
iterative method to solve the linear system, which is more efficient than using
the Gaussian elimination method. Finally, we present several examples to verify
our method, in which we make a comparison of our method with the finite element
method for solving a Riesz space fractional diffusion equation on a circular
domain. The numerical results demonstrate that our method can reduce CPU time
significantly while retaining the same accuracy and approximation property as
the finite element method. The numerical results also illustrate that our
method is effective and reliable and can be applied to problems on arbitrarily
shaped convex domains.Comment: 18 pages, 5 figures, 9 table
- …