Ministre de l’Enseignement Supérieur et de la Recherche scientifique

UNIVERSITÉ MOULOUD MAMMERI TIZI-OUZOU

Faculté du Génie de la Construction

Département de Génie Mécanique

Spécialité : Energétique

COMPTE RENDU
TP TTH
ECOULEMENT ET TRANSFER DE CHALEUR SUR UNE PLAQUE
PLANE

RÉALISÉ PAR : ENCADRÉ PAR :

 ZEGGAN WALID M. BENFERHAT.S

 SI HADJ MOHAND OUARDIA

Année universitaire : 2022/2023

1
  Introduction

Pour cette étude, nous allons nous concentrer sur la simulation d’une plaque
plane afin d'illustrer la solution d’un écoulement laminaire 2D sur cette dernière.
L’intérêt de l’étude est d’observer le développement de la couche limite
thermique, et d’observer également comment le nombre de Nusselt varie dont
le calcul s’effectue à partir du flux de chaleur `a la paroi.

Objectif :

L’objectif de ce TP est de s’initier à la manipulation des deux logiciels de CFD,

Gambit (le meilleur) et Fluent (le solveur et le post-processeur) à travers l’étude
de Écoulement et transfert de chaleur sur une plaque plane, Un écoulement de
fluide incompressible peut être permanent ou non (instationnaire). La distance au
cours de laquelle la couche limite thermique évolue correspond à la longueur de
la plaque. La notion d’écoulement permanent est une notion temporelle.

Brève présentation des logiciels :

Gambit et Fluent sont des logiciels, sous licence commerciale, permettant de

réaliser des simulations 2D ou 3D en mécanique des fluides et transfert de
chaleur allant de la construction du maillage avec Gambit à la résolution des
équations de Navier Stokes et au post-traitement avec Fluent. Largement
répandus dans l’industrie (automobile, aéronautique, espace, etc…) en raison de
leur interface graphique puissante et de l’abondance de leurs options, ils
permettent de réaliser des simulations sur tous types de géométries complexes
(fixes ou mobiles) associées à des maillages fixes ou adaptatifs et avec des
modèles physiques variés (diphasique, turbulent, etc…).

2
Interface du logiciel Gambit :

Figure-1 : Interface du logiciel Gambit.

Interface du logiciel Fluent :

On peut démarrer 4 versions de Fluent 2D, 3D, 2DDP, 3DDP ayants la même
interface. (Figure-2)

Figure-2 : Versions du logiciel Fluent.

 2D ( 2 Dimensions)
 2DDP ( 2 Dimensions Double Précision)
 3D ( 3 Dimensions)
 3DDP (3 Dimensions Double précision)

3
 Figure-3 : Interface du logiciel Fluent

Description du problème :

Dans ce problème, nous considérons une plaque plane de longueur L = 1 m

(Figure1). La vitesse du flux d’air est de sorte que le nombre de Reynolds Re =
1e+05. La plaque est maintenue à une température T = 400 K.

Figure-4 : Schéma du problème

4
  Les étapes qu’on a suivi durant le TP :
Ι. Gambit - Génération de la géométrie (plaque plane) et du domaine de
calcul.

Afin de générer un maillage structuré (rectangulaire), cette partie nous montre

comment résoudre un écoulement laminaire ou turbulent sur une plaque plane
par le biais du logiciel Fluent.

Figure-5 : Domaine du calcul

Pour créer une géométrie sur Gambit, il faut suivre les instructions suivantes :

1- Créer un nouveau dossier

2- Lancer Gambit (choisir le nom du fichier crée), puis cliquer sur
« browse ».
3- Choisir « commande de géométrie »
4- Dans la zone « Create Real Vertex-Global », créer le point (0,0,0) pour les
coordonnées x, y, z respectivement.
5- Créer les autres points avec les coordonnées correspondantes.
6- Créez des segments à partir de ces points pour définir le domaine de
calcul.
7- Générer une face à partir des segments créés :
8- Définir les points de nœuds le long des bords :

 Dans Opération, bouton de commande Mesh - bouton de

commande Edges.
 Sélectionner le segment gauche et droit (les deux ensembles)
  Changez l'option d'espacement d’ « Interval Size » à « Interval
Count ».
 Effectuer les commandes nécessaires…
 Fermer la fenêtre Mesh Edges.
9- Spécifier les types de frontières sur tous les bords tels spécifiés…
10- Générer le maillage sur la surface du domaine.
11- Enregistrer le maillage dans le format « .msh » ensuite quittez
Gambit…

𝚰𝚰. FLUENT
Préparation

1- Copiez les fichiers suivants dans votre dossier de travail :

• plate.msh
• point-x.jou
• iso-x.jou

2. Démarrer le solveur 2D double-précision (2ddp) de FLUENT

Configuration et solution

Etape 1 : Grille

1-Lire le fichier de
maillage plate.msh.
File → Read
→Case...
FLUENT va lire le
fichier de maillage
et rapporter la
progression dans la
console.

6
2- Vérifier la grille :
Grid →Check

3- Vérifiez l'échelle
de la grille.
Grid →Scale...

4- Affichez la grille
(Figure-6).
Display →Grid...

• Ouvrir la
case ‘’Colors ‘’
ensuite sélectionner
« Color by ID »
• Cliquez sur
« Display » et
fermez le panneau
Grid Display.

Figure-6 : Affichage de la grille

7
Etape 2 : Modèles

1- Conservez les
paramètres par défaut Et
du solveur. Define → ap
Models →Solver... e
3:
2- Activez l'équation M
d’énergie. Define → at
Models →Energy..
éri
au
x

Define →Materials...

• Conserver les
propriétés par défaut
de l'air.

Etape 4 : Conditions aux limites

Define →Boundary
Conditions…
1- Définirles conditions
aux limites…
• Entrer 1,4607 m/s pour
la grandeur de la vitesse.

8
2- Conserver les
conditions limites par
défaut pour la pression de
sortie ‘’PressureOutlet’’

3- Définir la condition aux

limites pour la paroi
‘’Wall’’ en entrant la
température à 400K.

Etape 5 : Solution

1- Set the solution controls. • ItemInitialize the flow.

Solve → Controls →Solution... Solve → Initialize →Initialize...

2- Activer le traçage (plot) des résidus pendant le calcul.

Solve → Monitors →Residuals...

9
 3- Enregistrer le
fichier de cas
(plate1.cas.gz). File →
Write →Case...

4. Commencez le
calcul pour 1000
itérations
Solve →Iterate..

Figure 7 : Scaled Residuals

5-Créer des iso-surfaces.

Surface →Iso-Surface.
• Entrer 0,1 pour « Iso-
values »
• Entrer 0,1 pour « New
surface name ».
Pour créer les 9 iso-
surfaces à des positions x
différentes ;
• Lire le fichier iso-x.jou. File
→ Read →Journal...
• Sélectionnez le fichier iso-
x.jou et cliquez sur OK.

10
6. Créez des points sur la
plaque correspondant aux
emplacements des iso-
surfaces.
Surface →Point...

(x=0.1 pour le point1)

(x=0.2 pour le point 2)…

• Lire le fichier point-x.jou.

File → Read →Journal...
• Sélectionner le fichier
point-x.jou et cliquez sur
OK.
7- Enregistrez les fichiers
case and data
(plate2.cas.gz et
plate2.dat.gz).
File → Write →Case &
Data...

Étape 6 : Post-traitement

1. Afficher les vecteurs de

vitesse (Figure 8).
Display→Vectors..
• Sélectionner « Velocity
dans vectors of »
• Sélectionner « Velocity
et « Velocity Magnitude »
dans « Color by » et
cliquer sur Display.
• Fermer le panneau
Vecteurs

11
 Figure 8 : Vecteurs de vitesses

2. Désactiver les titres

(titles et colormap).
Display →Options...

• Cliquer sur Appliquer et

fermez le Display Options
panel.
• Zoomer la vue près de la
plaque pour voir le profil
de vitesse comme indiqué
dans la Figure 9

12
 Figure 9 : Vue agrandie des vecteurs de vitesse

3- Afficher les contours de

la température (Figure 10).
Display→Contours...
• Sélectionner
« Température » et
« Température statique »
dans « Contours of »
• Activer l'option « Filled »
dans « options ».
• Cliquer sur Afficher.
• Fermer le panneau
Contours.

Figure 10 : Contours de la température statique

13
  Les calculs effectués
1- Calculer le nombre de Nusselt en utilisant la corrélation suivante :
Nu = 0.332 √𝑹𝒆Pr1/3

On prend : Le nombre de Prandtl Pr pour l'air est de 0,744176

On a : Le nombre de Reynolds Re est calculé en fonction de l'emplacement du

point sur la plaque
𝛒𝐕𝐋(𝐱)
𝐑𝐞 =
µ

Où : 𝛒 : Densité du fluide en [kg/m3]

𝐕 : Vitesse du débit en [m/s]

(𝐱) : Dimension linéaire caractéristique en [m]

µ : Viscosité dynamique du fluide en [Pl]

Le tableau de Nu(corrélation) du point 1 à 10.

Le La
point
3
position ρ[kg/m ] µ [Pl] Pr1/3 Re
√𝑅𝑒
Nu corrélation
[m]
1 0.1 1.225 1.7894*10-5 0.9061 9999.76249 99.99 30.08
2 0.2 1.225 1.7894*10-5 0.9061 19999.95249 141.41 42.53
3 0.3 1.225 1.7894*10-5 0.9061 29999.28747 173.20 52.10
4 0.4 1.225 1.7894*10-5 0.9061 39999.04996 199.99 60.16
5 0.5 1.225 1.7894*10-5 0.9061 49998.81245 223.60 67.26
6 0.6 1.225 1.7894*10-5 0.9061 59998.57494 244.94 73.68
7 0.7 1.225 1.7894*10-5 0.9061 69998.33743 264.57 79.58
8 0.8 1.225 1.7894*10-5 0.9061 79998.09992 282.83 85.08
9 0.9 1.225 1.7894*10-5 0.9061 89997.86241 299.99 90.24
10 1 1.225 1.7894*10-5 0.9061 99997.6249 316.22 95.12
Exemple de calcul :

Sachant que :
 ρ = 1.225 kg/m3
 V = 1.4607 m/s
 L(1) = 0.1 m
 µ = 1.7894 * 10-5 Pl

14
On aura :
(1.225∗1.4607∗0.1)
Re(1) = = 9999.76249
(1.7894∗10−5)

D’après ces résultats on obtient :

Nu(1)= 0.332 √9999.76249 (0,744176)1/3

Nu(1)= 30.0821

2-Calculer le nombre de Nusselt en utilisant les paramètres de FLUENT

suivants :

𝐪𝐩𝐗𝐩
𝑵𝒖 = (𝑻− 𝑻𝒑)𝝀

qp = Flux de chaleur au point p en [w/m2]

Xp = L'emplacement du point p, sur la plaque en [m]

T = Température de la plaque (400 K).

TP = Température du fluide au point p. en [K]

𝝀 = conductivité thermique du fluide en [w/ m-k]

Le tableau de Nu (fluent) du point 1à 10 :

Le point La position La température Le flux thermique de Nu (FLUENT)

[m] de position [K] position [w/m2]
1 0.1 300.16199 723.09045 29,92
2 0.2 300.22937 513.1322 42,5
3 0.3 300.2811 419.48712 52,15
4 0.4 300.32477 363.43494 60,26
5 0.5 300.36334 325.06693 67,48
6 0.6 300.39822 296.66501 73,84
7 0.7 300.43039 274.51437 79,75
8 0.8 300.46033 256.56198 85,21
9 0.9 300.48849 241.54576 90,27
10 1 300.5159 228.65506 94,96

15
 Exemple de calcul :

Sachant que :

 q1= 723.09045 w/m2

 L1=0.1m

 T= 400 K

 T1=300.16199 K
 λ = 0,0242 w/m-k
donc :
(723.09045∗0.1)
Nufluent(1)=(
400 – 300.16199)∗0.0242

Nu fluent(1) = 29.92

 Les graphiques :
• La courbe de Nu Fluent en fonction de la température T :

100
•
90
Nu(fluent)=F(Tb)
80
100.00
70
80.00
60
60.00
Nu(F)

50
40 40.00
Nu (FLUENT)
30 20.00
20
0.00
10 300.16 300.23 300.28 300.32 300.36 300.40 300.43 300.46 300.49

0 Tb[K]

Nu fleunt

16
• La courbe de Nu fluent en fonction de racine Re :

Nu=F(√Re)

100
100
80
90
60
Nu (c)

80
40
70
20
60
0
50 99,99 141,41 173,2 199,99 223,6 244,94 264,57 282,83 299,99
Nu (FLUENT)
40 √Re Nu (corrélation)
30
Nu corrélation
20

10

0
• La courbe de Nufluent et Nucorrélation en fonction de la position P :
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Nu(f),Nu(cor)=f(p)
100.00

80.00

60.00

40.00

20.00

0.00
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
position m

Nu fleunt Nu corrélation

17
Commentaires :

 Selon le premier graphique Nu fluent=f(TP), on remarque que :

• Le nombre de Nusselt varie en fonction de la Température c'est-à-dire
quand la température de position augmente, le nombre de Nusselt
augmente.

 Selon le deuxième graphique Nu fluent, Nucorrélation =f(P), on remarque :

•Nufluent et Nucorrélation augmentent en fonction de la position

•La courbe de Nufluent se superpose sur la courbe de Nucorrélation car les

valeurs théoriques et expérimentales sont presque les mêmes, ce qui nous
permet de Valider notre simulation sur fluent, et de dire que ce dernier est
un logiciel très fiable pour tout calcul numérique.

18
19