The Optimal <i>L</i><sup>2</sup>-Norm Error Estimate of a Weak Galerkin Finite Element Method for a Multi-Dimensional Evolution Equation with a Weakly Singular Kernel

This paper proposes a weak Galerkin (WG) finite element method for solving a multi-dimensional evolution equation with a weakly singular kernel. The temporal discretization employs the backward Euler scheme, while the fractional integral term is approximated via a piecewise constant function method....

Full description

Saved in:
Bibliographic Details
Main Authors: Haopan Zhou, Jun Zhou, Hongbin Chen
Format: Article
Language:English
Published: MDPI AG 2025-06-01
Series:Fractal and Fractional
Subjects:
Online Access:https://www.mdpi.com/2504-3110/9/6/368
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This paper proposes a weak Galerkin (WG) finite element method for solving a multi-dimensional evolution equation with a weakly singular kernel. The temporal discretization employs the backward Euler scheme, while the fractional integral term is approximated via a piecewise constant function method. A fully discrete scheme is constructed by integrating the WG finite element approach for spatial discretization. <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msup><mi>L</mi><mn>2</mn></msup></semantics></math></inline-formula>-norm stability and convergence analysis of the fully discrete scheme are rigorously established. Numerical experiments are conducted to validate the theoretical findings and demonstrate optimal convergence order in both spatial and temporal directions. The numerical results confirm that the proposed method achieves an accuracy of the order <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>O</mi><mfenced separators="" open="(" close=")"><mi>τ</mi><mo>+</mo><msup><mi>h</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup></mfenced></mrow></semantics></math></inline-formula>, where <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi>τ</mi></semantics></math></inline-formula> and <i>h</i> represent the time step and spatial mesh size, respectively. This work extends previous studies on one-dimensional problems to higher spatial dimensions, providing a robust framework for handling evolution equations with a weakly singular kernel.
ISSN:2504-3110