Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Importer
284 vues6 pages

TP mvf2

Transféré par

Andy Ben

Droits d'auteur :

© All Rights Reserved
Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
Télécharger au format pdf ou txt

TP mvf2

Transféré par

Andy Ben
284 vues6 pages

Titre original

tp-mvf2

Copyright

© © All Rights Reserved

Formats disponibles

PDF, TXT ou lisez en ligne sur Scribd

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

© All Rights Reserved
Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
Télécharger au format pdf ou txt
284 vues6 pages

TP mvf2

Transféré par

Andy Ben

Droits d'auteur :

© All Rights Reserved
Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
Télécharger au format pdf ou txt
Vous êtes sur la page 1/ 6

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

Vous aimerez peut-être aussi