La mécanique des fluides numérique (MFN), plus souvent désignée par le terme an
La mécanique des fluides numérique (MFN), plus souvent désignée par le terme anglais computational fluid dynamics (CFD), consiste à étudier les mouvements d'un fluide, ou leurs effets, par la résolution numérique des équations régissant le fluide. En fonction des approximations choisies, qui sont en général le résultat d'un compromis en termes de besoins de représentation physique par rapport aux ressources de calcul ou de modélisation disponibles, les équations résolues peuvent être les équations d'Euler, les équations de Navier-Stokes, etc. La MFN a grandi d'une curiosité mathématique pour devenir un outil essentiel dans pratiquement toutes les branches de la dynamique des fluides, de la propulsion aérospatiale aux prédictions météorologiques en passant par le dessin des coques de bateaux. Dans le domaine de la recherche, cette approche est l'objet d'un effort important, car elle permet l'accès à toutes les informations instantanées (vitesse, pression, concentration) pour chaque point du domaine de calcul, pour un coût global généralement modique par rapport aux expériences correspondantes. Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 1 sur 11 10/12/2019 à 16:40 Historique Méthodologie Méthodes de discrétisation Méthode des différences finies Méthode des volumes finis Méthode des éléments finis Résolution du problème discrétisé Forme générale des équations Discrétisation Maillage, symétries et périodicités Discrétisation spatiale et temporelle des équations Cas particuliers Diffusion numérique Propriétés des matériaux Conditions aux limites Modèles de turbulences Écoulement multiphasique Linéarisation Résidus Domaines d'application Méthodes de simulation de fluide Principaux logiciels Logiciels libres Logiciels commerciaux Références Voir aussi De manière générale, la résolution d'un problème de MFN passe par trois grandes phases : 1. la préparation du problème : ceci passe par la définition d'une géométrie, d'un maillage discrétisant le domaine de calcul, du choix des modèles et méthodes numériques employés ; 2. la résolution numérique du problème qui passe par l'exécution d'un programme informatique. Bien des problèmes suscitant un minimum d'intérêt nécessitent des ordinateurs aux très grandes capacités ; Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 2 sur 11 10/12/2019 à 16:40 3. l'exploitation des résultats : on vérifie d'abord leur cohérence, puis ils sont examinés afin d'apporter des réponses aux questions posées par le problème de MFN de départ. L'exploitation des résultats passe le plus souvent par des logiciels de post-traitement scientifique utilisés dans de nombreuses branches de la physique, ou bien par les modules de post-traitement disponibles dans certains logiciels de MFN commerciaux. La méthode des différences finies a une importance historique et est simple à programmer. Elle n'est actuellement utilisée que dans quelques codes spécialisés. La méthode des volumes finis est une approche commune utilisée dans les codes MFN. Les équations qui régissent le fluide sont résolues sur des volumes de contrôle discret. La méthode des éléments finis (MEF) est utilisée dans l'analyse structurale des solides, mais est également applicable aux fluides. Cependant, la formulation à éléments finis nécessite des soins particuliers pour assurer une solution conservative. La méthode consiste à résoudre des équations aux dérivées partielles (ÉDP) dites « équations de transport » ou « de conservation », et dont la forme générale est, pour une quantité scalaire φ donnée : où φ peut désigner : 1 (le nombre un) dans le cas de l'équation de continuité, Maillage d'une automobile. Méthodes de discrétisation Méthode des différences finies Méthode des volumes finis Méthode des éléments finis Forme générale des équations Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 3 sur 11 10/12/2019 à 16:40 les composantes du vecteur vitesse u, v et w, l'enthalpie h ; est le volume concerné et est la surface délimitant ce volume ; dV est le volume d'un élément infinitésimal de et le vecteur aire d'un élément infinitésimal de ; est l'opérateur nabla ; ρ est la masse volumique ; est le vecteur vitesse eulérienne d'une particule fluide ; Dφ est un coefficient de diffusion ; Sφ est un terme source, ou puits. Le maillage consiste à découper l'espace en cellules appelées « volumes de contrôle ». Le maillage est souvent plus délicat en mécanique des fluides que pour les éléments finis en résistance des matériaux : il faut mailler tout l'espace « vide » (la veine fluide), et les détails de surface ont de l'importance (puisqu'il génèrent des turbulences), on a donc souvent des maillages contenant de nombreuses mailles (typiquement plusieurs millions). Par ailleurs, alors qu'en résistance des matériaux, un maillage hexaédrique de l'espace est une stratégie intéressante (elle permet d'avoir la même qualité de résultats pour un nombre de nœuds plus faibles), en MFN, cela introduit des directions préférentielles qui peuvent avoir une influence sur le résultat (voir ci-après Diffusion numérique). Dans certains domaines, les pratiques professionnelles imposent parfois un maillage hexaédrique. Sinon, la stratégie retenue est souvent : à proximité des parois, avoir un maillage parallèle à la paroi — hexaédrique ou bien prismatique (hexaèdres coupés en deux) — et de faible épaisseur (trois couches de cellules ayant typiquement avec un facteur y+ de l'ordre de 1, voit ci-après Modèles de turbulences) ; au cœur du volume, avoir un maillage tétraédrique mais n'ayant pas de direction particulière (ce ne sont pas des cubes empilés coupés en deux). Une manière efficace de réduire la taille du modèle est d'utiliser les symétries et périodicités. Il faut pour cela que le système présente une symétrie ou une périodicité géométrique, et que les conditions limites présentent la même symétrie ou périodicité. Dans ce cas-là, les champs sortant sur une surface sont appliqués en entrée de la surface correspondante. Toutefois, il est délicat d'évaluer la pertinence de ces hypothèses. Par exemple, une condition de symétrie mal appliquée peut empêcher de voir un effet Coandă. La discrétisation spatiale consiste à remplacer les intégrales par des sommes sur des éléments de volume et de surface correspondant au maillage. Ainsi, pour chaque volume de contrôle (cellule), Discrétisation Maillage, symétries et périodicités Discrétisation spatiale et temporelle des équations Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 4 sur 11 10/12/2019 à 16:40 on peut écrire l'équation, en appliquant le théorème de flux-divergence : où Δ est l'opérateur laplacien ; ou bien, avec la notation nabla : La discrétisation temporelle consiste à faire des calculs à des instants déterminés, le résultat de la simulation à un instant ti étant la donnée d'entrée du calcul à l'instant ti + 1. Le pas de temps (ti + 1 − ti) peut être constant ou variable. La discrétion temporelle utilise typiquement la méthode des différences finies. Équation de continuité En prenant φ = 1, on trouve Équation de bilan de la quantité de mouvement En prenant φ = u (composante du vecteur vitesse sur l'axe x de vecteur directeur unitaire ), on trouve où p est la pression ; [τ] est le tenseur des contraintes dû à la viscosité ; est la flottabilité (poussée d'Archimède). Équation de la chaleur En prenant φ = htot (enthalpie totale), Cas particuliers Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 5 sur 11 10/12/2019 à 16:40 on trouve où p est la pression ; λ est la conductivité thermique ; T est la température absolue ; est la conduction thermique ; est le travail des forces visqueuses ; est la source d'enthalpie, typiquement la variation d'enthalpie due à un changement d'état (chaleur latente) ou à une réaction chimique (réaction exothermique ou endothermique). État stationnaire ou transitoire Un état stationnaire est une situation pour laquelle le terme transitoire est nul ; on néglige donc le premier terme des équations ci-dessus . Il s'agit d'un régime établi, d'un état de fonctionnement « longtemps » après le démarrage du système (ouverture du robinet, démarrage de la turbine) ; le résultat est indépendant du temps, il n'y a pas de discrétisation temporelle. Cette situation simplifie grandement les calculs. Cependant, cela ne permet pas de prendre en compte des fluctuations des conditions d'entrée ou de sortie, ou bien la description d'un remplissage. Dans ces cas-là, il faut faire un calcul transitoire (« instationnaire »), c'est-à-dire à plusieurs instants successifs. Si le pas de temps est trop grand, on introduit des erreurs numériques qui se propagent. Par ailleurs, cela peut rendre compliqué la convergence à chaque étape (voir ci-dessous Conditions aux limites). Mais plus le pas de temps est petit, plus le calcul est long et gourmand en ressources. Pour trouver un compromis, on s'attache en général à ce qu'une particule de fluide parcoure moins d'une maille entre chaque pas de calcul, c'est-à-dire à avoir un nombre de Courant, appelé dans le milieu CFL (condition de Courant, Friedrichs et Lewy), en général compris entre 0.1 et 0.6 (parfois plus dans le cas d'une résolution implicite, parfois moins selon l'application) : é é é é . Mécanique des fluides numérique — Wikipédia https://fr.wikipedia.org/wiki/Mécanique_des_fluides_numérique 6 sur 11 10/12/2019 à 16:40 La diffusion numérique est un artefact de calcul lié au maillage. Supposons qu'il y ait deux phases distinctes ; la frontière des phases passe au milieu de volumes de contrôle, la valeur des variables dans les-dits volumes est donc une « moyenne » des valeurs des deux phases. Il se crée donc un gradient artificiel local, puisque les valeurs sont potentiellement très uploads/Management/ mecanique-des-fluides-numerique.pdf
Documents similaires










-
51
-
0
-
0
Licence et utilisation
Gratuit pour un usage personnel Attribution requise- Détails
- Publié le Fev 03, 2021
- Catégorie Management
- Langue French
- Taille du fichier 0.4248MB