Université de Tunis El Manar

Département de Génie Civil

Projet d’analyse numérique 

Méthode de décomposition de
domaine en 1D par recouvrement

Elaboré par :
CHROUDI Eya
HALOUI Ines

Encadré par :

Mr. RJAIBI Bader Eddine

Classe : 1ère année Hydraulique

Année universitaire 2018/2019

MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Remerciement

Avant d’entamer ce rapport, nous profitons de l’occasion pour remercier nos

chers professeurs pour avoir créé cette occasion (le mini projet) et nous
permettre ainsi de voir nos acquis purement théoriques par le biais de la
programmation. On remercie ainsi notre encadreur Monsieur Rjaibi pour avoir
bien voulu encadrer ce travail ainsi que pour sa riche contribution et ses
remarques précieuses. Enfin on tient à remercier l’ensemble du corps enseignant
de l’ENIT.

2
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Sommaire

1. Remerciement…………………………………………………………………….....1
2. Introduction ………………………………………………………………………....2
3. Résolution avec condition aux limites de DIRICHLET…………………………....3
4. Méthode de NEWMANN…………………………………………………………...4
5. Méthode de recouvrement………………………………………………………….5
6. Schwarz Alternée…………...………………………………………………………6
7. Schwarz parallèle…………………………………………………………………...7
8. Conclusion………………………………………………………………………….8

3
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Introduction  :

Dans la nature, les systèmes et les phénomènes physiques les plus intéressants
sont aussi les plus complexes à étudier. Ils sont souvent régis par un grand
nombre de paramètres non-linéaires interagissant entre eux (transfert thermique,
météorologie, la turbulence des fluides...).
L'une des solutions est de recourir à une série d'expériences pour analyser les
paramètres et grandeurs du système. Mais les essais peuvent s'avérer très
coûteux (essais en vol, essais avec matériaux rares, instrumentations très
chères...) et ils peuvent être très dangereux (essais nucléaires, environnement
spatial...).

Enfin il peut être difficile de mesurer tous les paramètres : échelles du problème
trop petites (chimie du vivant, couche limite en fluide...) ou trop grandes
(astrophysique, météorologie, géophysique...). On peut aussi construire un
modèle mathématique permettant la représentation du phénomène physique. Ces
modèles utilisent très souvent des systèmes d'équations aux dérivées partielles
(EDP) non-linéaires dont on ne connait pas de solutions analytiques en général.
Il faut résoudre un problème physique avec un zone de recouvrement pour cela
on a utilisé la méthode de Schwartz mais la question qui se pose les quels des
conditions aux limites sont les meilleurs : celle de Dirichlet ou Neumann ?

4
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Résolution avec condition aux limites de

DIRICHLET

-Soit l`équation : -u``+ u = 10*sin (3x)

Cette méthode dépend des conditions au limites : u(0)=0 et u(1)=sin
(3)
-On cherche à résoudre ce système
-On compare la courbe du solution exacte avec celle du la solution
approchée

Résolution numérique :

La méthode des différence consiste à approcher la solution exacte sur

un intervalle [a,b]

aux n+1 points tel que

b−a
a=x 0< x1 … … …< x N =b , avec x i=ih et x i+1=h + x i où h est≤ pas de la discretisation , h=
N

err1=abs (U2sin (3*xm1))

Avec h= (b-a) /N : c’est le pas de discrétisation

5
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

AU N =b
( 2+ h2 ) −1 0 ⋯ 0 U1

1
h2
( −1
⋮
⋮
⋮
( 2+h2 )
⋱
⋱
…
⋱
⋱
⋱
…
⋱
⋱
⋱
−1
⋮
⋮
−1
( 2+h 2 )
) ()
U2
⋮
U N−2
U N−1
=
U0
f 1+

( )
h2
f2
⋮
f N −2
U
f N −1 + 2N
h

6
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

N=10

N=5

7
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

-On remarque que les courbes sont très proches donc l’erreur est presque
négligeable
-Ce qui est démontrée aussi par E=max |Uapp-U|

Passée en fenêtre de commande :

- On remarque que lorsqu’on augmente le nombre d’itération N, les 2 courbes

s’approchent de manière que l’erreur diminue.
-L’ordre de convergence représente la pente de la droite d’équation
p Log (H)+ Log(c)= Log(E) ,
Et celle peut être déterminée directement par la fonction polyfit

8
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Méthode de NEWMANN

- Dans cette méthode on cherche a résoudre la même équation :

- -u’’+ u=10*sin (3x)
- Avec les conditions : u’(0)=3 et u’(1)=3*cos(3)
- Dans cette méthode la résolution du système dépend des dérivées aux
points des extrémités

N=5 N=10

9
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

- On remarque que l’erreur est plus claire que celui de la méthode de

Dirichlet
- En exécutant l’expression : E=max |Uapp-U|

-L’ordre de convergence pour cette méthode déterminée directement par la

fonction polyfit est égale à 1.05

10
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Donc la méthode de Dirichlet converge plus rapidement que la méthode

de Newman

Méthode de recouvrement

-Le domaine est décomposé en 2 domaines avec une zone de recouvrement

qui représente la zone d’intersection .

-Puisque on a démontré que la méthode de Dirichlet converge la meilleur

donc on va l’utiliser dans cette approche pour résoudre les systèmes de Schwarz.
- soient x1 et x2 les points d’extrémités dans la zone de recouvrement

11
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

Schwarz Alternée

(-U1’’)ⁿ⁺₁ +U₁ⁿ⁺1 = 10*sin (3x)

U₁ⁿ⁺1(0)=0
U₁ⁿ⁺1(x1)=U₂ⁿ (x1)

-La décomposition en 2 domaines telle que la solution dans le 1 er domaine pour

l’itération n+1 en x1 dépend de la solution dans le 2eme domaine pour
l’itération n en x1.
(-U2’’)ⁿ⁺1+U₂ⁿ⁺1 = 10*sin (3x)
U₂ⁿ⁺1(xn)=sin(3)

12
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT
U₂ⁿ⁺1(x2)=U₁ⁿ⁺1 (x2)

En fait la décomposition en 2 domaines telle que la solution dans le 2eme

domaine d’itération n+1 en x2depend de la solution dans le 1 er domaine
d’itération n+1 en x2

Le graphique montre :

13
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

err1=abs (U2sin (3*xm1))

14
MÉTHODE Schwarz
1D PARParallèle
DE DÉCOMPOSITION DE DOMAINE EN RECOUVREMENT ENIT

(-U1’’)ⁿ⁺₁ +U₁ⁿ⁺1 = 10*sin(3x)

U₁ⁿ⁺1(0)=0
U₁ⁿ⁺1(x1)=U₂ⁿ (x1)
-La décomposition en 2 domaines tq la solution dans le 1er domaine pour
l’itération n+1 en x1 dépend de la solution dans le 2eme domaine pour
l’itération n pour en x1.

(-U2’’)ⁿ⁺1+U₂ⁿ⁺1 = 10*sin(3x)
U₂ⁿ⁺1(1)=sin(3)
U₂ⁿ⁺1(x2)=U₁ⁿ (x2)

En fait, la décomposition en 2 domaines tq la solution dans le 2éme domaine

pour l’itération n+1 en x2 dépend de la solution dans le 1ér domaine pour
l’itération n en x2.

15
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT

erreur=abs(U2-U1)

Si le nombre d’itérations augmente , les courbes s’approchent de la courbe

du solution exacte , de plus le choix des nœuds x1 et x2 intervient dans
l’approche.

16
MÉTHODE DE DÉCOMPOSITION DE DOMAINE EN 1D PAR RECOUVREMENT ENIT
Conclusion 

Ce projet illustre bien l’importance des méthodes numériques pour la résolution
des problèmes mathématique ainsi que leurs variétés, et permet de constater le
peu de différences concernant les solutions proposées par chaque méthode, d’où
leur efficacité. En effet , il nous a été très utile de travailler sur ce projet,
sachant que d’une part on a pu mieux concevoir l’idée de devoir résoudre un tel
problème mathématique, et d’autre part se familiariser davantage avec un outil
important pour un élève ingénieur et reconnaitre son utilité.

17