Heat Transfer 1. Conduction 2021
Heat Transfer 1. Conduction 2021
Heat Transfer 1. Conduction 2021
Volume 9
Heat Transfer 1
Conduction
Michel Ledoux
Abdelkhalak El Kami
Wiley
Heat Transfer 1
Mathematical and Mechanical Engineering Set
coordinated by
Abdelkhalak El Hami
Volume 9
Heat Transfer 1
Conduction
Michel Ledoux
Abdelkhalak El Hami
Wiley
First published 2021 in Great Britain and the United States by ISTE Ltd and John Wiley & Sons, Inc.
Apart from any fair dealing for the purposes of research or private study, or criticism or review, as
permitted under the Copyright, Designs and Patents Act 1988, this publication may only be reproduced,
stored or transmitted, in any form or by any means, with the prior permission in writing of the publishers,
or in the case of reprographic reproduction in accordance with the terms and licenses issued by the
CLA. Enquiries concerning reproduction outside these terms should be sent to the publishers at the
undermentioned address:
www.iste.co.uk www.wiley.com
Preface ................................................................................................. ix
Appendix 5. Table of erf (x), erfc (x) and ierfc (x) Functions .......... 295
Dividing the study of thermal science into three is the result of logic.
Presenting this work in three volumes is somewhat arbitrary; in our opinion,
x Heat Transfer 1
however, this split was necessary in order to keep the volumes in the
collection a reasonable size.
Within the Appendices, there are tables of error functions and their
offshoots: erf, erfc, ierfc, as well as reminders that are often essential for
hyperbolic functions. They also provide information on the notion of treating
certain non-stationary problems using Laplace transformations. Again, many
examples are given to illustrate another important aspect of conductive
transfers.
November 2020
Introduction
I.1. Preamble
Heating, lighting and operating the steam engines of the 19th century
were all very prosaic concerns. Yet this is where revolutions in the history of
physics began: the explosion of statistical thermodynamics driven by
Boltzmann’s genius, and quantum mechanics erupted with Planck, again
with Boltzmann’s invovlement.
through demonstration of black holes, the physics of which may still hold
further surprises for us.
The physics that we describe in this way, and to which we will perhaps
introduce some readers, is therefore related both to the pinnacles of
knowledge and the banality of our daily lives. Modestly, we will place our
ambition in this latter area.
The aim is to help learners who have not had high-level mathematical
training in their first years following the French Baccalaureate (therefore
accessible to apprentices), and pupils with more traditional profiles. At the
same time, we would like to show this broad audience the very new
possibilities in the field of digital processing of complex problems.
The revolution that digital tools has generated in the world of “science”
and “technical” fields, aside from the context of our daily lives, no longer
needs to be proved. We are a “has been” nowadays if we do not talk about
Industry 4.0. The “digital divide” is bigger than the social divide, unless it is
part of it...
Indeed, the memory of this revolution is now fading. Have students today
ever had a “slide rule” in their hands? Do they even know what it is? Yet, all
the physicists behind the laws of thermal science had only this tool in hand,
giving three significant figures (four with good visibility and tenacity),
leaving the user to find the power of ten of the result. It goes without saying
that a simple calculation of a reversible adiabatic expansion became an
ordeal, which played a part in degrading the already negative image of
thermodynamics held by the average student.
This reminder will seem useless to some; slide rules are at best sleeping
in drawers. But there is a moral to this story: no matter what type of
keyboard we type on, a calculator or a computer, our head must have control
over our fingers. This book has been written on the basis of this moral.
I.2. Introduction
Dividing the study of thermal science into three volumes is the result of
logic. Presenting this work in three volumes is somewhat arbitrary; in our
opinion, however, this split was necessary in order to keep the volumes in
the collection a reasonable size.
Before our readers immerse themselves in a text that, despite our best
efforts, remains intellectually demanding, we propose a short text that is a
little lighter.
I.3. Interlude
Let us entrust Ludwig Boltzmann to direct the film. Our B movie heroes
are getting agitated, delivering blows to one another. Each one of them has
moderate kinetic energy, distributed heterogeneously among them in the
room. For some reason, they get involved in a general brawl. Their average
kinetic energy becomes much greater.
E = 1 kT
C2
Using this model, we will return to the physical basis for all transport
phenomena.
We have just modeled the pressure, due to the transfer of the quantity of
movement on the surface, by the impact of molecules.
Let us now imagine that the altercation is initially located in the corner of
the room: a smaller group starts fighting between themselves.
From kicks to punches, after multiple impacts within the group and its
immediate neighbors, the agitation will spread: we have just seen the
mechanism of heat propagation by transfer of impacts.
Let us suppose that the great majority of the brawlers come from a ship
with a white uniform. Let us suppose that another boat in the port has
uniforms that are red. The red ones are initially all united. We will then
quickly see that the red mariners, as they receive and return blows, spread
out across the room. We have just shown the mechanism of diffusion of
matter, of a component within a mixture.
We can go further.
At the same time, these agitated clusters carry their disordered kinetic
energy, “thermal” agitation. We have just seen the mechanism of the thermal
boundary layer.
Finally, let us include a few red mariners in the crowd of white. They will
be carried with the clusters, and we have just invented the limit layer of
diffusion of a species.
All the individuals have a different speed, and the impacts are random, all
the bells start to jingle, each with a different frequency. The distribution of
frequencies will depend on the statistical distribution of speeds
(Boltzmann statistics), and the intensity of noise produced will depend on
the total agitation energy of the sailors.
NOTE.- the model is certainly simplistic. The emission comes from quantum
transitions in the gas atoms.
Here, we have already deviated from the pure substance of the book, but
we could go even further.
Let us suppose that our agitated sailors are in a room with one mobile
wall (a nightmare scenario frequently seen on the silver screen).
The incessant impacts of the fighters on this wall create a force that
pushes it. This force, reduced to a surface unit, explains the notion of
pressure.
By pushing against this wall, our crowd applies work that is greater than
the resistance.
We can see that the incidence of an average blow on the wall is rarely
normal.
Thus, we see that it will be impossible for the crowd (taken to mean a
gas) to give all its energy to a mobile wall.
book, will want to return to this text. They will then have understood, we
hope, the images that lead to the development of thermodynamics.
T=T(x,y,z,t)
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
2 Heat Transfer 1
The temperature field can only depend on a spatial variable. We say that
conduction is monodimensional or bidimensional.
This last remark divides the approach that we will adopt into two parts:
For many problems, the analytical approach is possible. This will be the
case, in particular, for the stationary or non-stationary problems with a single
dimension. In (most) other cases, a numerical approach is required.
1.2. Definitions
dT
dx
—— dT
grad T [1.2]
дУ
dT
dz
This property will be familiar to those who have studied electric fields.
4 Heat Transfer 1
ф=Q [1.3]
Density of the thermal flow ф, which can be even greater, is the flow
that passes through a unit of surface:
ф= Q
[1.4]
S tS
In this book, we will systematically note the flow with a capital Ф and
the flow densities with a lowercase ф.
These flows and flow densities may not be uniform through a surface,
which can also vary rapidly over time. We then deal with a definition on a
differential basis, through small surface areas and with small time intervals.
d Ф( t )= dQ [1.5]
dt
с!Ф d 2 Q
Ф = ^ = ^^ [1.6]
dS dt dS
The Problem of Thermal Conduction: General Comments 5
We have already pointed out that when we talk about thermal science,
thermodynamics is not far off. More precisely, it is its main principle.
1.3.1. Calorimetry
The minus sign indicates that one body is cooling and the other is
heating.
Q = mc A T [1.8]
6 Heat Transfer 1
Later on, Joule will demonstrate the equivalence of two forms of energy:
work and heat.
In principle, this is part of the repertoire of all high school pupils, but it is
not a bad thing to give a reminder of it.
becomes important. Moreover, we never recover all the heat that we provide
with the hot source in work. We always “throw” some into the cold source.
This will be the large Carnot work, among others.
Those who have closely read the text at the end of our preamble will
understand why. The blades of the mixer transform a “directed” energy into
a chaotic thermal agitation. Chaos will also never spontaneously return to
order.
All this comes from the genius (we defend our choice of word) of
Boltzmann, who will also be mentioned several times in the following
chapters.
To conclude, let us say that in all our transfers, we are dipping into the
realm of irreversibility and the permanent creation of entropy, but this will
be self-evident for the reader.
2
2.1. Introduction
2.2.1. Experiment
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
10 Heat Transfer 1
The whole apparatus is placed within a very thick layer of material that is
a very poor conductor of heat. Here, we make the hypothesis that the heat
dissipated by the electrical resistance passes all the way through the material
by thermal conduction.
The density of thermal flow is constant in the material, from left to right,
and is equal to:
P
Ф=- [2.1]
S
? = ATTT
—T [2.2]
e
Surface S
Thermocouples
1 = 1( x, y, z) [2.3]
1=1(T) [2.4]
1 = 0.12 W s-1K-1
Liquids are relatively poor conductors of heat. For example, for water at
20°C:
Gases are poor conductors of heat. For example, for air at 20°C :
REMARK.- Air is thus one of the best thermal insulators. However, between
two faces, separated by a layer of air of thickness of the order of 2 cm,
natural convection phenomena leads to a “parasitic” thermal transfer. Our
solution is thus to trap air in alveolar structures, such as networks of mineral
threads (glass fiber or rock “wool”). The solid material then leads to parasitic
thermal conduction, known as a “thermal bridge”, which explains how rock
wool presents a thermal conductibility of the order of h = 0.4 Ws-1K-1,
higher than that of air. If the faces are made of glass (double glazing), the
distance between the two panes must be small.
T1- T2
The law ф = h------ - implies that, in the case of the experiment, in other
heating and cooling plates. T(x) is then the temperature of a plane of the
material, perpendicular to Ox .
T1- T2
Indeed, the law ф = Л------- is valid for all distances e between two
planes perpendicular to Ox . This remains true for two planes that are
extremely close together. e= dx is then an infinitely small distance (as
small as we want), and the difference in temperature T1- T2 becomes a very
small difference dT. Fourier’s law is then written, in differential form, as:
dT
ф = л---- [2.5]
dx
Figure 2.3. Case of two infinitely close isotherms; relationship with the gradient.
For a color version of this figure, see www.iste.co.uk/ledoux/heat1.zip
The thermal flow that traverses dS1 and dS2 will be expressed as:
dT
d Ф = - 2---- d S1 [2.6]
dx
where x is an axis counted on the axis for n. Let us recall that the minus
sign indicates that heat propagates in the direction of decreasing
The Physics of Conduction 15
The same thermal heat flow evidently traverses dS', dS1 and dS2
dT dT
d Ф = - A---- dS1 ==-A----- dS cos в [2.8]
dx dx
REMARK.- Thus, the thermal flow that traverses a closed surface will be
written as:
Ф= jjAgradT.n dS [2.15]
S
a ) setting out a thermal balance. This balance can be written explicitly for
each problem, and this will be the case, in particular, in Chapter 4 concerning
quasi-stationary regimes.
We can also use an assessment written once and for all for a small region:
this will be the heat equation;
b ) writing special conditions for the problem, which will constitute the
conditions at the limits of the equation. This point will be examined in
greater detail in section 2.4, which follows.
- Cartesian coordinates:
T= T(x,y,z,t) [2.16a]
- cylindrical coordinates:
T = T (r ,0, z, t) [2.16b]
- spherical coordinates:
T = T (r ,0, Ф, t) [2.16c]
Often the device in question has symmetries that allow the problem to be
reduced to one spatial dimension. The temperature now only depends on two
variables, space and time. We can distinguish between:
- plane problems:
T= T(x,t) [2.17a]
- cylindrical symmetry:
T= T(r, t) [2.17b]
- spherical symmetry:
T= T(r, t) [2.17c]
For this reason, we will initially establish the heat equation in two
dimensions. This construction will be based on a direct expression of the
thermal balance.
For this purpose, we will establish the heat equation in four dimensions,
using the tools of vector geometry. For readers who are not accustomed to
these methods and to lighten the explanations a little, we have put this
establishment in the Appendix, simply giving its result below.
x + dx
Let us write
dT
ф( x ) = - A— [2.18]
dx
dT dT
— > 0 when the heat moves from right to left, and is negative
Л----
dx dx
when the heat moves in the direction of increasing x .
dT
The expression - A— is written at point x of the axis Ox . How is
dx
dT
- Л— written when we go from the x -axis plane to the x + dx -axis plane,
dx
with dx infinitely small?
dT
ф( x) = - Ad— is a specific f (x) function.
dT
Therefore, when we go from x to x + dx, -A— becomes
dx
ф( x) + d^xd=_ AdT+dL -AdT dx [2.19]
dx dx dx dx
20 Heat Transfer 1
We write:
dT d dT
^( x + dx ) = - A----- I----- ----- A---- dx [2.20]
dx dx dx
ddT d ( .TT 1
- A----- 1---- I -A----- I dx is positive when the heat moves in the direction
dx dx I dx )
of increasing x .
We will set up the balance of heat entering and exiting the section by
conduction.
When we set up this balance, we do not know1 what the direction of the
variation of T(x) is.
On the left-hand plane, if the heat goes from left to right, it is entering;
entering heat must be counted positively in the balance.
1 Since we are in the process of establishing the differential equation that will allow the
function T(x) itself to be determined!
The Physics of Conduction 21
dT
Therefore, this term will be written - 2— per m2 because it is positive
dx
if the heat moves in the direction of x (from left to right)
Through the plane x + dx, the entering heat must also be counted
positively in the balance.
Therefore, this is the opposite term, which will intervene in the balance,
in other words
d dT d L, dT L
2 +1 2I dx [2.21]
dx dx I dx j
REMARK.- The manipulation of the signs in front of the terms can surprise
or worry certain readers. It is fundamental to understand what is happening.
dT
It must be understood that the negative sign in the expression -2----
dx
indicates that the heat moves physically from hot to cold and that the flow is
positive if it moves in the direction of Ox , whereas the sign used in writing
a balance indicates whether the heat enters into the volume (necessarily a
positive term) or exits the volume (necessarily a negative term). Readers are
invited to take the time that is required to understand this correctly. This will
be of great benefit.
The overall heat balance in the small volume will therefore be:
ddT ddT d (d dT V d (d dT V
-2----- +- 2----- +---- 1 2 I dx = — I 2-----I dx [2.22]
dx dx dx ^ dx j dx ^ dx j
22 Heat Transfer 1
The term — I A^-T I dx is positive if more heat enters than exits. The
dx I dx )
small volume heats up.
— I A^- I dx is exchanged per unit of time. During the time dt, the
dx I dx j
quantity of heat exchanged with the outside is d I A^- I dxdt. It therefore
dx I dx j
modifies the temperature of the small volume of dT .
dQ = pcSdxdT [2.23]
Obviously,
d b
dQ = — ---- I1,
I AdT dx dt [2.24]
dx I dx j
S=1 [2.25b]
d T d (^ dT
Pc—= d x I d x [2.26]
dx
Let us note that we are now using d, because T = T (x, t) and the
derivatives are partial.
d T d2 T
---- = a----- [2.27]
dx d x2
2
a =---- [2.28]
Pc
± (2dT 1 = 0
[2.29]
dx ^ dx )
We note that T=T(x) and that we have therefore moved to the straight
'<” ’ that replaces the “ d ”, in particular partial derivatives.
24 Heat Transfer 1
д T д2 T
----= a----- [2.30]
д x д x2
dT A div gradrr
= ---- T [2.31a]
д t Pc
d T = a Дat
---- T [2.31b]
дt
In this case, the axis becomes an Oz axis and the heat (therefore the
variation in temperature) follows the vector radius r; positive orientation of
each radius from the cylinder axis to the outside of the cylinder. The
same reasoning given previously leads to the following conclusions.
д T 1 d {, д т
Pc a7= -----
r дr l
1 Ar —
дr
[2.32]
The Physics of Conduction 25
d T a д ( дT
---- =-------| r— [2.33]
дt r дr l дr
- stationary conduction:
d ( dT^ л
— [2.34]
drIl r dr)I = 0
As previously seen, here, we will provide the forms of the heat equation
in an axisymmetric coordinate system in a mono-dimensional case. To
26 Heat Transfer 1
lighten the text, we have included the details of the calculation in the
Appendix, based on the same approach as the one used above.
„ д T 1 д Г 2д T
p c---- =------- I A r ---- [2.35]
д t r2 д r l дr
д T a д ( 7d T
=I r 2 [2.36]
д t r2 д r l д r
- stationary conduction:
d ( dTT
--- I r 2---- =0 [2.37]
drl dr
The problems can also be dependent on time, in which case, they are
non-stationary. It is often the case that temperature is independent of time,
we then have a stationary problem.
The Physics of Conduction 27
In all cases, and no matter how sophisticated the problem is, there will be
two stages in resolving the problem:
1) determination of the temperature domain, generally T= T(x,y,z,t)
in the region D ;
2) determination of the flows at the interfaces, in general, from local
temperature gradients grad T (x,y,z,t) at these interfaces.
The heat equation can be set up by dividing the region into an infinite
number of small elements and by writing a heat balance for each of these
small elements. We will proceed in this way in the explanations given
immediately below. It will then be useful to set up this balance in two simple
cases, in the plane problem and in the axisymmetric problem.
28 Heat Transfer 1
A different method, faster and more elegant, consists of using the general
methods of differential geometry, which is associated with the previous
procedure in a way that is not immediately obvious.
Let us note that this heat equation in the most general case will include
time dependence. In addition, if necessary, it should take into account the
heterogeneity of the material in the region D .
This is particularly the case when the thermal exchange occurs with the
space outside the region by conduction, convection or radiation. In this way,
we can obtain various examples of the forms seen, since Te is a reference
temperature (distant temperature, external temperature of another material,
atmospheric temperature outside D , etc.):
,d T WT -eTe
T
Ф = ЛТ~ [2.38]
dn S R,uth
d> —
ф = Л- t
= £&\ [2.40]
dn S
The Physics of Conduction 29
We will see further examples of this later on. Some will be processed in
the section of the book on radiation, in particular, in Chapter 5.
NOTE.- The following examples have been placed at this point in this guide,
in line with the logical progression of the book and its presentations.
We have just examined the various ways of writing the heat equation.
The problems in this section are direct applications and solutions of this
equation. We show, in particular, the effect of a variation of the thermal
conductibility with temperature, the effects of the internal generation of heat,
and we examine some problems in spherical symmetry.
2 = 20 (1 + aT) [2.43]
Solution.-
—2—T = 0 [2.44]
dx dx
The Physics of Conduction 31
In this equation:
x =0 ; T=T1 [2.46a]
— A0 (1 + aT)—T = 0 [2.47]
dx dx
— (1 + aT )—T = 0 [2.48]
dx dx
(1 + aT )—T = C [2.49]
dx
d dT dT d T2 d ( T2A
(d 1 + aTd)— dT
T =—+ aT— = — + — a— =—I T + a— I [2.50]
dx dx dx dx dx 2 dx I 2 )
32 Heat Transfer 1
t2 1 ,
d I T + a— I = C,dx [2.51]
I 2 J 1
(_ T2 1 _
I T + a— I = C x + C7 [2.52]
I 2 J 1 2
( T2 1
x=0 ; T = T1 therefore I T1 + a— I = C1 *0 + C2 [2.53]
( T2 1
x=e ; T = T2 therefore I T2 + a— I = C1 * e + C2 [2.54]
It results that:
( T2 1
C7 = I T + a — I [2.55]
2 I 1 2 J
1 \( T2 1 ( T2 1
C, = — I T2 + a — I — I T + a — I [2.56]
1 e LL 2 2 JI 1 2 JJ
C2 = 38.06 [2.57a]
C1 = — 85.3 [2.57b]
Therefore,
ф = 2 (1 + ат )dx
—T = C1 2 [2.59]
т-т
Rth = T-T [2.61]
Rh = — = 5.05 m2 KW [2.62]
2.97
т-т
% = 2,-1---- 2 [2.63]
e
In an SI coordinate system:
^0 = 2.61W
34 Heat Transfer 1
The real flow is higher than this value. The real wall is less insulating
than the homogeneous wall.
T-T e
?с = [, with Rth0 = V [2.65]
Rth 0 Л)
The resistance of the real wall Rth = 5.05 m2 KW 1 comes out lower
than this value. We confirm that the real wall is less insulating than the
homogeneous wall.
The nuclear reactions that take place here produce a volumetric power q
expressed in qW.m-3 .
uranium melt in the center of the power plant, which is always best to
avoid?
Solution.—
2 d ( dT^
----- 1 r-----1 + q = 0 [2.68a]
r dr ^ dr )
R = — = 14.5 mm [2.68b]
2
d ( dT^ -qr
— 1 r---- 1 = ------ [2.70]
dr ^ dr ) 2
dT -qr2
r = + C1 [2.71]
dr 22
dT - qr C1
[2.72]
dr 22 r
The temperature must remain finite everywhere in the tube, which is not
the case in the center,
36 Heat Transfer 1
unless we use
C1 = 0 [2.74]
-qR 2
Te= +C2 [2.75]
42
qR 2
2e [2.76]
42
qR 2 2
T(r)-Te = 1 RJ [2.77]
42
qR2
t t t + q— [2.78]
max e 42
T = T + qmax R2 [2.79]
42
In other words
42(TF- T )
qmax [2.80]
The Physics of Conduction 37
4*27*(1132-200)
qmax ----- —i------- ------ -’ = 4.79.108 W.m-3 [2.81]
(14.5.10-3)2
EXAMPLE 2.3.- Maximum dissipation of an electrical resistance
We do not want the metal part of the resistance to exceed Tmax =150 °C.
The atmospheric temperature is Te =20 °C.
Solution.-
=1.075+2.12=3.197mKW-1 [2.83]
Д T 150 - 20
Ф =----=---------- = 40.66 W [2.84]
Rh 3.197
38 Heat Transfer 1
Solution.-
dT dT d dT7dT \
—2n 2— —2n 2----- l —2n 2 — I + 2nr dr q = 0 [2.87]
dr dr dr dr )
The Physics of Conduction 39
1 ddT q
r + =0 [2.88]
r dr dr X
1 ddT q
----- r---- + —= 0 [2.89]
rdr dr X
associated with a limit condition that we have imposed at the surface of the
wire:
ddT qr
— r----=------ [2.91]
dr dr X
dT qr2
r — = --—+ C1 [2.92]
dr 2X
dT _ qr C1
[2.93]
dr 22 r
At the surface,
r= R, we find T= Tw [2.95]
Therefore:
[2.96]
40 Heat Transfer 1
Finally:
T - T = qR- 1 - (r 1
[2.97]
w 4Л LRJ
For r = 0
Solution.-
1) The thermal power sought is equal to the incoming flow in the sphere
of radius r1 . This flow is equal to the flow traversing through each
The Physics of Conduction 41
For a value of r between r1 and r2 , the general expression of this flow is:
dT
Ф = 4п Ar2---- [2.101]
dr
dT
Ф = -4nAr 2---- = -4nAr 2 dT [2.102]
dr dr
. dT 2 . dT 2 d (. 2 dT ^ d (. 2 dT ^
-A—4nr2 + A—4nr2 + — \ A 4nr2 — I dr = — \ A4nr2 — I = 0 [2.103]
dr dr dx ^ dr J dr ^ dr J
d ( 2 dT
-^\ r -г =0 [2.104]
dr ^ dr
We note the term included below the derivative symbol. With this written
expression, we take into account the variation of the spherical surface
traversed by the flow when r varies from r to r+ dr.
d I 2 dT 1 [2.106]
dr 1^ r ----
— dr 1)= 0
2 dT r [2.107]
r T~ = C
dr
dT _C1
[2.108]
dr r2
C
T (r ) = - C1 + C2 [2.109]
r
C
r=r ; T = T1 ; T =---- + C2 [2.110]
r1
C
r = r2 ; T = T2 ; T2 =-------- -- + C2 [2.111]
r 2
T1 - T2 = - C1 I - - —1 [2.112a]
1 ri r2 )
c1 =-TT
1 -T\ [2.112b]
1 I111
I r1 r2 )
C
C2 = T + -L [2.113]
r1
C
—2 = T2 + _L [2.114]
r2
dr dr
Ф = - 4 nAr 2— r=r = -4 nAr 2 — 1 r=r1
[2.115]
dT d f- 1
dr 1 dr
Ф = -4nAr 2---- r=r = - 4nAr2---- I----L + C2 I r-r
I r 2) r=
-
= -4nAr 2— = -4 nAr 2 —1 [2.116]
r2 r=r1
dT -
Ф = -4 n Ar 2---- = = - 4nAr2r—
r=r2 2 r=rr2 = - 4n Ar2 —1
= [2.117]
dr
T-T
ф = 4 nr 2 T---- 2 [2.119]
Rth
44 Heat Transfer 1
By identification:
Rh =LT-1 [2.120]
^r1 r2
1 2
where Re is the Reynolds number based on the diameter of the sphere and
Д is the dynamic viscosity.
This equation will be valid for 3.5 <Re <80 000 and 0.7 <Pr <380
Solution.-
The second term gives the value of the Nusselt number when the
Reynolds numbers are zero. This corresponds to the absence of flow around
the sphere. Therefore, there is no longer any convective transport and term 2
corresponds to a purely conductive transfer, which in all cases coexists with
convection when the latter is present.
We will now write the balance for the space found between two spheres
of radii r and r+ dr:
d dT
Л—nr 2---- = 0 [2.122]
dr dr
d dT
л—nr 2---- = 0 [2.123]
dr dr
r ^ж ; T^T [2.124b]
0( r )= T - Te [2.125]
Defining:
00n = we
T -T [2.126]
d dT
--- r 2---- = 0 [2.127]
dr dr
r =R ; 0 = 00 [2.128a]
rd 00 = Cr
1---- [2.129]
dr
1 0
d 1 _ C
~11 [2.130]
dr r 1
-C
0( r ) = ——1 + C2 [2.131]
r
r =R ; =—
00,0R 1 + C22 [2.132]
r^^ ; 0 = C2 [2.133]
which gives:
C1=-R00 [2.134]
0, r )=R0 [2.135]
r
h= ^ = ^ = 2 = 22 [2.137]
Ta - T~ 0o~ R~ D
The Physics of Conduction 47
We deduce the classic expression for the Nusselt number from this:
Nu=D [2.138]
Nu = 2 [2.139]
D(t)=D0 (t)-Kt
where D0 is the initial diameter of the drop and K is a constant for which
the expression will be given.
This equation would be valid for 3.5 <Re <80 000 and 0,7 <Pr <380
Solution.—
1) The thermal flow ФW provided to the drop will, at instant t, have the
form:
m v = Pl ф w Lv [2.143]
4 3 D3
m (t )= —nR =n— [2.145]
36
Nu = 2 [2.146]
The Physics of Conduction 49
This result, which can be deduced directly (see previous exer cise), results
from a written theory in a stationary regime. When R (t) varies slowly with
temperature, we suppose that this Nusselt value remains valid.
Nu=D [2.147]
N Nu D
h = —— [2.148]
D
Finally:
ND
mv = PlLv 2nD*~D~(Ta - Tv ) = Pl Lv 2nDNu D(Ta - Tv) [2.149]
The balance for the mass of the drop is therefore written as:
dm d nD3
[2.150]
mv = dt ~ Pl dt 6
D* dD
vLL =- P
2 dt [2.151]
tm
v = PlLv 2nDNu D(Ta TV )= Pl 2 dt [2.152]
1 dD2
--------- = -L 4n N Л(T - T) [2.154]
2 dt Vr u \ a V)
where:
In other words:
dD = К
k [2.156]
dt
Writing:
K = Lv 8П Nu Л(Ta - T )
D2 =-Kt+C [2.157]
t=0;D(0)=D0 [2.158]
0 = D02 - К TV [2.160]
TV = D7 [2.161]
K
The Physics of Conduction 51
with
K = Lv 8П Nu Л(Ta - Tv ) [2.162]
This law is known under the name of the law expressed in D2. We
effectively note that the vaporization time is proportional to the square of the
initial diameter and proportional to the temperature differential between
atmosphere and saturation point.
When Nu = 2, we have:
D02
T [2.164]
V 16П Lv Л(Ta - Tv )
REMARK.- We have kept Nu until the end of the calculation because this
model can still be applied in the presence of wind blowing across the drop.
In this case, use the form proposed by Whittaker:
i 2
Nu=2+ 0.4 Re+0.06 Re p0.4 I Re [2.165]
r II Rw
We see that in this way, a vaporization drop continues to satisfy a law for
D 2. Convection obviously has the effect of increasing the constant K ,
therefore reducing the vaporization time. This result can be applied in
various practical situations: reactors in chemical engineering, agriculture
(spraying of pesticides), injection of petrol or diesel into a cylinder, rocket
engines, etc.
3
In fact, there are two thermal resistances that are very different:
a) the first is defined per unit of area of a wall. We will refer to it as
“thermal resistance”, which is denoted by Rth . We will use this almost
systematically in this chapter. We prefer to use this term because for a
homogeneous wall, this resistance is also a characteristic of the material it is
made of. It is generally related to the calculation of flow density.
b) the second is defined for a given surface of the wall. We will refer to it
as “global resistance”, which is denoted by RG .
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
54 Heat Transfer 1
[3.1]
e
If S is the surface area of the wall, the flow Ф that traverses it is then
T-T
ф = S ф = 2ST---- T— [3.2]
e
T1-T2 =Zr-Tt
ф= -1------- 2 [3.3]
e Rth
2
We have therefore defined the thermal resistance per unit of surface area,
Rh, or thermal resistance per m2:
e
Rth = [3.4]
2
U=R1h ■ 2 [3.5]
Conductance can also be used as commercial product data. For the sake
of simplicity, we will not use this parameter in examples proposed in this
book:
Conduction in a Stationary Regime 55
e -1
Rth = — is expressed in m2KW [3.6a]
U = — is expressed in W m 2 K 1 [3.6b]
e
T-T
ф = S у = 2ST-----T [3.7]
e
RG = 2S [3.9]
U = R- = — [3.10]
56 Heat Transfer 1
U = — is expressed in W mz K 1 [3.11b]
e
In this common case, the wall, for which we want to determine the
insulation performance, is made up of layers. We can still reason in terms of
density of thermal flow, because this is the same density that passes through
all materials.
The notion of thermal resistance per unit of surface area will therefore be
given preference. We can therefore obviously define a thermal resistance
e
Rh = 2 [ ]
for each component of the wall and a “resultant” thermal resistance Rth of
the wall as a whole.
Solution.—
We can set out the thermal resistance (per unit of surface area) for each
of the layers:
Rthi = eY [3.13a]
4
e
Rh 2 = ± [3.13b]
Z
e
Rh 3 = ± [3.13c]
Z3
We will denote the extreme temperatures applied to the two faces of the
wall as T1 and T4 , and the temperatures of the intermediate interfaces
between the layers of the wall as T2 and T3 .
The thermal flow density T1 can be written for each of the materials as:
71 - 72
9= [3.14]
Rh 1
72 — 73 [3.15]
9=
Rth 2
73 — 74 [3.16]
9=
Rth3
T1-T4 T-T
m =------------- 1--------- 4----------- = —---------- 4 [3.21]
Rth1+Rth2+Rth3Rth
Finally, we obtain:
T1-T4 T-T
m =------------- 1--------- 4----------- = —---------- 4 [3.22]
Rth1+Rth2+Rth3Rth
which leads us to the relationship between the overall thermal resistance Rth
and the thermal resistances of each wall:
We also note here a thermal analogy with Ohm’s law, assimilating the
density of thermal flow to an intensity and the temperatures to potentials.
b) Composite wall
In this case, the wall is no longer made up of layers that are simply
superimposed: it is heterogeneous. A problem of this kind, tackled from a
rigorous point of view, involves resorting to a numerical method. We can
find an approximate solution (more or less) to this problem using the notion
of overall thermal resistance.
Solution.—
We divide the composite wall into two plane walls made up of plane
layers:
- for the first wall, the concrete assembly and the two layers of insulation;
- for the second wall, the brick wall and its two layers of insulation.
We can define various overall thermal resistances for each layer of each
of these multilayer walls.
concrete:
ei [3.24]
RG =
XcS1
ei [3.25]
RGI =
XiS1
brick:
e1 [3.26]
G
RGB =
XbS2
60 Heat Transfer 1
R = 1 e [3.27]
GI2 0 c
AS
I 2
For each of the two sections of the composite wall, we therefore find:
1 -2T
ф = JT 2 T-T
2 2________
= ________ J1
[3.30]
RG1RGI1+RGC+RGI1
ф = JT T2
1 -2 =_______ J1T-2T2________
[3.31]
RG2RGI1+RGB+RGI1
The overall flow Ф traversing the complete wall is then given by:
R RG1 +Rg 2
[3.33]
G~ RG1RG2
№3
i■ 1I 2 Ii 3 ■3
Figure 3.1 . Thermal resistance. Electrical analogy for three superimposed walls.
For a color version of this figure, see www.iste.co.uk/ledoux/heat1.zip
The problem is a little more complex than the previous one. The
single-dimensional heat equation must be solved: we work in cylindrical
coordinates (r, в, z) and the temperatures vary only with the radius r. The
equation was given in Chapter 2, section 2.2.2:
d b
— dT I =I 0n
I Ar---- [3.34]
dr I dr )
d ( dT I л [3.35]
— I r---- I = 0
dr I dr )
This second-order equation must be solved with two limit conditions. The
temperatures on the two internal and external faces of the ring shape are
given by:
r=r; T=T
ee [3.36b]
first integration:
dT r
r—=C [3.37]
dr
Conduction in a Stationary Regime 63
dT C1
So — = — [3.38]
dr r
second integration:
T(r)=C1Lnr+C2 [3.39]
The two conditions at the limits will give the value of the constants:
Ti=C1Lnri+C2 [3.40]
Te=C1Lnre+C2 [3.41]
T-T T-T
=-i----rr
C11 e- = —i-----e- [3.42]
Ln — Ln —
reri
The value of C2 is not very important because to obtain the flow density,
we derive the flow per meter of length and the thermal resistance:
Ti- Te1
^ = - ^dT =—
C1
^ = - A----- - [3.43]
dr r r r
Ln—
ri
As we note, the flow density is not the same on both faces (the flow is the
same and the lateral surface is different):
internal face:
Ti- Te1
Pi = - e [3.44]
rr
Lnr ri
ri
external face:
T
i- e T 1
pe = - ^^—e~ [3.45]
rr
Ln r re
ri
64 Heat Transfer 1
The flow on a unit sleeve length will give (fortunately) the same result on
both faces:
T-T1 T-T
Ф i = 2п- yi
iii = 2i
nrik—rr = 2nk—r [3.46]
Ln- -i Ln—
rr
T-T1 T-T
= 2nr yi
Ф eei = 2nre
k—----
-- -— = 2nk—----
- - [3.47]
Ln— -e Ln —
r
Ln —
ri
R„,
th/m = [3.49]
2nk
T-T
Ф= L —---- r [3.50]
Rth/m
Conduction in a Stationary Regime 65
Vw = h (Tw - T) [3.51]
T-TT-T
Vw = h (Tw - Te ) = TVT = TT—- [3.53]
1 RthC
h
Rhe = h [3-54]
T-T T-T
Фw/m = 2nrh (Tw -Te) = 1 = -^ [3.55]
1 RthCL
2 nrh
66 Heat Transfer 1
1
RthCL = [3.56]
2 nr h
We note that the flow ФwL that escapes by convection (or that enters or
exits from the inside) from a length L of tube will thus be expressed by
T-T
ФwL = L -^-^
R
[3.57]
LU] ex I
1 +у +—L_
I rint J
Rth/m = [3.58]
2 nri ht f-gs 2 nX 2 nrehe
In this expression, we take into account the fact that the internal and
external convection coefficients can differ (possibly significantly if we have
a transfer to a liquid on the inside and to the air on the outside).
REMARK.- in the above formula, it goes without saying that, in the sum
of all the rings, we attribute the value of internal and external radii of each
element to re and ri (see examples for a better understanding).
T-T
Ф = L —---- e- [3.59]
Rth/m
This will obviously be the same on the internal and external faces.
Conduction in a Stationary Regime 67
In other words, Rth is the thermal resistance of the sleeve, resulting from
the resistance of the insulation and the resistance of convection.
We know that Rth has the expression for the case described here:
T fr1
Ln\ - I
r1
Rthth/. m =—— +--------- [3.60]
2nA 2nrh
and
T ' r1
_ h / m = JL
.d R . LnУ r-J I + 1 _ 1________ 1
dr dr 2nA 2nrh = 2nAr 2nr2 h
[3.61]
_
= 1 Г_1__ 1_’
2 nr A rh
68 Heat Transfer 1
In terms of signs:
1 1 1
sign = sign
2п2 2nrh 2 nr 2 ~~h
[3.62]
2
r-
rh-2 ___ h_
= sign Mi = sign 2r
r- rC
sign
2п2 2nrh = gn-jrC = sign (r -rC ) [3.63]
A polar explorer is invited by the Inuit into their igloo. The explorer, a
trained physicist then ponders the thermal behavior of this type of building.
The igloo is built of ice. The thickness of the walls is e= 40 cm. The
total surface area of the walls of the igloo is 5 = 25 m2
We will take into account a convection coefficient h for both inside and
outside the igloo.
Remark.- The Inuit coat themselves in seal fat for better thermal protection.
They also sleep fully clothed. Polar exploration means having to change
Western habits.
70 Heat Transfer 1
What thermal power P must the wood fire inside the igloo produce?
- all calculations assume that the walls have a practically plane geometry;
- losses through the ground (ice) and the smoke hole will be neglected
(the Inuit do not like to suffocate);
- the thermal conductibility of the ice is: AG = 2.1 W. m 1. K 1;
- the convection coefficient (both internal and external) is:
h= 10W.m-2.K-1.
Solution.-
The energetic power of the fire, that is, the energy that it liberates per unit
of time, must compensate for the thermal losses, meaning the thermal flow
resulting from the thermal resistance that has previously been calculated.
Thus:
25 8+40
P =Ф =S T -T [3.65]
Rth 0.39
P=3077W [3.66]
We give:
e =70.10-2m
2 = 3 W. m -2. K 1
e =25.10-2m
Solution.—
1) We must first calculate the overall thermal resistance Rth of the wall.
The resistance to transfer of each convection coefficient will be: Rth conv = —
h
e = 70.10-2 m
2=3W.m-2.K-1
e = 25.10-2 m
Using the additive nature of the resistances to the transfer, the resulting
thermal resistance is therefore:
3) The flow of outgoing heat for a surface area S of the wall can be
calculated by:
T-T
ф = S ф = S-----e- [3.70]
Rth
Solution.—
The thickness of each pane of the glass is 3.5 mm and the layer of air is
12 mm thick. The thermal conductivity of the glass is 0.7 W.m-1.K -1and for
the air, it is 0.024W.m-1.K-1.
Solution.—
Rth = J [3.74]
e = 3.5.10-3 m
J= 0.7W.m-2.K -1
- =1.2.10-2 m
J= 0.024W.m-2.K -1
The thermal losses per m2, that is, the density of flow of losses, is given
by:
ДT
Ю =----- [3.78]
Rh
ДT = 5C [3.79]
ю =— = —^ = 1000 W. m-2
[3.81]
Rthv 5.10—3
«—90-►
Figure 3.6. Dimensional drawing of the French window
Conduction in a Stationary Regime 77
Calculate the thermal flow Ф that passes all the way through the French
window.
Solution.—
Rth = J [3.82]
e= 6.10-3 m
J= 1.05W.m-2.K -1
-
Two glass panes A, with the thermal resistance RhA = j , and a layer of
e =1.5.10-2 m
2 = 0.026 W. m -2. K 1
Rr = - = 0.577K.m2W [3.84]
2
RthB=2*RthA+Rthair=0.588K.m2.W-1 [3.85]
Three glass panes A, with thermal resistance RthA , and two layers of air
with thermal resistance Rthair :
RthC=3*RthA+2*Rthair=1.171K.m2.W-1 [3.88]
2) Concrete wall
-
A concrete wall, with thermal resistance Rthconcrete = -j, a wall of rock
wool with thermal resistance RthLV , and a wall of BA13 plasterboard with
the thermal resistance:
RBA !3 = 2 [3.89]
Conduction in a Stationary Regime 79
e = 20.10-2 m
2 = 2 W. m -2. K 1
e =15.10-2m
2=0.041W.m-2.K-1
e =1.3.10-2m
2= 0.25W.m-2.K -1
RthB=Rthconcrete+RthLV+RthBA13=0.1+3.66+5.2.10-2=3.81K.m2.W-1 [3.92]
3) French window
This flow will be equal to the sum of the flows traversing the wood and
the window.
S=2.1*0.9=1.89m2 [3.93]
80 Heat Transfer 1
The surface area of the glass SV through which the flow passes is:
The surface area of wood SB through which the flow passes is equal to:
SB =1.89-1.05=0.84m2 [3.95]
The density of flow pBois passing through the wood will be:
with A T = 20 - 5 = 15 °C
and
e=4.10-2m
2 = 0.23 W. m K 1
Thus:
15
^wood = 0T77 = 0.8621W. m ’2 [3.98]
.
Vgla.. = P-99]
thglass
with
AT=20-5=15°C [3.100]
Conduction in a Stationary Regime 81
Thus:
15
<P^ = — = 12.81 W. m [3.102]
.
Slate
Batten
The temperature of the upper face of the slates is Te =30°C, and the
temperature of the attic is Ti =10°C.
82 Heat Transfer 1
Solution.—
eslate =3.10-3m
Aslate =2.2W.m-2.K-1
The slate, with thermal resistance Rtha = 1.363.10-3 K.m2.W 1 and the
batten with thermal resistance Rthvol :
e =1.5.10-2 m
A= 0.12W.m-2.K -1
Conduction in a Stationary Regime 83
S=2*(20*7)=280m2 [3.108]
AT
?b = T~
[3.109]
thB
A T = 30 -10 = 20 °C
Ф = SVb [3.111]
Solution.—
e = 20.10-2 m
2 = 1.4 W. m K '
e =15.10-2m
2LV = 0.045 W.m-2.K-1
The flow of the outgoing heat for a surface area S of the wall can be
calculated by:
T-T
Ф= S-i---- e- [3.117]
Rth
21 - 3
Ф= 102----- = 500.3 W [3.118]
3.67
P = ф= 500.3 W [3.119]
1) What is the thermal resistance per meter, Rth1 , of a rock wool sleeve
with an internal diameter di = 6cm and an external diameter de = 10cm.
The thermal conductibility of the rock wool is A = 0.045 W. m 1. K 1 .
2) We subject this to a temperature difference of A T = 25 °C and do not
take into account the convection. What is the flow Ф1 passing through a
length L = 15 m of the sleeve?
3) We take into account the two convection coefficients: an internal
coefficient hi = 5 W m-2K-1 and an external coefficient he = 8 Wm-2K-1.
Given these conditions, recalculate the flow Ф2 through 15 m of this
sleeve.
86 Heat Transfer 1
Solution.—
d i —10
Ln -e Ln
Rth1
dt 6 1.807mKW-1 [3.120]
2 nA 2 п* 0.045
т d-e
Ln
d 11
R th 2 =------
_ 1— +-------------- +
2 nA 2 nrtht 2п( rt + e) he
1 10
Ln
------- 6----+--------12------- + 1
6.316mKW-1 [3.122]
2п*0.045 2п*0.01*5 2 n*0.015*8
We want to reduce the flow of heat transferred by this tube to the room in
which it is found.
Solution.—
AT 200 - 20
Ф = =170W [3.125]
Rth1 1.061
REMARK.- We could have calculated this flow directly using the definition
of the convection coefficient. Calculating the flow through a lateral surface
of the tube of 1 m in length:
Ф1 = Sh A T = ndih AT
88 Heat Transfer 1
This obviously gives the same result, since the thermal resistance is
defined using this last expression!
d
re2 = 0.25 m
=— [3.126]
The thermal resistance takes into account the insulator and the
convection:
T (r ^ < 0.25 ^
Ln 1 — 1 Ln 1------ 1
„ IrJ 1 I 0.015 J 1
R^7 = + = + [3.127]
th2 2 nA 2 nreh 2 n *0.045 2 n 0.25*10
ДT
Ф2 = — = 10.01 = 17.98 W [3.129]
11
Rt,v = ---------- =------------------- = 6.37.10-2 mKW [3.130]
thconv 2nr h 2n* 0.25 *10
T-T T-T
Ф2 = TW------- e- = - WL-------- e = 17.98 w [3.131]
2 Rthconv 6.37.10-2
Therefore:
Moreover:
a) Calculate the thermal resistance per meter between the fluid and the
external atmosphere.
b) Calculate the thermal flow leaving the fluid and transferring to the
outside, for a length of 10 m of tube.
a) Calculate the new value of thermal resistance per meter between the
fluid and the external atmosphere.
b) Calculate the new value of the thermal flow leaving the fluid and
transferring to the outside over a length of 10 m of tube.
Numerical data
Solution.—
Rth [3.134]
2 пл
d
For the rock wool sleeve, the internal radius is r, = -2- and the external
radius is re + e:
. \ 3+5 1
Ln\ I
R^ ooi =----
thwool 2п 0.047- = 3.321 m2.KW
- ------- [3.135]
Therefore:
170 - 20
Ф = =10* = 451.7W [3.136]
Rth 3.321
Ln\ Ln (3 I
Rh = - ri =^-^ = 1.24.10-3 m2.KW-1
±±= [3.137]
th 2пХ 2*п*53
b) The thermal flow leaving the fluid and moving to the outside for a
10 m-long tube is calculated in the same way as above with the new thermal
resistance:
AT 170 - 20
Ф=----- = 10*----------- = 403.3 W [3.141]
Rh 3.72
2) Give the expression for the thermal resistance Rthtot of the assembly
tube + sleeve.
3) Determine the conditions for which the addition of a sleeve does
indeed lead to a reduction in thermal losses.
Data:
Solution.—
1) Let us calculate the thermal resistance of the tube alone. We will deal
with the problem for 1 m of tube.
The thermal resistance Rth per meter of a cylindrical sleeve, with internal
radius r1, external radius r2 and thermal conductibility 2, is calculated using
the expression:
Ln\
D _ I
\ r1 [3.142]
h = _ 0
2 П2
Ln\
Rth, (L) = thR.L = L [3.143]
2 ПЛ
Ln I Ln I -2-
-e- I
rr
Rh (L = 1 m ) =--------— =--------^-1- [3.144]
2 n2 2 n21
Ln I — I Ln I —
rr
R., (L = 1 m ) =-------- — =-------^2- [3.145]
hsleee ( ) 2 n2 2 n21
R =—1— [3.146]
thconv
2nrh
94 Heat Transfer 1
Ln \ — I Ln \ — I
,, l Г1 J lrJ 1 1
h = ~--------------------- +------------ +------------ +----------
Rthtot [3.147]
2 n11 2 n11 2 n r1 hi 2 n r3 —
Ln\r- I
+J [3.148]
2 n11 2 n r3 he
1
rc =± [3.149]
dR
- when r3 < rc, the derivative — is negative. An increase in r leads
dr
to a reduction in Rth . Increasing the thickness of insulation leads to an
increase in losses;
1
- when r3 becomes greater than the critical radius r3c = — , it becomes
r33c = —
6 = 1.67 cm [3.150]
We verify that r3c is indeed greater than r2. In addition, the increase in
the thickness of the sleeve insulation becomes effective when it is greater
than 1.67-1.5=0.17cm; in other words, 1.7 mm is low. We see that, in
practice, the condition is not very constraining.
Solution.—
1) Sleeve of D= 12cm
We calculate the thermal resistance Rth of the sides of the bottle, using
the equation:
rr
Ln — Ln—
r= r_ [3.151]
2n A 2n A ’
Di = 5.10-2 m
De =12.10-2 m
A=1.05 W.m-2.K-1
he = 8W.m-2.K -1
Ln D Ln 12
R = ------- i- =------------ 5 = 0.133 W. m 2.K [3.152]
th 2nA 2* n*1.05
L
Rtot
---- D- +--- 1 = 0.133 + 1 ,
2 nA 2 nrehe 2* n*6.102*8
=0.133+0.3316=0.464W.m-2.K-1 [3.153]
The flow of heat Ф exiting the system per meter for Tt = TB = 70 °C and
Te = 21°C is:
2) Sleeve of D= 26.3cm
We calculate the thermal resistance Rth of the sides of the bottle, using
the equation:
rr
Ln — Ln—
Rth = =----- - [3.155]
2nA 2nA
Di = 5.10-2 m
De = 26.3.10-2 m
A = 1.05 W.m-2kK4
he = 8W.m-2.K -1
LnD Ln —
Rh = ------- i- =-------- 5----- = 0.251 W.m-2kK [3.156]
th 2 nA 2* n*1.05
LnD
R = D. + 1 = 0.133 +------------- 1------------
tot 2nA 2nrehe 2*n*13.1510-2 *8
=0.251+0.151=0.402W.m-2.K-1 [3.157]
The flow of heat Ф exiting the system per meter for T = TB = 70 °C and
Te = 21°C is:
3) Sleeve of D= 36cm
We calculate the thermal resistance Rth of the bottle walls, using the
equation:
rr
Ln — Ln—
r =__ r_ [3.159]
2пЛ 2пЛ
Di = 5.10-2 m
De=36.10-2 m
Л= 1.05 W.m-2.K-1
he= 8W.m-2.K -1
Ln D Ln 36
Rth = ------- l- =-------5 = 0.299 W.mz.K [3.160]
th 2пЛ 2*п*1.05
LnD
1
Rtot ------ — +----------- = 0.299 +
2 пЛ 2 nrehe 2* п*18.10-2 *8
=0.299+0.1105=0.409W.m-2.K -1 [3.161]
The flow of heat Ф exiting the system per meter for Tt = TB = 70 °C and
Te = 21°C is:
*, = T^ = 70-21 = 119.8 W
[3.162]
0 tot
0.409
Conduction in a Stationary Regime 99
r = 18cm > rC ; Ф = 119.8W, the initial radius is larger than the critical
radius, r=13.15cm>rC , the thickness of glass increases, and the flow of
losses decreases.
Solution.—
b) This resistance is expressed per meter. The flow loss will therefore be:
_ , AT
Ф= L---- [3.167]
Conduction in a Stationary Regime 101
64
ф = 5000------ = 1.98.105 W [3.168]
1.613
2) We will therefore need to add the resistances per meter for the foam,
the steel (which is calculated using an identical equation) and the convection
coefficient.
r+e r
Ln —----- Ln —
riri1
Rth2 — + + [3.169]
2 п ЛmSi 2 п Л„ 2 nh
3) We will write a thermal balance between the input and the output
of the tube: the difference between the incoming heat flow and the outgoing
heat flow, with the flow, is equal to the thermal losses.
where p = 1000 kg.m-ъ is the density of water and c = 4180 J.kg 1 is its
specific heat capacity (values which should be known in principle). Tin and
Tout are, respectively, the input and output temperatures of water in the tube:
The loss of 17°C represents 24% of the input temperature. This is quite a
rough approximation.
4
Quasi-stationary Model
b) The temperature varies slowly over time: we then suppose that the
transfer to the exterior can be calculated using relationships taken from the
stationary theory. In addition, the thermal balance, as established, assumes
that the temperature is constant.
We can observe that these hypotheses are not sufficient, in several ways:
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
104 Heat Transfer 1
We will give the basis for the calculation here, by way of an example.
So, for a homogeneous body with volume VOL, density p and specific
heat capacity c, its mass will be:
m=pVOL [4.1]
During the time dt , the quantity of heat exchanged by the body and the
exterior is given by:
dQ=pcVOLdT [4.2]
This quantity of heat given or extracted from the body will be due to the
thermal transfer, which can generally be written as:
ф=SH-T
[4.3]
The balance will therefore be written as follows, taking into account the
above equations:
dQ =-Ф dt [4.4]
The minus sign is a result of the remark made above. We note that its
presence is a result of the choice made when writing dQ (choice of the
thermodynamic sign convention).
T(t)-T
dQ = p cVOL dT = -Ф dt = - -dt [4.5]
dT T(t)-T-
[4.6]
dt pcVOLRth
4.3. Resolution
O( t )= T (t)-Te [4.8]
d0 _ 0
[4.9]
dt PcVOLRth
d0 dt
------- = - adt [4.10]
~0 PScRth
with
1
a =------------ [4.11]
psc VOL Rth
Ln 0 = - at + LnC [4.12]
and
t= 0 T0 -Te =C [4.15]
which results in the final expression for the temperature of the bar:
T(t)-Te=(T0-Te)exp(-at) [4.16]
REMARK.- the radiation manifests itself here by an imposed heat flow term.
This example has a natural place in a chapter that is dedicated to conduction,
and not to radiation, where the latter field does not require any specialized
calculations.
3) How long will the temperature take to reach Ti =25°C? How long will
the temperature take to reach Ti =70°C?
4) Here, the model is very simplified. Under real conditions, do you think
that a temperature of 70°C can be reached in the room? Why?
Solution.—
the notation and values associated with this are classic pi = 1.3 kg.m-1,
c= 1000 J.kg-1.K-1
dT_ Ф S
[4.20]
dt pVOLc
2) Resolution is immediate.
ФS
T=--------------- t + Cst [4.22]
P VOL c
ФS
T (0) =-------- 0 + Cst = Tt [4.23]
p VOL c
Therefore:
t=T + Ф S t [4.24]
PVOLc
Numerically:
T =18+3.08.10-2 t [4.26]
4) The numbers are not realistic, above all the second, because at a
temperature of T=70°C, it will not be easy to neglect the thermal losses in
the building!
Where necessary, we will use the numerical data that is also given at the
end of the text.
The walls are made of concrete with a thickness of 20cm . The roof is
made of a concrete terrace area of the same thickness.
We insulate all the walls and the ceiling with a 25cm layer of expanded
polystyrene and BA13 plasterboard (13mm thick).
That is, hi and he are the internal and external convection coefficients
respectively (valid for the walls and ceiling).
3) Radiator breakdown.
The radiator breaks down. We presume that the thermal losses can still be
evaluated, as they can in the stationary regime, when the internal
temperature of the garden hut Ti (t) varies with time. The external
temperature Te remains constant.
We will call Rth the total thermal resistance of the garden hut, S the total
surface area of the walls and Vol the internal volume of the garden hut.
a) Write the expression for the instantaneous flow of thermal losses from
the garden hut.
b) Writing the thermal balance for the garden hut, give the differential
equation for the behavior of Ti (t). It will be useful to introduce the function
#( t )= T (t)-Te
4) Numerical application.
Te =6 °C, as above
T,0 =20°C
After the breakdown of the radiator, how long did the temperature take to
fall to T, =10 °C inside the garden hut? Give this time in hours and minutes.
b) Te and T,0 maintain the previous values of 6°C and 20°C. How long,
in hours and minutes, will it take for the temperature of the interior of the
garden hut to descend to T, = 10°C with this new insulation?
Additional data
hi =he =5W.m-2.K-1
Solution.—
e =20.10-2m
2 = 1.4 W. m к-1
e =25.10-2m
2 = 0.035 W. m K-1
e =1.3.10-2 m
2= 0.25W.m-2.K-1
The surface area which the leaks are distributed over is:
S=2*(7+5)*2.5+7*5=95m2 [4.32]
114 Heat Transfer 1
AT 14
Ф = S<p = S----- = 95* = 171.83 W [4.34]
R RhB 7,17
3) Radiator breakdown
a) We will give a general expression for the instantaneous flow of
thermal losses from the garden hut.
Ф=ST^
[4.35]
During a time period dt, the flow of heat exiting the garden hut Ф dt
leads to cooling dTi such that:
Ф dt = - pCVVOLdTi [4.36]
The minus sign expresses that the outgoing heat leads to the cooling of
the room.
Replacing Ф by its expression and introducing 0 , we obtain the
differential equation:
PairCvVoLdT = — S° [4.37]
dt R
= -а0 [4.39]
dt
with
S
а =------------------------------ [4.40]
PairCvVoLRthB
c) Calculation of Ti (t)
In other words:
00=Ti0-Te [4.41]
d0
= -a0 [4.43]
dt
d0
—=-adt [4.44]
0
dLn0 = -at + LnC [4.45]
0( t ) = C exp-at [4.46]
0( t ) = 00 exp-at [4.47]
116 Heat Transfer 1
4) Numerical application
Ti0 = 20°C
We are trying to find how long it takes after the radiator breakdown for
the temperature to drop to 10°C inside the garden hut.
Calculation of a gives
a=1.5.10-4 [4.49]
-14
t =---- Ln— = 8352 5 [4.51]
1.5.10-4 14
t=8352s=10h19mn12s
c) With this second solution, the cooling time is twice as weak as for the
initial insulation! Obviously, this solution is less satisfactory than the
previous one.
A steel block has a rectangular base with sides a and b and height H. Its
specific heat capacity is denoted c.
The cube can only transfer energy to the exterior via its other faces, for
which a convection coefficient h applies.
1) Give the expression of the flow density ф1, then of the flow Ф1
coming out of the block at the initial moment in time.
118 Heat Transfer 1
Deduce from this, by writing the instantaneous heat balance in the block,
the differential equation that describes the evolution of the temperature of
the block.
4) Numerical applications
We give
Density of steel:
h= 5W.m-2.K-1
Te =20°C
T0 =150°C
b) How long will it take for the difference between the temperature of the
block and the outside to descend to 10% of its initial value?
Quasi-stationary Model 119
Solution.—
Ф( t ) = Sh [ T (t)-Te ] [4.58]
We define a variable
0 = ( Tt - Te) [4.62]
do
= -a0 [4.63]
dt
Sh
a == P.e.Vo.c = 3225'10 [4.64]
120 Heat Transfer 1
d0
—= -a0 [4.65]
dt
d0
— = -adt [4.66]
0
Integrating the two terms:
0=Cexp-at [4.68]
t=0;0=00 [4.69]
4) Numerical applications
a) Values of ф1 and Ф1
0 = 00 exp - at [4.75]
00=150-20=130° [4.76]
a= 3.225.10-5 [4.77]
Quasi-stationary Model 121
t =-1ln— [4.78]
a —0
— = 0.9 [4.79]
—0
t =—ln — [4.80]
a—0
The total surface area of the walls (except the windows of course) is
covered with wood of thickness 8 cm.
-------------100 m
40 m
20 m
1) Calculate the thermal resistance (per m2) of the walls and windows
RthWall and RthWindows , respectively, taking into account the wooden
surrounding and the convection. Calculate the respective surface areas of the
wall, the terrace and the window.
a) Write the relationship between the instantaneous thermal loss rate and
the instantaneous internal and external temperatures Ti (t) and Te(t)
respectively.
Quasi-stationary Model 123
b) Write the thermal balance equation that connects the derivative of the
internal temperature with respect to time and this thermal loss.
d) How long will it take for the temperature to become 6°C inside the
castle? Give this time in hours and minutes.
It was cold in Versailles during Louis XIV’s era. Can you understand
him?
Density of air:
p = 1.3 kg. m -3
cv =718J.kg-1.K-1
Solution.—
=e_gasL + 11
+.1
tV A. h h [4.83]
glass i _
wall = 5400m2
S [4.85b]
ф= V T___
_J -TL + V T-T
_>___ L
wall glass [4.86]
thwall thglass
99
ф = 5400------ + 200-------- = 4.992.104 W [4.87]
1.179 0.207
Ф_ о T (t)- T. + о T (t)- Te
ф °wall p 1 о glass & [4.88]
thwall thglass
b) Balance equation between the outgoing flow and cooling of the air in
the castle:
Ф dt = - pCVVOLdTi [4.89]
The minus sign expresses that the outgoing heat leads to the cooling of
the room.
p cvdT=-о +v W—l
pair~^V O OL w wall glass [4.91]
dt Rthwall Rthglass
dT = d (T — Te) = de
[4.92]
dt dt dt
de
= -ав [4.93]
dt
with
оо
Owall + gljslass
a R,i ,,
thwall R„,
thglass
[4.94]
pair CV Vni
OL R.thh
126 Heat Transfer 1
a gives
SS
S wall + glass
We define
00=Ti0-Te [4.98]
d0
— -~a0 [4.99]
dt
d0
---- =-ddt [4.100]
0
dLn0=-at+LnC
t = tf; 0 = Of = Tf - Te = 6 - 5 = 1 [4.103]
-1 1
tf =---------- - Ln = 29572 5 [4.105]
f 7.43.10-5 9
t=8h12mn525 [4.106]
Taking into account the powers to be applied, the phenomenon must have
occurred often.
1) In the original construction of the farm, the walls of the farm are made
of granite of thickness e= 80 cm. The walls are not insulated.
a) Give the expression and the value of the thermal resistance per m2
Rth11 of the farm wall. Take into account an internal convection coefficient
b) To estimate the losses through the roof, we will consider that they are
added to the losses through the walls and are equal to 40% of the losses
through the vertical walls.
128 Heat Transfer 1
For the calculation, we will assimilate the door and windows (small
surface areas) to sections of wall.
3) The heating breaks down when it is Ti0 = 19°C inside the farm, and
the external temperature is still Te =5°C.
Give the expression for the total loss per unit of time Ф as a function of
the internal temperature difference Ti (t) and the external temperature Te .
b) Resolve this equation and give the expressions for в()) and T i ())
Numerical data:
Solution.—
R _ e 1 1
1 = \ + h + ~he [4.107]
S=184m2 [4.109]
T-T
Ф1 = S^----e- *1.4, [4.110]
where the coefficient 1,4 takes into account the losses through the roof.
That is:
Ф1 = 8196 W [4.111]
_ e 1 1 e2
1 = 1Q + h + h, +1 [4.112a]
T-T
Ф1 = S^---- — *1.4 [4.113]
1 Rth 2
Ф1 = 956.6 W [4.114]
t=0;Ti=Ti0=19°C: [4.116]
a) Introducing 0( t)
with:
[4.118]
Pair cV VOLRth 2
Quasi-stationary Model 131
Ln в = - at + LnC [4.120]
в = C exp (-at)
at
t= 0 в0 =C =14 [4.121]
Then,
в = 3° C [4.123]
We find:
a=1.63.10-4 [4.124]
-at=L^tT
[4.125]
T -T
i0 -
1 T-T
tS = aLnTi (t)-Te [4.126]
132 Heat Transfer 1
114
ts = Ln [4.127]
a 3
A DIY enthusiast has built an insulated box to protect their food for their
picnic.
The box has been filled in a house where the atmospheric temperature
was Tint0 =16°C.
This temperature is maintained in the box until it is taken out for the
picnic; it is forgotten outside during a family photo session; the atmospheric
temperature is Te =25°C. The external and internal convection coefficients
are he = 5W.m-2.K-1 and hi = 5W.m-2.K -1 respectively. In addition, no
thermal flow will pass through the base during the problem.
3) Since heat enters the box, the internal temperature Tint (t) increases.
What is the expression for the total instantaneous flow entering the box?
4) Writing the thermal balance for the air inside the box. Deduce the
differential equation from this, which gives the evolution of Tint (t).
Quasi-stationary Model 133
5) We consider that certain foods are in danger above 23°C. What is the
maximum authorized time for the photography session? Will the DIY
enthusiast be pleased by this first test?
6) The DIY enthusiast then uses their coolbox: they introduce a liter of
ice into the box.
- the internal temperature of the air in the box becomes Tint = 7°C ;
- this temperature is maintained as long as all the ice has not melted;
- all of the heat entering the box by conduction is used to melt the ice,
whose fusion temperature is L .
b) For how long can the temperature remain at Tint =7°C in the box?
We give:
Solution.—
The area which the leaks are spread out over is:
S=2*(0.5+0.3)*0.4+0.5*0.3=0.79m2 [4.130]
AT 9
ф = So = S---- = 0.79*------- = 3.648 W [4.133]
R Rh 1.949
ф=ST-u t)
[4.134]
Rth
Quasi-stationary Model 135
For a time dt, the flow of heat entering the box Ф dt leads to heating dTt
such that:
Ф dt = + pCVVOLdTi [4.135]
The plus sign expresses that the incoming heat leads to the heating of the
room.
PairCyVOJdTn- = - S — [4.136]
air V OL dt Rth
d—
= -a— [4.138]
dt
with
S
a = p CVR [4-139]
Pair V OL th
a gives
d0
= -a0 [4.142]
dt
d0 _
— =-a dt [4.143]
в
dLn0 = -at + LnC [4.144]
0( t) = C exp-at [4.145]
C = в0 [4.146]
0( t ) = 00 exp-at [4.147]
-12
t =-------- LLn— = 207s = 3mn 27s [4.151]
7.26.10-3 9
6) Use in an icebox
The thermal flow passing through the box is calculated for a constant
temperature:
T- T 25 7
Ф = £_j---- inL = 0.79*-------- = 7.296 W [4.152]
Rth 1.949
. = mL
ф t,fusion [4.154]
In other words:
This time is obviously more acceptable, which easily leaves time for
transportation and photographs!
Figure 4.2. The icebox, the ancestor of the fridge. For a color
version of this figure, see www.iste.co.uk/ledoux/heat1.zip
2) Taking this value of the flow, calculate a time tf for the fusion of the
ice cube.
3) Will this time be the observed fusion time? If not, will the exact result
be higher or lower than that seen here? We will not provide an answer by
calculation, but instead through reasoning.
Do you think that the order of magnitude found validates the concept of
an icebox?
Quasi-stationary Model 139
Solution.—
Ф tf = mLf [4.158]
where m = pa0 is the mass of the cube and Lf is the fusion heat of the ice.
This gives
pa03 Lf
tff =------- [4.159]
Ф -
917*(0.2)3*333.55.103
t, = ------ i---- ----------------- = 12.24.104 = 34 h [4.160]
f 20
3) This time cannot be exact. Effectively, the exchange surface between
the air and the ice reduces as the ice melts. The fusion time will therefore be
greater than the time found.
We can attempt to calculate this time; we will take the unknown a(t).
For a time dt, the flow of heat entering the ice cube Ф dt leads to a fUsion
of the mass dm
Ф dt = - dmLf [4.161]
dm is connected to da by:
m=pa3 [4.162]
140 Heat Transfer 1
dm = 2 pa2 da [4.163]
Ф dt = - 2 pa 2 daLf [4.164]
da _ -Ф
[4.165]
dt 2pa2 Lf
with
In other words:
And again
da
[4.168]
dt
with:
a (t) = a0 - at [4.171]
Quasi-stationary Model 141
tf =2.45.105s=67h58mn [4.173]
Numerical data:
Densities
Thermal conductibilities
EXAMPLE 4.8.- 1812. Retreat from Russia where “war is always cold”!
The walls and the roof of an isba are built from wooden logs. The walls
and roof here will be assimilated to homogeneous walls with a thickness of
e= 15 cm. The geometry and dimensions of this construction are given in the
figure below. From the point of view of thermal losses, we will assimilate
the doors and windows to walls. We will consider that no flow will escape
through the ground.
1) Calculate the flow Ф of thermal losses through the walls and the roof
of the isba. Deduce from this, the mass of wood m that must be burned each
hour to maintain the temperature in the room.
Solution.—
1) We must first calculate the overall thermal resistance Rth of the wall.
The wall is simply made of thick wood. The thermal resistance of each
1
[4.174]
h
—= 15.10-2 m
л= 0.15 W.m-2.K-1
The outgoing flow of heat for a surface area S of the isba will be
calculated by:
T-T
Ф = Sisba <P = STRTT [4.177]
Ti- T-
?= iR e [4.178]
Rth
Sfac =2.5*10
. =25m2 [4.179]
Isba =fac
S S +S
fond lat )+S
+(2*S roof
[4.184]
146 Heat Transfer 1
E=4802.5*3600W=1.73.107*
J [4.188]
1350
tS =------ = 1098 h [4.191]
S 1.23
in other words:
At this moment, the temperature of the Vodka in the bottle is Ti0 =12°C.
The Vodka, which contains 60% water is likely to freeze. Its volume
then increases, and the full bottle that contains it explodes.
Figure 4.4. The bottle of Vodka and its model. For a color
version of this figure, see www.iste.co.uk/ledoux/heat1.zip
1 We therefore neglect the losses through the neck and the cork here.
148 Heat Transfer 1
We will consider that the external wall of the bottle is at the external
atmospheric temperature.
2) Solving the differential equation. How long after the grenadier goes
out will the bottle explode?
How much longer will the grenadier who has stayed behind be able to
stay in the isba? Is his idea a good one?
Solution.—
a) We calculate the thermal resistance Rth of the walls of the bottle, via
the formula that is applicable to a homogeneous plane wall.
e= 5.10-3 m
2 = 1.05 W.mK-
Quasi-stationary Model 149
The internal surface of the bottle through which losses pass is composed:
S=4SC+SF=8.64.10-2m2 . [4.196]
Presuming a stationary process, the flow of heat Ф that exits the bottle
for a value Tt (t), for an area S = 40 m2 is calculated by:
T-T
p = Sp = S-e------- i [4.197]
b) For a time interval dt, the flow of heat exiting the bottle Ф dt leads to
heating dTi such that:
Ф dt = -pcVOLdTi [4.198]
The minus sign expresses that the outgoing heat leads to the cooling of
the Vodka.
в = Tt (t)- Te [4.199]
P-rCVo^T = - S в [4.200]
dt R
150 Heat Transfer 1
dT =- d (T — T)■ — de
[4.201]
dt dt dt
de
— -a0 [4.202]
dt
with
S
la [4.203]
pcVOLRth
a gives
a= 8.64.10-2 =3.54.10-3s-1
[4.205]
916*3500*1.6.10-3*4.76.10-3
de a
— -a0 [4.207]
dt
de „
— —-a dt [4.208]
e
dLnd — -at + LnC [4.209]
Quasi-stationary Model 151
0( t) = C exp-at [4.10]
6( t ) = 60 exp-at [4.212]
6=6
G0 exp-atG [4.214]
In other words:
tG = -1 Ln6 [4.215]
G a60
-113
tG =-------------------- L Ln [4.216]
G 3,538.10-3 42
EV =50*916*1.6.10-3*11.9.106=8.7.108J [4.219]
tL=50h26mn [4.221]
This is all right. He could have kept one aside to warm himself up in the
snow.
1) Calculate the flow of heat between the tube and the external
atmosphere for a tube length of 1 meter.
cV =716J.kg-1 [4.222b]
Solution.-
Ln— Ln—
Rh =—- + — = —— + — [4.223]
2nA 2nh 2nAS 2nh
We obtain:
The difference in temperature between the two sides of the sleeve is:
A T = TBe
- T = 70 - 21 = 49°C [4.225]
154 Heat Transfer 1
In other words, the thermal flow evacuated to the exterior per meter of
bar:
TT
Ф, = L----T [4.226]
Rh 1
with L= 1 m [4.227]
49
Ф=-------= 105.6 W [4.228]
0.464
We then obtain
Ф= 121.6 W [4.230]
We also find
Ф= 119.6 W [4.232]
We see that:
- from 12 cm to 26.3 cm, the diameter increases and the losses increase;
- from 26.3 cm to 36 cm, the diameter increases and the losses reduce.
Quasi-stationary Model 155
2 54
rc = - = — = 0.131 m = 13.1 cm [4.233]
Ch8
Therefore:
- for de < rc = 26.2cm (from 12cm to 26.3 cm), the diameter increases
and the losses increase;
- for de>rc=26.2cm (from 26.3cm to 36 cm ), the diameter increases
and the losses are reduced.
Part two. Steel bar and glass sleeve; evolution over time
We will look again at the steel bar from the previous example,
surrounded by the glass sleeve with diameter de2=2re2=26.3cm. The
external convection coefficient is still h= 8 Wm-2K-1.
In addition, we give:
We will neglect the thermal losses through the extremities of the bar and
the glass sleeve.
1) Write the relationship between the flow exiting the bar per meter of
length, Ф, T— (t) and Te
4) After a certain time period t0,5 , will the temperature of the bar have
reduced by half? Give this time in hours, minutes and seconds.
Solution.—
TT
ф= TB—— [4.234]
We had calculated:
dT T-T
PSC—B = —B------ e- [4.236]
dt Rth
t=0;TB=TBi=70°C [4.237]
The sign “ - ” in the equation indicates that the flow of heat (positive)
contributes to the reduction in the temperature; the derivative is therefore a
negative term.
d00
Psc OL dt ~ Rth [4.239]
(
VOL = nri2 = п 2.5.10-2 )2 = 1.96.10-3 m3 [4.240]
158 Heat Transfer 1
or:
do dt
— =---------- = - adt [4.241]
0 psCRth
1
with a = [4.242]
psc VOL Rth
О= Cexp(-at) [4.246]
at
TB(t)-T=(TBi-T)exp(-at) [4.248]
4) We calculate
in other words, TBi - T from TBi -T=49°C to TBi -T=35 -21 = 14°C
14 = 49exp(-atS) [4.250]
or even:
14 = 49exp(-atS) [4.251]
-114 1 49
tS =— Ln— = Ln— [4.252]
S a 49 3.47.10-4 14
tS =3610s=1hr10s [4.253]
a) Calculate the thermal resistance per meter of the glass cylinder, Rth1 ,
and deduce from this the flow Ф verre 1 that escapes from the reagents at the
initial instant of the reaction.
a) Give the expression2 for the resultant thermal resistance Rth2 for this
new reactor. Calculate its numerical value.
b) Establish that, taking into account this new value Rth2 , the flow
Фisol (t) exiting this new reactor configuration can be written as
Ф = в( T - Te); в is a coefficient whose numerical value will be given.
e) By how much will the temperature of the reagents reduce after half an
hour? Who is better, the physicist or the chemist?
2 Since the thickness is not known, this question will be treated literally up to e.
Quasi-stationary Model 161
Numerical data:
di = 5 cm
e = 6 cm
d
e= 3cm
L=1.5m
Tri =100°C
Te =19C
h= 5W.m-2.K-1
m= 15 kg
Solution.—
Rth [4.254]
2 nA
r 6
hr 1 Ln 5 1
r
1 " 2nA + 2nreh 2n3 2n 3.10-2, .5
= +
T-T 81
Ф = L^---- e- - = 1.5.-------- = 111.81Ж [4.258]
Rth 1 1.0867
Q = mc =15.4000=60,000J [4.259]
Lnr r+e 3
Ln—---- Ln Ln 6
r 1 =2.5 +---------- 3----- + 1
Rth1 = ri + +' e
2nA pyrex 2nALR 2nre h 2n*1.13 2n*0.045 2n*6.10-2*5
re'=6.10-2m [4.262]
b) The flow Фisol (t) exiting this new configuration of reactor can be
written as:
TT TT
Ф= L—--- e- = 1.5.—---- e- = 0.499(T -T) [4.264]
Rth1 3.0077
в = 0.499 [4.266]
c) Heat balance: by writing that the heat exiting per unit of time cools the
reagents, we draw up the differential equation for T(t)
Ф = в (T - T) [4.268]
with
в = 0.499 [4.269]
в = 0.499 [4.273]
164 Heat Transfer 1
0 = ( T - t. ) [4.274b]
dO
= -a0 [4.274c]
dt
в 0.499
a=в =---------- = 8.31.10-6 [4.275]
mc 15*4000
d0
= -adt [4.277]
0
Ln0= -at +Ln C [4.278]
0=Cexp-at [4.279]
e) By how much will the temperature of the reagents reduce after half an
hour.
00 =81° [4.282]
a= 8.31.10-6 [4.283]
0=81exp-8.31.10-6.1,800=0.985.81=79.98 [4.285]
A T = 1.2 ° [4.286]
On a hot summer day, a bottle of beer has been left next to a window.
The temperature of the beverage has reached T0 =25°C.
To drink a cold beer, at TF 1 =6°C, a tourist places the bottle in the main
body of the fridge, which is at a temperature of Te1 =4°C. We will consider
that the natural convection imposes a uniform convection coefficient on the
external surface of the bottle, he = 5W.m-2.K-1. We will also consider that
the temperature Ti (x) of the beer is uniform at all times in the bottle and
that, as a consequence, the cooling process is quasi-stationary.
In addition, we point out that the beer freezes at TG =-2°C and that, in
this case, the beer dilates and the bottle explodes.
a) Calculate the thermal resistance Rth per meter of length of the glass in
the bottle.
3 Here, we will neglect the thermal effects of the neck and the bottle cap.
166 Heat Transfer 1
b) Calculate the thermal flow Ф0 coming out of the bottle at the initial
moment in time; in other words, when the bottle is placed in the fridge.
c) How long will the tourist have to wait before drinking their cold
beer?
3) Since this time is too long, the tourist’s wife, who is not really a
physicist, places a second bottle of beer, also at T0 = 25°C in the freezer
compartment, where the temperature is Te2 =-15°C.
She thinks that the beer should be left for a shorter period in the freezer
and leaves it for a time t2 = 1 hour
Additional data:
Solution.—
rr
Lnre Lnr
__ T_ =__ [4.287]
2пЛ 2пЛ
ri = 3.10-2 m
re = 3.5.10-2 m
Di = 6.10-2 m
De = 7.10-2 m
Л=1.05 W.m-2.K
168 Heat Transfer 1
r
Ln Ln7
---- r- =------- 6— = 2.34.10-2 W.m-2.K [4.288]
2nA 2*n*1.05
r
Ln r
tot = --- rr +----1 = 2.34.10-2 +---------- 1----- ------
t 2nA 2nrh 2*n*3.3.10-2 *20
The flow of heat Ф exiting the bottle for an initial value Tt 25°C is
calculated for a cylinder length H= 26 cm, in other words:
T-T 25-4
Ф0 = H—---- e = 0.26*------- = 20.7W [4.290]
0 Rh 0.264
V(t)= hT Ro Te ( [4.291]
a) For a time dt , the heat flow exiting the bottle Ф dt leads to heating
dTi such that:
Ф dt = - pcVOLdTi [4.292]
The minus sign expresses that the outgoing heat leads to the cooling of
the liquid.
в = T. (t)- Te [4.293]
Quasi-stationary Model 169
de
=-ae [4.296]
dt
with
H
a =----------- [4.297]
PcVoR
(
VOL = n* 6.10 ) *26.10-2 = 7.35.10-4m3 [4.298]
a gives
t=0;e=e0=Ti0-Te=25-4=21°C [4.300]
de
=-ae [4.301]
dt
170 Heat Transfer 1
dO = -a ,
„ dt [4.302]
O
dLnO = -at + LnC [4.303]
O(t)=Cexp-at [4.304]
In other words:
t1 =-! Ln0-
[4.309]
1 aO0
-12
[4.310]
1 3.2.10-4 21
Of =12.64 [4.314]
-131
t =------- l Ln— = 3512 s = 58.5 mn [4.316]
exp 3.2.10-4 40
We are not far off, it is a shame. Disaster again, but only by a small
margin!
c) Unlucky
ri = 3.10-2 m
re = 3.1.10-2 m
Di = 6.10-2 m
De = 6.2.10-2 m
1 = 204 W.mK—
Lnr-e 6 6.2
Rh = 1= П ^2 = 2.56.10-5 W. m ~2. K4
—
[4.317]
2n1 2*n*204
172 Heat Transfer 1
r
Ln r
== rr + 1---- = 2.56.10-5 +---------- 1-----------
tot 2пА 2nrh 2* п*3.1.102*20
=2.56.10-5+0.257=0.257W.m-2.K-1 [4.318]
TT 25 4
Ф0 = Н-i------- e- = 0.3*-------- = 24.5 W [4.319]
0 R 0.257
V = п* (6.10
OL 4
) *0.3 = 8.48.10-4 m3 [4.320]
-:1---- Ln-2
t =------
[4.322]
1 3.29.10-4 21
COMMENT.-
Put more simply, we admit that this value of he is an “ad hoc” value,
which involves taking into account two phenomena that intervene in two
areas of the box or the bottle:
- natural convection through the top;
- direct contact of the bottle with the shelf, which necessarily imposes the
temperature Te =4°C on part of the external wall of the recipient and which
automatically eliminates all convection effects on this zone.
To heat the yurt, the head of the family has lit a wood fire that consumes
mkg of wood per hour. A kilogram of wood provides an energy Q by
combustion. We largely consider that only 50% of heat produced passes
through the felt by conduction. The remaining 50% exits directly with the
smoke and has no incidence on the temperature in the yurt.
We consider that when the fire is lit, we are in a stationary regime. The
internal temperature of the yurt Ti is constant. The outside temperature is Te .
174 Heat Transfer 1
4 This approximation is inspired by the construction of a “standard” yurt that is often carried
out using five plane vertical wall panels.
Quasi-stationary Model 175
2) Given the answer to the previous question, show that the correct
expression for the instantaneous flow exiting the yurt as a function of the
difference between the internal and external temperature takes the form
Ф = в(Tt - T) [4.324]
4) With no stock left, the head of the family went to get some wood in his
old Land Rover. He will take a whole day to do this. Will the family suffer
from the cold?
c) We consider that below 5°C, it is not possible to live in the yurt. How
long will it take for this temperature to be reached?
Will the head of the family come back in time to keep the family
comfortable?
A Mongolian sports club rented a yurt that is identical to the one seen in
previous questions and installed 10 home-trainers there. The club does not
176 Heat Transfer 1
have the means to heat the yurt and counts on the heat given off by the
cyclists to maintain an acceptable temperature in the tent.
Remark.-
Additional data:
D=6m
H=2.3m
e= 10mm
Te =-8°C
hi = 5W.m-2.K-1
he = 16W.m-2.K-1
Q= 14MJ.kg-1
Solution.—
a) Plane wall
„ nDD- .
5= + nDH = 28.27 + 43.35 = 71.62 m 2 [4.326]
AT 22
Ф = 5---- = 71.62--------= 2875 W [4.327]
Rth 0.548
b) Cylindrical wall
Ln 1 De-1 T (6.02^
Ln 1 1
Rth = +
1 Di J
2nXF
i 1
2nrihi
+
1
=
1 6 J
2nrehe 2n* 0.035
+
1 —
2n*3*5
+
1
2n*3.01*16
=1.5110-2+1.0610-2+3.3110-3=2.9.10-2W-1m2°K [4.328]
nD2
SRoof = = 28.27 m2 [4.329]
4
Д T = 22 ° [4.330]
ДT
Ф = 28.27 +2.3-------- - =1134.9+1744.8=2897.9W [4.331]
0.548 2.9.10-2
178 Heat Transfer 1
The form that takes account of the cylindrical geometry is more suitable.
Therefore:
AT
Ф = 28.27 +2.3---------r =в ( T - T) [4.333]
0.548 2.9.10-2
в=130.9 [4.334]
2898*3600*24=2.504.108J. [4.335]
5.077.108
It is therefore necessary to burn: = 35.77kg of wood every 24
14.106
hours.
Quasi-stationary Model 179
Ф = в (T — Te) [4.337]
в = 130.9 [4.338]
DD-
Vol = n—H = 65.03m3 [4.339]
4
0=(Ti-Te) [4.341]
d0
= -a0 [4.342]
dt
d0
=-a0 [4.345]
dt
180 Heat Transfer 1
d0 = -a
— „ л
dt [4.347]
о
Ln0 = - at + LnC [4.348]
0 = C exp - at [4.349]
13 =22exp-at [4.353]
113 1
t =—Ln _ =--------- - (-0.526) [4.354]
a 22 2.156.10-3
t= 244.01s [4.355]
b) It is a good thing that extra heating was not installed inside the yurt.
But after their efforts, the cyclists should not stay for too long inside the yurt
because it will cool down.
5
Non-stationary Conduction
дT
(
р c---- = div A grad T
дt
) [5.1]
In addition to this, generally, the physical properties of the material can vary.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
184 Heat Transfer 1
Remark.- Again we note that a parameter derived from the previous chapters
has appeared, and that this will have a pivotal role in all non-stationary
conduction problems, thermal diffusivity a, defined by:
a =— [5.2]
where T1 , T2 and T3 are temperatures, and which are not all necessarily
different, selected depending on the problem in question, and, in particular,
on its limit conditions.
dT d2 T
= a------ [5.4]
d t-----d x2
T(0,t)=Te [5.5]
T(x,0)=Ti [5.6]
REMARK.- We should recall that the heat equation results from a thermal
balance and the limit conditions from the specific conditions of an
experiment.
d0 d2 0
---- = a----- [5.9]
дt д X2
0(0,t)=0 [5.10]
0(x,0)=1 [5.11]
In other words,
REMARK.- This variable does not appear intuitively. We can look for its
expression in the form 0(x, t) = axm tn. Here, we do not give the
calculation and advise curious readers to work through it themselves.
The heat equation, which initially has two variables, will be transformed
into an equation with a single variable. We will therefore go from a written
expression in the form д to a written expression in d.
д0(x, t) _ d 0(^) д ^
[5.15]
дt ~ d^ дt
Non-stationary Conduction 187
Э0( x, t) _ d0(£)d d
[5.16]
dx d2 dx
With:
d_ _д[ x _ 1
[5.18]
d x dx 4 4at 41at
Э0(x, t) -d d 0(d)
[5.19]
dt ~ 21 d
d0( x, t) 1 d0(d)
[5.20]
dx 4 a at dd
1 1 d d0(d)
[5.21]
4 а at 4 aat d d dd
- d d0(d) _ 1 d 20(d)
[5.23]
21 dd 4 at d 222
188 Heat Transfer 1
d d°^
° = d{ [5.24]
°._ d2°(Я
° = d? [5.25]
or even:
integrates into:
Ln °1 = - ? + LnC1 [5.29]
°1 = C1 exp- ^2 [5.30]
° = C1 j^ exp - u2 du + C2 [5.31]
x
i = ,----- [5.34]
4 at
0(0) = 0 [5.35]
0(7 = 1 [5.36]
Hence:
Q
0=C1jO exp-u2du+C2 [5.37]
0
0 = C1 j^ exp- u2 du + C2 ^ C2 = 0 [5.38]
1 = C1 exp- u 2 du + C2 [5.39]
Hence:
( x \
T(x,t)=(Ti-Te)erf +Te [5.42]
This problem calls into question a wall that has a plane face. The physical
properties of the material are presumed invariable, that is, following the
usual notation of this book, p, c, 2.
Non-stationary Conduction 191
t=0
dT d2 T
---- = a------ [5.43]
dt d x2
X
a =----- [5.44]
p cP
The limit conditions express that the temperature away from the surface
of the wall is not affected by the modulation; the second expresses the
temperature condition imposed.
In other words:
T (~, t ) =T [5.45]
dT T [5.46]
- ^dx = ^°’
Here, we give the solution, which is, this time, in dimensional form:
у/
T ■ x
2 ^0 at1 t Zierfc
(x, 2t\) - T2 =------ x [5.47]
X ya4at7
The functions erf %, erf % and ierfc% are tabulated in Appendix 5. These
three functions are represented in Figure 5.3 as well as the Appendix:
This problem involves a wall with a plane face. The physical properties
of the material are presumed invariable, that is, following the usual notations
of this book, p, c, Л.
dT dT
---- = a----- [5.50]
dt d x2
A
a =----- [5.51]
p cp
The limit conditions express that the temperature away from the surface
of the wall is not affected by modulation; the second expresses the imposed
temperature condition.
In other words:
0=T-Ti [5.54]
d0 d2 0
— = a----- [5.55]
дt д X2
associated with
0(~, t ) = 0 [5.56]
We will use a complex formulation to solve this: we look for the reduced
temperature in a complex form; the solution will be the real value of this
complex function.
Therefore, we write:
0(x,t)=f(x)eiat [5.59]
196 Heat Transfer 1
в(,х, t) = 0 [5.60]
The solution will be the real part of the result that is found.
— = iafeiat [5.62]
dt J
( f"+l^L) ea = о
[5.65]
Leading to:
f"+iaL=о [5.66]
a
f (x)=Aexp [5.67]
REMARK.- iji is not necessarily familiar to all readers. We will easily find
its value by writing it in the form of a complex number where a and b are
real numbers to be found:
i/i = a + i b [5.68]
)
i = (y/i 2 = (a + ib)2 = a2 + b2 + 2i ab [5.69]
by identification:
a = b = —1=
22
Therefore:
г 1+i
i/i =---- [5.71]
2
at x = 0: [5.72]
f= 0 ^ A + B = T0 [5.73]
at infinity:
x ^^ ^ f^ 0 [5.74]
irn i® 1 + i I® 1 + i I®
f (x) = A exp x + B exp — x = A exp—— J —x + B exp---- =-. —x
a a a a
198 Heat Transfer 1
> ( ггт X ^
= A exp x +В exp 2a
a x 7 exp - i. .-x 7 [5.75]
V 7 V V
Therefore:
В = To [5.76]
Finally:
a .
e( x, t ) = f (x) eat = To exp - '— x expi at [5.77]
a
(
0( x, t) = T
jo at
V
(
= To exp - cos at - [5.79]
( a л
cos at - iaxJ [5.8o]
V
Non-stationary Conduction 199
( ( — 1
T(x,t)= Ti +T0 exp cos —t — 2 aa x [5.81]
к к 7
We see that the temperature at a point of the material is the sum of the
initial constant temperature and of a modulation.
T0 exp [5.82]
Ф= —x [5.83]
2 2a
These two materials are denoted by 1 and 2 ; they are initially at the
respective initial temperatures T1i and T2i. Their respective physical
properties, which are supposedly invariable, are p 1, cP 1, X 1 and p2, cP2, X2
when written using the usual notations of this book.
t=0
Therefore, contrary to the other problems seen here, the solution lies in
two temperature functions.
These two systems need to be written down for each environment, where
each differential equation is only valid for half of the space.
In environment 1 ,
д д T1 д T1
x<0 —= a [5.85]
дt ’ дx
In environment 2 ,
д дT
- д T2
x > 0 —-■ a. [5.86]
дt дx
A
a =----- [5.87]
p cp
We will define two temperature variables, which are not reduced here:
01 = T1 - T1 i and 02 = T2 - T2 i [5.88]
In environment 1 ,
д01 д01
x<0 —b = a —1 [5.89]
дt 1 дx
202 Heat Transfer 1
In environment 2,
д 02 д 02
x>0 —- = а^ —- [5.90]
дt 2 дx
Л
а =------ [5.91]
p cP
In environment 1 ,
x<0
01 (x,0)=0 [5.92]
01 (-ж, t ) = 0
In environment 2,
x>0
02 (x,0)=0 [5.93]
0i (+~, t ) = 0 [5.94]
x=0 Л^-T- = Л2 —
[5.95]
дx дx
d 61 d 62
x=0 4^ = ^23 [5.97]
dx dx
Here, we will not give the details of this quite complex solution, the
perspective being to provide an applicable relationship.
the variable
[5.99]
2jOt 4^ at
where
Here, we note the use of an absolute value. Indeed, in this environment, all
axes are negative.
204 Heat Transfer 1
T2 (x, t)- T 2 E x x
--------------- =---------- erfc .__ - [5.103]
T1 — T2 E + E2 4a^t
E = Xpc [5.104]
In other words:
E1 = 4 P1 c1 [5.105a]
E2 = A2 p2 c2 [5.105b]
How long will it take for the center of the concrete block to reach the
temperature of 40°C?
Non-stationary Conduction 205
Solution.—
Ti =20°C. [5.106]
g = T ( X. t)-Te [5.107]
T, - T
where
£ = 0.477 [5.113]
This value must be found in the middle of the block, in other words, at
x = 0.1m from the heating face.
xx
£( x, t )= =------ [5.114]
2^1at 4 a at
-1 f x Y
4 a ^%) [5.115]
2
1 f 0.1 Y
t= , ,L 1 [5.116]
4*6.93.10-7 ^0.477)
t=15855s=4h24mn15s [5.117]
We plunge a large cube of stainless steel into a bath of molten salts. The
initial temperature of the metal is 25°C. The temperature of the salt bath is
500°C. This temperature is imposed on the cube faces. Given the size of the
Non-stationary Conduction 207
cube, we estimate that the transfer of non-stationary heat for each face is
single-dimensional (the flow of heat is locally perpendicular to each face).
How long will it take for the temperature to reach 400°C at 2 cm from the
surface of each face of the cube?
Solution.—
xx
= = (x x, t ) = [5.119]
where erf (g) = -^^e u2 du is the error function, for which a table is
provided in Appendix 5.
400 - 500
о= = 0.21 [5.121]
25 - 500
208 Heat Transfer 1
% = 0.188 [5.122]
This value must be found at x = 2.10-2 m from the wall in contact with
the molten salt.
2
1 (x^
t =—1 - 1 [5.124]
4a L %J
2
-0^211 [5.125]
t= 1 .(
4*5.53.10-6 ^ 0.188 J
t=512s=8mn3s [5.126]
How long will it take for the temperature to fall to 5°C at a distance of
5 cm from the internal face of the concrete wall? Express this time in hours
and minutes.
Non-stationary Conduction 209
Solution.—
T (x,t)-T
в = -(x-)— e [5.127]
T-T
ie
xx
£ = £(x, t ) = [5.128]
2 f£ ,
where erf (g) = ^= £ e du is the error function, for which a table is
provided in Appendix 5.
5-0
в= = 0.333 [5.130]
15 - 0
This value must be found at x =20 -5 =15cm = 15.10-2 m from the face
of the wall in contact with the molten salt.
210 Heat Transfer 1
xx
£( x, t ) = —= = ^= [5.131]
lyj a t 4 a at
-1 f x Y
4 a ^%) [5.132]
1 f 0.15 У
t= , 7। । [5.133]
4*8.38.10-7 ^0.305)
We will also presume that over the course of the entire transfer, the water
will impose a temperature of 95°C onto each face of the cube.
Density:
p = 7816 kg.m3
Thermal conductibility:
c= 460 J.kg-1
Non-stationary Conduction 211
1) Evaluate in the context of this approximation how long it will take for
the temperature to reach 55°C at x = 4cm from each of the faces.
2) Evaluate by the same process (which then becomes quite large) how
long it will take for the temperature to reach 80°C at the center of the cube.
with
T (x,t)-T
в = -(-,-) [5.136]
T -T
ie
Prior to all problems of this type, we will calculate the thermal diffusivity
of the steel
we fix
the time, which is the value that we are looking for, will be deduced from the
value of giving в as prescribed.
In other words:
в = 55 95 = 0.54
[5.139]
21 - 95
212 Heat Transfer 1
^ 2 --------------- в"
t = ( ------2 [5.142]
^ 0.523 J 4*4.53.10-6
t=322.8s=5mn23s [5.143]
x=7cm=7.10-2 m [5.144]
0 = 80 95 = 0.2027 [5.145]
21 - 95
^ 2 -----------------
t = (-- ------2 [5.147]
I 0.182 J 4*4.53.10-6
t=8164s=2hr16mn4s [5.148]
Non-stationary Conduction 213
2 = 0.46 W. m-1.kg 1
Density:
p = 929.kg. m-i
c=1830J.kg-1
214 Heat Transfer 1
1) What is the flow density ф0 imposed by the iron at the top of the pile
of laundry?
2) Give the expression as a function of the time taken by the top of the
pile of laundry to heat up: T(x=0, t)-Ti .
How long will it take for the laundry at the top of the pile to begin to
melt?
SOLUTION.- 1) The density of flow ф0 imposed by the iron on the top of the
pile of laundry is obtained by dividing the flow P= 1000 W by the heating
surface S = 100 cm2 = 10-2 m2:
1000
P0 = =105W m-2 [5.149]
10-2
2) Expression as a function of the time taken for the top of the pile of
laundry to heat up:
xx
i = —/= = “/= [5.152]
2 at 4 a at
Non-stationary Conduction 215
In this case,
S=-0= = 0 [5.154]
2yJat
We calculate
a = — = 2.7.10 [5.157]
We will reach the temperature T(0,t)=220 °C at the top of the pile for
the time tF such that:
tF = 2.9 s [5.161]
Then:
5 10-2
c = —. • - = 0 878 [5.162]
2V2.7.10-7 *30
( )
в 5.10-2,30 = ^.10-J2.7.10-7 *30 * 0.073 [5.164b]
0.46
An egg, when it is laid, is made up of a shell, egg white and yolk. The
white of the egg is a part of the egg, made essentially of an albumin,
ovalbumin, which protects the yolk or the zygote.
In the egg yolk, the proteins also coagulate, but at a slightly higher
temperature of 65°C.
Non-stationary Conduction 217
1) A central cylindrical zone of radius 0.5 cm, known as the egg yolk.
2) An annular cylindrical zone surrounding the first, of thickness 5 mm,
called the egg white.
3) An extremely thin shell that will have no effect on the thermal
phenomena.
Water
We will consider that the white and the yolk of the egg have the same
thermal properties, which are very close to those of albumin. We will
therefore consider that the egg white and the yolk constitute a single
environment, which is also strictly immobile.
1) Initial calculation
A very thin plane wall, for which we will consider thermal resistance to
be zero, separates two spaces containing:
- On the left, water at 100°C, which will only have the effect of
maintaining this temperature at the separation for the duration of the
problem.
218 Heat Transfer 1
Water Albumin
Te = 100C Ti = 20C
Figure 5.8. Initial calculation
2 = 0.6 Wm-2 K -2
Density:
p = 1000 kg. m -3
c= 4180 Jkg-1K-1.
Remark.- Well-informed readers will certainly not miss the fact that these
are the thermal properties of water. It is not always easy to find precise data
in terms of the physics and chemistry of food items. Egg white contains 88%
water. Among the other constituents of albumin, we find 10.6% of globular
proteins, the main one being known as ovalbumin. For this simple model, we
have chosen data related to water, which are easier to access.
2) Cooking
We will admit that if the part at the center of the egg white has reached a
temperature of 57°C, the white will be coagulated and the yolk will remain
soft. We will have a soft-boiled egg.
a) Soft-boiled egg
b) Hard-boiled egg
In the same way, estimate the order of magnitude of the time required to
cook a hard-boiled egg.
220 Heat Transfer 1
c) Are these results coherent with your personal experience (if any)?
Solution.—
1) Initial calculation
—
a =--- = 1.43.10-7 m25-1 [5.167a]
Pc
a) We are looking for the time t required so that equal T1 reaches 57°C
in the albumin at a distance e= 0.5cm from the separation surface with
water at 100°C.
T-T
e=------ - [5.168]
Ti- Te
with
Ti = 20°C [5.169]
Te =100°C [5.170]
Non-stationary Conduction 221
0 = erf Q [5.171]
with:
x
Q = ^= [5.172]
4 a at
We write
T(x,t)=57°C [5.173]
then:
0= 0.537 [5.174]
Q= 0.519 [5.175]
Q= 5.10-3 = 6.61
[5.176]
^4*1.43.10-7 t tit
/= = 12.7 [5.177]
t=162s=2.7mn [5.178]
b) We now look for the moment when the temperature T2 reaches 65°C in
the albumin at a distance of e= 1 cm from the separation surface with the
water at 100°C.
T-T
0 =------- [5.179]
Ti- Te
222 Heat Transfer 1
i =20°C
T [5.180]
Te =100°C [5.181]
в = erf Q [5.182]
x
i = ,----- [5.183]
4 a at
T(x,t)=57°C [5.184]
в = 0.437 ; Q = 0.413
^= 10-2 = 13.22
[5.185]
^4*1.43.10-7 t t[t
4= = 32.01 [5.186]
t=1026s=17mn [5.187]
2) Cooking an egg
a) Soft-boiled egg
The preliminary calculations give us the result: using the modeling in the
albumin, looking for a temperature of 57°C at 5 mm from the shell, we
estimate the time taken for coagulation, therefore for cooking, as
2.7 minutes.
b) Hard-boiled egg
In the same way, using the modeling seen above in the albumin, looking
for a temperature of 65°C at 1 cm from the shell, we estimate the
coagulation time, therefore the cooking time, to be 17 minutes.
Over and above the approximation made on the thermal properties of the
white and the yolk, here, we use the results of a plane model in a revolving
system. The geometry selected here is particularly schematic.
Thermal conductivity:
Density:
c= 4000 Jkg-1K-1.
The slice of meat, which has just been taken out of the fridge, is at
Ti = 15°C. We will consider that the meat will be cooked “rare” at a
temperature of TCS = 50°C.
We will also consider that, although the slice has a finite thickness, the
evolution of temperatures within the meat is calculated as it is for a body of
infinite thickness.
224 Heat Transfer 1
We also consider that the physical properties of the meat are modified
very little by the cooking1 (which is a debatable approximation):
1) We are seeking to determine the thickness of the meat cooked after a
specific duration of time.
a) What is the density of the flow ф0 imposed by the plate on the meat?
b) What thickness of meat will be cooked after 2 mn?
2) After these two minutes, what will the temperature TS be at the surface
of the meat placed on the plate? All good cooks know that meat must be
“sealed” at the surface (meaning that the surface layer must be “pyrolyzed”.
Is this the case here?
3) A slice of meat with a thickness of 1.5 cm is cooked by placing it for
two minutes each side on the plancha.
Will the meat be:
- cooked “rare”?
- cooked “medium rare”?
- grilled all the way through?
Solution.-
1 We will admit it even for “pyrolyzed” meat at the surface (effect sought for “sealed” meat).
This approximation is obviously a subject of discussion. We will not discuss this here.
Non-stationary Conduction 225
A 0.47 -7 —2
a =---- =-----------------= 1.09.10 ms [5.191]
pc 1080*4000
xx
i= f~ = l----- [5.192]
2< lat 4 at
At time t=2 mn=120 s, we are looking for at what point the temperature
is equal to TCS =50°C, in other words,
35 * 0.47 2
ierfc (%) =-------------------------------- = = 9.47.10-2 [5.196]
2*2.4.104 * J1.09.10"7 *120
226 Heat Transfer 1
Hence,
4 = 0.782 [5.197]
2) What will the temperature TS be after these two minutes, at the surface
of the meat placed on the plate?
0 = TS - t [5.200]
t = 2 mn = 120 5 2 [5.187]
TS = 233°C [5.207]
Non-stationary Conduction 227
We will have “rare” cooked meat for 2 * 0.565 = 1.13cm and the center
will be very “raw” for 1.5-1.13=0.37cm =3.7 mm, with a well-sealed
surface.
We estimate that the welding is a very thin film of tin that does not act on
the temperature distribution in the copper block with sides a.
The fixing is sufficiently rapid for us to estimate that the upper face of the
cube is rapidly brought to a temperature of T0 =310°C, which it will
retain.
Density:
Thermal conductibility:
X = 386 kg.m-i Wm 1 K 1
c= 383 J.kg-1 .
Non-stationary Conduction 229
So, Ox is the vertical axis, with origin O at the center of the upper face
of the cube.
0
р
2) Tin melts at 232°C. How long after the cube and the plate are put into
contact will the cube B fall?
Solution.—
Defining:
T(x,t)-Te
О= [5.208]
230 Heat Transfer 1
we have:
where:
x [5.—10]
Q = I-----
4 a at
T(x,t)=(Ti-Te)erf(Q)+Te [5.—1—]
2) The tube will fall when the tin has melted. This will intervene at
time tF where the temperature will reach TF =232 °C in the plane that
a -2
separates the two half-cubes, in other words at x = — = 5 cm = 5.10 m .
Л 386 -a -I
[5.—13]
a =----= = 1.126.10 m 2s s
pc 8954*383
in other words:
x 5.10-2
Q (x, t) = —= = 0.244 = [5.216]
2Jat J4*1.126.10-4* t
Non-stationary Conduction 231
2
I 5.10-2 ] 1
t = I--------I ---------------- T [5.217]
^ 0.244 J 4*1.126.10-4
t= 93.23 s [5.218]
For the plate of zinc, we place it on the layer of tin. We presume that the
contact is perfect and that no layer of air can pass between the copper and
the tin.
On the upper part of the copper, we place a soldering iron tip with a
surface that is presumed to be much larger than the soldering zone. The
problem of conduction will thus remain purely single-dimensional.
x=0
Soldering plate
Tin layer
Base
x Axis
Figure 5.12. Diagram of a tin weld
232 Heat Transfer 1
We give:
Density:
р = 8954 kg.m-3
Thermal conductibility:
X = 386 Wm-1K-1
c= 383 J.kg-1
Density:
р= 7144 kg.m-3
Thermal conductibility:
X= 112 Wm-2K -2
c= 384 J.kg-1
a) The tin melts at TF =232°C. What is the minimum time that the iron tip
must be left on the plate of copper so that the weld can be performed?
After time t found in 1a, what will the temperature be in the tin?
What conclusions do you draw from this concerning the validity of our
previous calculation? In order to reason this answer as best possible, it will
be useful for us to be able to calculate the temperature of the upper surface
of the copper plate at time t calculated in 1a.
Solution.—
A 386 -4 -I
a =---- = = 1.125.10 m25 [5.219]
pc 8954*383
T-T
в =----------- e- [5.220]
T - Te
with
Ti =21°C [5.221a]
Te =400°C [5.221b]
в = erf Q [5.222]
with:
x
i = ,----- [5.223]
4 at
234 Heat Transfer 1
We want to reach
We have
0 = 0.443 [5.225]
% = 0.415 [5.226]
t -0.32s [5.229]
A 112
a —---- - =-------------- — 4.08.10-5 m25 [5.230]
Pc 7144*384
T-T
0 ------ e- [5.231]
T - Te
with
T. - 21°C [5.232]
Te - 400°C [5.233]
Non-stationary Conduction 235
0 = erf Q [5.234]
with:
x
Q = ,— [5.235]
4 at
We wish to reach
x =e=5 mm [5.237]
We have
0= 0.443 [5.238]
Q= 0.415 [5.239]
0.415 - 3.10-3 =
[5.240]
74*4.08.10-5t tit
t -0.32s [5.242]
We therefore find
е= 3.77°C [5.247]
which leads to
T =24.77°C [5.248]
which is obviously highly insufficient for a good weld. The welding time has
been significantly underestimated by the first hypothesis, which must
therefore be pushed back in time.
T =29.49°C [5.250]
6
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
238 Heat Transfer 1
(the wings occupy 2.4 cm of 3.14 cm of the outside), the exchange surface of
the eight wings will become 5' = 8*6.6.10-2 *1.10-2 = 5.28.10-3 m2.
Furthermore, we will see that the relation between the two flows defines
the efficiency n of the fin:
Ф
n = Ф0
A wing takes the form of an object, generally metallic, with one end
welded to a body that has a temperature T0 . This object will have a cross
section S that can vary as a function of the distance from one of its cross
sections with respect to the body.
Two main geometries can exist: a linear fin and a circular fin.
The basic problem that we have to resolve each time we study a wing is
as follows.
To resolve this problem, we will divide the fin into infinitely thin
sections, and write a thermal balance for each of these sections:
We will locate the various sections of the fin on the Ox axis that is
perpendicular to the body considered, where the origin O is on the body in
question. At this stage of the equation, the fin is not a parallelepiped, and
therefore S= S(x).
This hypothesis presumes that the lateral thermal gradients are negligible
in comparison with the longitudinal gradients. Physically, this comes down
to assuming that the lateral convective flows are weak in comparison with
the axial conductive flow.
We will write the flow. First, we note that the flows are directed from the
body towards the end of the wing; therefore, the temperature decreases
dT
according to Ox and---- < 0. The flow entering by conduction is given by:
dx
In this expression, p(x)dx is the lateral surface through which the flow
taken out by convection flows.
Ф e =Ф S +Ф Conv [6.4]
dT dT d ( dTV
- AS (x)--- = - AS (x)----- 1----- 1 - AS (x)---- I d x
ax 'dx dx ^ 'dxJ [6 5]
+ [p(x)dx]h(T(x) -T)
After simplification:
d AS (x) dT 1
-hp(x)(T(x)-Te)=0 [6.6]
dx dxJ
242 Heat Transfer 1
We note that in the most general problem, the term in S(x) remains in
the derivative, which can lead to equations that are complex to resolve.
6.1.3.1. Solution
To illustrate the type of solutions found, we consider here the
simplest case, that of a parallelepiped fin. We also provide a few other
solutions in Appendix 6.
First, a point about vocabulary and notation must be specified. We call S
the face of the wing that is parallel to the surface of the body (at the
extremity origin ( x = 0 )); a surface S is fixed to the body to be cooled.
S=le=Cte [6.7]
and
p =2(l+e)=Cte [6.8]
d2 T
(
''.Sdx - hp T (x)— Te = 0 ) [6.9]
where, again:
(
d2T - ^ T (X)— T = 0 ) [6.10]
Fin Theory: Notions and Examples 243
Here, this has led to the introduction of a parameter ю, a classic for these
fin problems. This parameter is defined by:
hp
Ю = —^~ [6.11]
AS
0(x)=T(x)-Te [6.12]
Then:
d0( x) dT (x)
[6.13]
dx dx
dО
------- Ю 0 = 0 [6.14]
d xг
-
eeT=0 [6.16]
x=0;0=T0-Te=00 [6.17]
x=L - A — = 0 ^ do = 0 [6.18]
dx dx
We can also take into account the flow at the extremity. Keeping the
same convection coefficient to simplify, we will have:
x=0 ; 0 = T0 - T = 00 [6.19]
x =L ; - A — = h (T - T) ^ — =-h 0 [6.20]
dx e dx A
d 20
— Q a2 0 = 0 [6.21]
dx 2
is known.
0=Ashax+Bchax [6.23]
This last form will most often be preferred for practical reasons.
Moreover, these two expressions are equivalent, since the sine and cosine
hyperbolics shax and chax are linear forms of the exponentials eax and
e-ax.
NOTE.- Readers who are not familiar with hyperbolic functions can refer to
Appendix 8, which contains a few reminders that are essential for these
functions.
de
— - &■e = о [6.24]
dx
x=0 ; e = Tо - Te = e0 [6.25]
x >° ; e ^ Te - Te = 0 [6.26]
In this specific case, we will use the solution of the sum of two
exponentials:
B = e0 [6.30]
We therefore obtain:
e=e0 e [6.31]
and then:
Ф = -ShdT
dx
(
=-SЛ(To -Te) -me-mx )x0 =
[6.33]
x=0
Ф = -S dT ( )
=—S Л( To — T) — me ~°,X x=0 [6.34]
x=0
Ф= S Лю( T0 - T ) = 00 S Лю [6.35]
m=. Ihp
[6.36]
ЛЛ8
= S лЛ ( T0 - Te )
Ф= SЛю(T0 -Te) [6.37]
We can compare this flow to Ф 0, the flow that we would have with an
ideal wing; in other words, the isotherm at the temperature of the body needs
to be cooled. This flow Ф0 is calculated by considering:
dв
——- a в = 0 [6.40]
dx
x=0 ; в = Tо - T = 00 [6.41]
d^ d0
x=L ; - 2---- = 0 ^---- = 0 [6.42]
dx dx
We obtain:
B = во [6.46]
We therefore obtain:
BashaL в0 shaL
A =------------- = - 0------------------ = -в0 thaL [6.47]
achaL------ chaL
and then:
- shaLshax + chaLchax
0( x )= 00 chaL
[6.49]
- shaLshax + chaLchax
chaL
We know that:
ch(x-y)=chxchy-shxshy [6.50]
Or even:
dT -mshm(L - x) 1
Ф = -S2----
dx
----ch------
aL
L\ + T. [6.53]
x=0 x=0
a=,\SS [6.55]
Fin Theory: Notions and Examples 249
In the above fin theory, the transfers involved are reduced to conduction
and convection. A wing can also have radiation. Examples of this case have
been attributed to the section dedicated to radiation in Chapter 5.
The convection coefficient between the fin and the ambient air is
h =10W.m-2.K-1.
The surface area of the fin (not including the part welded to the body to
be cooled) is S = 20 cm2.
1) If the entire surface of the wing was at T0, what would the flow Ф0
evacuated by the wing be?
2) The flow evacuated by the wing is, in fact, Ф = 0.9 W. What is the
efficiency n of the fin?
Solution.-
Solution.-
1) We define
The reduced variable:
Oo = To - T [6.62]
Fin Theory: Notions and Examples 251
pe=2(l+e) [6.63]
hp
O)2 = _^e. [6.64]
AS
We obtain:
pe=2(l+e)=2(5.10-2+3.10-3)=0.106m [6.67]
S=le=5.10-2*3.10-3=15.10-5m2 [6.68]
Ф = 7hpe AS 00 th ( aL)
3) The balance equation for the fin must be resolved with the
imposed limit condition. In other words:
d2T
ASdx; + hpe (T - Te ) = 0 [6.75]
0=T(x)-Te [6.76]
00=T0-Te [6.77]
a=_^e hp
[6.78]
AS
x=о ; e = t0 - Te = e0 [6.80]
x=L ; e = Te - Te = 0 [6.81]
x=L;e=0=AshaL+BchaL [6.84]
B = e0 [6.85]
A =---- — [6.87]
thaL
For the expression 0( x) and the flow Ф exiting the wing (calculated at
x = 0), we have:
e1
— =-------- sh ax + ch ax [6.88]
e0 thaL
ф =- S A— - S A—
dx dx
[6.89]
. ( a
=-s a— e0 I---------- shax + ch ax
d x 01 thaL
x=0
254 Heat Transfer 1
_ . 1 r—— 1
Ф = SЛО0 = JS2hpe00- [6.91]
thaL thaL
with:
90.92 - 25.72
=253% [6.94]
25.72
This is explained by the fact that the actual temperature found, presuming
that the flow at the end of the wing is zero, can be written as:
О
= -thaLshax + ch ax [6.95]
О0
This gives:
270
0( L )= T (L)-T =
e
270
ch 0.588
= 229.2 [6.98]
T(L)=229.2+30=259.2C [6.99]
Fin Theory: Notions and Examples 255
We note that this temperature is much higher than T0 . The flow of losses
from the edges is necessarily much lower than would be found if the
temperature descended to Te at the end of the wing.
Knowing that the circuit must not exceed the temperature TC =70 °C,
what is the maximum flow Ф1 that can be evacuated by natural convection?
6) Plate A is, by its very nature, connected to the upper part of the circuit
by a thin layer of x = 0.1 mm of a paste known as “thermal”, based on
silicon with thermal conductibility AP = 0.9 W.m-1.K-1. This paste avoids
the presence of air inclusions between the element and the plate which, even
if it is thin, would be a disturbing thermal resistance. We will denote TS2 as
the new temperature taken at the foot of each wing.
b) Deduce from this the flow evacuated by the radiator in the case where
the ventilator is used. Does this value pose a problem?
Solution.-
Ф1 = 0.9 W [6.102]
The distribution of the temperature along the fin will be given by:
e = e0 th ( aL) [6.103]
The flow at the origin of the fin (flow removed from the body to be
cooled) will be given by:
where:
e0 =TC-Te=45°C [6.105]
Perimeter:
pe =5+5+0.3+0.3=10.6cm=0.106m [6.106]
Surface area:
Parameter:
hp 8 8*0.106p
a = -^- =----------------------- — = 3.83 [6.108]
AAcs у 386*1.5.10-4
and
b) We will easily deduce from this the total flow Ф2 evacuated by the
radiator.
The flow to which the convection of the plane face must be added
becomes:
pe =5+5+0.3+0.3=10.6cm=0.106m [6.115]
260 Heat Transfer 1
Surface area:
Parameter:
hp 30*0.106 „
a= =J---------------- = = 7.41 [6.117]
A^ 3 386*1.5.10-4
and
The flow to which the convection of the plane face must be added
becomes:
4) The only acceptable solution is the one that evacuates more than
19W; that is, the one that includes the ventilator (question 3) is acceptable.
Fin Theory: Notions and Examples 261
19
Tm- = K1F4 + 25 = 975 °C [6'126]
Passing through the paste layer, the flow creates a difference in the
temperature TC- TS2 :
T-T 09
Ф5 = ' x = 10-4 = 9000 (TC — Ts 2) [6.127]
In addition, the relationship between the flow and Te , and TS2 remains
unmodified. We consider again the previous results:
b) The flow evacuated by the radiator when the ventilator is used is then
re-evaluated as:
So:
_ Ф
T - T =—5- [6.132]
C S2 9000
So:
Ф
T - TS 2 =—— [6.134]
es2 0.693
Ф5
(T - T„,)-(T - T„ ) = — [6.135]
Cs2es20.693 9000
1 1
TC - Te =Ф 5 =1.443Ф5 [6.136]
0.693 9000
‘c__ j 45
Ф5 = ----- = 31.2 W [6.137]
1.443 1.443
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
Appendix 1
Heat Equation of a
Three-dimensional System
Figure A1.1. Case of two infinitely close isotherms; relationship with the gradient.
For a color version of this figure, see www.iste.co.uk/ledoux/heat1.zip
The thermal flow that passes through dS1 and dS2 will be expressed by:
dT
dФ = -A d S1 [A1.1]
dx
where x is a measured axis on the axis bearing >7. We should recall that the
minus sign indicates that the heat propagates in the direction of decreasing
temperatures. In this expression, the sign implies a positive flow when it
flows in the direction of the normal.
The same flow of thermal heat obviously passes through dS', dS1 and
dS2 .
dT
grad7 T ) x
A----
dx
[A1.4]
Appendix 1 267
We note that:
In the end, the flow will take the form of a scalar product:
'| ( )
\rrad т||||и cos grad T, n' d S' [A1.7]
Taking into account the orientations of the normal that are specified
above:
Therefore, the flow of heat that passes through a closed surface will be
written as:
ф= jj A grad T. n dS [A1.9]
S
All outgoing flow will thus be counted positively in the integral, and all
flows entering will be counted negatively. The integral therefore represents a
resulting flow balance, positive if the volume contained in the surface loses
heat.
268 Heat Transfer 1
ф= jj2gradT.n dS . [A1.10]
S
This integral is positive if more heat exits than enters, negative if more
enters than exits.
In the former case, the flow will contribute, overall, to “cooling” the region
D. The sensitive quantity of heat integrated into D will diminish over a time
dt. In the latter case, the flow will contribute, overall, to heating the region
D. The sensitive quantity of heat integrated into D will increase over a
time dt.
Ф= jj2gradT.n dS [A1.11]
S
where p is the (local) density of the material and c is its (local) specific heat
capacity. The balance equation can therefore be written as:
dQ
= -Ф [A1.13]
dt
Appendix 1 269
or
Integrating the derivation with respect to the time in the integral and
transforming the integral of the surface into an integral for the volume, we
find
Rewritten as:
jjj
D Pc df-~ divAgradTda = 0 [A1.16]
pc dL = div(AgradT) [A1.17]
(
div Agrad T = 0 ) [A1.18]
The previous expressions are simplified and give forms in which the heat
equation is most generally known:
ДT = 0 [A 1.22]
d T = a at
-- ДT [A1.25]
dt
2
a =--- [A1.26]
Pc
is very practical. Writing equations therefore remains valid for all coordinate
systems. It is therefore necessary to use the known form of the Laplace
operator for each system, in order to have the heat equations in this coordinate
system.
The heat equation results from a balance between the heat that, by
conduction, enters, exits and accumulates in an elementary volume.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
274 Heat Transfer 1
Orthonormalized
Cylindrical system Spherical system
system
Reference and
Reference and coordinates Reference and coordinates
coordinates
Orthonormalized
Cylindrical system Spherical system
system
Elementary volume Elementary volume Elementary volume
In the most general case, the volumetric mass p, the specific heat
capacity °C and, in particular, the thermal conductibility A can vary with
temperature. The heat equation is then written as follows:
Appendix 2 275
^3T д К dT I д К dT I д К dT
pC=I AI + I AI+I A [A2.1]
dt dx 1 dxJ Эy 1 dyJ dz 1 dz
1 d П 2 - dT T | 1 d | д dT | d ( Q . d T |ГЛП/|-.
—I A r2sinФ— I + —---------- 1 A— I + 1 AsinФ I [A2.4]
r 2 sin Ф Э r 1 Э x J sin ФЭ01 d0 J ЭФ1 ЭФ J
a=A [A2.5]
pC
^T = a Д T [A2.6]
dt
дT a д 1 д2 T
д T +, ------------- ^дT
r 2sin Ф----- 1 ---- sin Ф----- [A2.12]
дt r 2sin Ф д r д x sin Ф д 02 дФ дФ
Appendix 2 277
дТ=о
[A2.13]
t
д ( дT i 1 д ( д дT i д Г . л дT
Ar 2sin Ф— I +------------ 1 A— I +------ 1 Asin Ф-----
=0 [A2.16]
дrl д x I sin Фд#1 д0 I дФ l дФ
In this case, the Laplace operator of the temperature is zero. The value of
the physical properties no longer intervenes in the general expression for the
temperature. Note that the thermal conductivity intervenes in the calculations
for flows
AT = 0 [A2.17]
278 Heat Transfer 1
A2.5. Writing the Fourier law. How can the expression for the
gradient be found?
It is therefore useful to know how to find the expression for the gradient
in various coordinate systems. To do so, we need to know:
a) How to write the expression for the components of the small
displacement vector dl due to an elementary modification of the
coordinates.
In the displacement dl, T therefore does not vary. The exact total
differential dT is zero.
дT
дx
-
дT
grad T [A2.26]
дУ
дT
дz
дT
дr
-
grad T 1 IT [A2.29]
r ди
дT
дz
дT 1 дT 1 дT ^ -
=---- dr +-------- rdФ+------------- rsinФdd = gradT. dl = 0 [A2.31]
дr r дФ rsinФ дz
дT
- i"t a
gradT - д— [A2.32]
1 дT
r sin Ф д 0
Appendix 3
In this case, we will no longer talk about an Ox axis, but we direct each
radius positively from the axis of the cylinder towards the outside of the
cylinder. Reasoning in the same way as before, we obtain the following
conclusions.
dT
- A--- 2nrL, with L = 1, that we will no longer write out in the following.
dr
1 Indeed, in the rest of the calculation, we see that the surface through which the conduction
takes place is proportional to r and is therefore going to remain within the derivative.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
284 Heat Transfer 1
dT < 0 when the heat goes from the axis towards the outside, and
dr
dT
- A---- is positive when the heat goes in the direction of increasing r.
dr
dT
— > 0 when the heat goes from the outside towards the axis, and
dx
dT
- A---- is negative when the heat goes in the direction of decreasing r.
dr
dT
Calculate the change in - A--- 2nr when we go from r to r + dr, with
dr
dr being infinitely small:
dT dT dd \^ dT L
- A- 2nr becomes A 2nr +----- 1 A2nr— I dr . [A3.1]
dr dr dr 1 dr I
-
dT
A A—2nr +-- 1 I -A2nrd------
2nr becomes -dT^d dTIU dr . [A3.2]
dr dr dr I dr I
We note that an “r” remains below the second derivative of the last
expression.
Hence the need to work on a linear meter of cylinder and not a m2.
dT d ( dT к
- A-2nr +----- 1 -A2nr---- I dr [A3.3a]
dr dr 1 dr I
Appendix 3 285
dT d
-A 2nr +— dr is positive when the heat goes in the
dr dr
direction of increasing r (s moving away from the axis).
dT d ( dT)
If -A----2nr +---- I -A2nr— I dr is negative, then the heat goes in the
dr dr 1 dr I
direction of the decreasing r (moving closer to the axis).
dTI dT d ( dT V
dQ = pc 2nr drdT = —A 2nr + A 2nr +----- I +A2nr I dr dt [A3.3b]
dr dr dr 1dr J
,
p. crdrdT d I +Ar dT I,
= — „
I dr dt [A3.4]
dr 1dr J
д T 1 d Г dT
---- =------ I Ar---- [A3.5]
дt r дr 1 dr
дT a д ( дT ^
-r- = -V-1 I [A3.6]
дt r дr 1 дr I
Stationary conduction:
d ( dT | ..
— I r---- I = 0 [A3.7]
dr1 dr I
In this case, we no longer talk about an axis Ox, but we positively direct
each radius of the axis of the cylinder to the outside of the sphere.
Reasoning in the same way as before, we reach the following conclusions.
The balance will be established for a space between two spheres that are
infinitely close to radius r and r+ dr.
V = 4nr3 = nd3
[A.9]
^ 6
dT
The flow through a sphere with radius r is then - 2----4nr2.
dr
dT dT
If —< 0, the heat goes from the origin towards the outside, and — 2----
dr dr
is positive when the heat goes in the direction of increasing r.
If dT > 0, the heat goes from the outside towards the axis, and — 2ddT- is
We note that an “r2” remains under the second derivative of the last
expression.
Let us now consider a sphere of radius r+ dr. The quantity of heat that
passes through it per meter of length of the cylinder is:
dT
- A----4nr2 is positive if the heat goes in the direction of increasing r.
dr
In this case, the heat enters the space and the term is counted positively in
the balance.
288 Heat Transfer 1
Through the plane r + dr, the incoming heat must also be counted
positively in the balance.
ddT d I _ dT )
Now -X—2nr +---- 1 -X2nr— I dr is positive when the heat goes in
dr dr 1dr J
the direction of increasing r, in other words when the heat exits the annular
space.
ddT . 2 ddT . 2 d I 2 dT | , d I 2 dT |
-X 4nr2 + X--- 4nr2 +— I X4nr2— I dr = — I X4nr2 r I [A3.12]
dr dr dx 1dr Jdr 1dr J
dT 2dT 2 d 2 dT \
dQ = pc 4nr 2 drdT = -X 4nr + X 4nr X4nr ---- I dr dt [A3.13]
dr dr dx dr J
d 2 dT L
pc 4nr 2 drdT = X4nr ---- I dr dt [A3.14]
dx dr J
д dT 1 d Г 2dT
pc=1 Xr2 [A3.15]
d t r2 d r 1 dr
dT a d ( 9dT^
-:;-=—I r I [A3.16]
d t r2 d r 1 dr J
Appendix 3 289
Stationary conduction:
d
[A3.17]
dr dr J
Hypotheses:
xx
% = ^= = ~n= [A4.1]
at 4a at
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
292 Heat Transfer 1
g = er (g)
f [A4.6]
g-T(x,t)-Ti [A4.7]
( Im } ( x A
0 =T i exP -J— x cos mt -J— x
I 2aa J I 2aa J
[A4.11]
amplitude:
( ,.
m
exp - — x [A4.1a]
I 2aa J
phase shift:
mx x
[A4.13]
a2 a
q
T(x, t)- Ti 2
[A4.i5]
2 - Ti1
Ti2-
Ei 4*i Pi C [A4.i6]
E2 /
— у ^2 P2 C2 [A4.i7]
E
Q = E1+E2 erfc(^)
Appendix 5
These three functions erf (x), erfc ( x) and ierfc(x), are fundamental in
processing single-dimensional non-stationary conduction problems.
x erf (x) erfc (x) ierfc (x) x erf (x) erfc (x) ierfc (x)
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
296 Heat Transfer 1
Complementary Information
Regarding Fins
T =T(x)
Te : ambient temperature
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
298 Heat Transfer 1
We define a parameter a:
hp>
a=,-?■ [A6.1]
AS
Relative temperatures:
0 = T ( x)-T [A6.2a]
00 = To - T [A6.2b]
Flow at the origin of the fin (flow taken from the body to be cooled):
0 = 0o th ( aL) [A6.5]
Flow at the origin of the fin (flow taken from the body to be cooled):
Flow at the origin of the fin (flow taken from the body to be cooled):
In problems that depend on time, the Laplace transform is often used with
respect to this variable.
In practice, we compare the result found with the Laplace transform with
a table of transforms that were previously established. We find many in the
literature.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
302 Heat Transfer 1
A7.1. Definition
The function f (t) must be real, locally integrable, defined for t > 0 . It
must be limited:
|f (t)|<Kea [A7.1a]
It is defined by:
The notation
^ F ( p ) = £ f (t), [A7.2]
df(t)
Transform of a derivative with respect to time-------- :
dt
df(t)
d
£ ( = -f (0) + pF (p) [A7.5]
We see the presence of the term - f (0). In many cases, f (0)= 0, but
we must always keep the existence of this term in mind in each new
problem, in the case where this condition of nullity at the origin of the signal
time is not fulfilled. We also see why it is important that the function is
limited, beyond the fact that this property guarantees the convergence of the
definition integral.
- e
-pt
TO
In the same way, all multiples of the Heaviside function, such that
t<0;f(t)=0 [A7.11]
304 Heat Transfer 1
-ae-pt a
f(t) ^F(p) = 10 ae-ptdt = [A7.13]
p p
exp-axdp
is transformed into F (p ) =------------ — [A7.14]
p
We see that F(p, x), the transform for f (t, x), responds to a differential
equation with one variable.
This time we are looking for the solution via a new reduced temperature:
T-T
в( x, t ) = T - T [A7.16]
ei
d f (x, t) Ъ 2 f (x, t)
[A7.17]
dt a d x2
в( x,0) = 0 [A7.18]
в( 0, t ) = 1 [A7.19]
0(™, t ) = 0 [A7.20]
в( x, t) is transformed into 0( x, p)
£в( x, t) = 0( x) [A7.21]
0(<™, t ) = 0 [A7.22]
£ d^ = p Q-в (t = 0) = p 0 [A7.24]
dt
306 Heat Transfer 1
0( 0, t ) = 1 gives:
£0 (0, t ) = 0( 0, t ) = - [A7.26]
p
0 (0, p) = —, therefore:
and:
B =- [A7.28]
p
fp
exp - x
0( x, p ) =
a [A7.29]
p
1
with a = —;= [A7.31]
fa
0( x, t ) = erfc x [A7.32]
4 at
Reminders Regarding
Hyperbolic Functions
exx - e--xx
hyperbolic sinus shx =---- ^---- [A8.1]
x -x
e +e [A8.2]
hyperbolic cosinus ch x =---- 2----
shx ex - e-x
hyperbolic tangent thx =----- =---------- [A8.3]
ch x ex + e-x
1 ch x ex + e-x
coth x =---- =------=---------- [A8.4]
th x sh x ex - e-x
It must be noted that the similarity stops at these relationships and that
the form of the hyperbolic functions is totally different from that of the
sinusoidal functions.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
310 Heat Transfer 1
sh(-x)=-shx [A8.6]
th (-x)=-th x [A8.8]
dshx
------- = chs [A8.9]
dx
dchx
-------- = shx [A8.10]
dx
dthx 1
[A8.11]
dx ch 2 x
Appendix 8 311
sh(x+y)=shxchy+chxshy [A8.12]
which implies:
sh(x-y)=shxchy-chxshy [A8.14]
[BYR 66] Byron Bird R., Stewart W.E., Lightfoot, E.N., Transport Phenomena,
John Wiley & Sons, New York, 1966.
[COS 97] Cosar P., Aide memoire du thermicien, Amicale des anciens eleves de 1’Ecole de
thermique (ed.), Elsevier, Paris, 1997.
[GOS 93] GOSSE J., Guide technique de thermique, Dunod, Paris, 1993.
[HOL 86] HOLMAN J.P., Heat Transfer, McGraw Hill, New York, 1986.
[JAN 16] JANNOT Y., MOYNE C., Transferts thermiques, Edilivres, Paris, 2016.
[KNU 58] KNUDSEN J.G., KATZ D.L., Fluid Dynamics and Heat Transfer, McGraw Hill,
New York, 1958.
[LEO 79] LEONTIEV A., Theorie des echanges de chaleur et de masse, Editions MIR,
Moscow, 1979.
[PER 97] PEREZ J.P., Thermodynamique, fondements et applications, Masson, Paris, 1997.
[ROH 98] Rohsenow, W.M., Hartnett, J.P., Cho, Y.I., Handbook of Heat Transfer,
McGraw Hill, New York, 1998.
[SAC 15] SACCADURA J.F., Transferts thermiques, initiation et approfondissement, Lavoisier
Editeur, Paris, 2015.
[VAN 76] Van Wylen G., Sonntag R.E., Fundamentals of Classical Thermodynamics,
SI Vers 2, John Wiley & Sons, New York, 1976.
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
Index
C, E F, G, H
conduction fin, 237-240, 242, 243, 245, 246,
non-sationary, 2, 183, 184, 190, 249-252, 255, 256, 258-260,
191, 195, 200, 225, 229, 295, 297-299
301, 304 flow, 1, 3-5, 10, 12-16, 18, 19, 21,
stationary, 2, 25, 26, 285, 289 27, 28, 30, 32-34, 37, 40-44, 46,
convection, 12, 27-29, 37, 44, 46, 51, 48, 53-60, 63-66, 70-72, 75, 77,
65-67, 69-71, 73, 74, 84, 85, 79-90, 92, 95-101, 107, 108, 111,
87-92, 95-101, 103, 105, 106, 110, 114, 117-119, 124, 125, 129, 130,
113, 117, 119, 122, 127-129, 132, 132-135, 137-139, 143-145, 148,
134, 138, 144, 152-155, 159, 162, 149, 152, 153, 156, 157, 159, 160,
165, 167, 168, 172, 173, 218, 223, 163, 166, 168, 172, 174, 175, 177,
237, 239-241, 244, 249, 252, 178, 180, 190, 191, 202, 207, 214,
255-261, 297, 299 224, 225, 233, 235, 237-241, 243,
coefficient, 5, 10, 11, 29, 37, 46, 244, 246, 248-250, 252-261,
66, 67, 69-71, 73, 74, 84, 85, 87, 266-268, 277, 283, 286, 292, 298,
89-92, 95-101, 103, 105, 106, 299
110, 113, 117, 119, 122, 127, density, 4, 5, 10, 19, 28, 30, 32, 33,
128, 132, 134, 138, 144, 152, 46, 53, 54, 56-58, 63, 65, 75, 80,
155, 159, 160, 165, 167, 168, 83, 117, 119, 190, 191, 214, 224,
172, 173, 175, 218, 237, 240, 225, 233, 235, 292
244, 249, 250, 255-257, 271, Fourier’s law, 9, 10, 13, 14, 16, 265,
297, 299 278, 301, 305
error function, 189, 203, 205, 207, gradient, 2, 3, 14, 15, 27, 28, 103,
209, 292 238, 240, 250, 265, 266, 278-281
Heat Transfer 1: Conduction, First Edition. Michel Ledoux and Abdelkhalak El Hami.
© ISTE Ltd 2021. Published by ISTE Ltd and John Wiley & Sons, Inc.
316 Heat Transfer 1
heat thermal
equation, 1, 2, 14, 16-18, 23-25, conductibility, 11, 12, 22-24, 29,
27-30, 35, 62, 104, 183, 185, 30, 33, 34, 37, 38, 41, 47, 54-56,
186, 191, 194, 201, 265, 58, 62, 64, 70, 77, 82, 85, 86, 90,
268-271, 273-277, 283, 304 91, 93, 99, 100, 112, 129, 142,
propagation, 14, 18, 266 143, 152, 153, 156, 161, 167,
176, 184, 205, 210, 213, 228,
I, R, T 232, 256, 257, 268, 271, 274,
275, 291, 297
isothermal, 2, 3, 9, 14, 18, 24, 44, 47,
conduction, 1, 2, 10, 12, 16, 27,
191, 194, 200, 238, 240, 246, 249,
239, 268
252, 265, 279, 280, 283, 286
diffusivity, 23, 184, 190, 192, 194,
surface, 2, 3, 14, 265
199, 201, 202, 206, 208, 211,
radiation, 12, 27-29, 105, 107, 249
220, 225, 230, 233, 270, 271,
temperature, 1-3, 5, 6, 9-17, 22, 24,
275, 291
26-31, 34-38, 40, 44, 47, 49, 51,
flow, 4, 10, 14, 16, 28, 46, 48,
54, 55, 57-62, 65, 69, 71, 73, 77,
54-58, 65, 70, 71, 77, 90, 92,
81, 84-89, 95, 99, 100-104,
124, 132-134, 137, 153, 154,
107-111, 114, 116-118, 122-125,
166, 191, 233, 246, 266
127, 128, 130-134, 137, 138, 143,
ladder, 190
144, 147, 148, 152, 153, 156-160,
resistance, 29, 30, 33, 34, 37, 40,
162, 164-166, 171, 173-176,
43, 53-59, 61-67, 69-79, 81-88,
178-180, 183-186, 190-195,
90-93, 96-98, 100, 103, 106,
199-211, 213-233, 236-241, 246,
110-113, 116, 122-124,
249, 250, 254-258, 261, 262, 265,
127-130, 132, 134, 144, 145,
266, 270, 274, 275, 277-279, 283,
148, 153, 157, 160-162, 165,
286, 291-293, 297-299, 304, 305
167, 168, 171, 172, 174, 177,
217, 257
Other titles from
in
Mechanical Engineering and Solid Mechanics
2021
DAHOO Pierre Richard, POUGNET Philippe, EL HAMI Abdelkhalak
Applications and Metrology at Nanometer Scale 1: Smart Materials,
Electromagnetic Waves and Uncertainties
(Reliability ofMultiphysical Systems Set - Volume 9)
2020
SALENQON Jean
Elastoplastic Modeling
2019
BAYLE Franck
Reliability ofMaintained Systems Subjected to Wear Failure Mechanisms:
Theory and Applications
(Reliability ofMultiphysical Systems Set - Volume 8)
BEN KAHLA Rabeb, BARKAOUI Abdelwahed, MERZOUKI Tarek
Finite Element Method and Medical Imaging Techniques in Bone
Biomechanics
(Mathematical and Mechanical Engineering Set - Volume 8)
IONESCU Ioan R., QUEYREAU Sylvain, PICU Catalin R., SALMAN Oguz Umut
Mechanics and Physics of Solids at Micro- and Nano-Scales
LE VAN Anh, BOUZIDI Rabah
Lagrangian Mechanics: An Advanced Analytical Approach
MICHELITSCH Thomas, PEREZ RIASCOS Alejandro, COLLET Bernard,
NOWAKOWSKI Andrzej, NICOLLEAU Franck
Fractional Dynamics on Networks and Lattices
SALENQON Jean
Viscoelastic Modeling for Structural Analysis
2018
BOREL Michel, VENIZELOS Georges
Movement Equations 4: Equilibriums and Small Movements
(Non-deformable Solid Mechanics Set - Volume 4)
FROSSARD Etienne
Granular Geomaterials Dissipative Mechanics: Theory and Applications in
Civil Engineering
RADI Bouchaib, EL HAMI Abdelkhalak
Advanced Numerical Methods with Matlab® 1: Function Approximation
and System Resolution
(Mathematical and Mechanical Engineering SET - Volume 6)
Advanced Numerical Methods with Matlab® 2: Resolution ofNonlinear,
Differential and Partial Differential Equations
(Mathematical and Mechanical Engineering SET - Volume 7)
SALENQON Jean
Virtual Work Approach to Mechanical Modeling
2017
BOREL Michel, VENIZELOS Georges
Movement Equations 2: Mathematical and Methodological Supplements
(Non-deformable Solid Mechanics Set - Volume 2)
Movement Equations 3: Dynamics and Fundamental Principle
(Non-deformable Solid Mechanics Set - Volume 3)
BOUVET Christophe
Mechanics ofAeronautical Solids, Materials and Structures
Mechanics ofAeronautical Composite Materials
BRANCHERIE Delphine, FEISSEL Pierre, BOUVIER Salima,
IBRAHIMBEGOVIC Adnan
From Microstructure Investigations to Multiscale Modeling:
Bridging the Gap
Chebel-Morello Brigitte, NICOD Jean-Marc, VARNIER Christophe
From Prognostics and Health Systems Management to Predictive
Maintenance 2: Knowledge, Traceability and Decision
(Reliability ofMultiphysical Systems Set - Volume 7)
EL HAMI Abdelkhalak, RADI Bouchaib
Dynamics ofLarge Structures and Inverse Problems
(Mathematical and Mechanical Engineering Set - Volume 5)
Fluid-Structure Interactions and Uncertainties: Ansys and Fluent Tools
(Reliability ofMultiphysical Systems Set - Volume 6)
KHARMANDA Ghias, EL HAMI Abdelkhalak
Biomechanics: Optimization, Uncertainties and Reliability
(Reliability ofMultiphysical Systems Set - Volume 5)
LEDOUX Michel, EL HAMI Abdelkhalak
Compressible Flow Propulsion and Digital Approaches in Fluid Mechanics
(Mathematical and Mechanical Engineering Set - Volume 4)
Fluid Mechanics: Analytical Methods
(Mathematical and Mechanical Engineering Set - Volume 3)
MORI Yvon
Mechanical Vibrations: Applications to Equipment
2016
BOREL Michel, VENIZELOS Georges
Movement Equations 1: Location, Kinematics and Kinetics
(Non-deformable Solid Mechanics Set - Volume 1)
BOYARD Nicolas
Heat Transfer in Polymer Composite Materials
CARDON Alain, ITMI Mhamed
New Autonomous Systems
(Reliability ofMultiphysical Systems Set - Volume 1)
DAHOO Pierre Richard, POUGNET Philippe, EL HAMI Abdelkhalak
Nanometer-scale Defect Detection Using Polarized Light
(Reliability ofMultiphysical Systems Set - Volume 2)
de Saxce Gery, Vallee Claude
Galilean Mechanics and Thermodynamics of Continua
DORMIEUX Luc, KONDO Djimedo
Micromechanics ofFracture and Damage
(Micromechanics Set - Volume 1)
EL HAMI Abdelkhalak, RADI Bouchaib
Stochastic Dynamics ofStructures
(Mathematical and Mechanical Engineering Set - Volume 2)
GOURIVEAU Rafael, MEDJAHER Kamal, ZERHOUNI Noureddine
From Prognostics and Health Systems Management to Predictive
Maintenance 1: Monitoring and Prognostics
(Reliability ofMultiphysical Systems Set - Volume 4)
KHARMANDA Ghias, EL HAMI Abdelkhalak
Reliability in Biomechanics
(Reliability ofMultiphysical Systems Set-Volume 3)
MOLIMARD Jerome
Experimental Mechanics ofSolids and Structures
RADI Bouchaib, EL HAMI Abdelkhalak
Material Forming Processes: Simulation, Drawing, Hydroforming and
Additive Manufacturing
(Mathematical and Mechanical Engineering Set - Volume 1)
2015
KARLICIC Danilo, MURMU Tony, ADHIKARI Sondipon, MCCARTHY Michael
Non-local Structural Mechanics
SAB Karam, LEBEE Arthur
Homogenization ofHeterogeneous Thin and Thick Plates
2014
ATANACKOVIC M. Teodor, PILIPOVIC Stevan, STANKOVIC Bogoljub,
ZORICA Dusan
Fractional Calculus with Applications in Mechanics: Vibrations and
Diffusion Processes
Fractional Calculus with Applications in Mechanics: Wave Propagation,
Impact and Variational Principles
CIBLAC Thierry, MOREL Jean-Claude
Sustainable Masonry: Stability and Behavior of Structures
ILANKO Sinniah, MONTERRUBIO Luis E., MOCHIDA Yusuke
The Rayleigh-Ritz Methodfor Structural Analysis
LALANNE Christian
Mechanical Vibration and Shock Analysis - 5-volume series - 3rd edition
Sinusoidal Vibration - Volume 1
Mechanical Shock - Volume 2
Random Vibration - Volume 3
Fatigue Damage - Volume 4
Specification Development - Volume 5
LEMAIRE Maurice
Uncertainty and Mechanics
2013
ADHIKARI Sondipon
Structural Dynamic Analysis with Generalized Damping Models: Analysis
ADHIKARI Sondipon
Structural Dynamic Analysis with Generalized Damping Models:
Identification
BAILLY Patrice
Materials and Structures under Shock and Impact
BASTIEN Jerome, BERNARDIN Frederic, LAMARQUE Claude-Henri
Non-smooth Deterministic or Stochastic Discrete Dynamical Systems:
Applications to Models with Friction or Impact
EL HAMI Abdelkhalak, RADI Bouchaib
Uncertainty and Optimization in Structural Mechanics
KIRILLOV Oleg N., PELINOVSKY Dmitry E.
Nonlinear Physical Systems: Spectral Analysis, Stability and Bifurcations
LUONGO Angelo, ZULLI Daniele
Mathematical Models ofBeams and Cables
SALENQON Jean
Yield Design
2012
DAVIM J. Paulo
Mechanical Engineering Education
DUPEUX Michel, BRACCINI Muriel
Mechanics of Solid Interfaces
ELISHAKOFF Isaac et al.
Carbon Nanotubes and Nanosensors: Vibration, Buckling
and Ballistic Impact
GREDIAC Michel, HILD Francois
Full-Field Measurements and Identification in Solid Mechanics
GROUS Ammar
Fracture Mechanics - 3-volume series
Analysis ofReliability and Quality Control - Volume 1
Applied Reliability - Volume 2
Applied Quality Control - Volume 3
RECHO Naman
Fracture Mechanics and Crack Growth
2011
KRYSINSKI Tomasz, MALBURET Francois
Mechanical Instability
SOUSTELLE Michel
An Introduction to Chemical Kinetics
2010
BREITKOPF Piotr, FILOMENO COELHO Rajan
Multidisciplinary Design Optimization in Computational Mechanics
DAVIM J. Paulo
Biotribolgy
PAULTRE Patrick
Dynamics ofStructures
SOUSTELLE Michel
Handbook ofHeterogenous Kinetics
2009
BERLIOZ Alain, TROMPETTE Philippe
Solid Mechanics using the Finite Element Method
LEMAIRE Maurice
Structural Reliability
2007
GIRARD Alain, ROY Nicolas
Structural Dynamics in Industry
Guinebretiere Rene
X-ray Diffraction by Polycrystalline Materials
KRYSINSKI Tomasz, MALBURET Francois
Mechanical Vibrations
KUNDU Tribikram
Advanced Ultrasonic Methods for Material and Structure Inspection
SIH George C. et al.
Particle and Continuum Aspects ofMesomechanics