WATER POLLUTION CONTROL RESEARCH SERIES • ORD-17O9OFTAO7/69
MATHEMATICAL MODEL OF
THE ELECTRODIALYSIS PROCESS
U.S. DEPARTMENT OF THE INTERIOR • FEDERAL WATER QUALITY ADMINISTRATION
-------
WAITER POLLUTION CONTROL RESEARCH SERIES
She Water Pollution Control Research Reports describe the results
and progress in the control and abatement of pollution in our
Nation's waters, ihey provide a central source of information on
the research, development, and demonstration activities in the
Federal Water Quality Administration, in the U. S. Department of
the Interior, through inhouse research and grants and contracts
with Federal, State, and local agencies, research institutions,
and industrial organizations.
A triplicate abstract card sheet is included in the report to
facilitate information retrieval. Space is provided on the card
for the user's accession number and for additional uniterms.
Inquiries pertaining to Water Pollution Control Research Reports
should be directed to the Head, Project Reports System, Planning
and Resources Office, Office of Itesearch and Development,
Department of the Interior, Federal Water Quality Administration,
Boom 1108, Washington, D. C. 20242.
-------
FWQ& Review Notice
Ihis report has been reviewed by the Federal Water
Quality Administration and approved for publication.
Approval does not signify that the contents neces-
sarily reflect the views and policies of the Federal
Water Quality Administration, nor does mention of
trade names or commercial products constitute
endorsement or recommendation for use.
-------
MODEL OF THE
ELECTRODIALYSIS PROCESS
by
Kenneth T. Pruyn
Joseph J. Harrington
J. Douglas Smith
Process Research, Incorporated
Cambridge, Massachsettts 02142
for the
FEDERAL WATER QUALITY ADMINISTRATION
DEPARTMENT OF THE INTERIOR
Program #17090 FTA
Contract #14-12-410
FW3A Project Officer, J. F. Roesler
Advanced Waste Treatment Research Laboratory
Cincinnati/ Ohio
July, 1969
For sale by the Superintendent of Documents, U.S. Government Printing Office
Washington, D.C. 20402 - Price 70cents
-------
Table of Contents
Page No.
Table of Contents ................................. ii
List of Figures ................................... iii
List of Tables ................ . ................... iv
Electrodialysis Process Description ............... 1
Electrodialysis Equipment ......................... 7
Mathematical Model ................................ 11
Electrodialysis Performance Equations ............. 16
Electrodialysis Process Cost Functions ............ 19
Capital Costs ..................................... 20
Operating Costs ................................... 27
Flow Sheet and Sample Outputs ........... . ......... 30
Recommendations and Suggestions ................... 39
Pertinent References ...................... , ....... 40
Appendix A - Derivation of Electrodialysis
Design Equations ........ 42
Appendix B
Program Description ............................. 47
Program Organization ............................ 48
Nomenclature .................................... 49
Listings ........................................ 55
-------
List of Figures
Figure No. Title Page No.
1 Basic Ion and Water Flow in Electro-
dialysis Stack 6
2 Expanded View of Electrodialysis Stack .. 8
3 Electrodialysis Process Flow Sheet 12
- iii -
-------
List of Tables
Table No. Title Page No.
1 Computer Output, 1 MGD Waste Treat-
ment Plant 35
2 Computer Output, 10 MGD Waste Treat-
ment Plant 36
3 Computer Output, 100 MGD Waste Treat-
ment Plant 37
- iv -
-------
Electrodialysis Process Description
The electrodialysis process removes ionized salts and
minerals from water by the migration of ions through plastic
ion-transfer membranes, under the influence of a D. C.
electrical field. Commercially available electrodialysis
membranes are thin sheets (4-40 mils thick) of cross-linked
organic polymeric materials. These membranes are similar
in composition to granular ion exchange resins, except for
their physical form. They commonly contain 40-50%, by
weight, of water which is distributed in extremely fine
pores honeycombing the plastic structure. The size of these
pores has been estimated at 10-100 A°. For the salt concen-
trations common to advanced waste treatment applications,
water transfer through these membranes will be virtually
zero. Two types of membranes are necessary to the process:
1) Cation permeable membranes, which will permit the
passage of cations but, because of the Donnan
exclusion principle, will prevent the passage of
anions; and
2) Anion permeable membranes, which for the same reason
will allow the passage of anions but not cations.
-------
If two metal electrodes are placed in an aqueous salt
solution and electrical current is passed between them, all
of the negative ions will move toward the positive electrode
and all of the positive ions will move toward the negative
electrode. If, now, alternative cation and anion permeable
membranes are interposed between the electrodes, the move-
ment of the ions in the electrical field will be seriously
restricted. Half of the ions will be able to move through
one type of membrane barrier and half will not; no ion
will be able to move through more than one barrier. This
is illustrated in Figure 1. In this figure, cations from
compartments whose left boundary is a cation membrane
(compartments 2, 4, 6, 8, 10, etc.) will move one compart-
ment to the left but must then stop because the second
membrane to be encountered will be an anion membrane, which
is impermeable to cations. Conversely, anions from compart-
ments whose right boundary is an anion membrane (compartments
2, 4, 6, 8, 10, etc.) will move one compartment to the
right, but must then stop because the second membrane they
encounter will be a cation membrane and impermeable to anions,
The net result, therefore, of the interposition of the
alternating cation and anion permeable membranes between
the electrodes is the creation of two types of compartments:
- 2 -
-------
diluting compartments (i.e., 2, 4, 6, 8, 10) which will lose
ions; and concentrating compartments (i.e., 3, 5, 7, 9)
which will gain ions.
The extent of salt removal, or demineralization,
obtained in the diluting compartments will depend upon the
amount of electrical current passed, the concentration of
ions present in the feed solution, and the flow rate of the
feed solution through the equipment. All ions are not, how-
ever, removed at identical rates, some moving preferentially
faster than others. This preferential, or selective, removal
of certain ions is described here by a Separation Factor. The
Separation Factor for each ion is the ratio of its percentage
reduction to the total concentration reduction through the
process.
In compartments 1 and 11, electrode reactions take place.
Hydrogen gas and hydroxide ions (OH") are produced at the
cathode, or negative, electrode; oxygen or chlorine gas and
hydrogen ions (H+) are produced at the anode, or positive
electrode. The effluent from these electrode compartments
may, in some situations, be useful for chemical operations,
being sources, respectively, of base and acid. The gases
produced at the electrodes may be used for disinfection, or
as a source of fuel. For purposes of this model, the cathode
and anode waste streams will be mixed with the concentrated
brine and be considered a part of the total waste stream.
- 3 -
-------
The mathematical model developed here for the electro-
dialysis process represents an attempt to describe the
present, moderately well documented state of the art of
commercially available equipment and operating practice.
Current state of the art has been established from site
visits and discussions with operating personnel at various
electrodialysis installations in this country, i.e., Lebanon,
Ohio; Pomona, California; Port Mansfield, Texas; etc.; from
the literature and from our own experience.
The modeling of presently available equipment in no way
implies that such equipment is optimally designed or operated
for waste water treatment applications. Our own independent
studies, in fact, indicate rather the opposite. Our design
optimization programs, for instance, suggest that modifications
of flow path geometry, i.e., thickness and/or width, could
result in greatly reduced membrane area and D.C. power
requirements. Our membrane fouling studies have indicated
modes of operation, or pretreatments, which would further
reduce membrane area and maintenance requirements.
The potential cost reductions, on the order of one-half
of present cost estimates, to be gained through these various
design and operating modifications are sufficient to justify
- 4 -
-------
their further examination and documentation on something
larger than a laboratory scale.
Such modifications have not been included in the model
developed here, because inadequate field experience has
been obtained to warrant their inclusion as present "state
of the art". An attempt, however, has been made to express
this model in sufficiently general terms as to readily
facilitate their later inclusion.
- 5 -
-------
Saline
Hater Saline Hater
to to cl
"-jj Ion Gaining Ion losing 2 ,
f; 2 Compartments Compartments or Oj
Cathode
*
- Feed
' 1
< 1 t
®
^Y
®
^^"^
0
i
^
J®
©x
©
i
i c
,@
* *
®k
^*. »
(s)
1
t
jjg,
[®J
(?)
i
1
r
' r
• i
X3>
* +
£L
"^ <
©
1
A | t
j@ l@
V ^
(•X, @\
"* 4
i
i
i
i
1
^®
®k
©
t
i
i
t <
l@
* *
\®\
^* t
\*«0
Anode
Feed
t
Jg)
©J
©
1
1
i
'
\
®
(ci^~*
(?)
» \~^
®
Cathode Concentrated Brine Product Anode
Waste or Haste Hater Waste
(7) Compartment Number (see text)
(N») Any Cation (Positive Ion) Like Sodium
(tr) Any Anion (Negative Ion) lake Chloride
C Cation Membrane
A Anion Membrane
FIGURE 1
BASIC ION AND HATER PLOW IN ELECTHODIALYSIS STACK
-------
Electrodialysis Equipment
The physical structure in which these demineralization
operations are carried out is called a membrane "stack".
Such a membrane stack is illustrated, in an expanded view,
in Figure 2. The stack consists of alternative anion and
cation permeable membranes, spacers, electrodes, end blocks
for making inlet and outlet manifold connections, heavy
steel end plates, and tie rods to provide uniform sealing
pressure.
The electrodes are thin metal sheets, stainless steel
or Hastelloy C being typical materials for the cathode.
Graphite, or platinum coated sheets of tantalum or titanium,
are typical materials for the anode.
The spacers are probably the most important component
in the electrodialysis stack, after the membranes them-
selves . They:
1) serve as gaskets, containing the water between the
membranes;
2) separate the thin, flexible membranes a constant
distance from each other; and
3) direct the flow of water across the surface of the
membranes.
- 7 -
-------
AJOUMULfll
00
I
Figure 2
Expanded view of Electrodialysis Stack
-------
A good spacer will be designed to provide an even distribution
of flow, and to promote sufficient turbulence at the membrane-
liquid interface to avoid the formation of polarization films
(to be discussed later). It has been suggested that by far
the greatest immediate process improvement, or cost reduction,
is to be found through optimal design of the geometry of
these spacers (see OSW R&D Report No. 325).
There are two general types of spacer designs used in
present commercially available equipment: the so-called
"sheet flow" design, as used in the Asahi Chemical electro-
dialysis plant at Webster, S. D.; and the so-called "tortuous-
path" design, as used in the Ionics, Inc., plant at Buckeye,
Arizona. The spacers shown in Figure 2 are of the tortuous-
path variety, the illustration being that of an earlv Ionics,
Inc. stack design.
Spacers are fabricated from sheets of electrical insulating,
plastic materials having good gasketing properties. Common
materials are plasticized polyvinyl chloride, polyethylene
or rubber-styrene compounds. In the tortuous-path type
spacer illustrated here, a flow path for the salt solution
is die-cut into one sheet of spacer material, along with
manifold holes and appropriate connections between manifold
holes and flow path. In the sheet-flow type spacer, a
-------
perimeter gasket containing the manifold holes and con-
nections to the flow path are die cut from sheets of gasket
material. Into the center of these perimeter gaskets, then,
are inserted plastic "screens", or mats, which provide the
necessary membrane support while still permitting hydraulic
flow through the compartments. Common insert support
materials have been expanded, or corrugated and perforated,
PVC or polypropylene.
Present commercially available electrodialysis stacks
contain approximately 500 membranes, or 250 anion-cation
membrane pairs. These stacks have hydraulic capacities of
the order of 250,000 gpd and conventionally provide a salt
reduction, or demineralization, or 40-60%. The concentrate,
or brine, stream is generally recirculated through the stack
at the same flow rate as the feed stream so as to minimize
pressure differentials across the somewhat fragile membranes.
- 10 -
-------
Mathematical^ Model
An electrodialysis process flow schematic, upon which
the mathematical model is based, is shown in Figure 3.
The development of the model follows the outline suggested
in Publication No. WP-20-9, "Preliminary Design and Simulation
of Conventional Waste Water Renovation Systems Using the
Digital Computer".
By writing material balances around the control volume
"A" in Figure 3, the following eguations can be derived
for flow and salt concentrations:
Q1 + Q9 = Q3 •»• Q? (1)
and:
Q1C1 * Q9C9 = Q3C3 + Q7C7
Assuming no water transport across the membranes, material
balances written around control volumes "B" and "C" will
establish'the following relationships:
Q2 * Q3 (3)
and:
Q5 « Q6. (4)
- 11 -
-------
FIGURE 3
Electrodialysis Process Flow Sheet
Fr-nri m
4«
r~
2
• --
*•
i
q
<
B,
1
1
^ np .... ,.,_ *
i
Uf^ J
T Cj
I ^
i
i
1
R). ' .8
V \
1
1
i
,»
ACID
A
i
1
^aste
DC - Dilute Compartments
CC - Concentrate Compartments
CR - Concentrate Recycle
ACID - Acid Feed
- 12 -
-------
If FRACT equals the fraction of total influent salt transferred
from diluting to concentrating compartments:
C3 = C2 d-FHACT)
Q2
C6 * C5 "*"
-------
and:
C6 = C? = Cg (12)
The above equations, then, can be combined to yield:
Q. Q-FRACT QQ
C SB C (—— + —- ) + C —— (13)
In this model, the stream vectors include stream identi-
fication numbers, volume flow rates (MGD), and concentrations
of the various contaminants pertinent to the electrodialysis
process. In the latter regard, the classification scheme
for contaminants is somewhat different from those generally
pertinent to other waste stream processes, in that primary
emphasis is placed on the composition of inorganic, rather
than organic, materials. The concentrations of each ionic
species present are specified in the stream vectors. The
electrodialysis stream vector is:
Contents of Stream Vector for Electrodialysis
Row Fortran Variable
No. Descriptive Parameter Names
2 Volume Flow, MGD SMATX( 2,I),Q
18 Concentration of Dissolved Sodium SMATX(18,I),PNA
19 Concentration of Dissolved Potassium SMATX(19,I),PK
20 Concentration of Dissolved Calcium SMATX(20,I),PCA
21 Concentration of Dissolved Magnesium SMATX(21,I),PMG
22 Concentration of Dissolved Ammonium SMATX(22,I),PNH4
23 Concentration of Dissolved Chloride SMATX(23,I),PCL
24 Concentration of Dissolved Bicarbonate SMATX(24,I),PHCO3
25 Concentration of Dissolved Sulfate SMATX(25,I),PS04
26 Concentration of Dissolved Nitrate SMATX(26,I),PN03
27 Concentration of Dissolved Phosphate SMATX(27,I),PPO4
- 14 -
-------
The decision vector for the electrodialysis process
contains parameters describing:
1) spacer design geometry, i.e., flow path thickness;
2) membrane properties, i.e., resistivity, current
efficiency and separation factors (see WP-20-AWTR-18)
for each of the ions present in the feed stream; and
3) mode of operation of the stack, i.e., current
density, product stream concentration, waste/
product ratio, waste stream pH, etc.
The electrodialysis decision vector is:
Contents of Decision Vector for Electrodialysis
Row
No. Descriptive Parameter
1 Total Dissolved Inorganic Solids, Product
Stream, ppm
2 Membrane Spacer Thickness, cm
3 Ratio of Spacer Thickness to Mesh Spacing
4 Input Stream Temperature, Deg-P
5 Cation Membrane Resistivity, ohm-sq.cm
6 Anion Membrane Resistivity, ohm-sq.cm
7 Current Efficiency
8 Limiting Current Density Parameter
9 Desired Waste/Product Stream Volume Ratio
10 Maximum Allowable Waste/Product Stream
Volume Ratio
11 Max Allowable Hydrogen Ion Concentration
in Waste Stream
12 Cosine of Mesh Angle, Theta
13 Separation Factor, Sodium
14 Separation Factor, Potassium
15 Separation Factor, Calcium
16 Excess Capacity Factor
17 Separation Factor, Magnesium
18 Separation Factor, Ammonium
19 Separation Factor, Chloride
20 Separation Factor, Bicarbonate
21 Separation Factor, Sulfate
22 Separation Factor, Nitrate
23 Separation Factor, Phosphate
Fortran Variable
Names
DMATX(
DMATX(
DMATX(
DMATX(
DMATX(
DMATX(
DMATX(
DMATX(
DMATX(
1
2
3
4
5
6
7
8
9
,N)
,N)
,N)
,N)
,N)
,N)
,N)
,N)
,N)
r
t
t
t
t
i
t
t
,
PPMOU
TKNS
TDEL
T
RCM
RAM
E
CURK
GAMMA
DMATX(10,N),GMAX
DMATX(11,N),HMAX
DMATX(12,N),COTHE
DMATX(13,N),ANA
OMATX(14,N),AK
DMATX(15,N),ACA
DMATX(16,N),ECF
DMATX(17,N),AMG
DMATX(18,N),ANH4
DMATX(19,N),ACL
DMATX(20,N),AHCO3
DMATX(21,N),AS04
DMATX(22,N),AN03
DMATX(23,N),AP04
- 15 -
-------
Electrodialysis Performance Equations
The influent salt concentration at Station 1, and the
required effluent salt concentration at Station 3, are
supplied to the program as input. The waste stream salt
concentration at Station 7 is calculated as the maximum
concentration that can be achieved without precipitating
calcium (as carbonate or sulfate) , or as the concentration
resulting from the lowest desirable Q7 to CU ratio supplied
as input, whichever is the smaller.
The repetitive physical unit of an electrodialysis
stack is called a cell-pair and consists of a diluting
compartment, concentrating compartment, and two membranes.
2
The average areal resistivity (ohm-cm ) of the cell pair, as
derived in Appendix A, is:
TKNS 71
RPA = RM
PRACT*CLAMCI U-FRACT) (C7/C1
The areal resistivity (ohm-cm ) at the effluent end of
the cell-pair iss
- 16 -
-------
The total electrical resistance of one cell-pair (ohm)
is the average areal resistivity divided by the effective
transfer area per membrane:
__ RPA / , f- »
(16)
The D.C. current (amperes) required to remove a given
amount of salt is obtained from Faraday's Law:
96,500*(C,-C,)*Q.*43.9
1 - - E*NCP - <17>
where E is the current efficiency of the process and NCP
is the number of cell-pairs per stack.
The D.C. power (kw) requirement then is obtained from
Ohm's Law:
EOPWR - , (18)
where EP is the applied cell-pair voltage.
The number of electrodialysis stacks required is also
obtained from Ohm's Law:
NSTK = (19)
There, in general , is a limiting value assigned to EP,
that is a complex function of the hydrodynamic design and
- 17 -
-------
operation of the particular electrodialysis equipment being
used. The expression for the polarization parameter, XICD,
is derived in Appendix A; the appropriate design value for
EP is obtained from Ohm's Law:
XICD*C *RPO
EP im— <20>
The pumping head loss (ft H20) through the electro-
dialysis stacks is calculated:
STKHD = 30 + PGRAD*LENGTH (21)
where 30 ft is the contribution due to piping and manifold
losses, PGRAD is derived in Appendix A, and LENGTH is the
flow path length defined as the ratio of membrane area to
flow path width.
- 18 -
-------
Electrodiajffsis Process Cost Functions
The goal of this effort was to develop simple relation-
ships that give a good approximation to the cost of com-
mercially available equipment. Data on cost of actual
operating -electrodialysis plants were taken from OSW R&D
Progress Report No. 134, "An Engineering Evaluation of the
Electrodialysis Process Adapted for Computer Methods for
Water Desalination Plants", Mason-Rust, Webster, South
Dakota (1965). Relationships among the various cost
items, particularly relative magnitudes, which were useful
in making functional approximations were given by OSW
R&D Progress Report No. 325, "Hydraulic Design Optimization
of the Electrodialysis Process", by Process Research,
Incorporated (1967). Some values used were in accordance
with conventions recommended by the Office of Saline Water
and are so noted where used in this development. Capital
costs are assumed to be referred to the Engineering News -
Record Cost Index January 1965.
- 19 -
-------
Capital Costs*
It was generally convenient to divide electrodialysis
capital cost items into three related groups and develop
separate cost equations for each group. The groups are
stack related equipment, D.C. power related equipment, and
auxiliary equipment.
Stack related costs include:
Electrodialysis Stacks
Stack Piping
Distribution Piping
Pumps, Dilute Stream and Concentrate Recycle.
Pumps for delivering feed water from the source to the
plant, and pumps for carrying off product water and waste
water which are usually included in electrodialysis plant
cost calculations were not included here.
D.C. power costs include:
Rectifiers
Process Electrical.
Auxiliary costs include:
Acid Storage Tank
Acid Feed Pumps
Instruments
Miscellaneous Equipment.
*Publisher's Note: Estimates of capital and operating and
maintenance cost presented in this report are low compared
to other authoritative sources. For the most reliable
cost information equipment manufacturers should be consulted.
- 20 -
-------
The cost of a building to house the plant would be in-
cluded in the cost of an electrodialysis plant operating on
its own. As a part of a waste water treatment facility the
building space would be shared. The building area required
by the ED process itself is calculated and passed along
to the executive program so that it may be combined with all
other building area requirements. Engineering cost, con-
tractor's profit, contingencies and omissions, and cost of
land are not included because they are handled by the
executive system.
Electrodialysis stack materials costs were based on
the following:
Membrane Cost
Membrane Size
Membrane Pairs/Stack
Spacers
Anode
Cathode
End Plate/End Block
Stack Hardware
Manufacturers Assembly Cost
Stack Erection Labor
Erection Labor Cost
Payroll Extras
General & Administrative
$1.00-3.00/sg. ft.
5-10 sq/ ft.
250-300
$.60-1.00 sq. ft.
$26.00/sq. ft.
$4.43/sq. ft.
$114.60/sq. ft.
$72.20/sq. ft.
$4.00/sq. ft. Membrane Area
160 man-hrs.
$4.50/hr.
15%
30%
- 21 -
-------
The last two values, payroll extras and general adminis-
trative overhead follow OSW conventions.
Combining typical values from the foregoing table gives
an erected stack cost as:
(2 x 275 x (1.00 + 3.00 + 4.00) + 26.00 + 4.43 4 2 x
114.60 + 72.20) x 5.0 + 160 x 4.50 x 1.15 x 1.30 -
$24,700/stack
Stack piping * $750/stack
Distribution piping * $200/stack.
Cost of pumps depends on the capacity of the pumps.
Flow in both the dilute and concentrate stream are essen-
tially equal and are given by:
Ql - Q7
The following values enter into pump costs:
Minimum number of concentrate pumps
Minimum number of dilute stream pumps
Base cost, stainless steel pump
Base cost, cast iron pump
Power cost, stainless steel pump
Power cost, cast iron pump
Stack piping head loss
Pump size scale factor
(22)
2
2
$1800
$1400
$52.50/HP
$30.30/HP
30 ft.
9.1
- 22 -
-------
Cost of pumps is given by:
(2 x 1400 + 30.30 x STKHD x 03/5.7) x (2 x 9.1/03)°'2 +
(2 x 1800 + 52.50 x STKHD x Q./5.7) x (2 x 9.1/O.j)0'2 «
[6400 + 14.5 x Q3 x STKHD] x (18.2/Q3)0'2
where the exponent 0.2 is a scale factor for pump costs,
and STKHD is the head loss through the stack and its associ-
ated piping.
A factor of 1.4 is applied to each item other than
stacks to account for erection labor. A factor of 1.05 is
applied to account for miscellaneous equipment. Combining
all of the foregoing costs gives
CSTK - 1.05 X 24700 x XNSTK + 1.05 x 1.4 x (750 + 200)
XNSTK + 1.05 x 1.4 x [6400 + 14.5 x 03 x STKHD] x
(18.2/Q3)°-2 (23)
where XNSTK = the number of stacks in the required plant, and
CSTK = the total capital cost of stack related items.
Equation (23) reduces to
CSTK = 27330 x XNSTK + [16800 + 38.08 x Q3 x
STKHD] x Q3~°'2 (24)
- 23 -
-------
D.C. power related items include rectifiers and other
process electrical equipment. Rectifier cost is made up
of
Rectifier Base $2700
Rectifier Power Cost $63/KW
Number of Stacks/Rectifier 4
So that rectifier cost is:
(2700/4) x XNSTK + 63 x EDPWR
where EDPWR is electrodialysis power in kilowatts.
Process electrical cost is given as
$600/stack + $40/KW of D.C. power
Combining D.C. power items gives
CDC = 1.05 x 1.4 x [(2700/4 + 600) x XNSTK + (63 + 40)
x EDPWR] (25)
which reduces to
CDC • 1874.25 x XNSTK + 151.41 x EDPWR (26)
Auxiliary equipment includes acid feed pumps, storage
for acid and process instrumentation. Data given by Mason-
Rust for the cost of acid storage tanks as a function of
tank capacity lead to the following relationship:
- 24 -
-------
Tank Cost = .0013 x (TANK SIZE)'715
The size required to hold a 14 day supply is
14 x Qg x 106
Combining these two gives:
.0013 x (14 x 09 x 106)*715
or
167.3 x Q9*715
Acid feed pump costs are based on the following values:
Number of Pumps 2
Pump Base Cost $640
Pump Capacity $101.90/GAL/DAY
so that cost is given by
2 x 640 + 101.9 x 106 x 09
Mason-Rust suggests the following function for process
instrumentation:
9600 + 750 x XNSTK x (XNSTK/4)'*17
- 25 -
-------
where $9600 is the base cost and $750 is the cost per stack
for a small plant with 4 stacks and
(XNSTK/4)"'17
is a scale factor to adjust down the unit cost of stack
instrumentation in very large plants.
Total auxiliary equipment cost, CAUX, is given by
combining the previous three costs
CAUX = 1.05 x 1.4 x [167.3 x Qg*715 + 1280 + 101.9
x 106 x Qg + 9600 + 750 x XNSTK x (XNSTK/4)~*17 (27)
which reduces to
CAUX = 15993.6 -I- 149.79 x 106 x Qg + 245.9 x Qg*715
+ 1102.5 x XNSTK x (XNSTK/4)"'17 (28)
The total electrodialysis plant capital cost, then, is
the sum of the stack, D.C. power and auxiliary equipment
related costs,
CCAP = CSTK + CDC + CAUX (29)
- 26 -
-------
Operating Costs
Operating costs for the electrodialysis process consist
of electrical power, membrane replacement, and other costs.
Pumping power is
STKHD x (Q- + Q., x CKWH x 3.15 x 36S/.827
J D
2780.5 x STKHD x Q3 x CKWH
where
and electrodialysis power is
EDPWR x CKWH x 24 X 36S/.9
so that total power cost, CPWR, is:
CPWR - (2780.5 x STKHD x Q3 + 9733.33 x EDPWR)x CKWH (30)
Membrane replacement and stack tear-down costs are based on
a membrane life of five years and stack disassembly frequency
of six times per year. In waste water treatment applications
operating conditions could be such that these figures would
be very much different than would actually be experienced.
Based on a labor rate of $3.50/hr. and a stack dis-
assembly time of 8 hrs./stack teardown cost is:
3.50 x 8 x 6 x 1.15 x 1.30 x XNSTK
- 27 -
-------
and membrane cost is
6.00 X 5 x 275 x .2 x XNSTK
so that membrane replace cost, CREPf is
CREP = 1900 x XNSTK (31)
Other costs include:
Chemicals (Sulfuric Acid)
Operating Labor
Maintenance Materials
Using an acid cost of $30/ton, the cost of acid is:
30 x 7619.7 x Qg x 365 - 83.44 x 106 x Q
Operating labor is
79«;
1.15 x 1.30 x 3.50 x (5.9 x Q3 ) * 365 = 11268 x Q3
where the exponent 0.725 is a scaling factor for large
sized plants.
Conforming approximately to OSW suggested costing pro-
cedures, maintenance materials and supplies are estimated
as 0.35% of plant cost, or in this case
.0035 x (CSTK + CDC + CAUX)
- 28 -
-------
Combining the preceding expressions gives total other
costs, COT, as:
COT = 83.44 x 106 x Og + 11268 x Q3'725 -I- .0035
x [CSTK + CDC + CAUX] (32)
Electrodialysis plant operating costs are the sum of
electrical power, membrane replacement, and other related
costs,
CTOP = CPWR + CREP + COT (33)
Building area, ABLDG, required to include stacks,
pumps and piping, instruments, membrane and parts storage
and chemicals is based upon the following:
Chemicals
37 + (Q3 x ID6)'33
Membrane and parts storage
55 + (Q3 x 106)'4
Instruments
(Q3 x 106)'33
Pumps and- Piping
60 + 20.36 x Q3
and Stacks
127.5 x XNSTK
Combining the above gives
ABLDG • 152 4- 127.5 x XNSTK + 20.36 x Q3 + 2 x
(Q3 x 106)'33 + (Q3 x lO6)'4 (34)
- 29 -
-------
Flow Sheet and Sample Outputs
A flow sheet for the computer program ELECT is
presented in the following figures, 4 through 6.
Sample outputs from the stand-alone version of the
program are given in Tables 1, 2, and 3, for typical
secondary effluent compositions and plant capacities of 1,
10 and 100 MGD, respectively.
- 30 -
-------
FLOW DIAGRAM: SUBROUTINE ELECT
Initialize
SMATX,
IERR
Pick up
Input
Values
Convert
Units, Find
Sums, Etc.
( Electroneutrality
I Satisfied?
>feo
Yes
T
Set:
IBRR(l) =
Negative Separation
Factors Proper?
Yes
Set:
IERR (2)
= 1
-©
c
Positive Separation
Factors Proper?
I Yes
Set:
IERR (3) •
= 1
- 31 -
-------
Find
Solubility
Parameters
c
Sulfate Concen-
tration too High?
res
Set:
IERR(4) •
= 1
No
/HC03 Concentration
V Essentially = 0?
Yes
No
FG - FG at
Desired
Gamma and
Acid Feed=0
Find FG
at HMAX
C
Square Root
Negative?
JNo
T
Set:
IERR(5) = 1
-©
C
Acid Feed Required
at HMAX Negative?
No
Set:
IERR(6)
= 1
- 32 -
-------
c
FG Less at HMAX
at Desired GAMMA?
No
Find FG at
Desired
GAMMA
r
Acid Feed < 0?
Find FG
at HMAX
-©
Find Acid
Feed at
Current FG
'Yes
Set Acid
Feed = 0
f Sulf ate Constraint
V Satisfied?
Yes
No
Reduce
Current
FG
- 33 -
-------
c
FG Too Small?
Yes
Set:
IERR(7) = 1
Find ED
Plant
Design
Values
Find
Cost
Values
Put Output
Values in
SMATX
Print
Sub-
routine
Name
"ELECT"
Acid Feed = 0?
Yes
'No
- 34 -
-------
TABLE 1
Computer Output, I MGD Waste Treatment Plant
CONCENTRATION SEPARATION
IN
130.
.. J 5.
60.
25.
20.
165.
300.
100.
10.
25.
850.
OUT
86.
8*
28.
15.
10.
77.
20" 1 .""
53.
5.
1.7.
500.
FACTOR
.79 SODIUM
1 .14 POTASSIUM
1 .25 CALCIUM
.98 MAGNESIUM
1.15 AMMONIUM
_1 .25 CHLORIDE
.77 BICARBONATE
1.11 SULFATE
1 .22 NITRATE
.72 PHOSPHATE
*100 SPACER_ THJCKNESS'Cil
.100 RATIO OF THICKNESS TO MESH SPACING
70..000 T.a WATER TEMPERATURE*_pEG F.
40.000 RCM=CATION" MEMBRANE RESISTIVITY*OHM-SQ.CM.
40 ..000. _ RAM=ANLQN MJENB RANE RESISTIVIT Y * OHM-SQ•CM •
.900 E=CURRENT EFFICIENCY
.707 COSINE OF.MESH ANGLE.THETA
.083 GAMMA=DESIRED RATIO Q7/Q3
2.000 GMAX=MAXJMgM RATIO Q7/Q3
.001 MAX HYDROGEN ION CONCENTRATION*WASTE STREAM
82.000 LIMITING CURRENT DENSITY PARAMETER*CURK
.010 COST OF ELECTRICITY»$/KWHR
ED POUER-KW
STACK HEAD LOSS-FT.H20
NUMBER OF STACKS
_NUMB£R_ OF. SECJT IF IERS_ ..
BUILDING AREA-SQ.FT.
FLOW PATH WIDTH-IN.
STACK COSTS-S
DC E9PT COSTS .-$.
AUX E9PT COSTS-S
TOTAL CAPITAL-^
POWER COST-CENTS/KGAL
._WEJfflRAN£_ R EP UAC.EMENT.rCENTS /KGA L
OTHER OPERATING COST-CENTS/KGAL
TOTAL OPERATING COST-CENTS/KGAL
INPUT STREAM-MGD
PRODUCT STREAM-MGD
WASTE STREAM-MGD
ACID FEED-GAL/KGAL WASTE
29.9472
40*3505
3602.5841
4
I
if 10. 0871
6.0005
12783LS3481
12031 .3024
24756.3443
164619.9948
1 .0823
2.0822
3.9093
7.0738
1 .0000
.9231
.0769
.3778
*EX1T*
- 35 _
-------
TABLE 2
Computer Output, 10 MGD Waste Treatment Plant
CONCENTRATION
IN
130.
15.
60.
25.
20.
165.
300.
100.
10.
25.
850.
OUT
86.
8.
28.
15.
10.
77.
201 .
53.
5.
17.
500.
SEPARATION
FACTOR
.79
1.14
1 .25
.98
1 .15
1 .25
.77
1 .11
1 .22
.72
SODIUM
POTASSIUM
CALCIUM
MAGNESIUM
AMMONIUM
CHLORIDE
BICARBONATE
SULFATE
NITRATE
PHOSPHATE
-.100 SPACER THICKNESS*CM
.100 RATIO OF THICKNESS TO MESH SPACING
70.000 T=WATER TEMPERATURE* DEG F.
40.000 RCM=CATION MEMBRANE RESISTIVITY.OHM-SQ.CM.
.40_.000 RAM=AN10N MEMBRANE RESISTiyiTY*OHM-SQ .CM.
•900 E=CURRENT EFFICIENCY
.•707 COSINE OF MESH ANGLE.THETA
•083 GAMMA=DESIRED RATIO Q7/Q3
2.000 GMAX=MAXIMUM RATIO Q7/Q3
•001 MAX HYDROGEN ION CONCENTRATION.WASTE STREAM
82.003 LIMITING CURRENT DENSITY PARAMETER*CURK
•010 COST OF ELECTRICITY.S/KWHR
ED POWER-KW
STACK HEAD LOSS-FT.H20
_TQTAL .TRANSFER AREA -SQ.FT. ,,. „
NUMBER OF STACKS
._NIJMBER OF RECTIFIERS
BUILDING AREA-SQ.FT.
FLOW PATH WIDTH-IN.
STACK COSTS-S
DC EQPT COSTS-$ _ _
"A"UX" EQPT COSTS -$"
TOTAL CAPITAL-S
POWER COST-CENTS/KGAL
_MEMBRANE REPLACEMENT-CENTS/KGAL
OTHER OPERATING COST-CENTS/KGAL
jroTAL OPERATING COST-CENTS/KGAL
INPUT STREAM-MGD
_PRQD_UCT STREAM-MGD _
WASTE STREAM-MGD
ACID -FJIED-GAL/KGAL WASTE
299.4718
40.3505
3602.5.8406.
40
10
6448.7172
6.0005
1113064.7588
J 2031 3. 0242
89335.6271
1322713.4101
1 .0823
2.0822
2.4300
__ 5.5945
10.0000
9.2308
.7692
.3778
- 36 -
-------
TABLE 3
Computer Output, 100 MGD Waste Treatment Plant
CONCENTRATION
IN
130.
15.
60*
25.
20.
165.
300.
100.
10.
25.
850.
OUT
86.
8.
28'
15«
10.
7.7.
201 .
53.
5.
1 7.
500 <
SEPARATION
FACTOR
.79 SODIUM
.14 POTASSIUM
.25 CALCIUM
.98 MAGNESIUM
.15 AMMONIUM
.25 CHLORIDE
.77 BICARBONATE
.11 SULFATE
.22 NITRATE
.72 PHOSPHATE
.100 SPACER THICKNESS'CM
.100 RATIO OF THICKNESS TO MESH SPACING
70.000 T=WATER TEMPERATURE* DEG F.
40.000 RCM=CATION MEMBRANE RESISTIVITY»OHM-SQ-CM.
40.000 RAM=ANION MEMBRANE RESISTI VITY»OHM-SQ .'CM.
.900 E=CURRENT EFFICIENCY
.707 COSINE OF MESH ANGLE*THETA
.083 GAMMA-DESIRED RATIO Q7/Q3
2.000 GMAX=MAXIMUM RATIO 97/93
.001 MAX HYDROGEN ION CONCENTRATION*WASTE STREAM
JB2.000 .LIMITING C.UBRENJ DENSITY PARAM£TER*CURK
.010 COST OF ELECTRICITY*S/KWHR
ED POW5R-KW
STACK HEAD LOSS-FT.H20
TOTAL TRANSFER AREA-SQ.FT.
NUMBER OF STACKS
NUMBER OF RECTIFIERS
BUILDING AREA-SQ.FT.
FLOW PATH WIDTH-IN.
STACK "COSTS -$
DC EQPT COSTS-S
AUX EQPT COSTS-S
TOTAL CAPITAL-S
POWER COST-CENTS/KGAL
MEMBRANE REPLACEMENT-CENTS/KGAL
OTHER OPERATING COST-CENTS/KGAL
TOTAL OPERATING COST-CENTS/KGAL
INPUT STREAM-MGD
PRODUCT STREAM-MGD
WASTE STREAM-MGD
_M!D_FEED-GAL/KGAL WASTE
2994.7179
40.3505
360258.4065
393
98
54524.1120
6.1074
10804862.9849
1190010.4918
649901.3938
12644774.8704
1 .0823
2.0458
1 .6556
4.7837
100.0000
92.3077
7.6923
.3778
- 37 -
-------
RECOMMENDATIONS AND SUGGESTIONS
It is our recommendation that a program be promptly
initiated to establish the effects of varying degrees of
pretreatment and operating conditions upon the process
design, and for that matter to establish the optimal design
of electrodialysis equipment for waste water treatment
applications. Such a program should follow the outline
suggested in the Conclusions of our final report under
FWPCA Contract No. 14-12-91, "Membrane Fouling Mechanisms
in the Electrodialysis Process".
Specifically, this program should include but not
necessarily be limited to the following:
1. Quantification of the effects of high temperature
(50°C) operation on both fouling rates and membrane
area reductions;
2. Further examination of the effects of competetive
surface area, in the form of entrained microscopic
air bubbles, on membrane fouling rates;
3. Determination of the effects of length and frequency
of resting periods, both with and without the
addition of enzyme active agents, on stack fouling
rates;
4. Determination of operating voltages (above the
point of polarization) necessary to prevent fouling
and the required additional acid feed to control
scale problems.
- 38 -
-------
Finally, the results of this program should be used as the
basis for electrodialysis design and operating improvements
to be incorporated into the mathematical model presented
here.
- 39 -
-------
PERTINENT REFERENCES
Calvit, B. W. and J. J. Sloan, "Operation Experience of the
Webster, South Dakota, Electrodialysis Plant", Pro-
ceedings, 1st Intl. Symp. on Water Desalination,
Washington, D. C. (1965).
Cooke, B. A., "Some Phenomena Associated with Concentration
Polarization in Electrodialysis", 1st Intl. Symp. on
Water Desalination, Washington, D. C. (1965).
Cowan, D. A. and J. H. Brown, "Effect of Turbulence on
Limiting Current in Electrodialysis Cells", Ind. and
Engr. Chem. 51, 1445 (1959).
Helfferich, Friedrich, Ion Exchange, McGraw-Hill, New York
(1962).
Mandersloot, W. G. B., "The Limitations of Electrodialytic
Water Desalting with Permselective Membranes and the
Requirements for Electrodialytic Equipment", 1st
Intl. Symp. on Water Desalination, Washington, D. C.
(1965).
Mason, E. A. and T. A. Kirkham, "Design of Electrodialysis
Equipment", Chem. Engr. Progr. Symp. Series 5_5, No. 25,
71 (1959).
Mason-Rust, "An Engineering Evaluation of the Electrodialysi
Process Adapted for Computer Methods for Water Desali-
nation Plants", Office of Saline Water, R&D Report No.
134 (1965).
Mintz, M. S., "Electrodialysis—Principles of Process Design",
I&EC 55_, 18-28 (1963).
Ments, M. V., "Water Desalination by Electrodialysis", I&EC
5_2, 148-152 (1960).
Rosenberg, N. W. and C. E. Tirrell, "Limiting Currents
in Membrane Cells", Ind. and Engr. Chem. 49, 780 (1957).
Smith, J. D. and J. L. Eisenmann, "Limiting Current Densities
in Wastewater Electrodialysis", Proc. Fifth Ind. Water
and Waste Conf., 141-164 (1965).
Smith, J. D. and J. L. Eisenmann, "Electrodialysis in
Advanced Waste Treatment", U. S. Department of the
Interior, Federal Water Pollution and Control Adminis-
tration, Publication No. WP-20-AWTR-18 (1967).
- 40 -
-------
Smith, J. D. and J. J. Harrington, "Hydraulic Design
Optimization of the Electrodialysis Process", Office
of Saline Water, R&D Report No. 325 (1968).
Smith, J. D., "Membrane Fouling Mechanisms in the Electro-
dialysis Process", Final Report under FWPCA Contract
No. 14-12-91 (1968).
Smith, J. D. and J. J. Harrington, "An Effective Model for
Limiting Current Density Characteristics of Electro-
dialysis Equipment", presented at 133rd National Meeting
of the Electrochemical Society, Boston, Mass. (May 5-9,
1968).
Speigler, K. S., Ed., Principles of Desalination, Academic
Press, New York (1966).
Solt, G. S., "Influence of Membrane Phenomena on Electro-
dialysis Operation", 1st Intl. Symp. on Water Desali-
nation, Washington, D. C. (1965).
Wilson, J. R., Ed., Demineralization byElectrodialysis,
Butterworths, London (1960).'
- 41 -
-------
APPENDIX A
DERIVATION OP ELECTRODIALYSIS DESIGN EQUATIONS
The areal resistivity (ohm-cm2) of one cell-pair in an
electrodialysis stack is obtained from the summation of the
resistivities of the two membranes and a dilute and concen-
trating compartment, i.e.,
RP = Rm + Rd
The resistivity contribution of the membrane-solution inter-
face is included in the membrane resistance term in this
model. The resistivities of the two solution compartments
are expressed in terms of solution concentrations and compart-
ment thicknesses:
Rc ' (A-3)
where:
t = compartment thickness , cm
A = equivalent conductance, mho-L/cm-equiv.
C = solution concentration, equiv/L
and d, c, are subscripts referring to diluting and concen-
trating compartments, respectively.
- 42 -
-------
Solution concentrations, at any point along the flow
path (con-current flow) are described as a function of the
fractional salt removed at that point:
Cd = C^l-f) (A-4)
Cc - C5 + C^f (A-5)
As a general practice, the two compartmental thicknesses
are maintained equal and the flow rates through the two
compartments are the same:
t - t = TKNS (A-6)
a c
Q2 - Q5 (A-7)
The concentrating stream concentration then can be described
in terms of the total salt reduction obtained through the
process (PRACT) and the outlet concentrate stream concen-
tration (C7) :
C = C_ - C,*FRACT + C,f (A-8)
c / i i
The cell pair resistivity then is obtained by substituting
equations A-8, A-6, A-4, A-5 and A-2 into A-l:
JL - PRACT •«- f
Cl
(Note: The term CLAM is used in the program for A)
- 43 -
-------
At the effluent end of the process, f = FRACT and C.(l-FRACT)
C3. The cell-pair resistivity at this point becomes simply:
"PC - \ *
The averaged resistivity through the entire process is
obtained from the integration:
FRACT
df
(A-ii)
FRACT
£ df
Performing this integration results in the expression:
- Rm
(1-FRACT)
Generalized models for predicting the mass (salt)
transfer characteristics and pumping energy requirements
for electrodialysis systems have been developed by Process
Research, Incorporated, under previous contract with the
Office of Saline Water, U. S. Department of the Interior.
These models are described in OSW R&D Progress Report No. 325
Limiting mass transfer rates are described through a
so-called polarization parameter, (i/C,) , which in turn
d lira
is expressed as a function of fundamental solution and
membrane properties and standard hydrodynamic parameters:
- 44 -
-------
/ 0.6
*•»
2/3
where:
2
i = current density, ma/cm
C, =* dilute stream solution concentration, eq/L
F = Faraday's constant, 96,500 coulombs/equivalent
k' a» empirical constant
2
v = kinematic viscosity, cm /sec
t = distance (cm) between supporting strands of
membrane support material
t , = flow path thickness , cm
z", » electrochemical valence of ionic species j
t = transference number of species j in the membrane
t = transference number of species j in the solution
S3 2
D . =* Fickian diffusion coefficient of species j , cm /sec
e * void fraction of spacer flow path
e = angle between superficial and actual flow directions
V = superficial flow velocity, cm/sec
n * number of cell pairs
a = flow path width, cm
V t
In this model, the term (evcose> is the msiSS transfer
Reynolds Number for the electrodialysis process. The mem-
brane and solution properties and the various constant
coefficients can be lumped together into one term:
- 45 -
-------
CURKV = - — - (A- 14)
For present commerciallv available membranes and tvpical
waste water compositions, CURKV is about 1.3.
The polarization parameter then is described in terms
of hydrodynamic properties only:
V t °'6
1 » l'l ( - ?__) (A- 15)
;
(In the program, the term XICD is used for (i/C^) . Also,
lim
since P_/nat, is used for V , CURKV is modified to:
CURK » (103)°'6CURKV - 82).
Pumping energy requirements are obtained from the
head loss gradient through the electrodialysis process:
hr u2 2V
IT ' ^-^775 - f
(±^.) cose
where:
ht
r= * pressure drop across electrodialysis spacer, cm H20
n « wetted surface area per unit of flow path volume,
cm /cm .
2V
In this model, the terms (flvcOg9) and (jp) are the
hydraulic Reynolds Number and hydraulic diameter, respec-
hL
tively. In the program, the term PGRAD is used for sp.
- 46 -
-------
APPENDIX B
PROGRAM DESCRIPTION
The computer program is for modeling the electrodialysis
process as a unit process in the simulation of advanced
waste treatment plants. The program is not intended to
design a completely self-contained electrodialysis plant.
There are two versions of the program. One is a sub-
routine for use with the "Executive Digital Computer Program
for Preliminary Design of Wastewater Treatment Systems'*.
The second version is an independent program which is
completely self-contained. The two versions are essentially
equivalent except for the handling of input and output.
The subroutine version gets input values from the stream
matrix SMATX, the decision matrix, DMATX, and from other
variables in COMMON. Output values are placed in the
stream matrix and in variables in COMMON.
The self-contained version accepts input values from
an input file, which could be a punched card deck, and puts
output values in an output file, which could be a line
printer or other device as appropriate.
Listings of both versions are given below.
- 47 -
-------
Program Organi 2ation
The program is composed of six parts. Ksch part of
the program i headed by comment cards.
Part 1 picks up input parameters from the stream matrix,
SMATX, and from the decision matrix, DMATX.
Part 2 performs unit conversions such as from ppm to
moles/liter, finds summations of concentrations, and finds
concentrations of each species in the product water stream.
Input stream validity checks, described in detail in the
section on error conditions, are also made.
Part 3 finds allowable waste stream concentration and
required waste stream acid feed based on bicarbonate and
sulfate solubility limitations. The resulting product
stream volume is also found.
Part 4 computes the electrodialysis system design para-
meters required.
Part 5 computes the capital cost values and operating
cost values.
Part 6 puts output stream values in stream matrix.
- 48 -
-------
NOTATIONS
Fortran Name
Description
ABLDG
ACA
ACL
AH 003
AK
AMG
ANA
ANH4
AN03
APO4
AS04
CAC
CADI
CAUX or CCOST(3,N)
CDC or CCOST(2,N)
GDI
CLAM or -A-
CLDI
CPWR or COSTO(1,3)
CREP or COSTO(2,N)
O03K1
COTHE
Building area required, sq ft
Separation factor, calcium
Separation factor, chloride
Separation factor, bicarbonate
Separation factor, potassium
Separation factor, magnesium
Separation factor, sodium
Separation factor, ammonium
Separation factor, nitrate
Separation factor, phosphate
Separation factor, sulfate
Acid feed required, equivalent acid/1 waste
Calcium concentration, mole/1
Auxiliary capital costs, $
D.C. related capital costs, $
Input stream concentration, equivalent/1
Specific solution conductance, mhos/cm per
equivalent/1
Chloride concentration, mole/1
Power cost, $/yr
Membrane replacement, $/yr
Ratio XCAC03/(XK1*XK2) where
XCACO3 = calcium carbonate solubility product
XK2 = second bicarbonate protolysis constant
Cosine of mesh angle, theta
- 49 -
-------
NOTATIONS
Fortran Name
Description
COT or COSTO(3,N)
CSTK or CCOST(1,N)
CURK
EDPWR
E
PG
FGMAX
FGMIN
FRACT
GAMMA
GMAX
HOO
H003D
HMAX
IERR
Other operating cost, $/yr
Stack related capital costs, $
Limiting current density parameter
D.C. power required for electrodialysis
Current efficiency
FRACT/GAMMA (Variable)
FRACT/desired GAMMA
FRACT/maximum GAMMA
Fraction demineralized
Volume ratio of waste stream to product
stream
Maximum allowable waste/product stream
volume ratio
Hydrogen ion concentration, waste stream,
mole/1
Bicarbonate concentration, mole/1
Max allowable hydrogen ion concentration
in waste stream, mole/1
Error indicator array; values defined below
The following table defines the meaning of the various error indications.
If the element listed in the table below is set = 1, the error condition
corresponding to that item exists
ELEMENT
IBRR(l)
Meaning if IERR(I)=1
The difference betwen SMP and SMN exceeds 5%
of SMP. Electroneutrality is violated.
- 50 -
-------
NOTATIONS
Fortran Name
Description
IERR(2)
IERR(3)
IBRR(4)
IERR(5)
IERR(6)
IERR(7)
The difference between SMN and SMCAN is not
within 5% of SMN. Negative separation factors
are improper.
The difference between SMP and SMCAP is not
within 5% of SMP. Positive separation factors
are improper.
Calcium sulfate solubility product is exceeded
by input stream concentrations.
Negative square root was attempted in quadratic
solution for acid feed rate at limiting waste
hydrogen ion concentration.
Minimum brine stream pH (HMAX) is violated even
at zero acid feed.
Calcium sulfate solubility is exceeded even at
the minimum brine stream concentration (GMAX).
NSTK
NRECT
PCA
PCL
PHCO3
PK
PMG
PNA
PNH4
PNO3
PO4DI
POCA
Number of stacks
Number of rectifiers
Input stream concentration of calcium, ppm
Input stream concentration of chloride, ppm
Input stream concentration of bicarbonate, ppm
Input stream concentration of potassium, ppm
Input stream concentration of magnesium, ppm
Input stream concentration of sodium, ppm
Input stream concentration of ammonium, ppm
Input stream concentration of nitrate, ppm
Phosphate concentration, mole/1
Product stream concentration of calcium, ppm
- 51 -
-------
NOTATIONS
Fortran Name
Description
POCL
POHCO
POK
POMG
PONA
PONH4
PONO3
POPO4
POS04
PPND4
PP04
PSO4
Ql
Q3
Q7
Q9
RAM
ROM
SMCAN
Product stream concentration of chloride, ppm
Product stream concentration of bicarbonate,
ppm
Product stream concentration of potassium,
ppm
Product stream concentration of magnesium,
ppm
Product stream concentration of sodium, ppm
Product stream concentration of ammonia, ppm
Product stream concentration of nitrate, ppm
Product stream concentration of phosphate,
ppm
Product stream concentration of sulfate, ppm
Total dissolved inorganic solids, product
stream, ppm
Input stream concentration of phosphate, ppm
Input stream concentration of sulfate, ppm
Input stream flow, mgd
Product water flow, mgd
Waste stream flow, mgd
Acid feed rate, gal. acid/k gal. waste
Anion membrane resistivity ohm-sq cm
Cation membrane resistivity ohm-sq cm
Sum of the products of the equivalents of each
negative species with its separation factor
- 52 -
-------
NOTATIONS
Fortran Name
Description
SMCAP
SMN
SMP
SO4DI
STKHD
T
TDBL
TKNS
TTA
VISC or
WIDTH
WMCA
WMCL
WMHCO
WMK
WMMG
WMNA
WMNH4
WMNO3
WMP04
WMSO4
XCASO
Sum of the products of the equivalents of each
positive species with its separation factor
Sum of the equivalents of each negative species
Sum of the equivalents of each positive species
Sulfate concentration, mole/1
Stack headless, ft - HO (includes manifold
loss) 2
Input stream temperature, F
Ratio of spacer thickness to mesh spacing
Membrane spacing thickness, cm
Total membrane transfer area required, sq ft
Kinematic viscosity of water, sq cm/sec
Flow path width, in.
Molecular weight, calcium
Molecular weight, chloride
Molecular weight, bicarbonate
Molecular weight, potassium
Molecular weight, magnesium
Molecular weight, sodium
Molecular weight, ammonium ion
Molecular weight, nitrate
Molecular weight, phosphate
Molecular weight, sulfate
Calcium sulfate solubility product
- 53 -
-------
NOTATIONS
Fortran Name
Description
XFTNA
XK1
XKDI
XKW
XMGDI
XNH4D
XNADI
XN03D
Flow parameter, I/sec/cm
First bicarbonate protolysis constant
Potassium concentration, mole/1
Dissociation constant for water
Magnesium concentration, raole/1
Ammonium concentration, mole/1
Sodium concentration, mole/1
Nitrate concentration, mole/1
See CLAM
See VISC
- 54 -
-------
Listings
Following are listings of the source code for both
versions of the program and a listing of the data file used
for the sample cases.
- 55 -
-------
. SUBROUTINE ELECT
COMMON/EDERR/ XCASO, SMP,SMN,SMCAP*SMCAN» I ERR ( 7)
COMMON/EDJP/ ABLOG,NSTK,NRECT*EDPWR,STKHD*TTA,WIDTH,Q9
00 5 1=1,7
5 IERRU>=0
00 6 1 = 1 * NROWS
SMATX=SMATX(I,IS1 >
SMATX<2,OS1 >=0.
c
C PART 1- PICK UP INPUT PARAMETERS
C
Q1=3MATX<2,IS1 )*DMATXU6,N)
CKWH=OMATX(1
PPMOU=DMATX(
TKNS=DMATX( /N)
TDEL=OMATX(
T=DMATX( ,
RCM=DMATX<
RAM=DMATX(
E=OMATX< *N
CURK=DMATX<
GAMMA=DMATX<
GMAX=DMATX<
HMAX=DMATX<
COTHE=DMATXC
ANA=DMATXC
AK=DMATX<
ACA=DMATX(
AMG=DMATX(
ACL=DMATX(
AHC03=DMATX<
AS04=DMATX( »N)
AN03=OMATX< >N>
AP04=DMATX( »N)
PNA=SMATX( ,131)
PK=SMATX<
PCA=SMATX<
PMG=SMATX( ,131)
PNH4=SMATX<
PCL=SMATX(
PHC03=SMATX( ,131)
PS04=SMATX< *IS1>
PN03=SMATX( ,IS1 )
PP04=SMATX< ,IS1 )
C
C PART 2- PERFORM UNIT CONVERSIONS, PPM TO MOL/LTR, ETC.
C SUrt CONCENTRATION VALUES* ETC.
C
WMNA=22.9898E3
WMK=39.102E3
WMMG=24*312E3
- 56 -
-------
WMCL=35.453E3
WMNH4=18.03858E3
WMHCO =61 .0094E3
WMS04=96.0616E3
WMN03=62»0049E3
WMP04=84»949E3
XNADi=PNA/WMNA . . . .
XKDI=PK/WMK
CADI=PCA/WMCA .
XMGDI=PMG/WMMG
CLDI=PCL/WMCL
XNH4D=PNH4/WMNH4
HC03D =PHC03/WMHC.O. .............. _
S04DI=PSQ4/WMS04
XN03D -PN03/WMN03 ............................ . . . . ..
P04DI=PP04/WMP04
SMP=XNAni+XKQI+2.*CADI+2-*XMGPJ _ + XNH4D
SMN=CLOI+HC03D *2.*S040I +XN03D + P04DI * 1.5
_PPMT*PNA+P.K+PQA+PMG+PC,LtPliC03+PS04+PN03 + PP04 +PNH4
sijMCI=SMP-SMN+HCd30
SMCApaANA*XNApI+AJK*XKJ)I*ACA*CApI*2.+AMG*XMQOI*2.+ANH4*XNH4D
SMCAN=ACL*CLDI+AS04*S04DI*2.+
1 ANp3*XN03p+HC030*AHC03*AP04*P040I*l .5
GDI =.5*
RAC1*AK_)_
.-FRACT*A'MG>
.-FRACT*ANH4)
PQHCO =PHC03*U .-FRACT*AHCQ3>
POS04*PS04*U .-FR ACT'* AS 04)
___
POCA=PCA*<1 .-FRACTfAC'A) "
P QP 04rPPQ4*l L^jiFR AC T*AP 04 >
TPdUT»P6NA+PbK+POCA+P6MG+PONH4+POCL+POHCO+POS04*PON03+POP04
C_ _ _____
C ERROR CHECK INPUT "VALUES
C _______ __ .
IF(ABS(SMP-SMN)-.05*SMP)832*832*830
830 CONTINUE
IERRC1)»l
JB32 CONTINUE _
IF(ABS
-------
C
C
C
7110
WASTE STREAM ACID FEED BASED ON SULFATE AND
BICARBONATE CONCENTRATION.
IF7110*71 10*7120
XKW = .13E-14*(2.3>**«T-32.)/18.)
GO TO 7150
7120 IFC113.-T) 7130*7138*7143
7130 XKW = 4.0E-14*<4«0)**«T-1 13.)/45.)
GO TO 7150
7140 XKW = 1 .47£-14*<2.72>**«T-86.>/27.)
7150 CONTINUE
IFCT-70.) 7210*7210*7220
7210 XCASO * /l 36 • >**2
GO TO 7270
7220 IF(T-103.>7230*7230*7240
7230 XCASO = 2.4E-4
GO TO 7270
7240 IFCT-130.) 7250*7250*7260
7250 XCASO » (2 . 1 5* <0 • 7953488 ) **« T- 100 • ) /30 . )/l 36 . ) **2
GO TO 7270
7260 XCASO = (1 .71*<0.7427)**<(T-130.)/30.)/136.)**2
7270 CONTINUE
C03K1 =10.** «12343.5/
C
C
._ CL...
-_.417
. 619
505
C
__C
C
__51fl_
. 20 „
C
CHECK SULFATE CONSTRAINT TO SEE IF IT MAY BE BINDING
IF=t
GO TO 1000
CONTINUE ... ...
FGMIN=FRACT/GMAX
IF < CADI *HCQ3D-1 .£-30)505*505*510
FGsFGMAX
IACID=2
.CHECK, TO J5.EE IF.HMAX ..LIMITS. H C 03 _C_ONCE NT RATION
C2F2=ACA+AHC03
HCO=HMAX* FGAM FROM EQN 2
S = l .-(C3F2*HMAX*(XK1+HMAX)) _________ _ _______
S=C2F2**2-4.*C1F2*S
- 58 -
-------
IF (S) 110* 130* 130
110 IERR<5)=1
GO TO 1000
130 FGs(-C2F2+SQRT(S)>/C2.*ClF2>
CAC=SUMCI+FG*SMCIA -XK1*HC03D *(!•+
1 FG*AHC03>/CXK1+HMAX)+HMAX-XKW/HMAX
C
C IF CAC IS NEGATIVE'ERROR
C
IF 140.155*155
140 IERR(6)=1
GO TO 1000
155 IF(FG-FGMAX>520*520*515
515 FG=FGMAX
C
C FIND H TO SATISFY HC03 FOR CURRENT FG
C (THIS IS RE-ENTRY FROM 671.) FIND CAC REQUIRED.
C
520 B*XK1
C=-(FG*FG*ACA*AHC03+FG*+1•)*CADI*HC03D /C03K1
HCO»<-B+SQRT(B*B-4.*C»/2.
IFCHCO-1.£-35)640*640*630
.C . .
C IF HCO IS 'ESSENTIALLY' ZERO THIS CONSTRAINT CAN'T BIND
C
630 CONTINUE
CAC«SUMCI+FG*SMCIA -XK1*HC03D *(1.+FG*AHC03>/(XK1+HCO)
1 +HCO-XKW/HCO
IF(CAC>640*640*650
640 CAC*0.
IACID=2
650 CONTINUE
.C
C NOW HAVE FG*CAC DESIRED FOR HC03 CAND HMAX).
C MAY HAVE TO ADJUST DOWN TO SATISFY S04-
C (THIS IS FG RE-ENTRY FOR CAC=0.J
C _ .
XLM3FG*FG*ACA*AS04*CADI*S04DI+FG*«ACA+AS04>*CADI*S04DI
1 +CADI*ACA*CAC>+CADI*S04DI+CADI*CAC-XCASO
IF(XLM)470*470*660
C
C IF XLM IS NOT>0 S04 CONSTRAINT IS SATISFIED.
_C
C TRY REDUCING FG BUT NOT BELOW FGMIN
C
660 IF
Q3=Q1-Q7
- 59 -
-------
Q9=Q7*CAC/34.7
C THE CONSTANT ABOVE* 34.7* IS BASED ON 93ZH2S04
C
C
C PART 4- COMPUTE ELECTROOIALYSIS SYSTEM DESIGN VALUES
C
FT = Q3 * 43.8079
C XFTNA=FLOW PARAMETER- FT/NA> LTR/SEC/CM
XFTNA=2.412E-3
CLAM = 0.1*EXP(.01 l*(T-68.»
VISC = .01*EXPC-.011*CT-68.»
RM = *«1 .-.393*
1TDEL)*COTHE>**.6)
RN=XFTNA*1000/((1.+TDEL)*VISC*COTHE>
HDIAM=TKNS*<1 -3.14159265*TDEL/8)/( 1+TDEL)
PGRAD=VISC**2*C.023*RN+6«98E-5*RN**2>/
XIGD=CUHLI /'TKNS*XFTNA**.6
RPO=RM+TKNS/
XKPP=9.65E+07*(FRACT*RM+TKNS/(CLAM*CDI>*
1 ALOG <1 ./«! .-FRACT)*U .-FRACT*SOL) » )/C E*RPO*< 1 .-FRACT»
XKP=96.5*CDI**2*FRACT*U «-FRACT)/E
XNAP= XKPP/XICD+FT
EDPWR=XKP*RPO*XICD*FT/t000.
PMH20=PGRAD*XNAP/FT*XFTNA/100.
STKHD=3 .281*PMH20+30.
TTA = XNAP * 1.0764E-3
CURDE =96500.*CDI*FRACT*FT/E/TTA
ELL=TKNS/TDEL*.3937
XNA=FT/XFTNA
NSTK=TTA/917.+.9
XNSTK=NSTK
NRECT=XNSTK/4.+.5
IF(NRECT)13»13*14
13 NRECT = 1
14 CONTINUE
WIDTH=XNA*.3937/CXNSTK*275.)
C
C PART 5- COMPUTE COST VALUES
C
CCOST(1*N)=27330.*XNSTK+(16800.+38.08*Q3*(STKHD))*Q3**<-.2)
CCOSTC2/N)=1874.25*XNSTK+151.41*EDPWR
CCOST(3*N>=15993.6+1 49.79E6*Q9+245.9*G>9**.71 5
1 +1102.5*XNSTK*(XNSTK/4.)**(-•17)
COSTO(1*N)=(2780.5*STKHD*Q3+9733.3333*EDPWR)*CKWH
COSTO(2»N)=1900.*XNSTK
COSTO(3*N)=83.44E6*Q9+11268.*Q1**.725+.0035*CCCOST<1*N)+
.1 CCOST<2*N)+CCOST(3*N)>
A8LDG=152.+127.5*XNSTK+20.36*Q3+2.*(Q3*I.E6)**.33+(Q3*i.£6)**.4
GAMMA=FRACT/FG
C
C PART 6- PUT OUTPUT STREAM VALUES IN OUTPUT MATRIX*ETC.
C
- 60 -
-------
SMATX< ,OS1>=PONA
SMATXC *OS1>=POK
SMATXC *OS1)=POCA
SMATXC *OS1>=POMG
SMATX< »OSl>sPONH4
SMATX( *OSD=POCL
SMATXC *OS1>=POHC03
SMATXC *OS1>=POS04
SMATXC *OS1)=PON03
SMATXC *OS1)=POP04
SMATXC2»OSl)=Q3
SMATXC ,OS2)aCPNA-PONA)/GAMMA
SMATXC »OS2)=CPK-POK)/GAMMA
SMATXC »OS2)sCPGA-POCA)/GAMMA
SMATXC *OS2)=CPMG-POMG)/GAMMA
SMATXC ,OS2>=CPNH4-PONH4)/GAMMA
SMATXC *OS2)=CPCL-POCL)/GAMMA
SMATXC ,OS2>=CPHC03-POHC03)/GAMMA
SMATXC »OS2>=CPS04-POS04>/GAMMA+CAC*WMS04/2.
SMATXC *OS2)=CPN03-PON03)/GAMMA
SMATXC ,OS2)=CPP04-POP04)/GAMMA
SMATXC2,OS2)=Q7+Q9
Q9=CQ9/Q7)*1000.
1000 CONTINUE
199 FORMATC5X»5HELECT>
PRINT 199
RETURN
END
- 61 -
-------
L = b
IN = 2
READC IN/b>PNA*PK*FCA*PMG*PNH4,PCL*PHCG3»PS04>PN03*FP04
IF(PI\1A>4, 7, 4
READC iN>b >ANA* AK* ACA» AMG* ANH4* ACL* AHC03»Aii04» ANOS* AP04
"~ ........ IF(Q1)6,3*6
6 READC IN*b)GAI*iMA.» GMAX*HMAX*C01HE»PHM
b FORMATC 10F7. 1 )
999 FORMATC1H1)
WMNA=22.9898E3
tot"iCA-:40.08L3
IWfiMG=24. 312E.3
MCL=3b.4b3E3
WMNH4=18.038b8E3
WhHCO =61 .0094E3
WMS04=96.0616E3
WMN03=62.0049L3
WMP04=84.949E3
XNADI=PNA/WMNA
XKDI = PK/l'.;MK
CADI=PCA/^MCA
HC03D =PHC03/WMHCO
XN03D_=PNg3/>MNQ3.
P04Di=PPG4/WMP04
.. SMP_=XNADIt.XKD!_+2...*CADI+2.*Xtt.GPI . + XNH4D
SMN=CLDI+HC03D +2 . *S04DI +XN03D + F04DI * 1.5
PPM7SPNA + PK + PCA + PMG+PCL + PHCP3 + PS04 + PN03 + .P.P04... +PNH.4
SUMCI=SMP-SMN+HC03D
'
CDI = . b* ( 'SMP + SKN )
CIAI = ANA*PNA + AK*PK + ACA*PCA + . AM.G*PMG. +ANH4*PNH4 + .......
1 ACL*PGL + AHC03*PHC03 + AS04*PS04 + AN03*PN03+APO 4* PP04
SMCIA =SMCAP-S^CANj-HCjp3.P..*AHC.03 ...... .................. __________ ......... _„..._. ........ ._.,
"FRACT = VPPMT - PPMOU >/CIAI
FGMAX=FRACT/GAMMA ______________________ _______ ___________ ........ ...._ ..... .. __ ....... _..
PCJNA=PNA*"( 1. -FRACT*ANA )
POK=PK*( 1 .-FRAC1*AK) _ _ ...... . ...... .. ....... _______
"" POM~G=PMG* ( 1 ."-FRAC1 *AMG)
POCL = PCL*( 1 .-FRACT*ACL) ...... . ... __ ......... . _ .. ...... _______________________________________________
"" PNH4*c"i .-FRAC1*ANH4)
-. =PHCQ3*C 1 . -FRAC1*AHC03.) _____ ........ __________ ....... ______________________ ..... _________
POS04=PS04*( 1 . -FRACT*AS04)
= PCA*( 1 .-FH ACT* AC A)
POP04 = PP04*( 1 . -FKACT*AP04)
- 62 -
-------
C
C" EkhOK CHECK INPU'l VALUES
C
I F (ABS ( SNP -SfrN > - • 0 b* Stop )8 32 • B 32, 8 30
830 CONTINUE
WKI1L(L.»S31 >
831 FOKMATC44H ELECT KONEU1KAL IT Y IS NOT WITHIN 5% AS GIVEN) ....
832 CONTINUE
IF(Abi>(i3MN-SMCAN>-.05*StoN>83Sj83b>833
'833 iukiTE!SMN,SMCAN
834 FOKMAK34H SEPARATION FACTORS IMPROPER.SMN= E 1 2 . A, 7H .SMCA.N = E 1 2 . 4)
" 8 3S "I F ( ABS C sKiP - ShC AP ) - . 05* SMP ) 8 38 * 8 38 * 8 3 6
836 WMIECLjSSTJbMPjbfiCAP
837 FOftMAT(34H SEPAKA1ION FACTORS I MPF.OPER. SKP= E12.4,7H SMCAP = E12-4)
838 CONTINUE
IF(T-8b.)7110,7110i7120
7110 XKh = .13E-14*(2.3>**C(1-32. )/lg. )
GO TO 7150
7120 IFC113.-1) 7130*7130*7140
7130 XKVv = 4.fe)E-1 4*( 4.0)**( (T-l 13. )/45. )
GO TO 7lb0
"7140 XKW = 1.47E-14*(2.72)**((T-86. )/27. )
7150 CONTINUE
IFCT-70.) 7210.7210.7220
7210 XCASO = « 1 .9 + 0.2*(T-40.)/30.)/136.)**2
GO TO 7270
7220 IFO-103. )7230. 7230> 7240
7230 XCASO = 2.4E-4
GO TO 7270
"7240 IFU-130.) 72b0* 72b0, 7260
72b0 XCASO - (2.1b*C0.79b3488)**((1-100.)/30.)/136.)**2
GO 10 7270
7260 XCASO = (1.71*(0.7427>**((1-130.)/30.)/136.)**2
7270 CONTINUE
C03K1 =10.** ((12343.5/CT+491.7)>-26.7bSSb+.023568*
-------
FG=FGMAX
GO TO 650
C
C CHECK TO SEE IF HMAX LIMITS HC03 CONCENTRATION
C
b!0 C1F2=ACA*AHC03
C2F2=ACA+AHCG3
20 C3F2=C03K1 /(CAUI*HC03D )
C
C SOLVE CASE 1
C HCO=HMAX, FLAM FROM EON 2
S=l . -(C3F2*HMAX*(XKH-HMAX) )
S=C2F2**2-4.*C1F2*S
IF (S) 1 10* 130; 130
110 MITE(L, 199)
WRITECL,MS)
lib FORMAT (31H NEGATIVE SQUARE ROOT IN CASE 1 )
GO TO 1000
1 30 FG=(-C2F2 + S(OKT( S ) ) / ( 2 . *C 1 F2 )
CAC=SUMCI+FG*SMCIA -XK1*HC03D *(!.+
1 FG*AHC03)/(XKl+HhAX)-i-HiviAX-XKVj/HMAX
C
C IF CAC IS NEGATIVE=EhKOK
C
IF (CAC) 140, 1 bb* 1 bb
140 U'kITE(L, 199)
IvRITECL* Ib0)
i b0 FORMAT OIH STATED PROBLEM is NOT FEASIBLE
1 41H MINIMUM bf>iINE STREAM FH IS VIOLATED
2 24H EVEN AT ZERO ACID FEED )
GO TO 1000
Ibb IFCFG-FGMAX)b20,b20,blb
bib FG=FGMAX
C
C FIND H TO SATISFY HC03 FOR CURRENT FG
C (THIS IS RE-ENTRY FROM 671.) FIND CAC REQUIRED.
C
b20 B=XK1
C = -(FG*FG*ACA*AHC03-i-FG*(ACA-fAHC03)-i-l . )*CADI*HC03D /C03K
)-l .£-35)640,640,630
C . . . ... .._. ._..
C IF HCO IS 'ESSENTIALLY' ZERO THIS CONSTRAINT CAN'T BIND
C . . . .
630 CONTINUE
_. _CAC.^SJJMCI+FG*SMC_IA.. -XK 1 .*HQPJ.D_.._*.( 1 • ,+ F.O>*AHCQ3,)./XXK 1+HCO. >
"7 +HCO-XKW/HCO
IF(CAC)640,640,650 .
640 CAC = 0.
IACID = 2
650 ~ CONTINUE
C _ .._
- 64 -
-------
C NOW HAVE FG*CAC DESIRED FOR HC03 (AND HMAX).
_C MAY HAVE TO ADJUST DOWN TO SATISFY SQ4.
C (THIS IS FG RE-ENTRY FOR CAC=0.)
_c ;
XLM=FG*FG*ACA*AS04*CADI*S04DI+FG*C(ACA+AS04)*CADI*S04DI
1 +CADI*ACA*CAC)-*-CADI*S04DI+CADI*CAC-XCASO
IF(XLM)470*470,660
_C
C IF XLM IS NOT>0 S04 CONSTRAINT IS SATISFIED.
C
C TRY REDUCING FG BUT NOT BELOW FGMIN
_C
660 IF(FG-FGMIN)665»665*670
665 WRITE(L*666)
666 FORMAT(53H SULFATE NOT SATISFIED AT MINIMUM WASTE CONCENTRATION)
GO TO 470
670 FG = FG*GMP
671 GO TO (520*650),IACID
470 CONTINUE
CMAX = GDI * ( l.+FG)
Q7=Q1/(FG/FRACT+1•)
Q3=Q1-Q7
Q9=Q7*CAC/34.7
C THE CONSTANT ABOVE* 34.7* IS BASED ON 93%H2S04
C
199 FORMAT(42H ERROR IN ELECTRODIALYSIS DESIGN PROBLEM )
FT = Q3 * 43.8079
_C XFTNA=FLOW PARAMETER- FT/NA, LTR/SEC/CM
XFTNA=2.412E-3
CLAM = 0.1*EXP(.011*(T-68.))
VISC = .01*EXP(-.01l*(T-68.>)
RM = (RAM»RCM)*EXP(-.01 l*(T-68. ))
SOL = CDI/CMAX
CURL I =_CU_RK_* T_PEL*_*. 4/ ( < VI SC** . 9 3) * ((1. -. 393*
1TDEL)*COTHE)**.6>
RN=XFTNA*1000/( C 1... -»-TDEL )*VI SC*COTHE )
HDIAM=TKNS*(1-3.14159265*TDEL/8)/(l+TDEL)~
PGRAD=VISC**2*( .023*RN-*-6.98E-5*RN**2)/(HDIAM**g.75»COTHE)
XICD=CURLI /TKNS*XFTNA**.6
-FRACT )+SOL)
= 9.65E+07*(FRACT*RM-»-TKNS/(CLAM*CDI)*
1 ALOG (1«/((1.-FRACT)*(1»-FRACT*SOL))))/(E*RPO*(^.-FRACT))
XKP=96.5*CDI**2*FRACT*(1.-FRACT)/E
XNAP= XKPP/XICD*FT
EDPWR=XKP*RPO*XICD*FT/1000.
PMH20=PGRAD*XNAP/FT*XFTNA/100 . ;
STKHD=3.281*PMH20+30.
TTA = XNAP * 1.0764E-3
CURDE =96500.*CDI*FRACT*FT/E/TTA
ELL=TKNS/TDEL*.3937 .
XNA=FT/XFTNA
NSTK=TTA/917.+.9
XNSTK=NSTK
NRECT=XNSTK/4.+.5
- 65 -
-------
13
14
IF(NRECT) 1 3* 13, 14
NRECT=1
CONTINUE
WIDTH=XNA* .3937/(XNSTK*275. )
CSTK-a7330« *XNSTK+( 1 6800 . +'JC) «08*Q3*( STKHD)) *( •'' * •* f -.2)
CDC=»°7/..2b*XNSTK+lbl .41*^ r >"'R
M 102. 5*XNSTK*(^'f TK/4. >**<-. 17)
CAPsCSTK+CDC+CAUX
CREP=1900.*XNSTK
COT=83. /i/'FM^o-* i j?68.*Ql**« 725+. 0035*
r
r
C
C
r
OPERATING COSTS IN CENTS/KGAL ARE. RA^n ov TN'pin; FLOW SO THAT.
COST IS PER I'IVTT n(r i-^FTE WATER TREATED. TO FIND COST PET: I 'N' T T
OF USABLE WATER PRODUCED REPLACE r' 1 T *• ^HF FOLLOWING EXPRESSIONS
COT=COT/Gl/3650.
" CREP=CREP/Q 1/3650.
ABLDG=152.+127.5*XNSTK+20.
1000 CONTINUE
Q9 = C&)9/Q7)*1000.
._IFCPRNT)30, 30*21
21 CONTINUE
950 FORMATC 30H CONCENTRATION
951
952
"953"
954
"955
956
957
958
"9b9
960
961
962
FORMATC
FORMATC
FORMATC
FORM AT C
FORMATC
FORMAT(
FORMATC
_FO,RMATC
FORMAT <
FORMATC
FORMATC
FORMATC
WRITECL
30H
F8
F8
F8
FB"
F8
F8
r.8.
F8
F8
F8
F8
* 9
9 9
•0'F.9
.0*F9
.0*F9
.0*F9
.0,F9
.0,F9
•0'F9
.0* F9
.0*F9
.0jF9
.0*F9
50)
IN
.0»F
i'0*F
.0*F
.0*F
.0,F
.0*F
.0*F,
.0, F
• 0*F
.0*F
•0.
1
0
•
10.
1
1
1
0
0
0
10
1
1
1
1
/
0
0
0
0
/
2*
2;
OUT
7H
10H
•_2*8H .
•
•
•
•
,
2*
2*
2>
2*
2,
10H
9.H
9H
J2H
8H
•2/8H
•
}
2>
10H
36*G3 + 2.*(Q3*1«E6>**.33+ 9 52 ) PNA* PONA, ANA
WRIjrECL*953>PKjPOK*AK .. _
WRITE"CL/954)PCA*POCA/ACA
__
WRI TEC L956 OPNH 4PONH4* ANH4
WRIJJLC Lj> 9 5 7 ) PCL ,, P.QCL *_A£L_
WRITE(L>958)PHC03*PbHCO>AHC03
WRITE(L,960)PN03*PON03>AN03
WKITE(L*961)P.P04»PP.PQ.4J!AP0.4_
WRITE(L,963)TKNS
- 66 -
-------
lvKITE(L,964)TDEL
WRITET
WRITE(L,966)HCM
lvkITE(L,?67)RAM
WRITECL,968)E
963
964
965
966
967
"968
969
970
971
'972
973
T74"
30
toKITEGAMMA
WRITECL,971 >GMAX
WRITECL,972)HMAX
WRITECL,973)CL'RK
toRITLCL,974>CKtoH
FORMATCF10.3,20H
FO'RMA'TC Fife). 3/35H "
FGKMATCFlia.3*28H
FGRMATCF10.3,43H
FORMAT(F10.3*43H
FORKATCF10.3,21H
FCRMAT
RCM=CATION MEMBRANE RESI STI V1T Y, OHM- SO . CM. )
RAM=ANIOl\ MEMBRANE RES I SI I VII Y, OHM-SQ . CM. )
E=CbKRENT EFFICIENCY )
COSINE OF MESH AN GLE, 1 HE! A- )
GAMMA = DESI RED RATIO Q7/Q3 )
GMAX=MAXIMUM RATIO Q7/Q3 )
"MAX HYDROGEN ION CONCENTRATION, WASTE STREAM
LIMITING CURRENT DENSITY PARAMETER, CURK )
COST OF ELECTRI CI TY, S/KVvHR //)
WRIJE904)NSTK
toRnENRECT
WRITE(L*920)ABLPG
WKl"TETL/906> WIDTH
UR1TE(L,908)CSTK
WRITE (1/909) CDC
WRnECL,910)CAUX
WR I_TE_(L * 9 1 1 ) CP WR
WRlfE(L,9l"2)CREP
WRITE(L.,913)COT
WRITYCL,919)OP
WRITE
-------
9f3 FORMAT<32H OTHER OPERATING COST-CENTS/KGAL F15.4)
91 ''
916
917
918
9 19
920
1001
7
FORMAT (32H
FORMAT ( 32H
FORMATC32H
FORMAT <32H
FORMATC32H
FORMAT ( 32H
FORMATC 32H
GO TO 8
CONTINUE
CALL "EXIT
END
INPUT STREAM-MGD
PRODUCT STREAM-MGD
WASTE STREAM-MGD
ACID FEED-GAL/KGAL WASTE
TOTAL CAPITAL- £
TOTAL OPERATING COST-CENTS/K GAL
BUILDING AREA-SO.FT.
Fl
Fl
Fl
Fl
Fl
Fl
Fl
b.
5.
b.
b.
5.
b.
5.
4)
4)
4)
4)
4)
4)
- 68 -
-------
THE FOLLOWING IS THE DATA FILE USED TO PRODUCE THE SAMPLE OUTPUTS*
130. 15.
• 79 1.14
10. 500.
.0833338.
100. S00.
.0833332*
1 . 500•
•0833332.
60.
1 .25
• 1
.001
.1
.001
.1
.001
25.
.98
.1
.707
.1
.707
• 1
.707
20.
1 .15
70.
1 .
70.
1 •
70.
1 .
165.
1 .25
40.
40.
40.
300.
.77
40.
40 .
40 <
100.
1 .11
.9
.9
.9
10.
1 .22
.01
.01
• 01
25.
.72
82.
82.
82.
- 69 -
-------
VARIABLE NAMES: DECISION MATRIX
THE FOLLOWING VARIABLES MUST BE PROVIDED AS INPUT TO THE
ELECTRODIALYSIS SUBROUTINE VIA THE DECISION MATRIX,DMATX.
THE NAMES LISTED HERE ARE THE FORTRAN VARIABLE NAMES USED WITHIN
THE SUBROUTINE.
TEXT PROGRAM
NAME NAME
PPMOU
TKNS TKNS
TDEL TDEL
T
RCM
RAM
E E
CURK CURK
GAMMA
GMAX
HMAX
COTHE
ANA
AK
ACA
AMG
ANH4
ACL
AH CO 3
AS04
AN03
AP04
NOMINAL
VALUE
500
. 1
. 1
70.
40 .
40.
.9
82.
.001
. 707
. 79
1.14
1.25
.98
1.15
1.25
.77
1.11
1 .22
.72
DESCRIPTION
TOTAL DISSOLVED INORGANIC SOLIDS*
PRODUCT STREAM.* PPM
MEMBRANE SPACER THICKNESS, CM.
RATIO OF SPACER THICKNESS TO MESH SPACING
INPUT STREAM TEMPERATURE, DEG-F .
CATION MEMBRANE RESI STI VI TY, OHM-SG . CM.
ANION MEMBRANE RESI STI V I TY, OHM- SQ . CM .
CURRENT EFFICIENCY
LIMITING CURRENT DENSITY PARAMETER
DESIRED WASTE/PRODUCT STREAM VOLUME RATIO
MAXIMUM ALLOWABLE WASTE/PRODUCT
STREAM VOLUME RATIO
MAX ALLOWABLE HYDROGEN ION CONCENTRATION
IN WASTE STREAM
COSINE OF MESH ANGLE,THETA
SEPARATION FACTOR, SODIUM
SEPARATION FACTOR* POTASSIUM
SEPARATION FACTOR, CALCI UM
SEPARATION FACTOR, MAGNESI UM
SEPARATION FACTOR, AMMONI UM
SEPARATION FACTOR, CHLORI DE
SEPARATION FACTOR, BI CARBONATE
SEPARATION FACTOR, SULFATE
SEPARATION FACTOR, NI TRATE
SEPARATION FACTOR, PHOSPHATE
t
- 70 -
U.S. GOVERNMENT PRINTING OFFICE:1970-0-408-309 516
------- |