Heat Equation with Non-Zero Temperature Boundaries — Question 6

PDF ↗

Question 6

A semi-infinite rod has ut=κuxxu_t=\kappa u_{xx} on x>0x>0, with κ,ω>0\kappa,\omega>0 and a maintained oscillating boundary u(0,t)=A+Bcos⁡(ωt)u(0,t)=A+B\cos(\omega t), where A>|B|>0A>|B|>0. Seek its periodic response of the form u(x,t)=A+Be−αxcos⁡(ωt−αx),α>0,u(x,t)=A+B e^{-\alpha x}\cos(\omega t-\alpha x),\qquad \alpha>0, with u(x,t)→Au(x,t)\to A as x→∞x\to\infty. The initial field is to be the trace of this periodic response, so no additional startup transient is requested. Use the rightward flux convention j=−κuxj=-\kappa u_x.

Tasks

  1. Determine α\alpha by direct substitution, and state the initial field that makes the response a solution of the full initial-boundary problem.

  2. Derive the amplitude attenuation and phase delay at distance xx. Find the depth where the oscillation amplitude falls to 1/e1/e of its boundary value, and express it in terms of the period P=2π/ωP=2\pi/\omega.

  3. Compute the boundary flux and its phase relative to the imposed temperature. Explain why constant far-field temperature does not imply zero instantaneous boundary flux.

  4. Prove positivity of the temperature and approach to the far-field value as x→∞x\to\infty, uniformly in time. For A=2,B=1,κ=1,ω=2A=2,B=1,\kappa=1,\omega=2, sketch profiles at phases ωt=0,π/2,π\omega t=0,\pi/2,\pi.

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

Question 6 – Solution

Strategy. A periodic boundary signal diffuses inward with both amplitude attenuation and a spatial phase shift.

Step 1: Match the PDE and the initial trace. Put θ=ωt−αx\theta=\omega t-\alpha x. Then ut=−ωBe−αxsin⁡θu_t=-\omega B e^{-\alpha x}\sin\theta and uxx=−2α2Be−αxsin⁡θu_{xx}=-2\alpha^2B e^{-\alpha x}\sin\theta. Thus α=ω/(2κ)\boxed{\alpha=\sqrt{\omega/(2\kappa)}}. The boundary trace is the imposed cosine, and the compatible initial field is f(x)=A+Be−αxcos⁡(αx)f(x)=A+B e^{-\alpha x}\cos(\alpha x). It is generally not the constant AA; a constant initial field would require an additional transient.

Step 2: Quantify attenuation and delay. The local oscillation amplitude is |B|e−αx|B|e^{-\alpha x} and its phase lag is αx\alpha x. Relative to the boundary waveform, the time delay is αx/ω\alpha x/\omega, understood modulo the period. The penetration depth is d=1/α=2κ/ω=κP/π.\boxed{d=1/\alpha=\sqrt{2\kappa/\omega}=\sqrt{\kappa P/\pi}.} Longer periods penetrate farther, with depth proportional to the square root of the period rather than the period itself.

Step 3: Determine the phase of the boundary flux. Differentiation gives ux=αBe−αx(−cos⁡θ+sin⁡θ)u_x=\alpha B e^{-\alpha x}(-\cos\theta+\sin\theta). Hence j(0,t)=Bκωcos⁡(ωt+π/4).\boxed{j(0,t)=B\sqrt{\kappa\omega}\cos(\omega t+\pi/4).} The flux leads the boundary temperature oscillation by π/4\pi/4 in phase and has zero average over a full period. Instantaneous storage and release near the driven boundary are compatible with a constant far-field limit.

Step 4: Verify bounds and interpret the profiles. For all x≥0,tx\ge 0,t, |u−A|≤|B|e−αx≤|B||u-A|\le|B|e^{-\alpha x}\le|B|, so u≥A−|B|>0u\ge A-|B|>0. Also sup⁡t|u(x,t)−A|≤|B|e−αx→0\sup_t|u(x,t)-A|\le|B|e^{-\alpha x}\to 0 as x→∞x\to\infty. This is a spatial far-field limit uniform in time, not decay of the maintained oscillation at a fixed point as time increases. For the plotted parameters, α=1\alpha=1 and the penetration depth is one.

See the diagram in the original worksheet below.

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

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