ENSAM-Meknès
Modélisation des systèmes
énergétiques
Pr. Chadia HAIDAR
Objectifs et compétences à acquérir
Connaitre les méthodes de résolution numérique : Différence finies,
volumes finies ….
Connaitre les techniques de simulation numérique pour des problèmes de
mécanique des fluides et transfert thermique,
Acquérir les compétences suivantes :
• Savoir résoudre des problèmes linéaires et non linéaires
• Etre capable de réaliser une simulation numérique par un code de calcul
Introduction
CFD (Computational Fluid Dynamics) : En français (dynamique des fluides numérique), c'est la
simulation numérique en mécanique des fluides.
• Ensemble de méthodes numériques permettant d'étudier les mouvements d'un fluide, ou leurs
effets.
CFD permet d'avoir des solutions très satisfaisantes
grâce au développement des méthodes numériques et à
des calculateurs de plus en plus puissants avec une
grande capacité de mémoire.
1. Applications de la CFD
Cette technique s'applique actuellement dans plusieurs applications :
• Aérodynamique
• Hydrodynamique,
• Energies renouvelables;
• Combustion,
• Génie biomédical
• Turbomachines,
• Météorologie,
• Sport
• Etc.
Champ de vitesse autour d’une maquette de voiture.
2. Logiciels CFD aujourd’hui
1. Intérêts de la CFD
Pour traiter un problème de mécanique des fluides, on peut le résoudre soit
par une approche :
• Analytique :
Le grand avantage de cette approche est qu'elle permet d'avoir des solutions
exactes. Cependant, elle est limitée pour des cas très simples ;
• Expérimentale : représente le mieux la réalité, souvent difficile à mettre en œuvre et
nécessite un certain temps pour résoudre tous les problèmes qu'on peut rencontrer, elle
peut devenir très coûteuse ;
1. Intérêts de la CFD
• Numérique : permet de faire moins d'hypothèses qu'une approche analytique, traite des
problèmes relativement complexes, moins coûteuse qu'une approche expérimentale. Par contre,
elle est limitée par l'ordre des méthodes numériques utilisées, la précision du modèle et les
moyens de calcul mis en œuvre
1. Intérêts de la CFD
La CFD, en pratique, est utilisée pour:
• Études de conception et optimisation de (nouveaux) produits;
• Études des performances de systèmes/installations;
• Fournir des données complémentaires aux données théoriques et expérimentales.
1. Intérêts de la CFD
• Les méthodes de discrétisation les plus connues sont :
• Méthodes des différences finies (MDF);
• Méthodes des éléments finis (MEF);
• Méthodes des volumes finis (MVF).
▪ Méthode des différences finies
• Historiquement, la plus ancienne des trois. (a Principe: Basée sur l’approximation
été utilisée par Richardson en 1910). des dérivées: Les dérivées première et
seconde sont approchées par des
Avantage: facile à mettre en œuvre. développements en série de Taylor. Le
Inconvénients: limité aux maillages simples, système algébrique résultant est résolu
et ne conserve pas la quantité de mouvement, numériquement.
l’énergie et la masse sur des maillages
grossiers.
1. Intérêts de la CFD
• Les méthodes de discrétisation les plus connues sont :
• Méthodes des différences finies (MDF);
• Méthodes des éléments finis (MEF);
• Méthodes des volumes finis (MVF).
Méthode des éléments finis
• La première utilisation: Courant (1943) pour résoudre un
problème de torsion.
• Clough (1960): a donné à la méthode son nom. Principe: Basée sur
l’approximation des solutions:
Avantage: précision maximale sur les maillages grossiers. Chercher des solutions
approchées sur des domaines
Excellente pour les problèmes dominés par la diffusion: e.g.,
décrets.
écoulements visqueux à surface libre.
Inconvénients: lent pour les problèmes de grandes
dimensions et peu adapté aux écoulements turbulents.
1. Intérêts de la CFD
• Les méthodes de discrétisation les plus connues sont :
• Méthodes des différences finies (MDF);
• Méthodes des éléments finis (MEF);
• Méthodes des volumes finis (MVF).
▪ Méthode des volumes finis
• La première utilisation bien documentée a été faite par
Evans et Harlow (1957)
• Patankar et Spalding (1972-1980): applications pour
les problèmes des écoulements.
Principe: Intégration des équations
Avantage: masse, quantité de mouvement, énergie sur des VCs. Intégration la forme
conservée même sur des maillages grossiers. Utilisation conservative des équations sur chaque
pratique pour des maillages non structurés; volume de contrôle pour obtenir les
Inconvénients: mal-estimation de la diffusion lorsque équations discrétisées qui conservent
des schémas simples sont utilisés. toutes les grandeurs physiques.
2. Étapes de la CFD
2. Étapes de la CFD
• 2.1 CAO
• Toute simulation CFD commence par la réalisation de la
géométrie en 2D ou 3D, soit avec un logiciel intégré au code
de calcul CFD soit à l'aide d'un logiciel de CAO (CATIA,
SolidWorks,...).
• Dans le deuxième cas, la géométrie doit être exportée en un
format lisible.
Les formats STEP (STandardized Exchange of Product) ou
IGES (Initial Graphics Exchange Specication) sont les plus
couramment utilisés
2. Étapes de la CFD
2.2 Maillage
• Le maillage est une partition de l'espace ou d'un domaine physique en un ensemble
d'entités élémentaires (cellules) ;
Il peut être constitué d'éléments linéiques, surfaciques et volumiques ;
2. Étapes de la CFD
Composants du maillage
• Le domaine de calcul est défini par un maillage qui représente le fluide et les faces
solides qui interviennent.
• "Cell" : volume de contrôle divisant la géométrie ;
• "Face" : frontière d'une "cell", où sont définies les conditions aux limites ;
• "Edge" : frontière d'une face ;
• "Node" : point de maillage ;
• "Zone" : groupe de "noeuds", "faces" et/ou "cells".
2. Étapes de la CFD
2.2 Maillage
La connectivité d'un maillage est la définition des liaisons entre les sommets de ses
éléments ;
La qualité du maillage a un sérieux impact sur la convergence, la précision de la solution et
le temps de calcul ;
Une bonne qualité de maillage repose sur la minimisation des éléments présentant des
"distorsions", sur une bonne "résolution" dans les régions présentant un fort gradient
(couches limites, ondes de choc,...)et un bon maillage doit également être suffisamment
lisse ;
2. Étapes de la CFD
• Principaux types d'éléments utilisés en 2D
2. Étapes de la CFD
• Principaux types d'éléments utilisés en 3D
2. Étapes de la CFD
2. Étapes de la CFD
Type de Maillage
Maillage non structuré (tri/tétra) : Les éléments de ce type de maillage sont générés
arbitrairement sans aucune contrainte quant-à leur disposition.
2. Étapes de la CFD
• Avantages :
• Peut être généré sur une géométrie complexe tout en gardant une bonne qualité des
éléments,
• Les algorithmes de génération de ce type de maillage (tri/tétra) sont très automatisés.
• Inconvénients :
• Très gourmand en nombre de mailles comparativement au maillage structuré.
• Impose une structure de données gourmande en capacités de stockage.
• Engendre des erreurs numériques qui peuvent être plus importantes si on le compare avec
le maillage structuré.
• Il est difficile de contrôler localement la densité des mailles.
2. Étapes de la CFD
• Maillage structuré (quadr/hexa) : Un maillage structuré est un maillage qui peut
être généré en reproduisant plusieurs fois une maille élémentaire.
2. Étapes de la CFD
Avantages :
• Économique en nombre d'éléments, présente un nombre inférieur de mailles par rapport à un maillage non
structuré équivalent.
• Lorsque l'écoulement moyen est aligné avec le maillage, un maillage structuré réduit les risques d'erreurs
numériques.
• L'utilisateur a le contrôle total de la qualité du maillage, contrairement aux maillages non structurés plus
dépendants de l'algorithme de maillage.
• Les calculs sont généralement plus rapides dans un maillage structuré que dans un maillage non structuré.
Inconvénients :
Difficile à générer dans le cas d'une géométrie complexe
2. Étapes de la CFD
• Maillage hybride : Maillage généré par un mélange d'éléments de différents types,
triangulaires ou quadrilatéraux en 2D, tétraédriques, prismatiques, ou pyramidaux en
3D. Il combine les avantages des maillages structurés et non structurés.
2. Étapes de la CFD
Techniques générales de génération du maillage
En pratique, il n'existe pas de règle précise pour la création d'un maillage valable, cependant il
existe différentes approches qui permettent d'obtenir une grille acceptable :
• Maintenir une bonne qualité des éléments,
• Assurer une bonne résolution dans les régions à fort gradient,
• Assurer un bon lissage dans les zones de transition entre les parties à maillage fin et les
parties à maillage grossier,
• Minimiser le nombre total des éléments (temps de calcul raisonnable).
2. Étapes de la CFD
Distorsion
• Le facteur de distorsion Fd (skewness) peut être
défini de deux façons différentes. Pour des éléments
triangulaires ou tétraédriques, le calcul est basé sur le
volume équilatéral :
L'autre formulation, applicable pour tout élément, est basée sur la déviation angulaire, avec θ les
angles de l'élément : (60 ou 90)
Table - Variation de la qualité des éléments de maillage en fonction de la valeur du Fd
2. Étapes de la CFD
• Raffinement de maillage
La notion de résolution concerne plus particulièrement les zones qui présentent un fort
gradient (gradient de pression, vitesse, température,...) ;
Une bonne résolution
un maillage localement plus fin- permet de mieux décrire les phénomènes physiques qui
existent dans ces zones ;
La plupart des mailleurs proposent des méthodes de maillage particulières pour réaliser
en proche paroi des maillages structurés dont on maîtrise la taille des éléments ;
2. Étapes de la CFD
• Lissage
Le changement dans la taille des éléments du maillage d'une zone maillée à une autre doit
être graduel, la variation de la taille des éléments de deux zones adjacentes ne doit
idéalement pas dépasser 20% à 30%.
2. Étapes de la CFD
• Nombre total d'éléments
Des techniques existent pour économiser un certain nombre d'éléments :
Utilisation des maillages non uniformes, en concentrant la bonne qualité du maillage
uniquement dans les zones où c'est nécessaire ;
Utilisation de la fonction adaptation de maillage pour raffiner uniquement sur des zones
bien précises ;
Utilisation des éléments de maillage hexaédriques dans les zones adéquates.
2. Étapes de la CFD
• Convergence en maillage
2. Étapes de la CFD
2.3 Simulation
Les principales étapes d'une simulation (Code FLUENT) sont dans l'ordre :
1. Paramétrage du solveur ;
2. Modélisation de la turbulence ;
3. Définition des caractéristiques du fluide
4. Paramétrage des conditions aux limites ;
5. Choix des critères de convergence ;
6. Initialisation du calcul ;
7. Paramétrage des sauvegardes automatiques en cours de simulation ;
8. Lancement de la simulation.
2. Étapes de la CFD
Choix des équations résolues par le solveur
• Dimension du problème, 2D, 3D, ou axisymétrique ;
• Caractère stationnaire ou instationnaire ;
• Caractère incompressible ou compressible de l'écoulement (et non du fluide) ;
• Besoin ou non de prendre en compte le bilan d'énergie ;
• Prise en compte de forces extérieures (gravité, force d'inertie liées à la rotation,...) ;
• Modélisation éventuelle de la turbulence ;
2. Étapes de la CFD
• Définition des caractéristiques du fluide
Les codes disposent de bibliothèques des fluides les plus courants.
• Conditions aux limites
On distingue des conditions d'entrée, de sortie, aux frontières solides, de périodicité ou de symétrie. Pour
le code FLUENT, on a le choix pour les conditions d'entrée entre :
Velocity Inlet : utilisée pour définir la vitesse et d'autres grandeurs scalaires de l'écoulement
d'entrée.
Pressure Inlet : utilisée pour définir la pression totale et d'autres grandeurs scalaires de
l'écoulement d'entrée.
2. Étapes de la CFD
• Mass Flow Inlet : utilisée dans les écoulements compressibles pour prescrire un débit massique à
l'entrée. Il n'est pas nécessaire d'utiliser des entrées de débit massique dans les écoulements
incompressibles car lorsque la densité est constante, les conditions aux limites d'entrée de vitesse
fixent le débit massique.
Inlet Vent : utilisée pour modéliser le vent d'entrée avec un coefficient des pertes de charge, une
direction d'écoulement, une température, et une pression totale.
Intake Fan : Il est utilisé pour modéliser un ventilateur d'admission externe avec une chute de
pression, une direction d'écoulement, une pression et une température d'admission spécifiés
2. Étapes de la CFD
pour les conditions de sortie :
Pressure Outlet : spécifie la pression statique de sortie et d’autres grandeurs scalaires.
Outflow : utilisée pour modéliser l’écoulement de sortie où la vitesse et la pression d’écoulement
ne sont pas connues. Cette condition n’est pas appropriée pour les calculs suivants :
- Si le problème possède une condition de type "Pressure Inlet" ;
- Si vous modélisez un écoulement compressible ;
- Si vous modélisez un écoulement instationnaire avec variation de la densité.
Pressure Far-field : uniquement en compressible, on modélise un écoulement libre de nombre de
Mach connu.
Outlet Vent : utilisée pour modéliser le vent de sortie avec un coefficient des pertes de charge, une
température, et une pression statique.
Outake Fan : Il est utilisé pour modéliser un ventilateur d’extraction externe avec une chute de
pression spécifié et une pression statique de refoulement.
2. Étapes de la CFD
Exemple
2. Étapes de la CFD
Post-traitement
C'est une rubrique dédiée au traitement et à l‘ affichage des résultats à savoir :
Domaine de calcul + maillage ;
Résultats (champ de vitesse, pression, température, turbulence,...) ; Courbes,... ;
Extraction des résultats sous forme fichiers
Animation des résultats dynamiques ;
Les logiciels de CFD fournissent souvent des outils intégrés de post-traitement.
Toutefois, rien n'interdit d'exporter les données dans divers formats afin de les traiter
sous d'autres logiciels de visualisation.
