Mcanique des Fluides Numrique
Suite ANSYS - FLUENT
Tutoriel N1: prise en main
Modlisation dun coulement dans une conduite rectiligne
de 10m
Objectif: calculer les pertes de charges liniques
F.OCANA - F.RIMBERT
septembre 2014
Tutorial 1_WB_V2 FRFO_V2014
Page 1
Modlisation d'un coulement simple
Nous cherchons modliser un coulement deau dans une conduite horizontale cylindrique. Ci-dessous
Ci
un
schma de cette conduite en coupe :
10m
0,1 m
y
Ecoulement deau C=2,1 m/s
Lobjectif est de retrouver les valeurs typiques des pertes de charge linique Jlinique de cette conduite.
Donnes techniques sur l'coulement et la conduite
C=2,1m/s
Conduite fonte DN100, diamtre 100 intrieur
Conduite
uite revtue intrieure par centrifugation, ou k =0,1mm
T de l'eau 10C
Masse volumique de l'eau 998,2 kg/m3
Viscosit cinmatique de l'eau 10C: =1,3 e-6 m2/s
Donnes typiques des pertes de charge liniques
lin
Jlinique en pa
par mtre de conduite
D'aprs Colebrook
D'aprs PAM
400
469
Rappel du calcul des pertes de charges linaires
L. C 2 moy
P = . .
en Pa
D
2
avec dtermin
termin partir du diagramme
de Colebrook-Moody
Re= 161500 , /d= 0,001 et 0,18
P=396,9 pa
Tutorial 1_WB_V2 FRFO_V2014
Page 2
Tutorial 1_WB_V2 FRFO_V2014
Page 3
Mise en situation : modlisation d'un coulement simple
Nous cherchons modliser lcoulement deau dans une conduite horizontale cylindrique. Ci-dessous
Ci
un
schma de cette conduite en coupe :
10m
0,1 m
y
Ecoulement deau C=2,1 m/s
Si on suppose le problme symtrique par rapport laxe central, il suffit de modliser la demi-conduite
demi
suprieure. Pour cet exercice, nous ne nous intresserons pas ce qui se passe dans la paroi de la conduite.
Partie modlise
Notre modle peut donc se rduire un simple rectangle 2D.
Donnes techniques sur le
e maillage
Le maillage doit tre structur, idalement de la forme :
wall
velocity
inlet
pressure
oulet
axis
Tutorial 1_WB_V2 FRFO_V2014
Page 4
Crer le projet sous Workbench
Dmarrer ANSYS Workbench.
Workbench Workbench est une plateforme sur laquelle est construite la suite
de produits de simulation numrique.
numrique
Enregistrez votre projet immdiatement (menu Fichier/Enregistrer sous ). Pensez
enregistrer rgulirement votre travail afin de ne rien perdre en cas de problme.
Glissez-dposez
dposez le systme danalyse Mcanique des fluides (FLUENT) depuis la bote
outils vers le schma de projet.
Renommer le projet en conduite de 10m (si ncessaire double-cliquer
cliquer sur le nom de projet
pour le modifier)
Tutorial 1_WB_V2 FRFO_V2014
Page 5
Crer la gomtrie sous Design Modeler
Pour crer le modle de gomtrie, double-cliquer sur Gomtrie . Attendre le dmarrage de
Design Modeler
Au dmarrage de Design Modeler, vrifier que le systme dunit de
longueur est bien rgl sur Mtre , puis valider par OK
Linterface se prsente comme ci-dessous :
Tutorial 1_WB_V2 FRFO_V2014
Page 6
Dans larborescence, slectionner le plan XY (important en 2D pour avoir accs toute les
options), puis cliquer sur le bouton nouvelle esquisse
Afin que le plan desquisse se confonde avec le plan de lcran, cliquer sur laxe Z en bas
droite de lespace de travail.
Tutorial 1_WB_V2 FRFO_V2014
Page 7
Dans longlet esquisse de la bote outils, activer loutil Dessiner/Rectangle puis tracer
un rectangle dans lespace de travail. Les dimensions de ce rectangle nont pas dimportance pour
linstant, mais il faut que le coin infrieur gauche soit confondu avec lorigine du repre.
2
3
1
En cas derreur, pour supprimer des lments tracs, il est possible dutiliser :
la fonction
, pour annuler la dernire fonction utilise
loutil Nouvelle slection
la touche Suppr du clavier.
, puis de slectionner les lments effacer, puis dappuyer sur
Pour dimensionner correctement notre rectangle, activer loutil Cotes/Gnral de la bote
outils. Cliquer sur le segment suprieur du rectangle, puis dposer la cote en cliquant lendroit
souhait. De mme, coter le segment de gauche.
Tutorial 1_WB_V2 FRFO_V2014
Page 8
2
3
6
5
Dans la Vue de dtails , modifier les dimensions V et H en 10m et 0,05m respectivement. Le
trac du rectangle est automatiquement mis jour avec ces dimensions.
Il est trs probable que le rectangle affich lcran devienne trs petit. Pour recadrer la vue, il suffit
de cliquer sur loutil Ajuster
dans la barre doutils.
Pour dplacer les cotes, il faut activer loutil Cotes/Dplacer dans la bote doutils.
Tutorial 1_WB_V2 FRFO_V2014
Page 9
Modle physique de lcoulement : Design Modeler (suite)
Cliquer sur longlet Modlisation , puis cliquer sur le menu Concept/Surfaces partir
desquisses .
Slectionner lesquisse dans larborescence, cliquer sur Appliquer dans le vue de dtails, et
enfin cliquer sur Gnrer
dans la barre doutils. Le rectangle se grise, signe que le
logiciel lidentifie comme un corps matriel surfacique.
Tutorial 1_WB_V2 FRFO_V2014
Page 10
2
Dvelopper 1 Pice, 1 Corps dans larborescence en cliquant sur le signe . Slectionner le
Corps surfacique et modifier sa proprit Fluide/Solide Fluide dans la vue de
dtails.
Tutorial 1_WB_V2 FRFO_V2014
Page 11
3
A prsent nous allons dfinir ce qui se passe aux frontires de notre corps surfacique. Activer le
filtre de slection Artes
dans la barre doutils.
Slectionner larte suprieure du rectangle. Puis cliquer sur le bouton droit et choisir Slection
nomme
dans le menu contextuel.
Dans la vue de dtails, modifier la proprit Slection nomme en wall . Le mot clef
wall sera plus tard reconnu par Fluent et interprt comme une paroi (la paroi de notre
conduite).
Tutorial 1_WB_V2 FRFO_V2014
Page 12
Cliquer ensuite sur 1 Arte dans la proprit Gomtrie et cliquer sur Appliquer .
Valider cette modification en cliquant sur
dans la barre doutils.
La Slection nomme est alors cre dans l'arborescence
Rpter les 3 tapes prcdentes, partir de la slection d'artes, et nommer les artes gauche,
droite et infrieure respectivement en velocity_inlet , pressure_outlet et axis (ces
noms sont galement constitus de mots clefs).
EN CAS DE BUG: slectionner dans l'arborescence, la slection nomme qui pose problme et
supprimer grce au menu contextuel (touche droite)
Notre modle est termin et prt tre maill. Vous pouvez quitter Design Modeler avec le menu
Fichier/Fermer DesignModeler . Le modle sera automatiquement enregistr, et une fois de retour
sur Workbench une coche
montre que la gomtrie est oprationnelle.
Tutorial 1_WB_V2 FRFO_V2014
Page 13
Crer un maillage sous Meshing
Dmarrer ANSYS MESH, double
ouble-cliquer sur Maillage . Attendre le dmarrage de Meshing
Raliser un bon maillage
Comme nous le verrons par la suite, le maillage doit tre structur afin de faciliter les calculs et permettre
perm
la convergence rapide de la solution.
Slectionner dans l'arboresence Maillage , puis insrer un Dimensionnement
. Activer le filtre de slection Artes
Tutorial 1_WB_V2 FRFO_V2014
dans la barre doutils.
Page 14
Caractristiques du maillage cot entre (velocity inlet)
Slectionner larte gauche ( l'entre en vitesse). Puis cliquer dans la vue de dtail sur
Appliquer
entrer les caractristiques suivantes
A l'aide du ZOOM, S'assurer que la taille de la premire maille est de l'ordre de 0,8mm (modifier le
facteur d'cart , ou la taille global des lments si ncessaire) . Passer le facteur d'cart 6
Caractristiques du maillage cot sortie (pressure outlet)
Tutorial 1_WB_V2 FRFO_V2014
Page 15
De la mme faon, slectionner l'autre arte et entrer les paramtres suivants
Caractristiques du maillage cot mur (wall) et axe (axis)
Slectionner les artes suprieure (wall) et infrieure (axis) avec la touche Ctrl, puis cliquer dans la
vue de dtail sur Appliquer
entrer les caractristiques suivantes
Tutorial 1_WB_V2 FRFO_V2014
Page 16
Gnrer le maillage
Slectionner l'icone Gnrer un maillage
cliquer sur l'icone maillage dans l'arbre , le maillage apparat
Constatation: le maillage n'est pas trs cohrent, il faut le Mapper r si la surface est mappable.
Mapper une surface
Reprer les surfaces mappables
Les faces mappables passent en vert
Tutorial 1_WB_V2 FRFO_V2014
Page 17
insrer alors "une face mappe"
Slectionner la surface mapper
Gnrer le maillage
cliquer sur l'icone
maillage dans
l'arbre , le maillage
apparat
Enregistrer le fichier
Tutorial 1_WB_V2 FRFO_V2014
Page 18
Qualit du maillage /statistiques
Remarque:
Ce maillage est d'assez bonne qualit.
L'influence du maillage est importante en particulier pour les cas plus complexes.
Nous verrons par la suite comment modifier le maillage
Tutorial 1_WB_V2 FRFO_V2014
Page 19
Configurer le problme sous FLUENT
Reprendre le maillage tabli avant les modifications.
A faire en priorit et dans l'ordre
2D Axisymetrique
Vrifier l'chelle
Checker le modle
Confirmer la qualit du maillage
Dsactiver l'affichage(display) du
maillage (trop gourmand en calculs)
Tutorial 1_WB_V2 FRFO_V2014
Page 20
Aide mmoire de configuration - ANSYS Fluent
TP 1: Ecoulement turbulent en conduite cylindrique
Memento
Menu de configuration (Fluent)
Problme axisymtrique 2D.
Modle visqueux : k-.
Utilisation de llment water-liquid de la bibliothque
standard de Fluent.
Le corps surfacique est de type fluide et constitu deau.
La vitesse dentre du fluide est de 2,1m/s ;
La pression en sortie est de 101300 Pa (Constante) ;
En entre et en sortie, lintensit de la turbulence est
de 3,57% et le diamtre hydraulique quivalent est de
0,1m ;
La rugosit de la paroi a une hauteur quivalente de
0,1mm.
La pression de rfrence est de 101300 Pa ;
Critre de convergence des rsidus : aucun : arrt la main
2000 itrations maximum
La maille en bord de paroi doit avoir une hauteur de environ H=0,8 mm
taille moyenne des mailles l=0,003m=3mm
Intensit turbulente I=3,57%
Tutorial 1_WB_V2 FRFO_V2014
Page 21
Complter le SETUP du problme, de haut en bas
Tutorial 1_WB_V2 FRFO_V2014
Page 22
Tutorial 1_WB_V2 FRFO_V2014
Page 23
Tutorial 1_WB_V2 FRFO_V2014
Page 24
Tutorial 1_WB_V2 FRFO_V2014
Page 25
Tutorial 1_WB_V2 FRFO_V2014
Page 26
Tutorial 1_WB_V2 FRFO_V2014
Page 27
Analyser les rsultats sous FLUENT
L'analyse des rsultats s'effectue sous Fluent dans la partie Results
L'objectif est ici de calculer et mettre en forme les donnes issues de la simulation:
Dtails des onglets :
L'onglet Graphics and Animations permet de crer
- des reprsentations graphiques (des champs des
vitesses , de pressions ) sous forme de cartes colores
- des reprsentations graphiques (des champs des
vitesses , de pressions ) sous forme de cartes vectorielles
- des animations graphiques des coulements sous forme
de video :
A utiliser avec prcaution, car cela ncessite une trs
grande puissance de calcul qui peut bloquer la RAM
L'onglet Plots permet de crer
- des graphiques X,Y des cartes vitesses, champ de pression,
valeurs des Y+
Tutorial 1_WB_V2 FRFO_V2014
Page 28
L'onglet Reports permet de calculer
- des flux, des efforts, des Cx, des pertes de
charges, des moyennes de vitesses ou de
pression
Tutorial 1_WB_V2 FRFO_V2014
Page 29
A/ Visualiser les champs de vitesses (ou de pression)
- "Tracer le champ des vitesses axiales dans la conduite ( Graphics / Contour)
B/ Tracer le graphique d'volution de la vitesse au niveau de l'axe de la conduite
Utiliser la fonction (Plot / XY plot)
- Conclusions:
Le champ des vitesses dans la conduite met une certaine distance avant de se stabiliser
Il n'est donc pas intressant de mesurer la perte de charge J entre l'entre et la sortie, car le champ de
vitesse influence le champ de pression
Il convient donc de dterminer la position dans la conduite partir de laquelle le champ de vitesse est
constant
Tutorial 1_WB_V2 FRFO_V2014
Page 30
C/ Tracer les cartes des vitesses diffrentes abscisses
- Tracer les 2 cartes de vitesses en entre et en sortie (Plot / XY plot)
- Crer 3 nouvelles "sections - line " 1m, 2,5m et 5m
en cliquant sur "New surface"
- Tracer les 5 cartes superposes.
- Conclusions:
Le champ des vitesses devient constant partir de l'abscisse 5m
Tutorial 1_WB_V2 FRFO_V2014
Page 31
D/ Calculer des pertes de charges dans la conduite
- Calculer les pressions statiques aux diffrentes abscisses (Reports / Surface Intgrals/Mass-Weighted
Average)
- Calculer le deltaP entre 5m et la sortie, en dduire la perte de charge linique. Comparer aux valeurs de
rfrences
Remarque: le calcul par Area-Weighted Average fourni le mme rsultat, ce qui n'est pas toujours le cas
Tutorial 1_WB_V2 FRFO_V2014
Page 32
E/ Calculer du dbit massique (Report / Fluxes /Mass Flow Rate)
Le dbit dtermin par simulation est de 16,46 kg/s,
Le dbit thorique est de 16,46 kg/s (qm=.C.S)
Remarque:
Il est important que les dbits d'entre et de sortie soit strictement identiques et exacts.
Ceci est un indice de la bonne qualit de la simulation
Conclusions
Les simulations Fluent donnent des valeurs cohrentes avec les valeurs de rfrence.
Nous verrons dans les prochains exercices que le choix du maillage et du modle de turbulence ont des
influences non ngligeables, qui ne sont pas ici dcelables (cas trop simple).
Tutorial 1_WB_V2 FRFO_V2014
Page 33
