Méthodes des volumes finis

Les méthodes des volumes finis discrétisent les lois de conservation en divisant le domaine en volumes de contrôle et en mettant à jour la valeur moyenne dans chacun à partir des flux traversant ses frontières, conservant ainsi la masse, la quantité de mouvement et l'énergie par construction.

Definition

Une méthode des volumes finis est une discrétisation d'une loi de conservation qui stocke la moyenne de la solution sur chaque volume de contrôle et fait évoluer ces moyennes en équilibrant les flux numériques à travers les frontières du volume, de sorte que le schéma discret est localement et globalement conservatif.

Scope

Ce sujet couvre la forme intégrale (de conservation) des équations aux dérivées partielles (EDP), les inconnues moyennées par cellule et les fonctions de flux numériques, la méthode de Godunov et les solveurs de Riemann approchés, les schémas à haute résolution avec des limiteurs de pente qui suppriment les oscillations parasites près des chocs, et le rôle des méthodes des volumes finis en dynamique des fluides numérique.

Core questions

• Pourquoi le fait de travailler à partir de la forme intégrale de conservation rend-il la méthode intrinsèquement conservative ?
• Comment les flux numériques aux interfaces des cellules sont-ils définis, et qu'est-ce qui rend un flux cohérent et stable ?
• Comment les schémas de type Godunov et les solveurs de Riemann capturent-ils les discontinuités telles que les chocs ?
• Comment les méthodes à haute résolution évitent-elles les oscillations des schémas d'ordre élevé près des discontinuités ?

Key theories

Conservation et flux numérique

En mettant à jour les moyennes des cellules à l'aide d'un seul flux numérique partagé entre les cellules adjacentes, la méthode conserve exactement la quantité sous-jacente ; la cohérence du flux avec le flux réel et une condition de stabilité appropriée conduisent à la convergence vers des solutions faibles de la loi de conservation.

Méthodes de Godunov et solveurs de Riemann

L'approche de Godunov traite chaque interface de cellule comme un problème de Riemann local dont la solution (exacte ou approchée) définit le flux, permettant au schéma de capturer les chocs et les discontinuités de contact de manière nette et correcte.

Schémas à haute résolution et limiteurs

Pour surmonter la précision de premier ordre des schémas de Godunov de base sans introduire d'oscillations parasites, les méthodes à haute résolution reconstruisent des états d'interface d'ordre supérieur et appliquent des limiteurs de pente ou de flux qui imposent une propriété de diminution de la variation totale près des discontinuités.

Mechanisms

L'intégration de la loi de conservation sur un volume de contrôle convertit les dérivées spatiales en flux de surface, de sorte que le taux de variation d'une moyenne de cellule est égal au flux net à travers ses faces. Étant donné que les cellules adjacentes partagent chaque flux de face, tout ce qui quitte une cellule entre dans sa voisine et la quantité totale est conservée exactement. Le flux d'interface est calculé en résolvant, exactement ou approximativement, le problème de Riemann posé par les états différents de part et d'autre ; cela permet de capturer la structure des ondes et les discontinuités. Les variantes à haute résolution reconstruisent d'abord un profil d'ordre supérieur à l'intérieur de chaque cellule et le limitent pour éviter de nouveaux extrema, atteignant ainsi une précision de second ordre dans les régions lisses tout en restant sans oscillation au niveau des chocs.

Clinical relevance

Les méthodes des volumes finis constituent la discrétisation standard en dynamique des fluides numérique et sont essentielles pour simuler les écoulements compressibles et incompressibles, l'aérodynamique, les phénomènes de choc et de détonation, les écoulements en eau peu profonde, atmosphériques et océaniques, ainsi que la simulation des milieux poreux et des réservoirs, précisément parce qu'elles respectent les lois de conservation physiques et gèrent les discontinuités de manière robuste.

History

L'approche conservative basée sur le problème de Riemann a vu le jour avec le schéma de Godunov en 1959 ; le développement de méthodes et de limiteurs à haute résolution, diminuant la variation totale, par van Leer, Harten et d'autres dans les années 1970 et 1980, a fait des méthodes des volumes finis le cadre dominant pour les écoulements compressibles et d'autres lois de conservation hyperboliques.

Key figures

• Sergei Godunov
• Peter Lax
• Bram van Leer
• Randall J. LeVeque
• Eleuterio Toro

Seminal works

• leveque2002
• toro2009

Frequently asked questions

Pourquoi les méthodes des volumes finis sont-elles conservatives ?

Elles mettent à jour les moyennes des cellules en utilisant des flux partagés entre les cellules voisines, de sorte que toute quantité quittant une cellule entre exactement dans la cellule adjacente. En sommant sur l'ensemble du domaine, la quantité totale ne change qu'à travers la frontière du domaine, ce qui reflète la loi de conservation physique.

Pourquoi les limiteurs sont-ils nécessaires près des chocs ?

Les schémas d'ordre élevé simples produisent des oscillations parasites (dépassements et sous-dépassements) près des discontinuités. Les limiteurs de pente et de flux détectent les gradients raides et réduisent localement l'ordre de reconstruction, maintenant la solution monotone et sans oscillation tout en préservant la précision dans les régions lisses.