République de Côte
d’Ivoire
Union-Discipline-Travail
UNIVERSITE
JEAN LOROUGNON GUEDE
UFR Environnement
TD D’ANALYSE NUMERIQUE
GENIE DE L’EAU ET DE L’ENVIRONNEMENT
ECUE 2 : RESOLUTION NUMERIQUE
DES EDP
Exercice 1
1) Trouver 𝛼, 𝛽 et 𝛾 tels que :
𝛼𝑓 (𝑥 + 2∆𝑥) − 𝛽𝑓 (𝑥 − ∆𝑥) + 𝛾𝑓(𝑥)
𝑓 ′ (𝑥 ) ≈
∆𝑥²
2) De quel ordre est ce schéma ?
Exercice 2 :
Soit l’équation aux dérivées partielles suivantes :
𝜕𝑢 𝜕2𝑢
(𝑥, 𝑡 ) − 2 (𝑥, 𝑡 ) = 0 pour (𝑥, 𝑡) ∈ ]0; 𝐿[ × ]0; 𝑇[
𝜕𝑡 𝜕𝑥
𝑢(0, 𝑡 ) = 𝑢(𝐿, 𝑡 ) = 0
{ 𝑢(𝑥, 0) = 𝑓(𝑥) pour 𝑥 ∈ ]0; 𝐿[
On notera ∆𝑡 le pas de temps, ∆𝑥 le pas d’espace, N le nombre de nœuds suivant
l’espace x et M le nombre de nœuds suivant le temps T.
1) a) En utilisant un schéma explicite, donner le problème discret du problème
continue ci – avant.
b) Ecrire un script scilab permettant d’approcher la valeur de la solution exacte aux
points de discrétisation à l’aide du schéma explicite.
3) a) En utilisant un schéma implicite, donner le problème discret du problème
continue ci – avant.
b) Pour N = 6 et M = 6, écrire le système matriciel associé au problème discret ci –
dessus.
Page 1 sur 3
Exercice 3 :
Résoudre l’équation de Laplace suivante
𝑢′′ = 1 − 𝑥 pour 𝑥 ∈ ]0,1[
{ 𝑢 (0) = 0
𝑢′(1) = 𝑥²
On prendra comme pas ∆𝑥 = 1/4.
N.B : On prendra soin de faire figurer sur sa copie, les différentes étapes de la
résolution et de préciser les méthodes utilisées.
Exercice 4:
Soit l’équation aux dérivées partielles suivantes :
𝜕𝑈 𝜕²𝑈
+𝐷 = 0 pour (𝑥, 𝑡) ∈ [−1,1] × [0; 𝑇]
(𝐼) 𝜕𝑡 𝜕𝑥²
𝑈(−1, 𝑡 ) = 𝑈𝐺 ; 𝑈(1, 𝑡 ) = 𝑈𝐷 ∀𝑡 ∈ [0, 𝑇]
{ 𝑈(𝑥, 0) = 𝑈0 (𝑥) ∀𝑥 ∈ [−1,1]
On notera ∆𝑡 le pas de temps, ∆𝑥 le pas d’espace, N le nombre de nœuds suivant
l’espace x, M le nombre de nœuds suivant le temps t.
1) En utilisant un schéma implicite en temps et un schéma avant en espace, discrétisez
l’équation (𝐼).
1
2) Pour ∆𝑥 = , écrire le système associé à l’équation discrétisée ci – avant.
4
3) Mettre ce système sous forme matricielle 𝐴𝑋 𝑛+1 = 𝑏 ; On précisera 𝐴, 𝑋 et 𝑏.
Exercice 5:
Soit l’équation (𝐸) suivante
𝜕𝑢 𝜕𝑢
𝑦 +𝑥 = 0 pour (𝑥, 𝑦) ∈ ]0,1[ × ]0,1[
𝜕𝑥 𝜕𝑦
(𝐸) ∶
𝑢(0, 𝑦) = 𝑢(1, 𝑦) = 𝑦 ∀𝑦 ∈ ]0,1[
{ 𝑢(𝑥, 0) = 𝑢(𝑥, 1) = 𝑥 ∀𝑥 ∈ ]0,1[
On considère également le quadrillage associé de pas ∆𝑦 et ∆𝑥 ci – dessous :
y
𝑢13 𝑢23 𝑢33
𝑦3
𝑢12 𝑢22 𝑢32
𝑦2
𝑦1 𝑢11 𝑢21 𝑢31
Page 2 sur 3
𝑥0 , 𝑦0
𝑥1 𝑥2 𝑥3 x
On prendra comme pas ∆𝑥 = ∆𝑦 = 1/3.
1. Quel type de condition avons-nous ? (préciser sa nature).
2. Combien y a t – il de nœuds inconnus ?
3. En utilisant la méthode des différences finies et un schéma décentré avant,
donner l’équation discrète associée à l’équation (𝐸).
4. Ecrire le système 𝐴𝑋 = 𝑏, associé à l’équation discrétisée ci – avant
Page 3 sur 3