Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Майнор / Книги / Chapter-5---Forced-convection_2021_Heat-Transfer-Engineering.pdf
Скачиваний:
0
Добавлен:
08.05.2026
Размер:
2.69 Mб
Скачать

5.6 Internal flows 157

We define such that = .

Tm Et mcvTm

Tm =

Ac

(ρu)(cvT )dAc

(5.217)

 

 

 

 

 

mc

v

 

 

 

 

For incompressible flow in a circular tube with constant cv , it follows that

Tm =

cv ρ r0 u T 2π r dr

(5.218)

ρ π r02 umean cv

 

0

 

Tm =

2

0r0 u T r dr

(5.219)

2

 

um r0

 

The bulk mean temperature, Tm is also called the mixing cup temperature.

Newton’s law of cooling

Once Tm is defined, we can use Newton’s law of cooling for heat transfer in this geometry (tube/duct flow)

q = h(Ts Tm )

(5.220)

where h is the local heat transfer coefficient, and Ts is the surface or wall temperature. There is an essential difference between external and internal flow. While in external flow Tis a constant, here Tm must vary in the flow direction. Stated explicitly, for heat transfer to occur, dTm /dx ≠ 0. dTm /dx is +ve, if the fluid is getting heated and

dTm /dx is -ve, if the fluid is getting cooled.

Fully developed conditions

In the hydrodynamic case, ∂u/ ∂x = 0 is valid for fully developed flow. But if heat transfer occurs, neither dTm /dx nor ∂T / ∂x equals zero at any radius r. T (r) changes continuously with x. However, this seeming paradox may be reconciled by working with a dimensionless temperature (TS (x)T (r, x)) / (TS (x)Tm (x)) and examining if this quantity becomes independent of x beyond a certain x. In fact, though T (r)

changes with x, (TS (x)T (r, x)) / (TS (x)Tm (x)) is not a function of x in the fully developed region. Hence, the mathematical condition for fully developed flow from

the viewpoint of convective heat transfer is

d

 

TS (x) T (r, x)

 

= 0

(5.221)

dx

 

T

(x) T

m

(x)

 

 

 

 

 

S

 

 

 

 

 

We can work a little with the dimensionless temperature considered above at r = r0, as follows.

 

 

 

 

 

 

 

 

T (r, x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

(x) T (r, x)

r

 

 

 

 

r r

f (x)

(5.222)

 

 

S

 

 

 

=

 

 

= 0

r

 

T

(x) T

(x)

 

T

(x) T

(x)

 

 

 

S

m

 

 

 

S

m

 

 

 

 

 

 

 

 

 

r=r0

 

 

 

 

158 CHAPTER 5  Forced convection

We know that

 

T

(5.223)

q = −k

 

 

r r=r0

 

or q / k (5.224)

= − T

r r =r0

Similarly, from Newton’s law of cooling, TS (x) Tm (x) = q/h. Substituting these results in, Eq. (5.222)

q/k

f (x)

(5.225)

q/h

 

 

 

h

f (x)

 

k

(5.226)

 

 

 

This is indicated in Fig. 5.10 for a constant k.

Hence, in the thermally fully developed flow of a fluid with constant properties, h f (x) and is a constant. One way of coming to peace with these “uneasy” arguments is to first state that from measurements, we know that Eq. (5.221) is valid, and using all our definitions and ideas of convection, we get to the fact that h f (x) in the fully developed region. If this too causes discomfiture, one can start with the argument that measurements show that h f (x) beyond a certain x and we would like to know what mathematical condition (Ts (x)Ts (r, x)) must satisfy in order for this to be true; this will logically lead us back to Eq. (5.221).

FIGURE 5.10

Variation of heat transfer coefficient, h, along the length of the tube.

5.6 Internal flows 159

Internal flow with constant heat flux, qw

From Newton’s law of cooling, q = h (Ts (x)Tm ) From the earlier discussion, we concluded that for fully developed internal flow, h is constant.

As q is constant in this case, the difference (Ts (x)Tm ) is also constant.

dxd (Ts (x)Tm (x)) = 0

dTs (x) = dTm (x) dx dx

Similarly, by performing differentiation on Eq. (5.221)

dxd (Ts (x) T (r, x)) = 0

dTs (x) = dT (r, x) dx dx

dTs (x) = dTm (x) = dT (r, x) dx dx dx

(5.227)

(5.228)

(5.229)

(5.230)

(5.231)

The key result we have now obtained is that the axial gradient of temperature is not a function of of the radial location under fully developed conditions for a constant heat flux.

From the energy balance, for ideal gases and incompressible liquids,

dqconv = mcpdTm = q p dx

Qconv = mcp (Tm,o Tm,i )

 

 

dTm

=

qp

=

p h (Ts (x)Tm (x))

 

dx

mc

p

 

mc

p

 

 

 

 

 

 

Integrating from x = 0 to x, we get

Tm (x) = Tm,i + qpx mcp

(5.232)

(5.233)

(5.234)

(5.235)

From Eq. (5.235), it is clear that the mean temperature variation is linear for constant heat flux. Tm (x) varies always linearly with x. Ts (x) varies linearly with x (for the constant heat flux case) in the fully developed region. These are qualitatively sketched for a typical case in Fig. 5.11.

Internal flow with constant wall temperature, Tw

or Ts

 

 

 

From Eq. (5.234), we know that dTm /dx = p.h(Ts (x)Tm (x))/mcp .

Let T = Ts (x)Tm (x), where Ts (x) is a constant

 

 

d ( T ) = −dTm

 

(5.236)

160 CHAPTER 5  Forced convection

FIGURE 5.11

Variation of surface temperature and mean temperature along the length of the tube.

dTm = −d ( T )

Substituting Eq. (5.237) in Eq. (5.234)

 

 

 

dTm

= −

d ( T )

=

 

p h T

 

 

 

dx

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

mc

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

To

d ( T )

 

 

= −

 

 

p

 

 

L h dx

 

Ti

T

 

 

 

 

 

 

m

cp

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρL

 

 

 

L

 

1

 

 

 

ln

To

=

 

 

 

 

h dx

 

 

 

Ti

 

m cp

 

0 L

 

 

 

 

T

 

= −

−pL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

o

 

 

 

 

 

h

 

 

 

 

 

Ti

m

cp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ts Tm,o

 

 

 

 

 

 

 

 

 

 

pLh

 

 

=

 

 

To = e mcp

 

 

Ts Tm,i

 

 

 

Ti

 

 

 

 

 

 

 

 

 

(5.237)

(5.238)

(5.239)

(5.240)

(5.241)

(5.242)

We could have also integrated from the tube inlet to some axial position x within the tube.

Ts Tm (x)

px h

 

 

 

mc

p

 

 

 

= e

 

 

 

 

(5.243)

Ts Tm,i

 

 

 

 

 

 

 

 

 

 

 

Qconv

 

(Tm,o Tm,i )

 

(5.244)

= mcp

 

 

Qconv = mcp

 

(Ts − Tm,i )(Ts − Tm,o )

(5.245)