TP N° 1 : Les ailettes

I. Répartition des températures le long des ailettes trapézoïdales et rectangulaires.

Considérons deux ailettes de formes distinctes, une trapézoïdale et l’autre rectangulaire. Les deux sont faites du
même matériau de conductivité thermique k=80,2 W/m.K et ayant chacune 20 mm d'épaisseur et 80 mm de
longueur.

L’installation des deux ailettes sur un mur refroidi à l’air a pour but d’augmenter la chaleur évacuée par ce mur.
La température du mur est Tb = 115°C et le coefficient de convection entre la surface du solide et l’air (T=15°C)
est h = 100 W/m2.K.

Figure 1.

1. En recourant au logiciel Matlab, on se propose d'étudier seulement la moitié de l’ailette. Justifier ce

choix :
…………………………………………………………………………………………………………………
…………………………………………………………………………………………………………………
…………………………………………………………………………

2. Ailette de section trapézoïdale

En supposant un transfert de chaleur en deux dimensions. Déterminer, par simulation numérique, les
températures à différentes positions 0  x  L, y  0  . Tracer sur le graphe de la figure 2 la répartition des
températures le long d'ailette.

X 0

T(x, y=0) 115

3. Ailette de section rectangulaire

En maintenant les mêmes hypothèses que dans l’ailette trapézoïdale, déterminer par simulation puis tracer
sur même graphe (figure 2) avec un motif différent la répartition des températures le long de l'ailette rectangulaire.

Comparer le comportement des deux ailettes.

X 0

T(x, y=0) 115

1 / 13
 Figure 2 : Répartition de la température en fonction de x.

II. Confrontation entre les résultats numériques et les résultats analytiques.

Nous proposons ici d'étudier par deux méthodes le comportement thermique d'ailette de section rectangulaire. Dans
cette étude, nous supposons que le flux à la surface extrême x=L est négligeable et égale à zéro.

1. Dans ses conditions, montrer que la solution analytique du problème s'écrit :

chm( L  x)
T ( x)  T  (Tb  T ) (Eq. 1)
chmL

Avec : m=h.P/(.S)

On donne : S=1*e et P=2*(e+1)

2. Déterminer par simulation la répartition des températures le long de l’ailette (rassembler les
résultats dans le tableau ci-dessous).
3. A partir de la solution analytique (Eq. 1), calculer les températures T(x).
X 0

T(°C) 115
analytique

T(°C) 115
numérique

2 / 13
 4. Tracer sur le même graphe ci-dessous (T(x) de la solution analytique tracée auparavant sur le
graphe) la répartition des températures le long d'ailette de la solution numérique. Comparer les
résultats des deux approches et conclure.

5. Déterminer par simulation la densité de flux (0,y) à la base de l’ailette.

Pour afficher les densités des flux, aller dans le menu Plot et sélectionner Paramètres. Dans la boite de
dialogue plot sélectionner apparait ; à la place de la température sélectionner heat flux. Pour plus de précision, il
est préférable de faire un zoom au voisinage de x=0. Pour cela, sélectionnez le menu Options puis l'item Axes
Limits.... Entrez les bornes du domaine du zoom, par exemple [0, 0.001] en x.

y 0 L

(0,y)

6. En utilisant des hypothèses simplificatrices, démontrer que la solution analytique du flux transmis
par l'ailette est égale à :
dT
  .S  P.h. .S (Tb  T ).thmL
dx x 0

Les hypothèses simplificatrices sont :

- La longueur dans le sens transversal est assez grande par rapport à l’épaisseur de l’ailette.

- La température du fluide est uniforme.

- La résistance de contact entre l’ailette et la paroi mère est nulle.

Calculer la densité du flux 0.

7. Tracer la variation de la densité de flux à la base de l'ailette en fonction de Y. Comparer avec le

résultat analytique.

3 / 13
 y(m)

8. Déterminer à partir du graphe la densité de flux moyenne moyenne.

e/ 2 e / 2
1 2
 moyenne    (0, y ) dy    (0, y ) dy
e e / 2 e 0

9. Déterminer la valeur de la constante a tel que :

 moyenne  a. (0,0) . a =……………..

III. Etude d'efficacité thermique E d'ailette

puissance évacuée avec l' ailette
L'efficacité E de l’ailette est donnée par : E =
puissance évacuée sans ailette

e / 2
2
e    0, y  .dy moyenne
E 0

h(Tb  T ) h.(Tb  T )

1. Ailette de section rectangulaire

Le choix d'une forme géométrique convenable d'une ailette exige un compromis entre le prix, le poids et
l'espace disponible. Pour une surface plane d'aire S, la résistance thermique est 1/ h*S. L'emploi des ailettes
augmente l'aire de la surface, mais en même temps, il introduit aussi une résistance thermique de conduction sur
cette partie de la surface initiale à laquelle sont liées les ailettes. Ainsi, la présence d’ailette n'aura toujours pas pour
effet une augmentation du flux de chaleur.

Par la suite, en suppose que la densité de flux moyenne est liée avec  (0, 0) (la densité de flux au point
(0.0)) par l’équation  moyenne  a. (0,0) avec a=1,077. Dans ces conditions, compléter le tableau ci-dessous. Pour
chaque valeur de L, simuler le problème par Matlab et calculer  (0,0) .

4 / 13
 L(m) m.L h(Tb-T) (0,0) moyenne E=moyenne/h(Tb-
T)

0,01 1,27 10000

0,09 11,45 10000

0,18 22,89 10000

0.27 34,34 10000

0.35 44,51 10000

0,44 55,96 10000

0,54 68,68 10000

0,60 76,31 10000

Tracer sur le graphe ci-dessus E en fonction de la variable m.L (grandeur sans dimension)

Interpréter l'allure du graphe. Que représente m.L=3. Conclure.

2. Mise en évidence de l'effet de fluide sur l'efficacité d'une ailette

Dans cette partie, on se propose de déterminer l'effet du fluide sur l'efficacité de l'ailette. Pour se faire,
plusieurs cas de fluides sont à étudier (voir tableau) dans les mêmes conditions que dans le paragraphe (III-1) et
pour une longueur L fixe égale à 0,08 m.

Pour chaque valeur de h du tableau ci-dessous, déterminer par simulation la densité du flux (0.0). Compléter le
tableau ci-dessous.

5 / 13
 Fluide h h(Tb-T) 0,25*P./S  (0,0) moyenne = E=

1,077*  (0,0) moyenne/h(Tb-T)

Air, convection naturelle 10 1000 2045

Eau, convection naturelle 100 10000 2045

Air, convection forcée 200 20000 2045

Huiles, convection forcée 3000 300000 2045

Eau, convection forcée 4000 400000 2045

Eau, ébullition 6000 600000 2045

Eau, condensation 10000 1000000 2045

Conclusion :.……………………………………………………………………………………………….…………
……………………………………………………………………………….………………….……………………
…………………………………………………………………….………………….………………………………
…………………………………...…………………..…………………………………………………………………
…………………………….…………………………….……

6 / 13
 TP2 : Conduction bidimensionnelle

en régime stationnaire avec terme de production.

Une paroi composée de 2 milieux homogènes de conductivité1 et 2, et d’épaisseur e1 et e2respectivement.

La première paroi est le siège d’une distribution de source volumiqueQv = 60 W/cm³.

Les conditions aux limites sont:

- Condition de Dirichlet à gauche T=T0

 hT  T f 
dT
- Condition de Newton (convection) à droite  
dx

Le contact est parfait entre les deux milieux.

On donne : 1  0,65 W / cm C ,  2  0,2 W / cm C , h  0,8 W / cm 2 C , T0  500 C , T f  150 C ,

e1  3 cm , e 2  2 cm , L  1 cm .

Les deux parois se comportent comme des Murs au sens thermique.

1. Donner les conditions aux limites des domaines y=0 et y=1cm.

2. Résoudre ce problème en utilisant pdetool de Matlab.
3. Tracer T(x)
4. Tracer (x)
5. Déterminer : T(x=e1), T(x=e2)
6. Déterminer la position et la valeur de Tmax.
7. Vérifier par analogie électrique les valeurs des températures T(x=e1) et T(x=e2)

7 / 13