Scribd
Importer
193 vues5 pages

Étude de la convergence en énergétique

Ce document présente une étude numérique de la résolution de l'équation de la chaleur en deux dimensions par la méthode de Gauss-Seidel. Différents paramètres sont étudiés comme la valeur du critère de convergence, la taille du maillage et la température initiale. Les résultats sont présentés sous forme de tableaux.

Transféré par

Youcef Bony
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
193 vues5 pages

Étude de la convergence en énergétique

Ce document présente une étude numérique de la résolution de l'équation de la chaleur en deux dimensions par la méthode de Gauss-Seidel. Différents paramètres sont étudiés comme la valeur du critère de convergence, la taille du maillage et la température initiale. Les résultats sont présentés sous forme de tableaux.

Transféré par

Youcef Bony
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

République Algérienne démocratique et

populaire
Ministère de l’enseignement supérieur et de la
recherche scientifique 
Université kasdi merbah Ouargla
Faculté des sciences technologies
Département de génie mécanique
Spécialité : énergétique

Presented by: Directed by:


 Ben mahcen Youcef. PR: Boubekri
 Ben mahcen Mouhamed Islam.

Group:
M1En1

Promotion:
2020/2021
A. Présentation du problème :
 le schéma de millage utilisé :

T¿

l2

l1


B. Traitements:
1. Paramètre d'influence sur la convergence:
 Etude l'effet de la valeur du critère de convergence « Epsilon » :

On a : T ¿= 100, et h = 1

Epsilon K

0,00001 142

0,0001 105

0,001 68

0,01 31

Le commentaire sue les résultats:


 Etude l’effet de la taille du millage « h » :

On a : T ¿= 100, et EPS = 0.0001

h K

0.5 465

1 142

1.5 74

2 42

Le commentaire sue les résultats:

 Etude l’effet du champ initial de la température :

On a : EPS = 0.0001, et h = 1

T¿ K

50 152

100 142

150 156

200 134

Le commentaire sue les résultats:


2.
C. Résultats en sortie formatées :

dimension t(100,100),

told(100,100)

real L1,L2

open (9,file='sortie.dat',status='unknown')

!lecture des données du probleme print*,'entrez les valeurs de L1 L2 et h et epsilon'

read(,)L1,L2,h,epsi print*,'entrez les valeurs de Tin T3 et Tdroite'

read(*,*)Tin,T3,Tdr

!Definition du maillage m=anint(L1/h)+1

n=anint(L2/h)+1

! introduction de champs initiaux

k=0

do 10

i=2,n-1

do 13

j=2,m-1

t(i,j)=Tin told(i,j)=Tin 13 continue 10 continue

! introduction de conditions aux limites

do 20 i=1,n t(i,1)=t3 20 t(i,m)=tdr do 21 j=1,m t(1,j)=t3 21 t(n,j)=t3

!clacul iteeratif des temperatures par la method de gauss seidel

88 k=k+1 do 100 i=2,n-1

do 100 j=2,m-1

t(i,j)=0.25*(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1)) 100 continue

!test de convergence

do i=2,n-1

do j=2,m-1 crc=abs((t(i,j)-told(i,j))/told(i,j))

if(crc.gt.epsi)

goto 101

end do

end do
goto 108 

do 102 i=2,n-1

do 102 j=2,m-1 told(i,j)=t(i,j) 102 continue

goto 88 108

write(,)k

write (9,*)'resaultat de'L1'L2'h',.epsi',Tin',T3',Tdr,k,'iteration'

stop

end

Vous aimerez peut-être aussi