Table of Contents
Fetching ...

Stabilizer-free Weak Galerkin Methods for Quad-Curl Problems on polyhedral Meshes without Convexity Assumptions

Chunmei Wang, Shangyou Zhang

TL;DR

This work develops a stabilizer-free weak Galerkin method for 3D quad-curl problems on polyhedral meshes, including non-convex elements, by employing bubble functions to eliminate stabilization terms. The proposed scheme preserves a symmetric, positive definite system, uses discrete weak gradient and curl-curl operators, and achieves optimal discrete-norm error with optimal-order $L^2$ error for $k>2$ (suboptimal for $k=2$). A complete theoretical framework is established, including well-posedness, error equations, energy and $L^2$ error estimates, and duality-based analyses, complemented by extensive numerical tests on convex and non-convex meshes. The method demonstrates robustness, efficiency, and accuracy in both the theoretical and practical settings, extending stabilizer-free WG methods to high-order 3D quad-curl problems on general polyhedral meshes.

Abstract

This paper introduces an efficient stabilizer-free weak Galerkin (WG) finite element method for solving the three-dimensional quad-curl problem. Leveraging bubble functions as a key analytical tool, the method extends the applicability of stabilizer-free WG approaches to non-convex elements in finite element partitions-a notable advancement over existing methods, which are restricted to convex elements. The proposed method maintains a simple, symmetric, and positive definite formulation. It achieves optimal error estimates for the exact solution in a discrete norm, as well as an optimal-order $L^2$ error estimate for $k>2$ and a sub-optimal order for the lowest order case $k=2$. Numerical experiments are presented to validate the method's efficiency and accuracy.

Stabilizer-free Weak Galerkin Methods for Quad-Curl Problems on polyhedral Meshes without Convexity Assumptions

TL;DR

This work develops a stabilizer-free weak Galerkin method for 3D quad-curl problems on polyhedral meshes, including non-convex elements, by employing bubble functions to eliminate stabilization terms. The proposed scheme preserves a symmetric, positive definite system, uses discrete weak gradient and curl-curl operators, and achieves optimal discrete-norm error with optimal-order error for (suboptimal for ). A complete theoretical framework is established, including well-posedness, error equations, energy and error estimates, and duality-based analyses, complemented by extensive numerical tests on convex and non-convex meshes. The method demonstrates robustness, efficiency, and accuracy in both the theoretical and practical settings, extending stabilizer-free WG methods to high-order 3D quad-curl problems on general polyhedral meshes.

Abstract

This paper introduces an efficient stabilizer-free weak Galerkin (WG) finite element method for solving the three-dimensional quad-curl problem. Leveraging bubble functions as a key analytical tool, the method extends the applicability of stabilizer-free WG approaches to non-convex elements in finite element partitions-a notable advancement over existing methods, which are restricted to convex elements. The proposed method maintains a simple, symmetric, and positive definite formulation. It achieves optimal error estimates for the exact solution in a discrete norm, as well as an optimal-order error estimate for and a sub-optimal order for the lowest order case . Numerical experiments are presented to validate the method's efficiency and accuracy.

Paper Structure

This paper contains 9 sections, 15 theorems, 147 equations, 4 figures, 7 tables, 1 algorithm.

Key Result

Theorem 2.1

sun Given $\bf\in H(div^0; \Omega)$, the problem weakform has a unique solution $({\mathbf{u}}; p)\in H_0(curl^2; \Omega)\times H_0^1(\Omega)$. Furthermore, $p=0$ and ${\mathbf{u}}$ satisfies

Figures (4)

  • Figure 9.1: The first three grids for the computation in Tables \ref{['t-1']}--\ref{['t-2']}.
  • Figure 9.2: The first three non-convex polygonal grids for the computation in Tables \ref{['t-3']}--\ref{['t-4']}.
  • Figure 9.3: The first three grids for the computation in Tables \ref{['t-5']}--\ref{['t-6']}.
  • Figure 9.4: The first three (non convex polyhedral) grids for the computation in Table \ref{['t-7']}.

Theorems & Definitions (31)

  • Theorem 2.1
  • Lemma 5.1
  • proof
  • Remark 5.1
  • Lemma 5.2
  • proof
  • Lemma 5.3
  • proof
  • Remark 5.2
  • Lemma 5.4
  • ...and 21 more