Utilisateur:Gigowatts/Schéma de Lax-Wendroff

Une page de Wikipédia, l'encyclopédie libre.
Comparaison de solutions numériques de l'équation d'advection linéaire en dimension 1 d'espace et temps obtenues par les schéma de Lax-Friedrichs, Lax-Wendroff, Leapfrog et Upwind. La courbe en noir représente la solution exacte.

Le schéma de Lax–Wendroff, d'après Peter Lax et Burton Wendroff, est défini en analyse numérique comme une technique de résolution numérique des équations aux dérivées partielles de type hyperbolique, basée sur la méthode des différences finies. Ce schéma est d'ordre 2 à la fois en espace et en temps.

Considérons l'équation d'advection en dimension 1 d'espace et temps:

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]

Ce schéma linéaire peut être étendu au cas général, non-linéaire de multiples façcons. L'une d'entre elle 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:

désigne les matrices jacobiennes évaluées pour les états .

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(xt)) 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:

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:

Alternativement, Première étape:

Deuxième étape:

Références[modifier | modifier le code]

  1. LeVeque, Randy J. Numerical Methods for Conservation Laws", Birkhauser Verlag, 1992, p. 125.
  • 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, 3rd, (ISBN 978-0-521-88068-8), « Section 20.1. Flux Conservative Initial Value Problems », p. 1040

Voir aussi[modifier | modifier le code]

{{Palette|Équations aux dérivées partielles - Méthodes numériques}} {{Portail|analyse}} [[Catégorie:Équations différentielles numériques]] {{ébauche}}