REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROON

Paix-Travail-Patrie Peace Work Fatherland

********** **********
MINISTERE DE L’ENSEIGNEMENT SUPERIEUR MINISTRY OF HIGHER EDUCATION
********** **********
UNIVERSITE DE NGAOUNDERE THE UNIVERSITY OF NGAOUNDERE
********** **********
ECOLE DE GENIE CHIMIQUE ET DES SCHOOL OF CHEMICAL ENGINEERING AND MINERAL
INDUSTRIES MINERALES INDUSTRIES SCHOOL
********** **********
B.P./P.O. Box: 454
Ngaoundéré

NOTES DE COURS DE METHODES NUMERIQUES

Code - Libellé Méthodes numériques

Nombre de crédit 2

TPE: Utiliser EXCEL pour créer une feuille de calculs avec des formules adapter à un cahier de
charge présenté.

TP: Présenter la modélisation et la simulation numérique de l’évolution d’une densité de

produit chimique qui sera présenté.
 CHAPITRE 0

1- Introduction
a- Définition 1:
𝝏𝒇 𝝏𝒇
𝑭 (𝒙𝟏 𝒙𝒏 ) 𝟎 (1)
𝝏𝒙𝟏 𝝏𝒙𝒏

- Une EDP (équation aux dérivées partielles) est une équation ans laquelle figure une fonction
de plusieurs variables indépendantes et des dérivées partielles de par rapport à ses
variables.
- Une telle équation est dite d’ordre quand elle contient au moins une dérivée d’ordre sans en
contenir d’autres d’ordre supérieur.
- Toute fonction qui satisfait identiquement à cette équation (1) [dans un sens à
préciser] appelée solution de (1).

b- Motivation :

Les EDP (Equations aux Dérivées Partielles) les plus courantes de la science proviennent de la
modélisation de quelques phénomènes physiques (étendus en chimie, ingénierie, biologie,…) comme :

 Le transport d’entités (convection de chaleur dans un liquide, convection d’un polluant dans
l’atmosphère) ;

 La diffusion (comme celle de la chaleur dans un solide) ;

 Les vibrations (son dans l’air ou dans les structures) ;

 L’équilibre (calcul de l’équilibre d’une structure soumise à des forces).

Exemples 0:

 Equation des ondes :

 Equation de la chaleur : ( )

 Eq. de Poisson : (Une EDP du 2nd ordre linéaire, non homogène)

 Une EDP du 1er ordre non linéaire : ( )

 Une EDP du 1er ordre linéaire mais non homogène: ( )

 Une EDP du 2nd ordre, linéaire et homogène: .

Dans ce premier chapitre, on présentera quelques exemples de techniques de résolution d’EDP.

A- Méthode de la séparation des variables pour des EDP du 2nd ordre et coefficients de
Fourier
a- Cas fondamental de l’équation de la chaleur sur un segment réel en espace
On cherche à résoudre l’équation de la chaleur dans { avec les

relations de compatibilité entre la condition initiale et les conditions aux limites

1ère étape : La méthode de séparation des variables consiste à poser dans l’équation et de
“séparer les variables”. L’équation devient On divise alors formellement par
pour avoir . Comme le membre de gauche ne dépend que de t et le membre de droite que

de x on en déduit qu’ils sont constants, c’est à dire qu’il existe avec { et on obtient bien deux
équations à variables séparées.

2ème étape : On cherche les solutions non nulles de l’équation en avec les conditions aux limites, soit
{ . Les solutions dépendent de la constante :

√ √
- Si alors les solutions de l’équation différentielle sont : . Si on cherche
maintenant à tenir compte des conditions aux limites, il vient que A = B = 0. Il n’y a donc pas de
solutions non nulles dans ce cas.
- Si alors , alors les conditions aux limites donnent A = B = 0. Il n’y a donc pas
aussi de solutions non nulles dans ce cas non plus.
- Si alors les solutions de l’équation différentielle sont :
( √ ) ( √ ) . Les conditions aux limites induisent que les solutions non
nulles forment une famille de fonctions ( ) pour , associées aux valeurs de

suivantes . Au total, on a obtenu une suite infinie de solutions associées chacune à une
valeur de . On appelle les solutions , les fonctions propres du problème et les n les valeurs propres
associées. La fonction propre, comme les vecteurs propre en algèbre linéaire, sont définis à un scalaire
multiplicatif prés.

3ème étape : On remarque que le produit scalaire ∫ , orthogonalise la suite des

, dans le sens où , pour . Ceci indique que la suite des est une base
sur laquelle on va pouvoir développer la solution en la projetant grâce au produit scalaire.

4ème étape : On résout l’équation en pour les valeurs de n trouvées précédemment et sans se
préoccuper de la condition initiale. On a à résoudre l’équation qui a pour solutions

Etant donné qu’il n’y a pas de condition initiale à cette équation différentielle d’ordre 1 on trouve un espace
vectoriel de dimension 1 de solutions. est une constante arbitraire pour le moment. A ce stade, les
fonctions sont solutions de l’E.D.P. et des conditions aux limites mais pas de la
condition initiale.

5ème étape : On L’équation étant linéaire, la somme de plusieurs solutions à l’équation est toujours solution
de l’équation. On écrit donc la solution comme somme de toutes les solutions élémentaires
∑ . Il faut maintenant déterminer les coefficients cn pour que la solution u(t, x) vérifie la
condition initiale. Cette condition initiale s’écrit ; ce qui donne ∑
∑ .

Ici, les peuvent s’interpréter directement comme étant les coordonnées de la décomposition de dans la
base des . Comme les sont orthogonales pour le produit scalaire défini précédemment, on a
∫ ( ) . On obtient donc la solution sous la forme d’une série. Ici il s’agit

d’une série de Fourier où seules les composantes en ( ) sont présentes, puisqu’on a imposé : u(0) =
u(L)=0.

B- Méthode des semi-groupes (pratiques)

CHAPITRE 1
 CHAPITRE 2

RESOLUTION NUMERIQUE DES SYSTEMES LINEAIRES

 CHAPITRE 3

METHODES ITERATIVES GENERALES

 TRAVAUX DIRIGES

Bibliographie
1- Mazen SAAD, ANALYSE NUMERIQUE, Ecole Centrale de Nantes D´ept. Info/Math Ann´ee universitaire 2011-
2012.
2- Takéo Takahashi Polycopié rédigé par Michel Pierre et Antoine Henrot, Analyse Numérique, Ecole des Mines
de Nancy, Cours électif CE33 Semestre S7 : 2013-2014.