Scribd
Importer
144 vues2 pages

Schémas numériques pour l'équation de transport

Ce document décrit plusieurs schémas numériques pour résoudre l'équation de transport, notamment les schémas aux différences finies explicites, le schéma de Lax-Wendroff et analyse leur précision et leur stabilité. Le document contient également des questions sur la dissipation et la dispersion numérique introduites par ces schémas.

Transféré par

Badreddine Essaidi
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
144 vues2 pages

Schémas numériques pour l'équation de transport

Ce document décrit plusieurs schémas numériques pour résoudre l'équation de transport, notamment les schémas aux différences finies explicites, le schéma de Lax-Wendroff et analyse leur précision et leur stabilité. Le document contient également des questions sur la dissipation et la dispersion numérique introduites par ces schémas.

Transféré par

Badreddine Essaidi
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

PC2 : Schémas numériques pour les équations

hyperboliques linéaires

05 Mai 2011

E XERCICE 1 (S CH ÉMAS NUM ÉRIQUES POUR L’ ÉQUATION DE TRANSPORT.)


On s’intéresse ici à la résolution numérique de l’équation de transport (c > 0) :

∂u ∂u


 +c = 0 x ∈ IR, t > 0,
∂t ∂x (1)

u(x, 0) = u0 (x) x ∈ IR.

Question 1. Etudier l’ordre et la stabilité, par la méthode de Fourier, des schémas


aux différences finies suivants :
un+1
j − ujn−1 unj+1 − unj−1
+c = 0, (2)
2∆t 2h
! !
1 un+1 n
j+1 − uj+1 un+1
j − unj c un+1 n+1
j+1 − uj unj+1 − unj
+ + + = 0. (3)
2 ∆t ∆t 2 h h

Question 2. Indiquer les difficultés spécifiques à l’utilisation pratique des deux schémas.

E XERCICE 2 (L E SCH ÉMA DE L AX -W ENDROFF POUR L’ ÉQUATION DE TRANSPORT )


Soit u une solution régulière de l’équation de transport (1) .

Question 1. Montrer que


∂u c2 ∆t2 ∂ 2 u
u xj , tn+1 = u (xj , tn ) − c∆t (xj , tn ) + (xj , tn ) + O(∆t3 ) .

(4)
∂x 2 ∂x2
Nous proposons la famille de schémas suivante (dépendant du paramètre µ) pour la
résolution numérique de l’équation du transport :

un+1
j − unj c  µ ∆t
unj+1 − unj−1 − 2 unj+1 − 2unj + unj−1 = 0 .

+ (5)
∆t 2h h
Obtenir la valeur de µ pour laquelle la méthode est consistante, au moins, à l’ordre
deux. Nous appellerons au schéma ainsi obtenu le schéma de Lax-Wendroff. Discu-
ter la précision du schéma lorsque c ∆t = h.

Question 2. Calculer le coefficient d’amplification du schéma de Lax-Wendroff et


en déduire sa stabilité sous condition CFL.

1
Question 3. On rappelle que la solution u de l’équation de transport peut s’écrire
en fonction de u0 : Z
u(x, t) = û0 (ξ)eı(x−ct)ξ dξ
R

Ceci montre que l’équation de transport n’est ni dissipative, les ondes planes eıξx ne
sont pas atténuées au cours du temps, et ni dispersive, les ondes planes eıξx se pro-
page toutes à la même vitesse c. Nous allons voir si le schéma numérique conserve
ces propriétés.

On pose Z
n n )ξ
n
uh (x, t ) = û0h (ξ)e−ah (ξ,∆t)t eı(x−ch (ξ,∆t)t dξ
R

Exprimez ah (ξ, ∆t) (la dissipation numérique) et ch (ξ, ∆t) (la dispersion numérique)
en fonction du coefficient d’amplification Ŝh (ξ, ∆t).

Question 4. On appelle taux de dissipation numérique la quantité

ra (ξ, ∆t) = ah (ξ, ∆t),

et erreur sur la vitesse de propagation des ondes planes la quantité

εc (ξ, ∆t) = c − ch (ξ, ∆t).

Étudier le comportement pour h petit (et CFL constante) de ces deux quantités. Com-
menter.

Vous aimerez peut-être aussi