Laplace's Equation — Question 1

PDF ↗

Question 1

A square plate occupies 0<x<10<x<1, 0<y<10<y<1. Its steady temperature satisfies uxx+uyy=0,u(0,y)=u(1,y)=u(x,0)=0,u(x,1)=sin⁡(πx)+12sin⁡(3πx).u_{xx}+u_{yy}=0,\qquad u(0,y)=u(1,y)=u(x,0)=0,\qquad u(x,1)=\sin(\pi x)+\tfrac 12\sin(3\pi x). The conductivity is one, so the outward heat flux is −∂nu-\partial_nu.

Tasks

  1. Construct the harmonic temperature by separation of variables and verify all four boundary values.

  2. At height y=1/2y=1/2, determine the attenuation of each boundary mode and their amplitude ratio. Explain which spatial feature penetrates farther into the plate.

  3. Compute the total outward heat flux on each edge and verify global balance. Use outward normals consistently.

  4. Sketch temperature profiles at y=1/4,1/2,1y=1/4,1/2,1. Prove that the temperature is strictly positive in the interior despite the sign changes of the third sine mode.

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

Question 1 – Solution

Strategy. Separate the two imposed spatial frequencies and track their normal derivatives as well as their temperatures.

Step 1: Extend each boundary mode. The zero side values give sin⁡(nπx)\sin(n\pi x); the corresponding vertical equation is Y″−(nπ)2Y=0Y''-(n\pi)^2Y=0. The bottom value selects sinh⁡(nπy)\sinh(n\pi y) and the top value normalizes it: u=sinh⁡(πy)sinh⁡πsin⁡(πx)+12sinh⁡(3πy)sinh⁡(3π)sin⁡(3πx).\boxed{u=\frac{\sinh(\pi y)}{\sinh\pi}\sin(\pi x) +\frac 12\frac{\sinh(3\pi y)}{\sinh(3\pi)}\sin(3\pi x).} Each term has zero Laplacian. Substitution gives all four traces, including compatible corner values. The maximum principle gives uniqueness among continuous harmonic extensions.

Step 2: Compare penetration of the modes. Since sinh⁡(a/2)/sinh⁡a=1/(2cosh⁡(a/2))\sinh(a/2)/\sinh a=1/(2\cosh(a/2)), the first and third modes are attenuated by 1/(2cosh⁡(π/2))1/(2\cosh(\pi/2)) and 1/(2cosh⁡(3π/2))1/(2\cosh(3\pi/2)). Their third-to-first amplitude ratio is cosh⁡(π/2)2cosh⁡(3π/2).\boxed{\frac{\cosh(\pi/2)}{2\cosh(3\pi/2)}.} It is smaller than the boundary ratio 1/21/2: the higher spatial frequency loses amplitude faster away from the driven edge.

Step 3: Balance the four edge fluxes. Integrating −∂nu-\partial_nu along each unit edge gives Ftop=−2coth⁡π−coth⁡(3π),Fbottom=2csch⁡π+csch⁡(3π),F_{\mathrm{top}}=-2\coth\pi-\coth(3\pi),\qquad F_{\mathrm{bottom}}=2\,\operatorname{csch}\pi+\operatorname{csch}(3\pi), Fleft=Fright=tanh⁡(π/2)+12tanh⁡(3π/2).F_{\mathrm{left}}=F_{\mathrm{right}} =\tanh(\pi/2)+\tfrac 12\tanh(3\pi/2). For example, at the bottom the outward normal is (0,−1)(0,-1), so the outward heat flux is uyu_y, not −uy-u_y. The identity coth⁡a−csch⁡a=tanh⁡(a/2)\coth a-\operatorname{csch}a=\tanh(a/2) proves that the four totals sum to zero.

Step 4: Interpret the profiles and their sign. Writing s=sin⁡(πx)∈[0,1]s=\sin(\pi x)\in[0,1], the top value is s+12(3s−4s3)=s(5/2−2s2)≥0s+\tfrac 12(3s-4s^3)=s(5/2-2s^2)\ge 0, and it is positive for 0<x<10<x<1. The weak maximum principle gives u≥0u\ge 0; the strong maximum principle gives u>0u>0 inside because the boundary data are not constant. Individual modes need not be positive.

See the diagram in the original worksheet below.

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

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