Etablissement des équations du flux de chaleur et de la température

Considérons une couche d'un matériau solide limitée par deux surfaces cylindriques coaxiales de rayons respectifs r0 et r1 et de longueur L supposée très grande par rapport à r1 (ceci permet de supposer que la propagation de la chaleur ne s'effectue que dans le sens radial et qu'il n'y a pas d'effet d'extrémité).

Soient θ0 et θ1 les températures respectives des deux surfaces (on considèrera θ0 > θ1). Supposons que la conductivité thermique est indépendante de la température. Pour des raisons de symétrie, les surfaces isothermes sont des surfaces cylindriques coaxiales.

En coordonnées cylindriques, à symétrie axiale avec L>>R, la variable z n'intervient pas.

En régime permanent, le terme d'accumulation est nul.

Sans génération, le terme de génération n'apparaît pas.

L'équation (30) se simplifie en :

d dr ( r d θ dr ) = 0 d over {dr} (r {{d%theta} over {dr}}) = 0 ..........................................................................(équation 44)

La résolution de cette équation permet d'établir les équations du flux de chaleur et de la température :

Q = 2 π λ L θ 0 θ 1 ln r 1 r 0 Q = 2 %pi %lambda L {{%theta_0} - {%theta_1}} over {ln{r_1 over r_0}} ........................................................(équation 45)

et

( θ θ 0 ) ( θ 1 θ 0 ) = ln r ln r 0 ln r 1 ln r 0 {({%theta} - {%theta_0})over ({%theta_1} - {%theta_0})} = {{ln r - ln r_0} over {ln r_1 - ln r_0}} ...............................................................(équation 46)

ComplémentEtablissement des équations (45) et (46)

En utilisant l'aire moyenne logarithmique A ml = A 1 A 0 ln A 1 A 0 A_{ml} = {A_1 - A_0} over {ln {A_1 over A_0}} , on peut écrire le flux sous la forme :

Q = θ 0 θ 1 e λ A ml Q = {{%theta_0} - {%theta_1}} over {e over {%lambda A_ml}} ......................................................................(équation 47)

Remarque

Les termes des dénominateurs expriment la résistance thermique de la couche cylindrique.

Si les 2 parois échangent par convection avec des milieux fluides (h coefficient de convection, θa températures des fluides, A surface des parois), on a alors :

Q = θ a 0 θ a 1 1 h 0 A 0 + e λ A ml + 1 h 1 A 1 Q = {{%theta_a0} - {%theta_a1}} over {{1 over {h_0 A_0}} + {e over {%lambda A_ml}} + {1 over {h_1 A_1}}} ............................................................(équation 48)

Cette formulation sera utile dans le cours de génie thermique pour calculer le coefficient d'échange global dans un échangeur.

.

Pour une couche cylindrique composite, le flux de chaleur s'écrit :

Q = θ a 0 θ an 1 h 0 A 0 + Σ e λ A mli + 1 h n A n Q = {{%theta_a0} - {%theta_an}} over {{1 over {h_0 A_0}} + {%SIGMA}{e over {%lambda A_mli}} + {1 over {h_n A_n}}} .................................................(équation 49)

où Ai = 2 π ri L et A mli = A i A i 1 ln A i A i 1 A_{mli} = {A_i - A_{i-1}} over {ln {A_i over A_{i-1}}} .

.