TD 1 : Différences Finies

TD 1 : Différences Finies
UF  Modélisation et calcul scientifique 
Formation ModIA  Modélisation et Intelligence Artificielle 

Nicolas Bertier ([Link]@[Link])

INSA, 4A

1/18
TD 1 : Différences Finies
Etude du schéma Leapfrog
Enoncé

Exercice 1 : Etude du schéma  Leapfrog  pour l’équation

d’advection

Equation d’advection discrétisée à l’aide du schéma Leapfrog :

uin+1 − uin−1 u n − ui−1
n
+ a i+1 =0
2∆t 2∆x

1 Quels sont les ordres en espace et en temps de ce schéma ?

2 Étudiez la stabilité de ce schéma.
3 Quels sont les avantages et inconvénients de ce schéma ?

2/18
TD 1 : Différences Finies
Etude du schéma Leapfrog
Analyse de stabilité

Analyse de stabilité
Ecriture explicite de uin+1 en fonction des autres états :
uin+1 = uin−1 − C ui+1
n n


− ui−1

Analyse de Von Neumann pour l’étude de stabilité :

• domaine périodique ;
• décomposition de la solution discrète en série de Fourier finie,
chaque terme étant de la forme An ejiϕ .

En réinjectant ce mode dans le schéma discret en espace et en temps on

obtient successivement :
 
An+1 ejiϕ = An−1 ejiϕ −C An ej(i+1)ϕ −An ej(i−1)ϕ
An+1 = An−1 − C An ejϕ −An e−jϕ

An+1 = An−1 − CAn ejϕ − e−jϕ

3/18
TD 1 : Différences Finies
Etude du schéma Leapfrog
Analyse de stabilité

Equation du facteur d’amplification

Par définition, on a :
An+1
G=
An
Or, le facteur d’amplification ne dépendant pas de n, on a également :
An An+1
G= =
An−1 An
n
En divisant l’expression précédente par A on obtient successivement :
An+1 An−1
− C ejϕ − e−jϕ


=
An An
1
− C ejϕ − e−jϕ


G =
G
En multipliant cette équation par G, on obtient l’équation du second
degré :
G 2 + G C ejϕ − e−jϕ −1 = 0


| {z }
B
4/18
TD 1 : Différences Finies
Etude du schéma Leapfrog
Analyse de stabilité

Résolution de l’équation du facteur d’amplification

Soit l’équation du second degré pour G :

G 2 + BG − 1 = 0

avec :
B = C ejϕ − e−jϕ = 2jC sin ϕ

Le discriminant ∆ de l’équation s’écrit :

∆ = B 2 + 4 = −4C2 sin2 ϕ + 4 = 4 1 − C2 sin2 ϕ

Les racines de l’équation sont alors :

q
1 − C2 sin2 ϕ


G1 = −jC sin ϕ +
q
1 − C2 sin2 ϕ


G2 = −jC sin ϕ −

5/18
TD 1 : Différences Finies
Etude du schéma Leapfrog
Analyse de stabilité

Stabilité du schéma selon C

On observe que ce schéma est conditionnellement stable :

• Si C > 1 alors les solutions sont purement imaginaires (le
discriminant est négatif) et de module supérieur à l’unité. Le schéma
est donc instable sous cette condition.
• Si C ≤ 1, le discriminant est positif et les solutions comportent alors
une partie imaginaire et une partie réelle. Le carré du module du
facteur d’amplification s’écrit alors (qu’il s’agisse de G1 ou G2 ) :

|G|2 = 1 − C2 sin2 ϕ + C2 sin2 ϕ = 1

|G| = 1

Le module du facteur d’amplification est donc toujours égal à l’unité,

quel que soit C ≤ 1. On parle alors dans ce cas de stabilité neutre.

6/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Exercice 2 : Etude d’une équation de diffusion discrète

Soit l’équation de diffusion :

∂u ∂2 u
−κ 2 =0
∂t ∂x

1 Proposez un schéma DF centré pour la dérivée seconde en espace de

u.
2 Utilisez le schéma d’intégration temporelle d’Euler explicite afin
d’écrire un schéma discret (en espace et en temps)
3 Quelles sont les condition de stabilité de ce schéma ?
4 Mettez en perspective les résultats obtenus par rapport à ce qui a
été étudié en cours sur l’équation d’advection.

7/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Rappel de trigonométrie

Rappel des formules d’Euler :

eiϕ + e−iϕ
cos ϕ =
2
eiϕ − e−iϕ
sin ϕ =
2

Ce qui permet d’écrire après quelques transformations :

1 − cos(2ϕ)
sin2 ϕ =
2

8/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 1

Question 1
Proposez un schéma DF pour la dérivée seconde en espace de u.

Principe de résolution
• De la même manière que l’on a exprimé les dérivées premières de u
par des différences entre des valeurs de u en différents points de
discrétisation, on va maintenant exprimer la dérivée seconde
comme une différence de dérivées premières.
• On choisit ici de prendre un décentrement dans le sens de
l’écoulement :
   
∂u 1 ∂u ∂u
≈ −
∂t ∆x ∂x i+1 ∂x i

9/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 1

Développement de Taylor à l’ordre trois de ui−1 au point i

! !
2 3
∂2 u ∂3 u
 
∂u (∆x) (∆x)
ui−1 = ui − ∆x + − + O(∆x)4
∂x i 2! ∂x 2 3! ∂x 3
i i

Développement de Taylor à l’ordre trois de ui au point i + 1

! !
2 3
∂2 u ∂3 u
 
∂u (∆x) (∆x)
ui = ui+1 −∆x + − +O(∆x)4
∂x i+1 2! ∂x 2 3! ∂x 3
i+1 i+1

10/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 1

Expression des dérivées premières aux points i et i + 1 :

! !
2
∆x ∂2 u ∂3 u
 
∂u ui − ui−1 (∆x)
= − + + O(∆x)3
∂x i ∆x 2! ∂x 2 3! ∂x 3
!i i
!
2
∆x ∂2 u ∂3 u
 
∂u ui+1 − ui (∆x)
= − + + O(∆x)3
∂x i+1 ∆x 2! ∂x 2 3! ∂x 3
i+1 i+1

11/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 1

On peut alors écrire la différence entre les expressions discrètes de ces
dérivées :
   
∂u ∂u ui+1 − 2ui + ui−1
− = + O(∆x)3
∂x i+1 ∂x i ∆x

Soit
   
1 ∂u ∂u ui+1 − 2ui + ui−1
− = + O(∆x)2
∆x ∂x i+1 ∂x i (∆x)2

On notera bien que ce schéma DF pour la dérivée seconde est d’ordre

deux : !
∂2 u ui+1 − 2ui + ui−1
= + O(∆x)2
∂x 2 (∆x)2
i

12/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 2

Question 2
Utilisez le schéma d’intégration temporelle d’Euler explicite afin d’écrire
un schéma discret (en espace et en temps)

Utiliser le schéma d’intégration d’Euler explicite implique :

• De discrétiser la dérivée partielle par rapport au temps comme :
 n
∂u u n+1 − uin
= i + O(∆t)
∂t i ∆t
• D’évaluer les états du second membre de l’équation à l’instant n.
En utilisant la discrétisation centrée de la dérivée seconde établie
précédemment et l’intégration Euler explicite, on obtient le schéma
discret en espace et en temps suivant :
!
n
uin+1 − uin ui+1 − 2uin + ui−1n
=κ 2
∆t (∆x)
13/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 3

Question 3
Quelles sont les condition de stabilité de ce schéma ?

Ecriture explicite de uin+1 en fonction des autres états :

κ∆t n
uin+1 = uin + u − 2uin + ui−1
n


∆x 2 i+1
On note un nouveau groupement combinant les pas d’espace et de
temps, que l’on note D :
κ∆t
D=
∆x 2
Soit :
uin+1 = uin + D ui+1
n
− 2uin + ui−1
n

14/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 3

Afin d’étudier la stabilité de ce schéma, on va mettre en

oeuvre l’analyse de Von Neumann, qui nécessite :
• de se placer sur un domaine périodique ;
• de décomposer la solution discrète en série de Fourier finie.

• On rappelle que chaque terme de la série de Fourier est de la forme

An ejiϕ .
• En réinjectant ce mode dans le schéma discret en espace et en
temps on obtient successivement :
 
An+1 ejiϕ = An ejiϕ +D An ej(i+1)ϕ −2An ejiϕ +An ej(i−1)ϕ
An+1 = An + D An ejϕ −2An + An e−jϕ

15/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 3

On rappelle que le facteur d’amplification G peut s’écrire :

An+1
G=
An
Soit successivement :

1 + D ejϕ −2 + e−jϕ


G =
= 1 + D (cos ϕ + j sin ϕ + cos ϕ − j sin ϕ − 2)
= 1 + 2D (cos ϕ − 1)

Or, on rappelle la formule de trigonométrie (directement issue de la

formule d’Euler) :
1 − cos(2ϕ)
sin2 ϕ =
2
Soit enfin :  ϕ  ϕ
G = 1 − 2D 2 sin2 = 1 − 4D sin2
2 2
16/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 3

On rappelle que le schéma numérique est stable ssi :

||G|| ≤ 1

Soit :
ϕ
|1 − 4D sin2 |≤1
2
Ce qui se traduit par la double inégalité :
ϕ
−1 ≤ 1 − 4D sin2 2 ≤1
−2 ≤ 4D sin2 ϕ2 ≤0
2≥ 4D sin2 ϕ2 ≥0

17/18
TD 1 : Différences Finies
Etude de l’équation de diffusion

Correction exercice 2 / question 3

Le schéma est donc conditionnellement stable
• κ > 0 : le coefficient de diffusion doit être positif ⇒ condition sur
la nature du phénomène physique étudié (pas une condition
numérique)
• D ≤ 12 : condition liée à la discrétisation du problème, qui est
l’équivalent pour une équation de diffusion de la condition sur C
pour une équation de convection.

Quelques remarques supplémentaires :

• On note que, contrairement au schéma centré pour la dérivée
spatiale dans l’équation de convection, la discrétisation centrée de la
dérivée seconde de l’équation est stable.
• Cela est lié à la manière dont l’information se propage avec les deux
types d’équations physiques.

18/18