Skip to content

Workshop 5: Steady heat conduction

Attempt these problems before opening the worked solutions.

Saturated steam at 130 ∘C130\,{}^\circ\mathrm{C} is flowing inside a steel pipe having an ID of 2.09 cm and an OD of 2.67 cm. The pipe is insulated with 3.8 cm of insulation on the outside. The convective coefficient for the inside steam surface of the pipe is estimated to be hih_i = 5680 W/m2.K\mathrm{W/m^2.K}, while on the outside of the insulation it is estimated to be hoh_o = 11.4 W/m2.K\mathrm{W/m^2.K}. The mean thermal conductivities of the metal and the insulation material are 45 W/m.K\mathrm{W/m.K} and 0.064 W/m.K\mathrm{W/m.K}, respectively. Calculate the heat loss for 1 m of pipe using the method of resistances if the surrounding air is at 26 ∘C26\,{}^\circ\mathrm{C}.

Worked solution

Heat flows radially from the steam, through the steel and insulation, and then into the air. The same heat rate passes through four resistances.

Cross-section of an insulated pipe. The inner and outer steel radii are r1 and r2; insulation extends to r3. Heat crosses the inner film, steel, insulation and outer film in series.
Figure 1: Cross-section of an insulated pipe. The inner and outer steel radii are r1 and r2; insulation extends to r3. Heat crosses the inner film, steel, insulation and outer film in series.

Data

QuantitySymbolValue
Pipe lengthLL1.0 m1.0\ \mathrm m
Inside diameterD1D_10.0209 m0.0209\ \mathrm m
Outside steel diameterD2D_20.0267 m0.0267\ \mathrm m
Radial insulation thicknesstt0.0380 m0.0380\ \mathrm m
Steel conductivityksk_s45 W m−1 K−145\ \mathrm{W\,m^{-1}\,K^{-1}}
Insulation conductivitykik_i0.064 W m−1 K−10.064\ \mathrm{W\,m^{-1}\,K^{-1}}
Steam and air temperaturesTsteam,TairT_{steam},T_{air}130,26 ∘C130,26\,{}^\circ\mathrm C
Inside and outside coefficientshi,hoh_i,h_o5680,11.4 W m−2 K−15680,11.4\ \mathrm{W\,m^{-2}\,K^{-1}}

We need the heat loss from the specified one-metre length.

Assumptions

  • Steady radial heat flow with negligible end losses.
  • Constant conductivities and heat transfer coefficients.
  • Perfect thermal contact between steel and insulation.
  • No separate radiation or fouling resistance.

Calculations

First convert diameters to radii. Insulation thickness is already a radial length, so it is added directly to the outside steel radius:

r1=D12=0.01045 m,r2=D22=0.01335 mr_1=\frac{D_1}{2}=0.01045\ \mathrm m,\qquad r_2=\frac{D_2}{2}=0.01335\ \mathrm mr3=r2+t=0.01335+0.0380=0.05135 mr_3=r_2+t=0.01335+0.0380=0.05135\ \mathrm m

Each film uses its own cylindrical area, 2πrL2\pi rL. Conduction through a cylindrical layer uses both bounding radii. The four resistances are

ElementResistance relationValue, K W−1\mathrm{K\,W^{-1}}
Inside steam filmRci=1/(hi2πr1L)R_{ci}=1/(h_i2\pi r_1L)0.002680.00268
Steel wallRs=ln⁡(r2/r1)/(2πksL)R_s=\ln(r_2/r_1)/(2\pi k_sL)0.0008660.000866
InsulationRi=ln⁡(r3/r2)/(2πkiL)R_i=\ln(r_3/r_2)/(2\pi k_iL)3.353.35
Outside air filmRco=1/(ho2πr3L)R_{co}=1/(h_o2\pi r_3L)0.2720.272

For example, the largest resistance is obtained from

Ri=ln⁡(0.051350.01335)2π(0.064)(1.0)≃3.35 K W−1R_i=\frac{\ln\left(\frac{0.05135}{0.01335}\right)}{2\pi(0.064)(1.0)} \simeq3.35\ \mathrm{K\,W^{-1}}

The elements are in series, so

Rtotal=Rci+Rs+Ri+Rco≃3.626 K W−1R_{total}=R_{ci}+R_s+R_i+R_{co}\simeq3.626\ \mathrm{K\,W^{-1}}

Use the steam-to-air temperature difference across the complete network:

Q˙=Tsteam−TairRtotal=130−263.626≃28.7 W\dot Q=\frac{T_{steam}-T_{air}}{R_{total}} =\frac{130-26}{3.626}\simeq28.7\ \mathrm W

The outer insulation temperature provides a useful direction check:

Ts,o=Tair+Q˙Rco≃26+28.7(0.272)=33.8 ∘CT_{s,o}=T_{air}+\dot QR_{co} \simeq26+28.7(0.272)=33.8\,{}^\circ\mathrm C

It lies above the air temperature and below the steam temperature. The insulation supplies about 92 percent of the total resistance. Its large temperature drop is consistent with the small heat loss.

Answer

The heat loss is approximately 28.7 W28.7\ \mathrm W from one metre of pipe, or 28.7 W m−128.7\ \mathrm{W\,m^{-1}}. This is a heat rate per pipe length, not a heat flux per surface area.

In a laboratory experiment, an electric current is passed through a slab of thickness 2δ2\delta. As a result, heat generation SeS_e takes place in the slab. The slab is insulated on both sides (thickness = t) to satisfy the safety considerations related to the danger stemming from the fact that an electric current is passing. In addition, the insulated slab is cooled on both sides by a fluid of temperature TcT_c, with convective coefficient of h. Determine the temperature at the surface of the insulation, assuming that the thermal conductivities for the slab and the insulation material are constants, denoted by ksk_s and kik_i, respectively.

Worked solution

Choose x=0x=0 at the slab centre and x=δx=\delta at its right face.

Symmetry gives dT/dx=0dT/dx=0 at the centre.

The slab energy balance is

ksd2Tdx2+Se=0.k_s\frac{d^2T}{dx^2}+S_e=0.

Integrating with the symmetry condition gives

dTdx=−Sexks.\frac{dT}{dx}=-\frac{S_ex}{k_s}.

The heat flux leaving either slab face, in W m−2\mathrm{W\,m^{-2}}, is therefore

q′′=−ks[dTdx]δ=Seδ.q''=-k_s\left[\frac{dT}{dx}\right]_{\delta}=S_e\delta.

Each face removes half the generation per unit area of the whole slab.

The insulation carries this constant flux, so

Seδ=h(To−Tc)=kit(Ti−To).S_e\delta=h(T_o-T_c) =\frac{k_i}{t}(T_i-T_o).

Hence the outer surface and interface temperatures are

To=Tc+Seδh,T_o=T_c+\frac{S_e\delta}{h},Ti=Tc+Seδh+Seδtki.T_i=T_c+\frac{S_e\delta}{h} +\frac{S_e\delta t}{k_i}.

Both expressions return kelvin.

The overall generation fixes the surface heat flux, so ksk_s affects the temperature gradient inside the slab rather than ToT_o.

For example the centre temperature is

Ti+Seδ22ks,T_i+\frac{S_e\delta^2}{2k_s},

also in kelvin.

A pipe designed to transport liquid metal is to be embedded in a 1.2 m thick wall. The inside surface of the wall has a constant surface temperature of 925 K, whereas its external surface is exposed to air with a convective heat transfer coefficient of 23 W/m2.K\mathrm{W/m^2.K} and 300 K. The thermal conductivity of the wall varies with temperature according to k=0.073+3.942×10−4Tk = 0.073 + 3.942\times10^{-4}T, where T is temperature in Kelvin, and k is in W/m.K\mathrm{W/m.K}. It is desired that the pipe is located at a point in the wall where the temperature is < 650 K. Calculate how far from the hot surface should the pipe be located and the heat flux for the wall.

Worked solution

Data

QuantitySymbolValue
Wall thicknessLL1.2 m1.2\ \mathrm m
Hot-face temperatureThT_h925 K925\ \mathrm K
Outside airT∞T_\infty300 K300\ \mathrm K
Outside coefficienthh23 W m−2 K−123\ \mathrm{W\,m^{-2}\,K^{-1}}
Conductivityk(T)k(T)0.073+3.942×10−4T W m−1 K−10.073+3.942\times10^{-4}T\ \mathrm{W\,m^{-1}\,K^{-1}}
Pipe temperature limitTlimitT_{limit}650 K650\ \mathrm K

We need the outward heat flux and the depth at which the wall reaches the pipe-temperature limit.

Assumptions

  • Steady one-dimensional conduction with no heat generation.
  • Negligible disturbance of the wall temperature field by the embedded pipe.
  • The given conductivity relation applies throughout the temperature range.

Calculations

Write k(T)=a+bTk(T)=a+bT, where a=0.073 W m−1 K−1a=0.073\,\mathrm{W\,m^{-1}\,K^{-1}} and b=3.942×10−4 W m−1 K−2b=3.942\times10^{-4}\,\mathrm{W\,m^{-1}\,K^{-2}}.

Let xx increase from the hot face towards the air, L=1.2 mL=1.2\,\mathrm{m}, and To=T(L)T_o=T(L).

With no generation the positive heat flux is constant:

q′′=−k(T)dTdx.q''=-k(T)\frac{dT}{dx}.

Integrating across the wall and applying convection gives

q′′L=a(925−To)+b2(9252−To2),q′′=23(To−300).\begin{aligned} q''L&=a(925-T_o) +\frac b2(925^2-T_o^2),\\ q''&=23(T_o-300). \end{aligned}(1)

Substituting the second relation into the first produces

1.971×10−4To2+27.673To−8516.17=0.1.971\times10^{-4}T_o^2+27.673T_o -8516.17=0.

The physical root is To≈307.07 KT_o\approx307.07\,\mathrm{K}.

The other root is negative and lies outside the specified temperature interval.

From Eq. 1, the outward heat flux is

q′′=23(To−300)≈162.6 W m−2.q''=23(T_o-300) \approx162.6\,\mathrm{W\,m^{-2}}.

Integrating from 925 K925\,\mathrm{K} to 650 K650\,\mathrm{K} gives

q′′x650=a(925−650)+b2(9252−6502),q''x_{650}=a(925-650) +\frac b2(925^2-650^2),x650≈0.648 m.x_{650}\approx0.648\,\mathrm{m}.

Since k(T)>0k(T)>0 throughout the wall, the temperature decreases with xx.

The pipe must lie farther than 0.648 m0.648\,\mathrm{m} from the hot face towards the cold face.

Equality gives 650 K650\,\mathrm{K} and does not meet the strict limit.

A finite pipe must have its entire occupied region on the cold side if every part must remain below the limit.

Answer

The outward heat flux is 163 W m−2163\ \mathrm{W\,m^{-2}}. The pipe must be placed more than 0.648 m0.648\ \mathrm m from the hot face towards the cold face to remain below 650 K650\ \mathrm K.