The Heat Equation — Question 7

PDF ↗

Question 7

A constant-area rod with constant k,C>0k,C>0 is prescribed the candidate field u(x,t)=Ta+at+bx(L−x),a,b,L>0,0≤x≤L.u(x,t)=T_a+at+b\,x(L-x),\qquad a,b,L>0,\quad 0\leq x\leq L. The model permits a volumetric source ss and obeys Cut=kuxx+sC u_t=k u_{xx}+s. Its lateral surface is insulated; end conditions must be determined consistently.

Tasks

  1. Determine the unique source field that makes the candidate satisfy the equation.

  2. State the initial and Dirichlet boundary data realized by the candidate. Decide whether it also satisfies insulated-end conditions.

  3. Compute both outward end heat fluxes and verify the integrated energy balance, using cross-sectional area AA.

  4. Suppose the source is removed while the same formula is retained but aa is allowed to be real. Find the required aa and explain why changing only the boundary values cannot fix a nonzero interior residual.

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

Question 7 – Solution

Strategy. A manufactured temperature field determines its source and boundary data; those data cannot be prescribed independently afterward.

Step 1: Compute the interior residual. The derivatives are ut=au_t=a, ux=b(L−2x)u_x=b(L-2x) and uxx=−2bu_{xx}=-2b. Substitution requires s=Ca+2kb.\boxed{s=Ca+2kb.} This is a positive spatially and temporally constant source, in W/m3\mathrm{W/m^3}. It supplies both the rising stored energy and heat lost at the ends.

Step 2: State the compatible data. At t=0t=0, u(x,0)=Ta+bx(L−x)u(x,0)=T_a+bx(L-x). Both endpoint temperatures equal Ta+atT_a+at. The endpoint derivatives are bLbL and −bL-bL, neither zero, so this field does not satisfy insulated-end conditions when b>0b>0. Prescribing insulation in addition to the realized end temperatures would contradict the candidate rather than add useful information.

Step 3: Verify total energy accounting. The outward fluxes are kux(0)=kbLku_x(0)=kbL and −kux(L)=kbL-ku_x(L)=kbL. The stored heat relative to TaT_a is E=CA∫0L(u−Ta)dx=CA(atL+bL36),E′=CAaL.E=CA\int_0^L(u-T_a)\,dx =CA\left(atL+\frac{bL^3}{6}\right),\qquad E'=CAaL. Total source power minus total outward power is AL(Ca+2kb)−2AkbL=CAaL=E′.\boxed{AL(Ca+2kb)-2AkbL=CAaL=E'.} This independently checks the source sign and both boundary signs.

Step 4: Remove the source consistently. For s=0s=0, the formula can satisfy the PDE only if a=−2kbC.\boxed{a=-\frac{2kb}{C}.} Its initially warm curved profile then decreases at a constant rate while both endpoint temperatures follow the same decreasing law. If aa retains a different value, the residual Ca+2kbCa+2kb is nonzero at every interior point. Boundary values alone cannot repair that local equation while keeping the same candidate field. Changing boundary data generally requires changing the field itself.

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

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