Laplace's Equation — Question 7

PDF ↗

Question 7

On the unit disk consider the Neumann problem Δu=0,∂nu=1+3cos⁡θ+4sin⁡(2θ)on r=1.\Delta u=0,\qquad \partial_nu=1+3\cos\theta+4\sin(2\theta) \quad\text{on }r=1. Here ∂nu\partial_nu is the outward normal derivative, not the outward conductive heat flux.

Tasks

  1. Decide whether a classical solution exists. Derive the necessary compatibility condition and evaluate it for these data.

  2. Subtract a constant cc from the prescribed derivative to make the data compatible. Find cc and construct a regular harmonic solution whose area average on the disk is zero.

  3. Prove uniqueness under that normalization, and state precisely the nonuniqueness without it.

  4. For the corrected solution, compute its boundary trace, Dirichlet energy and net conductive heat flux. Check the energy independently by integrating the Cartesian gradient.

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

Question 7 – Solution

Strategy. Test flux compatibility before separation; then divide each nonconstant disk mode by its normal-derivative multiplier.

Step 1: Test the solvability condition. The divergence theorem gives 0=∫diskΔudA=∫∂disk∂nuds.0=\int_{\mathrm{disk}}\Delta u\,dA =\int_{\partial{\mathrm{disk}}}\partial_nu\,ds. The prescribed derivative integrates to 2π2\pi, since both oscillatory terms integrate to zero. It violates a necessary condition, so the original problem has no classical solution.

Step 2: Correct the mean and construct the field. The corrected integral is 2π(1−c)2\pi(1-c), so the only constant correction is c=1c=1. Since ∂r(rn)=n\partial_r(r^n)=n at r=1r=1, the regular solution is u=3rcos⁡θ+2r2sin⁡(2θ)+C=3x+4xy+C.u=3r\cos\theta+2r^2\sin(2\theta)+C=3x+4xy+C. The nonconstant terms have zero area average by symmetry. Thus normalization gives C=0C=0, and u(x,y)=3x+4xy.\boxed{u(x,y)=3x+4xy.} Its Cartesian Laplacian vanishes, and its outward derivative is exactly 3cos⁡θ+4sin⁡(2θ)3\cos\theta+4\sin(2\theta).

Step 3: Establish uniqueness up to a constant. For two solutions, their difference ww has Δw=0\Delta w=0 and ∂nw=0\partial_nw=0. Green’s identity gives ∫|∇w|2=∫∂diskw∂nw=0\int|\nabla w|^2=\int_{\partial{\mathrm{disk}}}w\partial_nw=0. Hence ∇w=0\nabla w=0, so ww is constant on the connected disk. The zero-area-mean condition forces that constant to vanish. Without a normalization, every constant translate is another solution and these are all possibilities.

Step 4: Check energy in two coordinate systems. The boundary trace is 3cos⁡θ+2sin⁡(2θ)3\cos\theta+2\sin(2\theta). Orthogonality gives ℰ=∫02π[3cos⁡θ+2sin⁡(2θ)][3cos⁡θ+4sin⁡(2θ)]dθ=17π.\mathcal E=\int_0^{2\pi} [3\cos\theta+2\sin(2\theta)][3\cos\theta+4\sin(2\theta)]\,d\theta =\boxed{17\pi}. Independently, |∇u|2=(3+4y)2+16x2=9+24y+16r2|\nabla u|^2=(3+4y)^2+16x^2=9+24y+16r^2. Its disk integral is 9π+0+16(π/2)=17π9\pi+0+16(\pi/2)=17\pi. The corrected normal derivative integrates to zero, so its negative, the conductive flux, has zero total as required.

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

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