Scribd
Importer
111 vues7 pages

Rapport sur l'Équation de Chaleur 2D

Ce document présente un rapport sur Matlab pour la résolution numérique de l'équation de la chaleur en 2D. Il décrit la condition initiale, les conditions aux limites, les schémas d'approximation des dérivées et le schéma explicite découlant de l'équation.

Transféré par

Nassim Dadi
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
111 vues7 pages

Rapport sur l'Équation de Chaleur 2D

Ce document présente un rapport sur Matlab pour la résolution numérique de l'équation de la chaleur en 2D. Il décrit la condition initiale, les conditions aux limites, les schémas d'approximation des dérivées et le schéma explicite découlant de l'équation.

Transféré par

Nassim Dadi
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

Département Modélisation et informatique Scientifique

Rapport sur Matlab


TP2

Réalisé par : Encadré


par :
-ZDAIK kenza Mr. FIKEL
Najib
-LAFRAYAH Abdelmounim
-EL MALHI Abir
-DADI Nassim
Equation de chaleur en 2d :
Exo2
La condition initiale :

𝒖(𝒙, 𝒚, 𝟎) =200 𝒑𝒐𝒖𝒓 (𝒙, 𝒚) ∈]𝟎, 𝑳[ × ]𝟎, 𝒍[


Les conditions aux limites:
𝒖(𝟎, 𝒚, 𝒕) = 100 𝒆𝒕 𝒖(𝑳, 𝒚, 𝒕) = 100 𝒑𝒐𝒖𝒓 𝒕 ≥ 𝟎
𝒖(𝒙, 𝟎, 𝒕) = 100 𝒆𝒕 𝒖(𝒙, 𝒍, 𝒕) = 100 𝒑𝒐𝒖𝒓 𝒕 ≥ 𝟎
K = ∆𝒕 𝒍𝒆 𝒑𝒂𝒔 𝒅𝒆 𝒅𝒊𝒔𝒄𝒓é𝒕𝒊𝒔𝒂𝒕𝒊𝒐𝒏 𝒆𝒏 𝒕𝒆𝒎𝒑𝒔
𝒉 = ∆𝒙=∆y 𝒍𝒆 𝒑𝒂𝒔 𝒅𝒆 𝒅𝒊𝒔𝒄𝒓é𝒕𝒊𝒔𝒂𝒕𝒊𝒐𝒏 𝒆𝒏 𝒆𝒔𝒑𝒂𝒄𝒆

Approximation des dérivées :


∂u u ( xi , yj , tn+1 )−u(xi , yj, tn)
( xi , yj , tn ) =
∂t k

∂2 u u ( xi, yj+1 ,tn )−2 u ( xi, yj , tn ) +u(xi , yj−1 , tn)


( xi , yj ,tn )=
∂ y2 h
2

Et :
σ
u ( xi , yj+ 1 )−2 u ( xi , yj )+ u ( xi , yj−1 ) u ( xi+1 , yj ) −2u ( xi , yj )+u ( xi−1 , yj )
Δu ( xi , yj )=σ ( + )=¿
h2 h2

∂u
( xi , yj , tn ) =¿ σ Δu ( xi , yj )
∂t

En remplaçant par les termes d'approximation dans l'équation on obtient le


schéma explicite :
u ( xi , yj ,tn +1 )−u(xi , yj ,tn) u ( xi , yj+1 ,tn )−2 u ( xi, yj , tn ) +u ( xi , yj−1, tn ) u ( xi+1 , yj ,tn )−2 u ( xi, yj ) +u
=σ ( +
k h
2
h
2
u ( xi , yj ,tn+1 )−u(xi , yj ,tn ) u ( xi , yj+1 , tn ) −2u ( xi , yj, tn )+ u ( xi , yj−1 , tn ) u ( xi+1 , yj, tn )−2 u ( xi , yj ,t
=σ +σ
k h
2
h
2

σk σk
u ( xi, yj , tn+1 ) = 2
(u ( xi+1 , yj, tn )+ u ( xi−1, yj , tn ) +u ( xi , yj+1 , tn ) +u ( xi , yj−1 ,tn ) +(1−4 2 )u ( xi , yj, tn ))
h h
U(k+1)=AU(k)
A=I-kM

M= ( DN N
⋱
⋯ 0 0
N ¿ ¿
⋱ ¿ N ¿ ¿ ⋮¿ ¿ N ¿ ⋱ ¿ N ¿ 0 ¿ ¿ ⋯ ¿ N ¿ D ¿ )

( )
2 −1
⋯ 0 0
h2 h2 −1 −1 −1 −1 2
D= ⋱¿ 2 ¿
¿⋮ ¿ ¿ 2 ¿ ⋱ ¿ 2 ¿ 0 ¿ ¿ ⋯ ¿ 2 ¿ 2 ¿
−1 −1 h h h h h
2
⋱ 2
¿ ¿
h h

−1
N=- h 2 Id
Le programe sous matlab
Exo 1
Par analogie avec l’exercice précèdent, on a annulé
le second terme

Vous aimerez peut-être aussi