The Heat Equation — Question 9

PDF ↗

Question 9

Two homogeneous rods have lengths LL and 2L2L, the same k,C>0k,C>0, insulated sides and no sources. At both ends each rod exchanges heat with the same ambient TaT_a by an outward convection law. The first rod has coefficient h>0h>0. The two initial profiles have the same shape as functions of relative position. Set α=k/C\alpha=k/C and choose a nonzero temperature scale ΔT\Delta T.

Tasks

  1. Nondimensionalize the first problem using ξ=x/L\xi=x/L, η=αt/L2\eta=\alpha t/L^2 and θ=(u−Ta)/ΔT\theta=(u-T_a)/\Delta T. Identify the dimensionless boundary parameter.

  2. Verify the PDE and initial shape for the proposed rescaling u2(x,t)=u1(x/2,t/4)u_2(x,t)=u_1(x/2,t/4) on the longer rod.

  3. Derive exactly which convection coefficient h2h_2 makes this rescaling satisfy both boundary laws.

  4. Assess the claim that doubling length always quadruples every thermal response time when the same hh is used. State when the proposed scaling is valid, including the perfectly insulated case.

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

Question 9 – Solution

Strategy. Match the dimensionless boundary conditions as well as the interior diffusion time scale.

Step 1: Identify the dimensionless problem. The chain rule gives θη=θξξ\theta_\eta=\theta_{\xi\xi}. The outward convection laws become θξ(0,η)=Biθ(0,η),−θξ(1,η)=Biθ(1,η),Bi=hLk.\boxed{\theta_\xi(0,\eta)=\mathrm{Bi}\,\theta(0,\eta),\qquad -\theta_\xi(1,\eta)=\mathrm{Bi}\,\theta(1,\eta),\quad \mathrm{Bi}=\frac{hL}{k}.} This dimensionless coefficient is the Biot number for the chosen length scale. The interior time scale is L2/αL^2/\alpha, but the boundary parameter also affects the dimensionless response.

Step 2: Verify the rescaled interior field. For u2(x,t)=u1(x/2,t/4)u_2(x,t)=u_1(x/2,t/4), (u2)t=14(u1)t,(u2)xx=14(u1)xx.(u_2)_t=\tfrac 14(u_1)_t,\qquad (u_2)_{xx}=\tfrac 14(u_1)_{xx}. Thus the same diffusivity equation holds. At t=0t=0, relative positions agree: x/(2L)=(x/2)/Lx/(2L)=(x/2)/L, so the stated initial-shape requirement is met.

Step 3: Match both convection conditions. At the left endpoint, k(u2)x(0,t)=12k(u1)x(0,t/4)=h2(u2(0,t)−Ta).k(u_2)_x(0,t)=\tfrac 12 k(u_1)_x(0,t/4) =\frac h2\bigl(u_2(0,t)-T_a\bigr). The right outward derivative scales by the same factor. Therefore the coefficient making this a rescaling of the full problem is h2=h/2,h2(2L)/k=hL/k.\boxed{h_2=h/2,\qquad h_2(2L)/k=hL/k.} Both rods then have the same dimensionless boundary problem.

Step 4: Qualify the length-squared claim. If h2=hh_2=h, the longer rod instead has Biot number 2hL/k2hL/k. The proposed field generally fails the boundary law; the full response is not obtained merely by multiplying time by four. A special zero-excess equilibrium cannot establish a general scaling rule. The factor-four rescaling is valid when convection is also adjusted to h/2h/2 and the initial shapes match. If both rods are perfectly insulated (h=0h=0), the zero-flux conditions are preserved automatically. Interior dimensional analysis alone does not justify ignoring changed boundary parameters.

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

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