Méthodes des Différences Finies

2ème Année - Filière Génie Informatique

Département Mathématiques et Informatique

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 1 / 21
Définitions

Définition 1.1 (EDP) : Une équation aux dérivées partielles fait

intervenir plusieurs variables indépendantes (temps, espace ...), ainsi
que les dérivées partielles de la variable par rapport à ces variables
indépendantes. Par exemple, l’équation : ∂C ∂C
∂t + u ∂x = 0 est une EDP.
La variable dépendante est C, les variables indépendantes sont t et x.
u peut être fonction de t, x et C.
Définition 1.2 (ordre de l’EDP)
L’ordre d’une EDP est l’ordre le plus élevé parmi les dérivées partielles
apparaissant dans l’EDP.
Remarque : dans ce cours, on se limitera aux EDP d’ordre 2.

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 2 / 21
Définitions

Définition 1.3 (EDP)

On dit qu’une EDP est linéaire si elle ne fait intervenir que des
combinaisons linéaires des dérivées partielles de la variable dépendante. On
dit qu’une EDP est quasi-linéaire si elle est linéaire par rapport aux
dérivées d’ordre le plus élevé. Ainsi, l’EDP : a ∂u ∂u
∂t + b ∂y = c
est quasi-linéaire si a, b et c sont des fonctions réelles de x, y et u. Elle
serait linéaire si a, b et c ne dépendaient que de x et y, et à coefficients
constants si a, b et c ne dépendaient plus de x, y et u. De même,

∂2 u ∂2 u ∂2 u
a + b + c =d
∂x 2 ∂x ∂y ∂y 2
est quasi-linéaire si a, b, c et d sont des fonctions réelles de x, y , u, ux , uy .
En dehors de ces critères, l’EDP est non-linéaire.

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 3 / 21
Classification des E.D.P
u telle que u = u (x, y ) solution de :

∂2 u ∂2 u ∂2 u ∂u ∂u
a + 2b + c = f = f (x, y , u, , )
∂x 2 ∂x ∂y ∂y 2 ∂x ∂y

a, b, c ne dépendent que de x, y .

E.D.P : → Linéaire
: → 2ème ordre Mais ” général ”

b 2 − ac > 0 Hyperbolique

b 2 − ac = 0 Parabolique

Pr Hamid El Ouardi * 2GLD * 2 Méthodes des Différences Finies September 25, 2023 4 / 21
Exemples

∂2 u 2
2∂ u
− c =0
∂t 2 ∂x 2
Equation hyperbolique des ondes

∂2 u ∂2 u
+ = ∆u = 0
∂t 2 ∂y 2
Equation elliptique : Potentiel des ondes d’une plaque conductrice

∂u ∂2 u
−d 2 = 0
∂t ∂x
Equation parabolique de diffusion : Temp érature d’une barre métallique.
Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 5 / 21
Conditions aux limites

a) Dirichlet

u est connue sur la frontière

u = u0 sur ∂Ω

b ) Newmann
∂u
= f0 , f0 est connue .
∂n
c ) Mixte
∂u
αu + β = f1 , f1 est connue
∂n

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 6 / 21
Conditions initiales

(Problèmes d’évolution )

u (x, t0 ) = connue sur Ω

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 7 / 21
Un problème aux E.D.P nécessite la donnée de :

1 − Un domaine Ω;
2 − Une équation aux dérivées partielles;
3 − Des conditions aux limites;
4 − Des conditions initiales ( si problème d’évolutions)

L’approximation de l’E.D.P, c’est approcher :

1 − ” approcher ” Ω;
2 − ” Approximer ” les dérivées;
3 − ” Approcher ” les CL+ CI

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 8 / 21
La méthode des différences finies

(MDF) :

C ′ est comment ?
1 − comment approcher Ω;
2 − comment approximer les dérivées;
3 − comment approcher les CL+ CI;
4 − approximation est bonne ? Stabilité et Convergence.

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 9 / 21
Approximation des dérivées

On suppose satisfaites les conditions de validité des calculs qui suivent.

La formule de Taylor-Lagrange donne

h2 ′′ h3 hn
f (x + h ) = f (x ) + hf ′ (x ) + f (x ) + f ′′′ (x ) + f (n) (x ) + o (hn ).
2! 3! n!

∂f ∂f h 2 ∂2 f
f (x + h, y + k ) = f (x, y ) + h (x, y ) + k (x, y ) + (x, y ) +
∂x ∂y 2! ∂x 2
k 2 ∂2 f ∂2 f
( x, y ) + hk (x, y )+
2! ∂y 2 ∂x ∂y
  (n −1)
1 ∂f ∂f
.. h (x, y ) + k (x, y ) + o ((|h| + |k |)n ).
(n − 1)! ∂x ∂y

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 10 / 21
Par addition et soustraction nous obtenons
∼ f (x + h ) − f (x − h )
f ′ (x ) − ,
2h
∼ f (x + h ) − 2f (x ) + f (x − h )
f ′′ (x ) −
h2
∂f ∼ f (x + h, y ) − f (x − h, y )
(x, y ) − ,
∂x 2h
∂f ∼ f (x, y + k ) − f (x, y − k )
(x, y ) −
∂y 2k

∂2 f ∼ f (x + h, y ) − 2f (x, y ) + f (x − h, y )
2
(x, y ) −
∂x h2

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 11 / 21
 ∂2 f ∼ f (x, y + k ) − 2f (x, y ) + f (x, y − k )
2
(x, y ) −
∂y k2

∂2 f
∂x ∂y (x, y ) est approchée par

f (x +h,y +k )+f (x −h,y −k )−f (x −h,y +k )−f (x +h,y −k )

4hk

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 12 / 21
Condition de Dirichlet

On numérote ensuite les points de R et on sépare les numérotations sur le

bord de celle de l’intérieur.
Lorsque M décrit R, l’ensemble de toutes les équations précedentes ( il y
en a autant que d’élements de R ) est équivalent à un système linéaire de
la forme (3) AU = B. Le système linéaire (3) peut-être résolu par des
méthodes de Gauss s’il y a un nombre restreint d’inconnues. Si le nombre
d’inconnues est grand et la matrice est symétrique définie positive, on
pourra utiliser une méthode de Choleski, si la matrice n’est pas
symétrique, des méthodes itératives s’imposent.

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 13 / 21
Condition de Neumann

Lorsque les conditions sont de type (b), la valeur de l’inconnue n’est plus
donnée partout sur le bord mais seulement sur une partie.
∂u
∂n = f (s ), s ∈ ∂Ω

∆u = 0 dans Ω

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 14 / 21
 
u (C )−u (B )
∂u
∂n A = 2h / Ω!
+ o (h2 ), or C ∈

∂u h 2 ∂2 u
u (C ) = u (A) − h ( )A + ( 2 )A + o (h 3 )
∂n 2 ∂n
 2  2
Mais (∆u |A = 0 = ∂∂nu2 + ∂∂xu2
A A

∂2 u ∂2 u u (D ) − 2u (A) + u (E )
 
⇒ =− =−
∂n2 A ∂x 2 A h2

+o (h 2 )

∂u
d’où : ∂n = f (s ) est approchée par :

1
[4u (A) − 2u (B ) − u (D ) − u (E )]
2h

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 15 / 21
Etude de cas n◦ 1
A) Supposons que Ω = ]0, a[ × ]0, b [ et soient N et M deux entiers et
posons h = Na , k = M
b
, alors :

R = Rhk = [Mij : Mij = (ih, jk ); i = 0, ......N; j = 0, .....M ] .

Les points frontaliers du réseau sont donc les points

Mij ; i ∈ {0, N } ∀j; j ∈ {0, M } ∀i.

Si la condition aux frontières est u (x, y ) = g (x, y ); x, y ∈ Γ , on posera

alors ui,j = gi,j pour i ∈ {0, N } ∀j; j ∈ {0, M } ∀i.
Si la condition à la frontière est ∂u
∂n = g (0, y ) alors la dérivée normale est
la dérivée en x.
Nous avons donc :
∂u ∂2 u
u1,j = u0,j + h + h2 2∂x 2 + o (h3 ).
∂n

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 16 / 21
Les dérivées étant prises aux points (0; jk ). D’où (*)

∂2 u
 
2 ∂u
∂x 2
= 2 u1,j − u0,j −h + o (h ).
h ∂n

∂2 u
On peut approcher ∂x 2
au point (0; jk ) par

2
(u1,j − u0,j − hg0,j ) .
Pr Hamid El Ouardi * 2GLD *
h2 Méthodes des Différences Finies September 25, 2023 17 / 21
Si la condition à la frontière est : α ∂u
∂n (0, y ) + βu (0, y ) = g (0, y ) et en
2
reprenant (*), nous pouvons donc approcher ∂∂x 2u au point (0; jk ) par :
 
2 βh h
u1,j − (1 − )u0,j − g0,j .
h2 α α

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 17 / 21
Etude de cas n◦ 2

B) Considérons le problème suivant sur Ω = ]0, L[ × ]0, L[

 ∆u (x, y ) = f (x, y ); (x, y ) ∈ Ω



u (x, y ) = g (x, y ); (x, y ) ∈ ∂Ω



Soit un entier N, on obtient le réseau

 
L
R = Rhk = Mij : Mij = (ih, jk ); i, j = 0, ......N; h = .
N
Pour ne pas avoir des matrices trop grandes, nous prendrons N = 4.
Nous numérotons les points comme indiqué sur la figure 2.
Nous noterons fi la valeur de f au point numéroté i et gi la valeur de g au
point numéroté i (i = 10, 25) et Vi l’approximation au point i.

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 18 / 21
Nous avons l’approximation

u (x − h, y ) + u (x + h, y )+
∆u (x, y ) =
h2
u (x, y − h ) + u (x, y + h ) − 4u (x, y )
h2
= f (x, y )

Nous obtenons donc, en écrivant cette équation aux points

i = 1, .........., 9 :

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 19 / 21
 −4 1 0 1 0 0 0 0 0 V1
  

 1 −4 1 0 1 0 0 0 0 
 V2 


 0 1 −4 0 0 1 0 0 0 
 V2 


 1 0 1 −4 1 0 1 0 0 
 V3 


 0 1 0 1 −4 1 0 1 0 
 V4 


 0 0 1 0 1 −4 0 0 1 
 V5 


 0 0 0 1 0 0 −4 1 0 
 V6 

 0 0 0 0 1 0 1 −4 1  V7 
0 0 0 0 0 1 0 1 −4 V9

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 20 / 21
 h2 f1 − g11 − g25
 
 h2 f2 − g12 
2
 
 h f3 − g13 − g15 
h2 f4 − g24
 
 
 
=
 h 2 f5 

 2
h f6 − g16 
 
 2
h f7 − g21 − g23 
 
 h2 f8 − g20 
2
h f9 − g17 − g19

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 21 / 21
Ce qui est équivalent au système linéaire :

AV = B

Pr Hamid El Ouardi * 2GLD * Méthodes des Différences Finies September 25, 2023 21 / 21