---
title: Two-Grid Weak Galerkin Methods
url: https://www.emergentmind.com/topics/two-grid-weak-galerkin-method
type: topic
---

# Two-Grid Weak Galerkin Methods

The two-grid weak Galerkin method denotes a family of algorithms that combine weak Galerkin (WG) discretization with a coarse–fine decomposition in order to reduce the cost of nonlinear, eigenvalue, or linear algebraic computations while retaining the structural advantages of WG formulations. In the arXiv literature, the term covers at least three technically distinct constructions: a coarse nonlinear solve followed by a fine linearized WG solve for quasi-linear elliptic equations of non-monotone type, a coarse WG eigenproblem followed by a fine WG source problem and Rayleigh quotient for eigenvalue approximation, and a two-level auxiliary-space preconditioner for linear WG systems [2202.11164] [1708.08183] [2508.02065] [1405.7506] [1410.1012].

## 1. Conceptual scope

Within WG theory, “two-grid” does not refer to a single canonical algorithm. The common principle is that the expensive part of the computation is confined to a coarse space, while a fine-space correction recovers most of the accuracy otherwise obtained by a full fine-grid solve.

| Variant | Coarse component | Fine component |
|---|---|---|
| Nonlinear elliptic solver | Nonlinear WG solve on mesh $\mathcal{T}_\tau$ | Linearized WG solve on $\mathcal{T}_h$ |
| Eigenvalue acceleration | WG eigenproblem on mesh $\mathcal{T}_H$ | WG linear source problem on $\mathcal{T}_h$ plus Rayleigh quotient |
| Auxiliary-space two-level method | Conforming $P_1$ auxiliary problem | Smoothing and coarse correction for the fine WG system |

A recurrent misconception is that two-grid WG always means multigrid preconditioning. The literature shows a sharper distinction. In the nonlinear and eigenvalue settings, two-grid is an approximation strategy for the discrete PDE or eigenproblem itself; in the diffusion literature, it is a solver or preconditioner for the linear WG system generated on the fine mesh [2202.11164] [1708.08183] [2508.02065] [1405.7506] [1410.1012].

Another important distinction concerns the role of the coarse space. In the nonlinear and eigenvalue papers, the coarse space is associated with a genuinely coarser mesh with sizes $\tau>h$ or $H>h$. In the auxiliary-space preconditioning papers, the “coarse” space is often the conforming piecewise linear space on the same mesh, used as an auxiliary solver space rather than as a separate physical discretization level [1410.1012] [1405.7506].

## 2. Weak Galerkin structure

The WG framework replaces strong continuity requirements by weakly defined differential operators acting on pairs of interior and boundary unknowns. A generic weak function has the form $v=\{v_0,v_b\}$, where $v_0$ is an element-interior polynomial and $v_b$ is an edge or face polynomial. In stabilized polytopal formulations for quasi-linear elliptic and eigenvalue problems, one commonly uses spaces such as
\[
V_h=\{v=\{v_0,v_b\}: v_0|_T\in P_k(T),\ v_b|_e\in P_k(e)\ \text{or}\ P_{k-1}(e)\},
\]
with homogeneous boundary conditions enforced through the trace component. In RT/BDM-based diffusion formulations, the local choices are instead expressed through triples $(V(T),M(F),W(T))$, where the weak gradient lives in Raviart–Thomas or Brezzi–Douglas–Marini spaces [2202.11164] [1708.08183] [1405.7506].

The core operator is the discrete weak gradient. In several of the cited formulations it is defined locally by
\[
(\nabla_w v,\boldsymbol{\phi})_K
=
-(v_0,\nabla\cdot\boldsymbol{\phi})_K
+
\langle v_b,\boldsymbol{\phi}\cdot \boldsymbol{n}\rangle_{\partial K},
\]
for all local polynomial test vectors $\boldsymbol{\phi}$. For elasticity, the same pattern yields weak strain and weak divergence through $\varepsilon_w(v)=(\nabla_w v+(\nabla_w v)^T)/2$ and an elementwise weak divergence operator [2508.02065].

Stabilization is formulation-dependent. For the non-monotone quasi-linear elliptic problem, the stabilization is
\[
s_h(v_h,w_h)
=
\sum_{K\in \mathcal{T}_h}
h_K^{-1}\langle v_0-v_b,\ w_0-w_b\rangle_{\partial K},
\]
while the Laplace and elasticity eigenvalue papers use weighted stabilizers of the form $h_T^{-1+\varepsilon}$ or $h_T^{-1+\delta}$ on the boundary mismatch. By contrast, the auxiliary-space preconditioning paper for second-order diffusion studies a WG bilinear form without an additional stabilization term, namely
\[
a_h(u_h,v_h)=(\mathbb{A}\nabla_w u_h,\nabla_w v_h).
\]
This difference is structural rather than cosmetic: the particular stabilization enters both the spectral theory and the two-grid error analysis [2202.11164] [1708.08183] [2508.02065] [1410.1012].

Several norm equivalences underwrite the theory. In the nonlinear setting, the mesh-dependent seminorm
\[
|v_h|
=
\Big(
\sum_{K\in\mathcal{T}_h}
(
\|\nabla_w v_h\|_{L^2(K)}^2
+
h_K^{-1}\|v_0-v_b\|_{L^2(\partial K)}^2
)
\Big)^{1/2}
\]
is equivalent to the broken $H^1$-like norm $\|v_h\|_{1,h}$. In the diffusion-preconditioning literature, one similarly has $|v_h|_{1,h}\simeq \|\nabla_w v_h\|$, together with inverse and discrete Poincaré inequalities. These identities are the bridge between WG trace variables and classical energy estimates [2202.11164] [1410.1012].

## 3. Two-grid WG for quasi-linear elliptic problems of non-monotone type

For the quasi-linear model
\[
-\nabla\cdot(a(x,u)\nabla u)=f \quad \text{in }\Omega,
\qquad
u=g \quad \text{on }\partial\Omega,
\]
the coefficient satisfies $a(x,u)\in C^2(\overline{\Omega}\times\mathbb{R})$ with bounded derivatives up to second order and
\[
0<\alpha_0\le a(x,u)\le \alpha_1<\infty.
\]
The weak formulation is to find $u\in H_g^1(\Omega)$ such that
\[
\int_\Omega a(u)\nabla u\cdot \nabla v\,dx
=
\int_\Omega f\,v\,dx,
\qquad
\forall v\in H_0^1(\Omega).
\]
The corresponding WG scheme seeks $u_h=\{u_0,u_b\}\in V_h$ with $u_b=\mathcal{Q}_b g$ on $\partial\Omega$ and
\[
A_h(u_h;u_h,v_h)=(f,v_0),
\qquad
\forall v_h\in V_h^0,
\]
where
\[
A_h(u_h;v_h,w_h)
=
(a(u_0)\nabla_w v_h,\nabla_w w_h)_{\mathcal{T}_h}
+
s_h(v_h,w_h).
\]
The existence theory is based on Brouwer’s fixed point theorem applied to a map constructed from the discrete Fréchet derivative of the nonlinear operator. A discrete Gårding inequality for the linearized bilinear form $D_h$ provides the stability input, and the fixed-point map is shown to map a closed convex ball centered at $\mathcal{Q}_h u$ into itself; for sufficiently small $h$ it is also a contraction, which yields a unique WG solution in that ball [2202.11164].

The single-grid error analysis gives
\[
\|u-u_h\|_{1,h}\le C h^k |u|_{H^{k+1}(\Omega)},
\]
and
\[
\|u-u_0\|_{L^2(\Omega)}
\le
C\,(h^{k+1}+h^{2k}|u|_{H^{k+1}(\Omega)})\,|u|_{H^{k+1}(\Omega)}.
\]
These estimates use the projection commutativity relation $\nabla_w(\mathcal{Q}_h v)=\Pi_h(\nabla v)$, explicit consistency terms, and nonlinear remainder bounds derived from averaged derivatives $\widetilde a_u$ and $\widetilde a_{uu}$ [2202.11164].

The two-grid algorithm introduces a coarse mesh $\mathcal{T}_\tau$ and a fine mesh $\mathcal{T}_h$ with $\tau>h$. First, one solves the full nonlinear WG problem on the coarse grid:
\[
A_\tau(u_\tau;u_\tau,v)=(f,v_0),
\qquad
\forall v\in V_\tau^0.
\]
Second, one freezes the coefficient at the coarse solution and solves a single linear WG problem on the fine grid:
\[
A_h(u_\tau;u^h,v)=(f,v_0),
\qquad
\forall v\in V_h^0.
\]
The output is $u_h^{tg}:=u^h$. The main energy-norm estimate is
\[
\|u-u_h^{tg}\|_{1,h}
=
\|u-u^h\|_{1,h}
\le
C(u)\big(h^k+\tau^{k+1}\big),
\]
with $C(u)$ depending on $|u|_{H^{k+1}(\Omega)}$ and $|u|_{W^1_\infty(\Omega)}$, but not on $h$ or $\tau$. For $k=1$, the choice $\tau=\mathcal{O}(h^{1/2})$ yields the optimal rate in the energy-like norm [2202.11164].

The analysis is driven by a consistency identity comparing the fine-grid nonlinear solution $u_h$ and the two-grid solution $u^h$, together with a splitting
\[
a(u_{h,0})-a(u_{\tau,0})
=
[a(u_{h,0})-a(u)]
+
[a(u)-a(u_{\tau,0})].
\]
After Taylor expansion and use of the single-grid $L^2$ estimates, the coarse-grid error contributes at order $\tau^{k+1}$, while the genuine fine-grid discretization contributes at order $h^k$ [2202.11164].

## 4. Two-grid WG for eigenvalue problems

For eigenvalue problems, the two-grid WG method replaces a fine-grid generalized eigenproblem by a coarse-grid eigenproblem and a fine-grid linear correction. This strategy appears in both Laplace and linear elasticity settings.

For the Dirichlet Laplace eigenproblem, Zhai, Xie, Zhang, and Zhang consider a WG discretization with
\[
a_w(v,w)
=
\sum_{T\in\mathcal{T}_h}(\nabla_w v,\nabla_w w)_T
+
s(v,w),
\qquad
b_w(v,w)
=
\sum_{T\in\mathcal{T}_h}(v_0,w_0)_T,
\]
where
\[
s(v,w)
=
\sum_{T\in\mathcal{T}_h}
h_T^{-1+\varepsilon}
\langle Q_b v_0-v_b,\ Q_b w_0-w_b\rangle_{\partial T}.
\]
The standard WG eigenproblem seeks $(\lambda_h,u_h)$ with $\|u_h\|_b=1$ such that
\[
a_w(u_h,v_h)=\lambda_h b_w(u_h,v_h),
\qquad
\forall v_h\in V_h.
\]
Its single-grid approximation is asymptotically a lower bound: for sufficiently small $h$,
\[
\lambda_j-\lambda_{j,h}\in [C h^{2k},\, C h^{2k-2}],
\]
and
\[
\|u_j-u_{j,h}\|_V\le C h^{k-\varepsilon}\|u_j\|_{k+1},
\qquad
\|u_j-u_{j,h}\|_b\le C h^{k+1-\varepsilon}\|u_j\|_{k+1}.
\]
The two-grid algorithm computes $(\lambda_H,u_H)$ on $V_H$, solves
\[
a_w(\tilde u_h,v_h)=\lambda_H b_w(u_H,v_h),
\qquad
\forall v_h\in V_h,
\]
and then sets
\[
\tilde\lambda_h=\frac{a_w(\tilde u_h,\tilde u_h)}{b_w(\tilde u_h,\tilde u_h)}.
\]
With
\[
\bar k=\min\{2k-2,\ k+\gamma-\varepsilon\},
\]
one obtains
\[
\|Q_h u-\tilde u_h\|_V\lesssim H^{\bar k}+h^{k-\varepsilon},
\qquad
|\tilde\lambda_h-\lambda|
\lesssim
H^{2\bar k}+h^{2k-2}.
\]
For $k\ge 2$ and the favorable dual-regularity regime, this yields the characteristic doubling from a coarse-grid eigenvalue error of order $H^{2k-2}$ to a two-grid contribution of order $H^{4k-4}$. A sufficient condition for preserving the asymptotic lower-bound property is
\[
H^{2\bar k}\le C h^{2k+\delta}
\quad \text{for some } \delta>0,
\]
although the reported computations show that lower bounds may persist even when this condition is not strictly enforced [1708.08183].

The elasticity eigenvalue version developed by Lu and Zhai follows the same template but uses weak strain and weak divergence:
\[
a_w(v_h,w_h)
=
2\mu (\varepsilon_w(v_h),\varepsilon_w(w_h))
+
\lambda (\nabla_w\cdot v_h,\nabla_w\cdot w_h)
+
s(v_h,w_h),
\qquad
b_w(v_h,w_h)=(v_0,w_0),
\]
with
\[
s(v_h,w_h)=\sum_{T\in T_h} h_T^{-1+\delta}\langle v_0-v_b,\ w_0-w_b\rangle_{\partial T},
\qquad
\delta\in(0,1).
\]
After solving the coarse WG eigenproblem and the fine WG linear problem
\[
a_w(\tilde u_h,v_h)=\gamma_H b_w(u_H,v_h),
\qquad
\forall v_h\in V_h,
\]
the fine-grid Rayleigh quotient defines $\tilde\gamma_h$. Writing
\[
\bar k=\min\{2k-2\delta,\ k+2-2\delta\},
\]
the theory gives
\[
|||Q_h u-\tilde u_h|||
\lesssim
H^{\bar k}+h^{k-\delta/2},
\qquad
|\gamma-\tilde\gamma_h|
\lesssim
H^{2\bar k}+h^{2k-2\delta}.
\]
The scaling $h\lesssim H^2$, with the examples using $h\approx (2H)^2$, is the practical regime identified for the “doubling” effect. Lower-bound preservation is obtained if
\[
H^{2\bar k}\lesssim h^{2k+\delta_0}
\quad \text{for some } \delta_0>0,
\]
which ensures $\gamma-\tilde\gamma_h\ge 0$ for sufficiently small $h$ [2508.02065].

In both eigenvalue strands, the decisive observation is that the fine-grid computation is no longer a generalized eigenproblem but a symmetric positive definite linear solve driven by coarse spectral data. This is the principal source of the acceleration [1708.08183] [2508.02065].

## 5. Auxiliary-space two-level algorithms and preconditioners

A distinct branch of the literature interprets two-grid WG as a preconditioning framework for linear systems arising from WG discretizations of diffusion equations. Here the objective is not a new discrete approximation but rapid solution of the fine-grid WG algebraic system.

Li and Xie analyze a two-level algorithm for the WG discretization of
\[
-\operatorname{div}(a\nabla u)=f \quad \text{in }\Omega,
\qquad
u=0 \quad \text{on }\partial\Omega,
\]
on quasi-uniform simplicial meshes. The WG bilinear form is
\[
a_h((u_h,\lambda_h),(v_h,\mu_h))
=
(a\nabla_w(u_h,\lambda_h),\nabla_w(v_h,\mu_h))_\Omega,
\]
with block operator representation
\[
\mathcal A_h=
\begin{bmatrix}
C_h & B_h^t\\
B_h & D_h
\end{bmatrix}.
\]
They prove the norm equivalence
\[
\|(v_h,\mu_h)\|_h^2
\lesssim
a_h((v_h,\mu_h),(v_h,\mu_h))
\lesssim
h^{-2}\|(v_h,\mu_h)\|_h^2,
\]
and consequently
\[
\kappa(\mathcal A_h)=O(h^{-2}).
\]
The auxiliary space is the conforming piecewise linear space
\[
\widetilde V_h=\{v\in H_0^1(\Omega): v|_T\in P_1(T)\ \forall T\in T_h\},
\]
with prolongation $I_h:\widetilde V_h\to W_h$ defined by element and face moment matching. A two-level preconditioner $\mathcal B_h$ is constructed from pre-smoothing, left coarse-grid correction, right coarse-grid correction, and post-smoothing. Its analysis rests on an extended Xu–Zikatanov identity,
\[
(\mathcal B_h^{-1}w_h,w_h)_h
=
\inf_{v_h+I_h\tilde v_h=w_h}
\Big(
\|v_h+\mathcal R_h^t\mathcal A_h I_h\tilde v_h\|_{\overline{\mathcal R}_h^{-1}}^2
+
\|\tilde v_h\|_{\overline{\widetilde R}_h^{-1}}^2
\Big),
\]
which yields a mesh-independent contraction bound and does so without any elliptic regularity assumption [1405.7506].

Chen, Wang, Wang, and Ye develop an auxiliary-space multigrid preconditioner for a WG method for second-order diffusion equations. Their fine WG space is
\[
V_h=\{v=\{v_0,v_b\}: v_0|_{K_0}\in P_j(K_0),\ v_b|_e\in P_l(e)\},
\]
with $l=j$ or $l=j+1$, and the auxiliary space is the conforming piecewise linear space
\[
W_H=V_h^c=\{w\in H_0^1(\Omega): w|_K\in P_1(K)\}.
\]
The prolongation is the $L^2$ projection/inclusion $\Pi w=Q_h w=\{Q_0 w,Q_b w\}$, and the multiplicative two-grid preconditioner is characterized by
\[
I-BA=(I-RA)(I-\Pi B_H\Pi^t)(I-RA).
\]
Under smoother bounds, coarse-solver bounds, prolongation stability, and the existence of a stable quasi-interpolant $P:V_h\to W_H$, the preconditioned operator satisfies
\[
\kappa(BA)\lesssim C
\quad \text{independent of } h.
\]
The paper also derives a reduced edge/face system by eliminating interior unknowns and proves the same type of $h$-independent condition-number bound for the reduced operator $A_r$ [1410.1012].

The conceptual difference between these preconditioning papers and the approximation-oriented two-grid papers is substantial. In the preconditioning setting, the coarse space is part of an auxiliary-space decomposition of a single fine-grid operator; in the nonlinear and eigenvalue settings, coarse and fine meshes define distinct discrete problems whose solutions are combined to approximate the PDE or spectrum itself [1405.7506] [1410.1012].

## 6. Numerical behavior and methodological significance

The numerical literature consistently reports that the two-grid strategy preserves much of fine-grid WG accuracy while reducing either nonlinear/eigenvalue cost or Krylov iteration complexity.

| Setting | Representative observation | Source |
|---|---|---|
| Quasi-linear elliptic, $k=1$ | At $h=1/100$ on rectangular grids, single-grid WG took $\sim 272$ s with $\|u-u_h\|_{1,h}\approx 7.06\times 10^{-2}$; two-grid WG took $\sim 54$ s with $\|u-u^h\|_{1,h}\approx 8.29\times 10^{-2}$ | [2202.11164] |
| Auxiliary-space preconditioner | PCG iterations were essentially independent of $h$; in 2D unit-disk tests, the original WG system required about $13$ iterations and the reduced system about $8$ | [1410.1012] |
| Two-level diffusion solver | With a multigrid coarse solver, average error-reduction factors per iteration were about $0.56$, $0.36$, and $0.26$ for $m=1,2,3$ smoothing steps | [1405.7506] |
| Laplace eigenproblem | For $k=1$, the error $\lambda_1-\tilde\lambda_{1,h}$ decayed from $2.1554\times 10^{-1}$ to $1.3006\times 10^{-2}$ to $8.0627\times 10^{-4}$, with observed order about $4$ | [1708.08183] |
| Elasticity eigenproblem | For two-grid WG, reported eigenfrequency orders were around $3.6$–$3.7$ for $k=1$ and about $7.5$ for $k=2$, with lower bounds observed when $h=(2H)^2$ | [2508.02065] |

Several conclusions are common to these results. First, the WG weak-gradient machinery is compatible with coarse–fine acceleration without sacrificing the discrete structure that supports stability and error estimates. Second, the most favorable mesh couplings are problem-dependent: $\tau=\mathcal{O}(h^{1/2})$ is singled out for the $k=1$ quasi-linear elliptic case, whereas $h\approx H^2$ is the standard coupling in the eigenvalue literature [2202.11164] [1708.08183] [2508.02065]. Third, lower-bound preservation for eigenvalues is not automatic; it requires a quantitative relation between the coarse-grid correction error and the fine-grid lower-bound error, even though computations often remain lower-bounding outside the strict sufficient regime [1708.08183] [2508.02065].

A further methodological point is that WG formulations themselves are not uniform across applications. The cited papers use different polynomial pairings, different stabilizers, and in the diffusion-preconditioning setting even a formulation without an extra stabilization term. Two-grid analysis is therefore inseparable from the particular WG discretization chosen [1410.1012] [1405.7506] [2202.11164].

Taken together, these works establish the two-grid weak Galerkin method as a broad computational paradigm rather than a single algorithm. Its unifying feature is the use of WG trial spaces and weak differential operators together with a coarse–fine decomposition that transfers the dominant computational burden away from the full fine-grid problem, while rigorous theory quantifies when that transfer preserves optimal convergence, asymptotic lower bounds, or mesh-independent conditioning [2202.11164] [1708.08183] [2508.02065] [1405.7506] [1410.1012].

Source: https://www.emergentmind.com/topics/two-grid-weak-galerkin-method