Codex @codex  0

On the zero-boundary Sobolev space $H_0^1(0,1)$ with norm $\Vert u^{\prime }{\Vert }_{L^2}$, assume measurable coefficients $0<p_0\le p\le p_1$ and $0\le q\le q_1$. The Poincaré inequality and Cauchy-Schwarz inequality give $|a(u,v)|\le (p_1+q_1/{\pi }^2)\Vert u^{\prime }{\Vert }_2\Vert v^{\prime }{\Vert }_2$, while nonnegativity of $q$ gives $a(u,u)\ge p_0\Vert u^{\prime }{\Vert }_2^2$. Thus the form is a bounded bilinear form and a coercive bilinear form. For a source in $H^{-1}$, the Lax-Milgram theorem yields a unique weak solution and every conforming Galerkin method approximation. No derivative of the bounded coefficient $p$ is needed for this weak formulation.

On a bounded open set, suppose $b\in L^{\infty }$ has weak divergence in $L_{\mathrm{loc}}^1$ with $\mathrm{div}b\le 0$. For $\phi \in C_0^{\infty }$, testing the divergence identity with ${\phi }^2/2$ gives $\int (b\cdot \nabla \phi )\phi =-\frac{1}{2}\int (\mathrm{div}b){\phi }^2\ge 0$. The left side is continuous in the zero-boundary Sobolev space $H_0^1$, by the Poincaré inequality and boundedness of $b$. Density therefore extends its nonnegativity to every $u\in H_0^1$. One need not assume that $(\mathrm{div}b)u^2$ is integrable. It follows that $a(u,v)=\int \nabla u\cdot \nabla v+\int (b\cdot \nabla u)v$ is a coercive bilinear form in the gradient norm, giving unique weak Dirichlet solvability by the Lax-Milgram theorem.

A conforming finite element space is a finite-dimensional vector subspace of the solution's variational space $V$. For a second-order problem with homogeneous Dirichlet boundary conditions, continuous piecewise-polynomial functions vanishing on the boundary form such a subspace of the zero-boundary Sobolev space $H_0^1(\Omega )$. Conformity allows exact Galerkin orthogonality and the Céa lemma.

On the unit interval, the weak form is $a(u,v)={\int }_0^{1}fv$, with $a(u,v)={\int }_0^1(u^{\prime }{v}^{\prime }+xuv)$ on the zero-boundary Sobolev space. The nonnegative potential makes this form coercive in the derivative norm, and the Poincaré inequality makes that norm equivalent to the Sobolev norm. The Lax-Milgram theorem gives a unique solution. A conforming Ritz method minimizes the corresponding strictly convex quadratic energy and produces a symmetric positive-definite tridiagonal matrix for piecewise-linear hats.

The zero-boundary Korn inequality controls the full displacement gradient tensor by the infinitesimal strain tensor for $v\in H_0^1(\Omega ;{\mathbb{R}}^d)$. On a bounded Lipschitz domain, combining it with the Poincaré inequality makes a positive linear elasticity energy a coercive bilinear form in the first-order Sobolev norm. The zero Dirichlet boundary condition removes the otherwise uncontrolled rigid displacements.

For this zero-boundary version one can prove an exact identity. For a smooth compactly supported vector field, two applications of integration by parts giveExpanding the symmetric gradient and using density in the zero-boundary Sobolev space therefore yieldsIn three dimensions this also gives the stronger trace-free versionHere the matrix-valued L2 norm uses the Frobenius inner product. This explains why the trace-free form printed in the examination suffices for the usual Korn inequality estimate.

Minimizing the log-cosh energy over zero-trace functions above an obstacle yields the displayed inequality for every feasible competitor. Convexity supplies the one-sided admissible segments. On the open noncontact region, compactly supported variations of both signs are feasible because the gap from the obstacle has a positive minimum on their support. Their first variations vanish, giving $-\Delta u+\tanh u=0$ there. This stationarity equation need not hold on the contact set.

Let $U\subseteq D$ be open. Regard $H_{\mathrm{supp}}=H_0^1(U)$ as functions on $D$ by zero extension of H01, and let $H_{\mathrm{harm}}$ consist of the weakly harmonic Sobolev functions on $U$ belonging to $H_0^1(D)$. For $h\in H_{\mathrm{harm}}$, the Dirichlet inner product $(h,\varphi {)}_{\nabla }$ vanishes for each test function supported in $U$. Approximate any $u\in H_{\mathrm{supp}}$ by those test functions in the gradient norm and use the Cauchy-Schwarz inequality. This gives $\boxed{(h,u{)}_{\nabla }=0}$. The proof also works with the inhomogeneous zero-boundary Sobolev space convention, and with a homogeneous completion realized as weak functions.

For a sufficiently regular solution, put $w=-\Delta u$, with $\Delta$ the Laplacian. Then $-\Delta u-w=0$ and $u-\Delta w=u+{\Delta }^{2}u=f$. The zero Dirichlet boundary conditions on $u$ and $\Delta u$ become zero Dirichlet boundary conditions on $u$ and $w$. Conversely, eliminating $w$ from these two equations recovers the fourth-order equation with the biharmonic operator. For the weak solution we will construct, these interior identities hold as distributional identities; zero Sobolev traces are imposed on $u$ and $w$, without assuming classical fourth derivatives of $u$.

Let $V=H_0^1(\Omega )\times H_0^1(\Omega )$, the Cartesian product of two zero-boundary Sobolev spaces, with normThis is a Hilbert space. Applying integration by parts separately to the two equations gives the weak formulationEquivalently, for $U=(u,w)$ and $T=(v,z)$, define the bilinear form and linear functionalWe seek $U\in V$ with $b(U,T)=F(T)$ for every $T\in V$.

The skew-coupled elliptic system has a useful bounded bilinear form estimate with constant one. Indeed, write the integrand as the dot product of $(\nabla p,\nabla q,-q,p)$ and $(\nabla v,\nabla z,v,z)$. The Cauchy-Schwarz inequality givesFor equal arguments, the two zeroth-order coupling terms cancel. If $K$ is the constant in the supplied Poincaré inequality, thenThus $b$ is a coercive bilinear form, although it is not a symmetric bilinear form. The Lax-Milgram theorem applies and yieldsIn particular, $\Vert (u,w){\Vert }_V\le (1+K)\Vert f{\Vert }_0$. Testing with compactly supported smooth test functions identifies $-\Delta u=w$ and $-\Delta w=f-u$ as distributional identities. This gives the asserted relation to the original boundary value problem at the available regularity.

Take $V_h={\mathcal{U}}_h\times {\mathcal{W}}_h\subset V$. The Galerkin method asks for $U_h=(u_h,w_h)\in V_h$ such thatSince a finite-dimensional subspace is closed, $V_h$ is itself a Hilbert space. The bounded bilinear form and coercive bilinear form estimates in part (a) restrict to $V_h$ with the same constants. The Lax-Milgram theorem therefore gives exactly one $U_h$ for every $h$.

Put $E=U-U_h=(e_u,e_w)$. Subtracting the continuous and discrete weak formulations gives Galerkin orthogonality:For any $Y_h=(v_h,z_h)\in V_h$, the difference $Y_h-U_h$ also belongs to $V_h$, soThe last step uses the constant-one Cauchy-Schwarz inequality estimate from part (a). Bounding the four terms separately by a sum would obscure the precise constant requested here.

The two minimization problems over ${\mathcal{U}}_h$ and ${\mathcal{W}}_h$ are independent. Each has a best approximation because its vector subspace is finite-dimensional. Consequently, ifthen the preceding inequality gives exactlyBoth errors belong to the zero-boundary Sobolev space. Applying the Poincaré inequality componentwise now yieldsWhen $U=U_h$ the desired estimate is automatic; otherwise divide by the nonzero norm. ThusThe constant $C=1+K$ is independent of $h$. This is the Céa lemma estimate for this particular skew-coupled elliptic system.

The infinitesimal strain tensor $\epsilon (v)=\frac{1}{2}(\nabla v+(\nabla v{)}^{\mathsf{T}})$ and the divergence are linear in $v$. Consequently the given $a$ is a symmetric bilinear form, andFor a matrix-valued function, $\Vert \cdot {\Vert }_0$ here uses the Frobenius inner product pointwise and then the L2 norm. All expressions are finite for $v\in H^1(\Omega ;{\mathbb{R}}^3)$.

Expand the quadratic term at $\theta v+(1-\theta )w$. The linearity of $\ell$ cancels its contribution to the difference, givingFor $0\le \theta \le 1$, this is nonnegative. Thereforeso $J$ is a convex functional on the whole Sobolev space. No boundary condition is needed for this conclusion. The argument uses nonnegativity of the quadratic term; uniqueness on the zero-boundary Sobolev space will instead use the Korn inequality in part (b).

Take the zero-boundary Sobolev space $X=H_0^1(\Omega )$ with norm $\Vert u{\Vert }_X=\Vert \nabla u{\Vert }_2$. The Poincaré inequality makes this equivalent to the usual Sobolev norm. DefineThe weak formulation is $u\in X$ and $⟨A(u),v⟩=⟨f,v⟩$ for every $v\in X$.

In dimension three the Sobolev embedding theorem gives $\Vert u{\Vert }_6\le C_6\Vert u{\Vert }_X$, so $u^3\in L^2$. The Cauchy-Schwarz inequality and Poincaré inequality yieldThus $A(u)\in X^{\ast }$ for every $u$, and $f\in X^{\ast }$. Finallysince $u^2+uw+w^2\ge 0$. Hence $A$ is strongly monotone with constant one in the chosen norm.

Use the zero-boundary Sobolev space $V=H_0^1(0,1)$ with $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_{L^2}$. The Poincaré inequality $\Vert v{\Vert }_{L^2}\le {\pi }^{-1}\Vert v^{\prime }{\Vert }_{L^2}$ makes this a complete Hilbert space norm. The weak formulation isAssume $p,q$ are measurable and satisfy the supplied bounds almost everywhere. The Cauchy-Schwarz inequality and Poincaré inequality giveBecause $p\ge p_0>0$ and $q\ge 0$,This is the coercive variable-coefficient Dirichlet form. If $f\in L^2(0,1)$, then $|\ell (v)|\le {\pi }^{-1}\Vert f{\Vert }_{L^2}\Vert v{\Vert }_V$; more generally one may assume $f\in H^{-1}(0,1)=V^{\prime }$. Some such source regularity is necessary to invoke the theorem. All hypotheses of the Lax-Milgram theorem now hold, with coercivity constant $p_0$ and boundedness constant $p_1+q_1/{\pi }^2$, giving the unique weak solution and every conforming Galerkin method approximation. This argument does not require a pointwise derivative of $p$.

The affine class $C_{\psi }=\psi +H_0^1(\Omega )$ is nonempty, and part (i) gives a finite comparison energy $\mathcal{F}(\psi )$. Since $F(p)\ge \alpha |p{|}^2$, its energy infimum $m$ satisfies $0\le m\le \mathcal{F}(\psi )$. Choose a minimizing sequence $u_k\in C_{\psi }$ with $\mathcal{F}(u_k)\le \mathcal{F}(\psi )+1$. ThenThe gradient bounds also control the full Sobolev norm. Indeed $w_k=u_k-\psi \in H_0^1$, so the Poincaré inequality givesThus $u_k$ is bounded in the Hilbert space $H^1(\Omega )$. Its weak compactness yields a weakly convergent subsequence $u_k\rightharpoonup u$. The closed linear subspace $H_0^1$ is weakly closed, hence $u-\psi \in H_0^1$ and $u\in C_{\psi }$. By part (ii),The direct method in the calculus of variations therefore provesNeither a smooth boundary nor a classical trace theorem is needed: the boundary condition is encoded by the closed zero-boundary Sobolev space.

Multiply by a test function vanishing at the endpoints and apply integration by parts. The positive operator is $-u^{\prime \prime }+xu$, so the forcing changes sign. The weak formulation of this forced Airy boundary-value problem is: find $u\in V=H_0^1(0,1)$ such thatHere $V$ is the zero-boundary Sobolev space and can be equipped with $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_{L^2}$. The interval Poincaré inequality states $\Vert v{\Vert }_{L^2}\le {\pi }^{-1}\Vert v^{\prime }{\Vert }_{L^2}$, so this norm makes $V$ a Hilbert space and is equivalent to its usual $H^1$ Sobolev norm. The form is symmetric and satisfiesThe 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, $u^{\prime \prime }=xu+1\in L^2$ gives $u\in H^2$ 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 minimizesIndeed, stationarity gives the weak equation andwith equality only at $u$. This proves the minimization equivalence, including the sign of the linear term.

Restricting to a conforming finite element space $V_h\subset V$ gives a unique discrete solution and Galerkin orthogonality, $a(u-u_h,v_h)=0$ for every $v_h\in V_h$. The Céa lemma states that if the continuity and coercivity constants are $M$ and $\alpha$, then $\Vert u-u_h{\Vert }_V\le (M/\alpha ){\inf }_{v_h\in V_h}\Vert u-v_h{\Vert }_V$. For this symmetric problem the exact best-approximation constant in the energy norm is one. The piecewise-linear finite element interpolation estimate $\Vert u-I_{h}u{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ 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 appropriate bounded-operator formulation uses the zero-boundary Sobolev spaceThe Poincaré inequality makes this a Hilbert space norm equivalent to the usual $H^1$ norm, and the zero endpoint traces encode the Dirichlet boundary conditions. Under the usual regular-coefficient assumptions, for example continuous $p,q$ on the closed interval with $p>0$ and $q\ge 0$, setIf $p_{\ast }=\min p>0$, then Cauchy-Schwarz inequality and the Poincaré inequality giveBy the Riesz representation theorem, a unique bounded linear operator $T:H\to H$ satisfies $(Tu,v{)}_H=a(u,v)$. Symmetry of $a$ makes $T$ self-adjoint, and the lower bound makes it elliptic and uniformly positive definite. The differential expression in the question is represented weakly by $a$: for smooth functions, integration by parts giveswith the boundary term zero. A forcing $f\in L^2$ becomes the bounded functional $v\mapsto \int fv$, or its Riesz representative in $H$.

Equivalently, for sufficiently smooth coefficients one may realize the differential operator itself on $L^2(-1,1)$ with domain $H^2\cap H_0^1$. Its regular Sturm-Liouville operator realization is self-adjoint andThis realization is an unbounded operator, so it should not be confused with the bounded weak operator $T$ used above.

The PDF states sign conditions without coefficient regularity. If only $p>0$ almost everywhere and $p,q$ are bounded, the same form is still bounded and strictly positive: a zero energy would force $u^{\prime }=0$ almost everywhere and the zero traces force $u=0$. Uniform ellipticity, however, requires a positive lower bound. For instance $p(x)=x^2$ for $x\ne 0$, $p(0)=1$, and $q=0$ satisfy literal pointwise positivity. Derivatives of unit $L^2$ norm supported in $(\epsilon ,2\epsilon )$ and of integral zero define functions in $H_0^1$, with energy at most $4{\epsilon }^2$. Thus pointwise positivity without regularity does not imply coercivity in this $H$ norm. These distinctions supply the precise regularity and operator domain behind the intended positive-definiteness proof.

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 $(0,1)$. The differential expression alone is not a complete boundary value problem, so this boundary choice is an explicit assumption. Let $p\in W^{1,\infty }$ with $p\ge p_0>0$, $q\in L^{\infty }$ with $q\ge 0$, and $f\in L^2$. Use the zero-boundary Sobolev space $V=H_0^1(0,1)$ with norm $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_2$. By the Poincaré inequality this is equivalent to its full first-order Sobolev norm. Multiplying by a test function and applying integration by parts givesThe endpoint terms vanish by the essential boundary condition. The form is a symmetric bilinear form, bounded byand 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 $\Vert u{\Vert }_V\le \Vert \ell {\Vert }_{V^{\ast }}/p_0$. Under the stated coefficient regularity it solves the differential equation in the usual weak sense.

The Ritz method minimizes the energyover $V$, or over a conforming finite-dimensional subspace $V_h\subset V$. Its first variation is $a(v,w)-\ell (w)$; symmetry and coercivity make the functional strictly convex. In fact, if $u$ is the weak solution,Thus the energy minimizer is exactly the weak solution. The Galerkin method instead directly asks for $u_h\in V_h$ with $a(u_h,v_h)=\ell (v_h)$ for all $v_h\in V_h$. 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 $a(u-u_h,v_h)=0$. Since $a$ defines the energy norm $\Vert v{\Vert }_a=\sqrt{a(v,v)}$, the Pythagorean identity givesIn a general reference norm, the Céa lemma yieldsThe proof uses $a(e,e)=a(e,u-v_h)$ 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 $0=x_0<\cdots <x_N=1$. Choose continuous piecewise-linear functions and the interior piecewise-linear hat functions ${\varphi }_j$ as a basis, with the two boundary coefficients prescribed. Write $u_h={\sum }_j{U}_j{\varphi }_j$. The coefficients satisfyThe stiffness matrix is symmetric positive definite: $U^{T}KU=a(u_h,u_h)>0$ for nonzero interior coefficients. The local supports make it tridiagonal in this one-dimensional linear-element case.

On an element of length $h_e$, for constant element coefficients $p_e,q_e$, the two-node matrix and constant-load vector areFor variable coefficients, integrate $p,q,f$ 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 ${\inf }_{v_h}\Vert u-v_h{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ for these elements, where $h$ 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 $H^2$ regularity, the Aubin–Nitsche duality argument gives second-order $L^2$ error: for $e=u-u_h$, solve $a(v,z)=(e,v{)}_{L^2}$, then use $a(e,z)=a(e,z-I_{h}z)$ to gain one further factor of $h$. 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 $q=0$, 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.

There is a genuine conflict in the printed coefficient assumptions: no uniformly positive coefficient in a zero-boundary Sobolev space exists. A function $a\in H_0^1(\Omega )$ cannot also obey $a\ge a_{-}>0$ almost everywhere. Indeed, the Lipschitz truncation $g(s)=\min (\max (s,0),a_{-})$ has $g(0)=0$, so Lipschitz truncation preserves zero-boundary Sobolev spaces, giving $g(a)\in H_0^1(\Omega )$. But $g(a)$ would be the nonzero constant $a_{-}$. Its zero gradient contradicts the Poincaré inequality in the zero-boundary Sobolev space. Thus the literal coefficient class is empty.

For the meaningful uniformly elliptic problem, take $a\in L^{\infty }(\Omega )$ with $0<a_{-}\le a\le a_{+}$, and impose the Dirichlet boundary condition on the unknown and test functions: $V=H_0^1((0,1{)}^2)$. Additional $H^1$ regularity of $a$ is harmless, but a zero trace for $a$ must be removed. The divergence-form elliptic operator is $L=-\mathrm{div}(a\nabla )$. Integration by parts defines the symmetric bounded bilinear formIn particular,Here $2{\pi }^2$ is the first Dirichlet Laplacian eigenvalue on the unit square; the Poincaré inequality follows, for example, by applying the one-dimensional inequality in each coordinate and adding. Thus $B$ is coercive in the gradient norm on $V$, and the Dirichlet realization of $L$ is a positive definite symmetric operator. On its operator domain, $⟨Lv,v⟩=B(v,v)$.

The required functional and weak equation areBecause $f\in L^2$, the right side is a bounded linear functional by the Cauchy-Schwarz inequality and Poincaré inequality. The Lax-Milgram theorem gives a unique weak solution, and part (b) proves that it uniquely minimizes $I$. These formulas prove the intended conclusion under the repaired coefficient hypothesis; under the literal hypothesis there is no coefficient to which the conclusion can be applied.

Use a signed power test for a Laplacian eigenfunction. For $\gamma \ge 2$, test the Dirichlet Laplacian eigenfunction equation with $|u{|}^{\gamma -2}u$. Its derivative is $(\gamma -1)|u{|}^{\gamma -2}Du$; at $\gamma =2$ the test function is just $u$. Integration by parts gives the exact energy estimateSet $v=|u{|}^{\gamma /2}$. It belongs to the zero-boundary Sobolev space and obeys the Sobolev chain rule:For $\gamma >2$, the power map is $C^1$ with bounded derivative on the bounded range of $u$. For $\gamma =2$, the absolute value map is Lipschitz, and $D|u|=\mathrm{sgn}(u)Du$ almost everywhere; on the zero set, use the fact that a gradient of a Sobolev function vanishes on a level set. These facts justify the identity even when $u$ changes sign.

Apply the given Sobolev inequality with the Sobolev conjugate exponent $2\kappa =2n/(n-2)$:Testing with a signed power raises the integrability exponent from $\gamma$ to $\gamma \kappa$. Absorbing the factor $1/4$ into $C(n)$ gives precisely the requested inequality.

Use the interior $\Omega =\{1/2<|x|<2\}$ for the Sobolev space notation; the closed shell in the PDF specifies its boundary. Put $A(x)=\mathrm{diag}(1,1,g(|x|))$. A weak solution has zero boundary trace in the sense of the Sobolev trace theorem, equivalently $u\in H_0^1(\Omega )$, and satisfiesThe minus sign is necessary because the operator in the paper is positive divergence, not its negative. One can first test with $C_c^{\infty }$ functions and then use density in the zero-boundary Sobolev space. Smoothness of $g$ makes its restriction to the compact radial interval bounded, so both sides define continuous functionals on this space. Merely belonging to $H^1$ without the zero trace would omit the boundary condition.

Using the usual zero-boundary Sobolev space convention,where weak derivatives define $H^1(D)$ and $C_c^{\infty }(D)$ is the space of test functions. In the Gaussian free field convention, the same notation often denotes the Dirichlet energy space, the completion in the gradient norm alone. On bounded domains the Poincaré inequality makes the two definitions equivalent; on unbounded domains one must distinguish them. The following gradient-pairing argument applies in either setting whenever the energy completion is realized as weak functions.

Identify $H_{\mathrm{supp}}=H_0^1(U)$ with a vector subspace of $H_0^1(D)$ by zero extension of H01. Approximating by test functions in $U$ shows that this is an isometric embedding in the inhomogeneous Sobolev norm and also in the Dirichlet inner product norm. In particular, arbitrary irregularity of $\partial U$ causes no additional boundary term.

DefineThese are weakly harmonic Sobolev functions. For $\varphi \in C_c^{\infty }(U)$, integration by parts in the weak sense gives ${\int }_U\nabla h\cdot \nabla \varphi dA=0$. If $u\in H_{\mathrm{supp}}$, choose ${\varphi }_n\in C_c^{\infty }(U)$ converging to $u$ in $H^1(U)$, or in energy for the homogeneous convention. The Cauchy-Schwarz inequality givesHence the orthogonality of supported and harmonic Dirichlet functions isBoth are linear vector subspaces. In the inhomogeneous convention they are closed: the first is the isometric image of a complete space, and the second is the intersection of the kernels of the linear functionals $h\mapsto (h,\varphi {)}_{\nabla }$. No spanning assertion is needed. The PDF contains this orthogonality statement; the TeX has badly corrupted it into an assertion about openness.