Laplace's Equation — Question 6

PDF ↗

Question 6

A conducting annulus a<r<Ra<r<R has constant inner and outer temperatures Ti,ToT_i,T_o, where 0<a<R0<a<R and Ti≠ToT_i\ne T_o. Conductivity is one. Seek a radial harmonic field. Total outward heat flux means the integral of −∂nu-\partial_nu over the indicated boundary circle.

Tasks

  1. Derive the temperature and prove it is the unique continuous harmonic solution of these boundary data, including among nonradial competitors.

  2. Compute both total boundary fluxes with the annulus’s outward normals. Identify its conductance, defined as the magnitude of transmitted total heat flux divided by |Ti−To||T_i-T_o|.

  3. Compute the Dirichlet energy and describe its relation to the temperature difference and conductance.

  4. Hold R,Ti,ToR,T_i,T_o fixed and let a↓0a\downarrow 0. Determine the local limit of temperature away from the origin, the energy limit and the largest gradient. Explain the failure of uniform convergence to ToT_o. Sketch the profiles for R=1,Ti=1,To=0R=1,T_i=1,T_o=0 and a=0.1,0.3,0.5a=0.1,0.3,0.5.

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

Question 6 – Solution

Strategy. The logarithmic radial mode is admissible on an annulus, but its small-hole limit is not uniform near the shrinking inner boundary.

Step 1: Solve the radial equation. The equation (rur)′=0(ru_r)'=0 gives u=Alog⁡r+Bu=A\log r+B. Writing D=log⁡(R/a)>0D=\log(R/a)>0 and ΔT=Ti−To\Delta T=T_i-T_o, the boundary data give u(r)=To+ΔTlog⁡(R/r)D.\boxed{u(r)=T_o+\Delta T\,\frac{\log(R/r)}{D}.} Its radial Laplacian is zero and its two traces are correct. Any other continuous harmonic solution differs by a harmonic function with zero boundary data; the maximum principle on the annulus proves uniqueness.

Step 2: Respect the inner normal. Here ur=−ΔT/(rD)u_r=-\Delta T/(rD). At r=Rr=R the outward normal points toward increasing rr, whereas at r=ar=a it points toward decreasing rr. Thus Fouter=2πΔTD,Finner=−2πΔTD,G=2πD.\boxed{F_{\mathrm{outer}}=\frac{2\pi\Delta T}{D},\qquad F_{\mathrm{inner}}=-\frac{2\pi\Delta T}{D},\qquad G=\frac{2\pi}{D}.} The sum is zero. When ΔT>0\Delta T>0, heat enters the annulus at the inner boundary and leaves at the outer boundary.

Step 3: Integrate the energy. Direct radial integration yields ℰ=2π∫aR(ΔT)2r2D2rdr=2π(ΔT)2D=G(ΔT)2.\mathcal E=2\pi\int_a^R\frac{(\Delta T)^2}{r^2D^2}\,r\,dr =\boxed{\frac{2\pi(\Delta T)^2}{D}=G(\Delta T)^2}. Equivalently, Green’s identity gives ℰ=−TiFinner−ToFouter\mathcal E=-T_iF_{\mathrm{inner}}-T_oF_{\mathrm{outer}}, with the same sign and value.

Step 4: Separate local and uniform limits. As a↓0a\downarrow 0, D→∞D\to\infty. On each fixed r0≤r≤Rr_0\le r\le R the field tends uniformly to ToT_o, and ℰ→0\mathcal E\to 0. However sup⁡a≤r≤R|u−To|=|ΔT|\sup_{a\le r\le R}|u-T_o|=|\Delta T| for every aa. Also max⁡|∇u|=|ΔT|/[alog⁡(R/a)]→∞\max|\nabla u|=|\Delta T|/[a\log(R/a)]\to\infty at the inner circle. The concentration near a shrinking boundary is compatible with vanishing integrated energy.

See the diagram in the original worksheet below.

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

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