CT 2» 1371
WATER POLLUTION CONTROL RESEARCH SERIES • 16130EES11/70
RESEARCH ON
DRY - TYPE COOLING TOWERS
FOR THERMAL ELECTRIC
GENERATION
Part II
ENVIRONMENTAL PROTECTION AGENCY • WATER QUALITY OFFICE
-------
WATER POLLUTION CONTROL RESEARCH SERIES
The Water Pollution Control Research Series describes the
results and progress in the control and abatement of pollu-
tion of our Nation's waters. They provide a central source
of information on the research, development, and demon-
stration activities of the Water Quality Office, Environ-
mental Protection Agency, through inhouse research and grants
and contracts with Federal, State, and local agencies, re-
search institutions, and industrial organizations.
Inquiries pertaining to the Water Pollution Control Research
Reports should be directed to the Head, Project Reports
System, Office of Research and Development, Water Quality
Office, Environmental Protection Agency, Washington, B.C. 20242.
-------
RESEARCH ON DRY-TYPE COOLING TOWERS FOR
THERMAL ELECTRIC GENERATION: PART II
Computer Program Descriptions
by
R. W. Beck and Associates
600 Western Federal Savings Bldg.
Denver, Colorado 80202
for the
WATER QUALITY OFFICE
ENVIRONMENTAL PROTECTION AGENCY
Project # 16130 EES
Contract # 14-12-823
November 1970
For sale by the Superintendent of Documents, U.S. Government Printing office, Washington, D.C. 20402 - Price $1
-------
EPA Review Notice
This report has been reviewed by the Water Quality Office,
EPA, and approved for publication. Approval does not signi-
fy that the contents necessarily reflect the views and poli-
cies of the Environmental Protection Agency, nor does mention
of trade names or commercial products constitute endorsement
or recommendation for use.
-------
TABLE OF CONTENTS
Page
TITLE PAGE
TABLE OF CONTENTS '
GENERAL 1
ECONOMIC OPTIMIZATION PROGRAM "OPTDCT" 1
Description of Program "OPTDCT" 2
General Description of Subroutine "DATCHK" 4
General Description of Subroutine "TBLUQ" 5
Description of Data File "TURBIN" 5
Description of Data File "SITEXX" 6
Format of External Data Files 7
Output Data File 7
Operating Instruction for "OPTDCT" 7
"OPTDCT" Flow Chart Description 10
TOWER OPTIMIZATION PROGRAM "TOWSIZ" 18
Description of Program "TOWSIZ" 18
Description of Subroutine "TBLUQ" 19
General Description of Subroutine "TWRPIP" 19
General Description of Subroutine "PIPSIZ" 19
Description of Data File "SIZDAT" 20
Output Data 20
Operating Instruction for "TOWSIZ" 20
"TOWSIZ" Flow Chart Description 20
Subroutine "TBLUQ" Flow Chart Description 22
Subroutine "TWRPIP" Flow Chart Description 23
Subroutine "PIPSIZ" Flow Chart Description 24
APPENDIX A - Program Listings 25
APPENDIX B - Sample Data 42
APPENDIX C - Output Samples 49
APPENDIX D - Flow Charts 53
APPENDIX E - Glossary 83
-------
DRY COOLING TOWER OPTIMIZATION
COMPUTER PROGRAMS
GENERAL
The study of dry-type cooling towers as described in the report Research on
Dry-Type Cooling Towers for Thermal Electric Generation was facilitated by the
development of two computer programs to aid in the analysis of the large quantity
of data. The writeup herein describes the two programs in detail and provides in-
structions for their operation.
The computer used was a CDC 6400 in the timesharing system of United
Computing Systems, Inc., of Kansas City, Missouri. Input was by paper tape on a
remote teletype terminal. Remote and local batch output by terminal was possible;
however, remote batch was primarily used due to the large amount of output infor-
mation. The programs were coded in FORTRAN IV.
The physical dimensions of a natural-draft dry cooling tower and its capital
cost are evaluated by the program titled "TOWSIZ". The economically optimum
dry cooling system is determined by the program titled "OPTDCT".
ECONOMIC OPTIMIZATION PROGRAM "OPTDCT"
The purpose of the program OPTDCT is to determine, for a given set of con-
ditions at a particular location, the economically optimum dry-type cooling system
for a large thermal-electric generating plant.
The size of the dry cooling system is a function of the initial temperature
difference (ITD) which is the difference between the turbine exhaust steam tempera-
ture and the ambient air temperature. The economically optimum dry cooling
system for a specific set of conditions is that which results in the lowest total annual
cost.
The annual costs which are affected by the ITD of the dry cooling system
and which, therefore, must be considered in the economic optimization analysisare
the annual capital costs, operation and maintenance costs, total generating plant
fuel costs, auxiliary power costs and the cost of replacing generating capacity lost
at high ambient air temperatures.
The economic optimization program OPTDCT requires the input from two
external data files, TURBIN and SITEXX. All other data required for the analysis
are either in data statements within the program or are created by the program and
then stored in arrays.
-------
The data file TURBIN contains the following data pertaining to the turbine:
nameplate rating; type of fuel (fossil or nuclear); the design heat rejection; limiting
values for the operating back pressure; the turbine output for back pressure from 1 .0
inch Hg to 18.0 inches Hg in steps of 0.5 inches; and corresponding values of tur-
bine heat rate, station heat rate and heat rejection for the three turbine operating
conditions of full throttle, 3/4 load and 1/2 load.
The data file SITEXX contains the following site related information: site
name, elevation, reason of peak electrical demand, fuel costs, fixed-charge rates,
a construction cost multiplier, the capital cost per kw for peaking units, a cutoff
temperature used in determining when the peaking unit is operated, the capital cost
per kw for auxiliary power, operation and maintenance cost data and ambient air
temperature duration data for the site.
Description of Program "OPTDCT"
The following is a detailed section by section explanation of the computer
program. A glossary of terms, used in the program, can be found in Appendix E.
The first section of the program, ending at line 01080, contains declaration,
data and comment statements. The data within this section does not change with
the exception of the following: AUX100 data line no. 00360, KWAX data line no.
00440, CAPCST data lines 00960-01020, and the TOWER data line no. 01060.
These blocks of data have to correspond to the tower and turbine being analyzed.
A further explanation of the above mentioned data is in the section on operating
instructions.
The next section, line 01100 to line 01260, is the section that creates the
air temperature and back pressure arrays. The first part creates the 32 values of air
temperature starting at 117°F and decreasing by 5° increments to —38°F. The
second part creates the back pressure data. These values start at 1 .0 inch Hg and
go to 18.0 inches Hg by 0.5 inch increments.
In the section beginning on line 01260 and ending on line 01580 the turbine
data from file TURBIN is read into the program. Also within this section the nom-
inal generation for the given set of operating conditions is calculated. This is done
in the statements on lines 01320 and 01340.
After the first eight lines of data have been read into the program, the data
check subroutine (DATCHK) is called and the data is checked to see if the right
number of values have been read. The next fifteen lines of data are read and the
subroutine is called again. This is done until all of the turbine data has been read
into the program.
-------
The next section of the program reads in the site data. An explanation of
how the site data files are set up is contained in the section on data files. After the
site data is read, subroutine DATCHK is called again to verify that the correct num-
ber of values were read into the program.
Also within this section, an indicator (IWSP) is set up to direct the computer
to the different calculations which are involved with the different peaking periods.
This section also finds the air temperatures corresponding to the first and last non-
zero temperature durations.
The site identification data is written to the output file in the next section.
This starts at line 02520 and ends at line 03040.
The section starting with line 03060 sets up variables that are used in calcu-
lating operating costs. These variables are used later in the program, mostly in
table look-up calculations.
The next section of the program deals with capital cost. The computer picks
the capital cost for a specific ITD and elevation of the site being analyzed from a
table of capital costs.
The capital cost is then multiplied by the construction cost index for the site.
The next section, line 03580 through line 03820, initializes variables that
are accumulative. These variables are cumulative for each different ITD.
Starting at line 03900, the coefficient A is determined in the equation
ITD = A*HREJ**Z. The value of the exponent Z used in the program was .75 for a
natural-draft tower and .91 for a mechanical-draft tower.
The calculations to find the operating back pressure for the different ait-
temperatures and for the different turbine loadings begin at line 04000. The back
pressures are found by the method of successive approximations and are calculated
from line 04120 to line 04660. A limitation has been set on the number iterations
it takes to arrive at the turbine operating back pressure. If the limitation is ex-
ceeded then the operating back pressure is set equal to the allowable maximum. If
the back pressure is found before the limitation is reached, the computer then
checks to see if it is within the designated limit. If the back pressure is within the
maximum and minimum limitations, then that pressure is used to determine the power
of the turbine and also its heat rate ,
If the back pressure is not within the limiting values, the turbine operation
is adjusted to operate within the designated limits. These calculations are made
between lines 04680 and 04780.
-------
In the section beginning at line 04800 the annual fuel consumption (Btu) is
calculated. This value is calculated from the operating station heat rate and re-
flects the ambient air temperatures during the year and the generating plant loading
assumptions.
Starting at line 05000, the auxiliary power requirements for the varying
ITD's and air temperatures are calculated. The program is directed to the appro-
priate calculations depending upon the type of tower (natural draft or mechanical
draft). Auxiliary power and energy requirements for a mechanical-draft tower are
calculated between lines 05020 and 05440. Natural draft auxiliary power and
energy requirements are calculated between lines 05480 and 05540.
The next section beginning at line 05560 is the section in which the loss or
gain of capacity is calculated. If the operating back pressure is greater than 3-1/2
inches Hg, then there is a loss of capacity and the amount is calculated in the part
beginning at line 05760. If there is a gain, the amount is calculated in the part
beginning at line 05620.
The loss of capacity occurring at the temperature equalled or exceeded 10
hours per year is calculated on line 05960.
The next section of the program, beginning at line 06080, writes into the
output file the information identifying each run.
The section beginning at line 06240 writes into the output file the column
headings for the program results.
The last section of the program, starting at line 06560, calculates the total
annual cost. This annual cost data is then searched by the computer to find the
minimum annual cost and therefore the optimum ITD. After the optimum ITD is
found the results are placed in the output file.
The remainder of the program consists of two subroutines used by the main
program. These subroutines, DATCHK and TBLUQ, are described in the next two
sections.
General Description of Subroutine "DATCHK"
The purpose of this subroutine is to check if the correct number of values of
data have been read into the main program. This check is accomplished by testing
to see if the data check number, 1 .OE50, is in its proper place within the data file.
Data check numbers can be found on lines 00180, 00340, 00500, and 00660
in the turbine data file TURBIN and on line 250 in the site data file. After the
data preceding the data check number is read into the program, subroutine DATCHK
-------
is called and it tests to see if 1 .OE50 is in its proper place. If there are some data
missing or if there are extra data within the file, the subroutine will read some other
number than 1 .OE50. Thus, not finding the data check number it prints out an error
message, indicating that the data has an incorrect number of values. The message
printed out by the computer is the name of the data file it is reading and also the
line of data that it read instead of the data check number.
After the message is printed out, the data file is reread by the subroutine to
find the data check number. After finding this number, the data file has been
placed in the proper position so that the remaining data can be read without the
error being carried through the file.
General Description of Subroutine "TBLUQ"
This routine is used to interpolate values from tables created by the program
or tables entered in data statements. It uses the method of determinant solution of
two second order equations.
The first variable in the call is a reference to the first column of the table,
which is specified as the second argument. The subroutine first determines whether
the requested value lies within the limits of the specified table column. !f it does,
processing proceeds to the table search.
If the specified value is outside the column entry the program determines
whether it is above or below the table and sets indicators to extrapolate from the
appropriate end of the table. A return error flag is set and control is transferred to
the interpolation calculation.
The table search looks for successive elements of the reference column whose
values lie on either side of the requested value. Should an entry in the reference
column exactly equal the requested value, the corresponding entry from the look-up
column, the third argument, is returned as the subroutine value. Otherwise the
table search finds two bracketing values and checks whether they are the firstvalues
in the table. If they are, they are used in conjunction with the third variable for
interpolation. If they are not the first values, the bracketing values together with
the table entry preceding them are used for interpolation.
Description of Data File "TURBIN"
This external input file contains the data pertaining to the turbine perform-
ance for varying back pressures for program OPTDCT. Sample data for the fossil-
fueled and nuclear-fueled turbines are in Tables 5-B and 6-B of Appendix B.
-------
The first section of the data file TURBIN, starts with the nameplate capacity
of the generator (kw) and the type of fuel, specified as either FOSSIL or NUCLEAR.
Line 00110 is the nominal heat rejection of the turbine in thousands of Btu per hour.
Line 00120 contains the minimum and maximum values of back pressure (inches Hg)
at which the turbine is allowed to operate.
The next section, starting with line 00130, contains the 35 values of full
throttle generation (kw). These values correspond to back pressures from 1 .0 inch
Hg to 18.0 inches Hg in 0.5 inch increments. This set of data is followed by the
data check number on line 00180.
The next three sections, starting with lines 00190, 00350, and 00510 con-
stitute the main body of the file. These sections correspond to the turbine operating
conditions of full throttle, 3/4 load and 1/2 load, respectively. Each section ends
with the data check number 1 .OE50. This enables the program to check for the cor-
rect number of values.
Each section contains three sets, 35 values each of turbine heat rates (Btu/
kwh), plant heat rates (Btu/kwh), and heat rejection (10° Btu/hr) at various back
pressures. The back pressures correspond to the array of back pressures set up in the
main program.
Description of Data File "SITEXX" (1)
The SITEXX data file contains the data for each location that the economic
optimization program analyzes. Each file consists of 16 lines of data, the last line
of data being the data check number. A sample of a site data file is shown in Table
7-B of Appendix B.
The first three lines of data include the site name, elevation and the season
of peak electrical demands. The name of the site has been limited to thirty charac-
ters. The elevation is in feet above sea levei and the period of peak demand is
either summer or winter.
The next line of data, line 130, includes the number of base plant fuel costs
to be used in the analysis followed by the fuel costs in cents per million Btu. The
program is presently set up to take a maximum of three such fuel costs, but this
number can be increased by changing the array dimension.
Line 140 includes the number of fixed-charge rates to be used in the anal-
ysis followed by the fixed-charge rates. This line of data has been set for a maxi-
mum of five fixed-charge rates. It can be increased in the same manner as the
number of fuel costs. These values are entered as percentages.
(1) Where SITEXX is a file name of SITE01, SITE02, through SITE 27.
-------
The next two lines, 150 and 160, are the capital cost multiplier and the
capital cost of the peaking unit in dollars per kw respectively.
Line 170 contains the fixed-charge rates applicable to the peaking unit.
The same number of values are entered on this line as are entered on line 140.
These values are also entered as percentages.
The next four lines of data contain the peaking unit fuel costs in cents per
million Btu (line 180), the cutoff temperature for peaking generation in °F(line 190),
the capital cost of providing auxiliary power requirements in dollars per kw (line
200) and the operating and maintenance cost as a percentage of the capital cost (line
210). The number of peaking unit fuel costs corresponds to the number of base plant
fuel costs.
The next three lines of data (220, 230, 240) contain the 32 values of tem-
perature durations. These values are the percent of time during a year that the
corresponding air temperatures are expected to occur. The air temperatures corres-
pond to the array of air temperatures which are created in the main program.
The last line (250) contains the data check number 1 .OE50.
Format of External Data Files
Data input files TURBIN and SITEXX are desequenced files. These files can-
not have line numbers. The program reads each line starting with the first character.
The data files shown in Appendix B have line numbers. These are used only to enter
the data files in the computer and are used, in this case, only as reference numbers.
Once the files are entered in the computer they are desequenced.
Output Data File
Output file OPTOUT will contain the results of the computer runs. This file
must be created before the runs are made. To obtain the results from this file, list
file OPTOUT.
Operating Instructions for "OPTDCT"
Before running the computer program there are several items that need to be
checked to see if the data corresponds to the turbine and cooling system being ana-
lyzed. The items that need to be checked are:
1 . Auxiliary data
2 . Capital cost data
3. Type of tower
4. Turbine data
5. Tower characteristics equations
6. Capital cost multipliers
7. Capital cost of auxiliary power requirements
-------
The auxiliary power requirements data should correspond to the type of cool-
ing tower (natural or mechanical draft) and the type of fuel (fossil or nuclear). The
auxiliary power requirements for each of the towers and turbine units are shown in
Tables 1-B and 2-B of Appendix B. These values correspond to ITD's of 30°F through
80°F in 10° increments. The values of auxiliary power requirements for a mechani-
cal-draft tower are entered on line 00360 and the values for a natural-draft tower
are entered on line 00440.
The mechanical-draft cooling tower data are entered in the following form
using data for a fossil-fueled unit:
DATA AUX100/32000. ,26000. ,21000., 17000., 14000., 12000./
The natural-draft cooling tower data are entered in the following form, again
using data for a fossil-fueled unit:
DATA KWAX/12000.,9000.,7000.,6000.,5000.,4000./
The auxiliary data (AUXTMP) on lines 00380, 00400, and 00420 do not
change. The values in the program listing (Appendix A) are generalized due to their
proprietary nature. Contact the manufacturer for specific information. These values
represent the air temperatures associated with the minimum (53%) and maximum
(100%) proportion of full auxiliary power requirements corresponding to ITD values
of 30°F to 80 F in 10 increments for a mechanical-draft tower. They also corres-
pond to full throttle, 3/4 load and 1/2 load operations of the turbine.
The capital cost data, shown in Tables 1-B and 2-B of Appendix B, corres-
pond to the type of turbine and also to the type of cooling system. Depending upon
the type of tower and turbine being analyzed the corresponding capital cost must be
used. The data are entered on lines 00960, 00980, 01000, and 01020 in the follow-
ing form:
DATA CAPCST/32610000. ,22880000., 17200000., 14090000.,
+12040000., 10280000., 36150000.,24430000.,18220000.,
+14700000.,12500000., 10600000., 40500000.,2 6850000.,
+19770000., 15740000., 13260000.,11100000./
The data for sea level elevation is entered first, followed by the data for
elevation of 3,000 feet and 6,000 feet. These data represent three curves of capital
cost versus ITD for the three elevations.
-------
The type of tower being analyzed is entered on line 01060 in the following
forms:
DATATOWER/8HNATURALb,8HDRAFTbbb/
or
DATATOWER/8HMECHANIC,8HALbDRAFT/
These data are used later in the program in the output as well as in directing
the program to the auxiliary calculations that are associated with the type of tower
being analyzed.
The data file, TURBIN, containing the turbine data, corresponds to the type
of turbine that is being analyzed, either a fossil- or nuclear-fueled unit. These
data, along with the tower data determine the type of analysis that will be made.
The remaining data used in the calculations corresponds to these two sets of data.
The last item to be checked in the main program is the exponent in the tower
operating characteristics equations. The exponent must correspond to the type of
tower being used. The value of the exponent is 0.75 for natural draft and 0.91 for
mechanical draft. The three equations involving the exponent change are on lines
03680, 04140, and 04760. These equations are as follows for a natural-draft tower:
03680 A = THETA/HREJD**.75
04140 110 SATT = AIRT(NT) + A*HRJ**.75
04760 HREJMX = (PLITD*A)**(1 ./.75)
The last two items to be checked are in the site data. The first is the capital
cost multiplier and the second is the auxiliary capital cost.
The capital cost multiplier is to allow for variations in construction costs due
to local conditions. In areas subject to hurricane winds the construction cost will be
greater due to the requirement for stronger structures to resist the forces of these
winds. Therefore, a capital cost multiplier is used to provide for the resulting in-
crease in costs. The adjustment for hurricane winds is applied only to the analysis
of a natural-draft tower. The mechanical-draft towers are not sufficiently affected
by the high velocity winds to warrant an increase in the multiplier above the normal
construction cost index for that particular area.
The last item to be checked is the capital cost of auxiliary power. This
value depends on the type of unit (fossil- or nuclear-fueled). The capital cost of
auxiliary power is assumed to be $150/kw for a fossil-fueled unit and is assumed to
be $225/kw for a nuclear-fueled unit.
-------
After checking the above listed items, the program is ready to be run.
Each time that a different tower or turbine is used, the above checks must be re-
peated to insure that the correct data is in the program.
"OPTDCT" Flow Chart Description (1)
1 . Dimension variables for tables, temporary storage, etc. Define constants
and tables which are not dependent on the site or type of turbine.
Constants: Tables:
OPHRS AUXlOOvs. TBLITD
HRSPYR AUXTMP vs. TBLITD vs. PCTLOD
ASTER KWAX vs. TBLITD
TOWER P vs. T
SITE NO PCTTIMvs. PCTLD
CAPCST vs. TBLITD vs. TBLELV
Define ambient air temperatures for table of ambient air temperature vs.
annual duration, in percent, of that temperature (^2 F). AIRT ranges
from 117°F down to -38°F in steps of 5°F.
Define turbine exhaust pressures to associate with values at turbine heat
rate, plant heat rate, heat rejected by turbine, and full throttle kw which
are read from the turbine file. BP (back pressure) values range from 1 .0
inch Hg to 18.0 inches Hg in steps of 0.5 inch Hg.
2. Read the characteristics of the turbine used. This includes the rated kw,
whether fossil or nuclear, the design heat rejection and the minimum and
maximum back pressure allowable for operating the turbine. Also included
are values for tables of back pressure vs. full throttle kw, and values for
full throttle, 75% load and 50% load for tables of back pressure vs. net
heat rate, gross heat rate and heat rejection.
3. In several specific places the number 1 .OE50 must be put into the data file.
A subroutine DATCHK is called to check that these numbers are in the
correct place in the data file. If they are not in the correct place, there
are too many or too few values in the data file. In this case the subroutine
prints an error message indicating the data file and the location in the file
where the error occurred.
(1) Numbers on these pages correspond to numbers
on the flow chart blocks in Appendix D.
10
-------
4. The nominal annual mwh is the amount which could be produced if
operated at the given conditions of load and amount of time specified for
that load.
5. Start a loop to cycle through each of the different sites printing a complete
set of output for each .
6. Read the information from the appropriate site file. This should include all
information that might change for different locations. Included is the name
of the site, the elevation, whether winter or summer peak, the fuel costs,
fixed-charge rates, capital cost multiplier, cutoff temperature for peak
generation, operation and maintenance percentage, and the duration in
percent of a year associated with each of the temperatures in AIRT.
7. Set an indicator for applicable winter or summer peaking period.
8. Find the first and last non-zero temperature durations and set a pointer to
each so that later processing can be confined to meaningful values.
9. Call the subroutine DATCHK while reading the site file to verify that the
right number of values are on the file.
10. Output information which identifies the site being processed, the type of
turbine, the type of dry cooling tower, the turbine operating hours per
year and the three combinations of percent time at a given percent load
which describe the annual demand on the turbine. All output except for
error message is done indirectly; i .e., it is written onto a data file during
the execution of the program and then the data file is tested after the pro-
gram is finished.
11 .12. Compute several values for later use in interpolation of tables through the
13.14. use of subroutine TBLUQ. These include values of kw load for the maxi-
mum operating back pressure, kw load for the actual operating back pres-
sure, heat rejected by the turbine at the maximum operating back pressure,
the hours duration of the three highest temperatures, and the three highest
temperatures at that site. TBLUQ is a function subprogram that, given a
value in one column of a table, returns the corresponding value from the
second column of the table. If the given value is between values in the
table, a quadratic interpolation is done to arrive at the value to be re-
turned.
15. Find the air temperature which is equalled or exceeded 10 hours per year.
Use TBLUQ and the three valued table just computed of DUR vs. AIR.
11
-------
16. Set up a loop to perform the basic calculations with each ITD from 30°F to
80°F in 1 increments.
17. Convert the integer ITD into an equivalent real value THETA. The integer
I is the subscript which corresponds to each ITD. I ranges in value from 1
to 51 in increments of 1 .
18. Do two sets of table look-ups using TBLUQ to determine the capital cost of
the cooling system which is a function of both elevation and ITD. This
table is different for natural- and mechanical-draft towers and is deter-
mined outside of this program.
19. Multiply capital cost by capital cost multiplier.
20. The coefficient A must be calculated for each ITD. The given ITD and
design heat rejection are used to calculate A. Z is a given exponent and
is 0.75 for natural-draft towers and 0.91 for mechanical-draft towers.
Once A is computed this equation will give the ITD of a given tower for
all heat rejection values.
21 . Set up loop to cycle through each of the three loading conditions of per-
cent time vs. percent load.
22. Set up loop to cycle through each of the air temperatures between the
highest and lowest recorded at the given site.
23. Initialize iteration counter and two variables used for back pressures.
OPBP is initialized at BPMIN (minimum back pressure) because if the
operating BP is less than BPMIN it is possible to skip out of the loop be-
fore OPBT is computed. Skip to block 29 which calculates the heat re-
jected from the turbine.
24.- This group of blocks involves an iterative procedure to arrive at the actual
32 . operating back pressure. Essentially what must be found is the intersection
of the turbine characteristic curve for the given load and the cooling tower
characteristic curve for the given air temperature. Both the turbine and
cooling tower curves are in terms of heat rejection vs. back pressure. The
turbine curve is arrived at from a table of values from the turbine data file
and the cooling tower curve is arrived at by using the equation Part Load
ITD = A x (Heat Rejection) , the ambient air temperature and the rela-
tion between saturation temperature and saturation pressure.
The solution is arrived at by first assuming a back pressure (the minimum
back pressure is used for a starting value). For this back pressure and the
given load the turbine will produce a certain amount of heat rejection.
12
-------
This heat rejection is used to calculate the part load ITD of the cooling
system. Adding the ambient air temperature and the ITD gives the satu-
rated steam temperature . This temperature corresponds to a saturated steam
pressure which is exhaust back pressure of the turbine. If the preceding
back pressure and the newly derived back pressure are within 0.01 inch Hg
of each other then that is the back pressure which is used in the analysis.
If the two values differ by more than 0.01 inch Hg then the newly calcu-
lated back pressure is used as the starting point of the calculation . If after
15 iterations a solution has not been found, an error message is printed, the
last back pressure calculated is used and processing continues.
33.34. If the calculated back pressure is less than the minimum, make it equal to
the minimum .
35. In order to calculate energy requirements, it is necessary to calculate the
total number of hours of the year that the plant runs at the given load and
air temperature .
36. Check to see if the calculated back pressure is greater than the maximum
allowed for the turbine and if so, the turbine output must be decreased to
an allowable level . Transfer to block 39 to do this.
37. For the computer back pressure, find the maximum possible output in kw of
the turbine .
38. Compare the maximum possible output with that required and use the
smaller as the actual turbine output. Skip to block 41 .
39. If the calculated back pressure exceeds the maximum, then plant output is
limited by the amount of heat the cooling system can reject at the maxi-
mum back pressure. Calculate this heat rejection by using the maximum
back pressure to arrive at a corresponding saturation temperature; this tem-
perature minus the air temperature gives the maximum part load ITD which
is possible; from this part load ITD the corresponding heat rejection can be
determined using the equation:
Heat Rejection = Part Lo°d
= I
\
40. By doing a table look-up (with TBLUQ) in a table of heat rejection vs.
load at the maximum back pressure and using the calculated heat rejec-
tion as one value in the table, the maximum turbine output can be found.
13
-------
41 . The station heat rate may now be found from the table of station heat rate
vs. back pressure vs. load, by four calls to TBLUQ- three calls to build a
table of load vs. station heat rate for a constant back pressure and another
call to find the station heat rate for given load.
42. Calculate and accumulate the total annual fuel requirement (Btu). The Btu
requirements for each load and temperature combination are calculated
from the station heat rate multiplied by the total hours for that condition,
and then added to the previous subtotal.
43. The next few blocks calculate the auxiliary power and energy requirements.
These steps differ depending on whether the tower is natural or mechanical
draft. This block tests to see if the tov/er is natural draft and skips to 47 if
it is.
44. The maximum auxiliary power requirement for each ITD is found by inter-
polating a table of auxiliary power requirements vs. ITD
45. A linear interpolation of another table gives the percent of maximum auxil-
iary power requirement used for a given air temperature, ITD and load.
46. Compute the annual auxiliary kwh for mechanical draft by applying the
above percentage to the product of the maximum auxiliary power require-
ment computed in block 44 and the number of hours per year at this condi-
tion. Skip to block 49.
47. Determine the auxiliary power requirements for the natural-draft tower by
a call to TBLUQ using the ITD given and the table of ITD vs. KWAX.
48. Compute the auxiliary energy requirements by multiplying the auxiliary
power requirements by the number of hours at the given conditions.
49. Calculate and accumulate the total energy produced by the turbine. The
energy produced at the given conditions is computed by multiplying the kw
produced by the turbine times the number of hours at the given conditions.
50. If, at full throttle, the amount produced by the turbine is different than
the nameplate capacity, then a gain or loss of capacity and energy must
be computed.
51 .55. If there is a net gain, calculate that gain in kwh and Btu. The kwh gain
is calculated by taking the difference between the actual amount produced
and the amount required, times the number of hours at the given conditions.
The associated Btu requirement is computed by taking the kwh gain times
14
-------
the calculated station heat rate. A gain may only occur if the demand is
100% of nameplate and the turbine can produce more than that amount.
Skip to block 56.
52. If there is a loss, check whether the air temperature is above or below the
cutoff temperature.
53. If below the cutoff temperature, calculate the energy (kwh x hrs) and store
in KWHBLO. Skip to block 57.
54. If above the cutoff temperature, calculate the energy and store in KWHABV.
56. If processing one of the three highest temperatures, save the amount of
power loss or gain for later interpolation.
57.58. Go to next lower temperature and its duration, until all non-zero durations
have been used, then go back to block 23.
59.60. After making the calculations for all temperatures for a given load go to
the net load condition until calculations have been made for all load con-
ditions.
61 . The total energy produced in a year has now been calculated for the plant.
Save this value for each ITD.
62. Determine the capacity loss which occurs at the air temperature equalled
or exceeded 10 hours per year by interpolation of a table of the three
highest air temperatures and the corresponding kw losses or gains for those
temperatures.
63.64. All values which vary with ITD, for a given site and ITD, have now been
calculated and stored in arrays. Calculate the same values for the next
higher ITD until values have been computed for all ITD values up to and
including 80 F.
65. Set up counter to calculate annual costs using each of the fuel costs. Up
to three fuel costs may be used for each site.
66. Set up another counter to calculate the. annual costs for each of the five
different fixed-charge rates.
67. Write the headings on the output file for the page of outputs. Along with
the column headings, the capital cost factors, the fuel costs, and the
capital costs per kw for peaking and auxiliary power are outputted in order
to identify the case being run.
15
-------
68.69. Set up loop to calculate and output annual costs and intermediate results
for each ITD.
70. The annual cost which is associated with the capital cost of the cooling
system is equal to the capital cost multiplied by the fixed-charge rate.
71 . Calculate the operation and maintenance cost of the cooling system (which
is a percentage of the total capital cost) and add this cost to the total
annual cost.
72. Compute the total annual cost of fuel by taking total annual Btu times the
fuel cost per million Btu.
73.74. Check whether peak is in winter or summer. For a summer peak compute
75. the capital cost of replacing the lost capacity and a cost for the energy
lost. For a winter peak, compute only the penalty for the lost energy since
the peak system demand would not occur during the months when the capa-
city losses would occur. The capital cost of replacement capacity is com-
puted as the capacity loss which occurs at the air temperature equalled or
exceeded 10 hours per year times a cost in $/kw for peaking power to re-
place it, plus an operation and maintenance cost assumed as a $1 .20 per
kw loss.
The operating cost is determined by multiplying the energy lost by the
appropriate heat rate and the appropriate fuel cost. This cost is divided
into two parts, operating cost associated with energy lost above the cutoff
temperature and operating cost associated with energy lost below the cut-
off temperature. Above the cutoff temperature the lost energy was re-
placed by a peaking unit; thus a higher heat rate and a higher fuel cost
was used. Below the cutoff temperature the lost energy was assumed to be
replaced by another large base load unit; thus the base load unit heat rate
and fuel cost were used.
76. Add the capital and, operating costs for the penalty to get the total penalty
cost.
77. Add the kwh below the cutoff temperature to the kwh above the cutoff
temperature to get the total energy loss.
78. Capital cost of auxiliary power is computed by taking the maximum auxil-
iary power required times the $/kw cost of auxiliary power. The operating
cost of the auxiliaries is the auxiliary energy in kwh times the average
plant fuel cost in mills/kwh plus an operation and maintenance cost. The
total annual cost of the auxiliary power and energy is the sum of the
capital cost times a fixed-charge rate plus the operating cost.
16
-------
79. The annual cost of the cooling system is the summation of: annual capital
cost of the cooling system including an operation and maintenance charge,
the annual cost of the total plant fuel, the annual cost of replacing capa-
city and energy losses due to high turbine back pressure, the annual cost of
auxiliary power and energy, and the credit given for the generation of ex-
cess energy.
78.79. The optimum total annual cost is to be flagged in each case. Therefore, a
test must be made to save the lowest total annual cost and test with the
next one calculated. The ITD of the optimum value is also saved.
82. For each ITD the following values are saved on a temporary file for later
output:
a. ITD
b. Actual mwh produced
c. Credit energy in mwh
d. Penalty energy in mwh
e. Auxiliary energy in mwh
f. Capacity lost in kw
g. Auxiliary power in kw
h. Annual capital and operation and maintenance
cost of the dry cooling system
i . Annual cost of total plant fuel
j. Amount credited due to excess energy
k. Cost of replacing capacity and energy cost
I. Cost of auxiliaries
m. Total annual dollar cost
n. Total annual cost in mills/kwh
81 .82. Go back to block 68 and process with next ITD until ITD = 80°F; then
continue on .
85. Copy the temporary data file onto the permanent output file and look for
the ITD with the minimum annual cost. Before copying the line of data
associated with the optimum ITD, write the word "OPTIMUM" to indicate
that the next line is the optimum condition.
86.87. Calculate a whole new set of output values with a different fixed-charge
rate until all five fixed-charge rates have been used.
88.89. Calculate another set of five pages of output using the second set of fuel
costs, and the final set using the last set of fuel costs. This means that
there will be fifteen combinations of fixed-charge rates and fuel costs for
each site.
17
-------
90.91 . Now that all information for one site has been written onto the output file
go back and do the same thing with the next site until all 27 sites have
been processed. Once the program execution is finished the output file is
listed to obtain a printed copy of the printout.
TOWER OPTIMIZATION PROGRAM "TOWSIZ"
Program TOWSIZ computes the physical dimensions of a natural-draft dry-
type cooling tower and evaluates the capital cost of constructing the tower.
Six basic design parameters are required by the program as input. These
values are the ITD, range, heat rejection, water flow, ambient air temperature and
elevation. Using these parameters, the program evaluates a number of intermediate
variables which contribute to the computation of the tower size (height, top dia-
meter and bottom diameter).
After sizing the tower the program evaluates the capital cost of the cooling
system. The capital cost is divided into the areas of tower structure, condenser,
piping and valves and controls. The tower structure cost is the make-up of stack,
shed and coil costs. Piping and valves includes the cost for piping, valves, circu-
lating water pumps, cooling system filler pump and storage tank.
Within the subheading of structure costs, the dollar value of the tower it-
self is found by interpolation between curves for three evaluated tower diameters.
The shed cost is evaluated from the area difference between the base diameter and
top diameter charged at $7.60 per square foot. Coil cost is evaluated at $14,500
per coil.
Condenser cost is arrived at from a table of cost vs. heat rejection.
The costs of piping, valves, pumps, and storage tanks are computed in sub-
routine TWRPIP. They are developed by computing the water flow needed to reject
the amount of heat produced by the plant at its design point. From the GPM re-
quired, the pipe and valve sizes and quantity of each are determined. The number
of pumps to handle the flow is figured and the required storage capacity is calcu-
lated. From cost tables in the subroutine, the costs of the above items are figured
and then summed. $500,000 is included for controls.
Description of Program "TOWSIZ"
The first section of the program ending at line 02380 contains declaration,
data and comments statements. The data within this section do not change.
18
-------
Beginning at line 02400, the intermediate results for determining the tower
size is computed. After the intermediate results are calculated, the actual physical
size of the tower is computed. Results are the tower height, upper diameter, bottom
diameter and the water flow rate.
The next section beginning at line 03340 calculates the capital cost of the
cooling system. The stack, shed, and coil costs are calculated. These constitute
the structure costs. Then the condenser cost is calculated beginning at line 03540.
Line 03620 is the statement which calls routine TWRPIP. This routine cal-
culates the costs of pipes, valves, circulating water pumps, filler pumps and storage
tank costs.
The section beginning on line 03640 and ending on line 04140 prints out the
results of the tower sizing program. A sample of the printout is shown in AppendixC.
Each line of data in data file SIZDAT will result in a size and cost analysis of a
cooling system for the specified data.
Description of Subroutine "TBLUQ"
This subroutine is identical to the one in the economic optimization pro-
gram. See prior description beginning on page 8.
General Description of Subroutine "TWRPIP"
This routine is used to calculate the cost of piping, valves, pumps, and
storage facilities. This routine uses subroutine PIPSIZ to calculate the size of pipes
and valves.
Once the sizes of the pipe and valves are determined, subroutine TWRPIP
calculates the quantity of each that is required and their costs. After calculating
all the piping and valve costs, the pumping requirements are calculated. After
arriving at the pumping costs, the cost of water storage facilities are calculated.
These costs of piping, valves, pumps and storage are then summed and
added to the other cooling system costs.
General Description of Subroutine "PIPSIZ"
This routine calculates the size of the pipe required for a given water flow
rate. After the pipe size is calculated, the corresponding pipe and valve costs are
determined from a table of pipe and valve costs.
19
-------
Description of Data File "SIZDAT"
This external data file is used by program TOWSIZ. A sample is shown in
Appendix B. There are six values on each line. The first value is the initial tem-
perature difference followed by water cooling range, heat rejection, water flow,
ambient air temperature and elevation. Each line of the file has corresponding data.
There is no maximum limit set on the number of lines of data within the file.
This data file is a desequenced file. The file is created with line numbers.
When it is entered into the computer it is desequenced. The file is read line by line
until the last line of data is read. On completion of reading the last line, the exe-
cution of the program is terminated.
Output Data
A sample of the output data is shown in Appendix C. The listing shows in-
put data for the particular run, the physical dimensions of the cooling tower and its
costs, the condenser cost, piping facilities cost, the cost of controls and the total
cooling system cost, including contingencies.
Operating Instruction for "TOWSIZ"
The program is ready to run once it and the data file SIZDAT are in the
computer. Check to make sure the data file has been desequenced.
"TQWSIZ" Flow Chart Description (1)
1 .01 Data statements are used in the program to define data tables in the pro-
gram which will be used to perform interpolation for various computations.
Included are tables of stack diameter, stack height, stack cost, tempera-
ture ranges, ambient air temperature height adjustment factors, altitudes,
altitude height adjustment factors, air flow, heat rejection, pressure drop
across the coil and water flow.
1 .02 Values of entered heat rejection and pressure drop tables are scaled to
proper dimensions for program use.
1-04 The input variables, initial temperature difference, range, heat rejection,
G0 , ambient air temperature and elevation, are read for one case from
one line of the data file named "SIZDAT".
(1) Numbers on these pages correspond to the
numbers on flow chart blocks.
20
-------
1 .05 The number of cooling columns is calculated by dividing condenser flow by
G0 . A temporary variable is then set to the remainder of the number of
coils divided by four. If the remainder is zero the calculated number of
coils is used. If it is not zero the number of coils is increased to the next
number evenly divisible by four. This action is taken to ensure the number
of deltas, which is one half the number of columns will be an even number.
1 .06 Seven curves of heat rejection versus air flow have been entered for seven
different values of water flow. The program compares the entered value of
water flow with the acceptable values corresponding to the curves. If no
match is found an error message is printed in block 1 .07. A match branches
to block 1.10.
1 .07 The following message is printed: G0 = XXXXX NO CORRESPONDING
CURVE. Case aborted.
1 .08 If there is still data in the input file, the program initiates the new case.
When no data is left, the program terminates.
1 .10 Subroutine TBLUQ is called. If the error flag for data outside tables is set
in the subroutine, the main program will branch to block 1.11. A normal
subroutine execution will lead to block 1 .12.
1.11 The following message is printed: L0 OUTSIDE TABLE . Processing con-
tinues .
1 .12 Pressure drop across the coil is found using air flow looking into air flow
and pressure drop tables.
1.13, Velocities and losses are calculated.
1.14,
1.15
1.16 Ambient air temperature is used to reference a table of temperatures and
adjustment factors; the reference temperature of 50.0° is used in the same
way; the site elevation references a table of elevations and adjustment
factors. The result of these calls to TBLUQ is a set of factors which are
used to adjust the tower draft height according to varying air densities at
temperatures other than 50.0 and elevations other than sea level.
1.17 Apply the factors from block 1 .16 to calculate adjusted draft height and
add 80.0 feet to obtain total tower height.
1.18 The final air temperature is used to find an adjustment factor for the tower
upper area.
21
-------
1.19 Compute the tower upper area and diameter and lower diameter.
1 .20 Test to find if lower diameter is less than upper diameter. If not, go to
1 .22, otherwise go to 1 .21 .
1 .21 Set lower diameter equal to upper diameter. Continue at 1 .22.
1 .22 The program finds two successive entries in a table of stack diameters which
bracket the calculated upper diameter. When these are found, go to step
1 .24. If no such entries are found, move to 1 .23.
1 .23 The following error message is printed: UPPER DIAMETER XXXX OUTSIDE
OF TABLES. Control is sent to block 1 .08.
1 .24 Routine TBLUQ is called twice, once for each of the two diameters found
in step 1 .22, using the total calculated tower height to find a tower cost
for each of the two diameters.
1 .25 The calculated upper diameter is used to perform a linear interpolation be-
tween the calculated costs from step 1 .24.
1 .26 The roof cost is calculated as the difference between the cross-sectional
areas of the top and bottom stack diameters times $7.60 per square foot.
1 .27 Basic condenser costs are provided to the program in a table of water cool-
ing range vs. condenser cost in dollars per 10" Btu per hour. The range is
used through TBLUQ to find a basic condenser cost which is multiplied by
the heat rejection in 10 Btu per hour to get total condenser cost.
1 .28 The number of deltas is one half the number of cooling columns.
1 .29 Call Subroutine TWRPIP to compute all piping and associated costs.
1 .30 Print out summary for sizing and cost evaluation.
Subroutine "TBLUQ" Flow Chart Description
2.01 Checks to see if value is within the given table of values.
2.02 If value is not in table it checks to see if the value is less than the smallest
value in the table.
2.03 If value is less than the smallest value then it uses the first three values to
extrapolate for the desired value.
22
-------
2 .04 If value is not in the table it checks to see if the value is greater than the
largest table value.
2.05 If value is greater than the largest table value then it uses the last three
table values to extrapolate.
2.06 Sets an error flag for any value outside the table.
2.07 Checks to see if value is equal to a table entry.
2.08 If value equals table value then it returns corresponding entry.
2.09, Checks for exact location in the table.
2.10
2.11, Uses three table entries to evaluate the value being looked up. If value is
2.12, between the first and second table entry then it uses the first three entries
2.13, for the evaluations. If value is between last two table entries it uses last
2.14 three table entries. Anywhere else in the table it uses the two table values
preceding and the one following the value.
2.15 Value is evaluated for from appropriate table values and returned to pro-
gram .
Subroutine "TWRPIP" Flow Chart Description
3.01 Calculate the water flow rate required for the transfer of the heat produced
by the turbine. This flow rate establishes the size and number of supply
lines required to pipe the water to the cooling tower.
3.02, Subroutine "PIPSIZ" is called to calculate the diameters of the pipe and
3.05, also pick the corresponding pipe and valve costs from a table of pipe and
3.09, valve costs.
3.12,
3.16
3.03, Calculate the quantity of pipe and valves required and the cost.
3.06,
3.10,
3.13,
3.17
3.04 Calculate the water flow rate for the header lines.
3.07 Calculate the number of cooling deltas per tower sector.
23
-------
3.08 Calculate the water flow rate for the sector drain pipe.
3.11 Calculate the water flow rate for the filler water pipes.
3.14 Calculate the number of bypass valves and their cost.
3.15 Calculate the water flow rate for the emergency drain pipe.
3.18 Calculate the number of circulating water pump setups and their cost.
3.19 Calculate the cost of filler pumps.
3.20 Calculate storage tank costs from a curve of gpm vs. tank cost.
3.21 Sum up the costs of piping, valves, pumps and storage tanks.
Subroutine "PIPSIZ" Flow Chart Description
4.01 Calculate the pipe size from the water flow rate.
4.02 Check to see if the diameter is less than or equal to the maximum size in
the table of pipe size vs. cost. If it is less than or equal to the maximum,
the program proceeds to 4.04.
4.03 If diameter is greater than maximum, add on another line and proceed to
4.01 .
4.04 Find corresponding pipe and valve cost from table of pipe size vs. pipe
and valve cost.
24
-------
APPENDIX A
PROGRAM LISTINGS
25
-------
ECONOMIC OPTIMIZATION PROGRAM "0PTDCT"
00100 PROGRAM 0PTDCT(0UTPUT,TAPE 1,TAPE2,TAPE3,TAPE4)
00120 REAL KW*KWH*NHR*KWLOS(3),KWHABV<51>,KWHBL0(51>
REAL T0TCAP<5]),CRDBTUC51),T0TBTUC51),SITEN0C27>
REAL KWL0SSC51>,DURC3)*KWHL0S,CCFGTC5)
REAL MILPKW*MWHC51>*L0ADC3>,HREJT(3>
REAL PCTL0D(3),TBLITDC6),AUXKWHC51)
REAL AUXKW(51),L0AD1(3)
REAL KWHGANC51)*TEMPC3)*KWAX(6)
DIMENSION P(56),T<56),AIRT(32>*TDURC32>* SI TEC 5>*FULCSTC3>
DIMENSION BPC35)*FTKWC35),GHRC35,3),HREJ(35,3)*CCF(5>*AUXTEMC2)
00140
00160
00180
00200
00220
00240
00260
00280
00300
00320
DIMENSION PCTTIM(3)*PCTLDC3)*PFCST(3>*TBLELVC3)*CAPCST(6*3)
DIMENSION TOWER(2),AUX1OOC6),AUXTMPC2*6, 3>,BC3)*AIR(3)
00340** AUXILIARY DATA CORRESPONDS T0 FOSSIL OR NUCLEAR FUELED UNIT
00360 DATA AUX100 /48000., 39000.*32000., 26000.*22000.* 18000./
00380 DATA AUXTMP /55.* 65.*45.,55.,30.* 45., \5.,35.* 5.,25.* - 20., \ 5 . ,
00400+ 60.,70.,50.,60., 45.,55.,35.,45.,25.,40.,15.,30»*
00420+ 70.*75.* 65.>75.* 55.,65.,50.*60.*40.*55.*30.*50./
00440 DATA KWAX / 1 9000.* 14000.* 11000.* 9000.* 8000.* 7000./
00460 DATA SITEN0 / 6HSITEO1 *6HSITE02*6HSITE03*6HSITE04*6HSITE05*
00430+ 6HSITE06*6HSITE07*6HSITE08* 6HSITE09*6HSITE10,6HSITE 1 1 *
00500+ 6HSITE12* 6HSITE13* 6HSITE14* 6HSITE15*6HSITE16*6HSITE17*
00520+ 6HSITE18*6HSITE19* 6HSITE20*6HSITE21,6HSITE22*6HSITE23*
00540+ 6HSITE24*6HSITE?5*6HSITE26*6HSITE27 /
00560** STEAM TABLE DATA (SATURATION TEMPERATURES AND PRESSURES)
00580 DATA P / . 2* . 3* . 4, . 5* . 6* . 7* . 8* . 9* 1 . 0* 1 . 1 * 1. 2* 1 . 3* 1 . 4* 1 . 5*
00600+ 1.6*1.7*1.8*1.9*2.0^2.1*2.2*2.3*2.4*2.5*2.6*2.7*2.8*2.9*
00620+ 3.0*4.*5.*6.*7.*8.*9«* 10-* 11.* 12.* 13.* 14.* 15., 16.* 17-* 18.,
00640+ 19«*20.*21.*22.*23.,24.*25.*26.*27.,28.*29.*29.921 /
00660 DATA T / 34.57*44.96*52-64*58.80*63.96*68.41*72.32*75.84*79.
03;
00680+
00700+
00720+
00740+
00760+
00780+
00800
00820
00840
00860
00880
00900
00920
81. 9 6* 84. 6 4* 87. 17* 8Q- 5 1*91. 72*93.81* 95. 78*97. 65*99. 43*
101. 14* 102.77* 104. 33* 105.85* 107.30* 108. 71, 1 10.06* 1 1 1 -37*
06*125.43*133
28*172.78*176
68* 196* 90*
1 12.63* 1 13-36* 1 15
161.49* 165.54, 169
189-96*192.37*194
203.67,210.43*212.00 /
DATA. PCTTIM / 50.*50.*0. /
DATA PCTLD / 100.*75.*50. /
DATA PCTL0D / 100.*75.*50. /
DATA 0PHRS / 7500. /
DATA HRSPYR / 8760. /
DATA ASTER / 2H** /
DATA TBLITD / 30.,40-* 50.* 60
76,140.78,146.
05, 179. 14, 182.
199.03,201.09,203.
86, 152
05,184
08*205
.24, 157.09*
.82, 187.45,
00,206.87,
70.,80.
00940** CAPITAL COST DATA F0R VARIOUS ITD'S AND ELEVATIONS
00960 DATA CAPCST/50000000.,36080000.* 26920000•*21530000•* 17720000.*
00980+15870000.* 54800000.* 37840000.* 28210000.* 22340000-* 22340000.*
01000+18160000.* 16220000.*60000000.,41250000.* 29890000.,23670000.*
01020+19100000.,16810000./
01040 DATA TBLELV /O.* 3000.,6000. /
01060 DATA TOWER/8HNATURAL *8HDRAFT /
01080 CALL RETR (1,6H0PTOUT)
26
-------
01100** SET UP AIR TEMPERATURES IN ARRAY *AIRT*
0 1 120 D0 10 I = 1,32
01140 AIRTCI) = 122-1*5
01 160 10 C3NTINUE
01180** SET UP BACK PRESSURE ARRAY
01200 00 15 I = 1,35
01220 BPCI) = (1+1.)/2.
01240 15 CONTINUE
01260** READ TURBINE DATA FR0M FILE *TURBIN*
01280 CALL RETR C2,6HTURBIN)
01300 READ C2, ) KW,I TYPE,HREJD,BPMIN,BPMAX
01320 KWH = KW*0PHRS*CPCTTIM(1)*PCTLD(1)+PCTTIMC2)*PCTLDC2)+
01340+ PCTTIM(3)*PCTLD(3))/l.OE4
01360 TMWH = KWH/1000.
01380** READ FULL THF
-------
GAS TURBINE
DURATIONS
JRE DURATIONS
02380
02400
02420
02440
02460
02480
02500
IF
32
33
DO
IF
35
02100** READ CAPITAL COST FACTORS FOR
02120 READ (3,) ( CCF GT( I ) , I = 1 ,NCCF )
02140** READ IN PEAKING FUEL C0ST
02160 READ (3,) ( PF CST ( I ) „ I = 1 , NFJL)
02180**READ IN CUTOFF TEMPERATURE FOR PEAK GENERATION
02200 READ (3*) CUTOFF
02220** READ IN COST/KW FOR AUXILIARIES
02240 READ (3,) CPKWAX
02260** READ IN OPERATION AND MAINTENANCE PERCENTAGE
02280 READ (3*) OAMCT
02300** READ IN TABLE OF TEMPERATURE
02320 READ (3,) ( TDURC I ) * I = 1 , 32)
02340** FIND FIRST AND LAST N ON- ZERO
02360 DO 32 NT = 1,32
(TDUR(NT).NE.O. 0) GO TO 33
CONTINUE
NTL = NT
35 NT = 1, 32
(TDUR(NT).NE. 0.0) NTU = NT
CONTINUE
CALL DATCHKC3)
02520** WRITE OUT SITE IDENTIFYING DATA
02540 WRITE (1,360) SITE
360 FORMAT C//*SITE - *5A6)
WRITE (1,361) ELEV
361 FORMAT (*ELEVATI0N -*F5.0* FEET ABOVE SEA LEVEL*)
WRITE (1,362) PEAK
362 FORMAT (*PEAK PERIOD - *A6)
WRITE (1,364) KW,ITYPE
364 FORMAT (*TURBINE -*,F7-0,* KW *A7* UNIT*)
WRITE (1,365) 0PHRS, PCTTIM ( 1 ) , ( PCTLD( I ) , PCTTIM ( I ) , 1 = 2, 3)
365 FORMAT (* IN OPERATION -*F5.0* HRS/YR*/
* AT FULL THROTTLE*, F4. 0* 0/0 OF THE TIME*
02560
02580
02600
02620
02640
02660
02680
02700
02720
02740+
ENEFiL-.Y -*F
MV.'H NOMINAL*)
02820
02840
02860
02880
*2A8)
02760+ 2(/* AT*F4.0* 0/0 LOAD*F5.0* 0/0 OF THE TIME*))
02780 WRITE (1,363) TMWH
02800 363 FORMAT (*TUTAL ANNUAL
WRITE (1,368) TOWER
368 FORMAT (*T0WER -
WRITE (1,931) OAMCT
931 FORMAT(*COOLINb
02900+F5.2* 0/0*)
02920 WRITE (1,366) HREJD,ASTER
02940 366 FORMAT (*DESIGN HEAT
WRITE (1,367) CCM
367 FORMAT (^CAPITAL
WRITE (1,468)CUTOFF
468 FORMAT(*CUT0FF TEMPERATURE
F*)
SET
02960
02980
03000
03020
03040+*
03060**
SYSTEM OPERATION AND MAINTENANCE COST -*
REJECTION -*,F5.0,* X 10*A2*6 BTU/HR*)
COST MULTIPLIER -*F5.2)
FOR PEAKING GENERATION -*F4.0
UP VARIABLES USED IN CALCULATING OPERATING COSTS
03080 LOAD( 1 ) =T-BLUQ(BPMAX, BP,FTKW, 35)
28
-------
03100
03120
03140
03160
03180
03200
03220
03R40
03280
03300
03320
03340
03360
L0ADC2) =
L0ADC3) =
L0AD1 (2)
L0AD1 (3)
HREJT(l)
HREJTC2)
HREJTC3)
DIJRC1) -
DURC3) =
AIRC1) =
A I R C 2 ) =
AIRC3) =
AIR10 = T
KW*.75
KW*.50
L0ADC2)
L0ADC3)
TBLUQCBPMAX,BP,HREJ(1,1), 35)
TBLUO ( BPM AX* 8P, HREJ ( 1 , 2) , 35 )
TBLUOC3PMAX>BP,HREJ(1> 3), 35)
TD>'J/+TDURCNTL+2)*HRfPYR/100.
AIRT(NTL)
AIRTCNTL+1>
AIRTCNTL+2)
TBLUQC10.,DUR,AIR,3)
03380** CALCULATE CAPITAL C0ST F0R EACH THETA CORRESPONDING
03400** T0 SITE ELEVATION
03420 D0 200 ITD = 30,80
03440 I = ITD-29
03460 THETA = ITD
03480 DG 1501 IK = 1,3
03500 1501 B(IK) = TBLUO ( THETA, T BL ITD, CAPC^'T C ] , ! K ) , 6)
03520 TOTCAP(I) = TBLUQCELEV*TBLELV>B*3)
03540** ADJUST CAPITAL C0ST F0R SITE LOCATION
03560 TOTCAPCI) = TOTCAPCI)*CCM
03580** INITIALIZE ACCUMULATIVE VARIABLES
03600 TCTBTUCI) = 0.0
03620 CRDBTUCI) =0.0
03640 KWHBL0CI) =0.0
03660 KWHABV(I) = 0.0
03680 KWHGANCI ) = 0. 0
03700 AUXKW(I) = 0.0
03720 AUXKWH(I) = 0.0
03740 KWH = 0.0
03760 IT = 0
03780 KWL0SC 1 ) = 0«
03800 KWL0SC2) = 0.
03820 KWL0SC3) = 0.
03840** DETERMINE COEFFICIENT 0F "A" 0F EQUATION:
03860** EXPONENT "Z" CORRESPONDS T0 T0WER TYPE,
03880** NATURAL DRAFT Z = .75 MECHANICAL DRAFT
03900 A = THETA/HREJD**.75
03920** DETERMINE OPERATING BACK PRESSURE AS INTERSECTION 0F THE
03940** CURVES 0F HEAT REJECTION VS. BACK PRESSURE F0R THE
03960** TOWER AT A GIVEN AIR TEMPERATURE AND F0R THE TURBIN
03980** AT A GIVEN LOAD
04000 D0 150 IL = 1,3
04020 DO 145 NT = NTL,NTU
04040 ITER = 0
04060 BPT =0.0
04080 0PBP = BPMIN
THETA = A*HREJ**Z
Z = .91
29
-------
GIVEN L0AD AND TEMPERATURE
04100 G0 T0 117
04120 110 SATT = AIRTCNT)+A*HRJ**.75
04140 IF (SATT-750 120*115,115
04160 115 SATP = TBLUQ(SATT*T*P*56)
04180 0PBP = SATP
04200 IF (ABSC0PBP-BPT).LT.0.01) 60 T0 120
04220 BPT =' 0PBP
04240 117 TEMP(l) = TBLUGK0PBP*BP*HREJC1*1)*35)
04260 TEMP(2) = TBLUQC0PBP*BP*HREJC1*2)*35)
04280 TEMPO) = TBLUQC0PBP* BP* HRE J( 1* 3)* 35)
04300 HRJ = TBLUQCPCTLDCIL)*PCTL0D*TEMP*3)
04320 ITER = ITER+1
04340 IF (ITER.GT-15) PRINT**M0RE THAN 15 ITERATIONS T0 FIND**
04360+ * 0PERATING BACK PRESSURE **/* I TD* AI RTCNT) * 0PBP* HRJ
04380 IF (ITER.GT.15) G0 T0 120
04400 G0 T0 110
04420** MAKE OPERATING BACK PRESSURE N0T LESS THAN THE MINIMUM
04440** SPECIFIED IN THE FILE *TURBIN*
04460 120 0PBP = AMAX1(0PBP*BPMIN)
04480** CALCULATE THE NUMBER 0F H0URS AT THE
04500 HRS = 0PHRS*TDUR(NT)*PCTTIMCIL)/1.0E4
04520** DETERMINE KW LEVEL 0F 0PERATI0N P0SSIBLE
04540 BPT = AMIN1(0PBP*BPMAX)
04560 FTKWT = TBLUQCBPT*BP*FTKW,35)
04580 IF C0PBP.GT.BPMAX) G0 T0 122
04600 IF (PCTLD(IL).EQ.100.) TKW = FTKWT
04620 IF CPCTLDdD.LT. 100. ) TKW = PCTLDC IL)*KW/ 100.
04640 IF (TKW.GT.FTKWT) TKW = FTKWT
04660 G0 T0 126
04680** ADJUST KW LEVEL IF BACK PRESSURE GREATER THAN ALLOWABLE
04700 122 SATT = TBLUQCBPMAX*P*T*56)
04720 PLITD = SATT-AIRTCNT)
04740 HREJMX = (PLITD/A)**(1./-75)
04760 TKW = TBLUQ(HREJMX*HREJT,L0AD*3)
04780 0PBP = BPMAX
04800** CALCULATE THE T0TAL YEARLY BTU'S C0NSUMED F0R EACH 0F THE
04820** THREE L0ADINGS 0F 100* 75* AND 50 0/0
04840** DETERMINE GR0SS HEAT RATE F0R CALCULATED BACK PRESSURE
04860 126 TEMP(l) = TBLUGK0PBP, BP* GHR( 1* 1 ) , 35)
04880 TEMP(2) = TBLUGK0PBP, BP* GHRC 1 * 2) * 35)
04900 TEMPO) = TBLUGK0PBP* BP* GHRC 1 * 3) * 35 )
04920 L0AD1C1) = TBLUQC0PBP* BP* FTK W, 35)
04940 0PGHR = TBLUQCTKW* L0AD 1 * TEMP, 3)
04960 T0TBTU(I) = T0TBTUCI)+0PGHR*HRS*TKW
04980** DETERMINE AUXILIARY P0WER REQUIREMENTS F0R THE VARYING THETAS
05000 IF (T0WER(1).EQ.8HNATURAL ) G0 T0 1320
05020** AUXILIARY CALCULATI0NS F0R MECHANICAL DRAFT T0WER
05040 FULAUX = TBLUGKTHETA, TBLI TO* AUX 1 00* 6)
05060 D0 1210 IJ = 1*5
05080 IF ((THETA.GE.TBLITD(IJ)).AND.(THETA.LE.TBLITD(IJ+1))) G0 T0 1220
30
-------
05100 1210 CONTINUE
05120 PAUSE 4910
05140 1220 IJ1 = IJ+i
05160 IJ2 = IJ
05180 AUXTEMC1) = AUXTMPC 1 * IJ2, ID - ( ( ( THETA-TBL I TDC IJ) )/ 1 0- ) *
05200+ CAUXTMPC1,IJ2, ID-AUXTMPC 1, IJ 1 , ID >)
05220 AUXTEMC2) = AUXTMPC2, IJ2, ID -(( (THETA-TBL I TDCI J) )/ 1 0. )*
05240+ CAUXTMPC2,IJ2, ID - AUXTMPC 2* IJ \, ID ) )
05260 IF CAIRTCNT5.LE.AUXTEMC1)> G0 T0 1240
05280 IF CAIRTCNT).GE-AUXTEMC2)) G0 T0 1250
05300 PCT = 0.53+0.47*C(AIRT(NT)-AUXTEM(1))/(AUXTEMC2>-AUXTEMC1»>
05320 G0 T0 1260
05340 1240 PCT = 0.53
05360 G0 T0 1260
05380 1250 PCT = 1.00
05400 G0 T0 1260
05420 1260 AUXKWH(I) = AUXKWH(I) +FULAUX*PCT*HRS
05440 AUXKWCI) = FULAUX
05460 G0 T0 1290
05480** AUXILIARY CALCULATIONS F0R NATURAL DRAFT T0WER
05500 1320 AUXKWCI) = TBLUQCTHETA,TBLITD,KWAX, 6)
05520 AUXKWHCI) = AUXKWHCI)+AUXKWCI)*HRS
05540 1290 KWH = KWH + TKW*HRS
05560** CALCULATE L0SS 0R GAIN 0F CAPACITY
05580 CAPL0S = KW*PCTLDC ID/1 00.-TKW
05600 IF CCAPL0S.GT.0.0) G0 T0 130
05620** GAIN IN CAPACITY CALCULATIONS
05640 KWHGANCI) = KWHGANCI)-CAPL0S*HRS
05660 CRDBTUCI) = CRDBTUCI)-CAPL0S*HRS*0PGHR
05680 132 IT = IT+1
05700 IF CIT.GT.3) GO TO 145
05720 KWL0SCIT) = CAPL0S
05740 G0 TO 145
05760** LOSS 0F CAPACITY CALCULATIONS
05780 130 IFCAIRTCNT).GE.CUT0FF) GO T0 131
05800 KWHBL0CI) = KWHBLOCI) + CAPLOS*HRS
05820 G0 T0 132
05840 131 KWHABVCI) = KWHABVCI)+CAPL0S*HRS
05860 GO T0 132
05880 145 CONTINUE
05900 150 CONTINUE
05920 MWHCI) = KWH/1000.
05940** CALCULATE L0SS OF CAPACITY AT THE 10 HOUR AIR TEMPERATURE
05960 KWL0SSCI) = TBLUGCAIR10>AIR*KWL0S>3)
05980 IF CCKWL0SC2)+KWL0SC3)).EQ.O.O.AND.DURC1).LT.10.) KWL0SSCI) =
06000 200 CONTINUE
06020 D0 450 K = 1,NFUL
06040 D0 450 J = 1/NCCF
06060** WRITE OUT INFORMATION IDENTIFYING EACH RUN
06080 WRITE Cl*370) CCFCJ),CCFGTCJ)>CCFCJ),FULCSTCK),ASTER,PFCSTCK),
06100+ ASTER*PCCST,CPKWAX
31
0-0
-------
06120 370 F0RMAT(*CAPITAL C0ST FACTORS:*9X*PLANT -*F3.0* 0/0*3X*PE AK ING*
06140+* CAPACITY -*F3.0* 0/0*3X*AUXILIARIES -*F3.0* 0/0*/*PLANT FUEL *
06160+*C0ST -*F3.0* CENTS/10*A2*6 BTU*1OX*PEAKING FUEL C0ST -*F3«0
06180+* CENTS/10*A2*6 BTU*/*PEAKING CAPITAL C0ST -*F4.0* S/KW*16X
06200+*AUXILIARY CAPITAL COST ~*F4.0* $/KW*/>
06220** WRITE 0UT C0LUMN HEADINGS F0R THE DATA OUTPUT
06240 WRITE (1,380)
06260 380 F0RMATC68X*ANNUAL*/68X*CAPITAL*3X*ANNUAL*/17X*EXCESS*45X*AND *
06280+*0+M*4X*FUEL*
06300+33X*T0TAL ANNUAL*/*INIT.*3X*GR0SS*4X*ENERGY*47X*C0ST*5X*C0ST*3X
06320+*CREDIT*23X*CGST 0F C00LING*/*TEMP.*3X*ENERGY*3X*DUE T0*2X*CAPACI*
06340+*TY*21X*MAXIMUM*8X*0F DRY*5X*0F*5X*F0R*3X*CAPACITY*13X*SYSTEM AND*
06360+* T0TAL*/*DIFF.*3X*800 MW*4X*EXTRA*3X*PENALTY*IX*AUXILIARY*1X*L0S*
06380+*S 0F*2X*AUXILIARY*6X*C00LING*3X*800 MW*2X*EXCESS*2X*PENALTY*IX
06400+*AUXILIARY*6X*PLANT FUEL*/*(DEG)*4X*UNIT*3X*CAPACITY*2X*ENERGY*3X
06420+*ENERGY*3X*CAPACITY*3X*FOWER*9X*SYSTEM*4X*UNIT*3X*EN£RGY*
06440+3X*C0ST*5X
06460+*C0ST*5X*==================*/*( F )*,4C4X*CMWH)*), 5X*(KW)*6X*CKW)*
06480+ 1 1X*CS)*6X*(S)** 2C5X*( 3)=;-) , 6X* C 3) >=RX-• C S) ^^X-d ; ILL/KVJH) * )
06500 ITD = 29
06520 TANT = 1.0E50
06540 REWIND 4
06560 D0 425 I = 1,51
06580 ITD = ITD+1
06600** DETERMINE ANNUAL CAPITAL C0STS F0R EACH CAPITAL C0ST FACT0R
06620 ANNCAP = T0TCAPCI)*CCF(J)/100.
06640 0PMAT = T0TCAPCI)*0AMCT/100.
06660 ACCA0M = ANNCAP+0PMAT
06680** CALCULATE T0TAL FUEL C0ST FOR THE AM0UNT 0F BTU'S USED
06700** F0R THAT YEAR C0MBINED WITH THE DIFFERENT FUEL C0STS
06720 ANNFUL = T0TBTUCI)*FULCSTCK)/1.OE8
06740** CALCULATE CREDIT FOR EXCESS ENERGY PRODUCED
06760 CREDIT = CRDBTUCI)*FULCSTCK)/1.OE8
06780**CALCULATE CAPACITY L0SS CHARGES
06800 IF CI'VJSP.EQ. 1 > G0 T0 137
06820** SUMMER PEAKING CAPACITY L0SS CHARGES
06840 PENCAP = KWL0SSCI)*PCCST*CCFGT(J)/100.+1»2*KWL0SSCI)
06860 PEN0PR= KWHABVCI)*15.*PFCSTCK)/1.OE5+CKWHBL0(I)/1.OE4)*
06880+ (FULCSTCK)+1.)
06900 G0 T0 138
06920** WINTER PEAKING CAPACITY LOSS CHARGES
06940 137 PENCAP = 0.0
06960 PEN0PR = ((KWHABV*CFULCSTCK)+1.)
06980** SUMMATION 0F CAPACITY L0SS CHARGES
07000 138 PENLTY = PENCAP+PEN0PR
07020 KWHLOS = KWHABVCI>+KWHBL0(I)
07040**CALCULATE AUXILIARY POWER REQUIREMENT CHARGES
07060** ASSUME 0 AND M CHARGES ARE .10 MILLS/KWH
07080 390 0NM '=.10
32
-------
CORRESPONDING
07260 MILPKW =
07280** CHANGE
07300 KWHL0S =
07320 TAXKWH =
07340 TKWHGN =
07100 AUXCAP = AUXKWd)*CPKWAX
07120** CALCULATE MILLS/K.WH FIGURE F0R AUXILIARY KWH
07140 FCSTAX = ANNFUL/MWHCI)+0NM
07160 AUX0PR = AUXKWH(I)*FCSTAX/1000.
07180 AUXCST = AUXCAP*CCF(J)/100.+AUX0PR
07200** COMBINE OPERATING C0STS AND CAPITAL C0STS F0R
07220** THETAS
07240 TANCST = ACCAOM+ANNFUL+PENLTY+AUXCST-CREDIT
TANCST/TMWH
ENERGY VALUES T0 MWH
KWHL0S/1000.
AUXKWHCI)/1000«
KWHGANCI )/1000.
07360** DETERMINE THE OPTIMUM ITD
07380 IFCTANT.LT.TANCST) G0 T0 3
07400 TANT = TANCST
07420 ISAVE = ITD
07440** WRITE DATA IN 0UTPUT FILE
07460 3 WRITE(4,420)ITD,MWH(I),TKWHGN*KWHL0S*TAXKWH*KWL0SS(I)*
07480+ AUXKWd)* ACCA0M, ANN FUL* CREDIT* PENLTY* AUXCST* TANCST* MILPKW
07500 420 FORMAT ( I 4, F11.0,3F9.0, 2F10«0,F14.0,F9.0,F7.0*F9.0,
07520+ F10.0,F12.0,F8. 4)
07540 425 CONTINUE
07560 REWIND 4
07580 D0 1020 1=1*51
07600 READ ( 4, 420)ITD*MWHCI)*TKWHGN*KWHL0S*TAXKWH*KWL0SSCI ),
07620+ AUXKWd ), ACCA0M, ANNFUL*CREDIT*PENLTY*AUXCST,TANCST*MILPKW
07640 IF CITD.NE.ISAVE) G0 T0 1020
07660 WRITE (1,611)
07680 611 F0RMAT (*0PTIMUM:*)
07700 1020 WRITE (1, 420)ITD,MWHCI),TKWHGN,KWHL0S,TAXKWH,KWL0SSCI),
07720+ AUXKWd )*ACCA0M*ANNFUL,CREDIT*PENLTY* AUXCST* TANCST, MILPKW
07740 450 C0NTINUE
500 C0NTINUE
ENDFILE 1
REWIND 1
CALL REPL C1,6H0PT0UT)
ST0P
END
07880** DATA FILE CHECK SUBROUTINE
07900 SUBROUTINE DATCHK(NFIL)
DIMENSION LINEC24)
INTEGER FILEC3),ICNTC3)
DATA FILE /6HSPVSST,6HTURBIN,6HSITE /
READ (NFIL,) DATA
ICNTCNFIL) = ICNTCNFID+1
IF CDATA-1.OE50) 20,10,20
10 RETURN
20 BACKSPACE NFIL
READ (NFIL,30) LINE
07760
07780
07800
07820
07840
07860
07920
07940
07960
07980
08000
08020
08040
08060
08080
33
-------
- AB0VE LINE*
08100 30 FORMAT (S4A3)
08120 PRINT 40*LINE*FILE(NFIL)
08140 40 FORMAT (/24A3/*ERR0R IN READING FILE *A6*
08160+ * SH0ULD BE DATA CHECK NUMBER (1.0E50)*/)
08180 REWIND NFIL
08200 NCNT = ICNT(NFIL)
08220 D0 50 I = 1*NCNT
08240 45 READ (NFIL*) DATA
08260 IF (DATA-LOE50) 45*50*45
08280 50 CONTINUE
08300 RETURN
08320 END
08340** TABLE L00KUP SUBR0UTINE USING A CURVILINEAR INTERPOLATION
08360**TABLE L00KUP ROUTINE USING A DETERMINANT S0LUTI0N 0F TW0
08361** SEC0ND ORDER EQUATIONS
08380 DIMENSION XTC1 00)* YT ( 1 00)
((XT( 1)-X)*(XT(NVAL)-X)) 3*3*2
(ABSCX-XTC 1) >.LT- ABS CX-XKN VAL ) ) ) N = 1
(ABS(X-XT( 1)).GT. ABS(X-XT(NVAL))) N = NVAL-2
T0 25
NVAL-1
10 I = 1*N
IF(X.NE.XTCI)) G0 T0 5
) = YT(I)
IN
(CX.GT-XT(I)). AND. (X.LT.XT(H-l))) G0 T0 15
((X.LT-XTd )). AND. (X- GT.XT(I+1) )) GO T0 15
CONTINUE
08400
08420
08440
08460
08480
08500
08520
08540
08560
08580
08600
08620
08640
08660
08680
08700
08720
08740
08760
08780
08800
08820
08840
08860
08880
08900
08920
08940
08960
08980
09000
IF
2 I
IF
G0
3 N
D0
IFC
TBL
RET
5 I
IF
10
N =
G0
15
16
GO
20 I
25
X2
X3
Yl
Y2
Y3
X1S
X2S
X3S
DET
A =
B =
C =
NVAL-2
T0 25
IF (1-3)
N = 1
16* 20*20
TO 25
N = 1-1
XI = XT(N)
= XT(N+1)
= XKN + 2)
= YT(N)
: YT(N-H)
: YT(N+2)
= X1*X1
= X2*X2
= X3*X3
= X2S*X3-X3S*X2-X1S*X3+X3S*X1+X1S*X2-X2S*X1
(Y2*X3-X2*Y3-Y1*X3+Y3*X1+Y1*X2-Y2*X1)/DET
(X2S*Y3-X3S*Y2-X18*Y3+X3S*Y1+X1S*Y2-X2S*Y1)/DET
(Y1*(X2S*X3-X3S*X2)-Y2*(X1S*X3-X3S*X1)+Y3*(X1S*X2-X2S*X1)
09020+ )/DET
09040 TBLUQ = A*X*X+B*X+C
09060 RETURN
09080 END
34
-------
TOWER OPTIMIZATION PROGRAM "TOWSIZ"
00100 PROGRAM T0WSIZ(TAPE 1,OUTPUT)
00120C
00140C
00160C
00180C
00200C
00220C
00240C
00260C
00280C
00300C
00320C
00340C
00360C
OC380C
00400C
00420C
00440C
00460C
00480
00500C
00520C
00540C
00560C
C0580C
00600C
00620C
00640C
00660C
00680C
00700C
00720C
00740C
00760C
00780C
00800C
00820C
00840C
00860C
00880C
00900C
00920C
00940C
00960C
00980
0 1000
0 1020
0 1040
0 1060
0 1080
DRY C00LING T0WER SIZING AND C0ST EVALUATION PROGRAM
6/30/70 PJB
PROGRAM
NAME
ITD
RANGE
HTREJ
GO
AMAIRT
ELEV
CONDF
N
QO
CURVGO
LO
Q1TBL
L0TBL
DELTAP
DPTBL
AEXIT
EXITV
EXITLS
DRFTLS
TOTL0S
AIRHT
AIRDLT
FAIRTM
DRAFTH
ADRFTH
TOTALH
ADRT
ADRA
TEMPS
TFACTR
HEIGHT
AFACTR
UPAREA
UPDIAM
L0WDIA
PARAMETERS
USE
INITIAL TEMPERATURE DIFFERENCE
TEMP CHANGE OF C0NDENSATE
HEAT REJECTION
AMBIENT AIR TEMPERATURE
ELEVATION
C0NDENSATE FLOW
NUMBER OF COLUMNS
C7> TABLE 0F CONSTANT GO TO IDENTIFY CURVE 1
(31,7) TABULAR COORDINATES FROM CURVE 1 PLOT
31) LO FOR CURVE 1 AND DELTAP
COIL LOSS
(31) TABLE OF DELTAP VALUES
AIR EXIT VELOCITY
CONSTANT (20 FT/SEC) EXIT VELOCITY
EXIT LOSSES
DRAFT LOSSES
TOTAL LOSSES
AIR HEAT GAIN
AIR TEMPERATURE INCREASE
FINAL AIR TEMPERATURE
DRAFT HEIGHT-
ADJUSTED DRAFT HEIGHT
TOTAL DRAFT HEIGHT
AIR DENSITY RATIO - TEMPERATURE
AIR DENSITY RATIO - ELEVATION
( ) TEMPERATURE TABLE
( ) TEMP DENSITY RATIO TABLE
( ) ALTITUDE TABLE
( ) ALTITUDE DENSITY RATIO TABLE
AREA OF TOWER TOP
DIAM OF TOWER TOP
DIAM 0F TOWER BOTTOM
INTEGER STPTSC4)
REAL (M,LO
REAL LOTBL,L0WDIA
REAL ITD
DIMENSION CURVGO(7),01TBL(31* 7),LOTBL(31),DPTBL(31)
DIMENSION IENTR(7)
35
-------
01100 DIMENSION TEMPSC 1 8 ) , TFACTRC 1 8 ) * HEI GHT C 24) , AFACTRC 24)
0 1 120 DIMENSION STDIAM(4),STKHT<5, 4),STKCSTC5, 4) ,RTEMP(6), CNDCSTC6)
01140 DIMENSION DUMMYC31)
0 1 160 COMMON JJ
0 1 180 C0MM0N HIRE J, RANGE* DELTAS* LO WDI A, TPCST, TVCST, PC0ST,FPCS1
01200 COMMON STCST, T0TPCS, GPM
0 1220 EQUIVALENCECDUMMYC 1 ) * Q 1 TBL ( 1 , 7 ) )
01240 DATA STDIAM/21 7., 312. , 460. , 550. /
0 1260 DATA STKHT/ 1 20- » 240. , 3«7 . , 527- , 79 3. , 126.* 28 3. ,463., 6 03.*
01 28 0+800. ,126. ,28 3«, 463-, 603«, 813«, 200«, 400., 600., 800., 0«0/
0 1300 DATA STKCST/272466., 574247. ,897032.* 1470194. ,3111314. ,361950.,
0 1320+894091 ., 1565409. ,231 4192., 3971000. , 516062. , 1264230., 2067594.,
0 1 340+3326040. ,5360366. , 980000. , 2050000. , 3730000 ., 5980000., O./
01360 DAIA STPTS/5,5,5, 4/
013RO DATA RTEMP/5. 0, 1 0. 0, 20. 0, 30. 0, 40. 0, 50. O/
0 1 400 DATA CNDCST/273000. , 256000. , 22*000. , 208000. , 1 90000. , 1 67000. /
01420 DATA TEMPS /-20« , - 1 0. , 0. , 1 0. , 20. , 30. , 40. , 50. , 60. , 70. , 80. ,
0 1440+90., 100. , 1 10-, 120. , 1 30. , 140., 1 50 •/
01460 DATA TFACTR / 1 . 1 97, 1 . 1 75, 1 . 1 7 1 , 1 . 1 25, 1 . 1 02, 1 . 078, 1 . 057, 1 . 037,
0 1480+ 1 • 018, 1 .00, -986, .962,. 944, .933, .914, .899, .884, . 8 70/
01500 DATA HEIGHT /O. , 400. , 800. , 1 200. , 1 600. , 2000 ., 2400 ., 2800 •, 3200 •,
01520+ 3600. ,4000., 4400. ,48 00. ,5 200. ,56 00. ,6 000. ,6400. ,68 00. ,7200.,
0 1 5 40+ 7 6 00., 8 000., 8 400., 88 00., 9200. /
01560 DATA AFACTH /l.,.«85 , - «7 2, . ^58, . 944, - 93, . 9 1 7, . 904, . 89 1 , • 878,
0 1580+.864, .852, . 832, .826, -813, .802, . 788, • 776, . 755, . 754, . 743, . 73, .719,
0 1600+.70R/
01620 DATA IENTH / 3 1 , 3 1 , 3 1 , 3 1 , 27, 2 7, 3 1 /
01640 DATA. LOTBL /200. , 2 1 0 • , 220- , 230. , 240. , 250. , 260. , 270. , 280. , 290. ,
01660+ 300. ,310. ,320. ,330. ,340. ,350. ,360. ,370. ,380. ,390.,
0168 0+ 400. , 410. , 42 0., 4 30. ,440. ,450. ,460. ,4 70. ,480. ,490., 500. /
01700 DATA Q1TBL /35 . 6, 36. 6, 37 . 8, 38- 8, 3=> • ", A \ . , ;t? . l , /i3 • 2, 44. 3, 45 . 3,
0 1720+46. 4, 47. 4, 48. 4, 49. 2, 50. 1,51. ,51. 8, 52. 6, 53. 4, 54.2, 54. 9, 55. 6,
01740+56.
0 1760+34.
0 1780+45.
0 1800+52.
0 1820+33.
01840+43.
0 1860+51 .
0 1880+33.
0 1900+42.
0 1920+48.
0 1940+32.
0 1960+40.
0 1980+46.
02000+31 .
02020+38-
3,57.
6, 35.
3,46.
8,53.
7,34.
9, 44.
0,51.
1,34.
3, 43.
8, 49.
2, 33.
8, 41 .
5,46.
0, 31 .
2, 38.
, 57. 6, 58.2,58.8
5, 36.
1, 46.
4, 54.
7,35.
7,45.
6, 52.
0, 34.
0, 43.
3, 49.
0, 33.
4, 42.
9, 47-
7,32.
8, 39.
6,37.
8,47.
0,54.
7,36.
4, 46.
2,52.
9,35.
7,44.
8,50.
9,34.
0,42.
3,47.
4, 33.
3,39.
7, 38.
5, 48.
6, 55.
7,37.
1, 46.
7,53.
8, 36.
4, 45.
3,50.
8, 35.
6, 43.
6,48.
1,33.
8, 40.
,59. 4
7,39.
3,49.
2,55.
7,38.
8,47.
2,53.
7,37.
0,45.
8,51.
6,36 .
2,43.
0,48.
7, 34.
3, 40.
, 60. ,
8, 40.
0,49.
8, 56.
6, 39.
5,48.
8,54.
6,38.
6, 46.
3,51.
4,37.
8, 44.
4, 48-
4, 35.
7, 41 .
60. 6,
8,41 .
6,50.
3, 56.
6,40.
1, 48.
3,54.
4,39.
2,46.
7,52.
2, 38-
3,44.
8,49.
1,35.
1,41.
61 .2,
7, 42.
3, 51.
8, 57.
5,41.
7, 49.
8, 55.
3, 40.
7, 47.
1,52.
0, 38.
8,45.
1, 49.
8, 36.
5,41.
7,43.
0,51.
4,
4, 42..
3, 49.
3,
1, 40.
2,47.
5,
7, 39.
2, 45.
4,
4, 37.
9,42.
6, 44.
6, 52.
3, 43.
9, 50.
8, 41 .
7,48.
4,40.
6, 46.
0, 37.
3, 42.
5,
2,
1,
4,
6,
3,
2,
0,
6,
7,
02040+43. 1,43.4, 43.8,44.1,44.4,0.,0.,0.,O.X
02060 DATA DUMMY/ 25. 9, 26. 4, 25. Q, 2 7. 3, 27- 8, 28. 2, 28. 6, 2^. 0, 2^- 5, 29. 9, 30
02080+30. 6, 30. 9, 31. 2, 31. 5, 31. 7, 31. 9, 32. 2, 32. 4, 32. 5, 32. 7, 32.9,
02100+33. 1, 33. 2, 33. 4, 33. 6, 33. 7, 33. 8, 34. 0,35. 1,35.2 /
36
-------
176320-*154280-*132240.*
D0
10
20
02120 DATA DPI 8L / 1 8. 5* 1 8. 7* 1 8. 9* 1 9. 1 * 1 9 • 3* 1 9 • 5* 1 9 • 7* 1 9. 9* 20. 1
021 40+20. 3,20. 6, 20. K, 21 . 0, 21 . 2, 21 . 4, 21 . 6*21 .8*22.0*22- 3*22. 5*
02160+22.8* 23. 0*23.2, 23. 4, 23. 7,24.0, 24. 3,24. 6,24.8*25.2* 25. 5 /
02180 DATA CURVGO/264480.* 220400.* 1^8360.,
02200+88160./
02220 DATA EXITV/20.0/
DATA CONTRL/500000*/
CALL HE1RC 1, 6HSI ZDA'l )
DO 20 K=l,31
10 L=l,7
Q 1 T BL ( K , L ) = 0 H SL ( K , L ) * 1 0 0 0 .
DPTBL(K)=DPTBL(K)-16.0
INPUT BASIC SIZING AND COST EVALUATION DATA
READ (1,) I1D, RANGE, HTfi;:j, GO, AKATP.7*ELEV
C'JMPUTE INTERMEDIATE SIZING DATA
CONDF=HTREJ/RANGE
N=CONDF/GO
TEMP=AMUD(N,4.0)
IF (TEMP.ME.0.0) N=N-TEMP+4. 0
QO=CHTREJ/N)*0.252
Q1=00/(ITD*0.55555)
90 DO 100 1=1,7
IF (GO.EO. CURVGOd ))
100 CONTINUE
PRINT, *GO =*,GO,* NO
G0 TO 9000
110 ICURV=I
LO=TBLUQ(Q1,Q1TBL(1,ICURV)*LOTBL,TEN1R(ICURV),NFLAG)
IF (NFLAG.EO.2) PRINT**LO OUTSIDE TABLE*
DELTAP=TBLUQ(LO,LOTBL,DPTBL,31,NFLAG)/25,
AEXIT=EXI7V
EXITLS=((AEXIT/4005.)*60.0)**2
DRFTLS=0.15*EXITLS
TOTLOS=DELTAP+EXITLS+DRFTLS
AIRHT=GO/(LO*555.408)
AIRDLT=AIRHT/0.24
FAIRTM=50.0+AIRDLT
DRAFTH=TOTL3S/(7.659*C1.0/5!0.-!•0/CA!RDl
ADRT=TBLUQ(AMAIRT*TEMPS*TFACTR*18*NFLAG)
ADRT1=TBLUQ(50.0*TEMPS,TFACTH,18,NFLAG)
ADRA=TBLUO(ELEV* HE] i-iHl , AF ACT-<* 24,NFLAG)
ADRFTH = DRAFTH* ( (ADR'i 1 / ADRT ) **3) * ( (1-0/ADRA)**2)
TCTALH=ADRFTH+80.0
ADRT = TBL!JQ(FAIRTM,TEMPS, TFAC7R* 18)
UPAREA=N*LO/(ADRT*AEXIT)*8«174
UPDIAM=(JPAREA*!.2732)**0« 5
L0WDIA=3.28*(0.717+0«899*(N/4.0))
IF (LOWDIA.LT.UPDIAM) LOUDIA=UPDIAM
02240
02260
02280
02300
02320
02340
02360C
02380 50
02400C
02420
02440
02460
02480
02500
02520
02540
02560
02580
02600
02620
02640
02660
02680
02700
02720
02740
02760
02780
02800
02820
02840
02860
02880
02900
02920
02940
02960
02980
03000
03020
03040
03060
G0 TO 110
CORRESPONDING CURVE. CASE ABORTED.*
4
1 0. ) ) )
37
-------
T0WER SIZES
GO T0 2010
03080C COST EVALUATION SEGMENT
03100C VARIABLE USE
03120C STDIAM (3) TABLE OF C0ST EVALUATION
03140C STKHT (3*5) TABLE 0F STACK HEIGHTS
03160C STKCST (3*5) TABLE JF STACK COSTS
03180C STACK EVALUATEED COST 0F STACK
03200C DROOF EVALUATED COST 3F DELTA ROOF
03220C COILS COST OF COOLING COILS
03240C RTEMP (6) TEMP RANGE TABLE FOR CONDENSER C0STS
03260C CNDCST (6) COST TABLE FOR CONDENSER EVALUATION
03280C CPBTU CONDENSER COST PER BILLION BTU
03300C DNDNSR T0TAL CONDENSER COST
03320C TOWER STRUCTURE
03340 DO 2000 1=1*3
03360 IF ((UPDIAM.GE.STDIAM(I)).AND.(UPDIAM.LE.STDIAM(I+1)))
03380 2000 CONTINUE
03400 PRINT**UPPER DIAMETER = **UPDIAM** OUTSIDE 0F TABLES*
03420 GO TO 9000
03440 2010 A1=TBLUQ(T0TALH*STKHT(1,I)*STKCST(1*I)*STPTSCI>*NFLAG)
03460 A2 = TBLUQ(T:JTALH* STKHT ( 1, I + 1 ) , STKCST ( 1 > I + 1 ) , SI PTSCI + 1 )*NFLAG)
03480 STACK=A1+(A2-A1 )*< C UPDI AM-STDI AM ( I ) ) / C STDI AM (I + 1 )-STDIAM (I ) ) )
03500 DR00F=0.78539*(L3WDIA**2-UPDIAM**2)* 7.60
03520 C0ILS=14500.*N/2.
03540C CONDENSER COST
03560 CPBT!J=TBLUQ(RANGE*RTEMP*CNDCST*6*NFLAG)
03580 CNDNSR=CP3TU*(HTREJ/1.OE09)
03600 DELTAS=N/2.0
03620 CALL TWRPIP
03640 PRINT 1001* I TD, RANGE*HTREJ*GO*AMAIRT*ELEV,T0TALH* UPDIAM*L0WDIA*GPM
03660 1001 FCRMATC/////70C*-*>/*DRY CJQLING TOWER SIZING AND COST *
03680+*EVALUATION*///*DESIGN PARAMETERS*//9X*ITD = *F5«0*5X*RANGE = *
03700+F5.0/9X*HEAT REJECTION = *E9.1/9X*WATER FLOW PER HOUR = *E10«1/
03720+9X*AMBIENT AIR TEMP = *F5.0*4X*ELEVATI0N = *F7.0////*T0WER SIZING*
03740+//9X*TQWER HEIGHT = *F6.1/9X*UPPER DIAMETER = *F6.1/9X*BOTT0M *
03760+*DIAMETER = *F6.1/9X*GALL0NS PER MINUTE = *F9.0)
03780 TTOWER = STACK + DROOF + C'JILS
03800 PRINT 1 002* STACK, DH JJr , C JIL-S* ~ TOWER* CNDNSR
03820 1002 FORMAT C///*C'JST EVALUATI0N*/4X*TOWER STRUCTURE*//9X
03840+*STACK COST*8X*F11.0/9X*SHED COST*9X*F11.0/9X*C0IL C0ST*9X*
03860+F1 1.0//9X*TOTAL STRUCTURE*12X*F1 1.0///4X*CON DENSER*//
03880+9X*C0NDENSER COST*13X*F11.0)
03900 T0T=TOTPCS+TT0WER+CNDNSR+CONTHL
03920 PRINT 1 003* TPCST, TVCST* PC3ST* FPCST*'STCST, TOTPCS* C0NTRL* T0T
03940 1003 FORMAT C///4X*PIPING* VALVES* ETC.*//9X*PIPE COST*9X*F11.O/
03960+9X**VALVE C0ST*8X*F11.0/9X*PUMP C0ST*9X*F11.0/9X*FILLER PUMP C0ST*
03980+2X*F11.0/9X*ST0RAGE TANK COST *F11.0//9X*T0TAL PIPING FACIL*
04000 + *ITIES*4X*F11.0//4X*CONTROLS*//9X*CJNTR0L COST*15X*F 1 1.O///
04020+*COMPLETE TOWER FACILITIES*///9X*T0TAL T0WER C0ST*9X,F13.O//)
04040 T0T=TOT*1.25
04060 PRINT 1004*TOT
04080 1004 FORMAT (*TOTAL TOWER COST AND C0NTINGENCIES*F13.O///70(*-*))
38
-------
04100
04120
04140
04160
04180
04200
04220
04240
04260
04280
04300
04320
04340
04360
04380
04400
04420
04440
04460
04480
04500
04520
04540
04560
04580
04600
04620
04640
04660
04680
04700
04720
04740
04760
04780
04800
04820
04840
IF CENDrlLE
STOP
1) 9100,50
9000
9100
END
FUNCTION TBLUQ (X,XT,YT,NVAL,NFLAG)
DIMENSION XTC100),YTC100)
IF CCXTC1)-X)*CXTCNVAL)-X)) 3,3,2
2 IF CABSCX-XTC1)).LT-ABSCX-XT(NYAL))) N = 1
IF CABSCX-XTC1))•GT.ABSCX-XTCNYAL))) N = NVAL-?
NFLAG=2
G3 T0 25
3 N = NVAL-1
NFLAG=1
D0 10 I = 1,N
IF CX.NE-XTCI)) G0 T0 5
TBLUQ = YTCI)
RETURN
5 IF CCX.GT.XTCI)).AND.CX.LT.XTC1 + 1 ) )) G0 T0 15
IF CCX.LT-XTCI)).AND.CX.GT.XTCI+1))) G0 T0 15
10 C0NTINUE
N = NVAL-2
G0 T0 25
15 IF CI-3) 16,20,20
16 N = 1
G0 T0 25
20 N = I-1
25 XI = XTCN)
X2 = XTCN+1)
X3 = XTCN+2)
Yl = YTCN)
Y2 = YTCN+1)
Y3 = YTCN+2)
X1S = X1*X1
X2S = X2*X2
X3S = X3*X3
DET = X2S*X3-X3S*X2-X1S*X3+X3S*X1+X1S*X2-X2S*X1
A = CY2*X3-X2*Y3-Y1*X3+Y3*X1+Y1*X2-Y2*X1)/DET
B = CX2S*Y3-X3S*Y2-X1S*Y3+X3S*Y1+XIS*Y2-X2S*Y1)/DET
C = CY1*CX2S*X3-X3S*X2)-Y2*CX1S*X3-X3S*X1)+Y3*CX1S*X2-X2S*X1)
04860+ )/DET
04880 TBLUQ = A*X*X+B*X+C
04900 RETURN
04920 END
04940 SUBR0UTINE TURPIP
04960 DIMENSI0N N0DELC1),BTUPHRC1),RANGEC1), ITDC1)
04980 REAL ITD,N0DEL,N0LEV,NSECPL,NDELPL,NMSL,MPS,MPD,MVD, MPLEN
05000 REAL MPCST*MVCST,NMVAL,NHSL,NHV,NDV,NFV,NBPV,NEDL,NEDVAL
05020 C0MM0N JJ
05040 C0MM0N BTUPHR,RANGE,N0DEL,BDIAFT,TPCST,TVCST,PC0ST,FPCST,
05060+STCST,T0TAL*GPM
05080 DATA DELHI,PI/65.6,3.14159/
39
-------
05100
05120
05140
05 160
05180
05200
05220
05240
05260
05280
05300
05320
05340
05360
05380
05400
05420
05440
05460
05480
05500
05520
05540
05560
05580
05600
05620
05640
05660
05680
05700
05720
05740
05760
05780
05800
05820
05840
05860
05880
05900
05920
05940
05960
05980
06000
06020
06040
06060
06080
DATA GAPDEL,N0LEV>NSECPL/448.,2.>4./
D0 200 IR = 1*1
NDELPL = N0DELCIR)/N0LEV
JJ=1
GPM = BTUPHR(IR)/(RANGECIR)*8.337*60.
GPMPL = GPM/N0LEV
CALL PIPSIZ(GPMPL*NMSL,MPS,MPD,MVD>
NMSL*N0LEV
2.*400.*TNMSL+DELHI*(N0LEV-
MPLEN*MPD
2.*TNMSL
NMVAL*MVD
TNMSL
MPLEN = 2.*400«*TNMSL+DELHI*(N0LEV-1.)*NMSL
MPCST
NMVAL
MVCST
JJ = 2
HGPM = GPM/(N0LEV*NSECPL)
CALL PIPSIZ(HGPM,NHSL*HPS*HPD, HVD)
IF (NHSL.GT.1.)16* 20
16 PRINT,/>*T00 MANY HEADER LINES*
G0 T0 200
20 HPLEN = (BDIAFT + 3.55)*(1.+PI/2.)*N0LEV*NSECPL
HPCST = HPLEN*HPD
NHV = 4.*NSECPL*N0LEV
HVCST = NHV*HVD
DELPSE = N0DEL/
DGPM = DELPSE*GAPDEL
CALL PIPSIZ(DGPM,NDL,DPS,DPD, DVD)
DPLEN = BDIAFT*NSECPL*N0LEV
DPCST = DPLEN*DPD
NDV = 2.*NSECPL*N0LEV
DVCST = NDV*DVD
FGPM = DGPM/10.
CALL PIPSIZ(FGPM,NFL,FPS,FPD, FVD)
FPLEN = DPLEN
FPCST = FPLEN*FPD
NFV = NDV
FVCST = NFV*FVD
NBPV = NMVAL/2.
BPVALS = MPS
BPVALD = MVD
BPVCST = NBPV*BPVALD
EDGPM = GPM/TNMSL/4.
CALL PIPSIZCEDGPM,NEDL>EDPS, EDPD, EDVD)
NEDL = TNMSL
BDIAFT*.75*NEDL
EDPLEN*EDPD
EDPLEN :
EDPCST =
NEDVAL =
EDVCST :
IPSETS :
PC0ST =
NEDL
NEDVAL*EDVD
GPM/85000.+1.
IPSETS*300000.
IFPSET = FGPM/5500.+2.
FPCST = IF-PSET*20000.
40
-------
06 100
06120
06 140
06160
06180
06200
06220
06240
06260
06280
06300
06320
06340
06360
16912*GPM
12360.
STLBS*.5
MPCST+HPCST+DPCST+FPCST+EDPCST
MVCST+HVCST+DVCST+FVCST+EDVCST+BPVCST
TPCST+TVCST
TPVCST+PC0S7+FPCST+STCST
STLBS =
STCST =
TPCST =
TVCST =
TPVCST=
T0TAL =
200 CONTINUE
RETURN
END
SUBROUTINE PIPSIZCGPM,N0L,PDIA,DPRFT,DPRVAL)
REAL N0L,NSL
C0MM0N JJ
DIMENSI0N PSIZE(14),DPFTC2,14),DPVALC14)
DATA PS IZE/18«, 22.,24.,26.,32.,36.,42.,48 .,54.,60«, 72.,84.,96.,
06380+108./
06400 DATA DPFT/31., 45. ,32. ,49. ,33. ,51. ,36. ,57. ,40. ,64. ,42. ,69. ,52. ,84.,
06420+57.,93.,69., 109.,76., 121., 102., 155., 117., 178., 134.,203., 161-,238./
06440 DATA DPVAL/2960.,5350.,3720.,7330.,10480.,6380.,12670.,
06460+14520.,18650.,23650.,31280.,41370.,59340.,70920./
06480 NSL = 1.
06500 10 DIA = SQRTC3.55E-4*GPM/NSL)*12.
06520 IFCDIA.LE.108.)G0 T0 20
06540 NSL = NSL+1
06560 G0 T0 10
06580 20 D0 30 I = 1,14
06600 IFCDIA.LE.PSIZE CD) G0 T0 40
06620 30 CONTINUE
06640 I = 14
06660 40 N0L =NSL
06680 PDIA = PSIZECI)
06700 DPRFT = DPFTCJJ,!)
06720 DPRVAL = DPVALCI)
06740 RETURN
06760 END
41
-------
APPENDIX B
SAMPLE DATA
42
-------
TABLE 1-B
Auxiliary Power Requirements for 800-Mw Natural-Draft,
Dry-Type Cooling Tower
ITD (°F)
30
40
50
60
70
80
Fossil-Fueled Turbine (kw)
12,000
9,000
7,000
6,000
5,000
4,000
TABLE 2-B
Nuclear-Fueled Turbine (kw)
19,000
14,000
11,000
9,000
8,000
7,000
Auxiliary Power Requirements for 800-Mw Mechanical-Draft,
Dry-Type Cooling Tower
ITD (°F)
30
40
50
60
70
80
Fossil-Fueled Turbine (kw)
32,000
26,000
21,000
17,000
14,000
12,000
Nuclear-Fueled Turbine (kw)
48,000
39,000
32,000
26,000
22,000
18,000
43
-------
TABLE 3-B
ITD (°F)
30
40
50
60
70
80
ITD (°F)
30
40
50
60
70
80
pital Costs for 800-Mw Natural-Draft,
Dry-Type Cooling Tower
FOSSIL UNIT
Elevation
Sea Level
$32,610,000
22,880,000
17,200,000
14,090,000
12,040,000
10,280,000
NUCLEAR
Sea Level
$50,000,000
36,080,000
26,920,000
21,530,000
17,720,000
15,870,000
3,000 Ft.
$36,150,000
24,430,000
18,220,000
14,700,000
12,500,000
10,600,000
UNIT
Elevation
3,000 Ft.
$54,800,000
37,840,000
28,210,000
22,340,000
18,160,000
16,220,000
6,000 Ft.
$40,500,000
26,850,000
19,770,000
15,740,000
13,260,000
11,100,000
6,000 Ft.
$60,000,000
41,250,000
29,890,000
23,670,000
19,100,000
16,810,000
44
-------
TABLE 4-B
IIP (°F)
30
40
50
60
70
80
ITD (°F)
30
40
50
60
70
80
ral Costs for 800-Mw Mechanical-Draft,
Dry-Type Cooling Tower
FOSSIL UNIT
Elevation
Sea Level
$23,920,000
20,000,000
16,400,000
13,280,000
10,960,000
9,720,000
NUCLEAR
Sea Level
$35,280,000
29,480,000
24,120,000
19,640,000
16,160,000
14,240,000
3,000 Ft.
$24,460,000
20,450,000
16,770,000
13,580,000
11,210,000
9,940,000
UNIT
Elevation
3,000 Ft.
$36,070,000
30,140,000
24,660,000
20,080,000
16,520,000
14,560,000
6,000 Ft.
$25,120,000
21,000,000
17,220,000
13,940,000
11,510,000
10,210,000
6,000 Ft.
$37,040,000
30,950,000
25,330,000
20,620,000
16,970,000
14,950,000
45
-------
TABLE 5-B
"TURBIN" Data File
(Fossil Fuel)
00100
00110
00120
00130
001 40
0015 0
00160
00170
00180
00190
00200
00210
00220
00230
00240
00250
00260
00270
00280
00290
00300
00310
00320
00330
00340
00350
00360
00370
00380
00390
00400
0041 0
00420
00430
00440
0 0450
00460
00470
00480
00490
00500
005 1 0
00520
00530
00540
00550
00560
00570
00580
00590
00600
00610
00620
00630
00640
00650
00660
800000
4OOO.
1 . 0, 18
809975
790268
7*8560
7 1 3846
6841 50
1 .OE50
8 100.
3302.
8765.
9191.
9590.
9000.
9225.
9739.
10212
10655
3797.
3864.
4006.
41 25.
4226.
9
9
i
9
9
9
9
9
FOSSIL
.0
. , 809470. , 808966
. , 784068. , 77 7499
. , 743P49. , 738182
., 709 170.* 704555
.* 680209. * 676314
8105.
8368.
8827.
"252.
9645.
9006.
9298.
os 0*.
,
,
,
,
*
*
*
,
., 10280
•
9
9
9
9
9
,1071
3708.
3885.
4024.
4141-
4239.
7
9
9
9
9
9
81 10. ,81 30
8439. * 8509
8S~'8. * 8949
9312. , 9368
9701 . , 9756
901 l.,9034
9376. ,9455
0876., 9943
. , 1 0347. , 1
., 10779. , 1
3HOO. , 3S07
3907., 3929
4042. , 4059
/il 56. , 4! 7!
4253. y 4266
.,806917., 803960. ,800000. ,79 5 591.
., 771039., 765 1 3 b. , 759588., 754034.
., 733183., 728 ;\ bl., 72338 5., 718 58 4.
•>70 0377. ,696249. ,69 21 68., 688136.
. ,672533. ,668657. ,66 49 61. ,6611 70.
.,8161.
.,8575.
. ,9009.
. * 0423.
.* 0812.
. *9068.
. *9528.
. , 1001 0
04OO. , \
0840. , !
., 38 1 7.
. , 3950.
. , 4075.
. , '1 ! 8 b .
. , 4279.
*8201
*8638
*9070
*9/j79
* 9867
*91 12
*9597
.
•
.
.
.
.
.
,
,
,
,
,
,
,
. , 1 0078
0470-
0902.
,3831
, 39 68
, 4092
.-4199
, 429 1
,
,
.
.
•
•
.
1
1
,
,
*
*
,
8247.
8701.
9130.
9534.
9923.
9163.
9668.
. , 10145.
0532.,
0963.,
3846.
3987.
4108.
4212.
4304.
10594.
1 10?. 6.
1 . OEbO
8213.
8525.
8998.
9398.
9760.
9126.
9473.
9 9 9 H .
10442
1 0844
2 8 8 0 «
3067.
3351 .
3591 .
3808.
9
9
9
9
9
9
9
9
8218.
8600.
9057.
o 4 5 2 .
0809.
9132.
o 5 b 6 •
1 0063
3
9
9
9
9
9
9
•
. , 1 ObO?
. 9 1 0899
9
9
9
9
2883.
3112.
3386.
3623.
3838.
i
9
9
9
8238. , 827 1
8673. , 8743
9 1 15. * 9 1 74
9506. * 9560
0850. , 9908
9 I 53. ,9)90
o 6 3 7 . , o 1 \J>
, 1 0 ! 28 . , ! 0
. , !0b62. , 1
. , 10954. , 1
2895 . , PO i 5
3156., 3108
34?] . * 3456
3656. , 3688
3«67., 3897
.,831 1 .
., 881 1 .
.,9231.
. * 9607.
. * 0955.
. * 9 2 3 5 .
. , o 7 2659.
.,2832.
.,2985.
, °01 4
,9632
.,101
0493.
0876.
, ! 00 !
061 6.
1 ! 8 o .
1 65o.
2094.
, PP'if)
, 2488
,2685
,2855
, 3006
•
•
,
,
26
9
9
6
,
9
,
9
•
•
•
•
•
1
1
.
1
1
1
1
,
,
,
,
,
91 l&.
9709.
., 1019
0550.,
0929.,
, ! 0 ! P8
0703.,
1251 .,
' 723.,
21 44.,
2281 .
2518.
271 1 .
2878.
3026.
1.
10608.
10979.
.
1078H.
1 1323.
« 1 787.
12199.
1.0E50
46
-------
TABLE 6-B
"TURBIN" Data File
(Nuclear Fuel)
00100
00110
00120
00130
00140
00150
00160
00170
00180
00190
00200
00210
00220
00230
00240
00250
00260
00270
00280
00290
00300
00310
00320
00330
00340
00350
00360
00370
00380
00390
00400
00410
00^20
00430
00440
00450
00460
00470
00480
00^90
00500
00510
00520
00530
00540
00550
00560
00570
00580
00590
00600
00610
00620
00630
00640
00650
00660
800000. NUCLEAR
6000.
2. 0, 18. 0
81 36 87., 8 I 3687., R 136 87
778425., 768R44., 759427
7 17688., 7 10078., 704520
677 381. ,672321. ,667869
6 486 17., 64 49 68. ,641 410
.,811636.,807175.,80000P«, 789884.
., 750978. , 7427 1 5., 733924., 726030.
., 698465. , 692455., 68 7169., 68 19 62.
.,663475.,6b9710.,654910.,65174R.
.,63R377., 6353 73. ,63187]., 62 7937.
.OE50
0390.
0861 .
1780.
2481 .
3034.
0390.
CR61.
1780.
24R1 .
3034.
5677.,
5797.,
6005.,
6 1 42. ,
624Q. ,
, 10390. , 10390. , 1
, 10996
, 1 1906
, 12574
, 131 08
, 10390
, 10996
, 1 1906
.,11 132. , 1
• , 12000. , 1
. , 12658. , 1
• , 1 3180. , 1
., 1 0390. , 1
.,111 32. , 1
• , 12000. , 1
, 12574. , 1265R. , 1
, 131 OR
5677. ,
5830.,
6031 .,
6] 59.,
6253-,
., 1 3180. , 1
5677. , 5684
5P62., 589 1
6050. , 6070
61 75., 6190
6265., 6275
0416. ,
1257. ,
2104. ,
2742. ,
3243. ,
0416.,
1257. ,
2104. ,
2742.,
3243.,
. ,5699
. , 5919
.,6091
. ,6203
. ,6286
10474. ,
1 1 383. ,
12209. ,
12R15.,
1 3306. ,
10474. >
1 1383. ,
12209. ,
12815.,
1 3306. ,
.,5724.
. , 5949i
.,6109.
.,6219.
., 6298.
0568.,
1519.,
2303.,
2909.,
3379.,
0568.,
1 1519.,
12303.,
12909.,
13379.,
,5758.
, 5976.
,6127.
,6230.
, 631 1 .
10703.
1 1644.
12397.
2971.
3463.
0703.
1644.
2397.
2971.
3463.
1 .OE50
10535.
1 1 5.2 .
2093.
2761.
3264.
0535.
1 152.
12093.
12761 .
1 3264.
4273.,
4644 . ,
5208.,
5609.,
591 1 .,
.OE50
1 085.
2051.
31 49.
3924.
4490.
1085.
2051-
31 49.
3924.
4490.
3069.,
3455. ,
3894.,
4204. ,
4431 • »
1 . OE50
, 0535
, 1301
, 2215
, 2846
, 3330
, 0535
, 1 301
, 1 221 5
, 128^6
, 1 3330
4273.,
4733.,
5281 .,
5660. ,
5950.,
, 1 112R
, 12221
, 13293
, 1 4020
, 1 4564
, 1 1 128
, 12221
, 13293
, 14020
, 1 4564
3086.,
3523.,
3952. ,
4243. ,
4460. ,
. , 10552., 1
., 1 1 441 . ,
., 2307.,
., 2922.,
. , 3395. ,
., 0552.,
., 1 44] . ,
., 2307.,
., 2922.,
. , 33^5. ,
428/i. , 431 o
4R 1 7. , 4R92
5336. , 5398
57CS.,57b4
5989., 6022
. , 1 1 JPM . , 1
. , 12384.,
., 13398.,
. , 1 4103.,
. , 1 4637. ,
., 1 H81. ,
. , 12384. ,
., 1 339 H. , !
., 1 4103. , 1
. , 1 4637. , 1
31 07. , 31 44
3588. , 3648
3994. , 4042
4276. , 431 2
4490. > 451 3
0596.,
1 5fc>7 . ,
2409.,
3003. ,
3450. ,
C596. ,
1567.,
2409. ,
3003.,
3450. ,
. , 4352
. , ^970
. , 5458
., 579 1
. , 6051
1273.,
2532.,
3518.,
4193.,
4695. ,
1273. ,
2532. ,
3518. ,
4193.,
4695. ,
.,3210
. , 3708
. , 4089
. , 4343
. , 4534
10667. ,
1 1697.,
12510.,
1 3065. ,
1 3499. „
\ 0667. ,
1 1697.,
12510.,
1 3065.,
1 3499. ,
. , 442b.
.,5053.
.,5512.
.,5847.
. , 6093.
1 1437.,
12684. ,
13635.,
1 4269. ,
1 /i748. ,
1 1437.,
1 2684. ,
13635. ,
1 /I ,0 f; 9 . ,
] 474P. ,
.,3281.
.,3773.
. ,^129.
.,4382.
. , 4563.
10789.,
1 183b. ,
12600.,
13157.,
13567.,
10789.,
1 1835.,
12600. ,
13157.,
1356 7. ,
,4532.
,5129.
,5565.
,5879.
, 6134.
11615.,
12846. ,
3736. ,
4368.,
4821 .,
1615.,
2846.,
3736.,
4368.,
1 4*21 . ,
,3367.
,3832.
, 41 70.
, 440R.
,4593.
10967.
1 1961.
1268R.
1321 1 .
13637.
10967.
1 1961 .
12688.
1321 1.
13637.
1 1830.
12994.
13837.
14432.
14895.
1 1830.
12994.
13837.
14432.
14895.
47
-------
TABLE 7-B
"SITEXX" Data File
100 DENVER COLORADO
110 5300
120 SUMMER
130 3,20.,30.,35.
140 5,8.,10.,12.,15.,18.
150 .95
160 100.
170 8.,10.,12.,15.,18.
180 40.,40.,40.
190 0.0
200 150.
210 1.0
220 O.,0.,0.,.01,.11,1.18,2.69,3.79,4.98,6.26,7.80,8.93
230 8.34,7.73,8.03,7.89,8.18,8.22,6.31,4.09,2.46
240 1.36,.89,.41,.25,.07,.01,0.,.01,0.,0.,0.
250 1 .OE50
TABLE 8-B
"SIXDAT" Data File
57., 28.5, 5.0E09, 220400., 50., 600.
57., 28.5, 5.0E09, 220400., 75., 600.
57., 28.5, 5.0E09, 220400., 95., 600.
62., 31.0, 5.0E09, 220400., 50., 600.
62., 31.0, 5.0E09, 220400., 75., 600.
62., 31 .0, 5.0E09, 220400., 95., 600.
48
-------
APPENDIX C
OUTPUT SAMPLES
49
-------
TABLE 3-C
Tower Sizing Program "TOWSIZ"
Dry Cooling Tower Sizing and Cost Evaluation
Design Parameters
1TD = 40°F
Heat Rejection = 4.0E+09 Btu/hr.
Water Flow per Hour = 2.2E+05
Ambient Air Temp = 50°F
Range = 20°F
631.4ft.
421.8ft.
671.7ft.
399824
Tower Sizing
Tower Height =
Upper Diameter =
Bottom Diameter =
Gallons per Minute
Cost Evaluation
Tower Structure
Stack Cost
Shed Cost
Coil Cost
Total Structure
Condenser
Condenser Cost
Piping, Valves, etc.
Pipe Cost
Valve Cost
Pump Cost
Filler Pump Cost
Storage Tank Cost
Total Piping Facilities
Controls
Control Cost
Complete Tower Facilities
Total Tower Cost
Total Tower Cost and Contingencies
Elevation = 0 ft.
$3311729
1631071
6583000
$2383132
1309700
1800000
40000
39989
$11525799
$ 912000
$ 5362921
$ 500000
$18300721
$22875901
50
-------
TABLE 1-C
Program OPTDCT Output
SITE - DENVER, COLORADO
ELEVATION - 5300 FEET ABOVE SEA LEVEL
PEAK PERIOD - SUMMER
TURBINE - 800000 KW FOSSIL UNIT
IN OPERATION - 7500 MRS/YR
AT FULL THROTTLE 50 0/0 OF THE TIME
AT 75 0/0 LOAD 50 0/0 OF THE TIME
AT 50 0/0 LOAD 0 0/0 OF THE TIME
TOTAL ANNUAL ENERGY - 5250000 MWH NOMINAL
TOWER - NATURAL DRAFT
COOLING SYSTEM OPERATION AND MAINTENANCE COST - 1 .00 0/0
DESIGN HEAT REJECTION - 4000 X 10**6 BTU/HR
CAPITAL COST MULTIPLIER - .95
CUTOFF TEMPERATURE FOR PEAKING GENERATION -OF
51
-------
TABLE 2-C
Program OPTDCT Output
Ol
NO
CAPITAL COST FACTOP.SI PLANT - i ttt
IM.ANT FUEL COST - 20 CENTS/10"6 BTU
PEAKING
INIT.
TEMP.
OIFF.
« F >
30
31
33
3I|
35
36
37
!»
39
ill
42
43
1,5
47
48
49
51
53
51,
59
OPTJHuHl
56
58
59
60
62
64
66
67
68
69
70
71
72
73
74
75
76
77
7«
79
en
CAPITAL COST - 10
GROSS
ENERGY
800 NH
CHHHI
5281,268
5283562
5213164
52827"i6
5282280
5281763
5281208
5280616
5279287
5277700
5275892
S273797
5272630
5271398
5268699
5265621
5263957
5262208
5260 358
5256322
5254154
5251891
5246998
524161 "i
5235767
5233639
52293611
52221,3 1>
5218761
5211,939
5206822
5202557
5198139
5193561
5188826
5183939
5178903
EXCESS
ENERGY
DUE TO
EXTRA
IHHHI
3",351,
33739
33395
33051
32727
32360
31964
31535
30712
29797
28735
27801
27258
26677
2551,5
21,1,23
23793
23118
22563
21440
20804
20125
1901,7
17895
16665
16201
15671,
11,1,92
13881
13432
12353
11749
11106
10678
10189
9655
9080
0 1/KH
PEAKING CAPACITY -
PEAKING FUEL COST -
9 0/0 AUXILIARIES - « 0/0
40 CENTS/10"6 BTU
AUXILIARY CAPITAL COST - 150 I/KH
CAPACITY HAXIHUN
PENALTY AUXILIARY LOSS OF AUXILIARY
IHHHI
86
177
226
306
41,7
597
754
919
1425
2097
2843
4004
4628
5279
6847
8802
9836
10909
12205
15117
1665J
18234
20038
22049
26281
28499
30898
33561
36310
42858
45120
48493
55530
5919*
E2966
67117
71363
75716
80177
(HHHI
90624
85401
•2929
80550
78264
76071
73971
71964
68229
64666
61875
59256
58086
57009
"56025
54878
52724
51717
50756
49842
4815»
4/378
46650
45099
43656
42975
42321
41694
4109%
39975
39402
38847
38309
37788
37214
36798
36329
35«77
35442
35025
(KK)
6636
9199
106G5
12095
13677
15354
17108
18912
22742
26672
30696
34662
36641
38627
42689
47054
49297
51585
53899
58586
60941
63327
68286
73432
78818
81571
34409
90253
93236
96253
99114
101953
104879
107874
110906
113986
117161
120375
61569
85692
113110
128515
144837
161919
179503
217779
257967
299465
342909
387045
433822
485706
51Z55D
540037
569101
629692
660558
692010
725437
760526
833261
171119
910515
951824
994421
1082674
112(492
1176485
1223642
1271148
1320040
1370237
1423034
1476(50
1532178
1588507
III
318321
300008
283003
274991
267307
252921
239843
228077
217620
208476
200643
193223
185725
192223
178885
175712
169(61
167183
164669
161943
159312
154341
COST OF COOLING
PLANT FUEL
CO (NILL/KHH)
13128377 2.5006
12877647 2.4529
12653514 2.4102
12456902 2.3727
12286708 2.3403
12143357 2.3130
12025338 2.2905
11932196 2.2728
11964091 2.2598
11820324 2.2515
11785625 2.2449
11753228 2.2387
11742157 .2.2366
11734637 2.2352
11731217 2.2345
11735248 2.2353
11/42164 2. 2364.
11752521 2.2386
11770998 2.2421
117.99860 2.2476
149755 11(39467 2.2551
145554
141741
1397(9
137900
136074
134309
1J2606
130966
1293((
127(74
126421
125031
11891218 2.2650
11953792 2.2769
11983366 2.2(2*
12015346 2.2886
12047(25 2.2947
12079623 2.3009
1211417b 'J.iarb
12150619 2.3144
17189920 2.3219
12230734 2.3297
12273697 2.3379
1231(188 2.3463
-------
APPENDIX D
FLOW CHARTS
53
-------
FLOW CHART OF PROGRAM "OPTDCT1
START
OPTDCT
V
ASSIGN STORAGE
DEFINE CONSTANTS
AND TABLES
READ TURBINE
CHARACTERISTICS
FROM^FILE "TURBIN"
\(
CALL DATCHK
TO CHECK IF CORRECT
NUMBER OF DATA
VALUES READ
CALCULATE
NOMINAL ANNUAL
MWH
54
-------
ISITE 1
READ SITE INFORMATION
FROM FILE SITENO
(ISITE)
SET WINTER OR
SUMMER PEAK
FIND FIRST NCN-ZERO
AND LAST NCN-ZERO
TEMPERATURE DURATION
IN TABLE OF DURATIONS
TDUR (NTL> AND TDUR (NTU)
CALL DATCHK
55
-------
OUT-PUT HEADINGS
IDENTIFYING SITE, TYPE
OF TOvVER, TYPE OF
TURBINE AND
OTHER INFORMATION
CALCULATE THE K ,V LOAD
FOR FULL THROTTLE, 3/4 LOAD
AND l/'i LOAD AT THE
MAXIMUM OPERATING BACK
PRESSURE FOR LATER USE
IN THE INTERPOLATION
OF TABLES (LOAD)
11
w
CALCULATE THE KW LOAD
FOR 3/4 LOAD AND 1/2 LOAD
AT THE CALCULATED BACK
PRESSURE (LOAD1)
12
DETERMINE HEAT REJECTION
AT VARIOUS LOADS BY
INTERPOLATION
13
Q
56
-------
COMPUTE HOURS IN TH"EE
HIGHEST TEMPERATURE
DURATIONS
14
DETERMINE FIRST TEMPERATURE
FOR TEN HOURS DURATION
15
1TD = 30
0-
16
THETA ITD
I = ITD -Z9
DETERMINE CAPITAL COST
BY INTERPOLATION
18
MODIFY CAPITAL COST
BY CAPITAL COST FACTOR
19
57
-------
G>
CALCULATE A IN
ITD = A* (HEATREJ)2
u
SET LOOP CONTROL
FOR LOADS
V
SET LOOP CONTROL
FOR TEMPERATURES
20
21
INITIALIZE ITERATION COUNTER
(ITER = 0) AND BACK PRESSURES
(BPT = O, OPBP = BPMIN)
O
G>
CALCULATE COOLED WATER
TEMPERATURE
24
58
-------
FIND CORRESPONDING
PRESSURE FROM STEAM TABLES
INTERPOLATE TO FIND
HEAT REJECTION
AT BACK PRESSURE
ABS(OPBP-BPT) .LT. 0.01
59
-------
30
PRINT ERROR
MESSAGE
OPBP - BPM1N
34
YES
\
COMPUTE TOTAL ANNUAL
HOURS (HRS) AT GIVEN
LOAD AND TEMPERATURE
35
60
-------
INTERPOLATE TO FIND MAX.
KW AT CALCULATED BACK
PRESSURE
37
COMPUTE PLANT OUTPUT
SUBJECT TO MAX. CALCULATED
ABOVE
38
©*
©-
COMPUTE MAXIMUM HEAT
REJECTION
39
y
r^~*
8A
61
-------
CALCULATE REDUCED PLANT
OUTPUT BY TABLE LOCK-UP
40
INTERPOLATE TO FIND
STATION HEAT RATE
41
ACCUMULATE TOTAL BTU's
USED FOR A YEAR'S OPERATION
42
NATURAL OR
MECHANICAL DRAFT
TOiVER
TABLE LOOK UP TO FIND
MAXIMUM AUXILIARY
REQUIREMENTS .
62
-------
V
45
LINEAR INTERPOLATION TO GET
ACTUAL AUXILIARIES AT GIVEN
ITD AND AMBIENT AIR
TEMPERATURE
COMPUTE AUXILIARY KWH
(AUXKWHj)
46
DETERMINE AUXILIARY POWER
(AUXKWHj) REQUIREMENTS AS A
FUNCTION OF ITD BY TABLE
LOOK-UP
47
COMPUTE AUXILIARY ENERGY
REQUIREMENTS (AUXKvVHj)
48
49
ACCUMULATE TOTAL ENERGY
PRODUCED BY PLANT (KWH)
63
-------
COMPUTE. CAPACITY LOSS
OR GAIN (CAPLOS)
50
COMPUTE ENERGY GAIN IN
KWH AND BTU's
55
/ NO
AIRTNJ .GT. CUTOFF
ACCUMULATE ENERGY
LOST IN KWHBLOT
YES \
53
ACCUMULATE ENERGY
LOST IN KWHABV
STORE KW LOSS FOR THREE HIGHEST
AIR TEMPERATURES IN KWLOS
56
64
-------
NT
NT +1
57
YES
61
SAVE TOTAL ANNUAL ENERGY
PRODUCED IN MvVH,
62
DETERMINE CAPACITY LOSS
AT AIR TEMPERATURE EQUALLED
OR EXCEEDED 10 HRS-/.YR.
65
-------
WRITE HEADINGS ONTO
OUTPUT FILE
\
f
ITD = 29
TANT 1.OE50
I = 1
\
r
ITD = ITD + 1
68
69
66
-------
V
CALCULATE ANNUAL CAPITAL COST
AS A FUNCTION OF TOTAL
CAPITAL COST
70
CALCULATE ANNUAL OPERATION
AND MAINTENANCE COST. ADD TO
ANNCAP TO GET TOTAL
ANNUAL COST ON CAPITAL
INVESTMENT (ACCAOM)
71
CALCULATE ANNUAL
FUEL COST
72
WINTER OP
SUMMER PEAK
COMPUTE CAPITAL
AND OPERATING
PENALTY COSTS
COMPUTE OPERATING
COST PENALTY
COMPUTE TOTAL
PENALT-Y COST
76
67
-------
COMPUTE TOTAL
ENERGY LOSS
77
78
COMPUTE AUXILIARY CAPITAL
COST, FUEL COST AND
OPERATING COST, AND TOTAL
ANNUAL COST
79
SUM UP ALL
ANNUAL COSTS
WRITE TEMPORARY
OUTPUT FILE
68
-------
YES
COPY TEMPORARY OUTPUT
FILE ONTO ACTUAL
OUTPUT FILE
85
-------
70
-------
FLOW CHART OF PROGRAM "TOWSIZ"
DEFINITION OF CONSTANTS:
ALL, VALVES TO BE USED
FOR TABLE LOOK UP
PROCEDURES IN PROGRAM
EXECUTION
1. 01
SCALE TABLE
ENTRIES
ERMINATE\
RUN /
1. 00
(
_
1.02
TOWKR
SIZING
YES
-N
)
/
1. 03
1.04
1.08
NO
READ SITE DATA:
FROM ONE LINE OF FILE,
"SIZDAT" READ ITD,
RANGE HTREJ, G ,
AMAIRT AND ELEV
COMPUTE INTERMEDIATE SIZING:
CONDENSATE FLOW, NUMBER OF
COOLING COLUMNS (ROUND TO
NEXT INTEGER DIVISIBLE BY 4),
Q AND Q.
1.05
WRITE
ERROR
•1ESSAG
1-07 NO CURVE
CORRESPONDING
TO ENTERED
7
LOOP TO
DETERMINE WHICH
, CURVE IS BEING USED
CORRESPONDING
G0 CURVE FOUND
71
-------
PRINT
1.11
/Ql OUTSIDE
,/OF ENTERED
~Ql TABLES
1. 10
CALL ROUTINE "TBLUQ" TO
FIND LQ BY QUADRATIC
INTERPOLATION USING Qj
AND TABLES Qj AND LQ
VALUES.
v
NORMAL EXIT
1. 12
CALL "TBLUQ" TO FIND&P
BY QUADRATIC INTERPOLATION
IN ENTERED TABLES.
1.13
COMPUTE EXIT VELOCITY,
EXIT LOSSES, DRAFT
LOSSES, TOTAL LOSSES
1.14
COMPUTE AIR HEATING,
AIR TEMPERATURE
CHANGE AND FINAL
AIR TEMPERATURE
1.15
COMPUTE TOWER
DRAFT HEIGHT
72
-------
CALL "TBLUQ" TO FIND
HEIGHT ADJUSTMENT FACTORS
FOR AMBIENT TEMPERATURE
AND SITE ELEVATION
\
/
COMPUTE ADJUSTED DRAFT HEIGHT
FROM DRAFT HEIGHT AND
ADJUSTMENT FACTORS AND TOTAL
HEIGHT FROM ADJUSTED
HEIGHT PLUS 80 0 FEET
\
J
CALL "TBLUQ" FOR EXIT
VELOCITY ADJUSTMENT
IN UPPER DIAMETER
CALCULATION.
i
f
COMPUTE UPPER STACK AREA,
UPPER DIAMETER AND
LOWER DIAMETER
1.16
1.18
1.19
SET LOWER
DIAMETER EQUAL
TO UPPER
DIAMETER
1.31
l.ZO
LOWER DIAMETER
LESS THAN
UPPER DIAMETER?
73
-------
.1.23
NO SUCH
ENTRIES
FOUND
FIND
NTRIES IN TABL
F STACK DIAMETERS WHIC
LCULATED UPPER DIAMETE
FALLS BETWEEN
ENTRIES FOUND
CALL "TBLUQ" TO COMPUTE COST
OF STACK AT CALCULATED HEIGHT
FOR EACH OF TWO LIMITING
TABLE DIAMETERS
1.24
COMPUTE ACTUAL STACK
STRUCTURE COST BY
LINEAR INTERPOLATION OF
DIAMETER BETWEEN TWO
CALCULATED TABLE DIAMETERS
1.25
COMPUTE ROOF COST BY AREA
DIFFERENCE BETWEEN STACK
TOP AND BOTTOM. COMPUTE
TOTAL COIL COST
1.26
CALL"TBLUQ" TO FIND CONDENSE!
COST PER 109 BTU FROM
RANGE VALUE IN TABLE OF
CONDENSER COSTS. COMPUTE
TOTAL CONDENSER COST
1.27
74
-------
COMPUTE NUMBER
OF
DELTAS
\
p
CALL "TWRPIP" SUBROUTINE
TO COMPUTE ALL PIPING,
VALVE AND ASSOCIATED
COSTS
\
f
1.28
1.29
1.30
PRINT ONE PAGE
SIZE AND COST
SUMMARY
75
-------
FLOW CHART OF FUNCTION "TBLUQ'
START
TBLUQ
IS LOOK UP
VAL UE WITHIN LIMITS
OF GIVEN TABLE?
IS GIVEN
VALUE BELOW LOWEST
TABLE ENTRIES?
IS GIVEN
VALUE ABOVE HIGHEST
TABLE ENTRIES?-
SET RETURN
VALUE TO
CORRESPONDING
INTERPOLATION
TABLE ENTRY
IS ENTERED
VALUE EQUAL TO
THIS TABLE ENTRY?
ERMINAT
TBLUQ
;SET LOOP
THRU ALL
TABLE
ENTRIES
76
-------
YES
^. YES
IS ENTERED
VALUE GREATER THAN
THIS TABLE ENTRY AND
ESS THAN NEX
2.09
IS ENTERED
VALUE LESS THAN
THIS TABLE ENTRY AND
GREATER THAN
NEXT?
^"
NO
2. 10
WITHIN
FIRST TWO TABLE
ENTRIES ?
2. 12
2. 11
INTERPOLATE WITH PRECEDING,
SUCCEEDING AND SECOND
SUCCEEDING VALUES.
EVALUATE QUADRATIC
INTERPOLATION
EXPRESSION AND SET
RETURN VALUE
E.15
(TERMINATED
V TBLUQ J
77
-------
FLOW CHART OF SUBROUTINE "TWRPIP"
c
START TWRPIP
CALCULATE WATER FLOW RATE
TO TOWER AND WATEH FLOW
RATE PER LEVEL
3.01
CALL ROUTINE "PIPSIZ"
TO CALCULATE MAIN SUPPLY LINE
AND THE COST/FT PLUS
COST/VALVE
3.02
V
CALCULATE THE MAIN SUPPLY LINE
LENGTH AND COST, NUMBER OF
VALVES AND COST AND
NUMBER OF LINES
3.03
CALCULATE HEADER PIPE
WATER FLOW RATE
3.04
CALL ROUTINE "PIPSIZ" TO
CALCULATE HEADER PIPE SIZE,
ITS COST/FT AND CORRESPONDING
VALVE COST AND THE NUMBER
OF LINES
3.05
78
-------
0
CALCULATE PIPE LENGTH, PIPE
COST, NUMBER OF VALVES
AND VALVE COSTS
3.06
\f
3.07
CALCULATE THE NUMBER OF
DELTAS PER SECTOR
3. 08
CALCULATE DRAIN PIPE
WATER FLOW RATE
V
3.09
CALL PIPSIZ TO CALCULATE THE
DRAIN PIPE SIZE AND COST/FOOT.
THE NUMBER OF VALVES
AND COST PER VALVE
3.10
CALCULATE PIPE LENGTH AND
COST THE NUMBER OF VALVES
AND COST
79
-------
CALCULATE THE FILLER PIPE
WATER FLOW RATE
3. 11
CALL "PIPSIZ" TO CALCULATE THE
FILLER PIPE SIZE AND COST/FT
AND THE NUMBER OF VALVES AND
COST/VALVE
3. 12
V
CALCULATE THE LENGTH OF PIPE
AND ITS COST AND THE NUMBER
OF VALVES AND THEIR COSTS
3. 13
CALCULATE THE NUMBER OF
BY-PASS VALVES AND COSTS
3. 14
CALCULATE THE EMERGENCY
DRAIN PIPE WATER FLOW RATE
3.15
CALL PIPSIZ TO CALCULATE THE
EMERGENCY DRAIN PIPE SIZE AND
COST/FT AND THE COST PER VALVE
3.16
80
-------
3. 17
CALCULATE THE LENGTH OF
EMERGENCY DRAIN PIPE AND ITS
COST AND THE NUMBER OF VALVES
AND THEIR COSTS
CALCULATE THE NUMBER OF
CIRCULATING PUMPS AND THEIR
COST
3. 19
CALCULATE FILLER PUMP COST
3.ZO
CALCULATE STORAGE TANK COST
3.21
SUM UP PIPING, VALVES, PUMPS,
AND STORAGE TANK COSTS
V
c
TERMINATE TvVRPIP
81
-------
FLOW CHART OF SUBROUTINE "PIPSIZ"
CALCULATE PIPE SIZE
4.01
FIND THE CORRESPOND PIPE COST
AND VALVE COST
4.04
V
TERMINATE PIPSIZ
82
-------
APPENDIX E
GLOSSARY
83
-------
Glossary of Terms Used in the
Cooling Tower Optimization Program
"OPTDCT"
A - The coefficient of the equation
ITD = A x (Heat Rejection 106 Btu/hr)Z
which defines the operating characteristics of a dry cooling tower.
ABS - Absolute, Library function of the computer.
ACCAOM - Annual capital cost and operating and maintenance costs.
AIR - Variable used to store the three highest non-zero ambient air tem-
peratures for each site.
AIR10 - Variable used to store the air temperature which is equalled or
exceeded ten hours per year. (Sometimes called "the ten hour air
temperature".)
AIRT - An array that contains air temperatures corresponding to the table
of durations in the site data . Values from 117°F to — 38°F in 5*^
increments (117°, 112°, 107°, . . .-28°, -33°, -38°).
AMAXl - Maximum, Library function of the computer.
AMIN1 - Minimum, Library function of the computer.
ANNCAP - Annual capital cost determined by applying the appropriate (fixed-
charge rate) to the total capital cost for the given ITD.
ANNFUL - Annual fuel cost determined by applying the appropriate unit cost
of fuel to the total Btu's used for the given ITD.
ASTER - Variable used only in output format indicating a number raised to
a power.
AUX100 - Table of full auxiliary power requirements vs. ITD for mechanical-
draft towers.
AUXCAP - The capital cost of auxiliary power required in dollars.
84
-------
AUXCST
AUXKW
AUXKWH
AUXOPR
AUXTEM
AUXTMP
BP
BPMAX
BPMIN
BPT
CAPCST
CAPLOS
The total annual cost for the auxiliary power and energy necessary
to the cooling system.
The maximum power needed for the cooling system auxiliaries
(pumps and fans).
Variable used to store the energy value for the auxiliary require-
ments of the cooling system.
Annual operating cost incurred in supplying the necessary auxiliary
energy.
Variable used to store the maximum and minimum air temperature
used in calculating the percentage of full auxiliary power require-
ments for a mechanical-draft tower for a given ITD.
Table of air temperatures corresponding to the maximum (100%) and
minimum (53%) percentage of full auxiliary power requirements for
a mechanical-draft tower.
Variable used for temporary storage of capital costs for a given
ITD corresponding to elevations of sea level, 3,000 and 6,000 feet.
An array which contains a table of back pressures which cover the
range of operation from 1 .0 inch Hg to 18.0 inches Hg in steps of
.5 inch Hg. Corresponds to values of full throttle kw, station heat
rate, and heat rejected to the cooling system. From the turbine
data file.
Maximum allowable turbine back pressure (inches Hg).
Minimum allowable turbine back pressure (inches Hg).
Variable used for temporary storage of the calculated back pressure
used in the calculations for finding the operating back pressure of
the turbine.
Variable used for storing and calculating the capital costs of the
cooling system
Variable used in the calculation of the loss or gain of capacity
(kw) resulting from turbine operation above or below a back pres-
sure of 3-1/2 inches Hg.
85
-------
CCF
CCFGT
CCM
CPWAX
CRDBTU
CREDIT
CUTOFF
DUR
ELEV
FCSTAX
FTKW
FTKWT
FULAUX
FULCST
The capital cost factors (fixed-charge rates) applied to the cooling
system capital cost represents annual cost.
Capital cost factors (fixed-charge rates) applied to the capital cost
of peaking units used to replace capacity losses.
Capital cost multiplier. This variable is used to adjust for varia-
tions in construction cost.
Capital cost (in $/kw) of providing auxiliary power. ($150/kw for
fossil-fueled units and $225/kw for nuclear-fueled units).
Array used to retain the Btu of excess energy. Excess energy is
produced only under full throttle operation at back pressures below
3.5 inches Hg and in the energy resulting from power production in
excess of 800 mw. One value per given ITD.
Annual credit, in dollars, determined by applying the appropriate
fuel cost CRDBTU.
The air temperature at or above that at which the lost power and
energy is replaced by a peaking unit.
Variable used to retain the first three non-zero temperature dura-
tions corresponding to the three highest temperatures for each site.
Elevation (feet above sea level) of each site being studied.
Variable used for the incremental cost (mills/kwh) auxiliary energy.
This cost is the annual plant fuel in mills/l
-------
GHR
HREJ
HREJD
HREJ MX
HREJT
HRJ
MRS
HRSPYR
IJ, IJ2
Ul
IK
IL
I SAVE
Plant heat rate with a value corresponding to each value of back
pressure.
Heat rejected by the turbine with a value corresponding to each
value of GHR.
Design heat rejection of the turbine (4 x 10^ Btu/hr for a fossil unit
or 6 x ICr Btu/hr for a nuclear unit).
Maximum heat rejection obtainable when the turbine is backed
down to operate within its prescribed limits.
Heat rejection associated with the maximum allowable back pres-
sure for the three operating levels of full throttle, 3/4 load and
1/2 load.
Temporary value of heat rejected by the turbine for use in determin-
ing the operating back pressure.
Number of hours of operation for a given load (PCTTIM) at a given
air temperature (TDUR) based on the number of operating hours per
year (OPHRS).
The number of hours in a year. Assumed to be 8,760.
Variable subscript corresponding to values of ITD (1 = 1, ITD = 30;
I = 2, ITD = 31; etc.).
Variable subscript corresponding to curves of percent of full auxil-
iary power requirement vs. ambient air temperature for different
ITD for a mechanical-draft tower.
Variable subscript corresponding to the curve IJ + 1 of percent of
full auxiliary power requirement vs. air temperature for a mechan-
ical-draft tower.
Variable subscript associated with the capital cost curves corres-
ponding to elevations of sea level, 3,000 feet, and 6,000 feet.
Subscript used to cycle through the three combinations of percent
load (PCTLD) and percent time (PCTTIM).
Variable used to store the optimum ITD for each run.
87
-------
1SITE
IT
ITD,THETA -
ITER
ITYPE
IWSP
KW
KWAX
KWH
KWHABV -
KWHBLO -
Variable subscript used to denote the site on which the program is
operating.
Variable subscript used to associate the loss or gain of capacity with
the first three non-zero temperature durations.
Variable used interchangeably to indicate the initial temperature
difference which is the difference between the saturation tempera-
ture of the condensing steam and the ambient dry-bulb air temper-
ature .
Variable used as a counter to keep track of the number of iterations
needed to converge the operating back pressure.
The name of the basic type of fuel used for the turbine, either
fossil or nuclear.
Variable indicator relating to the peaking period for the site.
Peaking period being either winter (IWSP = 1) or summer (IWSP =
2).
Variable subscript used to denote the appropriate capital cost fac-
tor (fixed-charge rate) and capital cost factor for gas turbine asso-
ciated with each run.
Variable subscript used to denote the appropriate plant fuel cost
and peaking fuel cost associated with each run.
Turbine nameplate rating in kw.
Array of auxiliary power requirements for the natural-draft tower.
Six values corresponding to ITD values of 30°F, 40°F, 50°F, 60°F,
70° F and 80° F.
Variable used twice, first for calculating the nominal energy pro-
duced by the plant (kwh) and secondly for calculating the energy
produced under given conditions (kwh).
The total energy loss occurring when the ambient air temperature
is equal to or above the cutoff temperature. One value for each
ITD (kwh).
The total energy loss occurring when the ambient air temperature
is below the cutoff temperature. One value for each ITD (kwh).
88
-------
KWHGAN -
KWHLOS -
KWLOS
KWLOSS -
LOAD
LOAD1
MILPKW
MWH
NCCF
NFUL
NHR
NT
NTL
The excess energy produced when operating at full throttle and be-
low 3-1/2 inches Hg back pressure. One value for each ITD (kwh).
The energy lost when operating at full throttle and above 3-1/2
inches Hg back pressure. One value for each ITD (kwh).
Variable used in calculating the loss of capacity (KWLOSS) at the
ten hour air temperature. One value, in kw, for each of the thre,e
highest air temperatures.
The lost capacity associated with the ten hour air temperature. One
value, in kw, for each ITD.
The turbine output in kw corresponding to the maximum allowable
back pressure and the appropriate load condition (full throttle, 3/4
load and 1/2 load).
Three values of turbine output, the first value being the output at
the actual operating back pressure at full throttle and the remain-
ing two values corresponding to the output at the maximum allow-
able back pressure for 3/4 load and 1/2 load operation respectively.
Total annual cost of the cooling system and total plant fuel cost.
One value for each ITD (mills/kwh).
Variable used to retain the amount of energy actually produced by
the turbine for a given set of operating conditions.
The number of capital cost factors (fixed-charge rates) to be used
for each site.
The number of base fuel costs and the number of peaking fuel costs
to be used for each site.
Net heat rate of turbine. Not used in program calculations but
used to skip over net heat rate values on data file (Btu/kwh) in
order to leave data file general.
Variable subscript used to denote the 32 values of temperature dur-
ation and air temperatures.
Subscript of the first non-zero temperature duration in the array
TDUR.
89
-------
NTU
OAMCT
ONM
OPBP
OPGHR
OPHRS
OPMAT
OPTDCT
OPTOUT
P
PAUSE
PCCST
PCT
PCTLD
PCTLOD
PCTTIM
PEAK
Subscript of the last non-zero temperature duration in the array
TDUR
Operating and maintenance cost as a percentage of capital cost.
Variable operation and maintenance cost (mills/kwh).
Actual operating back pressure of the turbine for a given set of
conditions.
Actual operating station heat rate of the turbine for a given
conditions.
set of
The number of hours per year of plant operation.
Operation and maintenance cost.
Optimization of Dry Cooling Towers. The name of the program.
The name of the file in which the program results are stored.
An array which contains the saturation pressures for a table of
saturation temperatures (°F) vs. saturation pressures (inches Hg).
Library subroutine used if the program does not function correctly.
Peaking unit capital cost ($/kw).
Percent of full auxiliary power requirement for a mechanical-draft
tower for a given ITD, air temperature and load.
Percent load at which the turbine is operating; i.e., 100% for full
throttle, 75% for 3/4 load, and 50% for 1/2 load.
The three operating conditions of the turbine; 100% for full throttle,
75% for 3/4 load and 50% for 1/2 load.
The percent of time the turbine will be operating at a specific load
level.
Peaking period for each site (summer or winter).
90
-------
PENCAP - Annual capital cost of replacing the generating capacity lost at the
air temperature that is equalled or exceeded ten hours per year.
PENLTY - The annual cost of replacing capacity and energy losses.
PENLTY = PENCAP + PENOPR.
PENOPR - Annual operating cost related to the replacement of the energy lost
due to operation at back pressures in excess of 3-1/2 inches Hg.
PFCST - Peaking unit fuel cost (£/106 Btu).
PLITD - Part load ITD.
SATP - Saturation pressure (inches Hg).
SATT - Saturation temperature (°F).
SITE - Variable used to store the site names.
SITENO - An array that contains the names of the 27 data files containing
site information. SITE01, SITE02, , SITE 27.
T - An array which contains the saturation temperatures for a table of
saturation temperatures (°F) vs. saturation pressures (inches Hg).
TANCST - Total annual cost of the cooling system and total plant fuel ($).
TANT - Variable used in finding the optimum ITD.
TAXKWH - Auxiliary energy in mwh .
TBLELV - Table of elevations (0, 3,000, and 6,000 feet above sea level).
TBLITD - Table of ITD from 30°F to 80°F in 10° increments.
TBLUQ - Table look-up subroutine used in the program.
TDUR - Temperature durations in the site data. These durations are the
percent of time during a year that the ambient air temperature is
within a given 5° range.
TEMP - Temporary variable used for storage in several locations in the
program.
91
-------
TKW
TKWGAN -
TMWH
TOTBTU
TOTCAP -
TOWER
TURBIN
The maximum output of the turbine under a given set of conditions
within the limiting values.
Excess energy, in mwh, produced when operating at full throttle and
below 3-1/2 inches Hg back pressure.
Nominal annual energy, in mwh, produced by the 800-mw generat-
ing plant for the given loading assumptions.
The total annual fuel consumption, in Btu, for the given loading
assumptions.
Capital cost corresponding to site elevation and ITD.
Variable containing the type of tower used in the cooling system,
either "Natural Draft" or "Mechanical Draft".
Name of the data file containing the turbine data.
92
-------
Glossary of Terms Used in Subroutine
"DATCHK"
DATCHK - Data check, name of the subroutine.
FILE - Variable containing the names of the data files.
INCT - Variable used as a counter.
NFIL - The number of the tape the data file is on.
SITE - Name of the file containing the site data.
SPVSST - Variable used to hold space for a file name that could be inserted
at a later date.
TURBIN - Name of the file containing the turbine data.
93
-------
Glossary of Terms Used in Subroutine
"TBLUQ"
A - Coefficient of X2 in the equation TBLUQ = AX2 + BX + C.
ABS - Library function used to find the absolute value of a number or
function.
B - Coefficient of X in the equation TBLUQ = AX2 + BX -f C .
C - Constant in the equation TBLUQ = AX2 + BX + C .
DET - Variable used in the calculation of A, B, and C.
I - Subscript.
N - Subscript.
NVAL - Number of X and Y values in the table which the value is being
sought.
TBLUQ - Name of the subroutine and the value which is being sought.
X - The variable for a corresponding Y value.
XI ,X2,X3 - The three X values used in the interpolation from the table of
values being used.
X1S,X2S, - The XI, X2, and X3 values squared, respectively.
X3S
XT - The table of X values being used.
Yl ,Y2,Y3 - The three Y values used in the interpolation from the table of
values being used.
YT - The table of Y values being used.
94
-------
ADRA
ADRFTH
ADRT
AEXIT
AFACTR
AIRDLT
AIRHT
AMAIRT
CNDCST -
COILS
CONDF
CPBTU
CURVGO -
DELTAP
DNDNSR -
DPTBL
DRAFTH
DRAFTLS -
DROOF
ELEV
Glossary of Terms Used in the
Tower Optimization Program
"TOWSIZ"
Air density ratio - elevation.
Adjusted draft height.
Air density ratio - temperature.
Air exit velocity.
Altitude density ratio table.
Air temperature increase.
Air heat gain.
Ambient air temperature.
Table of condenser costs.
Cost of heat exchangers.
Condensate flow.
Condenser cost per billion Btu per hour.
Table of water flows.
Coil losses.
Table of condenser cost.
Table of coil losses.
Draft height
Draft losses.
Cost of delta roof.
Elevation.
95
-------
EX1TLS
EXITV
FAIRTN1
GO (G0)
HEIGHT
HTREJ
ITD
LO (L0)
LOTBL
LOWDIA
N
P
QO (Q0)
QI (QT)
Q1TBL
RANGE
RTEMP
STACK
STDIAM
STKCST
STKHT
TEMPS
Exit losses.
Exit velocity (20 fps).
Final air temperature.
Water flow (Ibs./hr.).
Table of altitudes.
Heat rejection.
Initial temperature difference.
Air flow (metric tons/hr.).
Table of values of air flow.
Bottom diameter of tower.
Number of heat exchanger columns.
Pressure drop across the coil.
Heat rejection per hour per column (Btu).
Heat rejection per hour per degree of ITD.
Table of values of heat rejection per hour per column per degree
ITD.
Temperature change of condensate .
Temperature range table for condenser cost.
Cost of stack.
Table of tower stack diameters.
Table of stack costs.
Table of stack heights.
Temperature table.
96
-------
TFACTR - Temperature density ratio table,
TOTALH - Total draft height.
TOTLOS - Total losses.
UPAREA - Tower area at top.
OPDIAM - Top diameter of tower.
97
-------
Glossary of Terms Used in Subroutine
"TWRPIP"
BDIAFT
BPVALD
BPVALS
BPVCST
BTUPHR
DELHI
DELPSE
DGPM
DPCST
DPD
DPFT
DPLEN
DPRFT
DPRVAL
DVCST
DVD
EDGPM
EDPCST
EDPLEN
E DVCST
Diameter of the base of the cooling tower.
Cost per valve for bypass valves.
Bypass valve size (inches).
Bypass valve cost ($).
Btu per hour.
Delta Height (ft.).
Number of deltas per sector.
Drain flow rate (gpm)
Drain pipe cost ($).
Unit cost of drain pipe ($/ft.).
Table of pipe cost ($/ft.).
Length of drain pipe (ft.).
Pipe cost corresponding to pipe size ($/ft.).
Valve cost corresponding to valve size ($/valve).
Drain valve cost ($).
Unit cost of drain valves ($/valve).
Rate at which water must be drained in case of an
emergency (gpm).
Cost of emergency drain pipe ($).
Emergency drain pipe length (ft.).
Emergency drain valve cost ($).
98
-------
EDVD - Cost per valve for emergency drain valves ($).
FGPM - Flow rate at which cooling system is filled (gpm).
FPCST - FiHer pipe cost ($).
FPD - Cost per foot for filler pipe ($/ft.).
FPLEN - Length of filler pipe (ft.).
FVCST - Filler valve cost ($).
GAPDEL - Gallons of water per delta.
GPM - Flow rate for cooling condensate (gpm).
GPMPL - Flow rate per level (gpm).
HGPM - Flow rate in header pipes (gpm).
HPCST - Header pipe cost ($).
HPD - Cost per foot of header pipe ($/ft.).
HPLEN - Length of header pipe (ft.).
HVD - Cost per valve for header valves ($).
HVCST - Cost of header valves ($).
IFPSET - Number of filler pumps required.
IPSETS - Number of circulating pumps required.
MPCST - Cost of main supply line pipe ($).
MPP - Cost per foot for main supply line ($/ft.).
MPLEN - Length of supply line (ft.).
MPS - Size of main supply line (inches).
MVCST - Cost of main supply line valves ($).
MVD - Cost per valve for main supply line valves ($).
99
-------
NBPV - Number of bypass valves.
NDELPL - Number of deltas per level.
NDV - Number of drain valves.
NEDL - Number of emergency drain lines.
NEDVAL - Number of emergency drain valves.
NFV - Number of filler line valves.
NHSL - Number of header supply lines.
NHV - Number of header valves.
NMSL - Number of main supply lines.
NMVAL - Number of main supply line valves.
NODEL - Number of cooling deltas.
NOL - Number of lines.
NOLEV - Number of levels of cooling deltas.
NSECPL - Number of sectors per level.
NSL - Number of supply lines.
PCOST - Circulating water pump cost ($).
PDIA - Diameter of the different pipe (inches).
PI - 3.14159.
PIPSIZ - Name of subroutine that calculates diameter.
PSIZE - Table of pipe sizes (inches).
STCST - Storage tank cost ($).
STLBS - Weight of storage tank (Ibs.).
TNMSL - Total number of main supply lines.
100
-------
TPCST - Total pipe cost ($).
TPVCST - Total pipe and valve cost ($).
TVCST - Total valve cost ($).
101
------- |