advanced30 min studyLesson 391 in Mathematics

Uniform ellipticity and principal symbols

Read the highest-order part of an operator and distinguish pointwise positivity from a uniform directional bound.

On this page

What you’ll understand

  • Compute a second-order scalar principal symbol with a stated sign convention.
  • Find ellipticity constants from a symmetric coefficient matrix.
  • Explain why pointwise positive coefficients need not give a uniform lower bound.

Elliptic equations describe spatial balance: a source acts throughout a region while boundary data constrain the response. Their highest-order coefficients decide whether every spatial direction receives genuine control. A matrix that looks positive at each sampled location can still lose that control near an unsampled boundary or interface. Uniform ellipticity records one positive lower bound that works everywhere in the chosen domain and in every direction. We begin with real scalar second-order operators; systems and higher-order equations need more specialized definitions.

Freeze the top-order coefficients

For a divergence-form operator Lu=−∂i(aij(x)∂ju)Lu=-\partial_i(a_{ij}(x)\partial_j u), the leading derivatives are described by the coefficient matrix A. When the coefficients are differentiable, expanding the product also produces first-order terms; these do not enter the second-order principal symbol. With the convention that differentiating a Fourier mode contributes i times its frequency, the two derivatives and the leading minus sign give the real quadratic symbol xi-transpose-A-xi. Stating this convention avoids confusing a positive elliptic operator with a negative Laplacian symbol.

The nondivergence operator −aij(x)∂i∂ju-a_{ij}(x)\partial_i\partial_j u has the same principal symbol but is not generally the same operator when A varies in space. Divergence form has a natural weak formulation requiring only bounded measurable coefficients and first weak derivatives of u. Nondivergence equations require a different interpretation at low regularity. Agreement of the principal symbols therefore identifies the leading directional behavior without erasing differences in lower-order terms, solution spaces, or theorems.

p2(x,ξ)=∑i,j=1naij(x)ξiξj,L=−∂i(aij∂j)p_2(x,\xi)=\sum_{i,j=1}^n a_{ij}(x)\xi_i\xi_j,\qquad L=-\partial_i(a_{ij}\partial_j)
Repeated indices are summed. This sign convention gives the positive symbol |ξ|² for −Δ. The symbol uses the symmetric part of a real scalar coefficient matrix.↔ Scroll wide equations horizontally.

One bound must survive all locations and directions

For a real symmetric measurable matrix field, assume positive constants lambda and Lambda bound the quadratic form below and above almost everywhere, for every real vector. This is the uniform ellipticity condition used throughout the module. It bounds the eigenvalues of A between those constants and makes the energy density comparable with the squared gradient. The ratio Lambda/lambda measures the permitted anisotropy, although estimates can also depend on coefficient regularity, domain geometry, and lower-order terms.

Pointwise positive definiteness only asks each location to have some positive smallest eigenvalue. Those location-dependent values may approach zero, leaving no common lambda. For example, the matrix diag(x squared,one) is positive definite for every x between zero and one, but its first-direction energy becomes arbitrarily small near zero. On a compact set, a continuous positive-definite matrix field has a uniform positive lower bound; on an open domain without control up to its boundary, that compactness argument is unavailable.

λ∣ξ∣2≤ξTA(x)ξ≤Λ∣ξ∣2,0<λ≤Λ<∞\lambda|\xi|^2\leq\xi^{\mathsf T}A(x)\xi\leq\Lambda|\xi|^2,\qquad 0<\lambda\leq\Lambda<\infty
The same constants work for almost every x and every ξ. Symmetric A and these bounds are assumed here; general nonsymmetric formulations separately control coefficient size and the coercive symmetric part.↔ Scroll wide equations horizontally.

Ellipticity is a starting hypothesis, not a full theorem

Uniform ellipticity does not by itself select a unique solution. A harmonic function can be changed by another harmonic function unless boundary conditions or a normalization remove that freedom. It also does not turn measurable coefficients into differentiable ones, so one cannot differentiate the equation and assume all resulting terms are classical functions. Existence, uniqueness, interior regularity, and boundary regularity each require an appropriate combination of data, domain, coefficient, and function-space assumptions.

The coefficient matrix can be interpreted as a directional conductivity in a simple diffusion model. That interpretation is useful for identifying which flux is continuous across an interface, but it does not justify treating every elliptic equation as a physical material model. In computations, test the smallest and largest eigenvalues across the represented coefficient family, and identify what happens between grid points. A finite list of positive eigenvalues is evidence about samples, whereas a uniform analytical bound controls the continuum model.

Read anisotropy without ignoring the cross term

For A=(3113)A=\begin{pmatrix}3&1\\1&3\end{pmatrix}, compute the symbol and tight ellipticity bounds. Evaluate it on ξ=(1,−1).

  1. The quadratic form is 3ξ12+2ξ1ξ2+3ξ223\xi_1^2+2\xi_1\xi_2+3\xi_2^2. The cross term has coefficient two because both off-diagonal entries contribute. The eigenvectors (1,1) and (1,−1) have eigenvalues four and two.

  2. Since A is symmetric, an orthonormal eigenbasis writes every directional energy as a weighted sum of squared components. Thus the tight constants are λ=2 and Λ=4. The matrix is uniformly elliptic on any domain because it is constant.

  3. At ξ=(1,−1), the symbol is 3−2+3=4 and the squared Euclidean norm is two. Their ratio is two, attaining the lower ellipticity bound. The most strongly controlled direction is instead parallel to (1,1).

The result

The cross term rotates the principal directions; looking only at the diagonal entries would miss the tight bounds two and four.

Explore the model

Interactive lab

Interactive lab · matrix

Turn geometric scaling into a positive energy

Pause, Reset, then set diagonal scale a=1 and shear b=0. Repeat from Reset with both a=0.25 and b=0. Predict singular values 0.25 and one in the second run. For the symmetric energy matrix G=TᵀT built from the displayed map T, those values square to ellipticity bounds 0.0625 and one. The triangular map itself is not being presented as a symmetric conductivity matrix or a PDE solver.

Turn geometric scaling into a positive energyThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA small singular direction produces a small quadratic-energy direction; a genuine uniform bound must control the whole allowed coefficient family.

Check your understanding

For A(x)=diag(x²,1) on 0<x<1, which statement is correct?

Try it yourself

For A(x)=diag(1+x²,2) on −1<x<1, find valid tight uniform constants. Why does the open endpoint not prevent an upper bound of two?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Bound each diagonal entry and remember that an inequality can use a supremum even when it is not attained.

Compare with the solution

The first eigenvalue lies between one and two, and the second is exactly two. Thus λ=1 and Λ=2 are valid and tight. The lower value one is attained at x=0; the upper value is attained by the second coordinate at every x. Even without that second coordinate, the supremum of 1+x² is two, and an upper inequality need not attain equality to remain valid.

Take it with you

  • A principal symbol isolates the leading directional action with an explicit sign convention.
  • Uniform ellipticity needs common bounds over both space and frequency directions.
  • Boundary conditions and coefficient regularity remain separate inputs to elliptic theory.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.