Numerische Lösung partieller Differentialgleichungen

Dieser Bereich entwickelt Methoden, die partielle Differentialgleichungen in Raum und Zeit diskretisieren, indem sie kontinuierliche Operatoren durch algebraische Systeme ersetzen, deren Lösungen das Verhalten von Feldern annähern, die physikalischen Gesetzen unterliegen.

Definition

Die numerische Lösung partieller Differentialgleichungen ist die Konstruktion und Analyse von Methoden, die die Lösungen von partiellen Differentialgleichungen durch Diskretisierung des räumlichen Bereichs (und der Zeit) annähern, wodurch endliche Systeme algebraischer Gleichungen entstehen.

Scope

Er umfasst die drei wichtigsten Diskretisierungsrahmen – Finite-Differenzen-, Finite-Elemente- und Finite-Volumen-Methoden –, angewendet auf elliptische, parabolische und hyperbolische Gleichungen; die Analyse von Konsistenz, Stabilität und Konvergenz (einschließlich des Lax-Äquivalenztheorems und der CFL-Bedingung); und die großen dünnbesetzten linearen und nichtlinearen Systeme, die durch die Diskretisierung entstehen.

Sub-topics

• Finite-Differenzen-Methoden
• Finite-Elemente-Methoden
• Finite-Volumen-Methoden

Core questions

• Wie werden Differentialoperatoren in Raum und Zeit in stabile, konvergente algebraische Systeme diskretisiert?
• Wie kombinieren sich Konsistenz und Stabilität, um Konvergenz zu gewährleisten, wie im Lax-Äquivalenztheorem?
• Wie bestimmt der Typ der partiellen Differentialgleichung – elliptisch, parabolisch oder hyperbolisch – die geeignete Methode und Stabilitätsbedingungen?
• Wie werden die resultierenden großen dünnbesetzten Systeme effizient gelöst?

Key theories

Lax-Äquivalenztheorem

Für eine konsistente Finite-Differenzen-Approximation eines wohlgestellten linearen Anfangswertproblems ist Stabilität notwendig und hinreichend für Konvergenz; dieses Theorem ist der Eckpfeiler, der den Konvergenzbeweis auf die Überprüfung von Konsistenz und Stabilität reduziert.

Stabilitätsbedingungen und die CFL-Zahl

Explizite Schemata für zeitabhängige partielle Differentialgleichungen sind nur unter Einschränkungen der Schrittweiten stabil; für hyperbolische Probleme erfordert die Courant-Friedrichs-Lewy-Bedingung, dass der numerische Abhängigkeitsbereich den physikalischen enthält, was die Zeitschrittweite relativ zum räumlichen Gitter begrenzt.

Variations- und Erhaltungsprinzipien

Finite-Elemente-Methoden basieren auf schwachen (variationalen) Formulierungen und der Galerkin-Projektion, während Finite-Volumen-Methoden diskrete Erhaltungsgesetze durchsetzen; jeder Rahmen bietet einen Weg zu konsistenten Diskretisierungen mit nachweisbaren Approximationseigenschaften.

Clinical relevance

Numerische PDE-Methoden sind die rechnerische Grundlage der Simulation in Ingenieurwissenschaften und den physikalischen Wissenschaften – Struktur- und Festkörpermechanik, Fluiddynamik und Aerodynamik, Wärmeübertragung, Elektromagnetismus, Geophysik, Wetter- und Klimamodellierung sowie medizinische Bildrekonstruktion – überall dort, wo kontinuierliche Feldgleichungen auf komplexen Geometrien gelöst werden müssen, die geschlossene Lösungen ausschließen.

History

Die Finite-Differenzen-Analyse von partiellen Differentialgleichungen begann mit der Courant-Friedrichs-Lewy-Arbeit von 1928; die Finite-Elemente-Methode entstand aus dem Bauingenieurwesen und der Variationsmathematik in den 1940er-60er Jahren, und Finite-Volumen-Methoden entwickelten sich parallel zur numerischen Strömungsmechanik, wobei das Lax-Äquivalenztheorem in den 1950er Jahren den vereinheitlichenden Konvergenzrahmen lieferte.

Key figures

• Richard Courant
• Peter Lax
• Olga Ladyzhenskaya
• Randall J. LeVeque

Seminal works

• morton2005
• leveque2007

Frequently asked questions

Warum gibt es drei verschiedene Diskretisierungsrahmen?

Finite Differenzen sind am einfachsten auf regulären Gittern, Finite Elemente behandeln komplexe Geometrien und Variationsprobleme auf natürliche Weise, und Finite Volumen erzwingen lokale Erhaltung, was sie ideal für Strömungen macht. Die Wahl hängt von der Geometrie, dem Gleichungstyp und den Eigenschaften ab, die erhalten bleiben müssen.

Was bedeutet die CFL-Bedingung?

Für explizite Schemata bei zeitabhängigen hyperbolischen Problemen begrenzt die Courant-Friedrichs-Lewy-Bedingung, wie groß der Zeitschritt relativ zum räumlichen Gitterabstand sein darf, um sicherzustellen, dass Informationen nicht mehr als eine Gitterzelle pro Schritt zurücklegen. Eine Verletzung führt zu Instabilität.