TP : La similitude numérique dans ANSYS Workbench:
Une approche pratique.
ABDELAZIZ EL HACEN
Maria Merouane
Introduction :
L’importance de la similitude numérique en mécanique des fluides est considérable. Elle permet de remédier à
plusieurs contraintes qu’on rencontre comme : réduire la complexité des problèmes tout en identifiant les nombres
sans dimension pertinents ce qui nous permettra de réduire le nombre de paramètres à étudier et à simplifier
l’analyse, généraliser les résultats obtenus pour un cas particulier à un étendu de classe plus large en utilisant les
nombres adimensionnels, optimiser les expériences physiques en réduisant les couts et les délais en concevant des
expériences à échelle réduite au lieu d’étudier un phénomène observé à grande échelle directement, valider les
modèles numériques en comparant les résultats numériques avec des expériences réalisées dans des conditions
similaires et enfin étudier des phénomènes difficiles à observer expérimentalement et avoir accès à des informations
détaillées sur l’écoulement.
En hydraulique, on constate que cette similitude joue un rôle crucial ; elle permet de modéliser des écoulements
complexes comme des écoulements turbulents, des écoulements à surface libre mais aussi des interactions,
optimiser les ouvrages hydrauliques avant même leur conception et parvenir à évaluer leur comportement dans
différentes conditions, étudier les impacts environnementaux et enfin réduire les coûts expérimentaux car la
simulation numérique permet de réaliser des expériences virtuelles à moindre coût.
Parmi les logiciels les plus utiliséss en hydraulique, on s'intéressera à ANSYS Workbench.
La similitude numérique dans ANSYS Workbench: Une approche
pratique.
La similitude numérique, en particulier dans un contexte d'ingénierie comme celui d'Ansys Workbench, consiste à
établir des correspondances entre un modèle numérique et un système réel. Cette correspondance se fait à travers
des nombres adimentionnels qui permettent de comparer des phénomènes physiques différents tout en conservant
les mêmes caractéristiques relatives.
La similitude dans ANSYS WorkBench permet de réduire la complexité tout en réduisant le temps de calcul,
généralisation des résultats, optimisation des simulations et enfin validation du modèle en comparant les résultats
numériques avec des données expérimentales pour enfin valider la fiabilité du modèle numérique.
Mise en œuvre dans Ansys Workbench
1. Définition du problème physique:
• Identifier les phénomènes physiques dominants(écoulement, transfert thermique, déformation, etc.).
• Définir les géométries, les matériaux et les conditions aux limites.
2. Choix des nombres adimentionnels:
• Sélectionner les nombres sans dimension pertinents en fonction du problème physique (Reynolds, Froude,
Mach, etc.).
3. Configuration du modèle numérique:
• Définir les propriétés des matériaux (densité, viscosité, conductivité thermique, etc.).
• Choisir les shémas de discrétisation spatiale et temporelle.
• Définir les conditions aux limites (pressr, vitesse, température, etc.).
4. Maillage:
• Créer un maillage adapté à la géométrie et aux phénomènes physiques étudiés.
• Raffinement du maillage dans les zones où les gradients dont importants.
5. Lancer la simulation:
• Vérifier la convergence de la solution.
• Analyser les résultats.
6. Validation et interprétation:
• Comparer les résultats numériques avec des données expérimentales ou des résultats analytiques.
• Interpréter les résultats en fonction des nombres adimentionnels.
Les conditions aux limites de notre modèle sur l'écoulement autour
d'un cylindre:
Sur le cylindre (paroi solide):
• Condition de non-glissement: La vitesse du fluide à la surface du cylindre est nulle.
À l'entrée du domaine de calcul :
• Profil de vitesse uniforme: On impose généralement un profil de vitesse uniforme à
l'entrée V= 1m/s
À la sortie du domaine de calcul:
• Condition de pression nulle ou de gradient de pression nul: On impose généralement
une condition de pression nulle ou de gradient de pression nul à la sortie.
Type d'écoulement:
• Écoulement incompressible et non turbulent.
Nature du liquide:
1/ Liquide visqueux.
2/ Liquide non visqueux.
1/Pour un liquide visqueux:
Voici les étapes à suivre pour ce modèle:
1•1/Géométrie:
Ouvrir la fenêtre de Ansys Workbench
Se diriger vers Système de composants
Translater la géométrie vers la fenêtre du shéma du projet
Cliquer sur la géométrie
Changer le type d'analyse en 2D.
Double cliquer sur géométrie et une fenêtre s'affiche géométrie-DesignModeler
Ajuster le plan en planXY
Cliquer sur esquisse
Cliquer sur réglages et cocher sur afficher la grille en 2D et et alligner sur la grille
Régler l'espassement principal de la grille
Régler l'échelle de 0 à 10m
Aller sur dessiner
Choisir rectangle, le dessiner sur le plan
Choisir cercle le dessiner sur le plan allant du centre de l'axe
Aller sur cotes et définir les dimentions du dessin
Cliquer sur général pour faire sortir le diamètre du cercle et aller dans la barre de vue en détail et définir le diamètre
D=D5=2m
Cliquer sur horizontale et faire sortire la longueur allant du centre jusqu'à la paroi verticale H2=H1=6m
Cliquer sur verticale et faire sortir les deux distances allant de l'axe des x aux parois verticales V4=V3=6m
Aller sur concept et choisir surfaces à partir d'esquisses
Cliquer sur l'une des parois et le modèle change de couleur
Générer le modèle et on aura une surface esquisse, une pièce et un corps
Enregister et fermer la fenêtre Géométrie- DesignModeler
1•2/Maillage:
Aller dans la fenêtre principale de workbench
Aller dans les systèmes de composants et faire translater Maillage vers le schéma du projet
Faire une translation entre géométrie de A et géométrie de B comme ça la géométrie qu'on a faite on la retrouve
dans notre fenêtre de maillage
On double clique sur Maillage de B
La fenêtre du maillage s'affiche B Maillage - meshing et on clique sur l'icône des arêtes
On définera les types de parois en sélectionnant chaque paroi et cliquant sur le bouton droit de la souris et créer une
sélection nommée
La paroi droite où l'eau entre on la nomme INLET
La paroi gauche on la nomme OUTLET
Les deux parois parallèle horizontales, on les nomme walls mais pour celà on crée une et puis on va dans sélections
nommées et on clique sur walls et on ajoute pour l'autre paroi ce même nom
On va dans modèle (B3)--» Maillage---» insertion---»Inflation
On régle les détails du maillage en cliquant sur maillage tout en définissant physique de préférence CFD et taille des
éléments 0,15 m et le format de l'export en Standard
On clique sur inflation et on clique sur notre modèle et on applique lavlimite sur ce modèle
Enfin on clique sur générer et on aura nitre modèle et notre maillage puis on enregistre notre maillage dans un
dossier
1•3/Fluent:
On va dans la fenêtre principale de workbench et on translate fluent dans le schéma de
projet
On translate Maillag de B dans configuration de C
On clique sur solution et on coche Double precision et do not show this panel again et on
choisit le 2 pour Solver processes et on double clique sur solution et une fenêtre s'affiche C:
Fluent Parallel Fluent
Aller dans Models---» viscous--» inviscid
Double cliquer sur Materials
Cliquer sur le bouton droit sur Fluid--» New--» Fluent database--» Fluent fluid Materials--» water-liquid(h2O «l»)--»
Density 998.2--» Copy--» Close
Supprimer air par le bouton droit
Aller dans cell Zone Conditions--» corps_surfacique--» Material Name--» Water liquid--» Apply
Aller dans cell zone boundary conditions--» inlet--» velocity magnitude 1 m/s--» Apply
Solution:
Methods--» scheme--» simple--» Apply
Metjods---» Residual--» changer les valeurs de continuity, velocity(x, y)--» fixer une valeur pour tous ces paramètres
0.000001--» Ok
Initializations--» Standard Initialization--» compute from Inlet--» initialize
Calculation acctivities--» Run Calculation--» Number if iterations--» 200--» calculate
Results--» bouton droit--» New--» Line/Rake...--» Name--» line x=0---» y0(0 m), y1(6 m)--» Save--» Close
Results--» graphics--» contours--» pressure--» static pressure; pour toutes les surfaces qu'a dessiné (on remarque
que la pression statique est grande sur le cercle dans le sens de l'écoulement et aux deux extrémités du haut et du
bas du cylindre la pression est tres petite sur la paroi du cylindre et en s'éloignant du cylindre la pression augmente)
Pour contours of velocity--» velocity magnitude(on constate que la vitesse est presque nulle sur la paroi du cylindre
et aux deux extrémités du haut et du bas la vitesse est très élèvée)
Pour Contours of Density(on constate que la densité est constante sur tout le modèle)=densité de l'eau)
Results--» graphics--» mesh--» on sélectionne les surface qu'on a--» save(on aura nitre maillage)
Cette fois on ne sélectionne pas le corps surfacique et son intérieur et on obtient ça:
Results--» graphics--» vectors--» velocity--» velocity magnitude--» pour toutes les surfaces sauf le corps surfacique
(la vitesse est représentée par des vecteurs)
Ici on ne sélectionne pas le corps surfacique et l'intérieur du corps
Graphics--» vectors--» color by pressure (la vitesse avec les couleurs de la pression)
Graphics--» pathlines--»--» cocher node values et write to file--» pressure--» static pressure--» pour toutes les
surfaces--» colore by pressure
Graphics--» pathlines--» colors by velocity magnitude--» pour toutes les surfaces
Graphics--» pathlines--» colors by density--» density
Results--» plots--» XY Plot--» node values--» write to file--» YAxis funcyion--» mesh--» Y-COordinate--» XAxis
Function--» Velocity---» X Velocity--» line x=0--» write--» on donne un nom au fichier et on l'enregistre
Ouvrir le fichier avec le bloc note et on aura ces valeurs suivantes qu'on va copier sur excel
Calcul de la solution exacte avec la formule 1+1/y^2
Figure 01: solution exacte et calculée pour un liquide non visqueux
On constate qu'à partir de cette courbe qu'elles sont presque les mêmes ce qui nous fait qire que notre modélisation
est bien réussit et que notre modèle est fiable
On calcule aussi le débit par la méthode des moindres carrés (1/2 *(Y(i+1)-Y(i))*(X(i+1)+X(i))) et on trouve qu'il est
presque le même de la valeur exacte qui est 6 m^3/s, il y a une très petite différence mais on peut constater que
notre modèle est fiable
Results--» Force--» direction vector X=1--» wall zones--» circle--» print
Results--» Surface integrals--»Report type---» Area--» circle--» print--» 6,266654 m^2
Results--» Surface integrals--»Report type---»integral--» Density--» circle--» compute---» 6266.654( kg/m^3)(m^2
Results--» Surface integrals--»Report type---»integral--» pressure--» static pressure--» circle---» compute--»-
2565.732(Pa)(m^2)
Results--» Surface integrals--»Report type---»integral--» velocity---» velocity magnitude---»compute--»
5.59002(m/s)(m^2)
Results--» Surface integrals--»Report type---»facet Average---» velocity--» velocity magnitude--» circle---» comoute--
»0.890543(m/s)
La partie de la viscosité laminar
On retourne à Setup et on change la viscosité :
Setup--->>viscous Laminar
Après on clique sur Boundary Conditions walls
On modifie Shear condition on valide Specified Shear apply
Et on change rien
Après on clique sur Initialization standard Initialization compute from « inlet» Initialization
Run calculation Calculation
Et on fait les mêmes étapes pour les conteurs et les vecteur
- Signiefie le sense inverse de l’ecoulement
Alors dans Results on trouve Forces
Double cliquer sur Forces
Cliquer sur circle write close
Et on change area to integral “pression”et “ velocity ” et “donsity”
Et on change integral to facet averge
Coclusion on conclut que notre simulation est fiable .
