Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

République Algérienne démocratique et populaire Ministère de l’enseignement sup

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é :  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: h l1 h T ¿ l2  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 uploads/Science et Technologie/ methode-numerique-tp-2.pdf