Université Sultan Moulay Slimane

Faculté Polydisciplinaire De Khouribga

MASTER MATÉRIAUX ET ÉNERGIES (MME)

Module : Méthodes Numériques

Rapport de TP : Méthodes Numériques pour les EDP

Réalisé par : Encadré par :

AFID Ikram Pr. Rachid Fakhar

Année universitaire 2025/2026

 Résumé
Ce rapport porte sur l’étude et l’implémentation de méthodes numériques pour la résolution
des équations différentielles et des équations aux dérivées partielles. Les méthodes des différences
finies et des éléments finis sont appliquées à des problèmes stationnaires et non stationnaires,
linéaires et non linéaires, en une et deux dimensions. Les résultats numériques, obtenus sous
MATLAB, mettent en évidence la stabilité, la convergence et la précision des schémas étudiés,
ainsi que l’efficacité des méthodes itératives pour les problèmes non linéaires.

-2-
 Abstract
This report focuses on the study and implementation of numerical methods for solving
ordinary and partial differential equations. Finite difference and finite element methods are
applied to stationary and time-dependent, linear and nonlinear problems in one and two
dimensions. Numerical simulations performed using MATLAB demonstrate the stability,
convergence, and accuracy of the proposed schemes, as well as the efficiency of iterative methods
for nonlinear problems.

-3-
Sommaire
Résumé...............................................................................................................................................2
Abstract ..............................................................................................................................................3
Introduction........................................................................................................................................6
Chapitre I: Résolution du problème stationnaire 1D avec conditions aux Limites de Dirichlet
homogène par MDF ........................................................................................................................8
I.1. Problème : ....................................................................................................................................8
I.2. Discrétisation du problème (méthode des différences finies) : ......................................................8
I.3. Approximation des opérateurs par développement de Taylor : .....................................................8
I.4. Schéma aux différences finies : .....................................................................................................8
I.5. Erreur :..........................................................................................................................................9
I.6. Résultats numériques sous MATLAB : .......................................................................................... 10
Chapitre II: Equation différentielle ordinaire d’ordre 2 à coefficients variables. ............................. 11
II.1. Problème : ................................................................................................................................. 11
II.2. Discrétisation du problème : ...................................................................................................... 11
II.3. Approximation des opérateurs par développement de Taylor : .................................................. 11
II.3.1. Approximation de la dérivée première :............................................................................... 12
II.3.2. Approximation de la dérivée seconde : ................................................................................ 12
II.4. Schéma aux différences finies :................................................................................................... 12
II.5. Résultats numériques sous MATLAB : ......................................................................................... 13
Chapitre III: Problème d'évolution 1D : .......................................................................................... 14
III.1. Problème : ................................................................................................................................ 14
III.2. Discrétisation par la méthode des différences finies : ................................................................ 14
III.3. Approximation de Taylor : ......................................................................................................... 15
III.4. Schémas aux différences finies : ................................................................................................ 15
III.4.1. Schéma explicite d’Euler : ................................................................................................... 15
III.4.2. Schéma implicite d’Euler :................................................................................................... 15
III.4.3. Schéma 𝜽 : ......................................................................................................................... 16
III.4.4. Erreur de consistance : ....................................................................................................... 16
III.5. Résultat numérique sous MATLAB : ........................................................................................... 16
Chapitre IV: Equation différentielle Non linéaire : ......................................................................... 18
IV.1. Problème : ................................................................................................................................ 18
IV.2. Résolution par la méthode de Newton : .................................................................................... 18
IV.3. Résultat numérique sous MATLAB :........................................................................................... 19
Chapitre V: MDF Bidimensionnelle : .............................................................................................. 20
V.1. Approximation par la méthode des différences finies : ............................................................... 20

-4-
 V.1.1. Équation modèle : ............................................................................................................... 20
V.1.2. Discrétisation spatiale : ....................................................................................................... 20
V.1.3. Approximation des dérivées : .............................................................................................. 21
V.1.4. Conditions aux limites : ....................................................................................................... 21
V.1.5. Forme matricielle : .............................................................................................................. 21
V.2. Problèmes stationnaires 2D : équation de Laplace : ................................................................... 22
V.3. Résultats numériques sous MATLAB :......................................................................................... 22
Chapitre VI: Modélisation par la MEF des équations de chaleur en 1D : ......................................... 24
VI.1. CAS 1 : Équation de la chaleur stationnaire (MEF 1D) : .............................................................. 24
VI.1.1. Problème : ......................................................................................................................... 24
VI.1.2. Discrétisation du problème : .............................................................................................. 24
VI.1.3. Approximation par la méthode des éléments finis : ............................................................ 24
VI.2. Équation de la chaleur non stationnaire (MEF 1D) : ................................................................... 25
VI.2.1. Problème : ......................................................................................................................... 25
VI.2.2. Discrétisation du problème : .............................................................................................. 25
VI.2.3. Approximation par éléments finis : ..................................................................................... 25
VI.3. Résultat numérique sous MATLAB :........................................................................................... 25
Chapitre VII: Problème Non Linéaire et Méthode des Eléments Finis : .......................................... 27
VII.1. Problème continu :................................................................................................................... 27
VII.2. Formulation variationnelle : ..................................................................................................... 27
VII.3. Approximation par la méthode des éléments finis : .................................................................. 27
VII.4. Écriture matricielle : ................................................................................................................. 27
VII.5. Résolution par la méthode de Newton–Raphson : .................................................................... 28
VII.6. Résultat Numérique sous MATLAB ........................................................................................... 28
Conclusion ........................................................................................................................................ 30

-5-
 Introduction

L’étude des méthodes numériques constitue un outil fondamental pour la résolution des
équations différentielles et des équations aux dérivées partielles rencontrées dans de
nombreux domaines de la physique et de l’ingénierie. En effet, la majorité des problèmes
réels issus de la diffusion thermique, de la mécanique des milieux continus, de
l’électromagnétisme ou encore des phénomènes de transport ne possèdent pas de solutions
analytiques explicites. Il devient alors indispensable de recourir à des méthodes numériques
afin d’obtenir des solutions approchées fiables et exploitables.

MATLAB (MATrix LABoratory) est un environnement de calcul scientifique largement

utilisé pour l’implémentation et l’analyse de ces méthodes numériques. Grâce à ses capacités
de manipulation matricielle, de résolution de systèmes linéaires et non linéaires, ainsi qu’à ses
outils de visualisation graphique, MATLAB constitue une plateforme privilégiée pour l’étude
numérique des équations différentielles ordinaires (EDO) et des équations aux dérivées
partielles (EDP).

Ce travail s’inscrit dans le cadre du module Méthodes Numériques et a pour objectif

principal d’étudier et de mettre en œuvre différentes techniques numériques pour la résolution
de problèmes stationnaires et non stationnaires, linéaires et non linéaires, en une et deux
dimensions. La démarche adoptée repose essentiellement sur la méthode des différences
finies (MDF) et la méthode des éléments finis (MEF), qui sont parmi les approches les plus
utilisées en calcul scientifique.

Dans un premier temps, nous étudions la résolution de problèmes stationnaires

unidimensionnels avec conditions aux limites de Dirichlet homogènes à l’aide de la méthode
des différences finies. Cette approche permet d’introduire les notions de discrétisation
spatiale, d’approximation des dérivées par développement de Taylor, de schémas aux
différences finies, ainsi que l’analyse des erreurs de consistance et de convergence.

Nous étendons ensuite cette méthode à des équations différentielles ordinaires d’ordre deux à
coefficients variables, ce qui permet d’illustrer l’adaptation des schémas numériques à des
opérateurs plus complexes. Par la suite, les problèmes d’évolution en une dimension,
notamment l’équation de la chaleur, sont étudiés en mettant en évidence le rôle du temps et en
comparant différents schémas temporels, tels que le schéma explicite d’Euler, le schéma
implicite d’Euler et le θ-schéma, avec une analyse de leur stabilité et de leur précision.

Une attention particulière est ensuite portée aux équations différentielles non linéaires, pour
lesquelles la résolution numérique nécessite l’utilisation de méthodes itératives, en particulier
la méthode de Newton–Raphson vectorielle, combinée à une discrétisation par différences
finies.

Le travail est ensuite élargi à la méthode des différences finies bidimensionnelle, appliquée
à des problèmes stationnaires en deux dimensions, notamment l’équation de Laplace. Cette
partie met en évidence la structure matricielle par blocs et l’impact des conditions aux limites
sur la formulation du système discret.

Enfin, les derniers chapitres sont consacrés à la méthode des éléments finis, appliquée à la
modélisation des équations de la chaleur en une dimension, aussi bien dans le cas stationnaire

-6-
que non stationnaire, ainsi qu’à la résolution de problèmes non linéaires par MEF couplée à la
méthode de Newton. Ces chapitres illustrent la puissance et la flexibilité de la MEF pour
traiter des problèmes plus généraux et physiquement réalistes.

L’ensemble des résultats numériques est implémenté et visualisé sous MATLAB, permettant
ainsi de valider les méthodes étudiées et de comparer les solutions numériques aux solutions
exactes lorsqu’elles sont disponibles.

-7-
Chapitre I: Résolution du problème stationnaire 1D avec conditions aux
Limites de Dirichlet homogène par MDF

I.1. Problème :

Considérons le problème aux limites suivant et on trouve u : [0,1] → ℝ tel que :

−𝑢′′ + 𝜋 2 𝑢(𝑥 ) = 2𝜋 2 𝑥𝑠𝑖𝑛(𝜋𝑥) − 2𝜋 cos(𝜋𝑥 ) 𝑑𝑎𝑛𝑠 ]𝑎, 𝑏[

{
𝑢 ( 0) = 0 , 𝑢 (1) = 0

La solution exacte de problème est : 𝒖𝒆𝒙 (𝒙) = 𝒙𝒔𝒊𝒏(𝝅𝒙)

I.2. Discrétisation du problème (méthode des différences finies) :

On commence par discrétiser l’intervalle [a,b] (ici [0,1]) en N+1 sous-intervalles de longueur :

𝑏−𝑎
ℎ=
𝑁+1
Les points du maillage sont : 𝒙𝟎 = 𝒂, 𝒙𝟏 … . . 𝒙𝒏 et 𝒙𝑵+𝟏 = 𝒃

I.3. Approximation des opérateurs par développement de Taylor :

On suppose que la solution u est suffisamment régulière (u ∈𝐶 4 ([𝑎, 𝑏]).

On utilise alors le développement de Taylor autour de x :

ℎ2 ℎ3 ℎ4
U(x+h)= u(x)+hu’(x)+ u’’(x)+ 𝑢3 (𝑥) + 𝑢4 (𝑥 + 𝜃ℎ)
2 6 24

ℎ2 ℎ3 ℎ4
U(x-h)= 𝑢(𝑥) − ℎ𝑢′ (𝑥) + 2
𝑢′′ (𝑥) − 6
𝑢3 (𝑥) + 24 𝑢4 (𝑥 − 𝜃ℎ)

En additionnant ces deux expressions :

𝑢(𝑥 + ℎ) − 2𝑢(𝑥) + 𝑢(𝑥 + ℎ)

𝑢′′ (𝑥) =
ℎ2

I.4. Schéma aux différences finies :

On remplace u’’ (𝑥𝑖 ) dans l’équation différentielle par son approximation discrète :
𝑢𝑖+1−2𝑢𝑖 +𝑢𝑖−1
- +𝛼𝑖 𝑢𝑖 =𝑓𝑖
ℎ2

Où :

o 𝛼𝑖 = 𝛼 (𝑥𝑖 ) = 𝜋 2
o 𝑓𝑖 = 𝑓 (𝑥𝑖 ) = 2𝜋 2 𝑥𝑖 sin(𝜋𝑥𝑖 ) − 2𝜋cos(𝜋𝑥𝑖 )

La forme matricielle : Au=B, où :

-8-
 𝛼1 0 𝑢
0 0 0 0 2 −1 0 0 0 0 𝑓(𝑥1 ) + 𝑎2
ℎ
⋮ 𝛼2 0 0 0 0 −1 2 ⋱ 0 0 0 𝑓(𝑥2 )
0 0 ⋱ 0 0 0
× 0 −1 ⋱ ⋱ ⋱ 0
1
A= + ℎ2 et B = ⋮
0 0 0 ⋱ 0 0 0 0 ⋱ ⋱ ⋱ 0
𝑓(𝑥𝑁−1 )
0 0 0 0 𝛼𝑁−1 0 0 0 0 ⋱ ⋱ −1 𝑢𝑏
[ 0 0 0 0 0 𝛼𝑁 ] [ 0 0 0 0 −1 2] (𝑓 (𝑥𝑁 ) + ℎ 2 )

I.5. Erreur :

On note 𝜋ℎ (𝑢𝑒𝑥 ) la projection de la solution exacte 𝑢𝑒𝑥 sur les points internes de la
subdivision spatiale. Cette projection est définie par :

𝑢𝑒𝑥 (𝑥1 )
𝜋ℎ (𝑢𝑒𝑥 ) = ( ⋮ ) ∈ 𝑅𝑁
𝑢𝑒𝑥 (𝑥𝑛 )

L’erreur de consistance est alors définie comme le vecteur de 𝑅𝑁 donné par :

𝑒ℎ (𝑢) = 𝐴𝜋ℎ (𝑢𝑒𝑥 ) − 𝐵

Un schéma numérique est dit consistant si l’erreur de consistance tend vers zéro lorsque le pas
de discrétisation h tend vers zéro, c’est-à-dire :

lim ∥ eh (u) ∥= 0
0

Plus précisément, il existe une constante C, indépendante de h, telle que :

∥ eh (u) ∥ ≤ Chq , où q désigne l’ordre de consistance du schéma.

Nous nous intéressons également à l’ordre de convergence, qui caractérise la vitesse à laquelle
la solution numérique 𝑢ℎ converge vers la solution exacte lorsque le pas h diminue. En
considérant deux pas de discrétisation ℎ1 et ℎ2 , et les erreurs correspondantes 𝑒1 𝑒𝑡 𝑒2 , on
obtient :
𝑞
𝑒1 ∥ eh1 (u) ∥ 𝐶ℎ1
= ≈ 𝑞
𝑒2 ∥ eh2 (u) ∥ 𝐶ℎ2

L’ordre de convergence peut alors être estimé par la relation suivante :

𝑒1
ln
𝑒2
𝑞≈ ℎ1
ln
ℎ2

Enfin, l’erreur de convergence est définie par le vecteur :

𝜀ℎ (𝑢) = 𝑢ℎ − 𝜋ℎ (𝑢𝑒𝑥 )

-9-
 qui mesure l’écart entre la solution numérique obtenue et la solution exacte projetée sur le
maillage.

I.6. Résultats numériques sous MATLAB :

 Résultat :

-10-
  Analyse : Les Figures (a) à (d) illustrent l’évolution de la solution numérique en
fonction du pas de discrétisation spatial. Pour des valeurs élevées de h,
l’approximation reste grossière, tandis que le raffinement progressif du maillage
1
permet une meilleure restitution de la solution exacte. Lorsque ℎ = la solution
32

numérique est pratiquement confondue avec la solution analytique, ce qui met en

évidence la convergence de la méthode des différences finies.

Chapitre II: Equation différentielle ordinaire d’ordre 2 à coefficients

variables.

II.1. Problème :

On considère le problème aux limites suivant : trouver une fonction u : [a,b] → ℝ

((1 + 𝑥 2 )𝑢′ (𝑥 ))′ = 2 + 6𝑥 2 + 2𝑥𝑐𝑜𝑠(𝑥 ) − (1 + 𝑥 2 )𝑠𝑖𝑛𝑥

{
𝑢(𝑎) = 1 , 𝑢(𝑏) = 2 + sin(1)

L’objectif est de déterminer numériquement la fonction u(x)u(x)u(x) satisfaisant cette équation

différentielle ainsi que les conditions aux limites associées.

II.2. Discrétisation du problème :

On subdivise l’intervalle [a,b] en N+1 sous-intervalles uniformes de longueur :

𝑏−𝑎
ℎ=
𝑁+1
On définit alors les points de la grille par : 𝒙𝟎 = 𝒂 , 𝒙𝑵+𝟏 = 𝒃 et 𝒙𝒊 = 𝒂 + 𝒊𝒉

II.3. Approximation des opérateurs par développement de Taylor :

-11-
II.3.1. Approximation de la dérivée première :

En supposant que la fonction u est suffisamment régulière, les développements de Taylor

autour du point x donnent :
ℎ2 ℎ3
U(x+h)= u(x)+hu’(x)+ u’’(x)+ 𝑢3 (𝑥) + 𝑂(ℎ3 )
2 6

ℎ2 ℎ3
U(x-h) = 𝑢(𝑥) − ℎ𝑢′ (𝑥) + 2
𝑢′′ (𝑥) − 6
𝑢3 (𝑥) + 𝑂(ℎ3 )

En combinant ces expressions, on obtient une approximation centrée de la dérivée première

d’ordre deux :
𝑢𝑖+1 − 𝑢𝑖−1
𝑢′ (𝑥𝑖 ) =
2ℎ

II.3.2. Approximation de la dérivée seconde :

De manière analogue, la dérivée seconde est approchée par la formule aux différences finies
centrées :

𝑢𝑖+1 − 2𝑢𝑖 + 𝑢𝑖−1

𝑢′′ (𝑥𝑖 ) =
ℎ2

II.4. Schéma aux différences finies :

L’équation différentielle peut s’écrire sous la forme équivalente :

((1 + 𝑥 2 )𝑢′ (𝑥))′ = 𝑓(𝑥)

En remplaçant les dérivées par leurs approximations discrètes aux points 𝑥𝑖 , on obtient le
schéma :

𝑢𝑖+1 − 2𝑢𝑖 + 𝑢𝑖−1 𝑢𝑖+1 − 𝑢𝑖−1

(1 + 𝑥 2 ) + 2𝑥 𝑖
ℎ2 2ℎ

En introduisant les notations :

𝝍𝒊 = 𝒙𝟐𝒊 ; 𝜶𝒊 = 𝟐𝒙𝒊

Formulation matricielle

Le système obtenu peut être exprimé sous la forme matricielle : A(x) u=B

Avec: A(x)=R+R(x)+S(x),

-12-
 2𝝍 𝟏 −𝝍𝟏 0 … 0
2 −1 0 … 0
−1 2 −1 ⋱ ⋮ −𝝍𝟐 2𝝍𝟐 − 𝝍𝟐 ⋱ ⋮
1 1
R= 2 0 −1 2 ⋱ 0 et R(x) = 2 0 −𝝍𝟑 2𝝍 𝟑 ⋱ 0
ℎ ℎ
⋮ ⋱ ⋱ ⋱ −1 ⋮ ⋱ ⋱ ⋱ −𝝍𝑵−𝟏
( 0 … 0 −1 2 ) ( 0 … 0 −𝝍𝑵 2𝝍𝑵 )

𝝍𝟏 𝛼1
𝑓 (𝑥1 ) − 𝑢𝑎 + 𝑢𝑎
ℎ2 2ℎ
𝑓(𝑥2 )
B= ⋮
𝑓(𝑥𝑁−1 )
𝝍𝑵 𝛼𝑁
(𝑓 (𝑥𝑁 ) − ℎ2 𝑢𝑏 − 2ℎ 𝑢𝑏 )

II.5. Résultats numériques sous MATLAB :

-13-
  Résultat :

 Analyse :

La solution numérique obtenue par la méthode des différences finies est régulière et
monotone sur l’intervalle étudié. L’absence d’oscillations et la continuité du profil confirment la
stabilité du schéma et la pertinence de la discrétisation spatiale.

Chapitre III: Problème d'évolution 1D :

III.1. Problème :

On considère l’équation de la chaleur unidimensionnelle définie sur l’intervalle spatial [0,1] et

sur l’intervalle de temps (0,T]. Le problème consiste à déterminer la fonction
u:[0,1] × [0,T]→R telle que :

𝝏𝒖 𝝏𝟐 𝒖
(𝒙, 𝒕) − 𝝂 𝟐 (𝒙, 𝒕) = 𝟎, 𝟎 < 𝑥 < 1, 0 < 𝑡 ≤ 𝑇,
𝝏𝒕 𝝏𝒙
𝒖(𝟎, 𝒕) = 𝒖(𝟏, 𝒕) = 𝟎, 𝟎 < 𝑡 ≤ 𝑇
{ 𝒖(𝒙, 𝟎) = 𝐬𝐢𝐧(𝝅𝒙) , 𝟎≤𝒙≤𝟏

La solution exacte associée à ce problème est connue et s’écrit :

2
𝑢𝑒𝑥 (𝑥, 𝑡) = 𝑒 −𝜋 𝑡 sin(𝜋𝑥)

III.2. Discrétisation par la méthode des différences finies :

La méthode des différences finies repose sur la discrétisation des variables spatiale et temporelle.

 L’intervalle temporel [0,T] est subdivisé en N sous-intervalles de pas

-14-
 𝑇
𝑘 = ∆𝑡 =
𝑁

 L’intervalle spatial [0,1] est subdivisé en M+1 sous-intervalles de pas

1
ℎ = ∆𝑥 =
𝑀+1

On définit alors les points de la grille :

𝑡𝑛 = 𝑛𝑘 , 𝑥𝑖 = 𝑖ℎ,

Et les approximations discrètes :

𝑢𝑖𝑛 ≈ 𝑢(𝑥𝑖 , 𝑡𝑛 ), 𝑓𝑖𝑛 = 𝑓(𝑥𝑖 , 𝑡𝑛 )

III.3. Approximation de Taylor :

 La dérivée temporelle est approchée par une différence finie d’ordre 1.

 La dérivée seconde spatiale est approchée par une différence centrée d’ordre 2 :

𝝏𝟐 𝒖 𝑢𝑛 − 2𝑢𝑛𝑖 + 𝑢𝑛𝑖−1
(𝑥𝑖 , 𝑡𝑛 ) ≈ 𝑖+1
ℎ2
𝝏𝒙 𝟐

III.4. Schémas aux différences finies :

III.4.1. Schéma explicite d’Euler :

Le schéma explicite d’Euler s’écrit :

𝑢𝑖𝑛+1 − 𝑢𝑖𝑛 𝑛
𝑢𝑖+1 − 2𝑢𝑖𝑛 + 𝑢𝑖−1
𝑛
−𝜈 2
= 𝑓𝑖𝑛
𝑘 ℎ
Sous forme matricielle, le système discret s’écrit :

1 1 𝜈
𝐼𝑀 𝑢𝑛+1 = 𝐵𝑛 + ( 𝐼𝑀 − 2 𝑅) 𝑢𝑛 ,
𝑘 𝑘 ℎ

Où :

 𝑰𝑴 : est la matrice identité de taille M×M.

 𝑹 ∶ est la matrice tri diagonale issue de la discrétisation spatiale.
 𝑩𝒏 : est le vecteur du second membre intégrant les conditions aux limites.

III.4.2. Schéma implicite d’Euler :

Le schéma implicite d’Euler est donné par :

-15-
 𝑢𝑖𝑛+1 − 𝑢𝑖𝑛 𝑛+1
𝑢𝑖+1 − 2𝑢𝑖𝑛+1 + 𝑢𝑖−1
𝑛+1
−𝜈 = 𝑓𝑖𝑛+1
𝑘 ℎ2
Le système matriciel correspondant est :

1 𝜈 1
( 𝐼𝑀 + 2 𝑅)𝑢𝑛+1 = 𝐵𝑛+1 + ( 𝐼𝑀 ) 𝑢𝑛 ,
𝑘 ℎ 𝑘

III.4.3. Schéma 𝜽 :

Le θ-schéma est une combinaison convexe des schémas explicite et implicite, paramétrée par
θ∈[0,1] :

𝑢𝑖𝑛+1 − 𝑢𝑖𝑛 𝑛
𝑢𝑖+1 − 2𝑢𝑖𝑛 + 𝑢𝑖−1
𝑛 𝑛+1
𝑢𝑖+1 − 2𝑢𝑖𝑛+1 + 𝑢𝑖−1
𝑛+1
− 𝜈 (1 − θ) − 𝜈θ = (1 − θ)𝑓𝑖𝑛 + θ𝑓𝑖𝑛+1
𝑘 ℎ2 ℎ2

 Pour 𝜃 = 0, le schéma est explicite

 Pour 𝜃 = 1, il est implicite
1
 Pour 𝜃 = , le schéma est Crank-Nicholson
2

Sous forme matricielle :

1 𝜈 1 𝜈
( 𝐼𝑀 + 𝜃 2 𝑅) 𝑢𝑛+1 − ( 𝐼𝑀 − (1 − 𝜃) 2 𝑅) 𝑢𝑛 = 𝜃𝐵𝑛+1 + (1 − 𝜃) 𝐵𝑛
𝑘 ℎ 𝑘 ℎ

III.4.4. Erreur de consistance :

L’erreur de consistance du θ-schéma est définie par :

1 𝜈 1 𝜈
𝑒ℎ𝑛 (𝑢) = (𝑘 𝐼𝑀 + 𝜃 2 𝑅) 𝜋ℎ𝑛+1 (𝑢𝑒𝑥 ) − ( 𝐼𝑀 − (1 − 𝜃) 2 𝑅) 𝜋ℎ𝑛 (𝑢𝑒𝑥 ) − 𝜃𝐵𝑛+1 -(1 − 𝜃) 𝐵𝑛
ℎ 𝑘 ℎ

Où 𝜋ℎ𝑛 (𝑢𝑒𝑥 ) désigne la projection de la solution exacte sur la grille spatiale au temps 𝑡𝑛 .

III.5. Résultat numérique sous MATLAB :

-16-
 1
 Tableau d’erreurs pour 𝜃 = , ∣𝑥𝑖 , 0.5) − 𝑢𝑖𝑁 ∣
4

-17-
Chapitre IV: Equation différentielle Non linéaire :

IV.1. Problème :

On considère le problème aux limites suivant :

−𝑢′′ (𝑥 ) + 𝑥𝑢(𝑥)3 = 𝑓 (𝑥 ), 0<𝑥<1

{ 𝑢 (0) = 0
𝑢 (1) = 1

Avec : f(x)=0

IV.2. Résolution par la méthode de Newton :

Après discrétisation par différences finies, on obtient :

𝐹1 (𝑢)
𝐹 (𝑢 ) = ( ⋮ )
𝐹𝑁 (𝑢)

Avec, pour i=1,…, N :

𝑢𝑖−1 − 2𝑢𝑖 + 𝑢𝑖+1

𝐹𝑖 (𝑢) = − + 𝑥𝑖 𝑢𝑖3 − 𝑓𝑖
ℎ2

En tenant compte des conditions aux limites 𝑢0 = 𝛼 et 𝑢𝑁+1 =𝛽, le vecteur s’écrit
explicitement :

-18-
 2𝑢1 − 𝑢2 𝛼
2
+ 𝑥1 𝑢13 − 2
ℎ ℎ
−𝑢1 + 2𝑢2 − 𝑢3
+ 𝑥2 𝑢23
ℎ2
𝐹 (𝑢 ) = ⋮
−𝑢𝑁−2 + 2𝑢𝑁−1 − 𝑢𝑁 3
+ +𝑥𝑁−1 𝑢𝑁−1
ℎ2
2𝑢𝑁 − 𝑢𝑁−1 3
𝛽
( + 𝑥 𝑁 𝑢 𝑁 − )
ℎ2 ℎ2

 La matrice Jacobienne est définie par :

𝜕𝐹𝑖
𝐽𝑖𝑗 =
𝜕𝑢𝑗

 Formule de Newton (vectorielle) :

À chaque itération k : J (𝑢𝑘 )𝛿𝑢𝑘 = −𝐹(𝑢𝑘 ) , 𝑢𝑘+1 = 𝑢𝑘 + 𝛿𝑢𝑘

Avec le critère d’arrêt :

∥𝑢𝑘+1 − 𝑢𝑘 ∥ < 𝑇𝑜𝑙é𝑟𝑎𝑛𝑐𝑒

IV.3. Résultat numérique sous MATLAB :

-19-
  Résultat :

Chapitre V: MDF Bidimensionnelle :

V.1. Approximation par la méthode des différences finies :

V.1.1. Équation modèle :

On considère le domaine bidimensionnel :

Ω = (𝑎𝑥 , 𝑏𝑥 ) × (𝑎𝑦 × 𝑏𝑦 )

De frontière Γ. On cherche une fonction :

𝑢∶ Ω →ℝ

Solution du problème aux limites suivant :

𝛿2 𝑢 𝛿2 𝑢
𝜈
𝛿𝑥 2
−𝜈 + 𝛼(𝑥, 𝑦)𝑢(𝑥, 𝑦) = 𝑓(𝑥, 𝑦) dans Ω
𝛿𝑦 2

Avec la condition aux bords :

u=g sur Γ

V.1.2. Discrétisation spatiale :

Le domaine est discrétisé par une grille régulière avec des pas :

𝑏𝑥−𝑎𝑥 𝑏𝑦 −𝑎𝑦
𝛿𝑥 = , 𝛿𝑦 =
𝑁𝑥 +1 𝑁𝑦 +1

-20-
 Les points de la grille sont notés :

𝑥𝑖 = 𝑎𝑥 + 𝑖𝛿𝑥 , 𝑦𝑗 = 𝑎𝑦 + 𝑗𝛿𝑦

Et l’approximation de la solution s’écrit :

𝑢(𝑥𝑖 𝑦𝑗 ) ≈ 𝑢𝑖𝑗 , i=1,……………,𝑁𝑥 , j=1,………., 𝑁𝑦

V.1.3. Approximation des dérivées :

Les dérivées secondes sont approchées par des différences finies centrées :

𝛿 2 𝑢 𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗

≈
𝛿𝑥 2 (𝛿𝑥 )2

𝛿 2 𝑢 𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1

≈
𝛿𝑦 2 (𝛿𝑦)2

On obtient alors le schéma aux différences finies :

𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗 𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1

−𝜈 −𝜈 + 𝛼𝑖,𝑗 𝑢𝑖,𝑗 = 𝑓𝑖,𝑗
(𝛿𝑥) 2 (𝛿𝑦)2

En posant :

𝜈 𝜈
𝛾𝑥 = (𝛿𝑥)2 , 𝛾𝑦 = (𝛿𝑦)2 , 𝛽 = 2(𝛾𝑥 + 𝛾𝑦 )

L’équation discrète devient :

−𝛾𝑥 𝑢𝑖−1,𝑗 − 𝛾𝑥 𝑢𝑖+1,𝑗 − 𝛾𝑦 𝑢𝑖,𝑗−1 , 𝛾𝑦 𝑢𝑖,𝑗+1 + (𝛽 + 𝛼𝑖,𝑗 )𝑢𝑖,𝑗 = 𝑓𝑖,𝑗

V.1.4. Conditions aux limites :

Les conditions aux bords sont imposées par :

𝑢0,𝑗 = 𝑔0,𝑗 , 𝑢𝑁𝑥 +1,𝑗 = 𝑔𝑁𝑥 +1,𝑗

𝑢𝑖,0 = 𝑔𝑖,0 , 𝑢𝑖,𝑁𝑦 +1 = 𝑔𝑖,𝑁𝑦+1

Les valeurs de bord apparaissent dans le second membre du système.

V.1.5. Forme matricielle :

Les inconnues 𝑢𝑖,𝑗 sont réorganisées dans un vecteur :

𝑢 = (𝑢1,1 , … … , 𝑢𝑁𝑥 ,1 , 𝑢1,2 , … . , 𝑢𝑁𝑥 𝑁𝑦 )𝑇

Le système discret s’écrit alors : Au=B

-21-
 Où A est une matrice tridiagonale par blocs :

𝐷 −𝐹 0 … 0
−𝐹 𝐷 −𝐹 … 0
𝐴 = 𝐴0 + 0 −𝐹 𝐷 ⋱ ⋮
⋮ ⋮ ⋱ 𝐷 −𝐹
( 0 0 … −𝐹 𝐷 )

D est donnée par :

𝛽 −𝛾𝑥 0 … 0
−𝛾𝑥 𝛽 −𝛾𝑥 … 0
𝐷= 0 −𝛾𝑥 𝛽 ⋱ ⋮
⋮ ⋮ ⋱ 𝛽 −𝛾𝑥
( 0 0 … −𝛾𝑥 𝛽 )

V.2. Problèmes stationnaires 2D : équation de Laplace :

Dans le domaine [0,1]×[0,1] on considère :

𝛿2 𝑢 𝛿2 𝑢
-
𝛿𝑥 2
− = 𝑓(𝑥, 𝑦), 𝑢 = 0 𝑠𝑢𝑟 𝛿Ω
𝛿𝑦 2

Pour un maillage uniforme h=δx=δy, le schéma à cinq points donne :

1
(−𝑢𝑖,𝑗−1 − 𝑢𝑖−1,𝑗 + 4𝑢𝑖,𝑗 − 𝑢𝑖+1,𝑗 − 𝑢𝑖,𝑗+1 ) = 𝑓𝑖,𝑗
ℎ2
En numérotant les inconnues par :

𝑘 = 𝑖 + (𝑗 − 1)𝑁𝑥

Le système s’écrit :

1
(−𝑢𝑘−𝑁𝑥 − 𝑢𝑘−1 + 4𝑢𝑘 − 𝑢𝑘+1 − 𝑢𝑘+𝑁𝑥 ) = 𝑓𝑘
ℎ2
La matrice associée possède une structure par blocs :

𝐷 𝐼 0 −4 1 0
(
𝐴= 𝐼 𝐷 𝐼) Avec : (
𝐷 = 1 −4 1)
0 𝐼 𝐷 0 1 −4

V.3. Résultats numériques sous MATLAB :

-22-
 Résultat :

 Les erreurs obtenues :

-23-
Chapitre VI: Modélisation par la MEF des équations de chaleur en 1D :

VI.1. CAS 1 : Équation de la chaleur stationnaire (MEF 1D) :

VI.1.1. Problème :

On considère l’équation de la chaleur stationnaire en une dimension spatiale, définie sur

l’intervalle [0,1] :

−𝑢′′ (𝑥 ) = 𝜋 2 𝑥𝑠𝑖𝑛(𝜋𝑥) − 2𝜋 cos(𝜋𝑥 ) , x ∈ [0,1]

Avec les conditions aux limites de Dirichlet :

𝑢(0) = 0, 𝑢(1) = 0

La solution exacte de ce problème est connue et donnée par :

𝑢(𝑥 ) = 𝑥𝑠𝑖𝑛(𝜋𝑥)

Ce problème est un cas particulier de l’équation de Poisson unidimensionnelle :

−𝑢′′ (𝑥 ) = 𝑓(𝑥)

VI.1.2. Discrétisation du problème :

L’intervalle [0,1] est discrétisé en N sous-intervalles de longueur uniforme h.

On introduit un espace d’éléments finis 𝑉ℎ constitué de fonctions de forme linéaires par morceaux
définies sur ce maillage.

La formulation variationnelle discrète consiste à chercher une solution approchée

𝑢ℎ ∈ 𝑉ℎ telle que :

𝑎(𝑢ℎ , 𝜈) = 𝐿(𝜈), ∀ 𝜈 ∈ 𝑉ℎ
1
Où : 𝑎(𝑢ℎ , 𝜈) = ∫0 𝑢′ ℎ (𝑥 )𝜈 ′ (𝑥)𝑑𝑥

1
𝐿(𝜈) = ∫ 𝜋 2 𝑥𝑠𝑖𝑛(𝜋𝑥 ) − 2𝜋 cos(𝜋𝑥) 𝜈 (𝑥 )𝑑𝑥
0

VI.1.3. Approximation par la méthode des éléments finis :

L’approximation par éléments finis conduit à un système linéaire matriciel de la forme :

𝑅𝑢 = 𝐵

Où :

 R est la matrice de rigidité, B est le vecteur du second membre.

-24-
 Les coefficients sont définis par :
1
𝑅𝑖𝑗 = ∫ 𝜑𝑖′ (𝑥 )𝜑𝑗′ (𝑥 )𝑑𝑥 ,
0

1
𝐵𝑖 = ∫ 𝑓(𝑥 )𝜑𝑖 (𝑥 )𝑑𝑥
0

Avec 𝜑𝑖 (𝑥 ) les fonctions de forme

VI.2. Équation de la chaleur non stationnaire (MEF 1D) :

VI.2.1. Problème :

On considère le problème d’évolution suivant :

𝜕𝑢 1 𝜕2𝑢
(𝑥, 𝑡) − (𝑥, 𝑡) = 0, x ∈ (0,1), 𝑡 ∈ (0, 𝑇)
𝜕𝑡 𝜋2 𝜕𝑥 2

Avec les conditions aux limites :

u(0,t)=u(1,t)=0

Et la condition initiale :

1
𝑢(𝑥, 0) = cos(𝜋 (𝑥 − ))
2

La solution exacte est donnée par :

1
𝑢𝑒𝑥 (𝑥, 𝑡) = 𝑒 −𝑡 cos(𝜋 (𝑥 − ))
2
VI.2.2. Discrétisation du problème :

L’intervalle spatial [0,1] est discrétisé par la méthode des éléments finis en N éléments.
L’intervalle temporel [0,T] est subdivisé en pas de temps de longueur Δt

Le schéma temporel utilisé est le schéma d’Euler implicite, assurant la stabilité du calcul.

VI.2.3. Approximation par éléments finis :

La formulation variationnelle discrète appliquée au problème non stationnaire conduit au

système matriciel suivant :

1
(𝑀 + 𝑅)𝑢𝑘+1 = 𝑀𝑢𝑘
Δt

VI.3. Résultat numérique sous MATLAB :

-25-
 Résultat :

-26-
Chapitre VII: Problème Non Linéaire et Méthode des Eléments Finis :

VII.1. Problème continu :

On cherche une fonction :

𝑢 ∶ [0,1] → ℝ

solution du problème non linéaire suivant :

−𝑢"(𝑥) + 𝑢3 (𝑥) + 3𝑢(𝑥) = 𝑓(𝑥), x ∈ ]0,1[

avec les conditions aux limites :

u(0)=u(1)=0

Ce problème est non linéaire à cause du terme 𝑢3

Il n’admet pas de solution analytique connue

VII.2. Formulation variationnelle :

On multiplie l’équation par une fonction test 𝜈 ∈ 𝐻01 et on intègre sur [0,1].

Après intégration par parties et en tenant compte des conditions aux limites, on obtient :

Trouver u ∈ 𝐻01 tel que :

1 1
∫ (𝑢′ (𝑥 )𝜈 ′ (𝑥) + 𝑢3 (𝑥 )𝜈(𝑥 ) + 3𝑢(𝑥 )𝜈(𝑥 ))𝑑𝑥 = ∫ 𝑓 (𝑥 )𝜈(𝑥 )𝑑𝑥
0 0

VII.3. Approximation par la méthode des éléments finis :

On discrétise l’intervalle [0,1] en N éléments finis et on introduit les fonctions de base {𝜑𝑖 }𝑁
𝑖=1

Les fonctions approchées s’écrivent :

𝑁 𝑁

𝑢 ℎ ( 𝑥 ) = ∑ 𝑢 𝑖 𝜑𝑖 ( 𝑥 ) , 𝜈ℎ (𝑥) = ∑ 𝜈𝑖 𝜑𝑖 (𝑥 )
𝑖=1 𝑖=1

En utilisant la méthode de Galerkin, on obtient le système discret :

1 1
∑𝑁 𝑁 3 𝑁 1
𝑖=1 𝑢𝑖 ∫0 𝜑′𝑖 (𝑥)𝜑′𝑗 (𝑥) 𝑑𝑥 + ∑𝑖=1 𝑢𝑖 𝜑𝑖 (𝑥)𝜑𝑗 (𝑥)𝑑𝑥+3 ∑𝑖=1 𝑢𝑖 𝜑𝑖 (𝑥)𝜑𝑗 (𝑥)𝑑𝑥 = ∫0 𝑓(𝑥)𝜑𝑗 (𝑥)𝑑𝑥

VII.4. Écriture matricielle :

Le système discret peut s’écrire sous la forme matricielle non linéaire :

(𝑅 + 3𝑀 )𝑢 + 𝑀𝑢3 = 𝐵

-27-
Où :

 R est la matrice de rigidité

 M est la matrice de masse
 B est le vecteur du second membre

On définit alors :

𝐹(𝑢) = (𝑅 + 3𝑀)𝑢 + 𝑀𝑢3 − 𝐵

VII.5. Résolution par la méthode de Newton–Raphson :

Le système non linéaire F(u)=0 est résolu par la méthode de Newton.

À chaque itération k, on calcule :

𝑢𝑘+1 = 𝑢𝑘 + 𝛿𝑢 avec 𝛿𝑢 = −𝐽(𝑢𝑘 )−1 𝐹(𝑢𝑘 )

La matrice jacobienne est donnée par :

(u) = R + 3M + K(u)

Avec :

Et :

VII.6. Résultat Numérique sous MATLAB

-28-
 Résultat

-29-
 Conclusion

À travers ce travail, nous avons étudié et implémenté différentes méthodes numériques pour
la résolution des équations différentielles et des équations aux dérivées partielles, en mettant
l’accent sur la méthode des différences finies et la méthode des éléments finis. Les résultats
obtenus montrent que ces méthodes permettent d’approcher efficacement les solutions exactes
lorsque celles-ci existent, et offrent des solutions numériques stables et convergentes dans le
cas général.

L’analyse des erreurs de consistance et de convergence a permis de mieux comprendre

l’influence du pas de discrétisation spatiale et temporelle sur la précision des solutions
numériques. L’étude comparative des schémas explicites, implicites et du θ-schéma pour les
problèmes d’évolution a également mis en évidence l’importance du choix du schéma
temporel en termes de stabilité.

Les problèmes non linéaires, résolus par la méthode de Newton–Raphson, ont montré la
nécessité d’utiliser des méthodes itératives adaptées et une formulation matricielle rigoureuse,
notamment dans le cadre de la méthode des éléments finis. De plus, l’extension aux
problèmes bidimensionnels a illustré la complexité croissante des systèmes numériques et
l’intérêt des structures matricielles par blocs.

En conclusion, ce travail met en évidence la complémentarité des méthodes numériques

étudiées et leur rôle essentiel dans la modélisation et la simulation de phénomènes physiques
complexes. L’utilisation de MATLAB a permis de valider efficacement ces approches et
constitue un outil indispensable pour l’apprentissage et l’application des méthodes
numériques en sciences et en ingénierie.

-30-
-31-