---
title: Reduced Oseen Discretization Techniques
url: https://www.emergentmind.com/topics/reduced-oseen-discretization
type: topic
---

# Reduced Oseen Discretization Techniques

Searching arXiv for recent work on reduced Oseen discretization and closely related Oseen elimination formulations.
Reduced Oseen discretization denotes a class of formulations for the steady Oseen problem in which pressure variables, or other local auxiliary variables attached to the PDE residual, are eliminated from the global algebraic system so that the final solve is posed in terms of velocity unknowns only. In the embedded Trefftz-DG construction for the Oseen problem, the reduced formulation is obtained from a full discontinuous Galerkin system by local inversion and static condensation; the resulting system is “posed in terms of the velocity unknown only,” and this is identified as “a crucial step in the analysis especially for the nonlinear Navier-Stokes problem in Part II” [2606.13229]. Closely related Oseen discretizations realize pressure elimination differently, including a pseudostress reformulation in mixed virtual elements and exactly divergence-free H(div)-conforming mixed methods whose velocity estimates are pressure-robust [2601.20050] [2511.21996].

## 1. Continuous Oseen setting

In the steady Oseen model underlying the reduced Trefftz-DG discretization, one seeks \((u,p)\) satisfying
\[
-\nu \Delta u + (w\cdot \nabla)u + \nabla p = f \quad \text{in }\Omega,
\]
\[
\nabla\cdot u = 0 \quad \text{in }\Omega,
\]
\[
u = 0 \quad \text{on }\partial\Omega,
\]
with \(\Omega\subset\mathbb{R}^d\) \((d=2,3)\), viscosity \(\nu>0\), given convective field \(w\), and body force \(f\). The weak form is: find \((u,p)\in [H_0^1(\Omega)]^d\times L_0^2(\Omega)\) such that for all \((v,q)\in [H_0^1(\Omega)]^d\times L_0^2(\Omega)\),
\[
a_0(u,v)+c(w;u,v)+b(v,p)=(f,v)_-,
\qquad
b(u,q)=0,
\]
where
\[
a_0(u,v)=\nu(\nabla u,\nabla v)_-,
\qquad
c(w;u,v)=((w\cdot\nabla)u,v)_-,
\qquad
b(v,p)=-(\nabla\cdot v,p)_- .
\]
This is the linearized incompressible flow problem to which the reduced discretization is applied [2606.13229].

Related Oseen formulations in the recent literature enlarge this baseline operator. One steady Oseen model adds an advection field \(b\) and a reaction coefficient \(c\),
\[
-\nu\Delta u + (b\cdot\nabla)u + c\,u + \nabla p = f,
\]
under the condition \(c(x)-\tfrac12\nabla\cdot b(x)\ge r_0>0\) on \(\Omega\) [2511.21996]. A two-dimensional generalized Oseen problem further includes a permeability term \(\kappa u\),
\[
-\nu\Delta u + (\beta\cdot\nabla)u + \kappa u + \nabla p = f,
\]
together with \(\mathrm{div}\,u=0\), homogeneous Dirichlet data, and \(\int_\Omega p=0\) [2601.20050]. These variants clarify that “reduced” refers to the discretization strategy rather than to a simplification of the continuous model.

## 2. Trefftz-DG spaces and local complement construction

The embedded Trefftz-DG method is built on a shape-regular simplicial mesh \(\mathcal{T}_h\) of \(\Omega\), with elementwise diameter \(h_T\), global \(h:=\max_T h_T\), and skeleton \(\mathcal{F}\). The ambient discrete space is
\[
X:= [P^k(\mathcal{T}_h)]^d \times P^{k-1}(\mathcal{T}_h)/\mathbb{R},
\]
while the test-side “PDE-residual” space is
\[
Q:= [P^{k-2}(\mathcal{T}_h)]^d \times P^{k-1}(\mathcal{T}_h).
\]
On each element \(T\), the local PDE operator \(B_w:X(T)\to Q(T)'\) is defined by
\[
\langle B_w(u,p),(r,s)\rangle_T
=
h_T\nu^{-1/2}(-\nu\Delta u + w\cdot\nabla u + \nabla p, r)_T
+
\nu^{1/2}(-\nabla\cdot u, s)_T .
\]
The Trefftz subspace is then
\[
T_w
=
\{(u_h,p_h)\in X:\ \langle B_w(u_h,p_h),(r,s)\rangle_T=0\ \forall (r,s)\in Q(T),\ \forall T\}.
\]
It is a kernel-in-a-relaxed-sense space: the Oseen operator and divergence constraint vanish against the residual test space on each cell [2606.13229].

A second ingredient is a local complement \(L(T)=L^u(T)\times L^p(T)\) with \(\dim L(T)=\dim Q(T)\). Its pressure part is
\[
L^p(T):=\{p_h\in P^{k-1}(T): \int_T p_h=0\}.
\]
Its velocity part is chosen as
\[
L^u(T)=\mathrm{Span}\Bigl\{ [B_0^{(k)}(T)]^d,\ \chi=x-x_T,\ n_{T,1},\ldots,n_{T,d}\Bigr\}
\]
subject to zero-mean normal flux on \(d\) of the \(d+1\) faces, where \([B_0^{(k)}(T)]^d := \{ P^k(T)\cap H_0^1(B_T)\}^d\). One checks
\[
X(T)=T_w(T)\oplus L(T),
\]
and, under a mild resolution assumption \(|w|_{h,d}\lesssim 1\), the local operator \(B_w\) is bijective on \(L(T)\) with stable inverse [2606.13229].

This decomposition is the structural basis of the reduction. The Trefftz part carries the globally coupled information, while the complement isolates those local modes on which the residual operator is stably invertible. A plausible implication is that the algebraic elimination is not an ad hoc post-processing step but is encoded in the space design itself.

## 3. Full embedded Trefftz-DG variational problem

The full formulation uses the usual SIPG-Stokes forms on \(X\). The viscous form is
\[
a_h(u_h,v_h)=\sum_T\Bigl[
\nu(\nabla u_h,\nabla v_h)_T
-
(\nu\partial_n u_h,\{v_h\})_{\partial T}
-
(\nu\partial_n v_h,\{u_h\})_{\partial T}
+
(\lambda \nu/h_T [u_h],[v_h])_{\partial T}
\Bigr],
\]
and the discrete divergence-pressure coupling is
\[
b_h(v_h,p_h)=\sum_T\Bigl[
-(\nabla\cdot v_h,p_h)_T + ([v_h\cdot n],\{p_h\})_{\partial T}
\Bigr].
\]
For convection one takes “an upwind- or Temam-modified convection trilinear \(c_h(w;u,v)\) satisfying continuity (2.21) and \(c_h(w;u_h,u_h)\ge 0\).” The global bilinear form is
\[
A_h^w((u_h,p_h),(v_h,q_h))
=
a_h(u_h,v_h)+b_h(v_h,p_h)+b_h(u_h,q_h)+c_h(w;u_h,v_h).
\]
The Trefftz-DG method reads: find \((u_h,p_h)\in X\) such that
\[
A_h^w((u_h,p_h),(v_h,q_h))
+
\langle B_w(u_h,p_h),(r_h,s_h)\rangle
=
(f,v_h)_- + (h_T\nu^{-1/2}f,r_h)
\]
for all \((v_h,q_h)\in T_w\) and \((r_h,s_h)\in Q\) [2606.13229].

Algebraically, this is “an upper-block-triangular system coupling the global Trefftz unknowns with local residual unknowns.” That structure is fundamental: the PDE-residual block is local, whereas the Trefftz block is the part that remains globally assembled. The reduced Oseen discretization is obtained by exploiting exactly this block form.

## 4. Pressure elimination and static condensation

The reduction proceeds from the splitting \(X=T_w\oplus L\) together with invertibility of \(B_w|_L\). On each element \(T\), one solves the local problem for \((u_L,p_L)\in L(T)\) satisfying
\[
B_w(u_L,p_L)=(h_T\nu^{-1/2}f,0)\quad\text{in }Q(T)' .
\]
After this local solve, one substitutes back into the global equation on \(T_w\) only. In the formulation summarized in the paper, this “removes all high-order pressure unknowns and the local flux degrees of freedom” [2606.13229].

In practice, the elementwise matrix representation of \(B_w|_{L(T)}\) is assembled and inverted, “e.g. via QR or SVD,” to build the static-condensation Schur complement. Writing
\[
K_T^{loc}=[B_w|_{L(T)}]^{-1}\circ ,
\]
the condensed contribution to the Trefftz block is
\[
k_w^T(u_T,v_T)
:=
h_{mg}^w((u_T,P_w u_T),(v_T,P_w v_T))
-
(\text{coupling through local solve})(u_T;v_T).
\]
Here \(P_w\) is the discrete pressure lifting \((-\Pi_\nabla A_h^w)\), so that \((u_T,P_w u_T)\in T_w\) [2606.13229].

The reduction therefore has two simultaneous meanings. First, it is a pressure elimination in the algebraic sense: the global system no longer carries the full pressure block. Second, it is a localization mechanism: the residual-driven complement variables are resolved cellwise and never appear as independent global unknowns.

## 5. Velocity-only reduced system and analytical properties

The reduced unknown is
\[
u_h\in T_w^u = \{u_h\in [P^k]^d:\ \exists p_h \text{ s.t. } (u_h,p_h)\in T_w\}.
\]
Imposing discrete \(\nabla\cdot u=0\) against \(P^0\) gives a subspace of \(T_w^u\). On that reduced space, the formulation establishes three key properties: coercivity,
\[
k_w(u_T,u_T)\ge \alpha \|u_T\|_{1,h}^2,
\]
continuity,
\[
k_w(u,v)\lesssim \|u\|_{1,h}\|v\|_{1,h},
\]
and inf-sup on the residual block \(K_w\), expressed as a stable pairing \(\|K_w u_L\|_{'}\gtrsim \|u_L\|_{1,h}\). Consequently, “the fully reduced bilinear form” on the reduced velocity space and local complement is inf-sup stable [2606.13229].

After condensation, the matrix system has the form
\[
S_{ij}u_T^j = F_i,
\]
where
\[
S_{ij} = \sum_T k_w^T(\phi_j,\phi_i)
\]
plus facet terms at \(\partial T\), and
\[
F_i = \sum_T (f,\phi_i)_- - \text{local solve contributions}.
\]
The paper characterizes this as “a global SPD-like system for the DOFs of \(u_T\)” [2606.13229].

The analytical results are stated in DG norms. With \(\|\cdot\|_X\) denoting the DG norm \((\|u\|_{1,h}+\nu^{-1}\|p\|_{0,h})\), Theorem 3.2 gives local invertibility: under \(|w|_{h,d}\le c_0\) small, \(B_w|_L(T)\) is bijective with
\[
\|B_w(\ell)\|_{Q'}\gtrsim \|\ell\|_X \qquad \forall \ell\in L(T).
\]
Theorem 4.3 states that the coupled system is inf-sup stable on the full space \(X=T_w\oplus L\). Theorem 5.1 gives quasi-optimal error in DG norm: if \((u,p)\) is the exact solution with \(u\in H^{k+1}\), \(p\in H^\ell\), and \(u_h\) is the reduced DG solution, then
\[
\|u-u_h\|_{1,h} + \nu^{-1/2}\|p-p_h\|_{0,h}
\lesssim
\inf_{v_h\in [P^k],\ \Pi^0 v_h=\Pi^0 u}\|u-v_h\|_{1,h}
+
\nu^{-1/2}h^\ell \|p\|_{H^\ell}.
\]
Using a BDM-interpolant, the corollary recovers the optimal \(O(h^k)\)–\(O(h^{k+1})\) convergence [2606.13229].

## 6. Assumptions, implementation, and related pressure-eliminating Oseen formulations

The reduced Trefftz-DG construction assumes a shape-regular simplicial mesh \(\mathcal{T}_h\), polynomial degree \(k\ge 1\), and a penalty parameter with \(\lambda\gtrsim k^2\) for coercivity of \(a_h\). The convection form may be “any standard upwind or Temam-modified DG form satisfying (2.21),(2.22),(2.23).” The convective field satisfies
\[
w\in [W^{1,4}(\Omega)]^d
\]
with the resolution bound
\[
|w|_{h,d}:=O(\nu^{-1}h^{1-d/4}\|w\|_{L^4}+\cdots)\le c_0.
\]
All constants in the estimates depend only on shape-regularity, \(k\), and \(\Omega\), not on \(h\) or \(\nu\). The implementation is correspondingly local-to-global: on each \(T\), assemble the small local matrix of \(B_w|_{L(T)}\), compute its inverse or a QR/SVD-based pseudoinverse, build the condensed element stiffness \(k_w^T\) on \(T_w^u\), assemble the global reduced system for \(u_h\in T_w^u\), and recover local complements \((u_L,p_L)\) by back-substitution if needed [2606.13229].

Other recent Oseen discretizations show that pressure elimination is not unique to Trefftz condensation. In a mixed virtual element method for the two-dimensional generalized Oseen problem, the additional unknown is the pseudostress
\[
\Sigma := \nu\nabla u - (u\otimes\beta) - pI,
\]
which yields
\[
-\mathrm{div}\,\Sigma + \kappa u = f,
\qquad
p = -\tfrac12\bigl(\mathrm{tr}\,\Sigma + \mathrm{tr}(u\otimes\beta)\bigr).
\]
The pressure is therefore eliminated from the first-order system and “recovered element-wise” after solving for \((\Sigma_h,u_h)\) [2601.20050]. In the H(div)-conforming finite element method for the steady Oseen equations, the discrete velocity is “exactly divergence-free pointwise,” the pressure enters “only via the constraint and not in the residual or stabilization,” and the velocity estimate contains “no factor \(1/\beta_h\) in front of the data error” [2511.21996]. This suggests that reduced Oseen discretization is a broader design principle: the global approximation may be organized around velocity-dominant unknowns even when the underlying formulations differ substantially.

A common misconception is that a reduced formulation discards pressure. The cited methods indicate a narrower and more precise statement. In the embedded Trefftz-DG setting, the global reduced system is velocity-only, but the pressure is still represented through the discrete pressure lifting and through local back-substitution. In the pseudostress VEM setting, the pressure is recovered from the trace relation. In the H(div) setting, pressure remains part of the mixed problem, but velocity approximation is decoupled from the discrete inf-sup constant in the stated pressure-robust estimate. The common theme is not the absence of pressure, but its displacement from the dominant global approximation mechanism.

Source: https://www.emergentmind.com/topics/reduced-oseen-discretization