TP2IFA

UNIVERSITE CHOUAIB DOUKKALI
ECOLE NATIONALE DES SCIENCES APPLIQUEES D’EL JADIDA
Filière : Génie Energétique et Electrique (G2E)
TP :
Modélisation numérique 2
Réalisé PAR :
-LOUMIR FADOUA
- MANSSAR AMINA Encadré PAR

-MELIANI ILHAM - ABALA ILHAM
Buts de ce TP :
Ce TP vise à étudier le transfert de chaleur convectif. Le premier exemple concerne le transfert
thermique au travers d’une fenêtre impliquant la convection naturelle et la convection externe. Pour la
résolution numérique du problème, des notions du cours Calcul Scientifique pour ingénieur seront
utilisées et implémentées sur le logiciel MATLAB. Le deuxième exemple concerne le design thermique
d’un séchoir à cheveux.
Mise en situation
Une fenêtre est un élément très important de l’enveloppe d’un bâtiment car une fraction importante des
déperditions thermiques y sont dues. La fenêtre étudiée comporte deux plaques de verre de 3 mm
d’épaisseur séparées d’une distance L, tel qu’illustré sur la Figure 1. Le gaz dans l’espacement entre les
fenêtres est de l’air à pression atmosphérique.
La largeur de la fenêtre (perpendiculaire à la figure) est de W=2 mètres.
PARTIE THÉORIQUE :
On retrouve dans une fenêtre un transfert de chaleur par conduction dans le verre et de la convection
dû au déplacement de l’air sur les parois intérieures et extérieures, ainsi que dans l’espace entre les deux
vitrages.
2
Bilan d’énergie sur les surfaces 1 et 4
Le bilan d’énergie sur une surface permet d’exprimer le flux de chaleur en fonction des températures en
jeu. Par exemple, pour la surface 4 (figure 1), soit la paroi extérieure de la fenêtre, nous obtenons
l’expression suivante :
Afin de déterminer le coefficient de convection moyen h , nous allons utiliser sur la surface 4 les
corrélations pour la convection externe présentées au chapitre 7. Celles-ci impliquent de calculer
d’abord le nombre de Reynolds donné par :
où la viscosité du fluide ν est évaluée à la température du film
Nous supposerons que le vent de vitesse U∞ est horizontal et parallèle à la

fenêtre, tel qu’illustré ci-dessous.
3
4
5
global Lv L kv Tint Text H g h1 h23 h4 Uext W;
%Initialisation des variables%

Lv=0.003;
L=0.01;
6
kv=1.4;
Tint=273+20;
Text=273+35;
H=1.5;
g=9.81;
Uext = 20;
W=2;
%Paramètres pour la méthode numérique sysnl

nmax=20;
eps=1.0e-6;
hi=1.0e-3;
%Vecteur de départ
Tini=[400 400 400 400];
T=sysnl('sys',Tini,nmax,hi,eps,'resul.dat');
T1 = T(size(T,1),1);
q1 = H*W*h1*(T1-Tint)
T4 = T(size(T,1),4);
q4 = H*W*h4*(Text-T4)
Fixons L=1cm et on fait varier Lv et on voit les valeurs des

températures.
- -pour Lv=3mm
T1 T2 T3 T4
400 400 400 400
298.7676 298.8294 307.3614 307.4232
297.9876 298.0401 307.4478 307.5012
297.9873 298.0407 307.4478 307.5013
Calculons
bilan d’énergie pour les surfaces 1 et 4
7
Pour les surfaces 2 et 3 :
8
- La courbe associée
9
- -pour Lv=0.01mm
T1 T2 T3 T4
400 400 400 400
298.8154 298.8156 307.4183 307.4185
298.0234 298.0236 307.4975 307.4977
298.0231 298.0233 307.4975 307.4977
- La courbe associée
10
- -pour Lv=6mm
T1 T2 T3 T4
400 400 400 400
298.7205 298.8431 307.3054 307.4279
297.9521 298.0583 307.3987 307.5048
297.9518 298.0580 307.3987 307.5048
- La courbe associée :
11
- Interprétation :
On remarque que les températures avec la variation de Lv les
températures varient d’une façon non considérable.
- MISE EN SITUATION 2
Le séchoir courant a une conduite de sortie ayant un diamètre interne de D = 4.6cm. L’élément
chauffant est composé d’un fil de nichrome et transmets à l’air ambiant une puissance de q = 1000 W. La
différence de potentiel appliquée à l’élément chauffant est de 110 volts. De plus, vous savez que
l’appareil aspire l’air ambiant à une température de Ti = 20°C et chauffe cet air jusqu’à une température
de 50°C.
12
 Pour D=0 .05 et Pour q variable
q(W) 800 900 1000

T(K) 664.67 730.9 799 .61
Dans ce cas on remarque lorsque on augmente q la température augmente aussi .
 Pour q=900 et D variable
q(W) 0 .06 0.05 0.04

T(K) 795.67 730.9 659.59
Dans ce cas on remarque que lorsque on diminue D la température diminue aussi
CONCLUSION :
À laide de ce TP on a pu bien comprendre le phénomène de convection de toutes ces formes
13

TP2IFA

Transféré par

Droits d'auteur :

Formats disponibles

TP2IFA

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TP2IFA

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITE CHOUAIB DOUKKALI

ECOLE NATIONALE DES SCIENCES APPLIQUEES D’EL JADIDA

Filière : Génie Energétique et Electrique (G2E)

- MANSSAR AMINA Encadré PAR

La largeur de la fenêtre (perpendiculaire à la figure) est de W=2 mètres.

où la viscosité du fluide ν est évaluée à la température du film

Nous supposerons que le vent de vitesse U∞ est horizontal et parallèle à la

%Initialisation des variables%

%Paramètres pour la méthode numérique sysnl

Fixons L=1cm et on fait varier Lv et on voit les valeurs des

q(W) 800 900 1000

Dans ce cas on remarque lorsque on augmente q la température augmente aussi .

 Pour q=900 et D variable

q(W) 0 .06 0.05 0.04

Dans ce cas on remarque que lorsque on diminue D la température diminue aussi

À laide de ce TP on a pu bien comprendre le phénomène de convection de toutes ces formes

Vous aimerez peut-être aussi