Ecole Nationale Supérieure des Travaux Publics

Département des sciences fondamentales

Cours de modélisation et simulation numérique

Niveau Master 1

RESOLUTION DE PROBLEMES PAR LA METHODE DE

DIFFERENCES FINIES IMPLICITES ET LA METHODE DE
ADAMS MOULTON D’ORDRE 4

PROJET DE SIMULATION NUMERIQUE EN MATLAB

Travail réalisé par : Chargé du cours :

DJIDDO AHMAT MACKI Dr. VOTSIA DJOKATA
DJEKOSSEME KEMLENGAR SYMPLICE
AHMAT ADANO ADOUDOU

ANNEE ACADEMIQUE 2023-2024

Table des matières
...................................................................................................................................................................
I. INTRODUCTION : ....................................................................................................................... 2
I. Partie I ............................................................................................................................................ 4
I . PRESENTATION DES PROBLEMES ABORDES ..................................................................... 4
I.1. Problème 1 : Equation de diffusion de réaction ....................................................................... 4
I.2. Problème 2 : Equation de transport .......................................................................................... 5
Partie II .................................................................................................................................................. 7
II. PRESENTATION DES METHODES NUMERIQUES UTILISEES POUR LA
RESOLUTION : .................................................................................................................................... 7
II.1. La Méthode des différences finies implicites : ........................................................................ 7
II.2. La Méthode ADAMS MOULTON d'ordre 4 : ........................................................................ 7
II. 3. Discrétisation et Résolution par la méthode des différences finies implicites :................... 8
3.1. Problème 1 : Equation de diffusion de réaction....................................................................... 8
3.2. Problème 2 : Equation de transport ....................................................................................... 12
III. Discrétisation et Résolution par la méthode de ADAMS MOULTON q=4 : .................. 15
IV. LES ALGORITHMES DE DEUX METHODES NUMERIQUES UTILISEES
PERMETTANT DE RESOUDRE LES DEUX PROBLEMES SOUS MATLAB : ....................... 17
III.1. Algorithme de la méthode des différences finies implicites : ............................................ 17
III.2. Algorithme de la méthode de ADAMS MOULTON ordre q=4 : ....................................... 18
V. LES FIGURES ISSUES DE LA RESOLUTION PAR LES DEUX METHODES
NUMERIQUES ................................................................................................................................... 18
IV.1. Les figures issues de la méthode des différences finies implicites : .................................... 19
IV.1.a. Courbe de l’équation de diffusion-réaction.................................................................. 19
IV.1.b. Courbe de l’équation de transport................................................................................ 19
IV.2. Les figures issues de la méthode de ADAMS MOULTON ordre q=4 ................................ 20
IV.1.a. Courbe de l’équation de diffusion-réaction...................................................................... 20
IV.1.b. Courbe de l’équation de transport................................................................................ 20
VI. COMPARAISON DE DEUX METHODES NUMERIQUES DE RESOLUTION ........... 21
VII. CONCLUSION ........................................................................................................................ 22
VIII. Annexes................................................................................................................................. 23
I. Les programmes MATLAB............................................................................................................. 23
I.1. Programme MATLAB méthode des différences finies .......................................................... 23
I.1.1. Programme ode45MDF équation de diffusion .................................................................... 23
I.1.2. Programme ode45MDF équation de transport ................................................................... 24
I.2. Programme MATLAB méthode de ADAMS MOULTON Q=4 ............................................ 24
I.2.1. Programme ode45 ADM4 équation de diffusion ................................................................. 24
I.2.2. Programme ode45 ADM équation de transport .................................................................. 26

1
 I. INTRODUCTION :

Dans le cadre de la formation académique du niveau master au sein de

l'Ecole Nationale Supérieure des Travaux Publics, les étudiants sont appelés
à traiter des projets de recherche relatifs aux matières enseignées. C'est dans
ce cadre qu’un projet de modélisation et simulation numérique nous ai
assigné en qualité d’étudiants de master 1 de génie civil. Ce projet doit d'être
traité par binôme.
De nos jours, l’enseignement des mathématiques appliquées ne peut se
concevoir sans l’expérimentation numérique. En effet, seule la mise en
œuvre d’une méthode sur un exemple concret permet d’en mieux
comprendre les qualités (précision et rapidité de calcul) et les défauts (coût
et limites d’utilisation). Bien sûr, ces éléments distinctifs peuvent être
appréhendés d’un point de vue théorique ; mais la complexité de
programmation et la variété des champs d’application, critères finalement
déterminants pour le succès d’une méthode, ne peuvent être mesurées que
par l’expérimentation.
Ce présent rapport réunit deux exemples d’application de l’analyse
numérique, mises en œuvre pour la résolution de problèmes (équation de
transport et celle de diffusion-réaction) concrets issus de la physique, la
mécanique, la chimie etc. Ces problèmes ont été choisis pour permettre aux
étudiants de mieux cerner la modélisation et la simulation numérique
considérée comme matière enseignée mais aussi comme outils permettant à
ceux derniers de maitriser la programmation sous Matlab
Au fait, le projet consiste précisément à résoudre deux problèmes (problème
de diffusion-réaction et problème de transport) en utilisant deux méthodes
numériques (méthode des différences finies implicites et la méthode de
Adams-Moulton d'ordre 4 sous MATLAB)
Dans ce travail, il sera question de présenter en premier lieu les problèmes
de transport et de diffusion-réaction qui ont été choisis pour cette étude,
expliquer en second lieu les principes et les équations mathématiques sous-
jacentes à chaque méthode numérique utilisée. En troisième approche, les
étapes de mise en œuvre de ces méthodes dans l'environnement MATLAB
seront mentionner. En quatrième et dernière position les analyses des
2
résultats obtenus pour chaque méthode et la comparaison entre les méthodes
utilisées offrant une précision, jugeant l'efficacité ainsi que la stabilité de
celles-ci se présentent en guise de conclusion.
Le cheminement suivi pour traiter ce travail se dessine en termes de résultats
attendus comme suit :
➢ Présentation de problèmes et méthodes numériques utilisées
➢ Les algorithmes pour la résolution
➢ Les programmes MATLAB
➢ Les figures issues de la résolution
➢ Comparaison de deux méthodes numériques utilisées
De façons structurelle, le présent rapport est structuré en deux parties en plus
de l’introduction et la conclusion.
La première partie porte sur la discrétisation et la résolution des problèmes ;
La deuxième partie présente les problèmes et méthodes numériques, les
algorithmes et programmes MATLAB et une comparaison de ces deux
méthodes.

3
 I. Partie I
I . PRESENTATION DES PROBLEMES ABORDES

De nos jours, la connaissance des méthodes numériques servant à modéliser

un problème physique s’avère indispensable pour un ingénieur. La
résolution des équations différentielles ordinaires mais surtout celle des
équations aux dérivées partielles à l’aide des méthodes numériques permet
d’anticiper via les modèles plus ou moins complexes la faisabilité, la
durabilité ainsi que l’efficacité d’un projet quel que soit sa nature. C’est
pourquoi, pour mieux cerner les méthodes et modèles numériques, nous, au
cours de ce travail, allons amorcer des exemples particuliers que sont les
équations de diffusion-réaction et de transport dans le but de traiter notre
travail. Notons que les équations de diffusion et de transport sont des types
courants d'équations utilisées pour modéliser divers phénomènes physiques
dans des domaines tels que la physique, la chimie, l'ingénierie et les sciences
appliquées. Chaque problème a un ordre, une classe, une forme générale
mais surtout chaque problème modélise un ou plusieurs phénomènes.

I.1. Problème 1 : Equation de diffusion de réaction

𝑢𝑡 = 𝜀𝑢𝑥𝑥 + 10𝑢 − 1 < 𝑥 < 1, 𝑡 > 1

{ 𝑢(−1, 𝑡) = 0, 𝑢(1, 𝑡) = 1
𝑢(0, 𝑥) = 0.53𝑥 + 0.47sin (−1.5𝜋)

De manière générale, un système de diffusion-réaction est un modèle

mathématique qui décrit l’évolution des concentrations d’une ou plusieurs
substances spatialement distribuées et soumises à deux processus : un
processus de réactions chimiques locales, dans lequel les différentes
substances se transforment, et un processus de diffusion qui provoque une
répartition de ces substances dans l'espace .

De façon spécifique, il s'agit d'une équation décrivant la diffusion d'une

quantité (u) dans une dimension spatiale (x) avec une réaction chimique
(10u). Les conditions aux limites sont données par : u (-1, t) = 0 et u (1, 0)
=1, et la condition initiale est u (0, x) = 0.53x + 0.47sin(-1.5π). La simulation
est réalisée pour ε = 10^ (-2).

4
 ✓ Ordre : 2 (dérivées d'ordre 2 en espace) ;
✓ Classe : Équation parabolique ;
✓ Forme générale : 𝒖𝒕 = 𝜺𝒖𝒙𝒙 + 𝐜𝒖 ,
Où u est une fonction de deux variables (x, t), 𝑢𝑥𝑥 est la dérivée
seconde de u par rapport à x, 𝑢𝑡 est la dérivée de u par rapport à t, et c
est une constante ;
✓ Phénomènes modélisés : Ce type d'équation de diffusion peut
modéliser la propagation de la chaleur, la diffusion d'une substance
dans un milieu, ou la répartition d'une quantité dans l'espace et le
temps.

I.2. Problème 2 : Equation de transport

𝑢𝑡 + 𝑢𝑥 = 0; 𝑥 𝜖 [−10,10], 𝑡 > 1
{𝑢(𝑥, 0) = 1 si 𝑥 𝜖 [−10,10], 0 𝑠𝑖𝑛𝑜𝑛
𝑢(−10, 𝑡) = 𝑢(10, 𝑡) = 𝑠𝑖𝑛𝑡

De façon générale, l’équation de transport appelée aussi équation de

continuité est une équation qui décrit le transport d’une certaine quantité.
Elle est particulièrement simple et puissante lorsqu’elle est appliquée à
une quantité conservée, mais elle peut être généralisée pour s’appliquer à
n’importe quelle quantité extensive. Étant donné que la masse, l’énergie, la
quantité de mouvement, la charge électrique et d’autres quantités naturelles
sont conservées dans leurs conditions appropriées respectives, une variété
de phénomènes physiques peuvent être décrits à l’aide d’équations de
transport.

De manière spécifique, il s'agit d'une équation qui décrit le transport d'une

quantité (u) dans une dimension spatiale (x) sans aucune réaction. Les
conditions aux limites sont données par : u (-10, t) = u (10, t) = sin(t), et la
condition initiale est u (x, 0) = 1 si x appartient à l'intervalle [-10, 10], sinon
u (x, 0) = 0.

✓ Ordre : 1 (dérivées d'ordre 1 en espace et en temps) ;

✓ Classe : Équation hyperbolique ;
✓ Forme générale : 𝑢𝑡 + c𝑢𝑥 = 0,

5
 où u est une fonction de deux variables (x, t), 𝑢𝑥 est la dérivée de u par
rapport à x, 𝑢𝑡 est la dérivée de u par rapport à t, et c est une constante ;

✓ Phénomènes modélisés : Ce type d'équation de transport peut

modéliser le transport d'une quantité ou d'une information à travers un
milieu. Par exemple, cela peut représenter le mouvement d'une
particule dans un fluide, la propagation d'une onde ou le transport de
quantités physiques telles que la masse, la vitesse ou la concentration.

Caractéristique Équation de transport Équation de diffusion

Nature de l'équation Équation hyperbolique Équation parabolique

Comportement général Propagation d'ondes et de Diffusion et lissage

discontinuités

Déplacement de quantités Transport des quantités sans Dispersion et mélange des quantités
modification de leur forme

Vitesse de propagation La vitesse de propagation est définie La diffusion se produit à une vitesse
par le terme convectif proportionnelle au gradient spatial

Conditions initiales Nécessite une condition initiale et des Nécessite une condition initiale et des
conditions aux limites spécifiques conditions aux limites spécifiques

Exemples d'applications Transport de polluants dans Diffusion thermique, diffusion de

l'atmosphère, modèles de trafic substances chimiques dans un milieu
routier, etc. poreux, etc.

Tableau comparatif des problèmes (de diffusion et de transport)

6
 Partie II

II. PRESENTATION DES METHODES NUMERIQUES

UTILISEES POUR LA RESOLUTION :

II.1. La Méthode des différences finies implicites :

La méthode des différences finies consiste à approximer les dérivées des
équations de la physique au moyen des développements de Taylor et se
déduit directement de la définition de la dérivée . C’est une méthode qui
discrétise l'espace et le temps en points et approxime les dérivées spatiales
et temporelles par des différences finies. Les équations résultantes sont
résolues implicitement à chaque pas de temps, ce qui permet de traiter
efficacement les problèmes avec des conditions aux limites complexes ou
des coefficients qui varient dans le temps. Le schéma implicite présente
l’avantage d’être toujours stable et donc de toujours converger vers la
solution. Exemple d’un schéma centré en espace et d’un schéma décentré en
temps.
𝐝𝐔 𝐔 𝐧+𝟏 − 𝐔 𝐧 𝐔𝐢+𝟏 − 𝟐𝐔𝐢 + 𝐔𝐢−𝟏
𝐔𝐭 = = ; 𝐔𝐱𝐱 =
𝐝𝐭 𝐝𝐭 (𝐝𝐱)𝟐

II.2. La Méthode ADAMS MOULTON d'ordre 4 :

La méthode de Adams Moulton q=4 est une méthode aux pas multiples qui
n’est pas autodéterminante. C’est-à-dire, pour déterminer les valeurs
initiales du problème, on utilise une parmi les méthodes à un pas (Euler,
Runge Kutta d’ordre 4, …) Elle permet d'obtenir une meilleure précision par
rapport à la méthode des différences finies implicites. Cependant, en
pratique, la méthode de Adams Moulton est utilisée conjointement avec la
méthode de Adams Bashforth de même ordre pour construire un schéma
prédicteur-correcteur. Le schéma de Adams Moulton q=4 est donné comme
suit :
Pas Ordre ADAMS MOULTON q=4

q=4 p=5 𝒉
𝒚𝒌+𝟏 = 𝒚𝒌 + (𝟐𝟓𝟏𝒇𝒌+𝟏 + 𝟔𝟒𝟔𝒇𝒌 − 𝟐𝟔𝟒𝒇𝒌−𝟏 + 𝟏𝟎𝟔𝒇𝒌−𝟐 − 𝟏𝟗𝒇𝒌−𝟑 )
𝟕𝟐𝟎

7
II. 3. Discrétisation et Résolution par la méthode des différences finies
implicites :

3.1. Problème 1 : Equation de diffusion de réaction

𝑢𝑡 = 𝜀𝑢𝑥𝑥 + 10𝑢
−1 < 𝑥 < 1, 𝑡 > 1
{ 𝑢(−1, 𝑡) = 0, 𝑢(1, 𝑡) = 1
𝑢(0, 𝑥) = 0.53𝑥 + 0.47sin (−1.5𝜋)

Discrétisation partielle en espace

−1 < 𝑥 < 1 est le domaine de l′espace

1−(−1)
 𝑑𝑥 =
𝑁𝑋

2
 𝑑𝑥 = 𝑝𝑎𝑠 𝑑𝑒 𝑙′𝑒𝑠𝑝𝑎𝑐𝑒
𝑁𝑋

𝑥𝑖 = 𝑖𝑑𝑥 𝑖 = 0, … , 𝑁𝑋

𝑈𝑖 = 𝑈(𝑥𝑖 , 𝑡) ;

𝑙𝑎 𝑑𝑖𝑠𝑐𝑟é𝑡𝑖𝑠𝑎𝑡𝑖𝑜𝑛 𝑑𝑒 𝑈𝑥𝑥 𝑝𝑎𝑟 𝑙𝑒𝑠 𝑑𝑖𝑓𝑓é𝑟𝑒𝑛𝑐𝑒𝑠 𝑓𝑖𝑛𝑖𝑒𝑠 𝑠𝑐ℎé𝑚𝑎 𝐷𝐴 𝑒𝑠𝑡:

𝑑𝑈 𝑈𝑖+1 −2𝑈𝑖 +𝑈𝑖−1
𝑢𝑡 = et 𝑈𝑥𝑥 =
𝑑𝑡 (𝑑𝑥)2

𝐸𝑛 𝑟𝑒𝑚𝑝𝑙𝑎ç𝑎𝑛𝑡 𝑢𝑡 𝑒𝑡 𝑢𝑥𝑥 𝑑𝑎𝑛𝑠 (1) 𝑜𝑛 𝑎𝑢𝑟𝑎 ∶

𝑑𝑈 𝜀
= (𝑈 − 2𝑈𝑖 + 𝑈𝑖−1 ) + 10𝑈𝑖
𝑑𝑡 (𝑑𝑥)2 𝑖+1

𝑑𝑈 𝜀 (𝑑𝑥)2
= (𝑈 − (2 − 10 )𝑈𝑖 + 𝑈𝑖−1 )
𝑑𝑡 (𝑑𝑥)2 𝑖+1 𝜀

8
 Itération

𝑑𝑈 𝜀 (𝑑𝑥)2
𝑖 =1 = (𝑑𝑥)2 (𝑈2 − (2 − 10 )𝑈1 + 𝑈0 )
𝑑𝑡 𝜀

𝑑𝑈 𝜀 (𝑑𝑥)2
𝑖 =2 = (𝑈 − (2 − 10 )𝑈2 + 𝑈1 )
𝑑𝑡 (𝑑𝑥)2 3 𝜀
𝑑𝑈 𝜀 (𝑑𝑥)2
𝑖 =3 = (𝑈 − (2 − 10 ) 𝑈3 + 𝑈2 )
𝑑𝑡 (𝑑𝑥)2 4 𝜀
⋮
𝑑𝑈 𝜀 (𝑑𝑥)2
𝑖 = 𝑁𝑋 − 1  = (𝑑𝑥)2 (𝑈𝑁𝑋 − (2 − 10 )𝑈𝑁𝑋−1 + 𝑈𝑁𝑋−2 )
𝑑𝑡 𝜀

(𝑑𝑥)2
𝑝𝑜𝑠𝑜𝑛𝑠 𝛼 = −2 + 10 , 𝑎𝑙𝑜𝑟𝑠 𝑜𝑛 𝑎 𝑙𝑎 𝑚𝑎𝑡𝑟𝑖𝑐𝑒
𝜀
𝛼 1 0 0… 𝑈1
1 𝛼 1 0… 𝑈1
0 1 𝛼 1 0… 𝑈1
𝑑𝑈 𝜀 ⋱ ⋱ ⋱
=
𝑑𝑡 (𝑑𝑥) 2 ⋮
⋮ ⋱ ⋱ ⋱ 1
0… 0 1 𝛼 𝑈𝑁𝑋−2
[ ] [ 𝑈𝑁𝑋−1 ]
𝛼 1 0 0… 𝑈1
1 𝛼 1 0… 𝑈1
0 1 𝛼 1 0… 𝑈1
⋱ ⋱ ⋱
𝑀= 𝑈=
⋮
⋮ ⋱ ⋱ ⋱ 1
0… 0 1 𝛼 𝑈𝑁𝑋−2
[ ] [ 𝑈𝑁𝑋−1 ]
𝑑𝑈 𝜀
= 𝑀𝑈
𝑑𝑡 (𝑑𝑥)2
𝑈(0) = 0.53𝑥 + 0.47sin (−1.5𝜋𝑥)

La résolution de ce problème est faite à l’aide de l’ODE45 sous MATLAB.

9
 Discrétisation en temps
0 < 𝑡 < T est le domaine temporel

𝑇−(0)
 𝑑𝑡 =
𝑁𝑇

𝑇
 𝑑𝑡 = 𝑝𝑎𝑠 𝑑𝑒 𝑙′𝑒𝑠𝑝𝑎𝑐𝑒
𝑁𝑇

𝑡𝑛 = 𝑛𝑑𝑡 , 𝑛 = 0, … , 𝑁𝑇 − 1

𝑈 𝑛 = 𝑈(𝑡𝑛 ) ;

𝑑𝑈 𝑈 𝑛+1 − 𝑈 𝑛
𝑈𝑡 = =
𝑑𝑡 𝑑𝑡
𝑑𝑈 𝜀 𝑈 𝑛+1 −𝑈 𝑛 𝜀
= (𝑑𝑥)2 𝑀𝑈  = (𝑑𝑥)2 𝑀𝑈
𝑑𝑡 𝑑𝑡

𝜀𝑑𝑡
 𝑈𝑛+1 = 𝑈 𝑛 + 2
𝑀𝑈 𝑛
(𝑑𝑥)

Itération en fonction de n :

𝜀𝑑𝑡
𝑛 = 0  𝑈1 = 𝑈 0 + 𝑀𝑈 0
(𝑑𝑥)2

𝜀𝑑𝑡
𝑛 = 1  𝑈 2 = 𝑈1 + 2
𝑀𝑈1
(𝑑𝑥)

𝜀𝑑𝑡
𝑛 = 𝑁𝑇 − 1  𝑈𝑁𝑇 = 𝑈 𝑁𝑇−1 + 2
𝑀𝑈 𝑁𝑇−1
(𝑑𝑥)

𝜀𝑑𝑡
𝑈 𝑁𝑇 = (𝐼𝑁 + 𝑀)𝑈 𝑁𝑇−1
(𝑑𝑥)2

La forme matricielle est donnée comme suit :

10
 𝜀𝑑𝑡
𝐼𝑁 0 … 𝑈1 𝑈0 + 𝑀𝑈 0
(𝑑𝑥 ) 2

−𝑁 𝐼𝑁 0 … 𝑈2 0

0−𝑁 𝐼𝑁 0 … 𝑈3
⋮ ⋱ ⋱ ⋱ ⋱ ⋱ ⋮
=

⋱ ⋱ 0
… 0 −𝑁 𝐼𝑁
[ ] [𝑈 𝑁𝑇−1 ] [ 𝐼𝑁 ]

𝐼𝑁 0 … 𝑈1

−𝑁 𝐼𝑁 0 … 𝑈2

0−𝑁 𝐼𝑁 0 … 𝑈3
⋮ ⋱ ⋱ ⋱ ⋱ ⋱
𝐴= ;𝑈 =
⋮

⋱ ⋱
… 0 −𝑁 𝐼𝑁
[ ] [𝑈 𝑁𝑇−1 ]
𝜀𝑑𝑡
𝑈0 + 2
𝑀𝑈 0
(𝑑𝑥 )
𝐵= 0
⋮
0
[ 𝐼𝑁 ]

𝐴𝑈 = 𝐵  𝑈 = 𝐴−1 𝐵

La résolution de ce problème est faite à l’aide de l’ODE45 sous MATLAB.

11
3.2. Problème 2 : Equation de transport

𝑢𝑡 + 𝑢𝑥 = 0
𝑥 𝜖 [−10,10], 𝑡 > 1
{
𝑢(𝑥, 0) = 1 si 𝑥 𝜖 [−10,10], 0 𝑠𝑖𝑛𝑜𝑛
𝑢(−10, 𝑡) = 𝑢(10, 𝑡) = 𝑠𝑖𝑛𝑡

Discrétisation partielle en espace

−10 < 𝑥 < 10 est le domaine de l′espace

10−(−10)
 𝑑𝑥 =
𝑁𝑋

20
 𝑑𝑥 = 𝑝𝑎𝑠 𝑑𝑒 𝑙′𝑒𝑠𝑝𝑎𝑐𝑒
𝑁𝑋

𝑥𝑖 = 𝑖𝑑𝑥 𝑖 = 0, … , 𝑁𝑋

𝑈𝑖 = 𝑈(𝑥𝑖 , 𝑡) ;

𝑙𝑎 𝑑𝑖𝑠𝑐𝑟é𝑡𝑖𝑠𝑎𝑡𝑖𝑜𝑛 𝑑𝑒 𝑈𝑥𝑥 𝑝𝑎𝑟 𝑙𝑒𝑠 𝑑𝑖𝑓𝑓é𝑟𝑒𝑛𝑐𝑒𝑠 𝑓𝑖𝑛𝑖𝑒𝑠 𝑠𝑐ℎé𝑚𝑎 𝐷𝐴 𝑒𝑠𝑡:

𝑑𝑈 𝑈𝑖 −𝑈𝑖−1
𝑢𝑡 = et 𝑈𝑥 =
𝑑𝑡 𝑑𝑥

𝐸𝑛 𝑟𝑒𝑚𝑝𝑙𝑎ç𝑎𝑛𝑡 𝑢𝑡 𝑒𝑡 𝑢𝑥𝑥 𝑑𝑎𝑛𝑠 (2) 𝑜𝑛 𝑎𝑢𝑟𝑎 ∶

𝑑𝑈 1
+ (𝑈 − 𝑈𝑖−1 ) = 0
𝑑𝑡 𝑑𝑥 𝑖
𝑑𝑈 1
= (𝑈 − 𝑈𝑖 )
𝑑𝑡 𝑑𝑥 𝑖−1

12
 Itération

𝑑𝑈 1
𝑖 =1 = (𝑈0 − 𝑈1 )
𝑑𝑡 𝑑𝑥

𝑑𝑈 1
𝑖 =2 = (𝑈1 − 𝑈2 )
𝑑𝑡 𝑑𝑥

𝑑𝑈 1
𝑖 =3 = (𝑈2 − 𝑈3 )
𝑑𝑡 𝑑𝑥

⋮
𝑑𝑈 1
𝑖 = 𝑁𝑋  = (𝑈𝑁𝑋−1 − 𝑈𝑁𝑋 )
𝑑𝑡 𝑑𝑥

−1 0 0 0 … 𝑈1 𝑈0
1 −1 0 0 … 𝑈2 0
𝑑𝑈 1 0 1 −1 0 … 𝑈3 1
 = +
𝑑𝑡 𝑑𝑥 ⋮ ⋱ ⋱ … ⋮ 𝑑𝑥 ⋮
0
[0 … 0 1 − 1] [𝑈𝑁𝑋−1 ] [𝑈𝑁𝑋 ]

= M U + V

𝑑𝑈 1 1
 = 𝑀𝑈 + 𝑉
𝑑𝑡 𝑑𝑥 𝑑𝑥

𝑈(−10, 𝑡) = 𝑈(10, 𝑡) = sin(𝑡)

𝑈(𝑁𝑇) = (𝑠𝑖𝑛𝑡1 𝑠𝑖𝑛𝑡2 … 𝑠𝑖𝑛𝑡𝑁𝑇 )’

𝑈 0 = 𝑈(0) = 1

13
 Discrétisation en temps

0 < 𝑡 < T est le domaine temporel

𝑇−(0)
 𝑑𝑡 =
𝑁𝑇

𝑇
 𝑑𝑡 = 𝑝𝑎𝑠 𝑑𝑒 𝑙′𝑒𝑠𝑝𝑎𝑐𝑒
𝑁𝑇

𝑡𝑛 = 𝑛𝑑𝑡 , 𝑛 = 0, … , 𝑁𝑇 − 1

𝑈 𝑛 = 𝑈(𝑡𝑛 ) ;

𝑑𝑈 𝑈 𝑛+1 − 𝑈 𝑛
𝑈𝑡 = =
𝑑𝑡 𝑑𝑡

𝑑𝑈 1 1 𝑈 𝑛+1 − 𝑈 𝑛 1 1
= 𝑀𝑈 + 𝑉  = 𝑀𝑈 + 𝑉
𝑑𝑡 𝑑𝑥 𝑑𝑥 𝑑𝑡 𝑑𝑥 𝑑𝑥
𝑑𝑡 𝑑𝑡
 𝑈𝑛+1 = 𝑈𝑛 + 𝑀𝑈 𝑛 + 𝑉
𝑑𝑥 𝑑𝑥

Itération en fonction de n :

𝑑𝑡 𝑑𝑡
𝑛 = 0  𝑈1 = 𝑈 0 + 𝑀𝑈 0 + 𝑉
𝑑𝑥 𝑑𝑥
𝑑𝑡 𝑑𝑡
𝑛 = 1  𝑈 2 = 𝑈1 + 𝑀𝑈1 + 𝑉
𝑑𝑥 𝑑𝑥
𝑑𝑡 𝑑𝑡
𝑛 = 2  𝑈3 = 𝑈2 + 𝑀𝑈 2 + 𝑉
𝑑𝑥 𝑑𝑥

𝑑𝑡 𝑑𝑡
𝑛 = 𝑁𝑇 − 1  𝑈 𝑁𝑇 = 𝑈 𝑁𝑇−1 + 𝑀𝑈 𝑁𝑇−1 + 𝑉
𝑑𝑥 𝑑𝑥
14
 𝑝𝑜𝑠𝑜𝑛𝑠 :

𝑑𝑡
𝑁 = 𝐼𝑁 + 𝑀 𝑜𝑛 𝑎𝑢𝑟𝑎:
𝑑𝑥
𝑑𝑡
𝑈0 + 𝑉
𝐼𝑁 0 0 0 … 𝑈1 𝑑𝑥
−𝑁 𝐼𝑁 0 0 … 𝑈2 𝑑𝑡
𝑉
0 − 𝑁 𝐼𝑁 0 … 𝑈3 𝑑𝑥
=
⋮ ⋱ ⋱ … ⋮ ⋮
[0 … 0 − 𝑁 𝐼𝑁 ] [𝑈 𝑁𝑇−1 ] 𝑑𝑡
𝑁𝑇
[−𝑈 + 𝑉
𝑑𝑥 ]
A U = B
 𝐴𝑈 = 𝐵  𝑈 = 𝐴−1 𝐵
Cette solution sera implémentée par ODE45 sous MATLAB.
III. Discrétisation et Résolution par la méthode de ADAMS MOULTON
q=4 :
Pour les deux problèmes nous allons utiliser les discrétisations précédentes
ensuite utiliser le schéma ADAMS MOULTON q=4 e suivant :
Cette méthode est utilisée comme prédicteur, avec une méthode implicite
d’ordre 5 comme correcteur.
Pas Ordre ADAMS MOULTON q=4

q=4 p=5 𝒉
𝒚𝒌+𝟏 = 𝒚𝒌 + (𝟐𝟓𝟏𝒇𝒌+𝟏 + 𝟔𝟒𝟔𝒇𝒌 − 𝟐𝟔𝟒𝒇𝒌−𝟏 + 𝟏𝟎𝟔𝒇𝒌−𝟐 − 𝟏𝟗𝒇𝒌−𝟑 )
𝟕𝟐𝟎

15
A la suite de chacune des discrétisation des problèmes 1 et 2 , on ajoute la
boucle qui permet une résolution par la méthode de ADAMS MOULTON
q=4 ci-dessous :

Boucle
𝒉
𝒚𝒌+𝟏 = 𝒚𝒌 + (𝟐𝟓𝟏𝒇𝒌+𝟏 + 𝟔𝟒𝟔𝒇𝒌 − 𝟐𝟔𝟒𝒇𝒌−𝟏 + 𝟏𝟎𝟔𝒇𝒌−𝟐 − 𝟏𝟗𝒇𝒌−𝟑 )
𝟕𝟐𝟎

U0=U ;
U1=U ;
U2=U ;
U3=U ;
U4=U ;
for i=1:NT
U5=U3+(h) / (720) *(251*U4+646*U3-264*U2+106*U1-19*U0) ;
U0=U1 ;
U1=U2 ;
U2=U3 ;
U3=U4 ;
U4=U5 ;

end

U=U0 ;

U0, U1, U2 et U3 sont les quatre premiers termes du schéma ADAMS

MOULTON q=4 et U est la solution obtenue par la méthode implicite qui est
dans notre cas la méthode dz RUNGE KUTTA d’ordre 4 ou méthode des
différences finies implicites.

16
IV. LES ALGORITHMES DE DEUX METHODES NUMERIQUES
UTILISEES PERMETTANT DE RESOUDRE LES DEUX
PROBLEMES SOUS MATLAB :

Les algorithmes sont adaptés en fonction des particularités du problème à

résoudre, notamment en ce qui concerne la discrétisation spatiale, les
conditions aux limites et les schémas de différences finies utilisés.

III.1. Algorithme de la méthode des différences finies implicites :

Paramètres d'entrée :
- Condition initiale u0
- Discrétisation spatiale dx considérée comme pas de grille de l’espace.
- Discrétisation temporelle dt considérée comme pas du temps.
- Nombre d'itérations temporelles N
1. Initialisation :
- Calculer le nombre de points spatiaux M = (L / dx) + 1, où L est la
longueur du domaine.
- Initialiser la grille spatiale x avec les valeurs xi = i * dx pour i = 0, ..,
M-1.
- Initialiser le tableau u avec la condition initiale u0.
2. Calcul d’autres paramètres :
- Calculer le coefficient de diffusion α en fonction des propriétés du
problème.
3. Boucle temporelle :
Pour n allant de 1 à N:
a. Construire la matrice A pour les coefficients des différences finies
implicites.
b. Construire le vecteur b pour les termes connus.
c. Résoudre le système linéaire A * u= b pour obtenir u en utilisant une
méthode de résolution appropriée.
d. Mettre à jour u .
4. Sortie :
- Retourner le tableau u contenant la solution discrétisée.

17
III.2. Algorithme de la méthode de ADAMS MOULTON ordre q=4 :

Paramètres d'entrée :
- Condition initiale u0
- Discrétisation spatiale dx considérée comme pas de grille de l’espace.
- Discrétisation temporelle dt considérée comme pas du temps.
- Nombre d'itérations temporelles N
1. Initialisation :
- Calculer le nombre de points spatiaux M = (L / dx) + 1, où L, est la
longueur du domaine.
- Initialiser la grille spatiale x avec les valeurs xi = i * dx pour i = 0, ...,
M-1.
- Initialiser le tableau u avec la condition initiale u0.
2. Boucle temporelle :
Pour n allant de 0 à N :
a. Si n <= 3 :
- Appliquer une étape de la méthode des différences finies explicite
pour obtenir yk de la méthode de Runge-Kutta d’ordre 4.
b. Sinon :
- Construire les termes de prédiction pour la méthode de Adams
Moulton q=4.
- Appliquer la méthode de Adams Moulton q=4 pour obtenir u.
c. Mettre à jour u = u.
3. Sortie :
- Retourner le tableau u contenant la solution discrétisée.

V. LES FIGURES ISSUES DE LA RESOLUTION PAR LES DEUX

METHODES NUMERIQUES

18
IV.1. Les figures issues de la méthode des différences finies implicites :

IV.1.a. Courbe de l’équation de diffusion-réaction

IV.1.b. Courbe de l’équation de transport

19
IV.2. Les figures issues de la méthode de ADAMS MOULTON ordre
q=4
IV.1.a. Courbe de l’équation de diffusion-réaction

IV.1.b. Courbe de l’équation de transport

20
VI. COMPARAISON DE DEUX METHODES NUMERIQUES DE
RESOLUTION

Il est important de noter que l’évaluation ou la comparaison des méthodes

numériques n’est pas aussi simple parce que leurs exécutions peuvent
dépendre de la caractéristique du problème posé. Il également à noter qu’il
y a autres facteurs à considérer, comme la stabilité, polyvalence, preuve
contre l’erreur d’exécution…. La comparaison des méthodes des
différences finies implicites et de Adams Moulton d'ordre 4 par rapport à la
consistance, la stabilité, la convergence.

Méthode des différences finies implicites sont généralement plus

consistantes, stables, donc convergentes. Comparativement, la méthode
d’Adams Moulton d’ordre 4, toute la difficulté d’adoption du pas de
discrétisation réside dans le changement de longueur du pas. De surcroit, la
stabilité de la méthode de Adams Moulton est cependant assez médiocre.

21
VII. CONCLUSION

Tout au long de ce travail, nous avons eu à relater de façon générique la

résolution des problèmes d’équation de diffusion-réaction et celle de
transport par l’utilisation de deux méthodes numériques qui sont
respectivement la méthode de différences finies implicites et la méthode de
Adams Moulton d’ordre 4. Les algorithmes de ces deux méthodes relatives
aux problèmes ainsi que leur programmation sous Matlab ont été le centre
d’intérêt de notre projet. Ainsi, les figures découlant de ces méthodes ont été
déterminant pour mieux simuler numériquement les problèmes abordés. En
somme, le présent projet était criblé de difficultés à notre égard mais
également une ouverture enrichissante dans le cadre de notre formation aussi
bien présente que future.

22
VIII. Annexes
I. Les programmes MATLAB
I.1. Programme MATLAB méthode des différences finies

I.1.1. Programme ode45MDF équation de diffusion

▪ Pardef

% Paramètres par défaut de l’EDO d’ordre

global T NT NX dt dx epsilon

% T=1 ; NX=10 ; NT=200 ; dt=T/NT ; dx=2/NX;

epsilon=0.01;

▪ Matrice
function M=matrice
%Matrice de l’EDO d’ordre 2 u_t=εu_xx+10u en utilisant un schéma centré en espace
global T NT NX dt dx epsilon
alpha=2-(10*(dx)^2)/(epsilon);
V = [-alpha 1 zeros(1,NX-3)];
M = toeplitz(V);
end

▪ Solution
function [U,t,x]= solution
% résoud l’EDO d’ordre 2 ut=εuxx+10u en utilisant un schéma
% centré en espace
epsilon=0.01;
global T NT NX dt dx epsilon
t=0:dt :T ; x=-1:dx:1 ;
U0=0.53*x'+0.47*sin(-1.5*pi*x') ; U0=U0(2 :end-1) ; U=U0; % condition initiale
M=matrice;
for i=1:NT
U0=U0+epsilon*dt/(dx)^2*M*U0 ; U=[U U0] ; % stockage de la solution
end
U=[zeros(1,NT+1) ; U ; zeros(1,NT+1)]; % ajout des conditions aux limites

end

▪ Courbe
% trace les courbes l’EDO d’ordre 2
% en utilisant une discrétisation totale
global T NT NX dt dx epsilon
[U,t,x]= solution;
[tt,xx]= meshgrid(t,x) ;
subplot (1,1,1), mesh(tt,xx,U),
xlabel('temps');
ylabel('espace');
zlabel('quantite');
title('courbe de l''equation de diffusion(différence finie implicite)')
legend('solution numérique pour T=1 NX=100 NT=200')

▪ Simulation
function simulation(T0,NX0,NT0)
global T NT NX dt dx epsilon
pardef
T= T0 ; NX=NX0 ; NT=NT0 ; dt=T/NT ; dx =2/NX ;
Courbe

23
I.1.2. Programme ode45MDF équation de transport

▪ Pardef
% Paramètres par défaut de l’EDO d’ordre 1
global T NT NX dt dx
T=1 ; NX=100 ; NT=100 ; dt=T/NT ; dx=20/NX;

▪ Matrice
function M = matrice
global T NT NX dt dx
V = [-1 1 zeros(1, NX-3)];
M = toeplitz(V);
for i=1:length(V)-1
M(i,i+1)=0;
end
M;
end

▪ Solution
function [U,t,x]= solution
% résoud l’EDO d’ordre 2 en utilisant un schéma centré en espace
global T NT NX dt dx
t=0:dt:T ; x=-10:(20/NX):10 ;
U0=ones(NT+1,1) ; U0=U0(2 :end-1) ; U=U0; % condition initiale
M=matrice;
for i=1:NT
U0=U0+(dt/dx)*M*U0 ; U=[U U0] ; % stockage de la solution
end
U=[sin(t) ; U ; sin(t)]; % ajout des conditions aux limites

end

▪ Courbe
% trace les courbes l’EDO d’ordre Ut+Ux=0
% en utilisant une discrétisation totale
global T NT NX dt dx
[U,t,x]= solution;
[tt,xx]= meshgrid(t,x) ;
subplot (2,1,1), mesh(tt,xx,U),
xlabel('temps');
ylabel('espace');
zlabel('quantite');
title('courbe de l''equation de transport(différence finie implicite)')
legend('solution numérique T=1 NX= 100 et NT=100')

▪ Simulation

function simulation(T0,NX0,NT0)
global T NT NX dt dx
pardef
T= T0 ; NX=NX0; NT=NT0; NT=NX0; dt=T/NT ; dx =2/NX ;
Courbe

I.2. Programme MATLAB méthode de ADAMS MOULTON Q=4

I.2.1. Programme ode45 ADM4 équation de diffusion

24
▪ Pardef
% Paramètres par défaut de l’EDO d’ordre
global T NT NX dt dx epsilon
% T=1 ; NX=100 ; NT=200 ; dt=T/NT ; dx=2/NX;
epsilon=0.01;

▪ Matrice
function M=matrice
%Matrice de l’EDO d’ordre 2 u_t=εu_xx+10u en utilisant un %schéma DA
global T NT NX dt dx epsilon
alpha=2-(10*(dx)^2)/(epsilon);
V = [-alpha 1 zeros(1,NX-3)];
M = toeplitz(V);
end

▪ Solution
function [U,t,x]= solution
% résoud l’EDO d’ordre 2 en utilisant un schéma DA
epsilon=0.01;
global T NT NX dt dx epsilon
t=0:dt :T ; x=-1:dx:1 ;
U0=0.53*x'+0.47*sin(-1.5*pi*x') ; U0=U0(2 :end-1) ; U=U0; % condition initiale
M=matrice;
for i=1:NT
U0=U0+epsilon*dt/(dx)^2*M*U0 ; U=[U U0] ; % stockage de la solution
end

%boucle de la méthode de ADAMS MOULTON d'ordre 4

Boucle
𝒉
𝒚𝒌+𝟏 = 𝒚𝒌 + (𝟐𝟓𝟏𝒇𝒌+𝟏 + 𝟔𝟒𝟔𝒇𝒌 − 𝟐𝟔𝟒𝒇𝒌−𝟏 + 𝟏𝟎𝟔𝒇𝒌−𝟐 − 𝟏𝟗𝒇𝒌−𝟑 )
𝟕𝟐𝟎

U0=U ;
U1=U ;
U2=U ;
U3=U ;
U4=U ;
for i=1:NT
U5=U3+(h) / (720) *(251*U4+646*U3-264*U2+106*U1-19*U0) ;
U0=U1 ;
U1=U2 ;
U2=U3 ;
U3=U4 ;
U4=U5 ;

end
U=U0;
U=[zeros(1,NT+1) ; U ; zeros(1,NT+1)]; % ajout des conditions aux limites

end

▪ Courbe
% trace les courbes l’EDO d’ordre 2
% en utilisant une discrétisation totale
global T NT NX dt dx epsilon
[U,t,x]= solution;
[tt,xx]= meshgrid(t,x) ;
subplot (1,1,1), mesh(tt,xx,U),
xlabel('temps');
ylabel('espace');
zlabel('quantite');
title('courbe de l''equation de diffusion(méthode de ADAMS M q=4)')

25
 legend('solution numérique pour T=1 NX=100 NT=100')

▪ Simulation
function simulation(T0,NX0,NT0)
global T NT NX dt dx epsilon
pardef
T= T0 ; NX=NX0 ; NT=NT0 ; dt=T/NT ; dx =2/NX ;
courbe

I.2.2. Programme ode45 ADM équation de transport

▪ Pardef
% Paramètres par défaut de l’EDO d’ordre
global T NT NX dt dx
T=1 ; NX=100 ; NT=200 ; dt=T/NT ; dx=20/NX;

▪ Matrice
function M = matrice
global T NT NX dt dx
V = [-1 1 zeros(1, NX-3)];
M = toeplitz(V);
for i=1:length(V)-1
M(i,i+1)=0;
end
M;
end

▪ Solution
function [U,t,x]= solution
% résoud l’EDO d’ordre 1 ut+ux=0 en utilisant un schéma
% centré en espace
global T NT NX dt dx
t=0:dt:T ; x=-10:(20/NX):10 ;
U0=ones(NT+1,1) ; U0=U0(2 :end-1) ; U=U0; % condition initiale
M=matrice;
for i=1:NT
U0=U0+(dt/dx)*M*U0 ; U=[U U0] ; % stockage de la solution
end

%boucle de la méthode de ADAMS MOULTON d'ordre 4

Boucle
𝒉
𝒚𝒌+𝟏 = 𝒚𝒌 + (𝟐𝟓𝟏𝒇𝒌+𝟏 + 𝟔𝟒𝟔𝒇𝒌 − 𝟐𝟔𝟒𝒇𝒌−𝟏 + 𝟏𝟎𝟔𝒇𝒌−𝟐 − 𝟏𝟗𝒇𝒌−𝟑 )
𝟕𝟐𝟎

U0=U ;
U1=U ;
U2=U ;
U3=U ;
U4=U ;
for i=1:NT
U5=U3+(h) / (720) *(251*U4+646*U3-264*U2+106*U1-19*U0) ;
U0=U1 ;
U1=U2 ;
U2=U3 ;
U3=U4 ;
U4=U5 ;
end
U=U0;
U=[sin(t) ; U ; sin(t)]; % ajout des conditions aux limites

end

26
▪ Courbe
% trace les courbes l’EDO d’ordre 1 ut+ux=0
% en utilisant une discrétisation totale
global T NT NX dt dx
[U,t,x]= solution;
[tt,xx]= meshgrid(t,x) ;
subplot (1,1,1), mesh(tt,xx,U),
xlabel('temps');
ylabel('espace');
zlabel('quantite');
title('courbe de l''equation de transport(méthode de ADAMS MOULTON q=4)')
legend('solution numérique T=1 NX= 100 et NT=100')

▪ Simulation
function simulation(T0,NX0,NT0)
global T NT NX dt dx
pardef
T= T0 ; NX=NX0; NT=NT0; NT=NX0; dt=T/NT ; dx =2/NX ;
courbe

27