1 subroutine transit(compo,dt0)
4 dimension tn(7),t(7),rf(7)
5 dimension tparoi(7),rparoi(7)
6 dimension maille(2),tflu(2,7)
29 c initialisation de temperature fluide a t=0
31 CALL cpeRE(compo,CP_TEMPS, ti, npas+1, 'tfluide', 6,
33 IF( info.NE. CPOK )GO TO 9000
35 c initialisation de la resistance thermique fluide a t=0
37 CALL cpeRE(compo,CP_TEMPS, ti, npas+1, 'rfluide', 6,
39 IF( info.NE. CPOK )GO TO 9000
41 c boucle temporelle jusqu'a 100.
43 do while( ti.lT.100. )
46 c lecture de la temperature de paroi entre ti et tf
48 CALL cplRE(compo,CP_TEMPS,ti, tf, npas+1,'tparoi', 6,
49 & nval, tparoi(2), info)
50 IF( info.NE. CPOK )GO TO 9000
52 c lecture de la resistance solide de bord entre ti et tf
54 CALL cplRE(compo,CP_TEMPS,ti, tf, npas+1, 'rparoi', 6,
55 & nval, rparoi(2), info)
56 IF( info.NE. CPOK )GO TO 9000
58 c calcul de la temperature a tf
61 smb=ro/dt*tn(i)+deb*t(i-1)+tparoi(i)/(rparoi(i)+rf(i))
62 t(i)=smb/(ro/dt+deb+1./(rparoi(i)+rf(i)))
65 write(6,*)'FLUID:temps=',tf,' temperature de sortie canal=',t(7)
68 c ecriture de la temperature fluide a tf
70 CALL cpeRE(compo,CP_TEMPS, tf, npas+1, 'tfluide', 6,
72 IF( info.NE.CPOK )GO TO 9000
74 c ecriture de la resistance thermique fluide a tf
76 CALL cpeRE(compo,CP_TEMPS, tf, npas+1, 'rfluide', 6,
78 IF( info.NE.CPOK )GO TO 9000
93 CALL cpfin(compo,CP_ARRET,info)
96 subroutine perma(compo)
99 dimension tparoi(7),rparoi(7)
110 c initialisation de temperature fluide a i=0
114 CALL cpeRE(compo,CP_ITERATION, ti, iter , 'tfi', 6,
116 IF( info.NE. CPOK )GO TO 9000
118 c boucle temporelle jusqu'a 100.
120 do while( iconv .EQ. 0)
122 c lecture de la temperature de paroi iteration iter
124 CALL cplRE(compo,CP_ITERATION,ti, tf, iter , 'tpi', 6,
125 & nval, tparoi(2), info)
126 IF( info.NE. CPOK )GO TO 9000
128 c calcul de la temperature
131 smb=deb*t(i-1)+tparoi(i)/(rparoi(i)+rf(i))
132 t(i)=smb/(deb+1./(rparoi(i)+rf(i)))
135 c ecriture de la temperature fluide a iter+1
137 CALL cpeRE(compo,CP_ITERATION,ti,iter+1, 'tfi', 6,
139 IF( info.NE. CPOK )GO TO 9000
142 write(6,*)'iter = ',iter,' temperature de sortie canal = ',t(7)
144 c lecture du flag de convergence iconv
146 CALL cplEN(compo,CP_ITERATION,ti, tf, iter , 'iconv', 1,
147 & nval, iconv , info)
148 write(6,*)"info:",info
149 write(6,*)"FLUIDE:",iter,iconv
152 IF( info.NE. CPOK )GO TO 9000
154 if(iconv.eq.1)go to 9000
159 CALL cpfin(compo,CP_ARRET,info)