Choose an upper ellipticity bound for the symmetric coefficient matrix, possible because its entries are bounded. Then
The Poincaré inequality on a bounded open set for zero-boundary functions states . Consequently
The symmetric positive bilinear form is therefore an inner product, and its square-root norm is equivalent to the usual Sobolev norm. Its triangle inequality follows from the Cauchy-Schwarz inequality for that inner product. completeness of in its original norm makes it a Hilbert space in the norm as well.
The linear functional is bounded, since
The Hilbert-space Riesz representation theorem states that every bounded linear functional is the inner product with a unique element. Applied in the inner product it gives a unique with for all , which is precisely the weak-solution identity. Alternatively uniqueness follows directly by testing the difference of two solutions with itself: . Thus .
If and , linearity of the weak equation makes solve the equation with forcing . Uniqueness therefore gives .
Test the equation for with itself. The Poincaré inequality yields
If the gradient norm is nonzero, divide by it; if it is zero, Poincare gives and the estimate is automatic. In either case
Hence is a bounded linear operator. Notice that the energy identity has a minus sign, so . Testing the equations for and against each other gives . Thus the bounded Dirichlet inverse is a negative self-adjoint operator on .
For fixed real , consider initial data and a solution continuous into at time zero, with homogeneous Dirichlet boundary conditions for positive times. Existence, uniqueness and continuous dependence in this norm give Hadamard well-posedness.
To construct a solution, put . Substitution gives
Let . The Fourier sine series solution is
Multiplication by is bounded on . Completeness of the sine basis gives convergence to in at time zero, while the exponential decay gives smoothness for positive times and the required equation and boundary values. For a classical solution continuous up to the endpoints at time zero, the usual initial-boundary compatibility conditions are additionally required.
For the energy method, integration by parts gives
Both boundary terms vanish. The same identity applies to the difference of two solutions, proving uniqueness and continuous dependence. The Poincaré inequality further gives
with the identical estimate for differences of initial data. This Dirichlet convection-diffusion contraction holds for every real , and its energy constant is independent of the drift strength.
Let and , the closure of smooth compactly supported functions in the Sobolev space . Its elements have zero boundary trace. Set
The form is symmetric and continuous by Cauchy-Schwarz. For a smooth function with zero boundary values, gives after integration. Density extends this Poincaré inequality to . Therefore the gradient norm makes complete and equivalent to its usual norm, with
All the hypotheses of part (a) hold. Green's identity identifies the weak solution of the Dirichlet problem as
For smooth this is with zero boundary values; in general it defines the unique weak solution without requiring pointwise second derivatives.
For the Ritz method, choose a finite-dimensional conforming space with basis and put . Minimize in the coefficients. Its derivative in each basis direction gives the finite element method equations
The stiffness matrix is symmetric positive definite: for a nonzero coefficient vector , the corresponding basis combination is nonzero and . Hence the Ritz system has a unique solution.
For example, triangulate the square and use continuous piecewise-affine nodal basis functions at the interior vertices, with boundary vertex values fixed to zero. On a triangle of area , let be its barycentric basis functions. If its counterclockwise vertices are , define and cyclically. Then , giving the directly assemblable element matrices and loads
If is constant on , ; otherwise integrate or quadrature the load. Assemble by adding element contributions with the same global vertex labels and eliminate prescribed boundary coefficients. Only overlapping basis supports interact, making sparse. This gives the Ritz-Galerkin equivalence for a symmetric coercive form in a practical finite element form.
Symmetry implies that all eigenvalues of lie between and , so . In particular the bilinear form from part (i) is bounded and symmetric under the printed assumptions. The symmetric elliptic Dirichlet energy with a nonpositive potential is
All its terms are finite. For an admissible variation , its first variation is . Thus stationarity is exactly the weak formulation of .
Here is the direct method in the calculus of variations. Put . Expanding the symmetric quadratic form gives
The Poincaré inequality controls by . Therefore is bounded below on the admissible class, and a minimizing sequence is bounded in . By weak compactness in this reflexive Banach space, a subsequence converges weakly to . The affine closed subspace is weakly closed, so remains admissible.
The symmetric square root is measurable and bounded. Since , the two bounded linear maps
take weak convergence to weak convergence. The weak lower semicontinuity of the Hilbert norm shows that is weakly lower semicontinuous. The term is weakly continuous. Hence and attains the minimum. Varying in either sign yields the weak equation.
The sign was used both to make the potential energy nonnegative, hence weakly lower semicontinuous by this argument, and to obtain coercivity. Finally for , since the first variation at vanishes. Therefore
There is a qualification in the printed hypotheses: quadratic ellipticity does not bound a nonsymmetric coefficient matrix. The intended standard existence proof needs , or more generally a bounded principal bilinear form. To see why this does not follow from the displayed inequalities, take and
Then for every . However, choose equal to , respectively, on a small rectangle about the origin. There , whose positive part has infinite integral. Thus the standard principal pairing is not a finite Lebesgue integral even for these smooth test functions. This demonstrates the missing form hypothesis, rather than assuming that all unbounded skew coefficients prevent solvability.
Under the intended bounded-coefficient hypothesis the proof is as follows. On use the Hilbert space norm , equivalent to the norm by the Poincaré inequality on a bounded domain. The form from part (i) is bounded, since
It is a coercive bilinear form even without symmetry:
Here the nonpositive sign of is essential. With , the required equation becomes
The Cauchy-Schwarz inequality, bounded coefficients and the Poincaré inequality show that is a bounded linear functional on . The Lax-Milgram theorem states that a bounded coercive bilinear form on a real Hilbert space represents each bounded linear functional in this way with a unique . It therefore yields with the required boundary condition and weak equation. For completeness, if are solutions, satisfies ; coercivity and the Poincaré inequality imply .
With the missing bounded-form hypothesis, the weak solution exists and is unique. For nonsymmetric matrices the printed quadratic inequalities omit the hypothesis needed for this standard formulation and proof. No maximum principle is used. In part (iii), symmetry makes the quadratic upper bound a bound on the whole matrix, so no such qualification is needed there.
The affine class is nonempty, and part (i) gives a finite comparison energy . Since , its energy infimum satisfies . Choose a minimizing sequence with . Then
The gradient bounds also control the full Sobolev norm. Indeed , so the Poincaré inequality gives
Thus is bounded in the Hilbert space . Its weak compactness yields a weakly convergent subsequence . The closed linear subspace is weakly closed, hence and . By part (ii),
The direct method in the calculus of variations therefore proves
Neither a smooth boundary nor a classical trace theorem is needed: the boundary condition is encoded by the closed zero-boundary Sobolev space.
Use a planar flux function and write
This automatically satisfies the solenoidal magnetic-field constraint. There is a sign convention to settle: with the usual definition , the magnetic vector potential is . The positive sign printed for instead gives the opposite field, or corresponds to defining . We use the standard curl convention and the flux oriented along the prescribed positive ; the flux evolution and all transport identities below are unchanged by reversing the overall potential convention.
For planar incompressible flow, substitute this representation into the resistive induction equation. The two in-plane components imply
An additive time-dependent constant in leaves the field unchanged and removes . Since , this gives
This is the scalar advection-diffusion equation underlying the planar anti-dynamo theorem, with the imposed-field source retained.
Put , and . Stationarity of the horizontal averages, , and incompressibility give
Here the horizontal derivative of the advective flux has zero average. Integration gives . Averaging once more, with , identifies . Thus
Decay or periodic averaging removes the endpoint contribution from the mean derivative. The periodic assumption used later replaces decay at vertical infinity; a nonzero periodic function cannot also decay there.
Subtract the averaged equation from the scalar equation and multiply by . The average of the time derivative is zero by stationarity; the advective term becomes a boundary flux and vanishes by incompressibility and the stated boundary averaging. Terms depending only on multiply . Integration by parts in the diffusion term therefore gives
Substituting and using yields the Zeldovich planar flux balance
Both terms on the right are dissipative. In a nontrivial stationary forced problem this requires transport opposing the imposed flux gradient.
For finite-energy localized fields, enforcing decay at both vertical infinities literally gives from the mean-flux boundary values; the corresponding integral balance then forces both dissipation terms to vanish. The later nonzero-transport, periodic case therefore replaces that localized boundary setting rather than satisfying both simultaneously.
Now assume and a positive turbulent magnetic diffusivity . Define
Since and , its spatial mean is . Hence
The Zeldovich planar flux balance becomes
Because , the fluctuation cannot vanish, and a nontrivial fluctuation has under the stated boundary or zero-mean periodic conditions. Therefore
If , then : the transport profile must be nearly constant in mean square. This does not, without additional regularity, imply a uniform pointwise bound on narrow spikes.
For the final estimate use the supplied Poincaré inequality in the periodic setting and the Cauchy-Schwarz inequality:
Insert the expression for , cancel , and multiply by . Setting gives
Consequently the variance bound for planar turbulent magnetic transport is
The factor comes from maximizing , so the weaker supplied Poincare constant already suffices. Under a zero vertical mean at each horizontal location, the usual period- Fourier bound is even stronger, with in place of .
Multiply by a test function vanishing at the endpoints and apply integration by parts. The positive operator is , so the forcing changes sign. The weak formulation of this forced Airy boundary-value problem is: find such that
Here is the zero-boundary Sobolev space and can be equipped with . The interval Poincaré inequality states , so this norm makes a Hilbert space and is equivalent to its usual Sobolev norm. The form is symmetric and satisfies
The Lax-Milgram theorem says that a bounded coercive bilinear form on a Hilbert space and a bounded linear functional determine exactly one weak solution. All its hypotheses have just been verified. Moreover, gives on the interval; the right-hand side is then continuous, so this weak solution is the classical solution with the prescribed endpoint values.
Since the form is symmetric, the equivalent Ritz method minimizes
Indeed, stationarity gives the weak equation and
with equality only at . This proves the minimization equivalence, including the sign of the linear term.
Restricting to a conforming finite element space gives a unique discrete solution and Galerkin orthogonality, for every . The Céa lemma states that if the continuity and coercivity constants are and , then . For this symmetric problem the exact best-approximation constant in the energy norm is one. The piecewise-linear finite element interpolation estimate on regular finite element meshes therefore gives convergence of a numerical method of the Ritz method, with first-order error in the energy norm.
The real Lax-Milgram theorem states that if is a Hilbert space, is a continuous bilinear form with for some , and is continuous, then there is a unique such that for all . Moreover .
Take the Sobolev space with norm . Since is bounded, the Poincaré inequality gives , making this norm equivalent to the full Sobolev norm. For the given bounded measurable coefficient define
Then , since , and . The Lax-Milgram theorem therefore gives a unique satisfying this identity, which is exactly the required weak solution. Strict positivity of the lower coefficient bound is unnecessary: the zero boundary condition and boundedness supply coercivity through the gradient term.
For a bounded connected domain with the interior cone property, and , the mean-zero form of the Poincaré inequality is
No zero boundary trace is imposed. Subtracting the mean is essential because nonzero constants have zero gradient.
Suppose the estimate failed. Subtract means and normalize to obtain with
This sequence is bounded in the Sobolev space . The allowed Rellich-Kondrachov compactness theorem on cone domains provides a subsequence converging strongly in to . Its mean is zero and its norm is one. For any compactly supported smooth test function ,
Thus every weak derivative of vanishes. To show this implies constancy, mollify on any ball compactly contained in . The resulting smooth function has zero gradient and is constant on the smaller ball. Passing to its limit makes constant almost everywhere on that ball. Overlapping balls have the same constant, and connectedness connects any two points by a chain of such balls. Hence is constant almost everywhere throughout . Its zero mean makes it zero, contradicting norm one. This proves the stated Poincaré inequality. The frequently used stronger Sobolev-Poincare formulation also follows: with , the cone-domain Sobolev embedding theorem gives . Insert and the estimate just proved to obtain
Thus both the same-exponent and critical-exponent mean-zero versions are covered. Connectedness cannot be dropped without subtracting a separate mean on every component. The usual meaning of the stated range is ; with is not part of this Sobolev-space formulation.
Testing a weak Neumann Poisson problem with the constant function one gives , proving necessity. For sufficiency, use the closed Hilbert space of mean-zero functions. The given Poincaré inequality makes an equivalent norm there. Set
The negative sign corresponds to . Poincare's inequality gives boundedness of and coercivity of . The Lax-Milgram theorem supplies a unique mean-zero solution. Every is the sum of its mean-zero part and a constant; since , the weak identity extends to every . This is precisely the weak equation with zero normal derivative as its natural boundary condition. Existence and uniqueness in the mean-zero space hold exactly when .
For , both and its weak gradient are square integrable. The assumed bounds give , so the potential term is integrable and is finite. Also
which is bounded below and coercive in by the Poincaré inequality.
For distinct and , convexity of and expansion of the squared gradient give
The last norm is nonzero: a zero weak gradient makes constant, and its zero trace forces that constant to vanish. Thus is strictly convex.
Weak lower semicontinuity means that in implies . A minimizing sequence is bounded by coercivity. Reflexivity of this Hilbert space gives a weakly convergent subsequence with limit still in . The assumed weak lower semicontinuity makes that limit a minimizer, and strict convexity makes it unique. Hence .
For a variation preserving both endpoint values and slopes, at the endpoints. The first variation is
Integrate the last term twice and the middle term once. All boundary terms vanish, giving
The fundamental lemma of the calculus of variations yields the higher-order Euler-Lagrange equation
For the present functional this becomes . Its general solution is . The left endpoint gives and . The right endpoint then gives
The coefficient determinant is , so , . The stationary function is .
For any admissible perturbation , exact quadratic expansion, followed by the stationarity equation and integration by parts, gives
Apply the supplied Poincaré inequality first to , then to , whose endpoint values also vanish:
Thus for every admissible perturbation, proving a global maximum through the clamped second-variation inequality.
The appropriate bounded-operator formulation uses the zero-boundary Sobolev space
The Poincaré inequality makes this a Hilbert space norm equivalent to the usual norm, and the zero endpoint traces encode the Dirichlet boundary conditions. Under the usual regular-coefficient assumptions, for example continuous on the closed interval with and , set
If , then Cauchy-Schwarz inequality and the Poincaré inequality give
By the Riesz representation theorem, a unique bounded linear operator satisfies . Symmetry of makes self-adjoint, and the lower bound makes it elliptic and uniformly positive definite. The differential expression in the question is represented weakly by : for smooth functions, integration by parts gives
with the boundary term zero. A forcing becomes the bounded functional , or its Riesz representative in .
Equivalently, for sufficiently smooth coefficients one may realize the differential operator itself on with domain . Its regular Sturm-Liouville operator realization is self-adjoint and
This realization is an unbounded operator, so it should not be confused with the bounded weak operator used above.
The PDF states sign conditions without coefficient regularity. If only almost everywhere and are bounded, the same form is still bounded and strictly positive: a zero energy would force almost everywhere and the zero traces force . Uniform ellipticity, however, requires a positive lower bound. For instance for , , and satisfy literal pointwise positivity. Derivatives of unit norm supported in and of integral zero define functions in , with energy at most . Thus pointwise positivity without regularity does not imply coercivity in this norm. These distinctions supply the precise regularity and operator domain behind the intended positive-definiteness proof.
For , two integrations by parts give the clamped Hessian identity
By density it remains valid on . Each has zero integral, first for compactly supported test functions and then by convergence. Applying the Neumann-Poincare inequality to each gives
The zero-boundary Poincaré inequality also gives . Hence, using a full-Hessian equivalent norm,
Thus is an inner product whose norm is equivalent to the complete norm on the clamped second-order Sobolev space. The functional is bounded for this norm by the Cauchy-Schwarz inequality and the displayed bound. Apply the Riesz representation theorem, or the Lax-Milgram theorem, to get a unique representing . The clamped biharmonic problem has a unique weak solution for every , with and no zero-integral compatibility condition.
The difference of two weak solutions lies in the clamped second-order Sobolev space. Testing with yields , so . Since , integration by parts gives
The Poincaré inequality for zero boundary values now implies . The clamped biharmonic problem has at most one weak solution. Unlike the Neumann Poisson problem, no additive constant is allowed by these boundary traces.
Take real and use the L2 norm on the spatial interval. Existence, uniqueness and continuous dependence are the three requirements of Hadamard well-posedness. An energy method supplies the decisive estimate. For a smooth solution with homogeneous Dirichlet boundary conditions, integration by parts gives
The drift contributes only a boundary term, which vanishes. The Poincaré inequality further gives
Apply the same argument to the difference of two solutions to obtain uniqueness and continuous dependence on the initial data.
For existence, use the Dirichlet gauge transform for constant drift: satisfies with zero boundary values. Expanding in its Fourier sine series gives
For the series defines a solution continuous in down to and smooth for positive time; multiplication by the fixed bounded exponentials preserves this interpretation. Its energy estimate follows by approximation with smooth initial data. For a classical solution at the initial corners, require the usual smoothness and boundary compatibility instead. The problem is well posed in , with a contraction estimate independent of the initial data.
The finite element method replaces an infinite-dimensional variational problem by a finite-dimensional one, usually using functions that are polynomial on small mesh elements. Ritz method and Galerkin method describe how the discrete equations are selected; neither intrinsically requires piecewise polynomials, though those spaces make assembly local and sparse.
For a concrete realization of the two-point problem, choose homogeneous Dirichlet boundary conditions at both endpoints of . The differential expression alone is not a complete boundary value problem, so this boundary choice is an explicit assumption. Let with , with , and . Use the zero-boundary Sobolev space with norm . By the Poincaré inequality this is equivalent to its full first-order Sobolev norm. Multiplying by a test function and applying integration by parts gives
The endpoint terms vanish by the essential boundary condition. The form is a symmetric bilinear form, bounded by
and it is a coercive bilinear form:
The load is a bounded linear functional. The Lax-Milgram theorem gives a unique weak solution and a bound . Under the stated coefficient regularity it solves the differential equation in the usual weak sense.
The Ritz method minimizes the energy
over , or over a conforming finite-dimensional subspace . Its first variation is ; symmetry and coercivity make the functional strictly convex. In fact, if is the weak solution,
Thus the energy minimizer is exactly the weak solution. The Galerkin method instead directly asks for with for all . For this symmetric coercive problem, Ritz-Galerkin equivalence for a symmetric coercive form shows that the discrete minimizer and Galerkin solution coincide. For a nonsymmetric problem, a Galerkin formulation can still apply but the same quadratic minimization generally cannot represent its full bilinear form; appropriate coercivity or inf-sup conditions must then justify the chosen spaces.
Subtract the continuous and discrete weak equations to obtain Galerkin orthogonality . Since defines the energy norm , the Pythagorean identity gives
In a general reference norm, the Céa lemma yields
The proof uses together with coercivity and continuity. Dense approximation spaces therefore imply convergence; the estimate reduces a numerical error problem to an approximation problem.
For implementation, partition the interval by nodes . Choose continuous piecewise-linear functions and the interior piecewise-linear hat functions as a basis, with the two boundary coefficients prescribed. Write . The coefficients satisfy
The stiffness matrix is symmetric positive definite: for nonzero interior coefficients. The local supports make it tridiagonal in this one-dimensional linear-element case.
On an element of length , for constant element coefficients , the two-node matrix and constant-load vector are
For variable coefficients, integrate against the same shape functions rather than silently replacing them by constants. Element contributions are added at shared nodes, and boundary values are eliminated or lifted into the right side. A sparse direct solve or a suitable iterative method gives the nodal coefficients. Numerical quadrature should preserve the required accuracy and, with positive coefficients and weights, the positive energy structure.
A finite element interpolation estimate gives for these elements, where is the largest element size and the solution has the stated regularity. The Céa lemma then gives first-order energy error. If the dual elliptic problem has regularity, the Aubin–Nitsche duality argument gives second-order error: for , solve , then use to gain one further factor of . Higher-order elements improve rates when the solution is sufficiently smooth.
Nonzero Dirichlet data are handled by a boundary lifting and an affine trial space. Neumann boundary conditions enter naturally through the integrated boundary term; Robin terms modify both the form and the load. With pure Neumann conditions and , constants lie in the kernel: existence needs the load/flux compatibility condition, uniqueness needs a mean constraint, and the preceding Dirichlet coercivity argument cannot simply be reused. Conforming approximation, boundary conditions and coercivity are the ingredients connecting the finite-element construction to a justified Ritz or Galerkin solution.
Use the Dirichlet Laplacian eigenvalues and eigenfunctions on :
An expansion reduces the wave equation to independent ordinary differential equations
When , put ; each coefficient is
For finite-energy data and , a uniform spatial bound follows from the energy method, rather than from unjustified absolute summation of the Fourier series. The conserved quantity is
The Poincaré inequality gives . Hence with for , and with for . Since the Dirichlet boundary condition gives , the Cauchy-Schwarz inequality yields
At , the first coefficient is , which is unbounded for some admissible data. For , that coefficient has an exponentially growing component for some admissible data. Thus all-time boundedness of a wave equation with a reaction term requires
At equality, boundedness for a particular data set requires ; the remaining modes are bounded by their spectral gap. Above the threshold, every unstable mode must have its growing component canceled, and any zero-frequency mode must have zero initial velocity. There is no condition on alone for arbitrary specially chosen data.
The printed reference to a limit needs qualification. Bounded oscillations generally have no limit as . If actual existence of that limit for every admissible initial datum is required, no real works: for every , choose with and nonzero pure oscillatory data in that eigenfunction. The boxed inequality answers the intended long-time boundedness question.