Heat Equation with Non-Zero Temperature Boundaries — Question 5

PDF ↗

Question 5

Let both endpoint temperatures rise at a constant rate r>0r>0: ut=κuxx,0<x<L,u(0,t)=A+rt,u(L,t)=B+rt,u_t=\kappa u_{xx},\quad 0<x<L,\qquad u(0,t)=A+rt,\quad u(L,t)=B+rt, where κ,L>0\kappa,L>0 and u(x,0)=ℓ(x):=A+(B−A)x/Lu(x,0)=\ell(x):=A+(B-A)x/L. You may use the sine expansion of x(L−x)x(L-x).

Tasks

  1. Subtract the instantaneous boundary line ℓ(x)+rt\ell(x)+rt and derive the PDE and initial data for the remaining field. Explain why that line alone is not a solution.

  2. Find a time-independent zero-endpoint correction P(x)P(x) so that ℓ+rt+P\ell+rt+P solves the PDE. Add the transient needed to recover the original initial field.

  3. Verify the reconstructed solution, justify the initial trace and positive-time differentiations, and determine its limiting offset relative to the moving boundary line.

  4. Compute the limiting midpoint lag, the limiting mean lag and the eventual total-heat growth rate. Explain why there is no stationary equilibrium even though the offset approaches a fixed shape.

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

Question 5 – Solution

Strategy. A boundary lifting with nonzero time derivative creates a source in the transformed equation; a persistent lag balances that source.

Step 1: Keep the source created by the lifting. For w=u−ℓ−rtw=u-\ell-rt, wt=κwxx−r,w(0,t)=w(L,t)=0,w(x,0)=0.w_t=\kappa w_{xx}-r,\qquad w(0,t)=w(L,t)=0,\qquad w(x,0)=0. The line ℓ+rt\ell+rt has time derivative rr and zero second spatial derivative, so it fails the original source-free PDE when r>0r>0.

Step 2: Find the lag profile and the initial correction. A fixed PP must satisfy κP″=r\kappa P''=r, P(0)=P(L)=0P(0)=P(L)=0, giving P=−rx(L−x)/(2κ)P=-r x(L-x)/(2\kappa). The remaining transient has initial data −P-P. Its odd sine coefficients are 4rL2/(κπ3n3)4rL^2/(\kappa\pi^3n^3). Thus u=ℓ+rt+P+4rL2κπ3∑n odde−κ(nπ/L)2tsin⁡(nπx/L)n3.\boxed{u=\ell+rt+P+\frac{4rL^2}{\kappa\pi^3} \sum_{n\text{ odd}}\frac{e^{-\kappa(n\pi/L)^2t}\sin(n\pi x/L)}{n^3}.}

Step 3: Verify all data and the moving-frame limit. The field ℓ+rt+P\ell+rt+P satisfies ut=κuxx=ru_t=\kappa u_{xx}=r; its added transient satisfies the homogeneous heat equation. All correction terms vanish at the endpoints. At t=0t=0, the summable sine expansion equals −P-P uniformly and recovers ℓ\ell. For t≥τ>0t\ge\tau>0, the exponential transient permits all termwise derivatives. Initial corner smoothness is not implied: the boundary time derivative is rr, whereas κℓ″=0\kappa\ell''=0. As t→∞t\to\infty, the transient tends uniformly to zero, so u−ℓ−rt→Pu-\ell-rt\to P uniformly.

Step 4: Quantify the persistent lag and heat input. The midpoint offset is P(L/2)=−rL2/(8κ)P(L/2)=-rL^2/(8\kappa), and the mean offset is L−1∫P=−rL2/(12κ)L^{-1}\int P=-rL^2/(12\kappa). Their negatives are the positive lags behind the boundary line. Since P′(L)−P′(0)=rL/κP'(L)-P'(0)=rL/\kappa, the eventual flux balance is H′=κ[ux]0L→rLH'=\kappa[u_x]_0^L\to rL, also obtained by differentiating the integrated solution. Boundary values grow indefinitely; it is the offset, not the absolute temperature, that approaches a fixed state.

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

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