The Heat Equation — Question 3

PDF ↗

Question 3

A rod has constant area, insulated sides, no internal source and conductivity k(x)=k0(1+x/L)k(x)=k_0(1+x/L), with k0,L>0k_0,L>0. Its end temperatures are u(0)=Thu(0)=T_h, u(L)=Tcu(L)=T_c, with Th>TcT_h>T_c. The steady equation is (k(x)ux)x=0(k(x)u_x)_x=0, and flux per area toward the right is q=−k(x)uxq=-k(x)u_x.

Tasks

  1. Derive the exact temperature and heat flux.

  2. Verify the equation and both boundary values, and locate the largest magnitude of the temperature gradient.

  3. Test the straight-line interpolant between the same endpoints in the correct variable-conductivity equation. Explain its physical failure.

  4. Define keffk_{\mathrm{eff}} by q=keff(Th−Tc)/Lq=k_{\mathrm{eff}}(T_h-T_c)/L. Compute it and compare it with the arithmetic spatial average of kk. Sketch the normalized exact and linear temperature profiles.

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

Question 3 – Solution

Strategy. Constant steady heat flux requires the gradient to vary inversely with conductivity.

Step 1: Integrate the flux law. Since (kux)x=0(ku_x)_x=0, the flux qq is constant. Integrating ux=−q/k(x)u_x=-q/k(x) gives Tc−Th=−qLk0log⁡2.T_c-T_h=-\frac{qL}{k_0}\log 2. Consequently q=k0(Th−Tc)Llog⁡2,u(x)=Th−(Th−Tc)log⁡(1+x/L)log⁡2.\boxed{q=\frac{k_0(T_h-T_c)}{L\log 2},\qquad u(x)=T_h-(T_h-T_c)\frac{\log(1+x/L)}{\log 2}.} The flux is positive: heat travels from the hotter left end toward the right.

Step 2: Check the profile and gradient. The logarithm is zero at 00 and log⁡2\log 2 at LL, giving the two endpoint values. Also ux=−Th−TcLlog⁡2(1+x/L),k(x)ux=−k0(Th−Tc)Llog⁡2.u_x=-\frac{T_h-T_c}{L\log 2\,(1+x/L)},\qquad k(x)u_x=-\frac{k_0(T_h-T_c)}{L\log 2}. Thus the product is constant as required. The gradient magnitude is largest at x=0x=0 and decreases by a factor of two from left to right.

Step 3: Diagnose the linear candidate. For ℓ=Th−(Th−Tc)x/L\ell=T_h-(T_h-T_c)x/L, one has ℓxx=0\ell_{xx}=0, but (kℓx)x=−k0(Th−Tc)L2≠0.\boxed{(k\ell_x)_x=-\frac{k_0(T_h-T_c)}{L^2}\ne 0.} Its rightward flux grows along the rod, so more heat leaves each small interval than enters. Without a source this profile would cool locally, not remain steady. Using kuxx=0ku_{xx}=0 discards the essential term k′uxk'u_x.

Step 4: Interpret the effective conductivity. The definition gives keff=k0/log⁡2k_{\mathrm{eff}}=k_0/\log 2, while the arithmetic average is 3k0/23k_0/2. In fact log⁡2>2/3\log 2>2/3: the function 1/t1/t is strictly convex and its average on [1,2][1,2] exceeds its midpoint value 2/32/3. Hence keff=k0log⁡2<3k02.\boxed{k_{\mathrm{eff}}=\frac{k_0}{\log 2}<\frac{3k_0}{2}.} Steady conduction uses the reciprocal average keff=L/∫0Lk−1dxk_{\mathrm{eff}}=L/\int_0^L k^{-1}\,dx. The graph uses ξ=x/L\xi=x/L and θ=(u−Tc)/(Th−Tc)\theta=(u-T_c)/(T_h-T_c).

See the diagram in the original worksheet below.

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

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