TP: Couche limite dynamique d’un écoulement laminaire sur une
plaque plane
I. Introduction
Lorsqu’un fluide visqueux s’écoule sur une paroi, les vitesses des particules fluides aux niveaux de
cette paroi s’annulent à cause de la condition d’adhérence. La vitesse augmente à mesure que l’on
s’éloigne de la paroi jusqu’à atteindre sa valeur à l’infini amont: cette zone s’appelle couche limite.
La vitesse de l’écoulement sur la plaque 𝑉̅→= 𝑢(𝑥, 𝑦)𝑒→𝑥 où u(x,y≥δ)=U0. Avec : U0 la vitesse initiale
de l’écoulement et δ l’épaisseur de la couche limite.
Des études analytiques montre que [1, 2] :
𝑢 0,32𝑦
= 𝑒𝑟𝑓( )
𝑈0 𝛿
𝜈𝑥
𝛿(𝑥) = 5,66 √
𝑈0
Avec : erf la fonction erreur et ʋ la viscosité cinématique du fluide.
II. Description du problème
Dans ce TP on s’intéresse à étudier la couche limite d’un écoulement 2D laminaire sur une plaque
plane. Cet écoulement va être modélisé en utilisant le système ANSYS-Fluent dans la suite ANSYS-
Workbench 15.0.
Les étapes suivies pour étudier ce problème sont :
La conception de la géométrie du domaine de fluide en utilisant le système « ANSYS-
DesignModeler »
La génération de maillage du domaine de fluide en utilisant le système « ANSYS-Meshing »
Lancer la simulation dans le système « ANSYS-Fluent » après la détermination des :
propriétés physiques de fluide, les conditions aux limites, les modèles numériques de
simulation….
Afficher les résultats (contours de pression, vitesse, …) par « ANSYS-Fluent » et/ou
« ANSYS-CFD-Post»
III. Modélisation numérique
Lancer la suite « ANSYS-Workbench » puis ouvrir le système d’analyse « Mécanique des Fluides
(FLUENT) »
Création du domaine de travail
=========== Dessiner les arêtes du domaine de travail et préciser les dimensions =========
1)
2)
3) et
4)
=========== Créer et diviser la surface (volume en 3D) interne du domaine de travail=========
1) surface à partir d’esquisses (sélectionner « Esquisse 1 » pour l’objet
de base et « ajouter un corps bloqué » pour l’opération)
2)
3)
(Dessiner une ligne verticale qui passe par x=0)
4) Ligne à partir d’esquisses (sélectionner « Esquisse 2 » pour l’objet de
base)
5) (sélectionner « Ligne 1 » pour l’arête et « corps surfacique »
pour la Cible)
=========== Définir la nature du domaine de travail et les conditions aux limites =========
1) (Cliquer sur le bouton droit de la sourie et désactiver le corps
filaire)
2)
3) Sélectionner l’arête d’entrée et cliquer sur le bouton droit :
Top
Inlet
outlet
Plaque
Symmetry
Refaire la même chose afin de définir les conditions aux limites pour les autres arêtes
Maillage du domaine de travail
1) Insertion
(Sélectionner les deux faces du domaine de travail)
2) Insertion
(Sélectionner l’arête d’entrée du domaine de travail)
(Refaire la même chose pour les arêtes verticales en tenant compte du sens d’accroissement d’écart
entre les nœuds)
3) Insertion
(Sélectionner l’arête qui est nommée « symmetry »)
(En tenant compte du sens d’accroissement d’écart entre les nœuds, refaire la même chose pour
l’arête horizontale en haut proche de l’entrée)
4) Insertion
(Sélectionner l’arête qui est nommée « plaque »)
(En tenant compte du sens d’accroissement d’écart entre les nœuds, refaire la même chose pour
l’arête horizontale en haut proche de la sortie)
5)
Maillage raffiné
dans la zone
proche parois
Simulation Numérique et Résultats
=========== Définir le modèle numérique, les conditions aux limites et simulation=========
1)
2)
3)
4)
5)
6)
7)
8) + + + +
9) (enregistrer le ficher XY)
Remplir le tableau suivant :
Dans notre cas : δ(x)th=0,967x0,510-2 m
x 0 0,5 1 1,5 2 2,5 3 3,5 4
δ(x)th 0
δ(x)num