GSIL1A, B, C A.

U 2023/2024 ENI Carthage

Série N°2Transfert Thermique

Ex 1: Conduction thermique dans un milieu uniforme (milieu mince)

Pour réchauffer du lait pour bébé, la mère verse le lait dans une mince paroi de verre dont le
diamètre est de 6 cm. La hauteur du lait dans le verre est de 7 cm. Elle place ensuite le verre
dans une grande casserole remplie d'eau chaude à 60°C. Le lait est agité en permanence, de
sorte que sa température est uniforme en tout temps. Si le coefficient de transfert de chaleur
entre l'eau et le verre est 120W/m2.°C, déterminer le temps qu’il faut pour que le lait se
réchauffe de 3°C à 38°C. En admet que les propriétés du lait sont les mêmes que ceux de
l'eau. Dans ces conditions, le lait peut être traité comme un milieu thermiquement mince,
Pourquoi ?
On donne : La conductivité thermique, densité et la chaleur spécifique de l’eau à 20°C sont
λ = 0.607 Wm-1.°C-1, ρ = 998 kg.m-3 et Cp = 4.182 kJ.kg-1.°C-1.
Solution :
Commençant par calculer le nombre de Biot
V π r2 L π 0.03 2 0.07 −2
LC = = 2
= 2
=1.05 10 m
S 2 π r 0 L+2 π r 0 2 π 0.03 0 0.07+ 2 π 0.030
h Lc (120)( 0.0105)
=2.076> 0.1 Le nombre de Biot dans ce cas est très important
Bi=
λ = 0.607
(beaucoup plus supérieur à 0.1). Toute fois le système globale d’un système uniforme est
encore applicable puisque le lait est constamment agité de telle façon à garder sa
température uniforme dans toute la bouteille.
T (t )−T ∞ hS h 38−60 120
=exp−( BiFo)=exp−( )t=exp−( )t ⇒ =exp−( )t
T i−T ∞ ρ Cp V ρ C p Lc 3−60 998 4182 0.0105
t = 348 s = 5min48s le temps nécessaire pour réchauffer le lait de 3°C à 38°C

Ex 2 : Conduction thermique dans un milieu épais (Méthode de séparation des

variable).
Soit un fil isolé de longueur L = 1m, de diffusivité thermique α = 3 10-3 m².s-1. Les deux
extrémités de ce fil sont plongées dans la glace. La température est initialement sous la
forme T(x,0) = 50x(1-x).
1- Déterminer l’expression de la température T(x,t) en utilisant la méthode de séparation de
variable pour un milieu thermique épais.

1
GSIL1A, B, C A.U 2023/2024 ENI Carthage

Utilisant la méthode de séparation des variables T(x,t) = f(x) g(t).

d 2 T 1 dT d 2 f f (x ) g (t)
L’équation de propagation est = α dt ⇒ g (t) = α
d x2 d x2 dt

f ' ' g' −α k t

2

=
2
=−k qui donne g(t )=C1 e et f ( x)=C 2 eikx +C3 e−ikx
f αg
Les conditions aux limites sont:T(0,t) = T(L,t) = 0 °C puisque les deux extimités du fil sont
plongées dans la glace.
ikx −ikx
T (0 , t )=0 ⇒ f (x=0)=C 2 +C3 =0 ⇒C 2=−C3 ⇒ f (x)=C(e −e )=C sin(kx )
nπ
T ( L ,t)=0 ⇒C sin(kL)=0 ⇒kL=n π ⇒ k=
L
2
nπ
−α ( )t nπ
Par la suite g(t )=C1 e L
et f (x)=C sin( x)
L
2
∞ nπ
−α ( )t
nπ
T (x , t)=∑ A n e L
sin( x)
n=1 L
∞
nπ 1π
T ( x , 0)=∑ A n sin ( x )=f (x)≈f max sin ( x)
n=1 L L
∞
nπ
On peut utiliser la décomposition en série de Fourier avec f ( x)=∑ A n sin( x )=50 x (1−x)
n=1 L
1 n 0 sin est paire
nπ x 200 (−1) 200
qui implique que An =2∫ 50 x (1−x) sin( )dx= 3 3 − { 400 }
0 L π n π 3 n3 3 3
si n est impaire
π n
Pour les valeurs de n supérieur à 1 les An sont négligeables
2- Déduire le temps nécessaire pour que la température de fil rechute de la moitié de la
température maximal.
A t = 0 T (x , t=0)=f ( x)=50 x (1−x )⇒ f ' ( x )=50 x −100⇒ f ' (x)=0 pour x=0.5 m
Ce qui implique que Tmax = (0.5, 0) = 12.5°C
la température de fil rechute de la moitié de la température maximal
T max
T (0.5 , t 1 )= =6.25° C
2
2
∞ nπ ∞
−α ( ) t1
nπ nπ
T (0.5 , t 1 )=∑ An e )=∑ A n e
L −α (n π ) t 2

sin ( sin( )=6.25 ° C

1

n=1 2 L n=1 2
2
−α (n π ) t1 nπ
T (0.5 , t 1 )=T max e sin( )=6.25 ° C
2

2
GSIL1A, B, C A.U 2023/2024 ENI Carthage

2
1π
T (0.5 , t 1 )=12.5 e−α (1 π ) t sin( 1
)=6.25° C
2
2
− α (1 π ) t 1 6.25 −log (0.5)
e = ⇒t 1= =23.41 s
12.5 απ2
Ex 3 :Milieu semi-infini : Densité de flux imposée
I-flux constant : On considère un mur, initialement à une température uniforme de Ti = 5
°C. Ce mur est exposé de 9:00 à midi à un éclairement solaire égal à 650 W.m 2 dont 40 %
est absorbé.
Φ0 Milieu semi-infini

Ti = T (x, t = 0 )
x

1-En utilisant la transformée déterminer l’expression de la température T(x,t).

T ( x , 0)=T i
2
∂ T 1 ∂T
L’équation de chaleur s’écrit : 2 = α avec les conditions aux limites T (∞ , 0)=T i
∂x ∂t ∂ T (0 , t)
−λ =Φ0
∂x
La dernière condition traduit la conservation du flux de chaleur au niveau de la surface du
milieu semi-infini. Posons θ ( x ,t )=T ( x , t)−T i les conditions aux limites deviennent :
θ (x , 0)=0
θ (∞ ,0)=0
∂ θ (0 , t )
−λ =Φ0
∂x
La transformée de Laplace de l’équation de propagation de chaleur est
∂ 2 θ = 1 ( P θ −θ ( x , 0))= P θ (x , P) puisque θ(x,0) = 0.
α α
∂ x2
θ ( x , 0)=0
La transformée de Laplace des conditions aux limites sera : θ (∞ , 0)=0
∂ θ (0 , P) Φ0
−λ =
∂x P
P
Posons k 2= α la solution de l’équation de chaleur est θ ( x , P)= A e kx + B e−kx. d’après les
conditions aux limites θ (∞ , P)=A e kx + B e−kx=0⇒ A=0 ⇒θ ( x , P)=B e−kx
∂ θ (0 , t) Φ Φ
=−kB= 0 ⇒ B= 0
∂x λP λkP

3
GSIL1A, B, C A.U 2023/2024 ENI Carthage

Φ 0 e− kx
Par la suite on aura θ (x , p)= λ .
kP

La transformée de Laplace inverse nous donne :θ ( x ,t )=

λ √
Φ0 4 α t
[ π exp(
−x2
4α t
)−x erfc (
x
2 √α t
)]

√
−qx 2
e αt −x x
On donne : L ⟨ ⟩=2 π exp( )−x erfc( )
pq 4α t 2 √α t
u
2 x
∫ e−t dt tel que u= 2
2

erfc (u)=1−erf (u)avec erf (u)= √ π

0 √α t
2-Déterminer la température de surface et la température à une profondeur de 4 cm.

√
2
-1 -1 −4 u
-7 2 -1
On donne : λ = 0.552 W.m .K et α = 8.26×10 m s et erf (u)≈ 1−exp ( π )
-2
t1 = 3(3600) = 10800s ; Φ0 = 0.4(650) = 260 W.m ( 40 % du flux solaire est absorbé)
- Calcul de la température à la surface

θ (0 , t 1)=
260
0.552
[ √
4×8.26×10−7 ×10800
π exp(
−0 2
4α t
)−0 erfc (
0
2 √α t
)]=50.21 °C

⇒T (0 , t 1)=θ (0 ,t 1)+T i=50.21+5=55.21 ° C

- Calcul de la température à une profondeur x = 0.04m

√
−7 2
260 4×8.26×10 ×10800 −0.04 0.04
θ (0.04 ,t 1)= [ π exp ( )−0.04 erfc ( )]
0.552 2 √ 8.26×10 t
−7 −7
4 8.26×10 t

√
−7 2
260 4×8.26×10 ×10800 −0.04 0.04
θ (0.04 ,t 1)= [ π exp ( )−0.04(1−erf ( ))]
0.552 2 √ 8.26×10 t
−7 −7
4 8.26×10 t

erf (
0.04
2 √ 8.26×10−7 t
)=erf (0.21) avec erf (0.21)= 1−exp( √ −4×0.212
π )=0.233522

θ (0.04 ,t 1)=33.56° C ⇒T (0.04 ,t 1)=θ (0.04 , t 1)+T i =33.56+5=38.56 ° C

II-Flux convectif φ=h[T∞ - T(0, t)] :

Φ = h ( T∞ - T(x=0,t)) Milieu semi-infini

Ti = T (x, t = 0 )
x

4
GSIL1A, B, C A.U 2023/2024 ENI Carthage

On soumet une plaque épaisse d’aluminium, à une condition de flux convectif a sa surface
dont la température du fluide et le coefficient d’échange par convection sont respectivement
T∞= 25 °C et h = 150 W.m-2.K-1. La température initiale était de 250 °C.
1- déterminer l’expression de la température T(x,t). On donne
−qx
e 1 x 1 x
+h √α t )
2
L⟨ ⟩= erf ( )− exp (hx +α t h )erfc(
p (q+ h) h 2 √α t h 2 √α t
Suivant le même raisonnement que la question précédente les conditions aux limites seront :
T ( x ,0)=T i θ ( x , 0)=0
T (∞ , 0)=T i siθ (x ,t )=T (x , t)−T i θ (∞ ,0)=0
∂ T (0 , t) ∂ θ (0 , t)
−λ =h⟦T ∞ −T ( x=0 , t)⟧ −λ =h⟦(T ∞ −T i )−θ ( x=0 , t)⟧
∂x ∂x
La transformée de Laplace des ces conditions aux limites
θ (x ,0)=0
θ (∞ , 0)=0
∂ θ (0 , P) h (T ∞ −T i)
−λ =−h θ (x=0 , P)+
∂x P
P
Posons k 2= α la solution de l’équation de chaleur est θ ( x , P)= A e kx + B e−kx. d’après les
conditions aux limites θ (∞ , P)=A e kx + B e−kx=0⇒ A=0 ⇒θ ( x , P)=B e−kx .
h(T ∞ −T i ) h
∂ θ (0 , P) h (T ∞ −T i) (T ∞−T i) P (T ∞ −T i)
−λ =−h θ ( x=0 , P)+ ⇒ λ Bk =−hB+ ⇒ B= =λ
∂x P P λ k+ h h
P( +k )
λ
e−kx h
Finalement θ ( x , p)=l(T ∞ −T i ) avec l= .
P(k +l) λ

La transformée de Laplace inverse nous donne

1 x 1 x
+l √α t))
2
θ ( x ,t )=(l(T ∞ −T i ))( erf ( )− exp(lx+ α t l ) erfc (
l 2 √α t l 2 √α t
2
x h h x h
T ( x , t)−T i=(T ∞ −T i)(erf ( )−exp ( x +α t ( ) )erfc( + √α t))
2 √α t λ λ 2 √α t λ

2- Quelle est alors la température de la plaque à une profondeur de 5 cm sous la surface

après 10 minutes.
On donne : λ = 230 W.m-1.K-1 et α = 9.465×10-5 m²s-1.
T (x , t)−T i x h h 2 x h
=(erf ( )−exp ( x +α t ( ) )erfc ( + √ α t))
T ∞−T i 2 √α t λ λ 2 √α t λ

5
GSIL1A, B, C A.U 2023/2024 ENI Carthage

T ( x , t)−T i h h 2
=erf (u)−exp( x + α t ( ) )erfc( v) avec
T ∞−T i λ λ

x
u=
T (x , t)−T i x h h
2
x h 2 √α t
=(erf ( )−exp ( x +α t ( ) )erfc ( + √α t))avec
T ∞−T i 2 √α t λ λ 2 √α t λ x h
v= + √α t
2 √α t λ
0.05
u= =0.105
2 √ 9.465 10−5 600
0.05 150 √ 9.46510−5 600
v= + =0.26
2 √ 9.465 10−5 600 230

√
erf (u)=erf (0.105)= 1−exp(
−4×0.105 2
π )=0.118

√
erf ( v)=erf (0.26)= 1−exp(
−4×0.26 2
π )=0.286
T ( x , t)−T i h h
2
=erf (u)−exp( x + α t ( ) )(1−erf (v ))
T ∞−T i λ λ
T (0.05,600)−25 150 −5 150
2
=0.118−exp( 0.05+ 9.46510 600 ( ) )(1−0.286) ⇒T (0.05,600)=221.58 °C
250−25 230 230