Galerkin finite element approximation of a stochastic semilinear fractional subdiffusion with fractionally integrated additive noise
AffiliationUniversity of Chester, Lvliang University, P. R. China, Indian Institute of Technology Bombay
MetadataShow full item record
AbstractA 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.
CitationWenyan 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
DescriptionThis 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/drab035
Showing items related by title, author, creator and subject.
A novel high-order algorithm for the numerical estimation of fractional differential equationsAsl, Mohammad S.; Javidi, Mohammad; Yan, Yubin; University of Tabriz; University of Chester (Elsevier, 2018-01-09)This 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 equationsYan, Yubin; Atallah, Samia A. (University of Chester, 2011-09)Fractional 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 equationsYan, Yubin; Nwajeri, Kizito U. (University of Chester, 2012-09)This 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 , predictor-corrector algorithm by Garrappa  and the convolution quadrature proposed by Lubich . 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.