arrow
Volume 13, Issue 4
$L2-1_\sigma$ Finite Element Method for Time-Fractional Diffusion Problems with Discontinuous Coefficients

Yanping Chen, Xuejiao Tan & Yunqing Huang

East Asian J. Appl. Math., 13 (2023), pp. 813-834.

Published online: 2023-10

Export citation
  • Abstract

A time-fractional diffusion equation with an interface problem caused by discontinuous coefficients is considered. To solve it, in the temporal direction Alikhanov’s $L2-1_σ$ method with graded mesh is presented to deal with the weak singularity at $t = 0,$ while in the spatial direction a finite element method with uniform mesh is employed to handle the discontinuous coefficients. Then, with the help of discrete fractional Grönwall inequality and the robustness theory of $α → 1^−,$ we show that the method has stable error bounds at $α → 1^−,$ the fully discrete schemes $L^2(Ω)$ norm and $H^1(Ω)$ semi-norm are unconditionally stable, and the optimal convergence order is $\mathscr{O}(h^2 + N^{−{\rm min}\{rα,2\}})$ and $\mathscr{O}(h + N^{−{\rm min}\{rα,2\}}),$ respectively, where, $h,$ $N,$ $α,$ $r$ is the total number of spatial parameter, the time-fractional order coefficient, and the time grid constant. Finally, three numerical examples are provided to illustrate our theoretical results.

  • AMS Subject Headings

65M10, 78A48

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{EAJAM-13-813, author = {Chen , YanpingTan , Xuejiao and Huang , Yunqing}, title = {$L2-1_\sigma$ Finite Element Method for Time-Fractional Diffusion Problems with Discontinuous Coefficients}, journal = {East Asian Journal on Applied Mathematics}, year = {2023}, volume = {13}, number = {4}, pages = {813--834}, abstract = {

A time-fractional diffusion equation with an interface problem caused by discontinuous coefficients is considered. To solve it, in the temporal direction Alikhanov’s $L2-1_σ$ method with graded mesh is presented to deal with the weak singularity at $t = 0,$ while in the spatial direction a finite element method with uniform mesh is employed to handle the discontinuous coefficients. Then, with the help of discrete fractional Grönwall inequality and the robustness theory of $α → 1^−,$ we show that the method has stable error bounds at $α → 1^−,$ the fully discrete schemes $L^2(Ω)$ norm and $H^1(Ω)$ semi-norm are unconditionally stable, and the optimal convergence order is $\mathscr{O}(h^2 + N^{−{\rm min}\{rα,2\}})$ and $\mathscr{O}(h + N^{−{\rm min}\{rα,2\}}),$ respectively, where, $h,$ $N,$ $α,$ $r$ is the total number of spatial parameter, the time-fractional order coefficient, and the time grid constant. Finally, three numerical examples are provided to illustrate our theoretical results.

}, issn = {2079-7370}, doi = {https://doi.org/10.4208/eajam.2022-178.101022 }, url = {http://global-sci.org/intro/article_detail/eajam/22064.html} }
TY - JOUR T1 - $L2-1_\sigma$ Finite Element Method for Time-Fractional Diffusion Problems with Discontinuous Coefficients AU - Chen , Yanping AU - Tan , Xuejiao AU - Huang , Yunqing JO - East Asian Journal on Applied Mathematics VL - 4 SP - 813 EP - 834 PY - 2023 DA - 2023/10 SN - 13 DO - http://doi.org/10.4208/eajam.2022-178.101022 UR - https://global-sci.org/intro/article_detail/eajam/22064.html KW - Time-fractional, interface problems, finite element, $L2-1_σ$ method, weak singularity. AB -

A time-fractional diffusion equation with an interface problem caused by discontinuous coefficients is considered. To solve it, in the temporal direction Alikhanov’s $L2-1_σ$ method with graded mesh is presented to deal with the weak singularity at $t = 0,$ while in the spatial direction a finite element method with uniform mesh is employed to handle the discontinuous coefficients. Then, with the help of discrete fractional Grönwall inequality and the robustness theory of $α → 1^−,$ we show that the method has stable error bounds at $α → 1^−,$ the fully discrete schemes $L^2(Ω)$ norm and $H^1(Ω)$ semi-norm are unconditionally stable, and the optimal convergence order is $\mathscr{O}(h^2 + N^{−{\rm min}\{rα,2\}})$ and $\mathscr{O}(h + N^{−{\rm min}\{rα,2\}}),$ respectively, where, $h,$ $N,$ $α,$ $r$ is the total number of spatial parameter, the time-fractional order coefficient, and the time grid constant. Finally, three numerical examples are provided to illustrate our theoretical results.

Chen , YanpingTan , Xuejiao and Huang , Yunqing. (2023). $L2-1_\sigma$ Finite Element Method for Time-Fractional Diffusion Problems with Discontinuous Coefficients. East Asian Journal on Applied Mathematics. 13 (4). 813-834. doi:10.4208/eajam.2022-178.101022
Copy to clipboard
The citation has been copied to your clipboard