10 REM ******** Program P4485 Rev. 6 Date: 30 Januarg 1987 ******** 12 M$-' EVALUATION OF THERMAL PERFORMANCE: WATER COOLING TOWER ' 14 L$=' *************************************a***************************' 16 N$=' TEST METHOD TO NS 4485: PART 2: 1987 ' 18 Q$=' INPUT FROM VDU. OUTPUT TO VDU ' 20 PRINT L$ 22 PRINT M$ 24 PRINT N$ 26 PRINT Q$ 28 PRINT L$ 30 PRINT 32 PRINT ' ********************** INPUT DESIGN DATA *********************' 34 PRINT 36 PRINT ' IS IT A MECHANICAL OR A NATURAL DRAUGHT TOWER ? (Input M or N)' 38 INPUT LINE I$ 40 IF SUB(I$,1)=’M’ THEN 46 42 IF SUB(I$,1)=’N’ THEN 46 44 GOTO 36 46 D$=SUS(I$,1) 48 PRINT DESIGN INLET WATER FLOW (m3/s) ?‘ 50 INPUT Q1 52 PRINT ‘ DESIGN HOT WATER TEMPERATURE (C) ?‘ 54 INPUT T1 56 PRINT ‘ DESIGN RECOOLED WATER TEMPERATURE (C) ?‘ 58 INPUT T2 60 IF D$=’M’ THEN 66 62 PRINT ‘ DESIGN DRY BULB TEMPERATURE OF INLET AIR (C) ?‘ 64 INPUT Dl 66 PRINT ‘ DESIGN WET BULB TEMPERATURE OF INLET AIR (C) ?‘ 68 INPUT W1 70 IF D$=’N’ THEN 76 72 PRINT ‘ DESIGN FAN POWER (kW) ?‘ 74 INPUT P3 76 PRINT ‘ IS THE DESIGN BASED ON ATMOS. PRESSURE. OR ALTITUDE? (Input P or A)’ 78 INPUT LINE P$ 80 IF P$=’P’ THEN 86 82 IF P$=’A’ THEN 92 84 GOTO 76 86 PRINT ‘ DESIGN ATMOSPHERIC PRESSURE (kPa) ?‘ 88 INPUT Pl 90 GOTO 96 92 PRINT ‘ DESIGN ALTITUDE (m above sea ievel) ?‘ 94 INPUT Al 96 PRINT ‘ DESIGN AIR/WATER RATIO L/G ?‘ 98 INPUT JI 100 PRINT ‘ "n" VALUE FOR TOWER CHARACTERISTIC CURVE KaV/L=C*(L/G)’n ?‘ 102 INPUT N1 104 N1=ABS(Ni) 106 PRINT 108 PRINT ' ********************** INPUT TEST DATA ***********************' 110 PRINT 112 PRINT ' TEST INLET WATER FLOW (m3/s) ?' 114 INPUT 02 116 PRINT ' TEST HOT WATER TEMPERATURE (C) ?' 118 INPUT T5 120 PRINT ' TEST RECOOLED WATER TEMPERATURE (C) ?' 122 INPUT T6 124 IF D$=’M’ THEN 130 126 PRINT ' TEST DRY BULB TEMPERATURE OF INLET AIR (C) ?' 128 INPUT D2 130 PRINT ' TEST WET BULB TEMPERATURE OF INLET AIR (C) ?' 132 INPUT W2 134 IF D$=’N’ THEN 140 136 PRINT ' TEST FAN POWER (kW) ?‘ 138 INPUT P4 140 PRINT / Calculating ... 142 REM SET LOOP TOLERANCES FOR KaV/L/ LIG (TEST)/ ENTHALPY BALANCE 144 REM AND EXPECTED WATER TEMPERATURE ( Ref. Line 7016 — 7018 ) 146 Y1=1.E-05 148 Y2=2.E-05 150 Y3=.0015 152 Y4=.0001 154 REM **************** READ DATA/ DEFINE FUNCTIONS ***************** 156 READ A, B, C, D, E, F, TO, T9, M1, M2, C1, C2, C3, LO, RO, P9, B1,B2 158 DIM N(4), W(4), A(4), D(4), T(6), E(2) 160 READ N(1), N(2), N(3), N(4), C0 162 GOSUB 3020 164 REM ************************ CALCULATIONS ************************ 166 IF D$='M' THEN GOSUB 1010 168 IF D$='N' THEN GOSUB 2010 170 IF F0+F1+F2+F3+F6+F7+F8 = 0 THEN GOSUB 2852 172 REM ********* PRINT OUTPUT (AND WARNINGS — IF APPLICABLE) ******** 174 GOSUB 4006 176 GOSUB 4306 178 END 1000 REM ************************************************************** 1002 REM SUBROUTINE: CALCULATION OF TOWER CAPABILITY — MECH. DRAUGHT 1004 REM ************************************************************** 1006 REM ********************** CHECK DATA INPUT ********************** 1008 REM *** SET ERROR FLAG F7 (INPUT DATA ERROR) 1010 F7=0 1012 REM *** SET ERROR FLAG F3 (NO CONVERGENCE) 1014 F3=0 1016 IF W1<3 THEN F7=1 1018 IF W2<3 THEN F7=1 1020 IF T2<=W1 THEN F7=1 1022 IF T6<=W2 THEN F7=1 1024 IF T1<=T2 THEN F7=1 1026 IF T5<=T6 THEN F7=1 1029 IF T1>90 THEN F7=1 1030 IF T5>90 THEN F7=1 1032 IF P$=/A' THEN P1=FNP(A1) 1034 REM * TEST ATMOSPHERIC PRESSURE IS ASSUMED EQUAL TO THE DESIGN VALUE * 1036 REM SEA LEVEL PRESSURE IS USED FOR ALTITUDES UP TO 300 m (BS4485: Pt 3) 1038 IF P1>=FNP(300) THEN IF P1O THEN IF D8<>SGN(K3-K4) THEN I8=.5*I8 1186 J3=J3+SGN(K3-K4)*IU 1188 IF J3<0 THEN J3=.01 1190 I3=I3+1 1192 GOTO 1144 1194 REM *********** CALCULATE PERCENTAGE TOWER CAPABILITY ************ 1196 E5=1OO*J3/J1 1198 RETURN 2000 REM ************************************************************** 2002 REM SUBROUTINE : CALCULATION OF TOWER CAPABILITY - NATURAL DRAUGHT 2004 REM ************************************************************** 2006 REM ********************** CHECK DATA INPUT ********************** 2008 REM ***SET ERROR FLAG F7 (INPUT DATA ERROR) 2010 F7=0 2012 *** SET ERROR FLAGS FO F1, F2, F3 (NO CONVERGENCE) 2014 FO=0 2016 F1=0 2018 F2=0 2020 F3=0 2022 IF W1<3 THEN F7=1 2024 IF W2<3 THEN F7=1 2026 IF T2<=W1 THEN F7=1 2028 IF T6<=W2 THEN F7=1 2030 IF T1<=T2 THEN F7=1 2032 IF T5<=T6 THEN F7=1 2034 IF T1>90 THEN F7=1 2036 IF T5>90 THEN F7=1 2038 IF D10 THEN IF P1=FNP(0) THEN R$='See NOTE (1)' 4058 PRINT ' Atmos. pressure (kPa)',FNB(P1), 'Design value',R$ 4060 IF P$='A' THEN PRINT ' Altitude (m asl) ',A1,'—' 4062 PRINT ' L/G water/air ratio',J1,FNB(J2) 4064 PRINT ' "n" For Kav/L=C*(L/Q)^n',N1,'Design value' 4066 PRINT ' Tower Kav/L ',FNB(K1),FNB(K2) 4068 PRINT 4070 PRINT ' TONER CAPABILITY (%)', '100',FNA(E5) 4072 IF F0+F1+F2+F3+F4+F6+F7+F8>0 THEN 4082 4074 PRINT ' EXPECTED RECOOLED WATER TEMPERATURE (C)',FNA(T6+D5) 4076 IF FNA(D5)<0 THEN PRINT ' WORSE THAN DESIGN BY',FNA(ABS(D5)):’(K)' 4078 IF FNA(D5)>0 THEN PRINT ' BETTER THAN DESIGN BY',FNA(D5):’(K)’ 4080 IF FNA(D5)=0 THEN PRINT ' EQUAL TO DESIGN’ 4082 PRINT 4084 IF R$='' THEN 4092 4086 PRINT ' NOTE (1): For altitudes O~3OO m above sea level, the test data' 4088 PRINT ' is evaluated at mean sea level pressure (BS 4485: Pt 3)' 4090 PRINT 4092 RETURN 4300 REM ************************************************************** 4302 REM SUBROUTINE: PRINT WARNINGS IF APPLICABLE 4304 REM ************************************************************** 4306 IF F0=1 THEN PRINT ' *** WARNING: NO CONVERGENCE (DESIGN EXIT AIR TEMP.)' 4308 IF F1=1 THEN PRINT ' *** WARNING: NO CONVERGENCE (TEST EXIT AIR TEMP.)' 4310 IF F2=1 THEN PRINT ' *** WARNING: NO CONVERGENCE (TEST L/G)' 4312 IF F3=1 THEN PRINT ' *** WARNING: NO CONVERGENCE (L/G FOR TOWER CAPABILITY)' 4314 IF F4=1 THEN PRINT ' *** WARNING: NO CONVERGENCE (EXPECTED COLD WATER TEMP.)' 4316 IF F6=1 THEN PRINT ' *** WARNING: NEGATIVE DRAUGHT' 4318 IF F7=1 THEN PRINT ' *** WARNING: INPUT DATA ERROR' 4320 IF F8=1 THEN PRINT ' *** WARNING: NEGATIVE ENTHALPY DRIVING FORCE' 4322 RETURN 5000 REM *********** DATA FOR FNV(X): A,B.C,D,E,F,G,TO,T9 ************ 5002 DATA —2948.997118,-2.1836674,—.000150474,—.0303738468,.00042873 5004 DATA 4.76955,25.83220018,273.15,273.16 5006 REM *************** DATA FOR FNW(W0): M1,M2,C1 ****************** 5008 DATA 18.01534,28.9645,.000666 5010 REM ***************** DATA FOR FNE(D0); C2,C3,L0 **************** 5012 DATA 1.00568,1.84598,2500.84 5014 REM ****************** DATA FOR FNS(W0): R0 ********************* 5016 DATA 8.31432 5018 REM ****************** DATA FOR FNP(X): P9,B1,B2 ************** 5020 DATA 101.325,-.0118917,4.94444E-07 5022 REM **** DATA FOR N(4): INTERVALS FOR NUMERICAL INTEGRATION ***** 5024 DATA .l,.4,.6,.9 5026 REM ********** SPECIFIC HEAT OF WATER (ReFerence 2); C0 ********** 5028 DATA 4.18684 7000 REM *********** Program P4485: Rev. 6: GENERAL NOTES ************* 7002 REM 7004 REM (A) The program uses constants of up to 10 digits for water vapour 7006 REM pressure correlation vs. temperature. This version of the program 7008 REM is also suitable for single precision BASIC with data storage 7010 REM limited to six digits, but such limitation mag affect the results 7012 REM slightly to an extent which is normallg negligible. 7014 REM 7016 REM (B) The tolerances for iteration loops were tested and are set 7018 REM wide to ensure convergence when using six—digit data storage. 7020 REM 7022 REM (C) The input data contains usual information For the design 7024 REM conditions and the averaged test results» as listed in lines 7026 REM 36-138. The 'n' value For the Fill characteristics also needs 7028 REM to be input: KaV/L=C*(L/G)^n 7030 REM 7032 REM (D) For altitudes From 0 to 300 m above sea level the program 7034 REM evaluates results at ambient pressure equal to that at sea level, 7036 REM in accordance with BS 4485: Part 3. For altitudes over 300 m 7038 REM above sea level the calculations are made at reduced ambient 7040 REM air pressure corresponding to the design altitude. 7042 REM (Reference: BS 2520, Barometer conventions and tables) 7044 REM 7046 REM (E) Four examples oF output are appended For information. 7048 REM 7050 REM (F) An error condition occurs if the limits of this BS Test 7052 REM Code are exceeded, or if there is no convergence after 50 7054 REM iterations. In such cases the flow of calculations stops and 7056 REM partial output is printed, which contains input data and 7058 REM results calculated prior to the occurence of this error. 7060 REM Noughts are displayed for the remaining data. This may help 7062 REM to trace the problem. Several diagnostic messages may appear 7064 REM in the output to warn the user oF possible input data errors, 7066 REM or ome unspecified abnormalitg which prevents the completion 7068 REM of calculations 7070 REM 7072 REM (G) The program may need adjustments to suit user's operating 7074 REM system, or the type of BASIC in use. An output file may be 7076 REM written by simple modiFication oF the subroutine 4000. 7078 REM 7080 REM ****** End of program P4485 Rev. 6 Date: 30 January 1987 ******