TP n°1 :
Prise en main de Pdetool de Matlab
2-
TP N° 2 : Les ailettes
I. Répartition des températures le long des ailettes triangulaire et rectangulaires.
Considérons deux ailettes de formes distinctes, une triangulaire et l’autre rectangulaire. Les deux sont faites du
même matériau de conductivité thermique (lamda) 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
1. En recourant au logiciel Matlab , on se propose d'étudier seulement la moitié de l’ailette.
Justifier ce choix :
Car l’axe de symétrie est un axe adiabatique, alors il va nous faciliter l’etude de l’ailette par étudier juste la
moitié.
2. Ailette de section triangulaire
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 0.01 0.02 0.03 0.04 0.05 0.06 0.07
T(x, y=0) 115 107 100 94 88 82 76 70
3. Ailette de section rectangulaire
En maintenant les mêmes hypothèses que dans l’ailette triangulaire, 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 0.01 0.02 0.03 0.04 0.05 0.06 0.07
T(x, y=0) 115 107 100 94 90 86 83 81
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 :
Avec : m=h.P/(.S)
On donne : S=1*e et P=2*(e+1)
A partir de la solution analytique (Eq. 1), calculer les températures T(x).
2. 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.
3. 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.
Conclusion :
On peut conclure d’après le tableau ci-dessus que l’effet des fluides influence sur l’efficacité
de l’ailette et laisser l’ailette d’avoir une même efficacité par rapport a n’importe quel effet
de fluide (air ,eau ,..).