Steady State Heat Transfer Notes

Steady State Heat Transfer

Conduction

  • Conduction: Heat transfer through a solid material due to a temperature gradient.

  • Thermal resistance circuit analogy.

  • Fourier's Law of Conduction: qxA=kdTdx\frac{q_x}{A} = -k \frac{dT}{dx} , where:

    • qxq_x is the heat transfer rate in the x-direction.

    • AA is the cross-sectional area perpendicular to the heat flow.

    • kk is the thermal conductivity of the material.

    • dTdx\frac{dT}{dx} is the temperature gradient in the x-direction.

  • Assumptions:

    • Steady-state conditions (temperature does not change with time).

    • One-dimensional heat transfer.

    • Constant thermal conductivity.

  • Boundary Conditions:

    • At x=x<em>1x = x<em>1, T=T</em>1T = T</em>1

    • At x=x<em>2x = x<em>2, T=T</em>2T = T</em>2

  • Integrating Fourier's Law:

    • <em>x</em>1x<em>2dx=kAq</em>x<em>T</em>1T2dT\int<em>{x</em>1}^{x<em>2} dx = -\frac{kA}{q</em>x} \int<em>{T</em>1}^{T_2} dT

    • q<em>x(x</em>2x<em>1)=kA(T</em>2T1)q<em>x (x</em>2 - x<em>1) = -kA (T</em>2 - T_1)

  • Temperature Distribution:

    • T(x)=T<em>1q</em>xkA(xx1)T(x) = T<em>1 - \frac{q</em>x}{kA} (x - x_1)

  • Thermal Resistance Concept:

    • R<em>t=T</em>1T<em>2q</em>x=x<em>2x</em>1kAR<em>t = \frac{T</em>1 - T<em>2}{q</em>x} = \frac{x<em>2 - x</em>1}{kA}

Example

  • Given a material with a temperature of 110°C on one side and 90°C on the other side, with a thickness of 1 cm, determine the temperature at 0.5 cm from the 110°C surface.

  • Calculations:

    • R<em>t=x</em>2x1kA=0.01m17W/(mK)1m2=5.88×104°C/WR<em>t = \frac{x</em>2 - x_1}{kA} = \frac{0.01 \, m}{17 \, W/(m \cdot K) \cdot 1 \, m^2} = 5.88 \times 10^{-4} \, °C/W

    • q<em>xA=T</em>1T<em>2R</em>t=110905.88×104=34013W/m2\frac{q<em>x}{A} = \frac{T</em>1 - T<em>2}{R</em>t} = \frac{110 - 90}{5.88 \times 10^{-4}} = 34013 \, W/m^2

    • Rt,0.5=0.005m17W/(mK)1m2=2.94×104°C/WR_{t,0.5} = \frac{0.005 \, m}{17 \, W/(m \cdot K) \cdot 1 \, m^2} = 2.94 \times 10^{-4} \, °C/W

    • T=T<em>1q</em>xARt,0.5=11034013×2.94×104=100°CT = T<em>1 - \frac{q</em>x}{A} R_{t,0.5} = 110 - 34013 \times 2.94 \times 10^{-4} = 100 \, °C

Cylindrical Coordinates

  • Heat transfer through a cylindrical geometry.

  • Fourier's Law in Cylindrical Coordinates: qr=kAdTdr=k(2πrL)dTdrq_r = -kA \frac{dT}{dr} = -k (2 \pi r L) \frac{dT}{dr}, where:

    • rr is the radial distance.

    • LL is the length of the cylinder.

  • Boundary Conditions:

    • At r=r<em>ir = r<em>i, T=T</em>iT = T</em>i

    • At r=r<em>or = r<em>o, T=T</em>oT = T</em>o

  • Temperature Distribution:

    • <em>r</em>ir<em>odrr=2πkLq</em>r<em>T</em>iTodT\int<em>{r</em>i}^{r<em>o} \frac{dr}{r} = -\frac{2 \pi k L}{q</em>r} \int<em>{T</em>i}^{T_o} dT

    • ln(r<em>or</em>i)=2πkLq<em>r(T</em>iTo)\ln \left( \frac{r<em>o}{r</em>i} \right) = \frac{2 \pi k L}{q<em>r} (T</em>i - T_o)

  • Thermal Resistance:

    • R<em>t=T</em>iT<em>oq</em>r=ln(r<em>o/r</em>i)2πkLR<em>t = \frac{T</em>i - T<em>o}{q</em>r} = \frac{\ln(r<em>o/r</em>i)}{2 \pi k L}

Example

  • Given a pipe with inner radius 0.03 m, outer radius 0.05 m, length 40 m, and thermal conductivity 43 W/m°C, with inner temperature 115°C and outer temperature 90°C, calculate the conductive heat loss.

  • Calculations:

    • R<em>t=ln(r</em>o/ri)2πkL=ln(0.05/0.03)2π×43×40=4.727×105°C/WR<em>t = \frac{\ln(r</em>o/r_i)}{2 \pi k L} = \frac{\ln(0.05/0.03)}{2 \pi \times 43 \times 40} = 4.727 \times 10^{-5} \, °C/W

    • q<em>r=T</em>iT<em>oR</em>t=115904.727×105=528903Wq<em>r = \frac{T</em>i - T<em>o}{R</em>t} = \frac{115 - 90}{4.727 \times 10^{-5}} = 528903 \, W

Multilayered Systems

  • Heat transfer through multiple layers of different materials.

  • Consider 3 layers: B, C and D having thickness δx<em>B\delta x<em>B, δx</em>C\delta x</em>C and δx<em>D\delta x<em>D and thermal conductivities K</em>BK</em>B, K<em>CK<em>C and K</em>DK</em>D respectively.

  • The total temperature difference is the sum of the temperature drops across each layer: ΔT=ΔT<em>B+ΔT</em>C+ΔTD\Delta T = \Delta T<em>B + \Delta T</em>C + \Delta T_D

  • Heat transfer through each layer:

    • q=K<em>BAΔT</em>BΔx<em>B=K</em>CAΔT<em>CΔx</em>C=K<em>DAΔT</em>DΔxDq = -K<em>B A \frac{\Delta T</em>B}{\Delta x<em>B} = -K</em>C A \frac{\Delta T<em>C}{\Delta x</em>C} = -K<em>D A \frac{\Delta T</em>D}{\Delta x_D}, where A is the area.

  • Total thermal resistance is the sum of individual resistances:

    • R<em>total=R</em>B+R<em>C+R</em>DR<em>{total} = R</em>B + R<em>C + R</em>D

  • Heat flux is calculated as follows:

    • q=T<em>1T</em>4δx<em>BK</em>BA+δx<em>CK</em>CA+δx<em>DK</em>DAq = \frac{T<em>1 - T</em>4}{\frac{\delta x<em>B}{K</em>B A} + \frac{\delta x<em>C}{K</em>C A} + \frac{\delta x<em>D}{K</em>D A}}

  • Two-Layer System:

    • R<em>t=R</em>tA+RtBR<em>t = R</em>{tA} + R_{tB}

    • Cylindrical: R<em>t=ln(r</em>2/r<em>1)2πLK</em>A+ln(r<em>3/r</em>2)2πLKBR<em>t = \frac{\ln(r</em>2/r<em>1)}{2 \pi L K</em>A} + \frac{\ln(r<em>3/r</em>2)}{2 \pi L K_B}

Example

  • Calculate the interface temperature between steel and insulation.

  • Given: Steel pipe with k<em>1=17W/(mK)k<em>1 = 17 \, W/(m \cdot K), insulation with k</em>2=0.035W/(mK)k</em>2 = 0.035 \, W/(m \cdot K), L=1mL = 1 \, m, T<em>=25°CT<em>\infty = 25 \, °C, T</em>1=130°CT</em>1 = 130 \, °C, r<em>1=0.04mr<em>1 = 0.04 \, m, r</em>2=0.06mr</em>2 = 0.06 \, m, r3=0.1mr_3 = 0.1 \, m

  • Calculations:

    • R<em>t1=ln(r</em>2/r<em>1)2πLk</em>1=ln(0.06/0.04)2π×1×17=0.0038°C/WR<em>{t1} = \frac{\ln(r</em>2/r<em>1)}{2 \pi L k</em>1} = \frac{\ln(0.06/0.04)}{2 \pi \times 1 \times 17} = 0.0038 \, °C/W

    • R<em>t2=ln(r</em>3/r<em>2)2πLk</em>2=ln(0.1/0.06)2π×1×0.035=2.3229°C/WR<em>{t2} = \frac{\ln(r</em>3/r<em>2)}{2 \pi L k</em>2} = \frac{\ln(0.1/0.06)}{2 \pi \times 1 \times 0.035} = 2.3229 \, °C/W

    • q=T<em>1T</em>R<em>t1+R</em>t2=130250.0038+2.3229=45.13Wq = \frac{T<em>1 - T</em>\infty}{R<em>{t1} + R</em>{t2}} = \frac{130 - 25}{0.0038 + 2.3229} = 45.13 \, W

    • T<em>i=T</em>1qRt1=13045.13×0.0038=129.83°CT<em>i = T</em>1 - q R_{t1} = 130 - 45.13 \times 0.0038 = 129.83 \, °C

Convection

  • Convection: Heat transfer between a surface and a moving fluid.

  • Convective Heat Transfer Coefficient (h): Energy transmitted per unit time from/to a surface of unit cross-sectional area.

  • Thermal resistance to convection: 1hA\frac{1}{hA}.

  • Nusselt Number (Nu): Enhancement in rate of heat transfer caused by convection over the conduction mode.

    • Nu=5Nu = 5: Rate of convective heat transfer is five times more than rate of heat transfer by conduction alone.

Dimensionless Numbers

  • Prandtl Number (Pr): Ratio of momentum diffusivity to thermal diffusivity.

    • Pr=να=μcpk\Pr = \frac{\nu}{\alpha} = \frac{\mu c_p}{k}

    • Indicates the relative thickness of the hydrodynamic and thermal boundary layers.

    • \Pr << 1: Molecular diffusivity of heat is larger than momentum diffusivity. Heat dissipates faster.

    • Gases: Pr0.7\Pr \approx 0.7

    • Water: Pr10\Pr \approx 10

Forced Convection

  • Forced Convection: Fluid motion is induced by external means (e.g., pump, fan).

  • Flow Regimes:

    • Laminar: Smooth, orderly fluid motion.

    • Transition: Intermediate regime.

    • Turbulent: Chaotic, disordered fluid motion.

  • Reynolds Number (Re): Dimensionless number indicating the flow regime.

    • Re=ρVDμRe = \frac{\rho V D}{\mu}

    • Laminar: Re < 2100

    • Transition: 2100 < Re < 10000

    • Turbulent: Re > 10000

Example 4.11

  • Water flowing at 0.02 kg/s is heated from 20°C to 60°C in a horizontal pipe (ID = 2.5 cm). The inside pipe surface temperature is 90°C. Estimate the convective heat transfer coefficient if the pipe is 1m long.

  • The solution is given by:

    • h=NNukD=(11.2)(0.633W/[m°C])(0.025m)=284W/(m2°C)h = \frac{NNuk}{D} = \frac{(11.2)(0.633 \, W/[m°C])}{(0.025 \, m)} = 284 \, W/(m^2 °C)

  • The convective heat-transfer coefficient is estimated to be 284 W/(m² °C).

Free Convection

  • Free Convection: Fluid motion is induced by density differences due to temperature gradients.

  • Rayleigh Number (Ra): Dimensionless number characterizing free convection.

  • Grashof Number (Gr): Ratio of buoyancy forces to viscous forces.

    • It's defined as: Gr=gβ(T<em>sT</em>)L3ν2Gr = \frac{g \beta (T<em>s - T</em>\infty) L^3}{\nu^2}

  • Film Temperature: Average of the surface and fluid temperatures.

    • T<em>f=T</em>s+T2T<em>f = \frac{T</em>s + T_\infty}{2}

  • Table 4.2: Coefficients for free convection calculations based on geometry and Rayleigh number range.

Example 4.15

  • Estimate the convective heat transfer coefficient for convective heat loss from a horizontal 10 cm diameter steam pipe. The surface temperature is 130°C, and the air temperature is 30°C.

  • Solution:

    1. Heat loss is by free convection.

    2. Film temperature: Tf=130+302=80°CT_f = \frac{130 + 30}{2} = 80 °C

    3. Properties of air at 80°C:

      • ρ=0.968kg/m3\rho = 0.968 \, kg/m^3

      • β=2.83×103K1\beta = 2.83 \times 10^{-3} \, K^{-1}

      • cp=1.019kJ/(kg°C)c_p = 1.019 \, kJ/(kg \cdot °C)

      • k=0.0293W/(m°C)k = 0.0293 \, W/(m \cdot °C)

      • μ=20.79×106Ns/m2\mu = 20.79 \times 10^{-6} \, N \cdot s/m^2

      • Pr=0.71Pr = 0.71

      • g=9.81m/s2g = 9.81 \, m/s^2

    4. Rayleigh number:

      • Gr=d3ρ2gβΔTμ2=(0.1)3(0.968)2(9.81)(2.83×103)(13030)(20.79×106)2=6.019×106Gr = \frac{d^3 \rho^2 g \beta \Delta T}{\mu^2} = \frac{(0.1)^3 (0.968)^2 (9.81) (2.83 \times 10^{-3}) (130-30)}{(20.79 \times 10^{-6})^2} = 6.019 \times 10^6

      • Ra=Gr×Pr=(6.019×106)(0.71)=4.27×106Ra = Gr \times Pr = (6.019 \times 10^6)(0.71) = 4.27 \times 10^6

    5. From Table 4.2, for a horizontal cylinder:

    6. Nusselt number:

      • Nu=0.6+0.387Ra1/6[1+(0.559/Pr)9/16]8/27=0.6+0.387(4.27×106)1/6[1+(0.559/0.71)9/16]8/27=22Nu = 0.6 + \frac{0.387 Ra^{1/6}}{\left[1 + (0.559/Pr)^{9/16}\right]^{8/27}} = 0.6 + \frac{0.387(4.27 \times 10^6)^{1/6}}{\left[1 + (0.559/0.71)^{9/16}\right]^{8/27}} = 22

    7. Convective heat transfer coefficient:

      • h=Nukd=220.02930.1=6.5W/(m2°C)h = \frac{Nu \cdot k}{d} = \frac{22 \cdot 0.0293}{0.1} = 6.5 \, W/(m^2 \cdot °C)

Overall Heat Transfer Coefficient

  • Combining convection and conduction resistances.

  • q=T<em>T</em>oRtotalq = \frac{T<em>\infty - T</em>o}{R_{total}}, where:

    • R<em>total=R</em>convection,inside+R<em>conduction+R</em>convection,outsideR<em>{total} = R</em>{convection,inside} + R<em>{conduction} + R</em>{convection,outside}

    • Rconvection=1hAR_{convection} = \frac{1}{hA}

    • R<em>conduction=ln(r</em>o/ri)2πkLR<em>{conduction} = \frac{\ln(r</em>o/r_i)}{2 \pi k L}

  • q=U<em>iA</em>i(T<em>iT</em>)q = U<em>i A</em>i (T<em>i - T</em>\infty), where UiU_i is the overall heat transfer coefficient based on the inside area.

    • 1U<em>iA</em>i=1h<em>iA</em>i+ln(r<em>o/r</em>i)2πkL+1h<em>oA</em>o\frac{1}{U<em>i A</em>i} = \frac{1}{h<em>i A</em>i} + \frac{\ln(r<em>o/r</em>i)}{2 \pi k L} + \frac{1}{h<em>o A</em>o}

Example 4.16

  • A 2.5 cm ID pipe conveys liquid food at 80°C. The inside convective heat transfer coefficient is 10 W/m²C. The pipe (0.5 cm thick) is steel (k = 43 W/mC). The outside ambient temperature is 20°C, and the outside convective heat-transfer coefficient is 100 W/m²C. Calculate the overall heat transfer coefficient and the heat loss from 1 m length of the pipe.

Heat Exchangers

  • Devices designed for efficient heat transfer between two fluids.

  • Assumptions:

    1. Steady state.

    2. Overall heat transfer coefficient is constant.

    3. No axial conduction.

    4. Well insulated.

  • dq=U<em>i(ΔT)</em>overalldAidq = U<em>i (\Delta T)</em>{overall} dA_i

  • dq=m<em>cc</em>pcdT<em>c=m</em>hc<em>phdT</em>h=U<em>i(T</em>hT<em>c)dA</em>idq = m<em>c c</em>{pc} dT<em>c = m</em>h c<em>{ph} dT</em>h = U<em>i (T</em>h - T<em>c) dA</em>i

  • Log Mean Temperature Difference (LMTD):

    • ΔT<em>lm=ΔT</em>1ΔT<em>2ln(ΔT</em>1/ΔT2)\Delta T<em>{lm} = \frac{\Delta T</em>1 - \Delta T<em>2}{\ln(\Delta T</em>1/\Delta T_2)}

  • q=UA(ΔT)lmq = U A (\Delta T)_{lm}

Example

  • Liquid food (inlet 20°C, exit 60°C, c<em>pc<em>p = 4 kJ/kg°C, 0.5 kg/s) is heated by hot water (inlet 90°C, c</em>pc</em>p = 4.18 kJ/kg°C, 1.0 kg/s).

  • Solution:

    1. Heat balance: q=m<em>cc</em>pcΔT=m<em>hc</em>phΔTq = m<em>c c</em>{pc} \Delta T = m<em>h c</em>{ph} \Delta T

      • (0.5)(4)(6020)=(1)(4.18)(90T2)(0.5)(4)(60-20) = (1)(4.18)(90-T_2)

      • T2=70.9°CT_2 = 70.9 °C

    2. The exit temperature of water is 70.9°C.

    3. Log-mean temperature difference:

      • ΔTlm=(9060)(70.920)ln((9060)/(70.920))=3050.9ln(30/50.9)=39.5°C\Delta T_{lm} = \frac{(90-60) - (70.9-20)}{\ln((90-60)/(70.9-20))} = \frac{30 - 50.9}{\ln(30/50.9)} = 39.5 °C

    4. The log-mean temperature difference is 39.5°C

    5. Heat transfer rate: q=(0.5)(4)(6020)=80kJ/sq = (0.5)(4)(60-20) = 80 kJ/s

    6. Heat exchanger length in the countercurrent case: L=(80kJ/s)(1000J/kJ)π(0.05m)(39.5°C)(2000W/[m2°C])=6.45mL = \frac{(80 \, kJ/s)(1000 \, J/kJ)}{\pi (0.05 \, m) (39.5°C) (2000 \, W/[m^2°C])} = 6.45 \, m

    7. For parallel-flow operation, assume the exit temperature will be the same as for counterflow, T2=70.9°CT_2 = 70.9°C.

    8. Log-mean temperature difference for parallel flow:

      • ΔTlm=(9020)(70.960)ln((9020)/(70.960))=31.8°C\Delta T_{lm} = \frac{(90-20) - (70.9-60)}{\ln((90-20)/(70.9-60))} = 31.8 °C

    9. The log-mean temperature difference for parallel flow is 31.8°C, about 8°C less than that for the countercurrent flow arrangement.

    10. Heat exchanger length in the parallel-current case: L=(80kJ/s)(1000J/kJ)π(0.05m)(31.8°C)(2000W/[m2°C])=8mL = \frac{(80 \, kJ/s)(1000 \, J/kJ)}{\pi (0.05 \, m) (31.8°C) (2000 \, W/[m^2°C])} = 8 \, m