Aubin--Nitsche-type estimates for space-time FOSLS for parabolic PDEs

Thomas Führer,Gregor Gantner
2024-12-06
Abstract:We develop Aubin--Nitsche-type estimates for recently proposed first-order system least-squares finite element methods (FOSLS) for the heat equation. Under certain assumptions, which are satisfied if the spatial domain is convex and the heat source and initial datum are sufficiently smooth, we prove that the $L^2$ error of approximations of the scalar field variable converges at a higher rate than the overall error. Furthermore, a higher-order conservation property is shown. In addition, we discuss quasi-optimality in weaker norms. Numerical experiments confirm our theoretical findings.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is: **Develop Aubin–Nitsche - type estimates to prove that under certain regularity assumptions, the error convergence rate in the \(L^2\) norm of the First - Order System Least - Squares finite element method (FOSLS) for the heat equation is higher than the overall error convergence rate**. Specifically, the goals of the paper include: 1. **Prove a higher convergence rate of the \(L^2(Q)\) error**: - In Theorem 13, the author proves that in a convex region and with appropriately smooth heat sources and initial conditions, the \(L^2(Q)\) error of the scalar field variable \(u\) converges faster than the overall error in the \(U\) norm. The specific formula is: \[ \|u - u_h\|_{L^2(Q)} \leq \gamma(h)\|u - u_h\|_U \] where \(0\leq\gamma(h)\to 0\) as \(h\to 0\), and \(h\) is the discretization parameter (mesh size). 2. **Prove high - order conservation properties**: - In Theorem 14, the author presents another main result, namely: \[ \|\Pi f-\text{div}_{t,x}u_h\|_{L^2(Q)} \leq \gamma(h)\|u - u_h\|_U \] This can be interpreted as high - order mass conservation. 3. **Discuss quasi - optimality in weak norms**: - The paper also explores quasi - optimality in weaker norms and points out that in some cases, quasi - optimality can be achieved by removing the divergence term. To achieve these goals, the paper utilizes the Aubin–Nitsche duality argument technique and considers the anisotropic regularity results of parabolic problems. In addition, the author also uses recent findings on the approximation and boundedness of interpolation operators in finite element spaces. Through numerical experiments, the paper verifies the validity of the theoretical results, indicating that the \(L^2\) error does converge at a higher rate.