MAILLOT Anne-Sophie 41003462 - NOURRY Auriana 40008619
Travaux pratique n°4:
Écoulement sur un plan incliné
Introduction
Au cours de cette séance, il vous est proposé de traiter un autre cas d’étude simple pouvant être
résolu analytiquement. Il s’agit de l’écoulement d’un film de liquide sur un plan incliné. Comme
la séance 1 et 2, il s’agit ici de comparer les résultats analytiques aux simulations numériques.
Présentation de l’étude
Nous étudions l’écoulement d’un fluide incompressible de masse volumique ρ, de viscosité η
avec une profondeur constante h dans un canal de largeur constante L, faisant avec l’horizontale
un angle α.
Approche théorique
1. Déterminer, en régime permanent, le champ de vitesses, de la forme ⃗v = v(x, y) ⃗ex et le
champ de pression p(x, z) en précisant les hypothèses employées.
Hypothèses :
Écoulement stationnaire : Les vitesses ne dépendent pas du temps.
Fluide incompressible : La divergence du champ de vitesses est nulle, c'est-à-dire ∇ · 𝑣 = 0
Écoulement unidirectionnel : Le fluide s'écoule uniquement dans la direction x le long du
→ →
canal, donc 𝑣 = 𝑣(𝑥, 𝑦) 𝑥, où 𝑣(𝑥, 𝑦) est la composante de vitesse dépendant de x et y.
Régime de cisaillement : Le fluide adhère aux parois inférieure et supérieure (condition de
non-glissement), donc la vitesse est nulle à ces limites.
Gravité : Le fluide s'écoule sous l'effet de la gravité, la force de gravité ayant une composante
dans la direction xxx (parallèle au canal), et une composante normale à la surface libre du fluide.
Equation de Navier-Stokes:
→ → 2→ →
ρ(𝑣 · ∇) 𝑣 =− ∇𝑃 + µ∇ 𝑣 + 𝑓
On a alors deux équations:
𝑑𝑃
Le profil de la pression le long de cet écoulement : 𝑑𝑦 =− 𝑔ρ 𝑐𝑜𝑠α → 𝑃(𝑦) = 𝑃0 − ρ𝑔ℎ 𝑐𝑜𝑠α
Le profil de vitesse le long de cet écoulement :
2
𝑑𝑃 𝑑𝑣 ρ𝑔 𝑠𝑖𝑛α 2 2
0 =− 𝑑𝑥
+ µ 2 + ρ𝑔 𝑠𝑖𝑛α → 𝑣(𝑦) = 2µ
(ℎ − 𝑦 )
𝑑𝑦
2. Quel est le débit volumique ? Définir et calculer une vitesse moyenne en fonction des
paramètres de problème en nommant L la largeur de l’écoulement.
Données glycérine :
3
ρ = 1260𝑘𝑔/𝑚
µ=1.49 Pa.s
Le débit volumique Q est défini comme l’intégrale de la vitesse sur toute la section transversale
du canal. Pour un canal de largeur L=10 mm=0.01m.
ℎ
2 ρ𝑔 𝑠𝑖𝑛α 3 −6
𝑄 = ∫ 𝑣(𝑦)𝑑𝑦 · 𝐿 = 3 µ
𝐿ℎ =1. 38 × 10 𝑘𝑔/𝑠
−ℎ
𝑄
𝑣𝑚𝑜𝑦 = 𝐿·ℎ
=0.138 m/s
Travaux de Simulation
Création du projet
— Ouvrir Workbench
— Dans la boîte à outils (Toolbox) à la section Systèmes d’analyse, double-cliquer sur l’onglet
Mécanique des fluides (Fluent)
Création de la géométrie
— Changer les unités pour "Millimètre"
— Dessiner dans le plan XY, un rectangle de hauteur 1 mm et de longueur 10 mm.
— Convertir l’esquisse en objet en 2D
Création du maillage
— Dimensionner le maillage en contrôlant la taille des éléments sur la face. Nous considérerons
une taille d’élément de 1.10−4m
— Nommer les frontières du domaine avec la nomenclature suivante : inlet, outlet, wall, free 4.4
Configuration du solveur Fluent
— General : Cocher la case Gravity et indiquer les valeurs des composantes du vecteur gravité
(9,81 m.s−2 ). On considère dans un premier temps un angle α de 30°.
— Models : viscous laminar
— Material : glycerin
— Boundary conditions :
— inlet composante selon x trouver à la Question 1. Sélectionner Components et expression à
l’aide du menu déroulant à côté du champ à compléter. Préciser les unités des grandeurs entre
crochets.
— outlet : pressure outlet (fixe la pression statique en ce point).
— free : la valeur de contrainte de cisaillement est fixée à 0.
— Modifier la valeur des résidus à 1e-5 pour un meilleur résultat
— Solution Initialization : Hybrid
— Run calculation avec 2000 itérations
Cas 1: Plan incliné à 30°
● Profil de vitesse le long de l’écoulement
● Profil de pression le long de l’écoulement
● Profil de vitesse sur l’épaisseur de l’écoulement
Refaire une simulation en entrant une composante de la vitesse selon x en entrée constante et
également à la vitesse moyenne. Comment évolue la vitesse et la pression le long de
l’écoulement ? Comparer les profils de vitesse sur l’épaisseur de l’écoulement au centre du
domaine de calcul avec celui obtenu précédemment.
Cas 2 :Vitesse constante imposée : 2mm/s
● Profil de vitesse le long de l’écoulement:
● Profil de pression le long de l’écoulement
● Profil de vitesse sur l’épaisseur de l’écoulement
Cas 3 : Vitesse moyenne = 0.138 m/s
● Profil de vitesse le long de l’écoulement
● Profil de pression le long de l’écoulement
● Profil de vitesse sur l’épaisseur de l’écoulement
Les deux graphiques montrent une distribution de vitesse similaire au début, mais avec des
différences importantes. Le cas 1 (vitesse théorique) présente une vitesse maximale initiale
d'environ 1.85 mm/s qui décroît rapidement, stabilisant vers 1.4 mm/s. Le cas 3 montre une
vitesse maximale plus élevée, autour de 2.3 mm/s, avant de se stabiliser à environ 2.15 mm/s. La
différence principale est que le cas 3 a une vitesse plus uniforme et élevée dans toute l'épaisseur,
tandis que le cas 1 présente un pic initial plus bas et une vitesse finale plus faible.
Refaire une simulation en considérant un angle α de 60°. Comparer les profils de vitesse dans
l’épaisseur de l’écoulement pour les deux angles. Est-ce en accord avec la théorie ?
Cas 4: Plan incliné à 60°
● Profil de vitesse le long de l’écoulement
● Profil de pression le long de l’écoulement
● Profil de vitesse sur l’épaisseur de l’écoulement
Ce qui est remarquable au niveau du centre du domaine de calcul, la vitesse augmente d’un coup
jusqu’à la vitesse maximale mais selon les cas on a deux profils :
Pour une vitesse constante, elle va tendre vers la moitié de la vitesse maximale et pour une
vitesse non constante, après le pic elle va décroitre rapidement vers la vitesse constante sur le
reste de la planche.
On reconnaît ici les profils laminaires des vitesses.
