Reference articles on history, science, culture and more
Encyclopedia

Trace operator

Boundary condition for generalized functions

In mathematical analysis, the trace operator extends the notion of the restriction of a function to the boundary of its domain to "generalized" functions in a Sobolev space. This is particularly important for the study of partial differential equations with prescribed boundary conditions (boundary value problems), where weak solutions may not be regular enough to satisfy the boundary conditions in the classical sense of functions.

01Motivation

On a bounded, smooth domain {\textstyle \Omega \subset \mathbb {R} ^{n}, consider the problem of solving Poisson's equation with inhomogeneous Dirichlet boundary conditions:

{\begin{alignedat}{2}-\Delta u&=f&\quad &{\text{in }}\Omega ,\\u&=g&&{\text{on }}\partial \Omega \end{alignedat}}

with given functions {\textstyle f and {\textstyle g with regularity discussed in the application section below. The weak solution {\textstyle u\in H^{1}(\Omega ) of this equation must satisfy

\int _{\Omega }\nabla u\cdot \nabla \varphi \,\mathrm {d} x=\int _{\Omega }f\varphi \,\mathrm {d} x for all {\textstyle \varphi \in H_{0}^{1}(\Omega ).

The {\textstyle H^{1}(\Omega )-regularity of {\textstyle u is sufficient for the well-definedness of this integral equation. It is not apparent, however, in which sense {\textstyle u can satisfy the boundary condition {\textstyle u=g on {\textstyle \partial \Omega: by definition, {\textstyle u\in H^{1}(\Omega )\subset L^{2}(\Omega ) is an equivalence class of functions which can have arbitrary values on {\textstyle \partial \Omega since this is a null set with respect to the n-dimensional Lebesgue measure.

If {\textstyle \Omega \subset \mathbb {R} ^{1} there holds {\textstyle H^{1}(\Omega )\hookrightarrow C^{0}({\bar {\Omega }}) by Sobolev's embedding theorem, such that {\textstyle u can satisfy the boundary condition in the classical sense, i.e. the restriction of {\textstyle u to {\textstyle \partial \Omega agrees with the function {\textstyle g (more precisely: there exists a representative of {\textstyle u in {\textstyle C({\bar {\Omega }}) with this property). For {\textstyle \Omega \subset \mathbb {R} ^{n} with {\textstyle n>1 such an embedding does not exist and the trace operator {\textstyle T presented here must be used to give meaning to {\textstyle u|_{\partial \Omega }. Then {\textstyle u\in H^{1}(\Omega ) with {\textstyle Tu=g is called a weak solution to the boundary value problem if the integral equation above is satisfied. For the definition of the trace operator to be reasonable, there must hold {\textstyle Tu=u|_{\partial \Omega } for sufficiently regular {\textstyle u.

A function defined on a rectangle (top figure, in red), and its trace (bottom figure, in red).
A function defined on a rectangle (top figure, in red), and its trace (bottom figure, in red).

02Trace theorem

The trace operator can be defined for functions in the Sobolev spaces {\textstyle W^{1,p}(\Omega ) with {\textstyle 1\leq p<\infty, see the section below for possible extensions of the trace to other spaces. Let {\textstyle \Omega \subset \mathbb {R} ^{n} for {\textstyle n\in \mathbb {N} be a bounded domain with Lipschitz boundary. Then there exists a bounded linear trace operator

T\colon W^{1,p}(\Omega )\to L^{p}(\partial \Omega )

such that {\textstyle T extends the classical trace, i.e.

Tu=u|_{\partial \Omega } for all {\textstyle u\in W^{1,p}(\Omega )\cap C({\bar {\Omega }}).

The continuity of {\textstyle T implies that

\|Tu\|_{L^{p}(\partial \Omega )}\leq C\|u\|_{W^{1,p}(\Omega )} for all {\textstyle u\in W^{1,p}(\Omega )

with constant only depending on {\textstyle p and {\textstyle \Omega. The function {\textstyle Tu is called trace of {\textstyle u and is often simply denoted by {\textstyle u|_{\partial \Omega }. Other common symbols for {\textstyle T include {\textstyle tr and {\textstyle \gamma.

Construction

This paragraph follows Evans, where more details can be found, and assumes that {\textstyle \Omega has a {\textstyle C^{1}-boundary. A proof (of a stronger version) of the trace theorem for Lipschitz domains can be found in Gagliardo. On a {\textstyle C^{1}-domain, the trace operator can be defined as continuous linear extension of the operator

T:C^{\infty }({\bar {\Omega }})\to L^{p}(\partial \Omega )

to the space {\textstyle W^{1,p}(\Omega ). By density of {\textstyle C^{\infty }({\bar {\Omega }}) in {\textstyle W^{1,p}(\Omega ) such an extension is possible if {\textstyle T is continuous with respect to the {\textstyle W^{1,p}(\Omega )-norm. The proof of this, i.e. that there exists {\textstyle C>0 (depending on {\textstyle \Omega and {\textstyle p) such that

\|Tu\|_{L^{p}(\partial \Omega )}\leq C\|u\|_{W^{1,p}(\Omega )} for all u\in C^{\infty }({\bar {\Omega }}),

is the central ingredient in the construction of the trace operator. A local variant of this estimate for {\textstyle C^{1}({\bar {\Omega }})-functions is first proven for a locally flat boundary using the divergence theorem. By transformation, a general {\textstyle C^{1}-boundary can be locally straightened to reduce to this case, where the {\textstyle C^{1}-regularity of the transformation requires that the local estimate holds for {\textstyle C^{1}({\bar {\Omega }})-functions.

With this continuity of the trace operator in {\textstyle C^{\infty }({\bar {\Omega }}) an extension to {\textstyle W^{1,p}(\Omega ) exists by abstract arguments and {\textstyle Tu for {\textstyle u\in W^{1,p}(\Omega ) can be characterized as follows. Let {\textstyle u_{k}\in C^{\infty }({\bar {\Omega }}) be a sequence approximating {\textstyle u\in W^{1,p}(\Omega ) by density. By the proven continuity of {\textstyle T in {\textstyle C^{\infty }({\bar {\Omega }}) the sequence {\textstyle u_{k}|_{\partial \Omega } is a Cauchy sequence in {\textstyle L^{p}(\partial \Omega ) and {\textstyle Tu=\lim _{k\to \infty }u_{k}|_{\partial \Omega } with limit taken in {\textstyle L^{p}(\partial \Omega ).

The extension property {\textstyle Tu=u|_{\partial \Omega } holds for {\textstyle u\in C^{\infty }({\bar {\Omega }}) by construction, but for any {\textstyle u\in W^{1,p}(\Omega )\cap C({\bar {\Omega }}) there exists a sequence {\textstyle u_{k}\in C^{\infty }({\bar {\Omega }}) which converges uniformly on {\textstyle {\bar {\Omega }} to {\textstyle u, verifying the extension property on the larger set {\textstyle W^{1,p}(\Omega )\cap C({\bar {\Omega }}).

The case p = ∞

If {\textstyle \Omega is bounded and has a {\textstyle C^{1}-boundary then by Morrey's inequality there exists a continuous embedding {\textstyle W^{1,\infty }(\Omega )\hookrightarrow C^{0,1}(\Omega ), where {\textstyle C^{0,1}(\Omega ) denotes the space of Lipschitz continuous functions. In particular, any function {\textstyle u\in W^{1,\infty }(\Omega ) has a classical trace {\textstyle u|_{\partial \Omega }\in C(\partial \Omega ) and there holds

\|u|_{\partial \Omega }\|_{C(\partial \Omega )}\leq \|u\|_{C^{0,1}(\Omega )}\leq C\|u\|_{W^{1,\infty }(\Omega )}.

03Functions with trace zero

The Sobolev spaces {\textstyle W_{0}^{1,p}(\Omega ) for {\textstyle 1\leq p<\infty are defined as the closure of the set of compactly supported test functions {\textstyle C_{c}^{\infty }(\Omega ) with respect to the {\textstyle W^{1,p}(\Omega )-norm. The following alternative characterization holds:

W_{0}^{1,p}(\Omega )=\{u\in W^{1,p}(\Omega )\mid Tu=0\}=\ker(T\colon W^{1,p}(\Omega )\to L^{p}(\partial \Omega )),

where {\textstyle \ker(T) is the kernel of {\textstyle T, i.e. {\textstyle W_{0}^{1,p}(\Omega ) is the subspace of functions in {\textstyle W^{1,p}(\Omega ) with trace zero.

04Image of the trace operator

For p > 1

The trace operator is not surjective onto {\textstyle L^{p}(\partial \Omega ) if {\textstyle p>1, i.e. not every function in {\textstyle L^{p}(\partial \Omega ) is the trace of a function in {\textstyle W^{1,p}(\Omega ). As elaborated below the image consists of functions which satisfy an {\textstyle L^{p}-version of Hölder continuity.

Abstract characterization

An abstract characterization of the image of {\textstyle T can be derived as follows. By the isomorphism theorems there holds

T(W^{1,p}(\Omega ))\cong W^{1,p}(\Omega )/\ker(T\colon W^{1,p}(\Omega )\to L^{p}(\partial \Omega ))=W^{1,p}(\Omega )/W_{0}^{1,p}(\Omega )

where {\textstyle X/N denotes the quotient space of the Banach space {\textstyle X by the subspace {\textstyle N\subset X and the last identity follows from the characterization of {\textstyle W_{0}^{1,p}(\Omega ) from above. Equipping the quotient space with the quotient norm defined by

\|u\|_{W^{1,p}(\Omega )/W_{0}^{1,p}(\Omega )}=\inf _{u_{0}\in W_{0}^{1,p}(\Omega )}\|u-u_{0}\|_{W^{1,p}(\Omega )}

the trace operator {\textstyle T is then a surjective, bounded linear operator

T\colon W^{1,p}(\Omega )\to W^{1,p}(\Omega )/W_{0}^{1,p}(\Omega ).

Characterization using Sobolev-Slobodeckij spaces

A more concrete representation of the image of {\textstyle T can be given using Sobolev-Slobodeckij spaces which generalize the concept of Hölder continuous functions to the {\textstyle L^{p}-setting. Since {\textstyle \partial \Omega is a (n-1)-dimensional Lipschitz manifold embedded into {\textstyle \mathbb {R} ^{n} an explicit characterization of these spaces is technically involved. For simplicity consider first a planar domain {\textstyle \Omega '\subset \mathbb {R} ^{n-1}. For {\textstyle v\in L^{p}(\Omega ') define the (possibly infinite) norm

\|v\|_{W^{1-1/p,p}(\Omega ')}=\left(\|v\|_{L^{p}(\Omega ')}^{p}+\int _{\Omega '\times \Omega '}{\frac {|v(x)-v(y)|^{p}}{|x-y|^{(1-1/p)p+(n-1)}}}\,\mathrm {d} (x,y)\right)^{1/p}

which generalizes the Hölder condition {\textstyle |v(x)-v(y)|\leq C|x-y|^{1-1/p}. Then

W^{1-1/p,p}(\Omega ')=\left\{v\in L^{p}(\Omega ')\;\mid \;\|v\|_{W^{1-1/p,p}(\Omega ')}<\infty \right\}

equipped with the previous norm is a Banach space (a general definition of {\textstyle W^{s,p}(\Omega ') for non-integer {\textstyle s>0 can be found in the article for Sobolev-Slobodeckij spaces). For the (n-1)-dimensional Lipschitz manifold {\textstyle \partial \Omega define {\textstyle W^{1-1/p,p}(\partial \Omega ) by locally straightening {\textstyle \partial \Omega and proceeding as in the definition of {\textstyle W^{1-1/p,p}(\Omega ').

The space {\textstyle W^{1-1/p,p}(\partial \Omega ) can then be identified as the image of the trace operator and there holds that

T\colon W^{1,p}(\Omega )\to W^{1-1/p,p}(\partial \Omega )

is a surjective, bounded linear operator.

For p = 1

For {\textstyle p=1 the image of the trace operator is {\textstyle L^{1}(\partial \Omega ) and there holds that

T\colon W^{1,1}(\Omega )\to L^{1}(\partial \Omega )

is a surjective, bounded linear operator.

05Right-inverse: trace extension operator

The trace operator is not injective since multiple functions in {\textstyle W^{1,p}(\Omega ) can have the same trace (or equivalently, {\textstyle W_{0}^{1,p}(\Omega )\neq 0). The trace operator has however a well-behaved right-inverse, which extends a function defined on the boundary to the whole domain. Specifically, for {\textstyle 1<p<\infty there exists a bounded, linear trace extension operator

E\colon W^{1-1/p,p}(\partial \Omega )\to W^{1,p}(\Omega ),

using the Sobolev-Slobodeckij characterization of the trace operator's image from the previous section, such that

T(Ev)=v for all {\textstyle v\in W^{1-1/p,p}(\partial \Omega )

and, by continuity, there exists {\textstyle C>0 with

\|Ev\|_{W^{1,p}(\Omega )}\leq C\|v\|_{W^{1-1/p,p}(\partial \Omega )}.

Notable is not the mere existence but the linearity and continuity of the right inverse. This trace extension operator must not be confused with the whole-space extension operators {\textstyle W^{1,p}(\Omega )\to W^{1,p}(\mathbb {R} ^{n}) which play a fundamental role in the theory of Sobolev spaces.

06Extension to other spaces

Higher derivatives

Many of the previous results can be extended to {\textstyle W^{m,p}(\Omega ) with higher differentiability {\textstyle m=2,3,\ldots if the domain is sufficiently regular. Let {\textstyle N denote the exterior unit normal field on {\textstyle \partial \Omega. Since {\textstyle u|_{\partial \Omega } can encode differentiability properties in tangential direction only the normal derivative {\textstyle \partial _{N}u|_{\partial \Omega } is of additional interest for the trace theory for {\textstyle m=2. Similar arguments apply to higher-order derivatives for {\textstyle m>2.

Let {\textstyle 1<p<\infty and {\textstyle \Omega \subset \mathbb {R} ^{n} be a bounded domain with {\textstyle C^{m,1}-boundary. Then there exists a surjective, bounded linear higher-order trace operator

T_{m}\colon W^{m,p}(\Omega )\to \prod _{l=0}^{m-1}W^{m-l-1/p,p}(\partial \Omega )

with Sobolev-Slobodeckij spaces {\textstyle W^{s,p}(\partial \Omega ) for non-integer {\textstyle s>0 defined on {\textstyle \partial \Omega through transformation to the planar case {\textstyle W^{s,p}(\Omega ') for {\textstyle \Omega '\subset \mathbb {R} ^{n-1}, whose definition is elaborated in the article on Sobolev-Slobodeckij spaces. The operator {\textstyle T_{m} extends the classical normal traces in the sense that

T_{m}u=\left(u|_{\partial \Omega },\partial _{N}u|_{\partial \Omega },\ldots ,\partial _{N}^{m-1}u|_{\partial \Omega }\right) for all {\textstyle u\in W^{m,p}(\Omega )\cap C^{m-1}({\bar {\Omega }}).

Furthermore, there exists a bounded, linear right-inverse of {\textstyle T_{m}, a higher-order trace extension operator

E_{m}\colon \prod _{l=0}^{m-1}W^{m-l-1/p,p}(\partial \Omega )\to W^{m,p}(\Omega ).

Finally, the spaces {\textstyle W_{0}^{m,p}(\Omega ), the completion of {\textstyle C_{c}^{\infty }(\Omega ) in the {\textstyle W^{m,p}(\Omega )-norm, can be characterized as the kernel of {\textstyle T_{m}, i.e.

W_{0}^{m,p}(\Omega )=\{u\in W^{m,p}(\Omega )\mid T_{m}u=0\}.

Less regular spaces

No trace in Lp

There is no sensible extension of the concept of traces to {\textstyle L^{p}(\Omega ) for {\textstyle 1\leq p<\infty since any bounded linear operator which extends the classical trace must be zero on the space of test functions {\textstyle C_{c}^{\infty }(\Omega ), which is a dense subset of {\textstyle L^{p}(\Omega ), implying that such an operator would be zero everywhere.

Generalized normal trace

Let {\textstyle \operatorname {div} v denote the distributional divergence of a vector field {\textstyle v. For {\textstyle 1<p<\infty and bounded Lipschitz domain {\textstyle \Omega \subset \mathbb {R} ^{n} define

E_{p}(\Omega )=\{v\in (L^{p}(\Omega ))^{n}\mid \operatorname {div} v\in L^{p}(\Omega )\}

which is a Banach space with norm

\|v\|_{E_{p}(\Omega )}=\left(\|v\|_{L^{p}(\Omega )}^{p}+\|\operatorname {div} v\|_{L^{p}(\Omega )}^{p}\right)^{1/p}.

Let {\textstyle N denote the exterior unit normal field on {\textstyle \partial \Omega. Then there exists a bounded linear operator

T_{N}\colon E_{p}(\Omega )\to (W^{1-1/q,q}(\partial \Omega ))',

where {\textstyle q=p/(p-1) is the conjugate exponent to {\textstyle p and {\textstyle X' denotes the continuous dual space to a Banach space {\textstyle X, such that {\textstyle T_{N} extends the normal trace {\textstyle (v\cdot N)|_{\partial \Omega } for {\textstyle v\in (C^{\infty }({\bar {\Omega }}))^{n} in the sense that

T_{N}v={\bigl \{}\varphi \in W^{1-1/q,q}(\partial \Omega )\mapsto \int _{\partial \Omega }\varphi v\cdot N\,\mathrm {d} S{\bigr \}}.

The value of the normal trace operator {\textstyle (T_{N}v)(\varphi ) for {\textstyle \varphi \in W^{1-1/q,q}(\partial \Omega ) is defined by application of the divergence theorem to the vector field {\textstyle w=E\varphi \,v where {\textstyle E is the trace extension operator from above.

Application. Any weak solution {\textstyle u\in H^{1}(\Omega ) to {\textstyle -\Delta u=f\in L^{2}(\Omega ) in a bounded Lipschitz domain {\textstyle \Omega \subset \mathbb {R} ^{n} has a normal derivative in the sense of {\textstyle T_{N}\nabla u\in (W^{1/2,2}(\partial \Omega ))^{*}. This follows as {\textstyle \nabla u\in E_{2}(\Omega ) since {\textstyle \nabla u\in L^{2}(\Omega ) and {\textstyle \operatorname {div} (\nabla u)=\Delta u=-f\in L^{2}(\Omega ). This result is notable since in Lipschitz domains in general {\textstyle u\not \in H^{2}(\Omega ), such that {\textstyle \nabla u may not lie in the domain of the trace operator {\textstyle T.

07Application

The theorems presented above allow a closer investigation of the boundary value problem

{\begin{alignedat}{2}-\Delta u&=f&\quad &{\text{in }}\Omega ,\\u&=g&&{\text{on }}\partial \Omega \end{alignedat}}

on a Lipschitz domain {\textstyle \Omega \subset \mathbb {R} ^{n} from the motivation. Since only the Hilbert space case {\textstyle p=2 is investigated here, the notation {\textstyle H^{1}(\Omega ) is used to denote {\textstyle W^{1,2}(\Omega ) etc. As stated in the motivation, a weak solution {\textstyle u\in H^{1}(\Omega ) to this equation must satisfy {\textstyle Tu=g and

\int _{\Omega }\nabla u\cdot \nabla \varphi \,\mathrm {d} x=\int _{\Omega }f\varphi \,\mathrm {d} x for all {\textstyle \varphi \in H_{0}^{1}(\Omega ),

where the right-hand side must be interpreted for {\textstyle f\in H^{-1}(\Omega )=(H_{0}^{1}(\Omega ))' as a duality product with the value {\textstyle f(\varphi ).

Existence and uniqueness of weak solutions

The characterization of the range of {\textstyle T implies that for {\textstyle Tu=g to hold the regularity {\textstyle g\in H^{1/2}(\partial \Omega ) is necessary. This regularity is also sufficient for the existence of a weak solution, which can be seen as follows. By the trace extension theorem there exists {\textstyle Eg\in H^{1}(\Omega ) such that {\textstyle T(Eg)=g. Defining {\textstyle u_{0} by {\textstyle u_{0}=u-Eg we have that {\textstyle Tu_{0}=Tu-T(Eg)=0 and thus {\textstyle u_{0}\in H_{0}^{1}(\Omega ) by the characterization of {\textstyle H_{0}^{1}(\Omega ) as space of trace zero. The function {\textstyle u_{0}\in H_{0}^{1}(\Omega ) then satisfies the integral equation

\int _{\Omega }\nabla u_{0}\cdot \nabla \varphi \,\mathrm {d} x=\int _{\Omega }\nabla (u-Eg)\cdot \nabla \varphi \,\mathrm {d} x=\int _{\Omega }f\varphi \,\mathrm {d} x-\int _{\Omega }\nabla Eg\cdot \nabla \varphi \,\mathrm {d} x for all {\textstyle \varphi \in H_{0}^{1}(\Omega ).

Thus the problem with inhomogeneous boundary values for {\textstyle u could be reduced to a problem with homogeneous boundary values for {\textstyle u_{0}, a technique which can be applied to any linear differential equation. By the Riesz representation theorem there exists a unique solution {\textstyle u_{0} to this problem. By uniqueness of the decomposition {\textstyle u=u_{0}+Eg, this is equivalent to the existence of a unique weak solution {\textstyle u to the inhomogeneous boundary value problem.

Continuous dependence on the data

It remains to investigate the dependence of {\textstyle u on {\textstyle f and {\textstyle g. Let {\textstyle c_{1},c_{2},\ldots >0 denote constants independent of {\textstyle f and {\textstyle g. By continuous dependence of {\textstyle u_{0} on the right-hand side of its integral equation, there holds

\|u_{0}\|_{H_{0}^{1}(\Omega )}\leq c_{1}\left(\|f\|_{H^{-1}(\Omega )}+\|Eg\|_{H^{1}(\Omega )}\right)

and thus, using that {\textstyle \|u_{0}\|_{H^{1}(\Omega )}\leq c_{2}\|u_{0}\|_{H_{0}^{1}(\Omega )} and {\textstyle \|Eg\|_{H^{1}(\Omega )}\leq c_{3}\|g\|_{H^{1/2}(\Omega )} by continuity of the trace extension operator, it follows that

{\begin{aligned}\|u\|_{H^{1}(\Omega )}&\leq \|u_{0}\|_{H^{1}(\Omega )}+\|Eg\|_{H^{1}(\Omega )}\leq c_{1}c_{2}\|f\|_{H^{-1}(\Omega )}+(c_{3}+c_{1}c_{2})\|Eg\|_{H^{1}(\Omega )}\\&\leq c_{4}\left(\|f\|_{H^{-1}(\Omega )}+\|g\|_{H^{1/2}(\partial \Omega )}\right)\end{aligned}}

and the solution map

H^{-1}(\Omega )\times H^{1/2}(\partial \Omega )\ni (f,g)\mapsto u\in H^{1}(\Omega )

is therefore continuous.

Watch videos about Trace operatorExplainers and documentaries on YouTube (opens in a new tab)

Sources and credits

This article is adapted from the Wikipedia article Trace operator, written by its contributors and licensed under CC BY-SA 4.0. Fathomly has changed the layout, removed citation markers, navigation and maintenance notices, and adjusted punctuation. This adapted version is shared under the same license. For references, see the original article.

Images, from Wikimedia Commons:

Fathomly is not affiliated with or endorsed by the Wikimedia Foundation. Spotted a problem? Tell us.