---
title: Non-Symmetric Nitsche Method
url: https://www.emergentmind.com/topics/non-symmetric-nitsche-method
type: topic
---

# Non-Symmetric Nitsche Method

Searching arXiv for recent and foundational papers on the non-symmetric Nitsche method.
The non-symmetric Nitsche method is a weak imposition technique for boundary or interface constraints in which the consistency term is retained but the adjoint or symmetry term is modified, sign-flipped, or omitted, so that the resulting discrete form is generally not symmetric. Across finite element, unfitted/interface, domain decomposition, contact, and boundary element settings, it appears in penalty-based, super-penalty, and penalty-free forms. Its recurring advantages are the elimination of Lagrange multiplier unknowns and, in several formulations, stability without a large penalty parameter; its classical drawback is loss of adjoint consistency and the resulting difficulty of proving optimal \(L^2\) estimates by standard duality arguments [2203.02603, 1102.0202, 2510.05597].

## 1. Algebraic structure and defining variants

A standard structural decomposition writes a Nitsche discretization as the sum of an interior bilinear form, a boundary or interface consistency term, a symmetry term, and a penalty term. In the abstract variational framework, the symmetric method is obtained by adding both the consistency contribution and its adjoint counterpart, whereas a non-symmetric method is obtained by modifying that last step. The tutorial literature makes this decomposition explicit and notes that the same construction steps can be used to build non-symmetric Nitsche methods [2203.02603].

For the Poisson problem with weak Dirichlet enforcement, a canonical non-symmetric form is
\[
a(u,v) := (\nabla u,\nabla v) - \langle\nabla u\cdot n,v\rangle_{\Gamma} + \langle u,\nabla v\cdot n\rangle_{\Gamma} + c_0\langle h^{-\alpha}u,v\rangle_{\Gamma},
\]
with \(c_0\ge 0\) and \(\alpha\ge 1\). The non-symmetry is concentrated in the boundary pair
\[
-\langle\nabla u\cdot n,v\rangle_{\Gamma} + \langle u,\nabla v\cdot n\rangle_{\Gamma},
\]
which is not invariant under interchange of trial and test functions. In the penalty-free case \(c_0=0\), this reduces to the classical continuous Galerkin non-symmetric Nitsche form studied without any boundary penalty term [2510.05597, 1106.5612].

Different sign conventions coexist in the literature. In the hypersingular boundary element formulation for domain decomposition, the interface bilinear form contains the switch parameter \(\sigma\in\{-1,1\}\); \(\sigma=1\) gives the symmetric variant, while \(\sigma=-1\) gives the non-symmetric or skew-symmetric variant. In that setting, the non-symmetric character is carried by the interface term
\[
\frac{\sigma}{2}\langle [u_h],\, T_1 v_h - T_2 v_h\rangle_{\gamma},
\]
and functions are allowed to be discontinuous across the interface \(\gamma\) [1102.0202].

## 2. Stability mechanisms: coercivity, inf–sup, and penalty-free formulations

Non-symmetric Nitsche methods are associated with two distinct stability mechanisms. In some formulations the bilinear form is directly coercive in a mesh-dependent energy norm; in others, stability is recovered through an inf–sup argument rather than through coercivity.

A basic coercive example is the hypersingular boundary element method with \(\sigma=-1\). There, discrete ellipticity holds for any \(\nu>0\), whereas the symmetric case \(\sigma=1\) requires \(\nu\ge C_1|\log \underline h|^3\). The same contrast appears in fitted and unfitted interface formulations: the non-symmetric choice trades symmetry for robustness with respect to the penalty parameter. In the unfitted Poisson interface method, the mesh-dependent norm
\[
\interleave v\interleave^2
:= \mu_1\|\nabla v_1\|_{\Omega_1}^2
 + \mu_2\|\nabla v_2\|_{\Omega_2}^2
 + c_0\|h^{-1/2}[\![v]\!]\|_\Gamma^2
\]
satisfies \(a(v_h,v_h)=\interleave v_h\interleave^2\), so the discrete problem is coercive on \(V_h\). The coefficient-weighted averages
\[
w_1=\frac{\mu_2}{\mu_1+\mu_2},\qquad
w_2=\frac{\mu_1}{\mu_1+\mu_2},\qquad
c_0=\frac{2\mu_1\mu_2}{\mu_1+\mu_2}
\]
are used precisely because they yield estimates independent of the coefficient jump [1102.0202, 2510.12151].

The penalty-free branch follows a different route. For the continuous Galerkin method without boundary penalty, testing with \(v_h=u_h\) only gives control of the interior gradient. The key observation is therefore that boundary control can be recovered through an inf–sup condition rather than via a coercivity estimate. In the Poisson case this is achieved by constructing special patch functions \(\varphi_r\) with prescribed average normal flux on boundary patches \(F_j\), and then testing with \(w_h=v_h+\eta\varphi_r\). Closely related patch constructions appear in penalty-free domain decomposition and in compressible and incompressible elasticity, where the non-symmetric form is shown to satisfy a Banach–Nečas–Babuška type stability estimate even though no term of the form \(\gamma h^{-1}\langle u_h,v_h\rangle\) is present on the boundary or interface [1106.5612, 1509.00542, 1407.2229].

This suggests a genuine bifurcation inside the subject. One branch keeps a penalty and exploits direct coercivity; the other removes the penalty and replaces coercivity by a carefully designed inf–sup argument. The non-symmetric method occurs in both branches.

## 3. Consistency, adjoint consistency, and error estimation

Primal consistency is a defining feature. In the hypersingular interface problem, the exact solution satisfies \([u]=0\) and the interface condition \(T_1u+T_2u=0\) on \(\gamma\), so the jump and interface terms are not required for consistency. The same pattern holds for the Poisson and elasticity formulations: the exact solution satisfies the discrete equations when inserted into the Nitsche bilinear form. What is typically lost in the non-symmetric variant is adjoint consistency [1102.0202, 1106.5612].

That loss dominated the classical \(L^2\) theory. For the penalty-free non-symmetric Poisson method, the proved energy estimate is optimal,
\[
\|u-u_h\|_{1,h}\le C h^k |u|_{H^{k+1}(\Omega)},
\]
but the \(L^2\) estimate is only
\[
\|u-u_h\|_{L^2(\Omega)}\le C h^{k+1/2}|u|_{H^{k+1}(\Omega)}.
\]
The same half-order loss appears in penalty-free domain decomposition for discontinuous material parameters and in penalty-free compressible elasticity. In the boundary element hypersingular setting, the main estimate is “almost quasi-optimal”:
\[
|u-u_h|_{H^{1/2}(\mathcal{T})}
\le C_2 |\log \underline h|^{3/2} h^{r-1/2}\|u\|_{H^r(\Gamma)},
\]
with the logarithmic perturbation attributed to the Nitsche coupling [1106.5612, 1509.00542, 1407.2229, 1102.0202].

Recent analyses changed this picture by modifying the dual problem instead of the primal method. For weak Dirichlet imposition, the adjoint bilinear form
\[
a^\star(u,v)
:= (\nabla u,\nabla v)
 + \langle\nabla u\cdot n,v\rangle_{\Gamma}
 - \langle u,\nabla v\cdot n\rangle_{\Gamma}
 + c_0\langle h^{-\alpha}u,v\rangle_{\Gamma}
\]
satisfies \(a^\star(u,v)=a(v,u)\). By defining the dual problem with \(a^\star\) rather than with the standard adjoint PDE, one restores adjoint consistency at the variational level. This yields optimal \(L^2\) convergence,
\[
\|u-u_h\|_0 \le C h^{k+1}\big(|u|_{k+1}+|u|_{k+1/2,\Gamma}\big),
\]
for super-penalty, standard-penalty, and penalty-free variants, including the practically important case \(\alpha=1,\ c_0>0\) on general shape-regular meshes. An analogous tailored dual problem yields optimal \(L^2\) estimates in unfitted interface FEM, together with the optimal energy estimate
\[
\interleave u-u_h\interleave
\le C h^k\big(\mu_1^{1/2}|u_1|_{k+1,\Omega_1}
+\mu_2^{1/2}|u_2|_{k+1,\Omega_2}\big).
\]
The modern theory therefore no longer treats suboptimal \(L^2\) convergence as intrinsic to non-symmetry; rather, it identifies the standard dual problem as the source of the gap [2510.05597, 2510.12151].

## 4. Interfaces, non-matching meshes, and unfitted discretizations

The non-symmetric Nitsche method is especially prominent in interface and domain decomposition problems because it weakly enforces continuity without introducing multiplier spaces. In the hypersingular boundary element setting, the surface is decomposed into two polygonal subdomains \(\Gamma_1\) and \(\Gamma_2\) with interface \(\gamma=\overline{\Gamma}_1\cap\overline{\Gamma}_2\). The discrete space
\[
X_h = X_{h,1}\times X_{h,2}
\]
contains functions that satisfy the homogeneous boundary condition on \(\partial\Gamma\) but are in general discontinuous across \(\gamma\). The Nitsche coupling acts directly on the jump \([v_h]\) and on interface operators \(T_i\), so non-matching meshes are allowed and no Lagrange multiplier is required. The paper explicitly lists the resulting advantages: no additional unknown on the interface, no inf–sup condition among interface spaces, elliptic discrete problems, and the possibility of standard linear solvers [1102.0202].

The same principle appears in fitted and unfitted Poisson interface problems with discontinuous coefficients. In the fitted penalty-free method, the interface bilinear form is
\[
A_h(u_h,v_h)
= \sum_{i=1}^2 (\mu_i\nabla u_h^i,\nabla v_h^i)_{\Omega_i}
- \langle \{\mu\nabla u_h\cdot n\},\llbracket v_h\rrbracket\rangle_\Gamma
+ \langle \{\mu\nabla v_h\cdot n\},\llbracket u_h\rrbracket\rangle_\Gamma,
\]
again with no interface penalty term. In the unfitted case the same interface form is combined with a ghost penalty acting on faces near the cut interface, not on \(\Gamma\) itself. This stabilization controls small cut elements while preserving the penalty-free character of the Nitsche coupling on the interface [1509.00542].

Unfitted interface FEM with piecewise polynomial spaces on \(\mathcal{T}_h^1\) and \(\mathcal{T}_h^2\) sharpens this picture. There the non-symmetric interface bilinear form is
\[
\begin{aligned}
a(u,v)
&:= \mu_1(\nabla u_1,\nabla v_1)_{\Omega_1}
 + \mu_2(\nabla u_2,\nabla v_2)_{\Omega_2} \\
&\quad - \langle\{\!\!\{\mu\nabla u\cdot\mathbf n\}\!\!\},[\![v]\!]\rangle_\Gamma
 + \langle[\![u]\!],\{\!\!\{\mu\nabla v\cdot\mathbf n\}\!\!\}\rangle_\Gamma
 + \langle c_0 h^{-1}[\![u]\!],[\![v]\!]\rangle_\Gamma .
\end{aligned}
\]
The special weights \(w_1,w_2,c_0\) are introduced to make the constants independent of the contrast, and a penalty-free variant adds a ghost penalty
\[
s_h(u_h,v_h)
= \gamma_g\sum_{i=1}^2\sum_{E\in\mathfrak E_\Gamma^i}\sum_{l=1}^k
\langle \mu_i h^{2l-1} [\![D_{n_E}^l u_h^i]\!], [\![D_{n_E}^l v_h^i]\!]\rangle_E .
\]
The interface literature thus uses non-symmetric Nitsche both with and without explicit penalty on the jump, and in both fitted and cut settings [2510.12151].

## 5. Extensions in mechanics and fluid dynamics

In contact mechanics, the non-symmetric Nitsche method is not limited to equality constraints. For the Signorini problem, the contact conditions are rewritten as
\[
u = -[P_\gamma(u)]_+,\qquad P_\gamma(u):=\gamma\partial_n u-u,
\]
and the resulting penalty-free non-symmetric formulation is
\[
(\nabla u,\nabla v)_\Omega
-\langle \partial_n u,v\rangle_{\Gamma_C}
+\langle u,\partial_n v\rangle_{\Gamma_C}
+\langle [P_\gamma(u)]_+,\partial_n v\rangle_{\Gamma_C}
=(f,v)_\Omega .
\]
The method uses nonconforming piecewise affine elements and achieves optimal error estimates in the energy norm [1609.03745].

Linear elasticity provides another major branch. For compressible elasticity, the nonsymmetric penalty-free boundary form is
\[
A_h(\boldsymbol u_h,\boldsymbol v_h)
= a(\boldsymbol u_h,\boldsymbol v_h)
- b(\boldsymbol u_h,\boldsymbol v_h)
+ b(\boldsymbol v_h,\boldsymbol u_h),
\]
with \(b(\boldsymbol u_h,\boldsymbol v_h)\) defined through the normal traction of the stress tensor. For incompressible elasticity, the same nonsymmetric Nitsche boundary treatment is coupled with Galerkin least-squares pressure stabilization, producing locking-free convergence in the mixed formulation [1407.2229].

In incompressible flow, a non-symmetric Nitsche method has been analyzed for the stationary Navier–Stokes equations with slip boundary conditions. There the normal component of the velocity on \(\Gamma_{\mathrm{Nav}}\) is enforced weakly through boundary consistency and penalty terms involving \(\mathbf n^T(\nu D(\mathbf u_h)-p_hI)\mathbf n\), while the tangential Navier law remains as a Robin-type term. By contrast, the later Stokes–Poisson–Boltzmann work with Navier slip explicitly uses the symmetric variant of Nitsche’s method, which isolates the non-symmetric method as one choice inside a broader weak-imposition framework rather than as the universal formulation [2307.03589, 2604.12396].

A general mechanics formulation pushes the same idea further: every Nitsche method is viewed as a stabilized Lagrange multiplier method in which the multiplier has been eliminated. In that framework, equality constraints arise by dropping the positive-part operator in the reduced energy, while inequality constraints retain the term
\[
\int_\Gamma \frac{\gamma(h,\kappa)}{2}
\Big(\lambda(u_h)-\tfrac{1}{\gamma(h,\kappa)}\beta(u_h)\Big)_+^2\,ds
-
\int_\Gamma \frac{\gamma(h,\kappa)}{2}\lambda(u_h)^2\,ds .
\]
This viewpoint does not privilege the symmetric method: it identifies symmetric and non-symmetric Nitsche variants as different eliminations of a stabilized saddle-point formulation [2603.05008].

## 6. Numerical behavior, parameter choices, and solver implications

The numerical record shows a persistent tension between old theory and observed performance. In the penalty-free Poisson method, the \(H^1\)-error is essentially insensitive to the penalty parameter \(\gamma\), and the \(L^2\)-error is often slightly smaller for moderate positive \(\gamma\), even though the theory with \(\gamma=0\) gives only \(h^{k+1/2}\). In fitted and unfitted domain decomposition, the proved rate is again \(h^{k+1/2}\) in \(L^2\), but the experiments with \(k=1\) show approximately \(\mathcal O(h^2)\). The penalty-free elasticity paper reports the same phenomenon: optimal \(H^1\) convergence in theory and numerically optimal \(L^2\) behavior in the tests [1106.5612, 1509.00542, 1407.2229].

Parameter sensitivity depends strongly on the chosen variant. In the hypersingular boundary element method, the skew-symmetric case \(\sigma=-1\) converges for all tested \(\nu>0\), whereas the symmetric case can fail to converge when \(\nu\) is too small and requires \(\nu\) scaling at least like \(|\log \underline h|^3\) in theory. In unfitted interface FEM, the standard non-symmetric penalty term \(c_0h^{-1}[\![u_h]\!]\) is not required to be “very large,” because stability already follows from the weighted flux structure and the choice of \(c_0\). In penalty-free variants, the practical stabilization task shifts from tuning a boundary penalty to designing the correct inf–sup or ghost-penalty mechanism [1102.0202, 2510.12151].

The algebraic consequences are equally consistent across the literature. Symmetric Nitsche preserves symmetry of the discrete operator for symmetric model problems; non-symmetric Nitsche yields non-symmetric matrices. In the hypersingular boundary element paper, the symmetric case admits conjugate gradients once the penalty is large enough, whereas the non-symmetric case calls for Krylov methods for non-symmetric systems such as GMRES or BiCGSTAB. The penalty-free Poisson formulation likewise produces a sparse non-symmetric matrix of the form \(A=K-B+B^T\), with no extra interface or multiplier unknowns [1102.0202, 1106.5612].

Recent numerical experiments finally align the \(L^2\) theory with longstanding computations. Three-dimensional tests for the weak Dirichlet problem confirm \(H^1\)-error \(\mathcal O(h^k)\) and \(L^2\)-error \(\mathcal O(h^{k+1})\) for standard penalty, super-penalty, and penalty-free non-symmetric Nitsche variants. In the unfitted interface setting, the tailored-dual analysis is also reported to match the observed optimal \(L^2\) rates. This suggests that the modern interpretation of non-symmetric Nitsche is no longer dominated by its historical half-order \(L^2\) loss, but by a more precise distinction between primal consistency, adjoint consistency, and the choice of dual problem in the analysis [2510.05597, 2510.12151].

Source: https://www.emergentmind.com/topics/non-symmetric-nitsche-method