1/7

I. Solution Analytique
𝜕𝑇 𝜕𝑇
= 𝛼 (1)
𝜕𝑡 𝜕𝑥 2

𝑇(0, 𝑡) = 100𝑜 𝑇(1, 𝑡) = 25𝑜

𝜋 𝜋
𝑇(𝑥, 0) = 25 sin ( 2 𝑥) + 100 𝑐𝑜𝑠 (2 𝑥)

Utilisons la méthode de changement des variables en posant

𝑇 (𝑥, 𝑡): 𝑓(𝑡) × 𝑔(𝑥) (2)
𝜕𝑓(𝑡) 𝑔(𝑥) 𝜕𝑓(𝑡)𝑔(𝑥)
En remplaçant 𝑇 (𝑥, 𝑡): par 𝑓(𝑡) × 𝑔(𝑥) , 𝑜𝑛 𝑜𝑏𝑡𝑖𝑒𝑛𝑡 = 𝛼
𝜕𝑡 𝜕𝑥 2

Si en plus on pose 𝛼 = 1, 𝑜𝑛 𝑜𝑏𝑡𝑖𝑒𝑛𝑡 𝑓 ′ (𝑡)𝑔 (𝑥) = 𝑓(𝑡)𝑔′′ (𝑥)

En regroupant les termes en fonction des variables, on obtient
𝑓 ′ (𝑡) 𝑔′′ (𝑥)
=
𝑓(𝑡) 𝑔 (𝑥)
Les variables étant séparées, l’égalité entre les deux expressions des variables
différentes est une constante donc :
1 1 1 𝑑𝑓 1 𝑑2 𝑔
× 𝑓 ′ (𝑥) = 𝑔 (𝑥) × 𝑔′′ (𝑥) = 𝑘 → × 𝑑𝑡 = 𝑔 =𝑘
𝑓(𝑡) 𝑓 𝑥2𝑑
𝑑𝑓 𝑓′
→ = 𝑘𝑑𝑡 → − 𝑘𝑓 = 0 (3)
𝑓 𝑔

𝑓 ′ − 𝑘𝑓 = 0
{
𝑔′′ − 𝑘𝑔 = 0
𝑓 ′ + 𝑤𝑓 = 0 (1)
Posons 𝑘 = − 𝑤 2 → {
𝑔′′ + 𝑤𝑔 = 0 (2)
D’après les formules de résolution des équations différentielles, 𝑓(𝑡) a pour solution
2𝑡
𝑓(𝑡) = 𝐴𝑒 −𝑤
Par contre, pour on a 𝑟 2 + 𝑤 2 𝑟 = 0 comme équation caractéristique.
𝑟1 = 0
Les valeurs probables de 𝑟 sont (𝑟 + 𝑤 2 ) = 0 → { →
𝑟2 + 𝑤 2 = 0
𝑟1 = 0
{
𝑟2 = −𝑤 2
𝑔(𝑥) = 𝐵1 𝑐𝑜𝑠 𝑤𝑥 + 𝐵2 sin 𝑤𝑥 donc 𝑇 (𝑥, 𝑡) = 𝑓(𝑡)𝑔(𝑥) devient :
2𝑡
𝑇 (𝑥, 𝑡) = 𝐴𝑒 −𝑤 (𝐵1 𝑐𝑜𝑠 𝑤𝑥 + 𝐵2 sin 𝑤𝑥)
Exploitons les paramètres initiaux pour déterminer 𝐴, 𝐵1 , 𝐵2 et 𝑤.
𝜋 𝜋
D’après l’énoncé 𝑇(0, 𝑡) = 100; 𝑇(1, 𝑡) = 25 et 𝑇(𝑥, 0) = 25 sin ( 2 𝑥) + 100 𝑐𝑜𝑠 ( 2 𝑥)
2𝑡
𝑇(0, 𝑡) = 100 → 𝐴𝑒 −𝑤 (𝐵1 𝑐𝑜𝑠 𝑤(0) + 𝐵2 sin 𝑤(0)) = 100𝑜
2𝑡
𝐴𝑒 −𝑤 (𝐵1 ) = 100𝑜 (1)
 2/7

2𝑡
𝑇(1, 𝑡) = 25 → 𝐴𝑒 −𝑤 (𝐵1 𝑐𝑜𝑠 𝑤(1) + 𝐵2 sin 𝑤(1)) = 100𝑜
2𝑡
→ 𝐴𝑒 −𝑤 (𝐵1 𝑐𝑜𝑠 𝑤 + 𝐵2 sin 𝑤) = 25𝑜
𝜋 𝜋 2(0)
𝑇(𝑥, 0) 25 sin ( 𝑥) + 100 𝑐𝑜𝑠 ( 𝑥) → 𝐴𝑒 −𝑤 (𝐵1 𝑐𝑜𝑠 𝑤𝑥 + 𝐵2 sin 𝑤𝑥)
2 2
𝜋 𝜋
= 25 sin ( 2 𝑥) + 100 𝑐𝑜𝑠 ( 2 𝑥)
𝜋 𝜋 𝜋
𝑇(𝑥, 0) 25 sin ( 2 𝑥) + 100 𝑐𝑜𝑠 ( 2 𝑥) → 𝐴(𝐵1 𝑐𝑜𝑠 𝑤𝑥 + 𝐵2 sin 𝑤𝑥) = 25 sin ( 2 𝑥) +
𝜋
100 𝑐𝑜𝑠 (2 𝑥)

𝜋 𝜋
→ 𝐴(𝐵1 𝑐𝑜𝑠 𝑤𝑥 + 𝐵2 sin 𝑤𝑥) = 25 (sin ( 𝑥) + 4 𝑐𝑜𝑠 ( 𝑥))
2 2
𝜋
𝐴 = 25 𝑤=2
Par identification, { {
𝐵1 = 4 𝐵2 = 1
𝜋2𝑡
𝜋 𝜋
𝑇 (𝑥, 𝑡) = 25𝑒 − 2 (4𝑐𝑜𝑠 𝑥 + sin 2 𝑥)
2
 3/7

II. Solution Numérique

I.1 Méthode Explicite
I.1.1 Discrétisation des membres de l’équation de la chaleur.
Rappel de la formule de Taylor pour une fonction 𝑓
𝒏
𝒇(𝒌) (𝒂) (1)
𝒇(𝒂 + 𝚫𝒙) = ∑ (𝚫𝒙)𝒌 + 𝑹𝒏(𝚫𝒙)
𝒌!
𝒌=𝟎

Avec a є D𝑓, et Δ𝑥 є R ∗
𝜕𝑇
I.1.1.1 Discrétisation du membre 𝜕𝑡

Par application de l’équation (1) au temps 𝑛, et à l’ordre 1, on a:

1
𝑛+1 (Δ𝑡)𝑛 𝜕 𝑘 𝑇
𝑇 𝑖
= 𝑇(𝑥𝑖 , 𝑡𝑛 + Δ𝑡) = ∑ + 𝑅𝑛(Δ𝑡) (2)
𝑘! 𝜕𝑡𝑘
𝑘=0

(Δ𝑡)0 𝜕0 𝑇 𝑛 (Δ𝑡)1 𝜕1 𝑇 𝑛 (3)

= ( 𝜕𝑡 0 ) 𝑖 + ( 𝜕𝑡 1 ) 𝑖 + 𝑅1(Δ𝑡)
0! 1!

𝜕𝑇 𝑛 (4)
≈ 𝑇 𝑛𝑖 + Δ𝑡 ( )
𝜕𝑡 𝑖

D’où:
𝒏+𝟏
𝝏𝑻 𝒏 𝑻 𝒊
− 𝑻 𝒏𝒊 (5)
( ) =
𝝏𝒕 𝒊 𝚫𝒕
𝛛𝟐 𝐓
I.1.1.1 Discrétisation du membre 𝛛𝐭 𝟐

De même, par application de la formule de Taylor au point 𝑥𝑖 , à l’ordre 2, on a:

2
𝑛 (Δ𝑥)𝑛 𝜕 𝑘 𝑇
𝑇 𝑖+1
= 𝑇(𝑥𝑖 + Δ𝑥, 𝑡𝑛 ) = ∑ + 𝑅𝑛(Δ𝑥) (6)
𝑘! 𝜕𝑥 𝑘
𝑘=0

(Δ𝑥)0 𝜕0 𝑇 𝑛 (Δ𝑥)1 𝜕1 𝑇 𝑛 (Δ𝑥)2 𝜕2 𝑇 𝑛 (7)

= (𝜕𝑥 0) 𝑖 + (𝜕𝑥 1 ) 𝑖 + (𝜕𝑥 2 ) 𝑖 + 𝑅2(Δ𝑡)
0! 1! 2!

𝑛 𝜕𝑇 (Δ𝑥)2 𝜕2 𝑇 𝑛
𝑇 𝑖+1
≈ 𝑇 𝑛𝑖 + Δ𝑡 (𝜕𝑥 ) 𝑛𝑖 + (𝜕𝑥 2 ) 𝑖
2 (8)
De la même façon, nous calculons:

2
𝑛 (−Δ𝑥)𝑛 𝜕 𝑘 𝑇
𝑇 𝑖−1
= 𝑇(𝑥𝑖 − Δ𝑥, 𝑡𝑛 ) = ∑ + 𝑅𝑛(Δ𝑥)
𝑘! 𝜕𝑥 𝑘
𝑘=0

(−Δ𝑥)0 𝜕0 𝑇 𝑛 (−Δ𝑥)1 𝜕1 𝑇 𝑛 (−Δ𝑥)2 𝜕2 𝑇 𝑛

= (𝜕𝑥 0 ) 𝑖 + (𝜕𝑥 1 ) 𝑖 + (𝜕𝑥 2) 𝑖 + 𝑅2(Δ𝑡)(9)
0! 1! 2!
 4/7

𝑛 𝜕𝑇 𝑛 (Δ𝑥)2 𝜕2 𝑇 𝑛
𝑇 𝑖−1 ≈ 𝑇 𝑛𝑖 − Δ𝑡 ( ) + (𝜕𝑥 2 ) 𝑖 (10)
𝜕𝑥 𝑖 2

En additionnant les équations (8), et (10), on obtient:

𝜕2 𝑇 𝑛 (11)
𝑛
𝑇 𝑖+1 𝑛
+ 𝑇 𝑖−1 ≈ 2𝑇 𝑛𝑖 − (Δ𝑥)2 ( )
𝜕𝑥 2 𝑖
D’où:
𝒏
𝝏𝟐 𝑻 𝒏 𝑻 𝒊−𝟏
− 𝟐𝑻 𝒏𝒊 + 𝑻 𝒏
𝒊+𝟏
(12)
𝜶 ( 𝟐) = 𝜶
𝝏𝒙 𝒊 (𝚫𝒙)𝟐
D’après les équations (5) et (12), on en déduit que:

𝜕𝑇 𝑛 𝑇
𝑛+1
𝑖
− 𝑇 𝑛𝑖 𝜕2 𝑇 𝑛 𝑇
𝑛
𝑖−1
𝑛
− 2𝑇 𝑖 + 𝑇
𝑛
𝑖+1
( ) = = 𝛼( ) =𝛼 (13)
𝜕𝑡 𝑖 Δ𝑡 𝜕𝑥2 𝑖 (Δ𝑥) 2

D’où il vient que:

𝚫𝒕
+ (𝟏 − 𝟐𝜶) (𝚫𝒙)𝟐 𝑻 𝒏𝒊 + 𝜶 (𝚫𝒙)𝟐 𝑻 𝒊+𝟏
𝚫𝒕𝒏 𝒏 𝚫𝒕
𝑻 𝒏+𝟏
𝒊
= 𝜶 (𝚫𝒙)𝟐 𝑻 𝒊−𝟏 (14)

Posons :
Δ𝑡
β = 𝛼 (Δ𝑥)2 (15)

Δ𝑡
λ = 𝛼(1 − 2𝛼) (Δ𝑥)2 (16)

D’où :

𝑻 𝒏+𝟏 𝒏 + 𝛌𝑻 𝒏 + 𝛃𝑻 𝒏
= 𝛃𝑻 𝒊−𝟏 (17)
𝒊 𝒊 𝒊+𝟏

Le problème qu’on cherche alors à résoudre se ramené a la recherche de la valeur

de la température au point 𝑥𝑖 de la tige, a l’instant 𝑛 +1, sachant que la température
au temps 𝑛, est connue, au point 𝑥𝑖 , 𝑥𝑖+1, et 𝑥𝑖−1..
Pour :
𝑖 = 1: 𝑇 𝑛+1
1
= β𝑇 𝑔𝑛 + λ𝑇 𝑛1 + β𝑇 𝑛2 (18)

𝑖 = 2: 𝑇 𝑛+1
2
= β𝑇 𝑛1 + λ𝑇 𝑛2 + β𝑇 𝑛3 (19)
𝑖 = 3: 𝑇 𝑛+1
3
= β𝑇 𝑛
2+ λ𝑇 𝑛3 + β𝑇 𝑛
4 (20)
 5/7

𝑖 = N − 1: 𝑛+1
𝑇 𝑁−1 𝑛 + λ𝑇 𝑛 + β𝑇
= β𝑇 𝑁−2 (21)
𝑁 𝑑

Définissons les matrices et vecteurs suivants :

La solution peut se mettre sous la forme matricielle ci-dessous:

A 𝑻 𝒏+𝟏
𝒊
= B 𝑻 𝒏𝒊 + 𝑻𝒄 (22)

Nous utilisons le logiciel Matlab pour faire les simulations (voir résultats de la
simulation dans le CD)

I.1 Schéma Implicite paramètre par Ө

La méthode de Euler, base sur l’application du théorème de Picard définie, pour
deux fonctions f et y continues, avec 𝑦 ′ = 𝑓(𝑥, 𝑦) 𝑒𝑡 𝑦(𝑥0) = 𝑦0 , l’égalité ci-
dessous :

𝑦 𝑛+1 = 𝑦 𝑛 + Δ𝑥[(1 − Ө) 𝑓(𝑥 𝑛 , 𝑦 𝑛 ) + Ө𝑓(𝑥 𝑛+1 , 𝑦 𝑛+1), 𝑛 = 0, … . , 𝑁 − (23)

1, 𝑒𝑡 Ө є [0, 1]
Posons :
(24)
𝑦 = 𝑇(𝑥, 𝑡)
𝜕 2 𝑇(𝑥, 𝑡)
𝑓(𝑥, 𝑡) = 𝛼 ( ) (25)
𝜕𝑥 2
D’après l’équation de la chaleur, et la relation (12), on a :

(26)
 6/7

𝑛
′
𝜕𝑇(𝑥, 𝑡) 𝜕 2 𝑇(𝑥, 𝑡) 𝑇 𝑖−1
− 2𝑇 𝑛𝑖 + 𝑇 𝑛
𝑖+1
𝑦 = = 𝛼( ) = 𝑓(𝑥, 𝑡) = 𝛼
𝜕𝑡 𝜕𝑥 2 (Δ𝑥)2
Par application de la méthode d’Euler, on peut écrire :

𝑛
𝑛+1 𝑇 𝑖−1 −2𝑇𝑛𝑖+𝑇 𝑖+1
𝑛
𝑇 𝑛+1 𝑛+1 𝑛+1
𝑖−1 −2𝑇 𝑖 +𝑇 𝑖+1
𝑇 𝑖
= 𝑇 𝑛𝑖 + Δ𝑡[(1-Ө)( 𝛼 (Δ𝑥)2
) + Ө(𝛼 (Δ𝑥)2
)] (27)

𝑇 𝑛+1 𝑛+1 𝑛+1

𝑖−1 −2𝑇 𝑖 +𝑇 𝑖+1
𝑛
𝑇 𝑖−1−2𝑇 𝑛𝑖+𝑇 𝑖+1
𝑛
(28)
𝑇 𝑛+1
𝑖
− ӨΔ𝑡(𝛼 (Δ𝑥)2
)= 𝑇 𝑛𝑖 + Δ𝑡[(1-Ө)( 𝛼 (Δ𝑥)2
)]

𝛼ӨΔ𝑡 2𝛼ӨΔ𝑡 𝛼ӨΔ𝑡 𝛼(1−Ө)Δ𝑡 2𝛼(1−Ө)Δ𝑡 𝛼(1−Ө)Δ𝑡

−( 𝑇 𝑛+1
𝑖−1
+ (1 + (Δ𝑥)2
) 𝑇 𝑛+1
𝑖
− (Δ𝑥)2 𝑇 𝑛+1
𝑖+1
= 𝑛
𝑇 𝑖−1+(1 − (Δ𝑥)2
) 𝑇 𝑛𝑖+ 𝑛
𝑇 𝑖+1 (29)
Δ𝑥)2 (Δ𝑥)2 (Δ𝑥)2

Δ𝑡
Posons : λ = 𝛼 (Δ𝑥)2 (30)

Alors:

−Өλ 𝑇 𝑛+1
𝑖−1
+ (1 + 2Өλ) 𝑇 𝑛+1
𝑖
− Өλ 𝑇 𝑛+1
𝑖+1
𝑛
= (1 − Ө)λ 𝑇 𝑖−1+(1 − 2(1 − Ө)λ ) 𝑇 𝑛𝑖+(1 − Ө)λ 𝑇 𝑖+1
𝑛 (31)

De même, cette dernière relation peut se mettre sous la forme matricielle ci-
dessous :

A 𝑻 𝒏+𝟏
𝒊
= B 𝑻 𝒏𝒊 + 𝑻 𝑪 (32)
Avec les matrices ci-dessous :
 7/7

Nous utilisons le logiciel Matlab pour faire les simulations avec Ө = 1 (on retrouve les
matrices vues en cours : voir résultats de la simulation dans le CD).