In
mathematics
Mathematics is an area of knowledge that includes the topics of numbers, formulas and related structures, shapes and the spaces in which they are contained, and quantities and their changes. These topics are represented in modern mathematics ...
, Sobolev spaces for planar domains are one of the principal techniques used in the theory of
partial differential equation
In mathematics, a partial differential equation (PDE) is an equation which imposes relations between the various partial derivatives of a Multivariable calculus, multivariable function.
The function is often thought of as an "unknown" to be sol ...
s for solving the
Dirichlet
Johann Peter Gustav Lejeune Dirichlet (; 13 February 1805 – 5 May 1859) was a German mathematician who made deep contributions to number theory (including creating the field of analytic number theory), and to the theory of Fourier series and ...
and
Neumann boundary value problems for the
Laplacian
In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols \nabla\cdot\nabla, \nabla^2 (where \nabla is the ...
in a bounded domain in the plane with smooth boundary. The methods use the theory of
bounded operator
In functional analysis and operator theory, a bounded linear operator is a linear transformation L : X \to Y between topological vector spaces (TVSs) X and Y that maps bounded subsets of X to bounded subsets of Y.
If X and Y are normed vector ...
s on
Hilbert space
In mathematics, Hilbert spaces (named after David Hilbert) allow generalizing the methods of linear algebra and calculus from (finite-dimensional) Euclidean vector spaces to spaces that may be infinite-dimensional. Hilbert spaces arise natural ...
. They can be used to deduce regularity properties of solutions and to solve the corresponding eigenvalue problems.
Sobolev spaces with boundary conditions
Let be a bounded domain with smooth boundary. Since is contained in a large square in , it can be regarded as a domain in by identifying opposite sides of the square. The theory of Sobolev spaces on can be found in , an account which is followed in several later textbooks such as and .
For an integer, the (restricted) Sobolev space is defined as the closure of in the standard
Sobolev space
In mathematics, a Sobolev space is a vector space of functions equipped with a norm that is a combination of ''Lp''-norms of the function together with its derivatives up to a given order. The derivatives are understood in a suitable weak sense t ...
.
* .
*Vanishing properties on boundary: For the elements of are referred to as " functions on which vanish with their first derivatives on ." In fact if agrees with a function in , then is in . Let be such that in the Sobolev norm, and set . Thus in . Hence for and ,
::
:By
Green's theorem
In vector calculus, Green's theorem relates a line integral around a simple closed curve to a double integral over the plane region bounded by . It is the two-dimensional special case of Stokes' theorem.
Theorem
Let be a positively orient ...
this implies
::
:where
::
:with the unit normal to the boundary. Since such form a dense subspace of , it follows that on .
*Support properties: Let be the complement of and define restricted Sobolev spaces analogously for . Both sets of spaces have a natural pairing with . The Sobolev space for is the annihilator in the Sobolev space for of and that for is the annihilator of . In fact this is proved by locally applying a small translation to move the domain inside itself and then smoothing by a smooth convolution operator.
:Suppose in annihilates . By compactness, there are finitely many open sets covering such that the closure of is disjoint from and each is an open disc about a boundary point such that in small translations in the direction of the normal vector carry into . Add an open with closure in to produce a cover of and let be a
partition of unity
In mathematics, a partition of unity of a topological space is a set of continuous functions from to the unit interval ,1such that for every point x\in X:
* there is a neighbourhood of where all but a finite number of the functions of are 0, ...
subordinate to this cover. If translation by is denoted by , then the functions
::
:tend to as decreases to and still lie in the annihilator, indeed they are in the annihilator for a larger domain than , the complement of which lies in . Convolving by smooth functions of small support produces smooth approximations in the annihilator of a slightly smaller domain still with complement in . These are necessarily smooth functions of compact support in .
*Further vanishing properties on the boundary: The characterization in terms of annihilators shows that lies in if (and only if) it and its derivatives of order less than vanish on . In fact can be extended to by setting it to be on . This extension defines an element in using the formula for the norm
::
:Moreover satisfies for ''g'' in .
*Duality: For , define to be the orthogonal complement of in . Let be the orthogonal projection onto , so that is the orthogonal projection onto . When , this just gives . If and , then
::
:This implies that under the pairing between and , and are each other's duals.
*Approximation by smooth functions: The image of is dense in for . This is obvious for since the sum + is dense in . Density for follows because the image of is dense in and annihilates .
*Canonical isometries: The operator gives an isometry of into and of onto . In fact the first statement follows because it is true on . That is an isometry on follows using the density of in : for we have:
::
:Since the adjoint map between the duals can by identified with this map, it follows that is a unitary map.
Application to Dirichlet problem
Invertibility of
The operator defines an isomorphism between and . In fact it is a Fredholm operator of index . The kernel of in consists of constant functions and none of these except zero vanish on the boundary of . Hence the kernel of is and is invertible.
In particular the equation has a unique solution in for in .
Eigenvalue problem
Let be the operator on defined by
:
where is the inclusion of in and of in , both compact operators by Rellich's theorem. The operator is compact and self-adjoint with for all . By the
spectral theorem
In mathematics, particularly linear algebra and functional analysis, a spectral theorem is a result about when a linear operator or matrix (mathematics), matrix can be Diagonalizable matrix, diagonalized (that is, represented as a diagonal matrix i ...
, there is a complete orthonormal set of eigenfunctions in with
:
Since , lies in . Setting , the are eigenfunctions of the Laplacian:
:
Sobolev spaces without boundary condition
To determine the regularity properties of the eigenfunctions and solutions of
:
enlargements of the Sobolev spaces have to be considered. Let be the space of smooth functions on which with their derivatives extend continuously to . By
Borel's lemma
In mathematics, Borel's lemma, named after Émile Borel, is an important result used in the theory of asymptotic expansions and partial differential equations.
Statement
Suppose ''U'' is an open set in the Euclidean space R''n'', and suppose that ...
, these are precisely the restrictions of smooth functions on . The Sobolev space is defined to the Hilbert space completion of this space for the norm
:
This norm agrees with the Sobolev norm on so that can be regarded as a closed subspace of . Unlike , is not naturally a subspace of , but the map restricting smooth functions from to is continuous for the Sobolev norm so extends by continuity to a map .
*Invariance under diffeomorphism: Any diffeomorphism between the closures of two smooth domains induces an isomorphism between the Sobolev space. This is a simple consequence of the chain rule for derivatives.
*Extension theorem: The restriction of to the orthogonal complement of its kernel defines an isomorphism onto . The extension map is defined to be the inverse of this map: it is an isomorphism (not necessarily norm preserving) of onto the orthogonal complement of such that . On , it agrees with the natural inclusion map. Bounded extension maps of this kind from to were constructed first constructed by Hestenes and Lions. For smooth curves the
Seeley extension theorem provides an extension which is continuous in all the Sobolev norms. A version of the extension which applies in the case where the boundary is just a Lipschitz curve was constructed by
Calderón
Calderón () is a Spanish and Sefardi occupational surname. It is derived from the Vulgar Latin "''caldaria''" ("cauldron") and refers to the occupation of tinker.
Calderón, or Calderon, may refer to:
* Alberto Calderón, Argentine mathematician ...
using
singular integral operators In mathematics, singular integrals are central to harmonic analysis and are intimately connected with the study of partial differential equations. Broadly speaking a singular integral is an integral operator
: T(f)(x) = \int K(x,y)f(y) \, dy,
w ...
and generalized by .
:It is sufficient to construct an extension for a neighbourhood of a closed annulus, since a collar around the boundary is diffeomorphic to an annulus with a closed interval in . Taking a smooth bump function with , equal to 1 near the boundary and 0 outside the collar, will provide an extension on . On the annulus, the problem reduces to finding an extension for in . Using a partition of unity the task of extending reduces to a neighbourhood of the end points of . Assuming 0 is the left end point, an extension is given locally by
::
:Matching the first derivatives of order ''k'' or less at 0, gives
::
:This matrix equation is solvable because the determinant is non-zero by
Vandermonde's formula. It is straightforward to check that the formula for , when appropriately modified with bump functions, leads to an extension which is continuous in the above Sobolev norm.
*Restriction theorem: The restriction map is surjective with . This is an immediate consequence of the extension theorem and the support properties for Sobolev spaces with boundary condition.
*Duality: is naturally the dual of H
−''k''0(Ω). Again this is an immediate consequence of the restriction theorem. Thus the Sobolev spaces form a chain:
::
:The differentiation operators carry each Sobolev space into the larger one with index 1 less.
*Sobolev embedding theorem: is contained in . This is an immediate consequence of the extension theorem and the Sobolev embedding theorem for .
*Characterization: consists of in such that all the derivatives ∂
α''f'' lie in for , α, ≤ ''k''. Here the derivatives are taken within the chain of Sobolev spaces above. Since is weakly dense in , this condition is equivalent to the existence of functions ''f''
α such that
::
:To prove the characterization, note that if is in , then lies in H
''k''−, α, (Ω) and hence in . Conversely the result is well known for the Sobolev spaces : the assumption implies that the is in and the corresponding condition on the Fourier coefficients of shows that lies in . Similarly the result can be proved directly for an annulus . In fact by the argument on the restriction of to any smaller annulus
ˆ’δ',δ'× T lies in : equivalently the restriction of the function lies in for . On the other hand in as , so that must lie in . The case for a general domain reduces to these two cases since can be written as with ψ a bump function supported in such that is supported in a collar of the boundary.
*Regularity theorem: If in has both derivatives and in then lies in . This is an immediate consequence of the characterization of above. In fact if this is true even when satisfied at the level of distributions: if there are functions ''g'', ''h'' in such that (''g'',φ) = (''f'', φ
''x'') and (''h'',φ) = (''f'',φ
''y'') for φ in , then is in .
*Rotations on an annulus: For an annulus , the extension map to is by construction equivariant with respect to rotations in the second variable,
::
:On it is known that if is in , then the ''difference quotient'' in ; if the difference quotients are bounded in ''H''
''k'' then ∂
''y''''f'' lies in . Both assertions are consequences of the formula:
::
:These results on imply analogous results on the annulus using the extension.
Regularity for Dirichlet problem
Regularity for dual Dirichlet problem
If with in and in with , then lies in .
Take a decomposition with supported in and supported in a collar of the boundary. Standard Sobolev theory for can be applied to : elliptic regularity implies that it lies in and hence . lies in of a collar, diffeomorphic to an annulus, so it suffices to prove the result with a collar and replaced by
:
The proof proceeds by induction on , proving simultaneously the inequality
:
for some constant depending only on . It is straightforward to establish this inequality for , where by density can be taken to be smooth of compact support in :
:
The collar is diffeomorphic to an annulus. The rotational flow on the annulus induces a flow on the collar with corresponding vector field . Thus corresponds to the vector field . The radial vector field on the annulus is a commuting vector field which on the collar gives a vector field proportional to the normal vector field. The vector fields and commute.
The difference quotients can be formed for the flow . The commutators are second order differential operators from to . Their operators norms are uniformly bounded for near ; for the computation can be carried out on the annulus where the commutator just replaces the coefficients of by their difference quotients composed with . On the other hand, lies in , so the inequalities for apply equally well for :
:
The uniform boundedness of the difference quotients implies that lies in with
:
It follows that lies in where is the vector field
:
Moreover, satisfies a similar inequality to .
:
Let be the orthogonal vector field
:
It can also be written as for some smooth nowhere vanishing function on a neighbourhood of the collar.
It suffices to show that lies in . For then
:
so that and lie in and must lie in .
To check the result on , it is enough to show that and lie in . Note that
:
are vector fields. But then
:
with all terms on the right hand side in . Moreover, the inequalities for show that
:
Hence
:
Smoothness of eigenfunctions
It follows by induction from the regularity theorem for the dual Dirichlet problem that the eigenfunctions of in lie in . Moreover, any solution of with in and in must have in . In both cases by the vanishing properties, the eigenfunctions and vanish on the boundary of .
Solving the Dirichlet problem
The dual Dirichlet problem can be used to solve the Dirichlet problem:
:
By Borel's lemma is the restriction of a function in . Let be the smooth solution of with on . Then solves the Dirichlet problem. By the
maximal principle, the solution is unique.
Application to smooth Riemann mapping theorem
The solution to the Dirichlet problem can be used to prove a strong form of the
Riemann mapping theorem
In complex analysis, the Riemann mapping theorem states that if ''U'' is a non-empty simply connected space, simply connected open set, open subset of the complex plane, complex number plane C which is not all of C, then there exists a biholomorphy ...
for simply connected domains with smooth boundary. The method also applies to a region diffeomorphic to an annulus. For multiply connected regions with smooth boundary have given a method for mapping the region onto a disc with circular holes. Their method involves solving the Dirichlet problem with a non-linear boundary condition. They construct a function such that:
* is harmonic in the interior of ;
* On we have: , where is the curvature of the boundary curve, is the derivative in the direction normal to and is constant on each boundary component.
gives a proof of the Riemann mapping theorem for a simply connected domain with smooth boundary. Translating if necessary, it can be assumed that . The solution of the Dirichlet problem shows that there is a unique smooth function on which is harmonic in and equals on . Define the
Green's function
In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.
This means that if \operatorname is the linear differenti ...
by . It vanishes on and is harmonic on away from . The
harmonic conjugate
In mathematics, a real-valued function u(x,y) defined on a connected open set \Omega \subset \R^2 is said to have a conjugate (function) v(x,y) if and only if they are respectively the real and imaginary parts of a holomorphic function f(z) of ...
of is the unique real function on such that is holomorphic. As such it must satisfy the
Cauchy–Riemann equations
In the field of complex analysis in mathematics, the Cauchy–Riemann equations, named after Augustin Cauchy and Bernhard Riemann, consist of a system of two partial differential equations which, together with certain continuity and differ ...
:
:
The solution is given by
:
where the integral is taken over any path in . It is easily verified that and exist and are given by the corresponding derivatives of . Thus is a smooth function on , vanishing at . By the Cauchy-Riemann is smooth on , holomorphic on and . The function is only defined up to multiples of , but the function
:
is a holomorphic on and smooth on . By construction, and for . Since has
winding number
In mathematics, the winding number or winding index of a closed curve in the plane around a given point is an integer representing the total number of times that curve travels counterclockwise around the point, i.e., the curve's number of turn ...
, so too does . On the other hand, only for where there is a simple zero. So by the
argument principle
In complex analysis, the argument principle (or Cauchy's argument principle) relates the difference between the number of zeros and poles of a meromorphic function to a contour integral of the function's logarithmic derivative.
Specifically, i ...
assumes every value in the unit disc, , exactly once and does not vanish inside . To check that the derivative on the boundary curve is non-zero amounts to computing the derivative of , i.e. the derivative of should not vanish on the boundary curve. By the Cauchy-Riemann equations these tangential derivative are up to a sign the
directional derivative
In mathematics, the directional derivative of a multivariable differentiable (scalar) function along a given vector v at a given point x intuitively represents the instantaneous rate of change of the function, moving through x with a velocity s ...
in the direction of the normal to the boundary. But vanishes on the boundary and is strictly negative in since . The
Hopf lemma In mathematics, the Hopf lemma, named after Eberhard Hopf, states that if a continuous real-valued function in a domain in Euclidean space with sufficiently smooth boundary is harmonic in the interior and the value of the function at a point on the ...
implies that the directional derivative of in the direction of the outward normal is strictly positive. So on the boundary curve, has nowhere vanishing derivative. Since the boundary curve has winding number one, defines a diffeomorphism of the boundary curve onto the unit circle. Accordingly, is a smooth diffeomorphism, which restricts to a holomorphic map and a smooth diffeomorphism between the boundaries.
Similar arguments can be applied to prove the Riemann mapping theorem for a doubly connected domain bounded by simple smooth curves (the inner curve) and (the outer curve). By translating we can assume 1 lies on the outer boundary. Let be the smooth solution of the Dirichlet problem with on the outer curve and on the inner curve. By the
maximum principle
In the mathematical fields of partial differential equations and geometric analysis, the maximum principle is any of a collection of results and techniques of fundamental importance in the study of elliptic and parabolic differential equations.
...
for in and so by the
Hopf lemma In mathematics, the Hopf lemma, named after Eberhard Hopf, states that if a continuous real-valued function in a domain in Euclidean space with sufficiently smooth boundary is harmonic in the interior and the value of the function at a point on the ...
the normal derivatives of are negative on the outer curve and positive on the inner curve. The integral of over the boundary is zero by Stoke's theorem so the contributions from the boundary curves cancel. On the other hand, on each boundary curve the contribution is the integral of the normal derivative along the boundary. So there is a constant such that satisfies
:
on each boundary curve. The harmonic conjugate of can again be defined by
:
and is well-defined up to multiples of . The function
:
is smooth on and holomorphic in . On the outer curve and on the inner curve . The tangential derivatives on the outer curves are nowhere vanishing by the Cauchy-Riemann equations, since the normal derivatives are nowhere vanishing. The normalization of the integrals implies that restricts to a diffeomorphism between the boundary curves and the two concentric circles. Since the images of outer and inner curve have winding number and about any point in the annulus, an application of the argument principle implies that assumes every value within the annulus exactly once; since that includes multiplicities, the complex derivative of is nowhere vanishing in . This is a smooth diffeomorphism of onto the closed annulus , restricting to a holomorphic map in the interior and a smooth diffeomorphism on both boundary curves.
Trace map
The restriction map extends to a continuous map for . In fact
:
so the
Cauchy–Schwarz inequality
The Cauchy–Schwarz inequality (also called Cauchy–Bunyakovsky–Schwarz inequality) is considered one of the most important and widely used inequalities in mathematics.
The inequality for sums was published by . The corresponding inequality fo ...
yields
:
where, by the
integral test
In mathematics, an integral assigns numbers to functions in a way that describes displacement, area, volume, and other concepts that arise by combining infinitesimal data. The process of finding integrals is called integration. Along with d ...
,
:
The map is onto since a continuous extension map can be constructed from to . In fact set
:
where
:
Thus . If ''g'' is smooth, then by construction ''Eg'' restricts to ''g'' on 1 × T. Moreover, ''E'' is a bounded linear map since
:
It follows that there is a trace map Ï„ of H
''k''(Ω) onto H
''k'' − ½(∂Ω). Indeed, take a tubular neighbourhood of the boundary and a smooth function ψ supported in the collar and equal to 1 near the boundary. Multiplication by ψ carries functions into H
''k'' of the collar, which can be identified with H
''k'' of an annulus for which there is a trace map. The invariance under diffeomorphisms (or coordinate change) of the half-integer Sobolev spaces on the circle follows from the fact that an equivalent norm on ''H''
''k'' + ½(T) is given by
:
It is also a consequence of the properties of Ï„ and ''E'' (the "trace theorem"). In fact any diffeomorphism ''f'' of T induces a diffeomorphism ''F'' of T
2 by acting only on the second factor. Invariance of H
''k''(T
2) under the induced map ''F''* therefore implies invariance of H
''k'' − ½(T) under ''f''*, since ''f''* = τ ∘ ''F''* ∘ ''E''.
Further consequences of the trace theorem are the two exact sequences
:
and
:
where the last map takes ''f'' in H
2(Ω) to ''f'',
∂Ω and ∂
''n''''f'',
∂Ω. There are generalizations of these sequences to H
''k''(Ω) involving higher powers of the normal derivative in the trace map:
:
The trace map to takes ''f'' to
Abstract formulation of boundary value problems
The Sobolev space approach to the Neumann problem cannot be phrased quite as directly as that for the Dirichlet problem. The main reason is that for a function in , the normal derivative cannot be a priori defined at the level of Sobolev spaces. Instead an alternative formulation of boundary value problems for the Laplacian on a bounded region in the plane is used. It employs
Dirichlet forms, sesqulinear bilinear forms on , or an intermediate closed subspace. Integration over the boundary is not involved in defining the Dirichlet form. Instead, if the Dirichlet form satisfies a certain positivity condition, termed coerciveness, solution can be shown to exist in a weak sense, so-called "weak solutions". A general regularity theorem than implies that the solutions of the boundary value problem must lie in , so that they are strong solutions and satisfy boundary conditions involving the restriction of a function and its normal derivative to the boundary. The Dirichlet problem can equally well be phrased in these terms, but because the trace map is already defined on , Dirichlet forms do not need to be mentioned explicitly and the operator formulation is more direct. A unified discussion is given in and briefly summarised below. It is explained how the Dirichlet problem, as discussed above, fits into this framework. Then a detailed treatment of the Neumann problem from this point of view is given following .
The Hilbert space formulation of boundary value problems for the Laplacian on a bounded region in the plane proceeds from the following data:
* A closed subspace .
* A Dirichlet form for given by a bounded Hermitian bilinear form defined for such that for .
* is coercive, i.e. there is a positive constant and a non-negative constant such that .
A weak solution of the boundary value problem given initial data in is a function ''u'' satisfying
:
for all ''g''.
For both the Dirichlet and Neumann problem
:
For the Dirichlet problem . In this case
:
By the trace theorem the solution satisfies in .
For the Neumann problem is taken to be .
Application to Neumann problem
The classical Neumann problem on consists in solving the boundary value problem
:
Green's theorem
In vector calculus, Green's theorem relates a line integral around a simple closed curve to a double integral over the plane region bounded by . It is the two-dimensional special case of Stokes' theorem.
Theorem
Let be a positively orient ...
implies that for
:
Thus if in and satisfies the Neumann boundary conditions, , and so is constant in .
Hence the Neumann problem has a unique solution up to adding constants.
Consider the Hermitian form on defined by
:
Since is in duality with , there is a unique element in such that
:
The map is an isometry of onto , so in particular is bounded.
In fact
:
So
:
On the other hand, any in defines a bounded conjugate-linear form on sending to . By the
Riesz–Fischer theorem
In mathematics, the Riesz–Fischer theorem in real analysis is any of a number of closely related results concerning the properties of the space ''L''2 of square integrable functions. The theorem was proven independently in 1907 by Frigyes Riesz ...
, there exists such that
:
Hence and so is surjective. Define a bounded linear operator on by
:
where is the map , a compact operator, and is the map , its adjoint, so also compact.
The operator has the following properties:
* is a contraction since it is a composition of contractions
* is compact, since and are compact by Rellich's theorem
* is self-adjoint, since if , they can be written with so
::
* has positive spectrum and kernel , for
::
:and implies and hence .
*There is a complete orthonormal basis of consisting of eigenfunctions of . Thus
::
:with and decreasing to .
*The eigenfunctions all lie in since the image of lies in .
*The are eigenfunctions of with
::
:Thus are non-negative and increase to .
*The eigenvalue occurs with multiplicity one and corresponds to the constant function. For if satisfies , then
::
:so is constant.
Regularity for Neumann problem
Weak solutions are strong solutions
The first main regularity result shows that a weak solution expressed in terms of the operator and the Dirichlet form is a strong solution in the classical sense, expressed in terms of the Laplacian and the Neumann boundary conditions. Thus if with , then , satisfies and . Moreover, for some constant independent of ,
:
Note that
:
since
:
Take a decomposition with supported in and supported in a collar of the boundary.
The operator is characterized by
:
Then
:
so that
:
The function and are treated separately, being essentially subject to usual elliptic regularity considerations for interior points while requires special treatment near the boundary using difference quotients. Once the strong properties are established in terms of and the Neumann boundary conditions, the "bootstrap" regularity results can be proved exactly as for the Dirichlet problem.
Interior estimates
The function lies in where is a region with closure in . If and
:
By continuity the same holds with replaced by and hence . So
:
Hence regarding as an element of , . Hence . Since for , we have . Moreover,
:
so that
:
Boundary estimates
The function is supported in a collar contained in a tubular neighbourhood of the boundary. The difference quotients can be formed for the flow and lie in , so the first inequality is applicable:
:
The commutators are uniformly bounded as operators from to . This is equivalent to checking the inequality
:
for , smooth functions on a collar. This can be checked directly on an annulus, using invariance of Sobolev spaces under dffeomorphisms and the fact that for the annulus the commutator of with a differential operator is obtained by applying the difference operator to the coefficients after having applied to the function:
:
Hence the difference quotients are uniformly bounded, and therefore with
:
Hence and satisfies a similar inequality to :
:
Let be the orthogonal vector field. As for the Dirichlet problem, to show that , it suffices to show that .
To check this, it is enough to show that . As before
:
are vector fields. On the other hand, for , so that and define the same distribution on . Hence
:
Since the terms on the right hand side are pairings with functions in , the regularity criterion shows that . Hence since both terms lie in and have the same inner products with 's.
Moreover, the inequalities for show that
:
Hence
:
It follows that . Moreover,
:
Neumann boundary conditions
Since , Green's theorem is applicable by continuity. Thus for ,
:
Hence the Neumann boundary conditions are satisfied:
:
where the left hand side is regarded as an element of and hence .
Regularity of strong solutions
The main result here states that if and , then and
:
for some constant independent of .
Like the corresponding result for the Dirichlet problem, this is proved by induction on . For , is also a weak solution of the Neumann problem so satisfies the estimate above for . The Neumann boundary condition can be written
:
Since commutes with the vector field corresponding to the period flow , the inductive method of proof used for the Dirichlet problem works equally well in this case: for the difference quotients preserve the boundary condition when expressed in terms of .
Smoothness of eigenfunctions
It follows by induction from the regularity theorem for the Neumann problem that the eigenfunctions of in lie in . Moreover, any solution of with in and in must have in . In both cases by the vanishing properties, the normal derivatives of the eigenfunctions and vanish on .
Solving the associated Neumann problem
The method above can be used to solve the associated Neumann boundary value problem:
:
By Borel's lemma is the restriction of a function . Let be a smooth function such that near the boundary. Let be the solution of with . Then solves the boundary value problem.
Notes
References
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*{{citation, title=An Introduction to Partial Differential Equations, volume =13, series= Texts in Applied Mathematics, first1=Michael, last1= Renardy, first2= Robert C., last2= Rogers, edition=2nd , publisher=Springer, year= 2004, isbn=0-387-00444-0
Partial differential equations
Harmonic analysis
Operator theory
Functional analysis