Misplaced Pages

Elliptic operator

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Elliptic partial differential operator) Type of differential operator
This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (September 2024) (Learn how and when to remove this message)
A solution to Laplace's equation defined on an annulus. The Laplace operator is the most famous example of an elliptic operator.

In the theory of partial differential equations, elliptic operators are differential operators that generalize the Laplace operator. They are defined by the condition that the coefficients of the highest-order derivatives be positive, which implies the key property that the principal symbol is invertible, or equivalently that there are no real characteristic directions.

Elliptic operators are typical of potential theory, and they appear frequently in electrostatics and continuum mechanics. Elliptic regularity implies that their solutions tend to be smooth functions (if the coefficients in the operator are smooth). Steady-state solutions to hyperbolic and parabolic equations generally solve elliptic equations.

Definitions

Let L {\displaystyle L} be a linear differential operator of order m on a domain Ω {\displaystyle \Omega } in R given by L u = | α | m a α ( x ) α u {\displaystyle Lu=\sum _{|\alpha |\leq m}a_{\alpha }(x)\partial ^{\alpha }u} where α = ( α 1 , , α n ) {\displaystyle \alpha =(\alpha _{1},\dots ,\alpha _{n})} denotes a multi-index, and α u = 1 α 1 n α n u {\displaystyle \partial ^{\alpha }u=\partial _{1}^{\alpha _{1}}\cdots \partial _{n}^{\alpha _{n}}u} denotes the partial derivative of order α i {\displaystyle \alpha _{i}} in x i {\displaystyle x_{i}} .

Then L {\displaystyle L} is called elliptic if for every x in Ω {\displaystyle \Omega } and every non-zero ξ {\displaystyle \xi } in R, | α | = m a α ( x ) ξ α 0 , {\displaystyle \sum _{|\alpha |=m}a_{\alpha }(x)\xi ^{\alpha }\neq 0,} where ξ α = ξ 1 α 1 ξ n α n {\displaystyle \xi ^{\alpha }=\xi _{1}^{\alpha _{1}}\cdots \xi _{n}^{\alpha _{n}}} .

In many applications, this condition is not strong enough, and instead a uniform ellipticity condition may be imposed for operators of order m = 2k: ( 1 ) k | α | = 2 k a α ( x ) ξ α > C | ξ | 2 k , {\displaystyle (-1)^{k}\sum _{|\alpha |=2k}a_{\alpha }(x)\xi ^{\alpha }>C|\xi |^{2k},} where C is a positive constant. Note that ellipticity only depends on the highest-order terms.

A nonlinear operator L ( u ) = F ( x , u , ( α u ) | α | m ) {\displaystyle L(u)=F\left(x,u,\left(\partial ^{\alpha }u\right)_{|\alpha |\leq m}\right)} is elliptic if its linearization is; i.e. the first-order Taylor expansion with respect to u and its derivatives about any point is an elliptic operator.

Example 1
The negative of the Laplacian in R given by Δ u = i = 1 d i 2 u {\displaystyle -\Delta u=-\sum _{i=1}^{d}\partial _{i}^{2}u} is a uniformly elliptic operator. The Laplace operator occurs frequently in electrostatics. If ρ is the charge density within some region Ω, the potential Φ must satisfy the equation Δ Φ = 4 π ρ . {\displaystyle -\Delta \Phi =4\pi \rho .}
Example 2
Given a matrix-valued function A(x) which uniformly positive definite for every x, having components a, the operator L u = i ( a i j ( x ) j u ) + b j ( x ) j u + c u {\displaystyle Lu=-\partial _{i}\left(a^{ij}(x)\partial _{j}u\right)+b^{j}(x)\partial _{j}u+cu} is elliptic. This is the most general form of a second-order divergence form linear elliptic differential operator. The Laplace operator is obtained by taking A = I. These operators also occur in electrostatics in polarized media.
Example 3
For p a non-negative number, the p-Laplacian is a nonlinear elliptic operator defined by L ( u ) = i = 1 d i ( | u | p 2 i u ) . {\displaystyle L(u)=-\sum _{i=1}^{d}\partial _{i}\left(|\nabla u|^{p-2}\partial _{i}u\right).} A similar nonlinear operator occurs in glacier mechanics. The Cauchy stress tensor of ice, according to Glen's flow law, is given by τ i j = B ( k , l = 1 3 ( l u k ) 2 ) 1 3 1 2 ( j u i + i u j ) {\displaystyle \tau _{ij}=B\left(\sum _{k,l=1}^{3}\left(\partial _{l}u_{k}\right)^{2}\right)^{-{\frac {1}{3}}}\cdot {\frac {1}{2}}\left(\partial _{j}u_{i}+\partial _{i}u_{j}\right)} for some constant B. The velocity of an ice sheet in steady state will then solve the nonlinear elliptic system j = 1 3 j τ i j + ρ g i i p = Q , {\displaystyle \sum _{j=1}^{3}\partial _{j}\tau _{ij}+\rho g_{i}-\partial _{i}p=Q,} where ρ is the ice density, g is the gravitational acceleration vector, p is the pressure and Q is a forcing term.

Elliptic regularity theorems

Let L be an elliptic operator of order 2k with coefficients having 2k continuous derivatives. The Dirichlet problem for L is to find a function u, given a function f and some appropriate boundary values, such that Lu = f and such that u has the appropriate boundary values and normal derivatives. The existence theory for elliptic operators, using Gårding's inequality, Lax–Milgram lemma and Fredholm alternative, states the sufficient condition for a weak solution u to exist in the Sobolev space H.

For example, for a Second-order Elliptic operator as in Example 2,

  • There is a number γ>0 such that for each μ>γ, each f L 2 ( U ) {\displaystyle f\in L^{2}(U)} , there exists a unique solution u H 0 1 ( U ) {\displaystyle u\in H_{0}^{1}(U)} of the boundary value problem
    L u + μ u = f  in  U , u = 0  on  U {\displaystyle Lu+\mu u=f{\text{ in }}U,u=0{\text{ on }}\partial U} , which is based on Lax-Milgram lemma.
  • Either (a) for any f L 2 ( U ) {\displaystyle f\in L^{2}(U)} , L u = f  in  U , u = 0  on  U {\displaystyle Lu=f{\text{ in }}U,u=0{\text{ on }}\partial U} (1) has a unique solution, or (b) L u = 0  in  U , u = 0  on  U {\displaystyle Lu=0{\text{ in }}U,u=0{\text{ on }}\partial U} has a solution u 0 {\displaystyle u\not \equiv 0} , which is based on the property of compact operators and Fredholm alternative.

This situation is ultimately unsatisfactory, as the weak solution u might not have enough derivatives for the expression Lu to be well-defined in the classical sense.

The elliptic regularity theorem guarantees that, provided f is square-integrable, u will in fact have 2k square-integrable weak derivatives. In particular, if f is infinitely-often differentiable, then so is u.

For L as in Example 2,

  • Interior regularity: If m is a natural number, a i j , b j , c C m + 1 ( U ) , f H m ( U ) {\displaystyle a^{ij},b^{j},c\in C^{m+1}(U),f\in H^{m}(U)} (2) , u H 0 1 ( U ) {\displaystyle u\in H_{0}^{1}(U)} is a weak solution to (1), then for any open set V in U with compact closure, u H m + 2 ( V ) C ( f H m ( U ) + u L 2 ( U ) ) {\displaystyle \|u\|_{H^{m+2}(V)}\leq C(\|f\|_{H^{m}(U)}+\|u\|_{L^{2}(U)})} (3), where C depends on U, V, L, m, per se u H l o c m + 2 ( U ) {\displaystyle u\in H_{loc}^{m+2}(U)} , which also holds if m is infinity by Sobolev embedding theorem.
  • Boundary regularity: (2) together with the assumption that U {\displaystyle \partial U} is C m + 2 {\displaystyle C^{m+2}} indicates that (3) still holds after replacing V with U, i.e. u H m + 2 ( U ) {\displaystyle u\in H^{m+2}(U)} , which also holds if m is infinity.

Any differential operator exhibiting this property is called a hypoelliptic operator; thus, every elliptic operator is hypoelliptic. The property also means that every fundamental solution of an elliptic operator is infinitely differentiable in any neighborhood not containing 0.

As an application, suppose a function f {\displaystyle f} satisfies the Cauchy–Riemann equations. Since the Cauchy-Riemann equations form an elliptic operator, it follows that f {\displaystyle f} is smooth.

Properties

For L as in Example 2 on U, which is an open domain with C boundary, then there is a number γ>0 such that for each μ>γ, L + μ I : H 0 1 ( U ) H 0 1 ( U ) {\displaystyle L+\mu I:H_{0}^{1}(U)\rightarrow H_{0}^{1}(U)} satisfies the assumptions of Lax-Milgram lemma.

  • Invertibility: For each μ>γ, L + μ I : L 2 ( U ) L 2 ( U ) {\displaystyle L+\mu I:L^{2}(U)\rightarrow L^{2}(U)} admits a compact inverse.
  • Eigenvalues and eigenvectors: If A is symmetric, b,c are zero, then (1) Eigenvalues of L, are real, positive, countable, unbounded (2) There is an orthonormal basis of L(U) composed of eigenvectors of L.
  • Generates a semigroup on L(U): -L generates a semigroup { S ( t ) ; t 0 } {\displaystyle \{S(t);t\geq 0\}} of bounded linear operators on L(U) s.t. d d t S ( t ) u 0 = L S ( t ) u 0 , S ( t ) e γ t {\displaystyle {\frac {d}{dt}}S(t)u_{0}=-LS(t)u_{0},\|S(t)\|\leq e^{\gamma t}} in the norm of L(U), for every u 0 L 2 ( U ) {\displaystyle u_{0}\in L^{2}(U)} , by Hille-Yosida theorem.

General definition

Let D {\displaystyle D} be a (possibly nonlinear) differential operator between vector bundles of any rank. Take its principal symbol σ ξ ( D ) {\displaystyle \sigma _{\xi }(D)} with respect to a one-form ξ {\displaystyle \xi } . (Basically, what we are doing is replacing the highest order covariant derivatives {\displaystyle \nabla } by vector fields ξ {\displaystyle \xi } .)

We say D {\displaystyle D} is weakly elliptic if σ ξ ( D ) {\displaystyle \sigma _{\xi }(D)} is a linear isomorphism for every non-zero ξ {\displaystyle \xi } .

We say D {\displaystyle D} is (uniformly) strongly elliptic if for some constant c > 0 {\displaystyle c>0} , ( [ σ ξ ( D ) ] ( v ) , v ) c v 2 {\displaystyle \left((v),v\right)\geq c\|v\|^{2}}

for all ξ = 1 {\displaystyle \|\xi \|=1} and all v {\displaystyle v} .

The definition of ellipticity in the previous part of the article is strong ellipticity. Here ( , ) {\displaystyle (\cdot ,\cdot )} is an inner product. Notice that the ξ {\displaystyle \xi } are covector fields or one-forms, but the v {\displaystyle v} are elements of the vector bundle upon which D {\displaystyle D} acts.

The quintessential example of a (strongly) elliptic operator is the Laplacian (or its negative, depending upon convention). It is not hard to see that D {\displaystyle D} needs to be of even order for strong ellipticity to even be an option. Otherwise, just consider plugging in both ξ {\displaystyle \xi } and its negative. On the other hand, a weakly elliptic first-order operator, such as the Dirac operator can square to become a strongly elliptic operator, such as the Laplacian. The composition of weakly elliptic operators is weakly elliptic.

Weak ellipticity is nevertheless strong enough for the Fredholm alternative, Schauder estimates, and the Atiyah–Singer index theorem. On the other hand, we need strong ellipticity for the maximum principle, and to guarantee that the eigenvalues are discrete, and their only limit point is infinity.

See also

Notes

  1. Note that this is sometimes called strict ellipticity, with uniform ellipticity being used to mean that an upper bound exists on the symbol of the operator as well. It is important to check the definitions the author is using, as conventions may differ. See, e.g., Evans, Chapter 6, for a use of the first definition, and Gilbarg and Trudinger, Chapter 3, for a use of the second.
  2. See Evans, Chapter 6-7, for details.

References

External links

Categories: