Équation différentielle raide

Une équation différentielle raide est une équation différentielle dont la sensibilité aux paramètres va rendre difficile la résolution par des méthodes numériques explicites. Plusieurs explications, aussi bien physiques que mathématiques, peuvent permettre d'appréhender la notion de raideur, qui reste difficilement formulable.

Définition

Il existe plusieurs définitions formelles de la raideur d'une équation différentielle. Une des plus simples est celle de Curtiss et Hirschfelder :

Modèle:Citation étrangère

Une formulation plus mathématique passe par le comportement des valeurs propres liés au système : Modèle:Énoncé où $\mathrm{Sp(M)}$ est le spectre de $M$.

La raideur vient donc d'une différence d'échelle entre deux phénomènes dynamiques : des phénomènes à évolution « rapide » vont tendre « lentement » vers une position d'équilibre.

Exemples

Un premier exemple simple est le système différentiel

$\left[\begin{array}{c}{x}^{\prime }(t)\\ y^{\prime }(t)\end{array}\right]=\left[\begin{array}{cc}-100& 0\\ 0& -1\end{array}\right]\cdot \left[\begin{array}{c}x(t)\\ y(t)\end{array}\right],$

dont les solutions sont de la forme :

$\left[\begin{array}{c}x(t)\\ y(t)\end{array}\right]=\left[\begin{array}{c}x(0){\mathrm{e}}^{-100t}\\ y(0){\mathrm{e}}^{-t}\end{array}\right],$

et la décroissance rapide de $x$ va imposer une forte contrainte sur la résolution de $y$.

On peut également citer les problèmes de Curtiss-Hirschfelder :

$y^{\prime }(t)=-k[y(t)-\varphi (t)],$

où $k$ est une constante positive réelle. Ce problème comporte lui aussi deux échelles de temps : la prédominance de $-ky$ pour les temps courts et de $+k\varphi$ sur les temps longs.

Les bilans de réactions chimiques peuvent également dans certains cas être décrites par des équations raides, comme la réaction de Belooussov-Jabotinski (réaction oscillante).

A-stabilité

Le comportement des méthodes numériques de résolution pour des problèmes raides peut être analysé en appliquant ces méthodes à l'équation différentielle test $y^{\prime }=ky$ avec pour condition initiale $y(0)=1$, et $k\in ℂ$. La solution de cette équation est $y(t)={\mathrm{e}}^{kt}$, et elle tend vers 0 pour $t\to \mathrm{\infty }$ dans le cas où $\mathrm{Re}(k)<0.$ Si la méthode numérique permet de retrouver cette propriété (pour un pas fixe), la méthode est dite A-stable.