Elliptic partial differential equation

In mathematics, an elliptic partial differential equation is a type of partial differential equation (PDE). In mathematical modeling, elliptic PDEs are frequently used to model steady states, unlike parabolic PDE and hyperbolic PDE which generally model phenomena that change in time. They are also important in pure mathematics, where they are fundamental to various fields of research such as differential geometry and optimal transport.

Definition

Elliptic differential equations appear in many different contexts and levels of generality.

First consider a second-order linear PDE in two variables, written in the form $$A{u}_{xx}+2B{u}_{xy}+C{u}_{yy}+D{u}_x+E{u}_y+Fu+G=0,$$ where Template:Math, Template:Math, Template:Math, Template:Math, Template:Math, Template:Math, and Template:Math are functions of Template:Math and Template:Math, using subscript notation for the partial derivatives. The PDE is called elliptic if $$B^2-AC<0,$$ with this naming convention inspired by the equation for a planar ellipse. Equations with $B^2-AC=0$ are termed parabolic while those with $B^2-AC>0$ are hyperbolic.

For a general linear second-order PDE, the "unknown" function Template:Mvar can be a function of any number Template:Math of independent variables; the equation is of the form $$\sum _{i=1}^n\sum _{j=1}^n{a}_{i,j}(x_1,\dots ,x_n)u_{x_i{x}_j}+\sum _{i=1}^n{b}_i(x_1,\dots ,x_n)u_{x_i}+c(x_1,\dots ,x_n)u=f(x_1,\dots ,x_n).$$ where Template:Math, Template:Math, Template:Mvar and Template:Mvar are functions defined on the domain subject to the symmetry Template:Math. This equation is called elliptic if, when Template:Mvar is viewed as a function on the domain valued in the space of Template:Math symmetric matrices, all of the eigenvalues are greater than some set positive number. Equivalently, this means that there is a positive number Template:Mvar such that $$\sum _{i=1}^n\sum _{j=1}^n{a}_{ij}(x_1,\dots ,x_n){\xi }_i{\xi }_j\ge \theta ({\xi }_1^2+\cdots +{\xi }_n^2),$$ for any point Template:Math in the domain and any real numbers Template:Math.Template:SfnmTemplate:Sfn

The simplest example of a second-order linear elliptic PDE is the Laplace equation, in which Template:Math is zero if Template:Math and is one otherwise, and where Template:Math. The Poisson equation is a slightly more general second-order linear elliptic PDE, in which Template:Mvar is not required to vanish. For both of these equations, the ellipticity constant Template:Mvar can be taken to be Template:Math.

The terminology elliptic partial differential equation is not used consistently throughout the literature. What is called "elliptic" by some authors is called strictly elliptic or uniformly elliptic by others.

Nonlinear and higher-order equations

Ellipticity can also be formulated for much more general classes of equations. For the most general second-order PDE, which is of the form

$F(D^{2}u,Du,u,x_1,\dots ,x_n)=0$

for some given function Template:Mvar, ellipticity is defined by linearizing the equation and applying the above linear definition. Since linearization is done at a particular function Template:Mvar, this means that ellipticity of a nonlinear second-order PDE depends not only on the equation itself but also on the solutions under consideration. For example, in the simplest kind of Monge–Ampère equation, the determinant of the hessian matrix of a function is prescribed:

$\det D^{2}u=f.$

As follows from Jacobi's formula for the derivative of a determinant, this equation is elliptic if Template:Mvar is a positive function and solutions satisfy the constraint of being uniformly convex.Template:Sfnm

There are also higher-order elliptic PDE, the simplest example being the fourth-order biharmonic equation.Template:Sfnm Even more generally, there is an important class of elliptic systems which consist of coupled partial differential equations for multiple 'unknown' functions.Template:Sfnm For example, the Cauchy–Riemann equations from complex analysis can be viewed as a first-order elliptic system for a pair of two-variable functions.Template:Sfnm

Moreover, the class of elliptic PDE (of any order, including systems) is subject to various notions of weak solutions, i.e., reformulating the above equations in such a way that allows for solutions to have various irregularities (e.g. non-differentiability, singularities or discontinuities) while still adhering to the laws of physics.Template:Sfnm Additionally, these type of solutions are also important in variational calculus, where the direct method often produces weak solutions of elliptic systems of Euler equations.Template:Sfnm

Canonical form

Consider a second-order elliptic partial differential equation

$A(x,y)u_{xx}+2B(x,y)u_{xy}+C(x,y)u_{yy}+f(u_x,u_y,u,x,y)=0$

for a two-variable function Template:Math. This equation is linear in the "leading-order terms" but allows nonlinear expressions involving the function values and their first derivatives; this is sometimes called a quasilinear equation.

A canonical form asks for a transformation Template:Math and Template:Math of the domain so that, when Template:Mvar is viewed as a function of Template:Mvar and Template:Mvar, the above equation takes the form

$u_{ww}+u_{zz}+F(u_w,u_z,u,w,z)=0$

for some new function Template:Mvar. The existence of such a transformation can be established locally if Template:Mvar, Template:Mvar, and Template:Mvar are real-analytic functions and, with more elaborate work, even if they are only continuously differentiable. Locality means that the necessary coordinate transformations may fail to be defined on the entire domain of Template:Mvar, although they can be established in some small region surrounding any particular point of the domain.Template:Sfnm

Formally establishing the existence of such transformations uses the existence of solutions to the Beltrami equation. From the perspective of differential geometry, the existence of a canonical form is equivalent to the existence of isothermal coordinates for the associated Riemannian metric

$A(x,y)d{x}^2+2B(x,y)dxdy+C(x,y)d{y}^2$

on the domain. (The ellipticity condition for the PDE, namely the positivity of the function Template:Math, is what ensures that either this tensor or its negation is indeed a Riemannian metric.) Generally, for second-order quasilinear elliptic partial differential equations for functions of more than two variables, a canonical form does not exist. This corresponds to the fact that, although isothermal coordinates generally exist for Riemannian metrics in two dimensions, they only exist for very particular Riemannian metrics in higher dimensions.Template:Sfnm

Characteristics and regularity

For the general second-order linear PDE, characteristics are defined as the null directions for the associated tensorTemplate:Sfn

$\sum _{i=1}^n\sum _{j=1}^n{a}_{i,j}(x_1,\dots ,x_n)d{x}^{i}d{x}^j,$

called the principal symbol. Using the technology of the wave front set, characteristics are significant in understanding how irregular points of Template:Mvar propagate to the solution Template:Mvar of the PDE. Informally, the wave front set of a function consists of the points of non-smoothness, in addition to the directions in frequency space causing the lack of smoothness. It is a fundamental fact that the application of a linear differential operator with smooth coefficients can only have the effect of removing points from the wave front set.Template:Sfnm However, all points of the original wave front set (and possibly more) are recovered by adding back in the (real) characteristic directions of the operator.Template:Sfnm

In the case of a linear elliptic operator Template:Mvar with smooth coefficients, the principal symbol is a Riemannian metric and there are no real characteristic directions. According to the previous paragraph, it follows that the wave front set of a solution Template:Mvar coincides exactly with that of Template:Math. This sets up a basic regularity theorem, which says that if Template:Math is smooth (so that its wave front set is empty) then the solution Template:Mvar is smooth as well. More generally, the points where Template:Mvar fails to be smooth coincide with the points where Template:Math is not smooth.Template:Sfnm This regularity phenomena is in sharp contrast with, for example, hyperbolic PDE in which discontinuities can form even when all the coefficients of an equation are smooth.

Solutions of elliptic PDEs are naturally associated with time-independent solutions of parabolic PDEs or hyperbolic PDEs. For example, a time-independent solution of the heat equation solves Laplace's equation. That is, if parabolic and hyperbolic PDEs are associated with modeling dynamical systems then the solutions of elliptic PDEs are associated with steady states. Informally, this is reflective of the above regularity theorem, as steady states are generally smoothed out versions of truly dynamical solutions. However, PDE used in modeling are often nonlinear and the above regularity theorem only applies to linear elliptic equations; moreover, the regularity theory for nonlinear elliptic equations is much more subtle, with solutions not always being smooth.

See also

• Elliptic boundary value problem
• Elliptic operator
• Hyperbolic partial differential equation
• Parabolic partial differential equation
• Maximum principle (property of solutions)
• Sobolev space

References

• Template:Cite book
• Template:Cite journal
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book

Further reading

• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book

External links

• Template:Springer
• Template:Springer
• Template:MathWorld