Modélisation et Calcul Scientifique 2 - MT 01
M1 SPMA/TI
Feuille de TD n◦ 4 : Equations différentielles. Problèmes aux limites
Exercice 1 : Problèmes aux limites
Dire si les équations différentielles suivantes sont linéaires ou non linéaires, et si elles sont avec
condition initiale ou avec condition aux bords :
1. y 00 + 4y 0 + 3y 2 = 0, y(0) = 1, y 0 (0) = 1,
2. y 00 + 2yy 0 = 0, y(0) = 1, y(1) = 0,
3. y 00 − y 0 + cos y = 0, y(0) = 1, y 0 (1) = 2,
4. y 00 + ty 0 + y = tet , y(0) = 1, y 0 (0) = −1.
Exercice 2 : Méthode des différences finies - Différentes conditions aux bords
On considère le problème vu en cours avec les conditions aux bords suivantes :
(1) Conditions aux bords de Dirichlet non homogènes :
−u00 (x) = f (x), x ∈ [0, 1], u(0) = a, u(1) = b,
(2) Conditions aux bords mixtes :
−u00 (x) = f (x), x ∈ [0, 1], u(0) = 0, u0 (1) = 0,
Dans chacun de ces deux cas :
1. Donner une solution exacte du problème considéré.
2. Ecrire le schéma aux différences finies pour approcher la solution de cette équation, en prenant
en compte les conditions aux bords ci-dessus.
3. La matrice du système linéaire obtenu est-elle inversible ?
Exercice 3 : Méthode des différences finies
On considère la généralisation suivante du problème vu en cours avec c(x) > 0 :
−u00 (x) + c(x)u(x) = f (x), x ∈ [0, 1], u(0) = 0, u(1) = 0.
1. Ecrire le schéma aux différences finies pour approcher la solution de cette équation.
2. La matrice du système linéaire obtenu est-elle inversible ?
3. Reprendre les preuves de la consistance, la stabilité, la convergence et l’ordre du schéma.
Exercice 4 : Méthode des différences finies - Conditions de Neumann homogènes
On considère le problème suivant :
−u00 (x) = f (x), u0 (0) = 0, u0 (1) = 0.
x ∈ [0, 1],
Z 1
1. Montrer que si ce problème admet une solution, alors f (t)dt = 0. On suppose que cette
0
condition est satisfaite, montrer que l’équation admet une solution à une constante près.
2. Ecrire le schéma aux différences finies pour approcher la solution de cette équation, en prenant
en compte les conditions aux bords ci-dessus.
3. La matrice du système linéaire obtenu est-elle inversible ? Donner un élément de son noyau.
Exercice 5 : Méthode des différences finies à 5 points
Montrer que la méthode des différences finies à 5 points définie par
−ui+2 + 16ui+1 − 30ui + 16ui−1 − ui−2
12h2
1
2
est d’ordre 4 pour approcher u00 .
Exercice 6 : Equation de la chaleur : introduction aux EDP
On cherche à trouver une valeur approchée de la fonction u(t, x) avec t > 0 et x ∈ [0, 1] vérifiant
l’équation de la chaleur suivante :
∂u ∂2u
− D 2 = 0.
∂t ∂x
On considère en espace les conditions aux bords
u(t, x = 0) = a, u(t, x = 1) = b, pour tout t > 0
et en temps la condition initiale
u(t = 0, x) = u0 (x), pour tout x ∈ [0, 1].
En temps, on considère le pas de temps ∆t > 0 et les temps discrets tn = n∆t, n ∈ N. En espace,
on considère le pas d’espace ∆x > 0 et les points discrets xi = i∆x, i ∈ {0, · · · , M + 1} tels que x0 = 0
et xM +1 = 1. On note uni une approximation de la valeur u(tn , xi ) de la fonction u évaluée au temps
tn et au point xi et on appelle U n le vecteur de taille M dont la i-ème composante est uni , 1 ≤ i ≤ M .
On utilise la méthode des différences finies pour discrétiser la dérivée en espace.
1. Utiliser la méthode d’Euler explicite pour la dérivée en temps et écrire le schéma ainsi obtenu
pour calculer le vecteur U n .
2. Comment calculer concrètement le vecteur U n ?
3. Utiliser la méthode d’Euler implicite pour la dérivée en temps et écrire le schéma ainsi obtenu
pour calculer le vecteur U n .
4. Comment calculer concrètement le vecteur U n ? La matrice du système linéiare ainsi obtenu
est-elle inversible ?
5. Estimer l’erreur de consistance de ces deux schémas. Quels sont leurs ordres en temps et en
espace ?
