On the torus , choose normalized Haar measure, so that , , is an orthonormal basis of . A periodic Sobolev space is defined, for any real , by
For negative , the coefficients describe distributions rather than necessarily functions. For a nonnegative integer , the weight is comparable to , so this norm is equivalent to the sum of the squared norms of all weak derivatives through order . Fourier truncation proves that smooth trigonometric polynomials are dense. The pairing identifies the continuous dual of with , giving Sobolev duality.
Two Fourier estimates explain the usefulness of these spaces. If , the Cauchy-Schwarz inequality gives
Uniform convergence of the differentiated Fourier series proves the Sobolev embedding . If , a bounded subset of has uniformly small high-frequency tails in :
Finite-dimensional Fourier truncation therefore proves the compact embedding , the periodic version of Rellich-Kondrachov compactness theorem. In addition, for multiplication is continuous , making it a Sobolev algebra. Indeed the weighted coefficient of a product is bounded by the sum of the convolutions of with and of with ; Cauchy-Schwarz puts unweighted coefficients in , and the case of Young's convolution inequality finishes the estimate.
For eigenfunction expansions, distinguish ellipticity from self-adjointness. An arbitrary non-self-adjoint elliptic operator need not have an orthonormal basis of eigenfunctions. A useful definite setting is
with smooth real periodic coefficients, , real , and uniform ellipticity
Take the self-adjoint realization in with domain . Its construction and spectral properties can be seen from the associated form on ,
For a sufficiently large constant , is a positive coercive form: , and . The Lax-Milgram theorem supplies the unique weak solution of for ; the form is Hermitian, which gives symmetry of its inverse on .
The regularity input is the periodic elliptic estimate
Here is why an elliptic operator gains two derivatives. With coefficients frozen at a point, its principal symbol dominates , and the estimate follows by multiplying Fourier coefficients by that symbol. On a sufficiently small coordinate patch, the difference of the leading coefficients from these constants has arbitrarily small supremum norm. In the integer-order differentiated estimate, this absorbs the highest derivatives into the left side; derivatives falling on coefficients contribute only lower-order terms. A partition of unity reduces the variable-coefficient problem to such patches. Its commutators with have order one, controlled by
which follows by splitting low and high Fourier modes. Absorbing these terms gives the displayed estimate. For a weak solution, the same localization with difference quotients first bounds its second weak derivatives for ; differentiating and repeating then gives for nonnegative integer . Interpolation and duality extend the scale to real orders. In particular is bounded.
Let . Compactness of makes a compact operator. Symmetry of the form makes self-adjoint, and positivity of makes positive and injective. Its range is dense because . The spectral theorem for compact self-adjoint operators consequently gives an orthonormal basis with
Elliptic regularity bootstraps every to a smooth function. Each eigenvalue has finite multiplicity. Thus a self-adjoint uniformly elliptic second-order operator on the torus has a complete smooth orthonormal eigenfunction expansion with discrete eigenvalues tending to infinity.
Writing , every has in , with . Iterating the elliptic estimate shows that the graph norm of is equivalent to the norm. Interpolation and duality give the more general spectral description
The comparison constants depend on the operator and . In particular smooth functions have coefficients decaying faster than every power of , and their expansions converge in every Sobolev norm and hence in every norm.
This gives a practical way to solve elliptic equations. The equation is solvable if and only if whenever ; then one may take on the other eigenspaces and add an arbitrary vector in . The Sobolev description proves the two-derivative gain and the Fredholm alternative, rather than leaving them as formal series manipulations. Similarly the solution of the heat equation is
For every the exponential decay at high eigenvalues places this solution in every Sobolev space, hence it is smooth.
For the Laplacian the eigenfunctions are precisely with eigenvalues . The number with is the number of lattice points in a radius- ball, asymptotic to its volume by comparison with unit cubes. For the general self-adjoint operator above, bounds on its form above and below by positive multiples of plus constants, together with the Courant–Fischer min-max principle, imply . Thus the growth reflects dimension and the two-derivative order. Compactness of the torus is essential: on , translated bumps prevent the compact embedding, and the Laplacian instead has continuous spectrum described by the Fourier transform.
The sobolev embedding theorem in three dimensions and the periodic elliptic estimate for the Stokes operator give
The last step uses the absence of the zero Fourier mode: on mean-zero periodic fields, the Poincaré inequality makes the homogeneous norm controlled by .
Take the inner product of the temperature equation with . The spectral projection disappears against , and part (a)(iii) cancels transport. Thus
The periodic Poincaré inequality, part (i), and the Cauchy-Schwarz inequality imply
The Gronwall inequality therefore gives, for ,
This defines a bound independent of , and part (i) then gives
Integrating the energy identity and using the same bound on its right-hand side gives
It remains to estimate the time derivative. For , the Fourier projection is a contraction in , and the skew identity from part (a) gives
The sobolev embedding theorem and the periodic elliptic estimate for the Stokes operator bound by . Moreover,
The already obtained bounds therefore imply
with independent of .