---
title: Discrete Stokes Polygonal Complex
url: https://www.emergentmind.com/topics/discrete-stokes-polygonal-complex
type: topic
---

# Discrete Stokes Polygonal Complex

Discrete Stokes polygonal complex most commonly denotes a compatible discretization of the two-dimensional Stokes sequence on polygonal meshes in which a scalar stream-function space, a velocity space, and a pressure space are linked by discrete curl and divergence operators in an exact sequence. The clearest polygonal realization in the cited literature is the Virtual Element Method sequence
\[
0 \xrightarrow[]{\,i\,} \Phi_h \xrightarrow[]{\,\curl\,} V_h \xrightarrow[]{\,\div\,} Q_h \xrightarrow[]{\,0\,} 0,
\]
which reproduces the continuous Stokes complex on simply connected polygonal domains and yields exactly divergence-free discrete velocities [1807.10650]. Closely related constructions appear on polygonal or polyhedral cell complexes through generalized finite element systems, Morley-type and Scott–Vogelius-type virtual spaces, weak Galerkin bases, compatible discrete operator schemes, and quadrilateral or rectangular nonconforming complexes [1605.08657].

## 1. Continuous sequence and the meaning of exactness

On a simply connected polygonal domain \(\Omega\subset \mathbb R^2\), the continuous Stokes complex used throughout this literature is
\[
0 \xrightarrow[]{\,i\,} H_0^2(\Omega) \xrightarrow[]{\,\curl\,} [H_0^1(\Omega)]^2 \xrightarrow[]{\,\div\,} L_0^2(\Omega) \xrightarrow[]{\,0\,} 0.
\]
With the convention
\[
\curl \psi := \left(\frac{\partial \psi}{\partial y}, -\frac{\partial \psi}{\partial x}\right),
\]
every divergence-free \(H^1\)-velocity is the curl of a scalar stream function, and every zero-mean pressure is the divergence of an \(H_0^1\)-field. Exactness therefore identifies the divergence-free kernel with the image of the scalar potential space rather than only with a weak constraint [1807.10650].

This exactness has several concrete consequences in the cited works. It yields pointwise or exact discrete incompressibility, supports equivalent stream-function/curl formulations of Stokes or Navier–Stokes, and separates velocity approximation from spurious pressure contamination. In the Morley-type stream-function formulation for Navier–Stokes, the same continuous identities are used in the form
\[
u=\operatorname{curl}\psi,\qquad \omega=\operatorname{rot}u=-\Delta\psi,
\]
so that the Stokes sequence becomes the structural device for recovering velocity, vorticity, and pressure from a scalar stream approximation [2212.02173].

## 2. Exact polygonal complexes in the Virtual Element Method

The most explicit polygonal discrete Stokes complex in the supplied literature is the two-dimensional Virtual Element construction of Beirão da Veiga, Russo, and Vacca. On general polygonal meshes satisfying the standard assumptions that each element \(E\) is star-shaped with respect to a ball of radius comparable to \(h_E\) and that the distance between any two vertices is bounded below by a multiple of \(h_E\), the pressure space is
\[
Q_h := \left\{ q \in L_0^2(\Omega) : q|_E \in P_{k-1}(E)\ \forall E\in\Omega_h \right\}.
\]
The velocity space \(V_h\subset [H_0^1(\Omega)]^2\) is an \(H^1\)-conforming vector VEM space with polynomial degree \(k\), and the new scalar stream-function space \(\Phi_h\subset H_0^2(\Omega)\) is a global \(H^2\)-conforming, \(C^1\)-virtual space on polygons. The local and global exact sequences are
\[
\mathbb R \xrightarrow[]{\,i\,} \Phi_h^E \xrightarrow[]{\,\curl\,} V_h^E \xrightarrow[]{\,\div\,} P_{k-1}(E) \xrightarrow[]{\,0\,} 0,
\]
and
\[
0 \xrightarrow[]{\,i\,} \Phi_h \xrightarrow[]{\,\curl\,} V_h \xrightarrow[]{\,\div\,} Q_h \xrightarrow[]{\,0\,} 0.
\]
The key identities are
\[
\curl \Phi_h^E = Z_h^E,\qquad \curl \Phi_h = Z_h,
\]
so every discrete divergence-free velocity in \(V_h\) is exactly the curl of a discrete stream function in \(\Phi_h\), while
\[
\div V_h \subseteq Q_h.
\]
The same paper also gives a reduced exact complex with piecewise constant pressure [1807.10650].

The local velocity space is defined through edge traces, divergence moments, and moments against \(\mathbf x^\perp P_{k-1}(E)\), while the local stream-function space is defined through polynomial traces, polynomial gradient traces, and a biharmonic residual in \(P_{k-1}(E)\). The degrees of freedom are chosen so that the canonical projections
\[
\Pi_k^{\nabla,E},\qquad \Pi_k^{0,E},\qquad \boldsymbol{\Pi}_{k-1}^{0,E}\nabla
\]
are computable for \(V_h^E\), and the stream-function degrees of freedom determine the velocity degrees of freedom of \(\curl\varphi_h\). This computability is essential because virtual functions are not known in closed form inside each polygon.

The exact sequence is also used to rewrite the discrete Navier–Stokes problem on the smaller scalar space. Since
\[
Z_h=\curl\Phi_h,
\]
the mixed formulation and the curl formulation produce the same discrete velocity \(u_h=\curl\psi_h\). Numerically, the curl formulation uses \(2(n_P-1)\) fewer global unknowns than the mixed formulation, but its condition number behaves like \(h^{-4}\), versus \(h^{-2}\) for the velocity-pressure formulation. The same study reports that the predicted convergence rates are confirmed and that the three trilinear form choices behave almost identically on the tested problems [1807.10650].

## 3. Polygonal VEM variants: divergence-free, Morley-type, and Scott–Vogelius-type constructions

A second major polygonal line is the divergence-free VEM for the Stokes problem on polygonal meshes. There the local velocity space is again a Stokes-type virtual space, but the decisive structural property is stated directly as
\[
\operatorname{div}\mathbf V_h = Q_h,\qquad \mathbf Z_h := \ker\bigl(b|_{\mathbf V_h\times Q_h}\bigr)\subseteq \mathbf Z.
\]
Because \(\operatorname{div}\mathbf u_h\in P_{k-1}(K)\) on each element and the discrete incompressibility equation is imposed against all of \(Q_h\), the final discrete velocity is pointwise divergence-free. The same paper proves that the full mixed problem is immediately equivalent to a reduced problem in which the velocity is unchanged and the reduced pressure is the elementwise constant projection of the full pressure; the saving is \(n_P\bigl((k+1)k-2\bigr)\) degrees of freedom [1510.01655].

The Morley-type virtual element method for Navier–Stokes uses a different polygonal pairing. Its stream-function space \(V_h\) is a nonconforming \(H^2\)-type polygonal analogue of the Morley element, its companion vector space \(\mathbf V_h\) is a Crouzeix–Raviart-type virtual element space, and the discrete Stokes-complex identity actually used is
\[
\operatorname{curl} V_h=\widehat{\mathbf V}_h,
\]
where \(\widehat{\mathbf V}_h\) is the subspace of \(\mathbf V_h\) with vanishing discrete divergence pairing against the piecewise constant pressure space \(Q_h\). In this setting the complex is the mechanism for pressure recovery: once \(\psi_h\in V_h\) is known, the recovered pressure \(p_h\in Q_h\) is obtained from an auxiliary mixed problem on \(\mathbf V_h\times Q_h\), and the paper proves optimal-order estimates for recovered velocity, vorticity, and pressure on general polygonal meshes, with domains not necessarily convex [2212.02173].

A different polygonal extension is the VEM generalization of the Scott–Vogelius method. Here the velocity space is obtained by applying the standard scalar conforming VEM space componentwise and pairing it with discontinuous elementwise polynomial pressures,
\[
Q_h=\mathbb P_{\underline k}(\mathcal T_h)\cap L_0^2(\Omega).
\]
This construction is explicitly presented as a polygonal generalization of Scott–Vogelius, but not as a full discrete Stokes complex. Its compatibility is expressed instead by the exact identity
\[
b_h(\mathbf v_h,q_h)=b(\mathbf v_h,q_h)\qquad \forall \mathbf v_h\in \mathbf V_h,\ \forall q_h\in Q_h,
\]
and by the projected divergence-free property
\[
\Pi_{\underline k}^0(\operatorname{div}\mathbf u_h)=0.
\]
For meshes satisfying the standard VEM assumptions, the paper proves a discrete inf-sup condition for \(k\ge 2\), while the lowest-order case on triangular and square meshes reproduces the classical unstable situations [2112.13292].

## 4. Generalized and compatible discrete frameworks on cell complexes

Beyond VEM, the supplied literature contains two broader structural frameworks. The first is the generalized finite element systems framework for smooth differential forms. It extends finite element systems by allowing pullback, full trace, and double-trace restrictions, and it is formulated on a cellular complex. In dimension two, it produces low-order and high-order exact complexes with regularity \((2,1,0)\) and \((2,1+,1)\). For example, on a Clough–Tocher split \(R\) of a triangle, one exact sequence is
\[
0 \to \mathbb{R}\to C^1\mathcal P^3\Lambda^0(R) \xrightarrow{d} C_d^0\mathcal P^2\Lambda^1(R) \xrightarrow{d} C^0\mathcal P^1\Lambda^2(R)\to 0,
\]
and another is
\[
0\to \mathbb{R}\to C^1\mathcal P^3\Lambda^0(R) \xrightarrow{d} C^0\mathcal P^2\Lambda^1(R) \xrightarrow{d} \mathcal P^1\Lambda^2(R)\to 0.
\]
The last two spaces in these complexes provide conforming Stokes pairs with continuous or discontinuous pressure. The framework also guarantees commuting interpolators and identifies some of the resulting spaces as minimal. The paper is explicit that the abstract framework applies to a cellular complex, but its concrete local constructions are composite polynomial spaces on refined simplices rather than arbitrary single polygonal cells [1605.08657].

The second framework is the compatible discrete operator approach on three-dimensional polyhedral meshes. There the discrete unknowns are organized as cochains on primal and dual complexes, the topological laws are discretized exactly by incidence matrices, and the constitutive relations are approximated by discrete Hodge operators. The discrete gradient, curl, and divergence operators satisfy the cochain identities
\[
CG=0,\qquad DC=0,
\]
and the schemes preserve local mass and momentum conservation. Two families are analyzed: one with pressure degrees of freedom at mesh vertices and one with pressure degrees of freedom at cells. Discrete stability is proved by new discrete Poincaré inequalities, and first-order error estimates are derived by commutator arguments for smooth solutions. This is a polyhedral, curl-based realization of the same compatible philosophy, but expressed in mimetic and cochain language rather than in virtual element exact sequences [1401.7842].

A third compatible route is the weak Galerkin divergence-free method on polygons and polyhedra. For the lowest-order weak Galerkin element,
\[
V_h=\left\{\{\mathbf v_0,\mathbf v_b\}: \{\mathbf v_0,\mathbf v_b\}|_T\in [P_1(T)]^d\times [P_0(e)]^d\right\},
\qquad
W_h=\{q\in L_0^2(\Omega): q|_T\in P_0(T)\},
\]
the paper defines the discrete divergence-free space
\[
D_h=\{\mathbf v\in V_h:\ b(\mathbf v,q)=0\ \forall q\in W_h\},
\]
and constructs explicit bases of \(D_h\) on general polygonal and polyhedral meshes. In two dimensions the basis decomposes into element-interior functions, edge-tangential functions, and vertex-hull functions; in three dimensions it decomposes into element-interior functions, face-tangential functions, and edge-solid-hull functions. The resulting reduced system on \(D_h\) is symmetric and positive definite [1602.08815].

## 5. Quadrilateral and rectangular discrete Stokes complexes

A separate line of work develops discrete Stokes complexes on quadrilateral and rectangular meshes using nonconforming finite elements made of piecewise polynomials directly on physical cells. On convex quadrilateral grids, Zhang constructs a stable nonconforming Stokes pair and proves the exact sequence
\[
0 \longrightarrow M_{h0}^Q \xrightarrow{\operatorname{curl}_h} V_{h0}^{QLTZ} \xrightarrow{\operatorname{div}_h} W_h^Q \longrightarrow 0.
\]
Here \(M_{h0}^Q\) is a quadrilateral Morley space,
\[
P_M^Q=P_2(Q)+\operatorname{span}\{\xi^3,\eta^3\},
\]
the scalar velocity space underlying \(V_{h0}^{QLTZ}\) is
\[
P_Q^{QLTZ}=P_1(Q)+\operatorname{span}\{\xi^2,\eta^2\},
\]
and the pressure space is discontinuous \(P_1\) on each quadrilateral. The paper also proves the commuting relations
\[
\operatorname{curl}_h \Pi_h^{QM} = \Pi_h^{QLTZ}\operatorname{curl}, \qquad
\operatorname{div}_h \Pi_h^{QLTZ} = \Pi_0 \operatorname{div},
\]
and extends the construction to mixed meshes consisting of triangles and convex quadrilaterals. In that precise sense, the quadrilateral complex becomes a polygonal-type complex on meshes built from two polygonal cell types [1310.4785].

On uniform rectangular meshes, a different nonconforming complex is
\[
0 \longrightarrow W_h \xrightarrow{\operatorname{curl}_h} \mathbf V_h \xrightarrow{\operatorname{div}_h} P_h \longrightarrow 0,
\]
with a commuting diagram
\[
\operatorname{curl}_h \mathcal I_h = \mathbf \Pi_h \operatorname{curl}, \qquad
\operatorname{div}_h \mathbf \Pi_h = \mathcal P_h \operatorname{div}.
\]
The scalar space \(W_h\) is a 12-DoF rectangular non-\(C^0\) plate element, \(\mathbf V_h\) is a 12-DoF vector space with weak continuity of normal moments and tangential averages, and \(P_h\) is piecewise constant pressure. Although the method is nonconforming and only rigorous for uniform rectangular partitions, it provides an exact commuting discrete Stokes complex and unexpectedly high accuracy: \(O(h^2)\) velocity error in the discrete \(H^1\)-norm, \(O(h^3)\) velocity error in \(L^2\) on convex domains, and \(O(h^2)\) postprocessed pressure accuracy [1812.05823].

## 6. Divergence-preserving reconstruction, pressure robustness, and three-dimensional directions

Several polygonal methods in the supplied literature show that an exact or projected divergence-free pair does not by itself settle the full structure of a discrete Stokes complex; the right-hand side and reconstruction operators must also preserve the relevant orthogonality. The clearest statement appears in the paper on “really pressure-robust” VEM. Starting from a divergence-free polygonal VEM, it argues that standard \(L^2\)-best approximation of virtual test functions destroys the orthogonality to gradient forces, and it introduces a local Raviart–Thomas reconstruction \(I_{\mathrm{RT}_m}\) on a polygonal subtriangulation with the commuting identity
\[
\operatorname{div}(I_{\mathrm{RT}_m}(\boldsymbol v_h)) = \pi_m(\operatorname{div}\boldsymbol v_h).
\]
For \(m=k-1\), divergence is preserved exactly. The reconstruction also satisfies
\[
\int_K (\boldsymbol v_h-I_{\mathrm{RT}_m}\boldsymbol v_h)\cdot \boldsymbol q_h\,dx=0
\qquad \forall \boldsymbol q_h\in \boldsymbol P_{m-1}(K),
\]
and it restores the pressure-robust annihilation of gradient forces. This provides a computable \(H(\mathrm{div})\)-conforming transfer from polygonal virtual spaces to Raviart–Thomas spaces on local simplicial refinements, which is a distinctly complex-like ingredient even though no full exact polygonal sequence is written down [2002.01830].

A related but different polygonal method is the lowest-order staggered discontinuous Galerkin discretization on general convex polygonal meshes. Its unknowns live on a primal polygonal partition plus a simplicial star-subdivision, and the right-hand side is modified by a divergence-preserving reconstruction \(\Pi^{RT}\) into a polygonal \(H(\mathrm{div})\) space. The key identity
\[
(\nabla p,\Pi^{RT}\boldsymbol v)=b_h^*(\pi_h p,\boldsymbol v)
\]
is the mechanism behind pressure robustness. The reconstructed discrete velocity \(\Pi^{RT}\boldsymbol u_h\) is exactly divergence free, and the method achieves optimal first-order convergence together with a superconvergent estimate
\[
\|\mathcal I_h\boldsymbol u-\boldsymbol u_h\|_0 \le C h^2\big(\|\boldsymbol u\|_2+\|\Delta \boldsymbol u\|_2\big)
\]
under stronger regularity [2007.00298].

Three-dimensional extensions remain materially harder. The Freudenthal-mesh paper states that conforming exactly incompressible discretizations are now well understood in two dimensions but remain poorly understood in three dimensions. It formulates two conjectures for the Scott–Vogelius pair on uniform meshes: inf-sup stability for velocity degree \(k\ge 4\), while the best result available in the literature is for \(k\ge 6\), and the existence of a stable space decomposition of the kernel of the divergence for \(k\ge 5\). It also presents numerical evidence supporting these conjectures [2211.05494]. This suggests that a fully satisfactory three-dimensional polygonal or polyhedral analogue of the two-dimensional exact Stokes polygonal complex remains an open direction rather than a closed theory.

Source: https://www.emergentmind.com/topics/discrete-stokes-polygonal-complex