Méthodes aux différences finies

Les méthodes aux différences finies approximent les dérivées par des quotients de différences sur une grille, transformant ainsi une équation différentielle en un système d'équations algébriques pour les valeurs de la solution aux points de la grille.

Definition

Une méthode aux différences finies est une discrétisation d'une équation différentielle dans laquelle les dérivées sont remplacées par des quotients de différences de l'inconnue évalués sur une grille structurée, produisant des équations algébriques dont la solution approche la solution de l'équation différentielle aux points de la grille.

Scope

Ce sujet couvre la construction d'approximations aux différences à partir de développements de Taylor, la discrétisation des EDP elliptiques, paraboliques et hyperboliques, les schémas d'intégration temporelle explicites et implicites (tels que Euler explicite, Euler implicite et Crank-Nicolson), l'analyse de stabilité de von Neumann, et le cadre de cohérence-stabilité-convergence spécialisé pour les schémas aux différences.

Core questions

• Comment les approximations aux différences précises des dérivées sont-elles dérivées et leur erreur de troncature quantifiée ?
• En quoi les schémas d'intégration temporelle explicites et implicites diffèrent-ils en termes de stabilité et de coût ?
• Comment l'analyse de von Neumann détermine-t-elle la stabilité d'un schéma aux différences ?
• Comment le type d'équation détermine-t-il le schéma approprié et les éventuelles restrictions sur la taille du pas ?

Key theories

Cohérence, stabilité et convergence

Un schéma aux différences est cohérent si son erreur de troncature tend vers zéro lorsque la grille est raffinée et stable si les erreurs ne croissent pas de manière illimitée ; selon le théorème d'équivalence de Lax, ces deux propriétés garantissent la convergence vers la vraie solution pour les problèmes linéaires bien posés.

Analyse de stabilité de von Neumann

La décomposition de l'erreur en modes de Fourier sur une grille uniforme réduit la stabilité à la limitation d'un facteur d'amplification pour chaque mode ; le schéma est stable lorsqu'aucun mode n'est amplifié, ce qui donne des conditions explicites sur la taille du pas, telles que les limites de diffusion et CFL.

Mechanisms

Les formules aux différences sont construites en combinant des développements de Taylor aux points de grille voisins pour annuler les termes d'ordre inférieur et isoler une dérivée, le terme résiduel principal donnant l'erreur de troncature et l'ordre de la méthode. Pour les problèmes dépendants du temps, les schémas explicites mettent à jour chaque nouvelle valeur directement à partir des anciennes mais doivent respecter une limite de stabilité (une borne de nombre de diffusion pour les équations paraboliques, la condition CFL pour les équations hyperboliques), tandis que les schémas implicites comme Crank-Nicolson couplent les nouvelles valeurs dans un système linéaire qui est inconditionnellement stable mais nécessite une résolution à chaque étape. L'analyse de von Neumann substitue des modes de Fourier pour tester la stabilité et déterminer ces limites.

Clinical relevance

Les méthodes aux différences finies sont largement utilisées pour des problèmes sur des domaines réguliers et des grilles structurées : conduction thermique et diffusion, propagation des ondes et modélisation sismique, électromagnétisme computationnel (la méthode des différences finies dans le domaine temporel), et l'évaluation d'options via l'équation de Black-Scholes ; leur simplicité et la facilité d'extension à des ordres élevés en font un premier choix lorsque la géométrie est simple.

History

Les fondements mathématiques ont été posés par l'article de Courant-Friedrichs-Lewy de 1928 sur les équations aux différences pour les EDP ; l'analyse de stabilité de von Neumann pendant la guerre et le théorème d'équivalence de Lax des années 1950 ont établi la théorie moderne, et les méthodes aux différences restent un pilier de la physique computationnelle et de l'ingénierie.

Key figures

• Richard Courant
• Kurt Friedrichs
• Hans Lewy
• John von Neumann
• Randall J. LeVeque

Seminal works

• leveque2007
• morton2005

Frequently asked questions

Quelle est la différence entre un schéma explicite et un schéma implicite ?

Un schéma explicite calcule chaque nouvelle valeur de la grille directement à partir des valeurs connues mais n'est stable que pour de petits pas de temps, tandis qu'un schéma implicite résout un système couplé pour toutes les nouvelles valeurs simultanément, permettant des pas de temps stables beaucoup plus grands au prix d'une résolution linéaire par étape.

Pourquoi les différences finies pourraient-elles être préférées aux éléments finis ?

Sur des géométries simples et régulières, les différences finies sont faciles à implémenter, peu coûteuses et simples à étendre à des ordres élevés. Les éléments finis deviennent avantageux principalement lorsque le domaine a une forme complexe ou que le problème possède une formulation variationnelle naturelle.