É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 :

« Stiff equations are equations where certain implicit methods, in particular BDF (en), perform better, usually tremendously better, than explicit ones. »

Une formulation plus mathématique passe par le comportement des valeurs propres liés au système :

Soit l'équation différentielle

${\frac {\mathrm {d} u}{\mathrm {d} t}}=F(u).$

L'équation est dite raide si elle vérifie :

$\exists (\lambda _{1},\lambda _{2})\in \operatorname {Sp} (\mathrm {D} F(u)),\,|\lambda _{1}|<<|\lambda _{2}|$

où $\operatorname {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

${\begin{bmatrix}x'(t)\\y'(t)\end{bmatrix}}={\begin{bmatrix}-100&0\\0&-1\end{bmatrix}}\cdot {\begin{bmatrix}x(t)\\y(t)\end{bmatrix}},$

dont les solutions sont de la forme :

${\begin{bmatrix}x(t)\\y(t)\end{bmatrix}}={\begin{bmatrix}x(0)\,e^{-100\,t}\\y(0)\,e^{-t}\end{bmatrix}},$

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'(t)=-k\,[y(t)-\phi (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\phi$ 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).