Résolution de L'équation Elliptique Avec La Méthode de SOR

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

Résolution de L'équation Elliptique Avec La Méthode de SOR

Transféré par

Droits d'auteur :

Formats disponibles

Résolution de L'équation Elliptique Avec La Méthode de SOR

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Résolution de L'équation Elliptique Avec La Méthode de SOR

Transféré par

Droits d'auteur :

Formats disponibles

Department de Physique

II- Application de l’équation elliptique

III- Discrétisation du domaine

IV- Résolution numérique de l’équation de

1) Condition aux limites (bc1) Condition de Dirichlet

3) Condition aux limites (bc2) Condition de Neumann

V- Méthode de sur relaxation successive SOR

En anglais, la méthode de relaxation successive (Successive Over Relaxation

Avec : 𝜔 : coefficient de sur - relaxation

Voici le programme en Fortran de SOR :

Nous obtenons que :𝜔opt = 1,329658 ;Avec une convergence rapide du

Vous aimerez peut-être aussi