Laplace's Equation — Question 8

PDF ↗

Question 8

A unit square has insulated vertical sides and prescribed horizontal temperatures: Δu=0,ux(0,y)=ux(1,y)=0,u(x,0)=0,u(x,1)=1+cos⁡(2πx).\Delta u=0,\qquad u_x(0,y)=u_x(1,y)=0,\qquad u(x,0)=0,\quad u(x,1)=1+\cos(2\pi x). Seek a solution smooth up to the boundary. Conductivity is one.

Tasks

  1. Construct the solution by separation, treating the zero spatial eigenvalue explicitly.

  2. Verify every datum and prove uniqueness for this mixed boundary problem using Green’s identity.

  3. Compute the horizontal average at each height and the integrated outward heat flux on each edge. Explain why zero side flux does not force constant temperature.

  4. Calculate the Dirichlet energy. Prove 0<u<20<u<2 at every interior point, justifying how the insulated sides are handled.

Original worksheet page 1: question and worked solution for 9-7-008
Show solutionHide solution

Question 8 – Solution

Strategy. Neumann side conditions select cosine modes, including a zero mode whose vertical factor is linear rather than hyperbolic.

Step 1: Keep the zero mode. The side eigenfunctions are 1,cos⁡(nπx)1,\cos(n\pi x) with eigenvalues 0,(nπ)20,(n\pi)^2. The mean mode satisfies Y″=0Y''=0 with bottom value zero and top value one, so it is yy. The second cosine mode gives u(x,y)=y+sinh⁡(2πy)sinh⁡(2π)cos⁡(2πx).\boxed{u(x,y)=y+\frac{\sinh(2\pi y)}{\sinh(2\pi)}\cos(2\pi x).} Dropping the zero mode would lose the entire prescribed mean temperature.

Step 2: Verify and prove mixed uniqueness. The displayed terms are harmonic. The sine factor in uxu_x vanishes at both vertical sides; the bottom and top traces are exactly zero and 1+cos⁡(2πx)1+\cos(2\pi x). For the difference ww of two solutions, Green’s identity yields ∫|∇w|2=∫∂Ωw∂nw=0\int|\nabla w|^2=\int_{\partial\Omega}w\partial_nw=0: w=0w=0 on the horizontal edges and ∂nw=0\partial_nw=0 on the vertical ones. Thus ww is constant, and its bottom value forces w=0w=0.

Step 3: Track the mean and edge fluxes. The cosine has zero horizontal integral, so u¯(y)=∫01u(x,y)dx=y\overline u(y)=\int_0^1u(x,y)\,dx=y. The integrated outward conductive fluxes are Ftop=−1,Fbottom=1,Fleft=Fright=0.\boxed{F_{\mathrm{top}}=-1,\qquad F_{\mathrm{bottom}}=1,\qquad F_{\mathrm{left}}=F_{\mathrm{right}}=0.} They balance. Insulation forbids heat crossing the vertical sides; it does not forbid a vertical gradient or heat transfer between the horizontal edges.

Step 4: Compute energy and bound the field. Only the top contributes to ∫∂Ωu∂nu\int_{\partial\Omega}u\partial_nu. At the top, uy=1+2πcoth⁡(2π)cos⁡(2πx)u_y=1+2\pi\coth(2\pi)\cos(2\pi x), so orthogonality gives ℰ=1+πcoth⁡(2π).\boxed{\mathcal E=1+\pi\coth(2\pi).} For 0<y<10<y<1, strict convexity of sinh⁡(2πy)\sinh(2\pi y) gives 0<sinh⁡(2πy)/sinh⁡(2π)<y0<\sinh(2\pi y)/\sinh(2\pi)<y. Therefore 0<y−sinh⁡(2πy)/sinh⁡(2π)≤u≤y+sinh⁡(2πy)/sinh⁡(2π)<2y<20<y-\sinh(2\pi y)/\sinh(2\pi)\le u \le y+\sinh(2\pi y)/\sinh(2\pi)<2y<2. This direct estimate avoids incorrectly treating the unspecified side temperatures as Dirichlet data.

Original worksheet page 2: question and worked solution for 9-7-008

Original worksheet layout. Use Enlarge or open the PDF for a closer view.