Solving the Heat Equation — Question 8

PDF ↗

Question 8

For α>0\alpha>0, solve the forced heat problem ut=uxx+e−αtsin⁡x,0<x<π,t>0,u(0,t)=u(π,t)=0,u(x,0)=0.u_t=u_{xx}+e^{-\alpha t}\sin x,\qquad 0<x<\pi,\ t>0, \qquad u(0,t)=u(\pi,t)=0,\quad u(x,0)=0. Use a field of the form u(x,t)=qα(t)sin⁡xu(x,t)=q_\alpha(t)\sin x. This is interior forcing; the boundary temperatures remain zero.

Tasks

  1. Derive and solve the scalar initial-value problem for qαq_\alpha, treating α=1\alpha=1 separately.

  2. Verify the PDE and all data, and prove that the constructed amplitude is positive for t>0t>0. Explain why uniqueness holds in the class of smooth solutions continuous into L2L^2 initially.

  3. Determine the unique time of maximum amplitude for α≠1\alpha\ne 1 and its limit as α→1\alpha\to 1. Check the maximum directly for α=1\alpha=1.

  4. Derive the total-heat balance and compare matching temporal decay rates here with unbounded oscillatory resonance. Does the factor tt in the exceptional solution cause unbounded temperature?

Original worksheet page 1: question and worked solution for 9-5-008
Show solutionHide solution

Question 8 – Solution

Strategy. The source has one spatial eigenfunction, so the full problem reduces to a forced scalar decay equation.

Step 1: Solve the amplitude equation including the exceptional rate. Substitution gives qα′+qα=e−αtq_\alpha'+q_\alpha=e^{-\alpha t}, qα(0)=0q_\alpha(0)=0. Multiplying by ete^t and integrating yields qα(t)={(e−αt−e−t)/(1−α),α≠1,te−t,α=1.\boxed{q_\alpha(t)=\begin{cases} (e^{-\alpha t}-e^{-t})/(1-\alpha),&\alpha\ne 1,\\ te^{-t},&\alpha=1. \end{cases}} The exceptional formula follows from the same integrating-factor integral; it is not obtained by substituting into a zero denominator.

Step 2: Verify existence, sign and uniqueness. The scalar ODE verifies the PDE, while the sine factor and qα(0)=0q_\alpha(0)=0 verify endpoints and initial data. The integral representation qα(t)=∫0te−(t−s)e−αsdsq_\alpha(t)=\int_0^t e^{-(t-s)}e^{-\alpha s}\,ds is strictly positive for t>0t>0. The difference of two solutions with the same source solves the homogeneous zero-data problem. Its squared L2L^2 norm is nonincreasing by integration by parts; starting at a positive time and using the zero initial norm proves uniqueness in the stated class.

Step 3: Locate the unique maximum. For α≠1\alpha\ne 1, qα′=(e−t−αe−αt)/(1−α)q_\alpha'=(e^{-t}-\alpha e^{-\alpha t})/(1-\alpha). The equation qα′=0q_\alpha'=0 is e(α−1)t=αe^{(\alpha-1)t}=\alpha, giving tα=log⁡αα−1>0.\boxed{t_\alpha=\frac{\log\alpha}{\alpha-1}>0.} There is exactly one root; qα′(0)=1>0q_\alpha'(0)=1>0 and the derivative is negative for sufficiently large time, proving it is the unique maximum. As α→1\alpha\to 1, tα→1t_\alpha\to 1. Directly (te−t)′=e−t(1−t)(te^{-t})'=e^{-t}(1-t), with maximum 1/e1/e at t=1t=1.

Step 4: Interpret the forced heat balance. The total heat is H=2qαH=2q_\alpha. Since [ux]0π=−2qα[u_x]_0^\pi=-2q_\alpha and the source integral is 2e−αt2e^{-\alpha t}, H′=−2qα+2e−αt,H'=-2q_\alpha+2e^{-\alpha t}, which agrees with the scalar ODE. At α=1\alpha=1, matching the source decay rate to the homogeneous decay produces te−tte^{-t}. This field is bounded and tends to zero; its time factor does not imply unbounded oscillatory resonance. For every fixed α>0\alpha>0, the source and the solution ultimately decay.

Original worksheet page 2: question and worked solution for 9-5-008

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