Équations différentielles raides et stabilité

Les équations différentielles raides contiennent des processus évoluant sur des échelles de temps très différentes, de sorte que les méthodes explicites sont contraintes de prendre des pas de temps d'une petitesse impraticable pour assurer la stabilité ; leur résolution efficace nécessite des méthodes implicites dotées de fortes propriétés de stabilité.

Definition

Une équation différentielle est dite raide lorsqu'elle admet des composantes de solution qui décroissent sur des échelles de temps très différentes, de sorte que la stabilité numérique, plutôt que la précision, dicte la taille du pas ; la théorie de la stabilité analyse quelles méthodes peuvent prendre de grands pas sans croissance d'erreur.

Scope

Ce sujet aborde le phénomène et la définition informelle de la raideur, l'équation test linéaire et la région de stabilité absolue, les concepts de A-stabilité, A(alpha)-stabilité et L-stabilité, les raisons de l'échec des méthodes explicites sur les problèmes raides, ainsi que les méthodes implicites — Runge-Kutta implicites et formules de différenciation arrière (backward differentiation formulas) — qui permettent de les résoudre.

Core questions

• Qu'est-ce qui rend un problème raide, et pourquoi cela met-il en échec les méthodes explicites ?
• Comment la région de stabilité absolue est-elle définie à travers l'équation test linéaire ?
• Que requièrent la A-stabilité et la L-stabilité, et pourquoi sont-elles importantes pour les problèmes raides ?
• Quelles méthodes offrent la stabilité nécessaire pour les systèmes raides et différentiels-algébriques ?

Key theories

Stabilité absolue et l'équation test

L'application d'une méthode à l'équation test linéaire scalaire produit un facteur d'amplification ; l'ensemble des produits (taille du pas multipliée par la valeur propre) pour lesquels ce facteur a une magnitude au plus égale à un est la région de stabilité absolue de la méthode, laquelle doit contenir les valeurs propres raides du problème pour permettre de grands pas.

A-stabilité et L-stabilité

Une méthode est A-stable si sa région de stabilité contient tout le demi-plan gauche, elle est donc stable pour tous les modes décroissants quelle que soit la taille du pas, et L-stable si elle amortit de plus complètement les modes très raides ; ces propriétés distinguent les méthodes implicites adaptées aux problèmes raides.

Mechanisms

Dans un problème raide, le mode à décroissance la plus rapide possède une grande valeur propre négative ; la région de stabilité bornée d'une méthode explicite contraint la taille du pas à résoudre ce mode même longtemps après qu'il ait physiquement disparu, rendant le calcul désespérément lent. Les méthodes implicites telles que la méthode d'Euler arrière (backward Euler method), les schémas de Runge-Kutta implicites et les formules de différenciation arrière (backward differentiation formulas) ont des régions de stabilité couvrant le demi-plan gauche (ou la majeure partie de celui-ci), elles restent donc stables avec de grands pas et permettent de choisir la taille du pas uniquement en fonction de la précision. Chaque pas nécessite alors la résolution d'un système algébrique (généralement non linéaire), typiquement par une itération de Newton utilisant le Jacobien.

Clinical relevance

La raideur est omniprésente dans les réseaux de réactions chimiques, la combustion, les circuits électriques, les systèmes de contrôle et les discrétisations par la méthode des lignes (method-of-lines) des équations aux dérivées partielles paraboliques ; reconnaître la raideur et sélectionner un solveur implicite suffisamment stable est essentiel pour obtenir des résultats dans un délai raisonnable, et la plupart des logiciels de production d'EDO incluent une détection automatique de la raideur et une commutation de méthode.

History

La notion de raideur a été identifiée par Curtiss et Hirschfelder en 1952, et la théorie de la stabilité sous-jacente — A-stabilité et les barrières d'ordre — a été développée par Dahlquist ; les codes de formules de différenciation arrière (backward differentiation formula) de Gear et, plus tard, les méthodes de Runge-Kutta implicites d'ordre élevé ont établi la boîte à outils pratique pour les problèmes raides et différentiels-algébriques.

Key figures

• Germund Dahlquist
• C. William Gear
• Ernst Hairer
• Gerhard Wanner

Seminal works

• hairer1996
• iserles2008

Frequently asked questions

Qu'est-ce qui rend précisément une EDO raide ?

La raideur apparaît lorsque le système possède des composantes qui décroissent beaucoup plus rapidement que la solution d'intérêt n'évolue. Il n'existe pas de définition unique et rigoureuse, mais la caractéristique pratique est que les méthodes explicites sont contraintes d'utiliser de très petits pas pour la stabilité, même lorsque la précision permettrait d'en prendre de grands.

Pourquoi les problèmes raides nécessitent-ils des méthodes implicites ?

Les méthodes implicites peuvent avoir des régions de stabilité couvrant tout le demi-plan gauche (A-stabilité), elles restent donc stables avec de grands pas pour les modes à décroissance rapide. Les méthodes explicites ont des régions de stabilité bornées, ce qui contraint à de très petits pas et les rend impraticables pour les problèmes raides.