A Systematic Study on Weak Galerkin Finite Element Method for Second Order Parabolic Problems

TL;DR

This paper provides a comprehensive numerical analysis of the weak Galerkin finite element method for second order linear parabolic problems, establishing convergence and demonstrating robustness through numerical experiments.

Contribution

It extends the analysis of WG methods from elliptic to parabolic problems, allowing various polynomial degrees and establishing convergence in multiple norms.

Findings

Convergence in $L^{ abla} (L^2)$ and $L^{ abla} (H^1)$ norms for semidiscrete and fully discrete solutions.

Validation of robustness, reliability, and accuracy through numerical experiments.

Extension of WG method analysis to time-dependent parabolic problems.

Abstract

A systematic numerical study on weak Galerkin (WG) finite element method for second order linear parabolic problems is presented by allowing polynomial approximations with various degrees for each local element. Convergence of both semidiscrete and fully discrete WG solutions are established in $L^{\infty}(L^2)$ and $L^{\infty}(H^1)$ norms for a general WG element $({\cal P}_{k}(K),\;{\cal P}_{j}(\partial K),\;\big[{\cal P}_{l}(K)\big]^2)$, where $k\ge 1$, $j\ge 0$ and $l\ge 0$ are arbitrary integers. The fully discrete space-time discretization is based on a first order in time Euler scheme. Our results are intended to extend the numerical analysis of WG methods for elliptic problems [J. Sci. Comput., 74 (2018), 1369-1396] to parabolic problems. Numerical experiments are reported to justify the robustness, reliability and accuracy of the WG finite element method.