File: HEATEX.FC of Tape: Various/ETH/fc1
(.FC,.FO Focal format converted to 'WRITE' listing) 
C-HEATEX.FC SEP-76

01.01 C HEAT EXCHANGER TOY-PROGRAM
01.02 S LQ=0;T!,%6.03,"UNITS:  CGS, CAL, ATM !";A "LIQUID ?",DU,!;I (DU,=,0YES)1.03;S LQ=1
01.03 T !,"DOWNFLOWING WARM GAS IS CALLED '1', UPFLOWING COLD GAS '2'",!
01.04 T !,:30,"GAS,TUBE #1",:50,"GAS,TUBE #2"
01.05 A !,"VISCOSITY AT 0C (*1E6)",:30,E0(1),:50,E0(2)
01.06 A !,"Cv IN UNITS OF R/2",:30,C0(1),:50,C0(2)
01.07 A !,"MOLAR WEIGHT",:30,MW(1),:50,MW(2)
01.08 S CV(1)=C0(1)/MW(1);S CP(1)=CV(1)+2/MW(1);S CV(2)=C0(2)/MW(2);S CP(2)=CV(2)+2/MW(2)
01.09 A !,"# OF  TUBES",:25,"INNER",N(1),:45,"OUTER",N(2)
01.11 A !,"INNER 0 OF TUBES",:30,D(1),:50,D(2)
01.12 A !,"WALL THICKNESS OF INNER TUBE",:40,TX
01.13 S PE(1)=PI*D(1);S PE(2)=PI*D(2)+N(1)*PI*(D(1)+TX)
01.14 S AR(1)=PI*D(1)^2/4;S AR(2)=PI*D(2)^2/4-N(1)*PI*(D(1)+TX)^2/4
01.15 S PH(1)=PE(1);S PH(2)=N(1)*PI*(D(1)+TX)
01.20 S DH(1)=4*AR(1)/PE(1);S DH(2)=4*AR(2)/PE(2)
01.30 S DE(1)=4*AR(1)/PH(1);S DE(2)=4*AR(2)/PH(2)
01.35 A !,"TEMPERATURE IN K,  SIDE:",:26,"WARM",TH(1),:46,"COLD",TC(2)
01.40 A !,"MASS FLOW RATIO OF GASES 1/2",:40,RM
01.42 I (RM,<=,1) 1.45;S RM=1.001
01.45 A !,"MASS FLOW OF GAS #1",:30,M;S M(2)=M/RM/N(2);S M(1)=M/N(1)
01.55 S PX=PI*(DI+TX/2)
01.60 A !,"CONTACT WALL CONDUCTIVITY",:40,LX
01.65 A !,"AVERAGE PRESSURE",:30,P(1),:50,P(2)
01.90 S TH=TH(1)-TC(2);S T=TC(2)+TH/2;F I=1,2;D 20
01.95 S RA=MC(1)*N(1)/MC(2)/N(2)
01.96 I (LQ,=,1) 1.97;S RA=1.001
01.97 S T0=TH(1)-TH/RA

02.05 D 21
02.10 D 22;F EF=.99,-.01,.79;D 10;D 11;D 23
02.20 G 1.09

10.10 S TC(1)=TH(1)-EF*(TH(1)-T0)
10.20 S AL=1/(HT(1)*PH(1)*N(1))+TX/(PX*LX*N(1))+1/(HT(2)*PH(2)*N(2))
10.25 S LE=AL*MC(1)*RA*FLOG((TC(1)-TC(2))/(TH-RA*(TH(1)-TC(1))))/(RA-1)
10.30 S LE(1)=LE;S LE(2)=LE

11.10 S DP(1)=PS(1)*LE(1)*G(1)^2/(2*981000*DH(1)*RO(1))
11.20 S DP(2)=PS(2)*LE(2)*G(2)^2/(2*981000*DH(2)*RO(2))

20.10 S TR=T/273
20.15 S ET(I)=FEXP(.75*FLOG(TR))*E0(I)*1E-6
20.20 S LA(I)=2*ET(I)*CV(I)
20.25 S G(I)=M(I)/AR(I)
20.30 S PR(I)=ET(I)*CP(I)/LA(I)
20.35 S RE(I)=G(I)*DH(I)/ET(I)
20.40 S NU(I)=.023*FEXP(.8*FLOG(RE(I)))*FEXP(.4*FLOG(PR(I)))
20.45 S HT(I)=NU(I)*LA(I)/DE(I)
20.46 I (LQ,=,1) 20.5;S HT(2)=10000
20.50 S PS(I)=64/RE(I)
20.55 I (RE(I),>,2000)20.6;S PS(I)=.316*FEXP(-.25*FLOG(RE(I)))
20.60 S MC(I)=M(I)*CP(I)
20.65 S RO(I)=.0118*P(I)*MW(I)/T

21.05 T %6.04
21.10 T !,"CV",:30,CV(1),:50,CV(2)
21.15 T !,"CP",:30,CP(1),:50,CP(2)
21.20 T !,"VISCOSITY*1E6",:30,ET(1)*1E6,:50,ET(2)*1E6
21.25 T !,"THERM CONDUCTIVITY*1E6",:30,LA(1)*1E6,:50,LA(2)*1E6
21.30 T !,"PRANDTL",:30,PR(1),:50,PR(2)
21.35 T !,"REYNOLDS",:30,RE(1),:50,RE(2)
21.37 T !,"NUSSELT",:30,NU(1),:50,NU(2)
21.40 T !,"HEAT TRANSF COEFF",:30,HT(1),:50,HT(2)
21.45 T !,"LAM.TURB. PSI FACTOR",:30,PS(1),:50,PS(2)
21.50 T !,"DELTA T MAX",:40,TH
21.53 T !,"LOWEST ATTAINABLE TEMP",:40,T0
21.55 T !,"MASS FLOW OF GAS #2",:40,M(2)
21.60 T !,"AVERAGE DENSITY",:30,RO(1),:50,RO(2)

22.10 T !!,"EFFICIENCY",:20,"END TEMP",:40,"LENGTH",:60,"DELTA P 1",:80,"DELTA P 2"
22.20 T !;F I=1,90;T "-"

23.10 T !,%4.02,EF,:20,%5.01,TC(1),:40,%6.04,LE,:60,%4.02,DP(1),:80,DP(2)