En analyse numérique, la consistance d’un schéma numérique aux différences finies est une propriété locale de l’algorithme qui en découle. Elle concerne essentiellement la capacité du schéma à représenter une solution régulière satisfaisant localement les équations aux dérivées partielles, ceci lorsque les pas de discrétisation (discrétisation en temps, en espace, ) tendent tous vers 0. Plus précisément, si les données d’une étape du traitement algorithmique sont issues d’une solution exacte, les résultats de ce traitement tendent vers cette solution.

La consistance est une propriété distincte de la convergence, cette dernière étant de portée globale. Sous certaines hypothèses (et en particulier celle de la consistance), le théorème de Lax montre que la stabilité est une condition nécessaire et suffisante pour assurer la convergence.

Définition pour une EDO

On considère une équation différentielle ordinaire (EDO) de la forme :

qu'on essaie de résoudre numériquement par un schéma de la forme

avec une subdivision de et .

L’erreur de consistance relative à une solution exacte de l'équation est donnée par :

La consistance dans le théorème de Lax

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 , soit ) ; des équations aux dérivées partielles du type où 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 ) pour établir une approximation numérique de l’état d’origine. Discrétisation de la variable temporelle (sur avec un pas ) pour entreprendre un processus se déroulant par étapes successives au cours desquelles l’état numérique se transforme.

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

La consistance exige que, pour une fonction régulière, alors soit une bonne approximation de , ceci au sens de la définie par .

Références

Voir aussi

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

Catégorie:Analyse numérique