Méthodes numériques pour intégrer des équations différentielles ordinaires et des flux, conçues pour préserver des structures géométriques qualitatives du système continu (par exemple la forme symplectique, les intégrales premières, le volume ou la structure de groupe de Lie) plutôt que de ne viser que la minimisation de l'erreur locale de troncature.

Principe

Construire des schémas d'intégration qui respectent les objets géométriques invariants de la dynamique continue (formes symplectiques, quantités conservées, symétries de groupe, réversibilité) afin que la trajectoire numérique reflète le comportement qualitatif à long terme.

Démonstration

Un intégrateur symplectique comme le point milieu implicite ou la méthode Stoermer–Verlet appliqué à un système hamiltonien préserve la forme symplectique et présente une erreur d'énergie bornée sur des durées exponentiellement plus longues que les méthodes de Runge–Kutta classiques de même ordre.

Mauvaise application

Utiliser une méthode choisie uniquement pour son ordre local élevé (par ex. un Runge–Kutta explicite d'ordre élevé) sur un système quasi-hamiltonien en espérant une fidélité énergétique à long terme ; ou imposer un schéma symplectique à un problème fortement dissipatif où préserver une forme symplectique déforme la dissipation physique.

Conséquence

Lorsque les structures géométriques sont préservées, les solutions numériques conservent des invariants qualitatifs (dérive d'énergie bornée, topologie correcte de l'espace des phases, moments conservés), permettant des simulations de longue durée fiables et de meilleures prédictions qualitatives.

Inversion

Les intégrateurs traditionnels non-structurants privilégient l'erreur locale de troncature et peuvent présenter une dérive séculaire des invariants, un amortissement ou une croissance artificielle, et un comportement asymptotique qualitativement incorrect.

Limite

S'applique principalement aux flux d'EDO déterministes de dimension finie et aux champs vectoriels réguliers ; les extensions aux problèmes raides, aux EDP ou aux systèmes stochastiques exigent des techniques supplémentaires ou des cibles de préservation différentes et peuvent ne pas préserver tous les invariants continus.

Tension sémantique

Tension avec la recherche d'un ordre local élevé et d'une erreur par pas minimale : un intégrateur géométrique d'ordre inférieur peut surpasser un méthode non géométrique d'ordre supérieur en fidélité qualitative sur le long terme ; tension également avec les intégrateurs strictement conçus pour la stabilité des problèmes raides.

Synthèse

L'intégration géométrique regroupe la conception et l'analyse d'intégrateurs numériques qui sacrifient l'erreur locale pure pour préserver des invariants géométriques continus, garantissant un comportement qualitatif correct à long terme pour les systèmes dynamiques.