Galerkin finite element approximation of a stochastic semilinear fractional subdiffusion with fractionally integrated additive noise
Affiliation
University of Chester, Lvliang University, P. R. China, Indian Institute of Technology BombayMetadata
Show full item recordAbstract
A Galerkin finite element method is applied to approximate the solution of a semilinear stochastic space and time fractional subdiffusion problem with the Caputo fractional derivative of the order $ \alpha \in (0, 1)$, driven by fractionally integrated additive noise. After discussing the existence, uniqueness and regularity results, we approximate the noise with the piecewise constant function in time in order to obtain a regularized stochastic fractional subdiffusion problem. The regularized problem is then approximated by using the finite element method in spatial direction. The mean squared errors are proved based on the sharp estimates of the various Mittag-Leffler functions involved in the integrals. Numerical experiments are conducted to show that the numerical results are consistent with the theoretical findings.Citation
Wenyan Kang, Bernard A Egwu, Yubin Yan, Amiya K Pani, Galerkin finite element approximation of a stochastic semilinear fractional subdiffusion with fractionally integrated additive noise, IMA Journal of Numerical Analysis, 2021;, drab035Publisher
Oxford AcademicAdditional Links
https://academic.oup.com/imajna/advance-article-abstract/doi/10.1093/imanum/drab035/6278467?redirectedFrom=fulltextType
ArticleDescription
This is a pre-copyedited, author-produced PDF of an article accepted for publication in IMA Journal of Numerical Analysis following peer review. The version of record Wenyan Kang, Bernard A Egwu, Yubin Yan, Amiya K Pani, Galerkin finite element approximation of a stochastic semilinear fractional subdiffusion with fractionally integrated additive noise, IMA Journal of Numerical Analysis, 2021;, drab035 is available online at: https://doi.org/10.1093/imanum/drab035EISSN
1464-3642Collections
Related items
Showing items related by title, author, creator and subject.
-
A novel high-order algorithm for the numerical estimation of fractional differential equationsThis paper uses polynomial interpolation to design a novel high-order algorithm for the numerical estimation of fractional differential equations. The Riemann-Liouville fractional derivative is expressed by using the Hadamard finite-part integral and the piecewise cubic interpolation polynomial is utilized to approximate the integral. The detailed error analysis is presented and it is established that the convergence order of the algorithm is O(h4−a). Asymptotic expansion of the error for the presented algorithm is also investigated. Some numerical examples are provided and compared with the exact solution to show that the numerical results are in well agreement with the theoretical ones and also to illustrate the accuracy and efficiency of the proposed algorithm.
-
A finite element method for time fractional partial differential equationsFractional differential equations, particularly fractional partial differential equations (FPDEs) have many applications in areas such as diffusion processes, electromagnetics, electrochemistry, material science and turbulent flow. There are lots of work for the existence and uniqueness of the solutions for fractional partial differential equations. In recent years, people start to consider the numerical methods for solving fractional partial differential equation. The numerical methods include finite difference method, finite element method and the spectral method. In this dissertation, we mainly consider the finite element method, for the time fractional partial differential equation. We consider both time discretization and space discretization. We obtain the optimal error estimates both in time and space. The numerical examples demonstrate that the numerical results are consistent with the theoretical results.
-
Stability regions of numerical methods for solving fractional differential equationsThis dissertation deals with proper consideration of stability regions of well known numerical methods for solving fractional differential equations. It is based on the algorithm by Diethelm [15], predictor-corrector algorithm by Garrappa [31] and the convolution quadrature proposed by Lubich [3]. Initially, we considered the stability regions of numerical methods for solving ordinary differential equation using boundary locus method as a stepping stone of understanding the subject matter in Chapter 4. We extend the idea to the fractional differential equation in the following chapter and conclude that each stability regions of the numerical methods differs because of their differences in weights. They are illustrated by a number of graphs.
