Consistance (mathématiques)

On considère une équation différentielle ordinaire (EDO) de la forme : ${\displaystyle y'(t)=f(t,y(t)),\ \forall t\in [0,T],}$ qu'on essaie de résoudre numériquement par un schéma de la forme ${\displaystyle y_{n+1}=y_{n}+\Delta t_{n}\Phi (t_{n},y_{n},\Delta t_{n}),\forall n\in [0,N].}$ avec t₀,..., t_N+1 une subdivision de [0 , T] et Δt_n=t_n+1-t_n.

L’erreur de consistance relative à une solution exacte z de l'équation est donnée par : ${\displaystyle e_{n}=z(t_{n+1})-y_{n+1}=z(t_{n+1})-z(t_{n})-\Delta t_{n}\Phi (t_{n},z(t_{n}),\Delta t_{n}),\forall n\in [0,N].}$

Considérons un problème supposé être bien posé qui modélise un système évolutif caractérisé par :

• une condition initiale précisant son état d’origine (variables spatiales en t = 0, soit U (x,0) = U₀(x)) ;
• des équations aux dérivées partielles du type

${\displaystyle \partial _{t}U(x,\,t)=D\,U(x,\,t)}$ où D est un opérateur différentiel relatif aux variables spatiales ;

• des conditions de bord auxquelles est soumis l’état du système au cours de son évolution.

Dans ce contexte, un schéma numérique procède de la manière suivante :

• Discrétisation des variables spatiales (pas Δx) pour établir une approximation numérique de l’état d’origine.
• Discrétisation de la variable temporelle (sur [0 , T] avec un pas Δt) pour entreprendre un processus se déroulant par étapes successives au cours desquelles l’état numérique se transforme.

Notons C(Δt) l’opérateur de modification de l’état discret au cours d’une étape, ceci en supposant une relation liant Δx à Δt qui contraint Δx à converger vers 0 lorsque Δt fait de même.

La consistance exige que, pour une fonction u(x) régulière, alors C(Δt)u – u soit une bonne approximation de Δt Du, ceci au sens de la ${\displaystyle \|.\|_{\infty }}$ définie par ${\displaystyle \|v\|_{\infty }={\text{sup}}\,v(x)}$.

La consistance est définie par la propriété suivante [1] :

Pour toute solution régulière u(x,t), alors

${\displaystyle \lim _{\Delta t\to 0}\left\|\left({\frac {C(\Delta t)\,-\,I}{\Delta t}}\,-\,D\right)\,u(.,\,t)\right\|_{\infty }=0}$

et la convergence est uniforme sur ${\displaystyle 0\leqslant t\leqslant T}$.

• Théorème de Lax
• Stabilité d'un schéma numérique