Université de Annaba–“Departement de Genie de Procedes”
Master 1-Genie Chimique 2019–2020
Analyse Numerique
Les reponses doivent etre envoyées a l’adresse email: [Link]@[Link]
Remarques Importantes:
• Les reponses doivent etre envoyées au plus tard Mercredi 23 Sept à 9h matin a l’adresse email suivante: abdallah.
bradji@[Link]
• Aucune reponse ne sera acceptée apres Mercredi 23 Sept à 9h matin.
• Si vous avez des questions, vous m’ecrivez via l’adresse email [Link]@[Link]
• La presentation de la réponse est comptée sur 3pts
Exercice 1. Questions de cours.
1. Pour resoudre un probleme en Maths Appliquées, on suit quelques etapes en commençant de Phenomene Physique.
Citer au moins trois etapes.
2. Quels sont les principes de la methode de differences finies.
Exercice 2. Questionnaire à choix multiples (QCM).
1. L’equation ut − uxx = f represente l’un des phynomenes physiques (Choisir la bonne reponse):
1-Diffusion de Chaleur
2-Vibration d’Onde
2. Phynomene de Vibration peut etre modelisé par l’une des equations suivantes (Choisir la bonne reponse):
√
1-iut − uxx = f avec = −1.
2-utt − uxx = f
Exercice 3. Construire des approximations de Differences Finies pour les problemes suivants:
1.
−u00 (x) = f (x), x ∈ (0, 1) (1)
avec u(0) = 1 et u(1) = −1.
2.
−u00 (x) + u(x) = f (x), x ∈ (0, 1) (2)
avec u(0) = u(1) = 0.
Corrigé Examen
Barème
1. Exercice 1 : 3pt +3pt.
2. Exercice 2 : 1,5pt +1,5pt
3. Exercice 3 : 4pt+4pt
4. Presentation de copie: 3points
Solution Exercice 1.
1. Les etapes a suivre pour resoudre un probleme en Mathematiques Appliquees en commencant par le phynomene
physique : Phynomene Physique, Modelisation, L’etude Mathematique, L’approximation Numerique, Programma-
tion, Validation des Resultats
2. Les principes de la methode de Differences Finies : Subdivision de domaine par des points et l’approximation des
derivees dans ces points de maillage.
Solution Exercice 2.
1. La reponse juste est le choix 1 : L’equation ut − uxx = f respresente le phynomene de diffusion de chaleur.
2. La reponse juste est le choix 2 : Phynomene de vibration peut etre modelisé par l’equation utt − uxx = f .
Solution Exercice 3.
Probleme 1. I will consider here the uniform mesh h = 1/(M + 1), with M ∈ N? . The mesh points are denoted by,
for i = 0, . . . , M + 1
xi = ih (3)
Replacing now x by xi in (1) to get
−uxx (xi ) = f (xi ). (4)
The term uxx (xi ) can be approximated using the three point central scheme
u(xi+1 ) − 2u(xi ) + u(xi−1 )
. (5)
h2
We will have then the scheme:
un n n
i+1 − 2ui + ui−1
− = f (xi ), i = 1, . . . , M (6)
h2
with
u0 = 1 and uM +1 = −1.
The scheme (6) leads to the following linear system:
1
AU = F, (7)
h2
where
U = (u1 , u2 , . . . , uM )T , (8)
2 −1 0 ... 0
−1 2 −1 ... 0
A= . (9)
. ..
. .
0 ... 0 −1 2
and
1 1
F = (f (x1 ) + , f (x2 ), f (x3 ), . . . , f (xM −1 ), f (xM ) − 2 )T . (10)
h2 h
Probleme 2. I will consider here the uniform mesh h = 1/(M + 1), with M ∈ N? . The mesh points are denoted by,
for i = 0, . . . , M + 1
xi = ih (11)
Replacing now x by xi in (2) to get
−u00 (xi ) + u(xi ) = f (xi ). (12)
The term uxx (xi ) can be approximated using the three point central scheme (5).
We will have then the scheme:
un n n
i+1 − 2ui + ui−1
− + ui = f (xi ), i = 1, . . . , M (13)
h2
with u0 = u1 and uM +1 = 0.
The scheme (13) leads to the following linear system:
1
( A + I)U = F, (14)
h2
where
U = (u1 , u2 , . . . , uM )T , (15)
2 −1 0 ... 0
−1 2 −1 ... 0
A= . (16)
. ..
. .
0 ... 0 −1 2
I is the identity matrix
1 0 0 ... 0
0 1 0 ... 0
0 0 1 ... 0
I= (17)
.. ..
. .
0 ... 0 0 1
and
F = (f (x1 ), f (x2 ), f (x3 ), . . . , f (xM ))T . (18)