Fact-checked by Grok 2 weeks ago

Elliptic partial differential equation

An elliptic partial differential equation (PDE) is a type of second-order linear PDE of the form \sum_{i,j=1}^n a_{ij}(x) \frac{\partial^2 u}{\partial x_i \partial x_j} + \sum_{i=1}^n b_i(x) \frac{\partial u}{\partial x_i} + c(x) u = f(x), where the coefficient matrix (a_{ij}) is symmetric and uniformly elliptic, meaning there exists \theta > 0 such that \sum_{i,j} a_{ij} \xi_i \xi_j \geq \theta |\xi|^2 for all \xi \in \mathbb{R}^n.^[1] These equations model steady-state or equilibrium phenomena in physical systems, such as electrostatic potentials or steady heat distributions, and are distinguished from parabolic and hyperbolic PDEs by the lack of real characteristics, leading to well-posed boundary value problems rather than initial value problems.^[1] The classification of second-order PDEs into elliptic, parabolic, and hyperbolic types originates from the discriminant of the principal part, analogous to conic sections; for the equation a u_{xx} + 2b u_{xy} + c u_{yy} + \ lower\ order\ terms = 0 in two variables, it is elliptic if b^2 - ac < 0.^[2] This classification was formalized by Jacques Hadamard in 1923, building on earlier work, and emphasizes the elliptic case's connection to the Laplace operator, which has no real characteristics except the trivial one.^[2] Prominent examples include the Laplace equation \Delta u = 0, which describes harmonic functions and arises in potential theory for irrotational, incompressible fluid flow or gravitational fields, and the Poisson equation \Delta u = f, a inhomogeneous variant modeling sources like charge distributions in electrostatics.^[1] More generally, divergence-form elliptic equations like -\nabla \cdot (A(x) \nabla u) = f appear in diffusion processes at equilibrium, with uniform ellipticity ensuring the matrix A(x) is positive definite.^[1] Key properties of solutions to elliptic PDEs include the maximum principle, which states that for the Laplace equation on a bounded domain, the maximum of a non-constant harmonic function is attained on the boundary, implying no interior maxima and aiding in uniqueness proofs.^[1] Regularity theory further asserts that weak solutions in Sobolev spaces H^1(\Omega) are classically smooth interior to the domain if coefficients are sufficiently regular, with interior C^{2,\alpha} estimates from the De Giorgi-Nash-Moser theory developed in the 1950s-1960s.^[1] These features underpin existence and uniqueness via variational methods, such as the Lax-Milgram theorem for Dirichlet problems.^[1] Historically, elliptic PDEs trace to Pierre-Simon Laplace's 1782 work on gravitational potentials, with the Dirichlet problem posed by Bernhard Riemann in 1857 and solved variationally by Henri Poincaré in 1890; David Hilbert revived the Dirichlet principle in 1900, spurring modern developments like Sergei Bernstein's 1904 analyticity results for two-dimensional nonlinear cases.^[2] Applications span physics (e.g., quantum mechanics via the Schrödinger equation in steady state), geometry (e.g., minimal surfaces), and engineering (e.g., steady-state heat conduction), with ongoing research in nonlinear and fully nonlinear variants like the Monge-Ampère equation.^[2]

Introduction

Overview and motivation

Elliptic partial differential equations (PDEs) form a class of equations that model steady-state or equilibrium phenomena in physical systems, where the solution describes a time-independent balance of forces or quantities. Unlike hyperbolic PDEs, which govern wave propagation with finite speed and sharp fronts, or parabolic PDEs, which describe diffusive processes evolving over time, elliptic PDEs lack a preferred direction of information propagation and typically yield smooth solutions throughout the domain. This classification arises from the mathematical structure of the equations, particularly for second-order linear cases, and reflects their role in applications like potential theory and stationary flows.^[3] Prominent examples include Laplace's equation, \nabla^2 u = 0, which arises in steady-state heat conduction without sources, where u represents the temperature in a homogeneous medium at equilibrium. Another key instance is Poisson's equation, \nabla^2 u = -\rho / \epsilon_0, modeling the electrostatic potential u due to a charge distribution \rho in vacuum, with \epsilon_0 as the permittivity of free space. These equations illustrate how elliptic PDEs capture balanced, non-transient behaviors in physics, such as thermal equilibrium or irrotational fluid flow.^[4]^[3] The general form of a second-order linear elliptic PDE in two variables is a u_{xx} + 2b u_{xy} + c u_{yy} + d u_x + e u_y + f u = g(x,y), where the principal part satisfies the ellipticity condition b^2 - ac < 0. This discriminant condition ensures the equation's elliptic nature, distinguishing it from parabolic (b^2 - ac = 0) and hyperbolic (b^2 - ac > 0) cases by analogy to conic sections. Lower-order terms and the right-hand side g influence the solution but do not alter the classification.^[5] A brief historical note traces the origins of this classification to early 19th-century developments, with Joseph Fourier's work on the heat equation in his 1822 treatise Théorie analytique de la chaleur exemplifying the parabolic type and inspiring the broader categorization of PDEs into elliptic, parabolic, and hyperbolic classes based on their physical and mathematical behaviors.^[2]

Historical context

The origins of elliptic partial differential equations trace to 18th- and 19th-century developments in physics and mathematics, particularly through equations modeling steady-state phenomena. The classification of second-order PDEs into elliptic, parabolic, and hyperbolic types, analogous to conic sections, was formalized by Jacques Hadamard in 1923.^[6] In the early 19th century, Joseph Fourier advanced the field through his 1822 treatise Théorie Analytique de la Chaleur, where his analysis of heat conduction in steady states led to the formulation of Laplace's equation as a prototypical elliptic PDE.^[7] Building on this, Carl Friedrich Gauss contributed to potential theory in the 1830s, particularly with his 1839 work Allgemeine Theorie des Erdmagnetismus, which applied elliptic equations to model magnetic fields via scalar potentials.^[8] Concurrently, George Green independently developed key aspects of potential theory in his 1828 self-published essay An Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism, introducing Green's functions and integral representations essential for solving elliptic boundary value problems.^[9] The 20th century saw significant progress in the analytic properties of solutions to elliptic PDEs, spurred by David Hilbert's 19th problem posed at the 1900 International Congress of Mathematicians in Paris, which asked whether solutions to certain elliptic equations with analytic coefficients are themselves analytic.^[10] Sergei Bernstein provided a partial affirmative solution in 1904, proving analyticity for twice continuously differentiable solutions of two-dimensional elliptic equations in his doctoral thesis published in Mathematische Annalen. Later advancements in elliptic regularity theory were driven by Sergei Sobolev's introduction of Sobolev spaces in the 1930s, which enabled weak formulations and embedding theorems crucial for establishing higher regularity of solutions to elliptic problems. Lars Hörmander further refined interior regularity results in the mid-20th century, notably in his 1958 paper demonstrating that solutions to elliptic equations with sufficiently smooth coefficients inherit the regularity of the data.^[11]

Classification and definition

Linear second-order equations

Linear second-order partial differential equations (PDEs) in n variables take the general form

\sum_{i,j=1}^n a_{ij}(x) \frac{\partial^2 u}{\partial x_i \partial x_j} + \sum_{i=1}^n b_i(x) \frac{\partial u}{\partial x_i} + c(x) u = f(x),

where the coefficients a_{ij}, b_i, c, and f are given functions defined on a domain \Omega \subset \mathbb{R}^n, and the matrix (a_{ij}) is symmetric, i.e., a_{ij} = a_{ji}.^[4] This equation models steady-state phenomena, such as electrostatic potentials or equilibrium temperatures, through a balance of diffusion-like second-order terms and lower-order effects.^[12] The principal part of the operator is the second-order term \sum_{i,j=1}^n a_{ij}(x) \frac{\partial^2 u}{\partial x_i \partial x_j}, which governs the highest-order behavior and determines the type of the PDE.^[4] The associated principal symbol is the quadratic form \sum_{i,j=1}^n a_{ij}(x) \xi_i \xi_j, where \xi = (\xi_1, \dots, \xi_n) \in \mathbb{R}^n is the dual variable.^[12] A linear second-order PDE is elliptic at a point x \in \Omega if the symmetric matrix A(x) = (a_{ij}(x)) is positive definite, meaning all its eigenvalues are positive or, equivalently, \sum_{i,j=1}^n a_{ij}(x) \xi_i \xi_j > 0 for all \xi \neq 0.^[4]^[12] In two dimensions, for the equation a(x,y) u_{xx} + 2b(x,y) u_{xy} + c(x,y) u_{yy} + \ lower\ terms = f, the matrix is A = \begin{pmatrix} a & b \\ b & c \end{pmatrix}, and positive definiteness holds if the trace a + c > 0 and the determinant ac - b^2 > 0, which simplifies to the discriminant condition b^2 - ac < 0.^[4] For example, Laplace's equation \Delta u = u_{xx} + u_{yy} = 0 has a = 1, b = 0, c = 1, so the discriminant is $0 - 1 \cdot 1 = -1 < 0, confirming it is elliptic.^[4] In contrast, the one-dimensional wave equation u_{tt} - u_{xx} = 0 (viewed in space-time variables t, x) has coefficients a = -1 for x and c = 1 for t with b = 0, yielding discriminant $0 - (-1) \cdot 1 = 1 > 0, making it hyperbolic.^[4] A stronger condition, uniform ellipticity, requires the existence of constants \lambda, \Lambda > 0 (independent of x) such that

\lambda |\xi|^2 \leq \sum_{i,j=1}^n a_{ij}(x) \xi_i \xi_j \leq \Lambda |\xi|^2

for all x \in \Omega and all \xi \in \mathbb{R}^n \setminus \{0\}.^[4] This bounds the eigenvalues of A(x) away from zero and infinity, ensuring consistent elliptic behavior across the domain and facilitating proofs of existence, uniqueness, and regularity for solutions.^[12] For Laplace's equation, uniform ellipticity holds with \lambda = \Lambda = 1.^[4]

General and nonlinear cases

The classification of partial differential equations (PDEs) as elliptic extends beyond the linear second-order case to more general settings, including nonlinear equations and higher-order operators. For nonlinear PDEs of the form F(x, u, Du, D^2 u) = 0, where F is a smooth function depending on the position x, the solution u, its first derivatives Du, and second derivatives D^2 u, the equation is considered elliptic at a solution u if the associated linearized operator is elliptic. Specifically, the linearization L_v = \sum F_{u_i} v_{x_i} + \sum F_{u_{ij}} v_{x_i x_j} + \text{lower-order terms}, where the coefficients are evaluated along the solution, must satisfy the uniform ellipticity condition that there exist positive constants \lambda, \Lambda > 0 such that \lambda |\xi|^2 \leq \sum a_{ij}(x) \xi_i \xi_j \leq \Lambda |\xi|^2 for all x in the domain and all real vectors \xi, with a_{ij} = F_{u_{ij}}. A prominent subclass is the quasilinear elliptic PDE, typically written as \operatorname{div}(A(x, u, Du) Du) = f(x, u, Du), where A is a matrix-valued function. This equation is elliptic if A(x, u, p) is uniformly positive definite for all arguments in the relevant domain, meaning its eigenvalues are bounded below by a positive constant \lambda > 0, ensuring the principal part behaves like a uniformly elliptic operator. This condition guarantees that the equation inherits key analytic properties from its linear counterparts, such as regularity of solutions under suitable assumptions. For linear PDEs of higher even order, say order $2m, the operator P = \sum_{|\alpha| \leq 2m} a_\alpha(x) D^\alpha is elliptic if its principal symbol p_{2m}(x, \xi) = \sum_{|\alpha| = 2m} a_\alpha(x) (i \xi)^\alpha has no real zeros except at \xi = 0, or more precisely, p_{2m}(x, \xi) \neq 0 for all \xi \in \mathbb{R}^n \setminus \{0\} and x in the domain. This generalizes the second-order case, where the symbol is a quadratic form without real characteristics. A canonical example is the bi-Laplace equation \Delta^2 u = 0, a fourth-order (m=2) linear elliptic PDE whose principal symbol is |\xi|^4, which vanishes only at \xi = 0 and thus satisfies the ellipticity condition uniformly.

Fundamental properties

Maximum principle

The maximum principle is a cornerstone of elliptic partial differential equation theory, asserting that solutions to certain elliptic inequalities attain their extrema on the boundary of the domain rather than in the interior. This property, which holds under suitable ellipticity conditions, provides essential bounds for solutions and underpins uniqueness results in boundary value problems. For the Laplace equation \Delta u = 0 in a bounded domain \Omega \subset \mathbb{R}^n, the weak maximum principle states that a continuous solution u satisfies \max_{\overline{\Omega}} u = \max_{\partial \Omega} u and \min_{\overline{\Omega}} u = \min_{\partial \Omega} u. This extends to subharmonic functions satisfying \Delta u \geq 0, where the maximum is still attained on the boundary, reflecting the "mean value property" that subharmonic functions lie below their averages over balls. The strong maximum principle strengthens this by asserting that if \Delta u \geq 0 and u achieves its maximum at an interior point of \Omega, then u must be constant throughout \Omega. These principles generalize to linear second-order elliptic operators of the form

Lu = \sum_{i,j=1}^n a_{ij} \frac{\partial^2 u}{\partial x_i \partial x_j} + \sum_{i=1}^n b_i \frac{\partial u}{\partial x_i} + c u,

where the matrix (a_{ij}) is uniformly elliptic (i.e., there exist positive constants \lambda, \Lambda such that \lambda |\xi|^2 \leq \sum a_{ij} \xi_i \xi_j \leq \Lambda |\xi|^2 for all \xi \in \mathbb{R}^n) and c \leq 0. In this setting, if Lu \geq 0 in \Omega, the weak maximum principle holds: \max_{\overline{\Omega}} u = \max_{\partial \Omega} u, assuming u is continuous up to the boundary. The strong version similarly prohibits non-constant solutions from attaining an interior maximum. A key refinement is the Hopf boundary point lemma, which applies when u achieves its maximum at a boundary point x_0 \in \partial \Omega with \Omega lying on one side of the boundary (e.g., C^2 boundary). Under the conditions Lu \geq 0, c \leq 0, and non-constancy of u, there exists a direction (the outward normal) in which the directional derivative at x_0 is strictly positive, implying u increases away from the boundary near x_0. This lemma, originally due to Hopf, ensures strict boundary control and is crucial for handling cases where equality might otherwise hold. Proofs of these results typically rely on the mean value property for harmonic functions, extended via perturbation or barrier arguments for general elliptic operators. For harmonic functions (\Delta u = 0), the mean value property states that u(x) = \frac{1}{|B_r(x)|} \int_{B_r(x)} u(y) \, dy for balls B_r(x) \subset \Omega, implying that an interior maximum would require u to be constant by Jensen's inequality or direct averaging. For general Lu \geq 0, one constructs auxiliary functions (e.g., v = u + \epsilon |x|^2) to apply a contradiction argument: assuming an interior maximum leads to Lu(v) > 0 at that point, violating the inequality unless u is constant. The Hopf lemma follows from a local barrier construction using radial solutions near the boundary point.

Uniqueness and existence theorems

For linear second-order elliptic equations, uniqueness of solutions to the Dirichlet boundary value problem is established using the maximum principle and energy methods. Specifically, if Lu = 0 in a bounded domain \Omega \subset \mathbb{R}^n with u = 0 on \partial \Omega, where L is a uniformly elliptic operator with continuous coefficients, the maximum principle implies that u \equiv 0 in \Omega, as any non-trivial solution would attain a non-zero maximum or minimum interior to \Omega, contradicting the principle. For self-adjoint elliptic operators in divergence form (e.g., the Laplacian), energy methods provide an alternative proof: multiplying the equation by u and integrating by parts yields \int_\Omega |\nabla u|^2 \, dx \leq 0 (or =0 if c=0), implying u=0 under coercivity. For general operators, uniqueness relies primarily on the maximum principle.^[13] Another fundamental property is Harnack's inequality, which for positive solutions to Lu=0 states that \sup_B u / \inf_B u \leq C for balls B \subset \Omega, providing quantitative control on oscillations and aiding in regularity proofs.^[13] Existence of solutions for the Dirichlet problem for the Laplace equation \Delta u = 0 in \Omega with continuous boundary data u = \phi on \partial \Omega is guaranteed by Perron's method, which constructs the solution as the infimum of the class of superharmonic functions dominating \phi on the boundary. This method yields a harmonic function that attains the boundary values continuously under mild conditions on \Omega, such as the Wiener criterion for regularity. For more general linear elliptic boundary value problems, the Fredholm alternative applies when the associated operator is compact on appropriate Sobolev spaces, such as H^1(\Omega). In this framework, the problem Lu = f with homogeneous Dirichlet conditions has a unique solution if and only if f is orthogonal to the kernel of the adjoint operator; otherwise, solutions exist but are non-unique, with solvability determined by a finite-dimensional compatibility condition.^[14] In the nonlinear case, existence for semilinear elliptic equations like \Delta u + f(u) = 0 with Dirichlet data is often proved using Schauder estimates combined with the Schauder fixed-point theorem. These estimates provide Hölder continuity of solutions and their derivatives, allowing the nonlinear operator to be mapped into a compact subset of a Hölder space, where a fixed point yields the solution. A representative example is the Poisson equation \Delta u = g in \Omega with u = \phi on \partial \Omega, where g \in L^2(\Omega) and \phi \in H^{1/2}(\partial \Omega). Uniqueness follows from the maximum principle applied to v = u - u_0, where u_0 is a particular solution, and existence is obtained via the Lax-Milgram theorem in the weak formulation, ensuring a unique solution in H^1_0(\Omega).

Canonical forms

Transformation to standard form

For linear second-order partial differential equations with constant coefficients in two dimensions, the principal part is a quadratic form associated with the symmetric matrix corresponding to the coefficients of the second derivatives. Since the equation is elliptic, the discriminant is negative, ensuring both eigenvalues of this matrix are of the same sign and nonzero, allowing a linear change of variables—via rotation to align with eigenvectors followed by scaling—to diagonalize the form and reduce it to the Laplace operator plus lower-order terms:

\frac{\partial^2 u}{\partial \xi^2} + \frac{\partial^2 u}{\partial \eta^2} + \text{lower-order terms} = 0.

This canonical form simplifies analysis by aligning the equation with the well-understood Laplace equation. A concrete example illustrates this process for the equation u_{xx} + 5 u_{yy} = 0. Here, the matrix is diagonal with eigenvalues 1 and 5. Introduce new variables \xi = x and \eta = y / \sqrt{5}. Then, \partial^2 u / \partial y^2 = 5 \partial^2 u / \partial \eta^2, so the equation becomes u_{\xi\xi} + u_{\eta\eta} = 0, the standard Laplace equation. In the general linear case with variable coefficients, where the equation is uniformly elliptic, local diffeomorphisms can be employed to simplify the operator. Unlike hyperbolic or parabolic equations, elliptic equations lack real characteristics, so no "straightening" of characteristic curves is needed; instead, local coordinate changes near a point freeze the coefficients and reduce the principal part to the Laplacian, yielding a form \Delta u + \text{lower-order terms} = 0 in suitable coordinates. This local reduction relies on the ellipticity condition ensuring the principal symbol is invertible. However, variable coefficients generally prevent a global transformation to the standard form across the entire domain, as the required coordinate changes may not extend consistently without singularities or distortions.

Role of characteristics

In the theory of partial differential equations (PDEs), characteristics are defined as the curves or surfaces in the domain along which the principal part of the PDE governs the propagation of singularities or information. For a linear second-order PDE of the form \sum_{i,j=1}^n a_{ij}(x) \partial_i \partial_j u + \ lower\ order\ terms = 0, the principal symbol is the quadratic form p(x, \xi) = \sum_{i,j=1}^n a_{ij}(x) \xi_i \xi_j, where \xi \in \mathbb{R}^n \setminus \{0\}. The PDE is elliptic if this symbol never vanishes for \xi \neq 0, meaning p(x, \xi) \neq 0 for all real \xi \neq 0, or equivalently, the quadratic form is definite (positive or negative). This condition ensures there are no real directions \xi where the principal symbol vanishes, hence no real characteristics exist. In contrast, hyperbolic PDEs possess real characteristics, which are curves along which discontinuities or singularities in the solution and its derivatives propagate at finite speeds. For example, in the wave equation \partial_t^2 u - \Delta u = 0, the characteristics are light cones, and singularities travel along these cones without diffusing. The absence of real characteristics in elliptic PDEs implies that information does not propagate along specific real paths; instead, the domain of dependence for any interior point is the entire spatial domain, making the solution at that point dependent on boundary data everywhere. This global dependence leads to well-posed boundary value problems where smooth boundary data typically yield smooth solutions throughout the domain, without localized propagation of irregularities. The characteristics in the elliptic case are inherently complex, arising from the complex roots of the symbol equation. These complex characteristics play a role in analytic continuation properties, allowing solutions to elliptic PDEs with analytic coefficients to be analytically continued across the domain in suitable complex directions, which underpins regularity results. This differs from the real canonical transformations used to reduce elliptic equations to standard forms like the Laplace equation, but highlights the fundamentally nonlocal nature of elliptic problems.

Boundary value problems

Dirichlet and Neumann problems

The classical Dirichlet boundary value problem for a linear second-order elliptic partial differential equation Lu = f in a bounded domain \Omega \subset \mathbb{R}^n with smooth boundary \partial \Omega requires prescribing the solution values u = g on \partial \Omega, where L is an elliptic operator of the form L u = a_{ij} \partial_i \partial_j u + b_i \partial_i u + c u with uniform ellipticity.^[15] For sufficiently smooth data f \in C^\infty(\Omega) and g \in C^\infty(\partial \Omega), and smooth bounded domains, the problem is well-posed: a unique classical solution u \in C^2(\Omega) \cap C^0(\bar{\Omega}) exists, as established by existence and uniqueness theorems relying on potential theory and integral representations.^[16] This well-posedness extends the classical results for the Laplace equation \Delta u = 0 with u = g on \partial \Omega, where the solution is unique by the maximum principle.^[15] The Neumann boundary value problem specifies the normal derivative \frac{\partial u}{\partial \nu} = h on \partial \Omega for Lu = f in \Omega, where \nu is the outward unit normal.^[17] For smooth bounded domains and regular data, solvability requires a compatibility condition derived from the divergence theorem: \int_\Omega f \, [dx](/page/DX) = \int_{\partial \Omega} h \, d\sigma, ensuring consistency with the elliptic operator's structure (for the prototype -\Delta u = f, this is -\int_\Omega f = \int_{\partial \Omega} h).^[15] Under this condition, solutions exist but are unique only up to an additive constant, with the problem well-posed in appropriate function spaces like H^1(\Omega) for weak formulations.^[16] A representative example of the Dirichlet problem is the Poisson equation \Delta u = f in the unit disk \Omega = \{ (x,y) \in \mathbb{R}^2 : x^2 + y^2 < 1 \} with u = g on the unit circle, solved explicitly using the Poisson kernel:

u(r, \theta) = \frac{1}{2\pi} \int_0^{2\pi} \frac{1 - r^2}{1 - 2r \cos(\theta - \phi) + r^2} g(\phi) \, d\phi + \int_\Omega P(r, \theta; y) f(y) \, dy,

where P is the appropriate kernel extension, yielding a harmonic function adjusted for the source term.^[17] For the Neumann problem, consider steady-state incompressible fluid flow modeled by \Delta \phi = 0 in \Omega, with \frac{\partial \phi}{\partial \nu} = h on \partial \Omega specifying the normal velocity component; compatibility \int_{\partial \Omega} h \, d\sigma = 0 ensures conservation of flux, and the velocity potential \phi is determined up to a constant.^[17] Solutions to these problems can be represented using Green's functions G(x,y), which satisfy L G = \delta_y in \Omega with appropriate boundary conditions. For the Dirichlet problem, G = 0 on \partial \Omega, yielding

u(x) = \int_\Omega G(x,y) f(y) \, dy + \int_{\partial \Omega} \frac{\partial G}{\partial \nu_y}(x,y) g(y) \, d\sigma(y).

For the Neumann problem, \frac{\partial G}{\partial \nu} = \frac{1}{|\partial \Omega|} on \partial \Omega to fix the constant, giving

u(x) = \int_\Omega G(x,y) f(y) \, dy + \int_{\partial \Omega} G(x,y) h(y) \, d\sigma(y) + C,

where C is chosen arbitrarily; the fundamental solution in \mathbb{R}^n (e.g., \frac{1}{2\pi} \log |x-y| for n=2) is modified by a regular part to enforce boundary conditions.^[15]

Mixed boundary conditions

Mixed boundary conditions for elliptic partial differential equations combine Dirichlet conditions, where the solution u is prescribed as u = g on a portion \Gamma_D of the boundary \partial \Omega, and Neumann conditions, where the normal derivative \partial u / \partial n = h on the complementary portion \Gamma_N, such that \partial \Omega = \Gamma_D \cup \Gamma_N with \Gamma_D \cap \Gamma_N = \emptyset. This setup arises in modeling scenarios where different physical constraints apply to distinct boundary segments, such as insulated and prescribed-temperature surfaces in heat conduction. Well-posedness of these mixed boundary value problems hinges on geometric conditions at the interfaces between \Gamma_D and \Gamma_N, particularly the angles formed there. Singularities can develop at these junction points unless the interior angle satisfies specific constraints, such as being less than or equal to \pi radians for the Laplace equation, to ensure finite energy solutions and avoid unbounded gradients. Pierre Grisvard's analysis establishes that under such angle conditions, the problem admits a unique solution in appropriate Sobolev spaces, with regularity determined by the domain's smoothness away from the interfaces. The variational formulation seeks weak solutions in the Sobolev space H^1(\Omega) with zero trace on \Gamma_D, defined by integrating the elliptic operator against test functions vanishing on \Gamma_D, incorporating the Neumann data via the boundary integral. The Lax-Milgram theorem guarantees existence and uniqueness in this setting when the bilinear form is coercive and continuous, provided \Gamma_D has positive measure to ensure stability.

Advanced topics

Higher-order elliptic equations

Higher-order elliptic partial differential equations generalize the second-order case to operators of order m \geq 2, where the principal symbol determines the elliptic nature. A linear partial differential operator P of order m in n variables is locally expressed as Pu = \sum_{|\alpha| \leq m} a_\alpha(x) D^\alpha u, with the principal symbol given by the homogeneous polynomial p(x, \xi) = \sum_{|\alpha| = m} a_\alpha(x) \xi^\alpha, where \xi \in \mathbb{R}^n and D^\alpha denotes partial derivatives.^[18] The operator is elliptic if p(x, \xi) \neq 0 for all x in the domain and all \xi \neq 0, ensuring the symbol is invertible away from the origin and implying no real characteristics.^[18] A canonical example is the biharmonic equation \Delta^2 u = f, where \Delta is the Laplacian, which is elliptic of order 4 with principal symbol |\xi|^4. This equation arises in the modeling of thin elastic plates under transverse loading, where solutions represent the deflection of the plate.^[19] For smooth coefficients, elliptic operators of any order satisfy hypoellipticity: if Pu is smooth (or C^\infty) in an open set, then u is also smooth there, extending the regularity properties beyond the order of the operator.^[18] The theory of higher-order elliptic equations relies on parametrix constructions to analyze invertibility and regularity. For an elliptic operator P of order m on a compact manifold, a parametrix is a pseudodifferential operator Q of order -m such that PQ = I - R and QP = I - S, where I is the identity and R, S are smoothing operators (vanishing to infinite order). This approximate inverse facilitates proofs of Fredholm properties and elliptic estimates.^[18]

Nonlinear elliptic equations

Nonlinear elliptic partial differential equations (PDEs) generalize linear elliptic PDEs by incorporating nonlinearity in the highest-order derivatives, leading to more complex analytical behavior while preserving the elliptic structure through uniform ellipticity conditions. Quasilinear elliptic PDEs are nonlinear in the lower-order derivatives but linear in the second-order terms, typically of the form a^{ij}(x, u, Du) u_{ij} + b(x, u, Du) = 0, where the matrix (a^{ij}) is uniformly elliptic. A canonical example is the minimal surface equation, which describes surfaces of least area and takes the form

\operatorname{div}\left( \frac{\nabla u}{\sqrt{1 + |\nabla u|^2}} \right) = 0

in a domain \Omega \subset \mathbb{R}^n. This equation arises in the calculus of variations for minimizing the area functional \int_\Omega \sqrt{1 + |\nabla u|^2} \, dx subject to boundary conditions, and its solutions exhibit graph representations of minimal hypersurfaces. Fully nonlinear elliptic PDEs involve nonlinearity directly in the second derivatives, often expressed as F(D^2 u, Du, u, x) = 0, where F satisfies uniform ellipticity, meaning there exist constants $0 < \lambda \leq \Lambda < \infty such that \lambda I \leq D^2_{p p} F(M + P) \leq \Lambda I for symmetric matrices M, P with M \geq 0. A prominent example is the Monge-Ampère equation \det(D^2 u) = 1 in \Omega, which plays a central role in optimal transport theory, where it characterizes the optimal mapping between probability measures via the Brenier potential. Existence of smooth convex solutions to the Dirichlet problem for this equation was established under strict convexity of the domain, with interior C^{2,\alpha} regularity following from the convexity of the solution. Challenges in nonlinear elliptic PDEs include potential non-uniqueness of solutions without additional structural assumptions, such as convexity of the admissible set or the operator. For the Monge-Ampère equation, non-convex solutions may fail to exist or be unique, as the equation degenerates outside the convex regime, complicating boundary value problems. In control theory, Bellman equations of the form \sup_{a \in A} \left( - \operatorname{tr}(A^a(x) D^2 u) + H(x, Du, a) \right) = 0 arise as fully nonlinear elliptic PDEs for value functions in stochastic optimal control, where viscosity solution theory ensures uniqueness under monotonicity and continuity of the Hamiltonian, but multiple subsolutions may exist without proper discounting. The De Giorgi-Nash theory establishes Hölder continuity of weak solutions to quasilinear elliptic equations with bounded measurable coefficients, extending the interior regularity results for linear cases to structures like the p-Laplace equation \operatorname{div}(|\nabla u|^{p-2} \nabla u) = 0 for $1 < p < \infty. This theory relies on iterative estimates from integral inequalities, yielding C^{0,\alpha} bounds independent of the solution's magnitude, which is crucial for bootstrapping higher regularity in nonlinear settings.

Regularity theory

Interior estimates

Interior estimates provide a priori bounds on the regularity of solutions to elliptic partial differential equations within compact subsets of the domain, independent of boundary behavior. These estimates are crucial for establishing higher-order smoothness and Hölder continuity for solutions of linear uniformly elliptic equations of the form Lu = f, where L is a second-order operator with bounded measurable coefficients and f belongs to appropriate Hölder or Lebesgue spaces. Such bounds rely on the ellipticity condition and enable bootstrapping arguments to achieve classical regularity.^[20] A foundational result in this context is the De Giorgi-Nash-Moser theorem, which establishes Hölder continuity for weak solutions of uniformly elliptic equations with measurable coefficients. For the divergence-form equation -\operatorname{div}(A(x) \nabla u) = 0 with A uniformly elliptic (bounded measurable entries), solutions u \in W^{1,2}_{\mathrm{loc}}(\Omega) satisfy u \in C^\alpha_{\mathrm{loc}}(\Omega) for some \alpha > 0 depending on dimension n and ellipticity constants, with explicit estimates like

\osc_{B_r(x_0)} u \leq C \left( \frac{r}{R} \right)^\alpha \osc_{B_R(x_0)} u

for concentric balls B_r \subset B_R \subset \subset \Omega. The Nash-Moser iteration provides Harnack inequalities for positive solutions, while De Giorgi's geometric approach yields the full Hölder result; this enables bootstrapping to higher regularity when coefficients are smoother.^[20] Schauder estimates constitute a cornerstone of interior regularity theory, quantifying the Hölder continuity of second derivatives in terms of the data. For a solution u to Lu = f in a domain \Omega \subset \mathbb{R}^n, where L = a_{ij}(x) \partial_{ij} + b_i(x) \partial_i + c(x) satisfies uniform ellipticity with constants \lambda, \Lambda > 0, the interior Schauder estimate asserts that for any compact set K \subset \Omega and $0 < \alpha < 1,

\|u\|_{C^{2,\alpha}(K)} \leq C \left( \|f\|_{C^{\alpha}(K)} + \|u\|_{L^\infty(\Omega)} \right),

where C depends on n, \alpha, \lambda, \Lambda, and the bounds on the coefficients. This result, originally due to Schauder, extends to more general elliptic operators and underpins global regularity via covering arguments.^[20] In the L^p framework, Calderón-Zygmund estimates address integrability properties through the decomposition of solutions into singular integral operators. For divergence-form elliptic equations -\operatorname{div}(A(x) \nabla u) = \operatorname{div} \mathbf{F} with A uniformly elliptic, these estimates yield

\|\nabla^2 u\|_{L^p(B_r)} \leq C \left( \|\mathbf{F}\|_{L^p(B_{2r})} + \|u\|_{L^p(B_{2r})} / r \right)

for balls B_r interior to the domain and $1 < p < \infty, with C independent of the specific ball. This theory, building on singular integral decompositions, applies to nondivergence forms via perturbation and provides W^{2,p} bounds essential for Sobolev regularity.^[20]^[21] Moser's Harnack inequality offers a pointwise control for positive solutions, linking values across subdomains. For nonnegative solutions u to Lu = 0 in B_{2r} \subset \Omega, it states

\sup_{B_r} u \leq C \inf_{B_r} u,

where C depends only on n, \lambda, \Lambda, and the coefficient bounds; this holds under the strong maximum principle, which ensures nonconstant positive solutions cannot attain interior minima. The inequality facilitates oscillation decay and Hölder continuity via iteration.^[20]^[22] Proofs of these estimates often invoke potential theory to represent solutions via Newtonian or Riesz potentials, followed by Campanato-Morrey characterizations of Hölder spaces. For Schauder estimates, the second derivatives are estimated by differentiating the potential integral and applying Schauder integrability conditions on the kernel. Iteration techniques, such as Moser's method of testing with subsolutions and supersolutions raised to powers, bound oscillations by applying the maximum principle repeatedly to auxiliary functions, yielding exponential convergence to the mean value.^[20]^[23]

Boundary regularity results

Boundary regularity results extend the interior Schauder estimates to regions near the boundary of the domain, providing Hölder continuity and higher-order derivatives for solutions of linear elliptic equations up to the boundary, assuming suitable regularity of the boundary data and the domain itself. These estimates are crucial for establishing global regularity in bounded domains and solving boundary value problems classically. For flat boundaries, such as in the half-space \mathbb{R}^n_+ = \{x = (x', x_n) \in \mathbb{R}^n : x_n > 0\}, Schauder boundary estimates yield bounds on the C^{2,\alpha} norm of the solution u near the boundary in terms of the C^\alpha norm of the right-hand side f and the boundary data. Specifically, for a solution u to Lu = f in B_1^+ \cap \mathbb{R}^n_+ with L a uniformly elliptic operator with C^\alpha coefficients, the estimate takes the form

\|u\|_{C^{2,\alpha}(\overline{B_r^+})} \leq C \left( \|f\|_{C^\alpha(B_1^+)} + \|\phi\|_{C^\alpha(B_1 \cap \partial \mathbb{R}^n_+)} + \|u\|_{L^\infty(B_1^+)} \right),

where \phi is the Dirichlet boundary data on the flat boundary \{x_n = 0\} \cap B_1, B_r^+ is a ball of radius r intersecting the boundary, and C depends on n, \alpha, and ellipticity constants. This result, derived using reflection principles and potential theory, holds for Dirichlet boundary conditions and extends to other linear conditions under appropriate assumptions.^[20] In curved domains \Omega with C^1 boundary, global boundary regularity is obtained by locally flattening the boundary via C^1 diffeomorphic charts that map portions of \partial \Omega to flat hyperplanes, preserving the ellipticity of the operator up to lower-order terms. Under the assumption that \partial \Omega \cap B_r is C^1 and the coefficients and data are C^\alpha, the Schauder estimates localize to yield

\|u\|_{C^{2,\alpha}(\overline{\Omega} \cap B_r)} \leq C \left( \|Lu\|_{C^\alpha(\Omega \cap B_1)} + \|\text{boundary data}\|_{C^\alpha(\partial \Omega \cap B_1)} + \|u\|_{L^\infty(\Omega \cap B_1)} \right),

with C independent of the local geometry beyond the C^1 assumption. This flattening technique ensures the estimates hold uniformly near the boundary, enabling global C^{2,\alpha} regularity when combined with interior estimates.^[20] For oblique derivative problems, where the boundary condition is \beta \cdot Du = g on \partial \Omega with \beta a non-tangential vector field, boundary estimates require the oblique direction to satisfy a non-degeneracy condition, such as |\beta \cdot \nu| \geq \theta > 0 for the outward normal \nu. In this setting, Agmon-Douglis-Nirenberg established C^{2,\alpha} estimates near the boundary for solutions in domains with C^{1,1} boundary, generalizing the flat case and yielding

\|u\|_{C^{2,\alpha}(\overline{\Omega} \cap B_r)} \leq C \left( \|Lu\|_{C^\alpha} + \|g\|_{C^{1,\alpha}} + \|u\|_{C^\alpha} \right),

provided the operator and boundary data satisfy the necessary Hölder conditions. These results apply to general linear boundary conditions, including mixed types, as long as the conormal or oblique directions avoid tangential degeneracy.^[20] A representative example is the Dirichlet problem for the Laplace equation -\Delta u = 0 in a C^{1,\alpha} domain [\Omega](/page/Omega) with u = [\phi](/page/Phi) on [\partial \Omega](/page/boundary), where [\phi](/page/Phi) \in C^\alpha([\partial \Omega](/page/boundary)). The solution u satisfies Hölder continuity up to the boundary, with \|u\|_{C^\alpha(\overline{[\Omega](/page/Omega)})} \leq C \|[\phi](/page/Phi)\|_{C^\alpha([\partial \Omega](/page/boundary))}, ensuring classical solvability and paving the way for higher regularity when [\phi](/page/Phi) \in C^{k,\alpha}. This follows directly from the Schauder boundary estimates applied iteratively.^[20]

Applications

Physical models

Elliptic partial differential equations commonly arise in physical models describing steady-state phenomena and equilibrium configurations, where time derivatives vanish and the systems reach a balance without transient effects. In electrostatics, the electric potential \phi in a region with charge density \rho satisfies Poisson's equation \nabla^2 \phi = -\rho / \epsilon_0, where \epsilon_0 is the permittivity of free space; this elliptic equation determines the potential from which the electric field \mathbf{E} = -\nabla \phi is derived, modeling the equilibrium distribution of electrostatic forces.^[24] In the absence of charges (\rho = 0), the equation reduces to Laplace's equation \nabla^2 \phi = 0, applicable to charge-free regions like conductors or insulators in equilibrium.^[24] Steady-state heat flow in regions without internal heat sources is governed by Laplace's equation \nabla^2 u = 0 for the temperature field u, representing thermal equilibrium in insulated materials where heat flux balances across boundaries; this models constant-temperature distributions in homogeneous media under fixed boundary conditions.^[25] With distributed sources, such as uniform heating, the equation becomes Poisson's form \nabla^2 u = -f / k, where f is the source term and k the thermal conductivity, capturing balanced heat generation and conduction.^[25] In fluid dynamics, the Stokes equations describe slow, viscous incompressible flows at low Reynolds numbers, where inertial effects are negligible; the system consists of -\Delta \mathbf{u} + \nabla p = \mathbf{f} for the velocity \mathbf{u} and pressure p, coupled with the incompressibility condition \nabla \cdot \mathbf{u} = 0, forming a linear elliptic PDE system that models creeping flows around obstacles, such as sediment particles or biological swimmers in viscous media.^[26] Newtonian gravitational theory employs Poisson's equation \nabla^2 \Phi = 4\pi G \rho for the gravitational potential \Phi, with mass density \rho and gravitational constant G; the acceleration \mathbf{g} = -\nabla \Phi follows, representing the equilibrium potential due to a static mass distribution in the weak-field, non-relativistic limit.^[27] For vacuum regions (\rho = 0), it simplifies to Laplace's equation, applicable to gravitational fields outside massive bodies.^[27]

Geometric and analytical uses

Elliptic partial differential equations play a central role in geometric analysis, particularly in the study of flows that evolve submanifolds to minimize area. Stationary points of the mean curvature flow, which describes the motion of hypersurfaces by their mean curvature vector, are minimal hypersurfaces satisfying the quasilinear elliptic minimal surface equation \operatorname{div} \left( \frac{\nabla u}{\sqrt{1 + |\nabla u|^2}} \right) = 0 for graphs.^[28] These solutions represent critical points of the area functional and exhibit regularity properties derived from the ellipticity of the operator.^[29] In higher dimensions, such stationary configurations inform the understanding of singularity formation and long-time behavior in the flow.^[30] In complex analysis, elliptic PDEs underpin conformal mapping results through the Dirichlet problem for harmonic functions. The Riemann mapping theorem asserts that any simply connected domain in the complex plane, distinct from the whole plane, is conformally equivalent to the unit disk, with the mapping obtained as the real part of a holomorphic function whose imaginary part solves a Dirichlet boundary value problem for the Laplace equation.^[31] This construction relies on the solvability and uniqueness of the elliptic Dirichlet problem, ensuring the harmonic extension matches prescribed boundary data continuously.^[32] The theorem's proof via this method highlights the interplay between elliptic regularity and geometric uniformity in the plane. Variational problems in analysis often lead to elliptic PDEs as Euler-Lagrange equations for energy functionals. For the Dirichlet energy functional \int_\Omega |\nabla u|^2 \, dx, minimizers satisfy the linear elliptic equation -\Delta u = 0, corresponding to harmonic functions that achieve balance between boundary conditions and smoothness.^[33] This principle extends to more general functionals, where the resulting second-order equations inherit ellipticity from the positive definiteness of the quadratic form in the first variation, enabling existence via direct methods in the calculus of variations. Spectral theory of elliptic operators provides tools for analyzing manifold geometry through eigenvalues and eigenfunctions. On a compact Riemannian manifold, the Laplace-Beltrami operator -\Delta_g is a self-adjoint elliptic operator with discrete spectrum \{ \lambda_k \}_{k=1}^\infty of non-negative eigenvalues accumulating at infinity, where \lambda_1 = [0](/page/0) corresponds to constant functions and higher eigenvalues encode geometric invariants like diameter and volume. Estimates on eigenvalue sums and gaps, such as Weyl's law N(\lambda) \sim c \operatorname{Vol}(M) \lambda^{n/2} for dimension n, relate spectral data to asymptotic manifold properties.^[34] These results facilitate applications in shape optimization and inverse problems on manifolds.

References

[1]
[PDF] Chapter 4: Elliptic PDEs - UC Davis Math
There is often considerable freedom in how one defines a weak solution of a. PDE; for example, the function space to which a solution is required to belong is.
[2]
None
Summary of each segment:
[3]
Elliptic Equation - an overview | ScienceDirect Topics
Second-order partial differential equations can be classified into elliptic, parabolic, and hyperbolic forms. Elliptic equations are encountered in steady state ...
[4]
[PDF] Notes on Partial Differential Equations John K. Hunter - UC Davis Math
These are notes from a two-quarter class on PDEs that are heavily based on the book Partial Differential Equations by L. C. Evans, together with other sources ...
[5]
1.3 Classification
The classification for these equations is: $b^2 - ac > 0$ : hyperbolic; $b^2 - ac = 0$ : parabolic; $b^2 - ac < 0$ : elliptic. Example. Question: Classify the ...
[6]
[PDF] Partial Differential Equations
Jan 6, 2024 · ... Classification of equations . . . . . . . . . . . . . . . . . . . 9 ... (d'Alembertian or d'Alembert operator). It appears in elasticity ...
[7]
[PDF] s heat conduction equation: History, influence, and connections
Fourier's heat conduction equation, formulated in the 19th century, describes heat conduction in solids and is used for diffusion problems in many fields.Missing: classification | Show results with:classification
[8]
[PDF] Carl Friedrich Gauss – General Theory of Terrestrial Magnetism
Feb 5, 2014 · Gauss' Theory also helped to resolve several issues of discussion. For ex- ample, the introduction of a magnetic potential and its spher- ical ...
[9]
An essay on the application of mathematical analysis to the theories ...
Jun 21, 2022 · An essay on the application of mathematical analysis to the theories of electricity and magnetism ; Publication date: 1828 ; Topics: Potential ...
[10]
[PDF] Mathematical Problems
Mathematical Problems. Lecture delivered before the International Congress of. Mathematicians at Paris in 1900. By Professor David Hilbert ∗. Wednesday, August ...
[11]
On interior regularity of the solutions of partial differential equations
On interior regularity of the solutions of partial differential equations. Lars Hörmander, ... elliptic differential equations, Comm. Pure Appl. Math., Vol ...
[12]
[PDF] Partial Differential Equations - UC Berkeley math
This article is an extremely rapid survey of the modern theory of partial differential equations. (PDEs). Sources of PDEs are legion: mathemat-.
[13]
[PDF] Boundary value problems for elliptic partial differential equations - HAL
... Fredholm's alternative theorem is useful tool to solve the equation u−Au = f. (2.7). Theorem 2.5 says that we have the following alternative. • For any f ...
[14]
Methods of Mathematical Physics: Partial Differential Equations
Apr 19, 1989 · Courant and Hilbert's treatment restores the historically deep connections between physical intuition and mathematical development ...
[15]
[PDF] notes on elliptic boundary value problems for the laplace operator
The classical problems of this type are the Dirichlet and Neumann problems for the Laplace equation on a smoothly bounded domain, Ω ⊆ Rn : Dirichlet Problem.
[16]
[PDF] Boundary Value Problems
We begin by considering an abstract “weak” problem (2.1) motivated by certain carefully chosen formulations of the Dirichlet and Neumann problems for the ...
[17]
[PDF] Chapter 6 Partial Differential Equations
PARTIAL DIFFERENTIAL EQUATIONS. A hyperbolic second-order differential equation Du = 0 can therefore be written in either of two ways: α ∂∂x + β ∂∂y U1 + F1 = 0 ...
[18]
Regularity of a transmission problem and periodic homogenization
This paper is concerned with the regularity theory of a transmission problem arising in composite materials.
[19]
[PDF] Elliptic Partial Differential Operators
Aug 10, 2023 · We define the symbol of a partial differential operator, we define what it means for an operator to be elliptic, and we prove the existence ...
[20]
[PDF] Fourth-order compact finite difference schemes for solving ... - arXiv
Sep 4, 2024 · The biharmonic equation, a fourth order partial differential equation(PDE), plays a pivotal role in various branches of applied mathematics ...
[21]
[PDF] Elliptic Partial Differential Equations of Second Order
David Gilbarg • Neil S.Trudinger. Elliptic Partial. Differential Equations of Second Order. Reprint of the 1998 Edition. Springer. Page 2. Table of Contents.Missing: book excerpt
[22]
A Geometric Approach to the Calderón-Zygmund Estimates
In this note, we will show a new proof of the classical Calderón-Zygmund estimates established in [7] 1952. The estimates are among the fundamental estimates ...
[23]
On Harnack's theorem for elliptic differential equations
On Harnack's theorem for elliptic differential equations ; First published: August 1961 ; Citations · 570 ; This paper represents results obtained under Contract ...
[24]
[PDF] PDE II – Schauder estimates
Schauder estimates for PDEs, like the interior estimate, involve a constant C. The fundamental estimate for the Laplacian on Rn is |D2u|α ≤ C|4u|α.Missing: seminal paper
[25]
5.15: Poisson's and Laplace's Equations - Engineering LibreTexts
Sep 12, 2022 · Note that Poisson's Equation is a partial differential equation, and therefore can be solved using well-known techniques already established ...
[26]
4.9: Steady State Temperature and the Laplacian - Math LibreTexts
Feb 23, 2025 · Solutions to the Laplace equation are called harmonic functions and have many nice properties and applications far beyond the steady state heat ...
[27]
[PDF] 7 The Navier-Stokes Equations - MIT Mathematics
For low Reynolds number it may be possible to ignore the inertial terms in the Navier-Stokes equations and obtain the so-called slow (or creeping) flow ...Missing: textbook reference
[28]
[PDF] 1 Foundations of Newtonian gravity - UF Physics
Dec 16, 2013 · ... gravitational field is Poisson's equation ... Foundations of Newtonian gravity one of the key defining equations for the Newtonian potential.
[29]
[PDF] arXiv:math/0411589v1 [math.AP] 26 Nov 2004
He shows that the mean curvature flow (the minus gradient flow. 5. Page 6. Luca Martinazzi of the area functional) of the initial graph Gψ (now ψ is intended ...
[30]
The problem of dirichlet for quasilinear elliptic differential equations ...
This paper is concerned with the existence of solutions of the Dirichlet problem for quasilinear elliptic partial differential equations of second order.
[31]
[PDF] research statement - theodora bourni - Mathematics
There is an obvious connection between minimal surfaces and mean curvature flow which is the fact that the former are stationary solutions of the latter. A ...
[32]
[PDF] arXiv:1604.04071v2 [math.CV] 13 Dec 2016
Dec 13, 2016 · The Riemann Mapping Theorem is one of the most remarkable re- sults of nineteenth century mathematics. Even today, more than a.Missing: scholarly | Show results with:scholarly<|separator|>
[33]
The Riemann mapping theorem from Riemann's viewpoint
Jan 16, 2017 · This article presents a rigorous proof of the Riemann mapping theorem via Riemann's method, uncompromised by any appeals to topological intuition.
[34]
Dirichlet Energy Integral and Laplace Equation - SpringerLink
The Euler–Lagrange equation of this problem is the q-Laplace equation,div(\mid\bigtriangledown u \mid^{q-2}\bigtriangledown u )\, = \,0 , which has to be ...
[35]
Spectral Properties of the Laplace-Beltrami Operator and Applications
... eigenvalues of the Laplace-Beltrami operator associated to a closed Riemannian manifold (M ... eigenvalues, in “Spectral Theory and Geometry” (edited by B.
[36]
On sums of eigenvalues of elliptic operators on manifolds - arXiv
Jul 9, 2015 · Among the operators we consider are the Laplace-Beltrami operator on compact subdomains of manifolds. These estimates become more explicit ...