Design og modellering af metanolanlæg til VEnzin-visionen Bilag

methanol.for
d:/DTU/Eksamensprojekt/bilag/
elseif(abs(domega).gt.abs(omega_old*.0025)) then
if(abs(dtau).gt.abs(tau_old*.0025)) tau_new=tau_old
$ +dtau/abs(dtau)*tau_old*.0025*abs(1.d0
$ −tau_new/tau_old)/abs(1.d0−tau_new/tau_old)
$ /factor
omega_new=omega_old+domega/abs(domega)*omega_old
$ *.0025/factor
endif
else
if ((abs(dtau)).gt.(abs(tau_old*.05))) tau_new
$ =tau_old+dtau/abs(dtau)*tau_old*.05
if ((abs(domega)).gt.(abs(omega_old*.05)))
$ omega_new=omega_old+domega/abs(domega)*omega_old*.05
endif
else
if ((abs(dtau)).gt.(abs(tau_old*.2/factor))) tau_new
$ =tau_old+dtau/abs(dtau)*tau_old*.2/factor
if (domega.lt.(−(omega_old*.9/factor))) then
omega_new=omega_old−omega_old*.9/factor
elseif (domega.gt.(omega_old*9/factor)) then
omega_new=omega_old+omega_old*9/factor
else
omega_new=omega_old+domega
endif
endif
call calculate(f1_new_omega,tau_old,omega_new,in1)
call calculate(f1_new_tau,tau_new,omega_old,in1)
call calculate(f2_new_omega,tau_old,omega_new,in2)
call calculate(f2_new_tau,tau_new,omega_old,in2)
f1_old=f1_new
f2_old=f2_new
i=i+1
goto 5
else
goto 10
endif
10 if (.not.found) then
error=5
c if (fiters) print*
c $ ,’Warning: No convergence in iterer_omega_tau’,in1,f1
c $ ,f1_new,abs(df1/f1),in2,f2,f2_new,abs(df2/f2),1/(omega_old
c $ *rho_ref*M)*1.d−6,1/(omega_new*rho_ref*M)*1.d−6,T_ref
c $ /tau_old−T_0C,T_ref/tau_new−T_0C,factor,i
if ((abs(df1/f1).gt.0.01).or.(abs(df2/f2).gt.0.01)) print*
$ ,’Warning: No convergence in iterer_omega_tau’,abs(df1/f1)
$ ,abs(df2/f2)
endif
tau=tau_new
omega=omega_new
c if (fiters) print*,’iterer_omega_tau. Antal iterationer=’,i
return
end
subroutine iterer_omega(f,tau,omega,in,error)
integer in,i,error,IMAX
double precision f,tau,omega,omega_old,omega_new,domega,f_new
$ ,f_old,df,dfdomega,EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M
$ ,T_0C
logical found,saturation,fiters
common/constants/EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C,IMAX
$ ,fiters
common/switches/iterate_rho_l
c omega_old=omega
omega_old=0.032148/rho_ref
call calculate(f_old,tau,omega_old,in)
omega_new=omega_old*(1d0+1d−3)
if (in.eq.1) then
saturation=.true.
else
saturation=.false.
endif
i=0
found=.false.
5 if ((i .lt. IMAX) .and.(.not.found)) then
if (saturation) call test_saturation(tau,omega_new,error,f)
call calculate(f_new,tau,omega_new,in)
if (error.eq.2) then
8/19
19−03−2007