Ecole Nationale d’Ingénieurs de Tunis
Unité Pédagogique de Mathématiques
Devoir Surveillé
Classes : 1ère année GC–GHE–GI–GM–MINDS Date : 7 Mars 2019
Matière : Mathématiques de l’Ingénieur II Documents : non autorisés
Enseignants : [Link], H. Chaker, N. Chouaieb, Y. Daoud, [Link] Durée : 2h
Exercice 1
On considère l’équation d’advection linéaire :
∂u ∂u
∗
+c = 0, (x, t) ∈ R × IR+
∂t ∂x (1)
x∈R
u(x, 0) = u0 (x),
où c > 0 et u0 est une fonction C 4 (R) donnée. On suppose dans toute la suite que l’équation (1) admet une
solution unique u(x, t) de classe C 4 (R × R+ ).
Pour approcher la solution u, on propose le schéma suivant : Pour (β, θ) ∈ [0, 1]2 ,
n+1
vj+1 n
− vj+1 vjn+1 − vjn n+1
vj+1 − vjn+1 n
vj+1 − vjn
β + (1 − β) + c(θ + (1 − θ) )=0 (2)
∆t ∆t ∆x ∆x
où vjn désigne une approximation de u(xj , tn ), avec xj = j∆x et tn = n∆t, où j ∈ Z et n ≥ 0. On suppose que
∆t
= λ, avec λ ∈ R∗+ fixé.
∆x
1. Soit εnj l’erreur commise en remplaçant dans l’équation (2) vjn par unj = u(j∆x, n∆t). Montrer que
n n
∂2u ∂2u
1 1
εnj = ∆t ( − θ) + c∆x ( − β) + O(∆t2 ) + O(∆x2 ).
∂t2 j 2 ∂x2 j 2
2. En déduire que le schéma (2) est consistant avec l’équation (1).
3. Donner l’ordre en temps et en espace du schéma (2).
4. Etudier la L2 - stabilité du schéma pour les deux cas suivant :
1
(a) β = 0 et θ = .
2
1
(b) β = θ = .
2
Que peut-on dire sur la convergence pour ces deux cas?
Exercice 2
On considère l’équation de la chaleur:
∂u(x, t) ∂ 2 u(x, t)
− = 0, ∀x ∈ R, ∀t > 0. (3)
∂t ∂x2
Partie I
Vérifiez les propriètés suivantes des solutions de l’équation (3).
1
1. Si S1 (x, t) et S2 (x, t) sont deux solutions de (3) alors aS1 (x, t) + bS2 (x, t) est aussi solution de (3), ∀a, b ∈ R.
2. Soit y ∈ R. Si S(x, t) est une solution de (3) alors S(x − y, t) est aussi solution de (3).
3. Soit g une fonction localement intégrable
Z ∞sur R.
si S(x, t) est une solution de (3) alors S(x − y, t)g(y)dy est aussi solution de (3).
−∞
√
4. Si S(x, t) est une solution de (3) alors S( ax, at) est aussi solution de l’équation (3), ∀a ∈ R+ .
Partie II
Résolution d’un problème auxiliaire
∂Q(x, t) ∂ 2 Q(x, t)
− = 0, ∀x ∈ R, ∀t > 0,
∂x2
∂t
( (4)
1, ∀x ≥ 0,
Q(x, 0) =
0, ∀x < 0.
On considère le changement de variables suivant :
x
Q(x, t) = H( √ )
t
1. Montrer que H statisfait l’équation différentielle suivante :
00 p 0 x
H (p) + H (p) = 0, où p = √ . (5)
2 t
0
2. En posant G(p) = H (p), écrire l’équation satisfaite par G et donner son expression.
∂Q 1 x
3. Vérifier que = √ G( √ ), ∀x ∈ R, ∀t > 0.
∂x t t
∂Q
4. Soit S(x, t) = (x, t), vérifier que S est la solution du problème de Cauchy suivant :
∂x
∂S(x, t) ∂ 2 S(x, t)
− = 0 ∀x ∈ R, ∀t > 0,
∂t ∂x2 (6)
S(x, 0) = δ(x), ∀x ∈ R.
(
0, si x 6= 0,
avec δ(x) =
1, si x = 0.
5. On considère maintenant le probème de Cauchy relatif à l’equation de chaleur suivant :
∂u(x, t) ∂ 2 u(x, t)
− = 0 ∀x ∈ R, ∀t > 0,
∂t ∂x2 (7)
u(x, 0) = u (x),
0
avec u0 une fonction localement intégrable sur R.
Z +∞
Montrer que u(x, t) = S(x − y, t) u0 (y)dy est une solution de (7).
−∞