Department de Physique
Master
Modelisation et Simulation en Mecanique (MSM) &
Mecanique et Ingenierie (MI)
Mini Projet
Résolution de l’équation
elliptique avec la méthode
de SOR
Réalisé par : Demander par :
SAIH LAMYAA Pr. SEHAQUI RACHID
SAIH FATIMA ZAHRA
2023/2024
Table des matières
I- Définition : ......................................................................... 1
II- Application de l’équation elliptique ................................ 1
III- Discrétisation du domaine ............................................... 1
IV- Résolution numérique de l’équation de laplace ............... 2
1) Condition aux limites (bc1) Condition de Dirichlet ......... 6
2) Condition aux limites (bc2) Condition de Neumann 1: ... 7
3) Condition aux limites (bc2) Condition de Neumann 2 : .. 8
V- Méthode de sur relaxation successive SOR ..................... 9
I- Définition :
Une Équation aux Dérivées Partielles (EDP) est une expression mathématique
qui établit une relation fonctionnelle impliquant des dérivées partielles. En
général, si u représente une fonction à valeurs scalaires dépendant des variables
x et y, une EDP prend la forme :
𝜕𝑣 𝜕𝑣 𝜕𝑣
𝑓(𝑣, 𝑥, 𝑦, 𝑧, , , )=0
𝜕𝑥 𝜕𝑦 𝜕𝑧
II- Application de l’équation elliptique
Equation de laplace ∆𝑇 = 0 :
𝜕²𝑇 𝜕²𝑇
+ =0
𝜕𝑥² 𝜕𝑦²
III- Discrétisation du domaine
Une plaque à 2 dimensions d’hauteur H=10 cm et de largeur L=20 cm :
Avec : 𝑛𝑥 = 9 : nombre des nœuds selon x
𝑛𝑦 = 5 : nombre des nœuds selon y
𝐿
∆𝑥 =
𝑛𝑥 − 1
𝐻
∆𝑦 =
𝑛𝑦 − 1
Et : ∆𝑥 = 2,5 𝑐𝑚
∆𝑦 = 2,5 𝑐𝑚
1
H = 10 cm
L = 20 cm
IV- Résolution numérique de l’équation de
laplace
On a l’équation de laplace qui s’écrit sous la forme :
𝜕²𝑇 𝜕²𝑇
+ =0
𝜕𝑥² 𝜕𝑦²
Alors , d’une manière discret :
𝑇𝑖+1,𝑗 + 2𝑇𝑖,𝑗 + 𝑇𝑖−1,𝑗 𝑇𝑖,𝑗+1 + 2𝑇𝑖,𝑗 + 𝑇𝑖,𝑗−1
+ =0
(∆𝑥)² (∆𝑦)²
Trouvons 𝑇𝑖,𝑗 :
2 2 1 1
(− − )𝑇𝑖,𝑗 + (𝑇𝑖+1,𝑗 + 𝑇𝑖−1,𝑗 ) + (𝑇 + 𝑇𝑖,𝑗−1) = 0
(∆𝑥)² (∆𝑦)² (∆𝑥)² (∆𝑦)² 𝑖,𝑗+1
2((∆𝑥)² + (∆𝑦)²) 1 1
(− )𝑇𝑖,𝑗 + (𝑇𝑖+1,𝑗 + 𝑇𝑖−1,𝑗 ) + (𝑇 + 𝑇𝑖,𝑗−1)
(∆𝑥)²(∆𝑦)² (∆𝑥)² (∆𝑦)² 𝑖,𝑗+1
=0
2
(∆𝑥)² (∆𝑦)²
𝑇𝑖,𝑗 = (𝑇𝑖,𝑗+1 + 𝑇𝑖,𝑗−1) + (𝑇 + 𝑇𝑖−1,𝑗 )
2((∆𝑥)² + (∆𝑦)²) 2((∆𝑥)² + (∆𝑦)²) 𝑖+1,𝑗
(∆𝑥)² (∆𝑦)²
Avec : 𝑃 = et 𝑄 =
2((∆𝑥)²+(∆𝑦)²) 2((∆𝑥)²+(∆𝑦)²)
Donc :
𝑇𝑖,𝑗 = 𝑃 ∗ (𝑇𝑖,𝑗+1 + 𝑇𝑖,𝑗−1) + 𝑄 ∗ (𝑇𝑖+1,𝑗 + 𝑇𝑖−1,𝑗 )
Pour résolu le système il devient il nous faut des conditions sur T au bord du
domaine. Les exemples types de conditions aux limites qui seront étudié sont :
Condition de Dirichlet
Condition de Neumann
Voici le programme en fortran :
3
4
5
Le programme se fait avec un test de convergence donné par :
𝑘+1 𝑘
|∑ 𝑇𝑖,𝑗 − ∑ 𝑇𝑖,𝑗 |
𝑘+1 ≤𝜀
∑ 𝑇𝑖,𝑗
Avec : 𝜀 = 10−6
Le calcul s’arrête à une valeur de k=34
1) Condition aux limites (bc1) Condition de Dirichlet
𝑇(𝑥, 0) = 0
𝑇(𝑥, 𝐻 ) = 0
𝑇(0, 𝑦) = 0
𝑇(𝐿, 𝑦) = 100
0
100
6
2) Condition aux limites (bc2) Condition de
Neumann 1:
100
𝜕𝑇
=0 𝜕𝑇
𝜕𝑥 =0
𝜕𝑥
0
𝑇(𝑖, 𝑛𝑦) = 100
𝑇 (𝑖, 1) = 0
𝜕𝑇 −3𝑇𝑖,𝑗 + 4𝑇𝑖+1,𝑗 − 𝑇𝑖+2,𝑗
=
𝜕𝑥 2∆𝑥
Pour i=1 :
4 1
𝑇(1, 𝑗) = 𝑇(2, 𝑗) + 𝑇(3, 𝑗)
3 3
7
𝜕𝑇 3𝑇𝑖,𝑗 − 4𝑇𝑖−1,𝑗 + 𝑇𝑖−2,𝑗
=
𝜕𝑥 2∆𝑥
Pour i=nx :
4 1
𝑇(𝑛𝑥, 𝑗) = 𝑇(𝑛𝑥 − 1, 𝑗) + 𝑇(𝑛𝑥 − 2, 𝑗)
3 3
3) Condition aux limites (bc2) Condition de Neumann
2:
𝜕𝑇
=0
𝜕𝑦
0 100
𝜕𝑇
=0
𝜕𝑦
𝑇 (1, 𝑗) = 0
𝑇(𝑛𝑥, 𝑗) = 100
8
𝜕𝑇 −3𝑇𝑖,𝑗 + 4𝑇𝑖,𝑗+1 − 𝑇𝑖,𝑗+2
=
𝜕𝑦 2∆𝑦
Pour j=1 :
4 1
𝑇(𝑖, 1) = 𝑇(𝑖, 2) + 𝑇(𝑖, 3)
3 3
𝜕𝑇 3𝑇𝑖,𝑗 − 4𝑇𝑖,𝑗−1 + 𝑇𝑖,𝑗−2
=
𝜕𝑦 2∆𝑦
Pour j=ny :
4 1
𝑇(𝑖, 𝑛𝑦) = 𝑇(𝑖, 𝑛𝑦 − 1) + 𝑇(𝑖, 𝑛𝑦 − 2)
3 3
V- Méthode de sur relaxation successive SOR
En anglais, la méthode de relaxation successive (Successive Over Relaxation
Method), abrégée en SOR, représente une variation de la méthode de Gauss-
Seidel utilisée pour résoudre un système d'équations linéaires. Généralement,
l'algorithme SOR converge de manière plus rapide, Gardons le langage Fortran
précédent et l’expression de 𝑇𝑖,𝑗 devient :
𝑘
𝑇𝑖,𝑗 = (1 − 𝜔)𝑇𝑖,𝑗 + 𝜔[𝑃 ∗ (𝑇𝑖,𝑗+1 + 𝑇𝑖,𝑗−1 ) + 𝑄 ∗ (𝑇𝑖+1,𝑗 + 𝑇𝑖−1,𝑗 )]
Avec : 𝜔 : coefficient de sur - relaxation
Si 𝜔 = 1 : méthode de Gauss – Seidel
𝜔 Optimum est donnée par :
2(1 − √1 − 𝜆)
𝜔𝑜𝑝𝑡 =
𝜆
9
Avec :
𝜋 ∆𝑦 𝜋
𝑐𝑜𝑠( ) + ( )²𝑐𝑜𝑠( )
𝑛𝑥 − 1 ∆𝑥 𝑛𝑦 − 1
𝜆=( )²
∆𝑦
1 + ( )²
∆𝑥
Voici le programme en Fortran de SOR :
Nous obtenons que :𝜔opt = 1,329658 ;Avec une convergence rapide du
programme, et le nombre d'itérations et donné par k=16 pour les conditions de
Neumann (1) et k=15 pour les conditions de Dirichlet et Neumann (2)
10
11
12
13