Université Mohammed V, Rabat ENSAM, Rabat
Département de Mathématiques Appliquées & G. Informatique Année : 2024– 2025
Module : Méthodes Numériques S1 F. Ing. -GM, I2A
TD N°4
Exercice 1 :
Considérons le problème de convection thermique unidimensionnel le long d’une barre cylindrique de longueur L
modélisé par l’équation aux dérivées partielles suivante :
𝑑 𝑑 𝑇(𝑥)
𝑑𝑥
(𝛾𝐴 𝑑𝑥 ) − ℎ𝑃(𝑇(𝑥) − 𝑇∞) = 0 x]0, L[
𝑇(0) = 𝑇𝑏 𝑝𝑜𝑢𝑟 𝑥 = 0, (*)
𝑑𝑇
{ (𝐿)
= 0 𝑝𝑜𝑢𝑟 𝑥 = 𝐿
𝑑𝑥
𝑇∞ est la température externe
ℎ𝑃
1) En Supposant que 𝛾𝐴 et hP sont des constantes et en posant 𝑛2 = 𝛾𝐴 , donner la solution exacte de (*)
2) En notant ∆x le pas de la discrétisation uniforme de [0, L], donner le schéma des volumes finis sur le volume
de contrôle ]𝑥𝑖−−1 , 𝑥𝑖+−1 , [
2 2
3) donner l'approximation des volumes finis
a) au nœud 1, b ) au nœud N
4) Donner le schéma numérique des volumes finis associés au problème (*) à {T1; T2, T3, ...., TN}
5) Montrer que le problème discrétisé se met sous la forme de A.X= F et donner la matrice A et les vecteurs X et F.
Exercice 2
Considérons le problème de conduction thermique unidimensionnel le long d’une barre cylindrique de longueur L=1m
modélisé par l’équation aux dérivées partielles suivante :
𝑑 𝑇(𝑥,𝑡) 𝑑 2 𝑇(𝑥,𝑡)
𝑑𝑡
= 𝛼 𝑑𝑥 2
x]0, 1[
{𝑇(0, 𝑡) = 𝑇𝑔 𝑝𝑜𝑢𝑟 𝑥 = 0, 𝑇(1, 𝑡) = 𝑇𝑑 𝑝𝑜𝑢𝑟 𝑥 = 1 (*)
𝑇(𝑥, 0) = 𝑇0 𝑝𝑜𝑢𝑟 𝑡 = 0
𝛼 est le coefficient de diffusivité thermique que l’on suppose égale à 1, (𝛼 = 1)
1) On note Δx le pas de la discrétisation uniforme de [0, L] et Δt le pas de la discrétisation du temps, t n = n. ∆t ). En
utilisant le schéma d’Euler explicite et le schéma des volumes finis sur le volume de contrôle ]𝑥𝑖−1 , 𝑥𝑖+1 [,
2 2
montrer que
Tn+1
i − Tn
i ∂T n ∂T n
Δx Δt
= [ ( ∂x ) − ( ∂x ) ]
x 1 x 1
i+ i−
2 2
∂T n ∂T n
2) Aux nœuds 1 et N, donner les approximations de ( ) et ( )
∂x x1 ∂x x 1
N+
2 2
∆t
3) On pose β = ∆x2, et le vecteur V = { T1, T2, …., 𝑇𝑁 }𝑡 , V n = V n (tn). tn = n. ∆t
Donner le schéma numérique des volumes finis et temporel V n+1 = A. V n + F et donner la matrice A et le
vecteur F associés.
