Stokes operator

The Stokes operator, named after George Gabriel Stokes, is an unbounded linear operator used in the theory of partial differential equations, specifically in the fields of fluid dynamics and electromagnetics.

Definition

If we define [math]\displaystyle{ P_\sigma }[/math] as the Leray projection onto divergence free vector fields, then the Stokes Operator [math]\displaystyle{ A }[/math] is defined by

[math]\displaystyle{ A:=-P_\sigma\Delta, }[/math]

where [math]\displaystyle{ \Delta\equiv\nabla^2 }[/math] is the Laplacian. Since [math]\displaystyle{ A }[/math] is unbounded, we must also give its domain of definition, which is defined as [math]\displaystyle{ \mathcal{D}(A)=H^2\cap V }[/math], where [math]\displaystyle{ V=\{\vec{u}\in (H^1_0(\Omega))^n|\operatorname{div}\,\vec{u}=0\} }[/math]. Here, [math]\displaystyle{ \Omega }[/math] is a bounded open set in [math]\displaystyle{ \mathbb{R}^n }[/math] (usually n = 2 or 3), [math]\displaystyle{ H^2(\Omega) }[/math] and [math]\displaystyle{ H^1_0(\Omega) }[/math] are the standard Sobolev spaces, and the divergence of [math]\displaystyle{ \vec{u} }[/math] is taken in the distribution sense.

Properties

For a given domain [math]\displaystyle{ \Omega }[/math] which is open, bounded, and has [math]\displaystyle{ C^2 }[/math] boundary, the Stokes operator [math]\displaystyle{ A }[/math] is a self-adjoint positive-definite operator with respect to the [math]\displaystyle{ L^2 }[/math] inner product. It has an orthonormal basis of eigenfunctions [math]\displaystyle{ \{w_k\}_{k=1}^\infty }[/math] corresponding to eigenvalues [math]\displaystyle{ \{\lambda_k\}_{k=1}^\infty }[/math] which satisfy

[math]\displaystyle{ 0\lt \lambda_1\lt \lambda_2\leq\lambda_3\cdots\leq\lambda_k\leq\cdots }[/math]

and [math]\displaystyle{ \lambda_k\rightarrow\infty }[/math] as [math]\displaystyle{ k\rightarrow\infty }[/math]. Note that the smallest eigenvalue is unique and non-zero. These properties allow one to define powers of the Stokes operator. Let [math]\displaystyle{ \alpha\gt 0 }[/math] be a real number. We define [math]\displaystyle{ A^\alpha }[/math] by its action on [math]\displaystyle{ \vec{u}\in \mathcal{D}(A) }[/math]:

[math]\displaystyle{ A^\alpha \vec{u}=\sum_{k=1}^\infty \lambda_k^{\alpha} u_k\vec{w_k} }[/math]

where [math]\displaystyle{ u_k:=(\vec{u},\vec{w_k}) }[/math] and [math]\displaystyle{ (\cdot,\cdot) }[/math] is the [math]\displaystyle{ L^2(\Omega) }[/math] inner product.

The inverse [math]\displaystyle{ A^{-1} }[/math] of the Stokes operator is a bounded, compact, self-adjoint operator in the space [math]\displaystyle{ H:=\{\vec{u}\in (L^2(\Omega))^n| \operatorname{div}\,\vec{u}=0 \text{ and }\gamma(\vec{u})=0\} }[/math], where [math]\displaystyle{ \gamma }[/math] is the trace operator. Furthermore, [math]\displaystyle{ A^{-1}:H\rightarrow V }[/math] is injective.

References

• Temam, Roger (2001), Navier-Stokes Equations: Theory and Numerical Analysis, AMS Chelsea Publishing, ISBN 0-8218-2737-5 
• Constantin, Peter and Foias, Ciprian. Navier-Stokes Equations, University of Chicago Press, (1988)