Burt Rutan's White Knight and SpaceShip One, Photo Courtesy of Scaled Composites
Thermodynamics and Propulsion


17.2 Combined Conduction and Convection

We can now analyze problems in which both conduction and convection occur, starting with a wall cooled by flowing fluid on each side. As discussed, a description of the convective heat transfer can be given explicitly as

$\displaystyle \frac{\dot{Q}}{A}=\dot{q} = h(T_w-T_\infty).$

This could represent a model of a turbine blade with internal cooling. Figure 17.6 shows the configuration.

Figure 17.6: Conducting wall with convective heat transfer
Image fig9WallCondConv_web

The heat transfer in fluid 1 is given by

$\displaystyle \frac{\dot{Q}}{A} = h_1(T_{w1}-T_1).$

which is the heat transfer per unit area to the fluid. The heat transfer in fluid 2 is similarly given by

$\displaystyle \frac{\dot{Q}}{A} = h_2 (T_2 -T_{w2}).$

Across the wall, we have

$\displaystyle \frac{\dot{Q}}{A} = \frac{k}{L}(T_{w2} -T_{w1}).$

The quantity $ \dot{Q}/A$ is the same in all of these expressions. Putting them all together to write the known overall temperature drop yields a relation between heat transfer and overall temperature drop, $ T_2 - T_1$ :

$\displaystyle T_2-T_1 = (T_2-T_{w2})+(T_{w2}-T_{w1})+(T_{w1}-T_1) =\frac{\dot{Q}}{A}\left[\frac{1}{h_1}+\frac{L}{k}+\frac{1}{h_2}\right].$ (17..20)

We can define a thermal resistance, $ R$ , as before, such that

$\displaystyle \dot{Q} =\frac{(T_2 -T_1)}{R},$

where $ R$ is given by

$\displaystyle R = \frac{1}{h_1A} + \frac{L}{Ak} + \frac{1}{h_2A}.$ (17..21)

Equation (17.21) is the thermal resistance for a solid wall with convection heat transfer on each side.

For a turbine blade in a gas turbine engine, cooling is a critical consideration. In terms of Figure 17.6, $ T_2$ is the combustor exit (turbine inlet) temperature and $ T_1$ is the temperature at the compressor exit. We wish to find $ T_{w2}$ because this is the highest metal temperature. From (17.20), the wall temperature can be written as

$\displaystyle T_{w2} = T_2 - \frac{\dot{Q}}{Ah_2}= T_2 - \frac{T_2 - T_1}{R}\frac{1}{Ah_2}.$ (17..22)

Using the expression for the thermal resistance, the wall temperatures can be expressed in terms of heat transfer coefficients and wall properties as

$\displaystyle T_{w2} = T_2 - \cfrac{T_2 - T_1}{\cfrac{h_2}{h_1}+\cfrac{Lh_2}{k}+1}.$ (17..23)

Equation (17.23) provides some basic design guidelines. The goal is to have a low value of $ T_{w2}$ . This means $ h_1/h_2$ should be large, $ k$ should be large (but we may not have much flexibility in choice of material) and $ L$ should be small. One way to achieve the first of these is to have $ h_2$ low (for example, to flow cooling air out as in Figure 17.1 to shield the surface).

A second example of combined conduction and convection is given by a cylinder exposed to a flowing fluid. The geometry is shown in Figure 17.7.

Figure 17.7: Cylinder in a flowing fluid
Image fig9CylinderInFluid_web

For the cylinder the heat flux at the outer surface is given by

$\displaystyle \dot{q} = \frac{\dot{Q}}{A} = h(T_w-T_\infty)\quad \textrm{ at
}r=r_2.$

The boundary condition at the inner surface could be either a heat flux condition or a temperature specification; we use the latter to simplify the algebra. Thus, $ T = T_1$ at $ r = r_1$ . This is a model for the heat transfer in a pipe of radius $ r_1$ surrounded by insulation of thickness $ r_2 - r_1$ . The solution for a cylindrical region was given in Section 16.5.1 as

$\displaystyle T(r) =a\ln\left(\frac{r}{r_1}\right) +b.$

Use of the boundary condition $ T(r_1) = T_1$ yields $ b = T_1$ .

At the interface between the cylinder and the fluid, $ r = r_2$ , the temperature and the heat flow are continuous. (Question: Why is this? How would you argue the point?)

$\displaystyle \dot{q} = \underbrace{-k\frac{dT}{dr}}_{\substack{\textrm{heat fl...
...2}{r_1}\right)+T_1\right)-T_\infty\right]} _\textrm{surface heat flux to fluid}$ (17..24)

Plugging the form of the temperature distribution in the cylinder into Equation (17.24) yields

$\displaystyle -a\left(\frac{k}{r_2}+h\ln\left(\frac{r_2}{r_1}\right)\right)=h(T_1-T_\infty).$

The constant of integration, $ a$ , is

$\displaystyle a=\cfrac{-h(T_1-T_\infty)}{\cfrac{k}{r_2}+h\ln\left(\cfrac{r_2}{r...
...t)}
=
-\cfrac{(T_1-T_\infty)}{\cfrac{k}{hr_2}+\ln\left(\cfrac{r_2}{r_1}\right)}$

and the expression for the temperature is, in normalized non-dimensional form,

$\displaystyle \frac{T_1-T}{T_1-T_\infty} = \cfrac{\ln(r/r_1)}{\cfrac{k}{hr_2}+\ln(r_2/r_1)}.$ (17..25)

The heat flow per unit length, $ \dot{Q}$ , is given by

$\displaystyle \dot{Q}=\cfrac{2\pi(T_1-T_\infty)k}{\cfrac{k}{hr_2}+\ln(r_2/r_1)}.$ (17..26)

The units in Equation (17.26) are W/m-s.

A problem of interest is choosing the thickness of insulation to minimize the heat loss for a fixed temperature difference $ T_1
-T_\infty$ between the inside of the pipe and the flowing fluid far away from the pipe. ( $ T_1
-T_\infty$ is the driving temperature distribution for the pipe.) To understand the behavior of the heat transfer we examine the denominator in Equation (17.26) as $ r_2$ varies. The thickness of insulation that gives maximum heat transfer is given by

$\displaystyle \frac{d}{dr_2}\left(\frac{k}{hr_2}+\ln\left(\frac{r_2}{r_1}\right)\right)=0.$ (17..27)

(Question: How do we know this is a maximum?)

From Equation (17.27), the value of $ r_2$ for maximum $ \dot{Q}$ is thus

$\displaystyle (r_2)_{\textrm{maximum heat transfer}} = \frac{k}{h}.$ (17..28)

If $ r_2$ is less than this, we can add insulation and increase heat loss. To understand why this occurs, consider Figure 17.8, which shows a schematic of the thermal resistance and the heat transfer. As $ r_2$ increases from a value less than $ r_2 = k/h$ , two effects take place. First, the thickness of the insulation increases, tending to drop the heat transfer because the temperature gradient decreases. Secondly, the area of the outside surface of the insulation increases, tending to increase the heat transfer. The second of these is (loosely) associated with the $ k/hr_2$ term, the first with the $ \ln(r_2/r_1)$ term. There are thus two competing effects which combine to give a maximum $ \dot{Q}$ at $ r_2 = k/h$ .

Figure 17.8: Critical radius of insulation
Image fig9CylinderInFluidCriticalRadius_web





Muddy Points

In the expression $ 1/hA$ , what is $ A$ ? (MP 17.4)

It seems that we have simplified convection a lot. Is finding the heat transfer coefficient, $ h$ , really difficult? (MP 17.5)

What does the ``K'' in the contact resistance formula stand for? (MP 17.6)

In the equation for the temperature in a cylinder (17.25), what is ``r?'' (MP 17.7)

UnifiedTP