Vous pouvez partager vos connaissances en l’améliorant (comment ?) selon les recommandations des projets correspondants.
Comparaison de solutions numériques de l'équation d'advection linéaire en dimension 1 d'espace et temps obtenues par les schémas de Lax-Friedrichs, Lax-Wendroff, Leapfrog et Upwind. La courbe en noir représente la solution exacte.
Considérons l'équation d'advection en dimension 1 d'espace et temps :
où x et t sont des variables indépendantes, et u(x, 0) est l'état initial connu.
Dans le cas linéaire, on a par définition f(u) = Au , avec A une constante[1],
Avec , le facteur de discrétisation.
Ce schéma linéaire peut être étendu au cas général, non linéaire de multiples façons. L'une d'entre elles consiste à ré-écrire A(u) sous la forme
On obtient alors la formulation conservative du schéma de Lax-Wendroff dans le cas général non-linéaire :
où désigne les matrices jacobiennes évaluées pour les états .
Méthode de Richtmyer
Afin d'éviter le calcul de cette matrice jacobienne, on peut procéder en deux étapes. Dans ce qui suit, on présente la méthode dite de Richtmyer, constituée de deux étapes. Dans un premier temps, on calcule f(u(x, t)) pour un demi pas de temps, tn + 1/2, et en prenant les points xi + 1/2 situés aux interfaces entre deux cellules du maillage. Dans un deuxième temps, les valeurs au temps tn + 1 sont calculées en utilisant celles aux temps tn et tn + 1/2.
Première étape (Lax) :
Deuxième étape :
Méthode de MacCormack
Une autre méthode voisine a été proposée par MacCormack, consistant à utiliser d'abord une différence finie décentrée amont, puis une autre décentrée aval :
Première étape :
Deuxième étape :
Alternative
Première étape :
Deuxième étape :
Références
↑(en) Randy J. LeVeque, Numerical Methods for Conservation Laws, Birkhauser Verlag, , p. 125
(en) P.D Lax et B. Wendroff, « Systems of conservation laws », Commun. Pure Appl Math., vol. 13, no 2, , p. 217–237 (DOI 10.1002/cpa.3160130205)
Michael J. Thompson, An Introduction to Astrophysical Fluid Dynamics, Imperial College Press, London, 2006.
(en) WH Press, SA Teukolsky, WT Vetterling et BP Flannery, Numerical Recipes : The Art of Scientific Computing, New York, Cambridge University Press, , 3e éd., 1235 p. (ISBN978-0-521-88068-8, lire en ligne), « Section 20.1. Flux Conservative Initial Value Problems », p. 1040