1-Introduction:
Ce travail appliqué concerne le transfert de chaleur par conduction thermique avec terme Source en
régime permanent et unidimensionnelle dans une plaque.
• L’épaisseur plaque L = 2 cm
• Conductivité thermique constante λ=0.5W/MK
• Ta=100°C et Tb=200°C.
• Génération de chaleur q = 1000 kW/m3.
• dx=L/n=0.004.
2-rappel théorique :
Équation de discrétisation pour le nœud frontière 1 :
TP:MVF 1
Équation de discrétisation pour les points nodaux 2, 3,4,
Équation de discrétisation pour le nœud frontière 5
Regroupement des résultats :
nœud Aw Ae Sp Ap=aw+ae-sp Su
1 0 125 -250 357 29000
2 125 125 0 250 4000
3 125 125 0 250 4000
4 125 125 0 250 4000
5 125 0 -250 375 54000
TP:MVF 2
-introduire les conditions aux limites dans le programme
!****************** Program avec_terme_source_positif*************************
dimension t(500),x(500),aw(500),ae(500),ap(500),Sp(500),Su(500),su0(500),ap0(500)
real L,k
open(0,file="input.dat")
!The entries are in the following order
!L ta tb n k q s
open(1,file="output.dat")
read(0,*)L,ta,tb,n,k,q,s
!boundary node 1
aw(1)=0
ae(1)=s*(k/(L/n))
Sp(1)=-2*s*(k/(L/n))
ap(1)=aw(1)+ae(1)-Sp(1)
Su(1)=q*s*(L/n)+2*ta*s*(k/(L/n))
!boundary node n
aw(n)=s*(k/(L/n))
ae(n)=0
Sp(n)=-2*s*(k/(L/n))
ap(n)=aw(n)+ae(n)-Sp(n)
Su(n)=q*s*(L/n)+2*tb*s*(k/(L/n))
x(1)=(L/n)/2
!boundary node 2:n-1
do i=2,n-1
x(i)=x(i-1)+(L/n)
aw(i)=s*(k/(L/n))
ae(i)=s*(k/(L/n))
TP:MVF 3
Sp(i)=0
ap(i)=aw(i)+ae(i)-Sp(i)
Su(i)=q*s*(L/n)
enddo
x(n)=x(n-1)+(L/n)
write(1,*)" i aw ae Su Sp ap"
do i=1,n
write(1,*)i,aw(i),ae(i),Su(i),Sp(i),ap(i)
enddo
do i=1,n
ap0(i)=ap(i)
Su0(i)=Su(i)
enddo
!TDMA
do i=2,n
ap0(i)=ap0(i)+ae(i-1)*(-aw(i)/ap0(i-1))
Su0(i)=Su0(i)-Su0(i-1)*(-aw(i)/ap0(i-1))
enddo
t(n)=Su0(n)/ap0(n)
do i=n-1,1,-1
t(i)=(Su0(i)+ae(i)*t(i+1))/ap0(i)
enddo
write(1,*)" i x t"
do i=1,n
write(1,*)i,x(i),t(i)
enddo
end
TP:MVF 4
-Comparer le résultat avec la solution analytique :
TP:MVF 5
Conclusion :
Le résultat correspond à la solution analytique, La solution exacte est un polynôme du
deuxième degré T(X)=-10^6X^2+25000X+104.
TP:MVF 6
