Méthodes des EDP en physique computationnelle

Les équations de champ de la physique, de la diffusion et des ondes à l'électrostatique, sont des équations aux dérivées partielles, et leur résolution numérique implique de discrétiser l'espace et le temps en une grille et d'y propager ou de relaxer le champ.

Definition

Les méthodes des EDP en physique computationnelle sont des schémas numériques qui approximent la solution d'équations aux dérivées partielles sur une grille discrète, en remplaçant les dérivées spatiales et temporelles par des différences finies ou des opérateurs apparentés.

Scope

Ce sujet couvre la discrétisation par différences finies des classes canoniques d'EDP, elliptiques, paraboliques et hyperboliques, ainsi que les méthodes de pas de temps explicites et implicites, de relaxation et multigrilles pour les problèmes aux limites, et les critères de stabilité qui les régissent. Les approches par éléments finis et spectrales sont considérées comme des méthodes connexes.

Core questions

• Comment les EDP elliptiques, paraboliques et hyperboliques sont-elles discrétisées et résolues différemment ?
• Qu'est-ce que la condition de Courant-Friedrichs-Lewy et pourquoi limite-t-elle les pas de temps explicites ?
• Comment les méthodes de relaxation et multigrilles résolvent-elles efficacement les problèmes aux limites de grande taille ?
• Quand un schéma implicite justifie-t-il son coût supplémentaire par rapport à un schéma explicite ?

Key theories

Discrétisation par différences finies

Les dérivées spatiales et temporelles sont remplacées par des quotients de différences sur une grille, convertissant une EDP en un grand système d'équations algébriques dont la précision est déterminée par l'espacement de la grille et l'ordre du stencil.

Condition de stabilité CFL

Pour les schémas explicites résolvant des équations hyperboliques et paraboliques, la condition de Courant-Friedrichs-Lewy borne le pas de temps par rapport à l'espacement de la grille et à la vitesse de propagation, au-delà de laquelle la solution numérique diverge.

Relaxation et multigrille

Les problèmes aux limites elliptiques, tels que l'équation de Poisson, sont résolus par relaxation itérative, les méthodes multigrilles accélérant la convergence en corrigeant les erreurs à travers une hiérarchie de résolutions de grille.

Clinical relevance

Les solveurs d'EDP calculent les champs électrostatiques et magnétostatiques, la conduction thermique et la diffusion, la propagation des ondes et l'équation de Schrödinger, constituant l'épine dorsale de l'électromagnétisme computationnel, de la dynamique des fluides et de la simulation en physique des milieux continus.

History

La théorie systématique des solutions par différences finies des EDP a débuté avec l'article de Courant-Friedrichs-Lewy de 1928 sur la stabilité, s'est considérablement développée avec les ordinateurs au milieu du XXe siècle, et a été rendue efficace pour les problèmes de grande taille par le développement des méthodes multigrilles dans les années 1970.

Key figures

• Richard Courant
• Kurt Friedrichs
• Randall J. LeVeque

Seminal works

• leveque2007
• press2007

Frequently asked questions

Quelle est la différence entre les pas de temps explicites et implicites ?

Les schémas explicites calculent le niveau de temps suivant directement à partir du niveau actuel et sont peu coûteux par pas, mais limités par une condition de stabilité sur la taille du pas. Les schémas implicites résolvent un système couplé à chaque pas, coûtant plus cher par pas mais restant stables pour des pas beaucoup plus grands, ce qui est avantageux pour les problèmes raides ou diffusifs.

Pourquoi classer les EDP en elliptiques, paraboliques ou hyperboliques ?

La classification reflète la manière dont l'information se propage : les équations elliptiques décrivent des champs d'équilibre avec un couplage global, les équations paraboliques décrivent une diffusion lissante dans le temps, et les équations hyperboliques décrivent des ondes se propageant à vitesse finie. Chaque classe requiert des stratégies de discrétisation et de stabilité différentes.