3. Analysis of Shell and Tube Heat Exchanger 705
c
data pi/3.1415927/
open(5,file='hx1.in')
open(6,file='hx1.out')
read(5,*) dmi,dmo
read(5,*) Ti_in,To_in
read(5,*) dip,dop,aL,aN
read(5,*) aks
read(5,*) fi,fo
c
di=dip/12.00
do=dop/12.00
fai=pi*di*di/4.
fao=pi*do*do/4.
Ao=pi*do*aL*aN
c Only for the first trial, find properties based on the inlet temps.
iter=0
Ti_out=Ti_in
To_out=To_in
100 continue
iter=iter+1
Ti_avg=0.5*(Ti_in+Ti_out)
To_avg=0.5*(To_in+To_out)
Ts=0.5*(Ti_avg+To_avg)
call intrpl(Ti_avg,cpfi,cpgi,amufi,amug,akfi,akg,Prfi,Prg,
1 sigf,betaf,rofi,rog,anuf,anug,vf,vfg,vg)
call intrpl(To_avg,cpfo,cpgo,amufo,amug,akfo,akg,Prfo,Prg,
1 sigf,betaf,rofo,rog,anuf,anug,vf,vfg,vg)
call htc(Ti_avg,Ts,dmi,aN,di,fai,hi)
call htc(To_avg,Ts,dmo,aN,do,fao,ho)
c
Call Uover(di,do,hi,ho,fi,fo,aks,Uo)
Call NTU(dmi,cpfi,dmo,cpfo,Ao,Uo,Cmin,Cr,aNTU,eff)
Qdot=eff*Cmin*abs(Ti_in-To_in)
If(Ti_in.lt.To_in) go to 1
Ti_out=Ti_in-(Qdot/(dmi*cpfi))
To_out=To_in+(Qdot/(dmo*cpfo))
go to 2
1 continue
Ti_out=Ti_in+(Qdot/(dmi*cpfi))
To_out=To_in-(Qdot/(dmo*cpfo))
2 continue
DTlmtd=Qdot/(Uo*Ao)
if(iter.eq.1) go to 100
Tsn=0.5*(Ti_avg+To_avg)
eps=abs(Tsn-Ts)/Tsn
if(eps.le.1.e-6) go to 101
if(iter.lt.30) go to 100
print *,'Steady-State Iteration Did Not Converge'
stop
101 continue
write(6,3) dip,dop,aL,aN,dmi,dmo,Ti_in,To_in,Ti_out,To_out,Ts,
1 DTlmtd,hi,ho,Qdot,Uo,aNTU,eff
3 format(
1' Tube inside diameter (in):.....',f8.3,5x,
1' Tube outside diameter (in):....',f8.3,5x,/,
1' Tube total length (ft):........',f8.1,5x,
1' Total number of tubes (-):.....',f8.1,5x,/,
1' Tube mass flow rate (lbm/h):...',e8.0,5x,
1' Shell mass flow rate (lbm/h):..',e8.0,5x,/,
1' Tube inlet temperature (F):....',f8.0,5x,
1' Shell inlet temperature (F):...',f8.0,5x,/,
1' Tube outlet temperature (F):...',f8.1,5x,
1' Shell outlet temperature (F):..',f8.1,5x,/,
1' Tube surface temperature (F):..',f8.1,5x,
1' Log mean temp. difference (F):.',f8.1,5x,/,
1' Tube-side HTC (Btu/ft2 h F):...',f8.1,5x,