Sobolev space

In mathematics, a Sobolev space is a vector space of functions equipped with a norm that is a combination of "L^p" norms of the function itself as well as its derivatives up to a given order. The derivatives are understood in a suitable weak sense to make the space complete, thus a Banach space. Intuitively, a Sobolev space is a Banach space or Hilbert space of functions with sufficiently many derivatives for some application domain, such as partial differential equations, and equipped with a norm that measures both the size and smoothness of a function.

Sobolev spaces are named after the Russian mathematician Sergei L. Sobolev. Their importance lies in the fact that solutions of partial differential equations are naturally in Sobolev spaces rather than in the classical spaces of continuous functions and with the derivatives understood in the classical sense.

Introduction

There are many criteria for smoothness of mathematical functions. The most basic criterion may be that of continuity. A considerably stronger notion of smoothness is that of differentiability (because functions that are differentiable are also continuous) and a yet stronger notion of smoothness is that the derivative also be continuous (these functions are said to be of class "C"¹ — see smooth function). Differentiable functions are important in many areas, and in particular for differential equations. In the twentieth century, however, it was observed that the space "C"¹ (or "C"², etc.) was not exactly the right space to study solutions of differential equations.

The Sobolev spaces are the modern replacement for these spaces in which to look for solutions of partial differential equations.

Sobolev spaces on the unit circle

We start by introducing Sobolev spaces in the simplest settings, the one-dimensional case on the unit circle. In this case the Sobolev space "W"^k,p is defined to be the subset of "L"^p such that function "f" and its weak derivatives up to some order "k" have a finite "L"^p norm, for given "p" ≥ 1. Some care must be taken to define derivatives in the proper sense. In the one-dimensional problem it is enough to assume that the ("k"-1)-th derivative of function "f", "f"^(k-1), is differentiable almost everywhere and is equal almost everywhere to the Lebesgue integral of its derivative (this gets rid of examples such as Cantor's function which are irrelevant to what the definition is trying to accomplish).

With this definition, the Sobolev spaces admit a natural norm,

:$|f|\_\{k,p\}=Big(sum\_\{i=0\}^k\; |f^\{(i)\}|\_p^pBig)^\{1/p\}\; =\; Big(sum\_\{i=0\}^k\; int\; |f^\{(i)\}(t)|^p,dt\; Big)^\{1/p\}.$

Space "W"^k,p equipped with the norm $|cdot|\_\{k,p\}$ is a Banach space. It turns out that it is enough to take only the first and last in the sequence, i.e., the norm defined by

:$|f^\{(k)\}|\_p\; +\; |f|\_p$

is equivalent to the norm above.

=The case "p" = 2=

Sobolev spaces with "p" = 2 are especially important because of their connection with Fourier series and because they form a Hilbert space. A special notation has arisen to cover this case::$,H^k\; =\; W^\{k,2\}.$

The space $H^k$ can be defined naturally in terms of Fourier series whose coefficients decay sufficiently rapidly, namely,

:

where $widehat\{f\}$ is the Fourier series of $f$. As above, one can use the equivalent norm

:$|f|^2=sum\_\{n=-infty\}^infty\; (1\; +\; |n|^\{k\})^2\; |widehat\{f\}(n)|^2.$

Both representations follow easily from Parseval's theorem and the fact that differentiation is equivalent to multiplying the Fourier coefficient by "in".

Furthermore, the space "H"^"k" admits an inner product, like the space "H"⁰ = "L"². In fact, the "H"^"k" inner product is defined in terms of the "L"² inner product::$langle\; u,v\; angle\_\{H^k\}=sum\_\{i=0\}^klangle\; D^i\; u,D^i\; v\; angle\_\{L\_2\}.$The space "H"^"k" becomes a Hilbert space with this inner product.

Other examples

Some other Sobolev spaces permit a simpler description. For example, $W^\{1,1\}(0,1)$ is the space of absolutely continuous functions on $(0,1)$, while "W"^1,∞("I") is the space of Lipschitz functions on $I$, for every interval $I$.All spaces "W"^k,∞ are (normed) algebras, i.e. the product of two elements is once again a function of this Sobolev space, which is not the case for "p" < ∞. (E.g., functions behaving like |"x"|^−1/3 at the origin are in "L"², but the product of two such functions is not in "L"²).

obolev spaces with non-integer "k"

To prevent confusion, when talking about "k" which is not integer we will usually denote it by "s", i.e. $W^\{s,p\}$ or $H^s.$

= The case "p" = 2 =

The case "p" = 2 is the easiest since the Fourier description is straightforward to generalize. We define the norm

:$|f|^2\_\{s,2\}=sum\; (1+n^2)^s|widehat\{f\}(n)|^2$

and the Sobolev space $H^s$ as the space of all functions with finite norm.

Fractional order differentiation

A similar approach can be used if "p" is different from 2. In this case Parseval's theorem no longer holds, but differentiation still corresponds to multiplication in the Fourier domain and can be generalized to non-integer orders. Therefore we define an operator of "fractional order differentiation" of order "s" by

:$F^s(f)=sum\_\{n=-infty\}^infty\; (in)^swidehat\{f\}(n)e^\{int\}$

or in other words, taking Fourier transform, multiplying by $(in)^s$ and then taking inverse Fourier transform (operators defined by Fourier-multiplication-inverse Fourier are called multipliers and are a topic of research in their own right). This allows to define the Sobolev norm of $s,p$ by

:$|f|\_\{s,p\}=|f|\_p+|F^s(f)|\_p$

and, as usual, the Sobolev space is the space of functions with finite Sobolev norm.

Complex interpolation

Another way of obtaining the "fractional Sobolev spaces" is given by complex interpolation. Complex interpolation is a general technique: for any 0 ≤ t ≤ 1 and "X" and "Y" Banach spaces that are continuously included in some larger Banach space we may create "intermediate space" denoted ["X","Y"] _"t". (below we discuss a different method, the so-called real interpolation method, which is essential in the Sobolev theory for the characterization of traces).

Such spaces "X" and "Y" are called interpolation pairs.

We mention a couple of useful theorems about complex interpolation:

"Theorem (reinterpolation): " [ ["X,Y"] _"a" , ["X,Y"] _"b" ] _"c" = ["X,Y"] _{"cb"+(1-"c")"a"}.

"Theorem (interpolation of operators): if "{"X,Y"}" and "{"A,B"}" are interpolation pairs, and if T is a linear map defined on X"+"Y into A"+"B so that T is continuous from X to A and from Y to B then T is continuous from " ["X,Y"] "_t to " ["A,B"] "_t. and we have the interpolation inequality:"

$|T|\_\{\; [X,Y]\; \_t\; o\; [A,B]\; \_t\}leq\; C|T|\_\{X\; o\; A\}^\{1-t\}|T|\_\{Y\; o\; B\}^t.$

See also: Riesz-Thorin theorem.

Returning to Sobolev spaces, we want to get $W^\{s,p\}$ for non-integer "s" by interpolating between $W^\{k,p\}$-s. The first thing is of course to see that this gives consistent results, and indeed we have

"Theorem: $left\; [W^\{0,p\},W^\{m,p\}\; ight]\; \_t=W^\{n,p\}$ if n is an integer such that n=tm."

Hence, complex interpolation is a consistent way to get a continuum of spaces $W^\{s,p\}$ between the $W^\{k,p\}$. Further, it gives the same spaces as fractional order differentiation does (but see extension operators below for a twist).

Multiple dimensions

We now turn to the case of Sobolev spaces in R^"n" and subsets of R^"n". The change from the circle to the line only entails technical changes in the Fourier formulas — basically a change of Fourier series to Fourier transform and sums to integrals. The transition to multiple dimensions brings more difficulties, starting from the very definition. The requirement that "f"^("k"−1) is the integral of "f"^("k") does not generalize, and the simplest solution is to consider derivatives in the sense of distribution theory.

A formal definition now follows. Let "D" be an open set in Rⁿ, let "k" be a natural number and let 1 ≤ "p" ≤ +∞. The Sobolev space "W"^"k","p"("D") is defined to be the set of all functions "f" defined on "D" such that for every multi-index "α" with |"α"| ≤ "k", the mixed partial derivative

:$f^\{(alpha)\}\; =\; frac\{partial^\; f\}\{partial\; x\_\{1\}^\{alpha\_\{1\; dots\; partial\; x\_\{n\}^\{alpha\_\{n\}$

is both locally integrable and in "L"^"p"("D"), i.e.

:

There are several choices of norm for "W"^"k","p"("D"). The following two are common, and are equivalent in the sense of equivalence of norms:

:

and

:

With respect to either of these norms, "W"^"k","p"("D") is a Banach space. For finite "p", "W"^"k","p"("D") is also a separable space. As noted above, it is conventional to denote "W"^"k",2("D") by "H"^"k"("D").

The fractional order Sobolev spaces "H"^"s"(R^"n"), "s" ≥ 0, can be defined using the Fourier transform (using the fact that the Fourier Transform is a unitary transformation) as before:

:

However, if "D" is not a periodic domain like R^"n" or the torus T^"n", this definition is insufficient, since the Fourier transform of a function defined on an aperiodic domain is difficult to define. Fortunately, there is an intrinsic characterization of fractional order Sobolev spaces using what is essentially the "L"² analogue of Hölder continuity: an equivalent inner product for "H"^"s"("D") is given by

:

where "s" = "k" + "t", "k" an integer and 0 < "t" < 1. Note that the dimension of the domain, "n", appears in the above formula for the inner product.

Examples

In higher dimensions, it is no longer true that, for example, "W"^1,1 contains only continuous functions. For example, 1/|"x"| belongs to "W"^1,1(B³) where B³ is the unit ball in three dimensions. For "k" > "n"/"p" the space "W"^"k","p"("D") will contain only continuous functions, but for which "k" this is already true depends both on "p" and on the dimension. For example, as can be easily checked using spherical polar coordinates, the function "f" : B^"n" → R ∪ {+∞} defined on the "n"-dimensional ball and given by

:$f(x)\; =\; frac1\{|\; x\; |^\{alpha$

lies in "W"^"k","p"(B^"n") if and only if

:

Intuitively, the blow-up of "f" at 0 "counts for less" when "n" is large since the unit ball is "smaller" in higher dimensions.

ACL characterization of Sobolev functions

Let Ω be an open set in Rⁿ and 1≤p≤∞. If a function is in "W"^1,p(Ω), then the restriction to every line parallel to the coordinate directions in Rⁿ is absolutely continuous. Conversely, if the restriction of "f" to every line parallel to the coordinate directions is absolutely continuous, then the pointwise gradient ∇"f" exists almost everywhere, and "f" is in "W"^1,p(Ω) provided "f" and |∇"f"| are both in "L"^p(Ω). In particular, in this case the weak partial derivatives of "f" and pointwise partial derivatives of "f" agree almost everywhere.

obolev embedding

Write $W^\{k,p\}$ for the Sobolev space of some compact Riemannian manifold of dimension "n".Here "k" can be any real number, and 1 &le; "p" &le; &infin;. (For "p" = &infin; the Sobolev space $W^\{k,infty\}$ is defined to be the Hölder space "C"^"n",α where "k" = "n" + α and 0 < α &le; 1.) The Sobolev embedding theorem states that if "k" &ge; ℓ and "k" − "n"/"p" &ge; ℓ − "n"/"q" then

:$W^\{k,p\}subseteq\; W^\{ell,q\}$

and the embedding is continuous. Moreover if "k" > ℓ and "k" − "n"/"p" > ℓ −"n"/"q"then the embedding is completely continuous (this is sometimes called Kondrakov's theorem). Functions in $W^\{ell,infty\}$ have all derivatives of order less than ℓ continuous, so in particular this gives conditions on Sobolev spaces for various derivatives to be continuous. Informally these embeddings say that to convert an "L"^"p" estimate to a boundedness estimate costs 1/"p" derivatives per dimension.

There are similar variations of the embedding theorem for non-compact manifolds such as R^"n" harv|Stein|1970:

Traces

:"Main article Trace operator."

Let "s" > &frac12;. If "X" is an open set such that its boundary "G" is "sufficiently smooth", then we may define the "trace" (that is, "restriction") map "P" by

:$Pu=u|\_G,$

i.e. "u" restricted to "G". A simple smoothness condition is uniformly $C^m$, "m" &ge; "s". (There is no connection here to trace of a matrix.)

This trace map "P" as defined has domain $H^s(X)$, and its image is precisely $H^\{s-1/2\}(G)$. To be completely formal, "P" is first defined for infinitely differentiable functions and is extended by continuity to $H^s(X)$. Note that we 'lose half a derivative' in taking this trace.

Identifying the image of the trace map for $W^\{s,p\}$ is considerably more difficult and demands the tool of real interpolation. The resulting spaces are the Besov spaces. It turns out that in the case of the $W^\{s,p\}$ spaces, we don't lose half a derivative; rather, we lose 1/"p" of a derivative.

Extension operators

If "X" is an open domain whose boundary is not too poorly behaved (e.g., if its boundary is a manifold, or satisfies the more permissive but more obscure "cone condition") then there is an operator "A" mapping functions of "X" to functions of R^"n" such that:

# "Au"("x") = "u"("x") for almost every "x" in "X" and
# "A" is continuous from $W^\{k,p\}(X)$ to $W^\{k,p\}(\{mathbb\; R\}^n)$, for any 1 &le; "p" &le; &infin; and integer "k".

We will call such an operator "A" an extension operator for "X".

Extension operators are the most natural way to define $H^s(X)$ for non-integer "s" (we cannot work directly on "X" since taking Fourier transform is a global operation). We define $H^s(X)$ by saying that "u" is in $H^s(X)$ if and only if "Au" is in $H^s(mathbb\; R^n)$. Equivalently, complex interpolation yields the same $H^s(X)$ spaces so long as "X" has an extension operator. If "X" does not have an extension operator, complex interpolation is the only way to obtain the $H^s(X)$ spaces.

As a result, the interpolation inequality still holds.

Extension by zero

We define $H^s\_0(X)$ to be the closure in $H^s(X)$ of the space $C^infty\_c(X)$ of infinitely differentiable compactly supported functions. Given the definition of a trace, above, we may state the following

"Theorem: Let X be uniformly C^m regular, m &ge; s and let P be the linear map sending u in $H^s(X)$ to "

:$left.left(u,frac\{du\}\{dn\},\; dots,\; frac\{d^k\; u\}\{dn^k\}\; ight)\; ight|\_G$

"where d/dn is the derivative normal to G, and k is the largest integer less than s. Then $H^s\_0$ is precisely the kernel of P."

If $uin\; H^s\_0(X)$ we may define its extension by zero $ilde\; u\; in\; L^2(\{mathbb\; R\}^n)$ in the natural way, namely

:$ilde\; u(x)=u(x)\; ;\; extrm\{\; if\; \}\; ;\; x\; in\; X,\; 0\; ;\; extrm\{\; otherwise.\}$

"Theorem: Let s > &frac12;. The map taking u to $ilde\; u$ is continuous into $H^s(\{mathbb\; R\}^n)$ if and only if s is not of the form n + &frac12; for n an integer."

References

* R.A. Adams, J.J.F. Fournier, 2003. "Sobolev Spaces". Academic Press.
* L.C. Evans, 1998. "Partial Differential Equations". American Mathematical Society.
*springer|id=i/i050230|title=Imbedding theorems|first=S.M.|last= Nikol'skii
*springer|id=S/s085980|title=Sobolev space|first=S.M.|last= Nikol'skii
*S.L. Sobolev, "On a theorem of functional analysis" Transl. Amer. Math. Soc. (2) , 34 (1963) pp. 39–68 Mat. Sb. , 4 (1938) pp. 471–497
*S.L. Sobolev, "Some applications of functional analysis in mathematical physics" , Amer. Math. Soc. (1963)
*citation|last=Stein|first=E|title= Singular Integrals and Differentiability Properties of Functions, |publisher=Princeton Univ. Press|year=1970| ISBN= 0-691-08079-8