In numerical partial differential equations, the Ladyzhenskaya–Babuška–Brezzi (LBB) condition is a sufficient condition for a saddle point problem to have a unique solution that depends continuously on the input data. Saddle point problems arise in the discretization of Stokes flow and in the mixed finite element discretization of Poisson's equation. For positive-definite problems, like the unmixed formulation of the Poisson equation, most discretization schemes will converge to the true solution in the limit as the mesh is refined. For saddle point problems, however, many discretizations are unstable, giving rise to artifacts such as spurious oscillations. The LBB condition gives criteria for when a discretization of a saddle point problem is stable. The condition is variously referred to as the LBB condition, the Babuška–Brezzi condition, or the "inf-sup" condition.
Saddle point problems The abstract form of a saddle point problem can be expressed in terms of Hilbert spaces and bilinear forms. Let V {\displaystyle V} and Q {\displaystyle Q} be Hilbert spaces, and let a : V × V → R {\displaystyle a:V\times V\to \mathbb {R} } , b : V × Q → R {\displaystyle b:V\times Q\to \mathbb {R} } be bilinear forms. Let f ∈ V ∗ {\displaystyle f\in V^{*}} , g ∈ Q ∗ {\displaystyle g\in Q^{*}} where V ∗ {\displaystyle V^{*}} , Q ∗ {\displaystyle Q^{*}} are the dual spaces. The saddle-point problem for the pair a {\displaystyle a} , b {\displaystyle b} is to find a pair of fields u {\displaystyle u} in V {\displaystyle V} , p {\displaystyle p} in Q {\displaystyle Q} such that, for all v {\displaystyle v} in V {\displaystyle V} and q {\displaystyle q} in Q {\displaystyle Q} ,
a ( u , v ) + b ( v , p ) = ⟨ f , v ⟩ b ( u , q ) = ⟨ g , q ⟩ . {\displaystyle {\begin{aligned}a(u,v)+b(v,p)&=\langle f,v\rangle \\b(u,q)&=\langle g,q\rangle .\end{aligned}}}
For example, for the Stokes equations on a d {\displaystyle d} -dimensional domain Ω {\displaystyle \Omega } , the fields are the velocity u {\displaystyle u} and pressure p {\displaystyle p} , which live in respectively the Sobolev space H 1 ( Ω ) d {\displaystyle H^{1}(\Omega )^{d}} and the Lebesgue space L 2 ( Ω ) {\displaystyle L^{2}(\Omega )} . The bilinear forms for this problem are
a ( u , v ) = ∫ Ω μ ∇ u : ∇ v d x b ( u , q ) = ∫ Ω ( ∇ ⋅ u ) q d x , {\displaystyle {\begin{aligned}a(u,v)&=\int _{\Omega }\mu \nabla u:\nabla v\,dx\\b(u,q)&=\int _{\Omega }(\nabla \cdot u)q\,dx,\end{aligned}}}
… excerpt ends here. Continue reading the full article.
