Heat Equation with Non-Zero Temperature Boundaries — Question 7

PDF ↗

Question 7

On 0<x<10<x<1, let ut=uxxu_t=u_{xx} and u(x,0)=Au(x,0)=A. The right endpoint stays at AA. The left is held at A+BA+B for 0<t<τ0<t<\tau and returns to AA for t≥τt\ge\tau, where B,τ>0B,\tau>0. You may use the unit boundary-step response S(x,t)=1−x−2π∑n≥1e−n2π2tsin⁡(nπx)n,t>0,S(x,t)=1-x-\frac 2\pi\sum_{n\ge 1}\frac{e^{-n^2\pi^2t}\sin(n\pi x)}n, \qquad t>0, which has zero initial interior trace, left boundary one and right boundary zero. Set S(x,0)=0S(x,0)=0 for 0<x<10<x<1 when discussing the switching time.

Tasks

  1. Construct the solution before and after the switch by superposing boundary step responses. Verify the boundary traces on both open time intervals.

  2. Derive a homogeneous-boundary sine series for t>τt>\tau, with coefficients at time τ\tau. Explain continuity at the switch for interior points and why continuity at the switching boundary corner is impossible.

  3. Compute the total excess heat after the switch and prove that it strictly decreases. Determine the long-time equilibrium.

  4. For A=B=1A=B=1, τ=0.1\tau=0.1, sketch profiles at t=0.03,0.08,0.13,0.3t=0.03,0.08,0.13,0.3. Explain why the rod does not return instantly to its initial field when the boundary pulse ends.

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

Question 7 – Solution

Strategy. Switching off a boundary step is a delayed negative step; its instantaneous boundary effect does not erase the stored interior heat.

Step 1: Superpose the two boundary events. Linearity gives u=A+BS(x,t)(0<t<τ),u=A+B[S(x,t)−S(x,t−τ)](t>τ).\boxed{u=A+B S(x,t)\quad(0<t<\tau),\qquad u=A+B[S(x,t)-S(x,t-\tau)]\quad(t>\tau).} Each term solves the heat equation in its time interval. Before the switch, the left trace is A+BA+B; afterward the two unit left traces cancel and leave AA. Both right traces are always AA. The initial interior trace is AA.

Step 2: Identify the post-switch initial data. The stationary lines cancel after the switch, leaving u−A=2Bπ∑n≥11−e−n2π2τne−n2π2(t−τ)sin⁡(nπx).u-A=\frac{2B}{\pi}\sum_{n\ge 1} \frac{1-e^{-n^2\pi^2\tau}}n e^{-n^2\pi^2(t-\tau)}\sin(n\pi x). These are the sine coefficients of BS(x,τ)B S(x,\tau) on the open interval. For every interior xx, S(x,t−τ)→0S(x,t-\tau)\to 0 as t↓τt\downarrow\tau, so the interior trace agrees with the pre-switch value. The left boundary itself jumps from A+BA+B to AA; no single continuous trace at that corner can match both. For times strictly after the switch, all differentiated series converge uniformly by Gaussian decay. The switching trace is also recovered in L2L^2.

Step 3: Track the retained heat. For t>τt>\tau, integration yields Hex(t)=4Bπ2∑n odd1−e−n2π2τn2e−n2π2(t−τ).H_{\mathrm{ex}}(t)=\frac{4B}{\pi^2}\sum_{n\text{ odd}} \frac{1-e^{-n^2\pi^2\tau}}{n^2}e^{-n^2\pi^2(t-\tau)}. Its derivative is −4B∑n odd(1−e−n2π2τ)e−n2π2(t−τ)<0-4B\sum_{n\text{ odd}}(1-e^{-n^2\pi^2\tau}) e^{-n^2\pi^2(t-\tau)}<0. This equals [ux]01[u_x]_0^1 after the switch. At any later fixed starting time the exponentially damped coefficients are summable, which proves uniform convergence to the equilibrium u=Au=A as t→∞t\to\infty.

Step 4: Interpret the finite pulse. The profiles before the switch have left value two; those after it have left value one. The right value is one throughout. At the switch the excess heat is B∫01S(x,τ)dx>0B\int_0^1 S(x,\tau)\,dx>0, as follows from the positive derivative of step-response heat and its zero initial value. It must subsequently diffuse out. A changed boundary value is not an instantaneous reset of the interior.

See the diagram in the original worksheet below.

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

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