EPA-600/2-75-062
December 1975
Environmental Protection Technology Series
psKs^s^Sf^!)^!??*^^
^m'^^
^^
:OMBINED SEWER OVER FLOW REGULATOR
unicipal Environmental Research Laboratory
-------
• lRe^search''regbrts of the office of Research an3 Development,
tf.i;S."'Environmental'/?^ SgencyJ "Gave Been grouped" Jn"t0"
" five' "series' These five broad"' categorTeis were e"s taoTSCined to
facilitate further development and application of environmental
techno logy"!1' '" 'El^ina'tion' 0*'"f" 'TtraSit" iona 1 grouping" was consciously
planned to foster technology transfer and a maximum interface in
related fields. The five series are:
1. Environmental Health Effects Research
2, Environmental Protection Technology
3. Ecological Research
4. Environmental Monitoring
5. Socioeconomic Environmental Studies
This report has been assigned to the ENVIRONMENTAL PROTECTION
TECHNOLOGY series. This series describes research performed to
develop and demonstrate instrumentation, equipment and methodology
to repair or prevent environmental degradation from point and non-
point sources of pollution. This work provides the new or improved
technology required for the control and treatment of pollution
sources to meet environmental quality standards.
iiwill in1" IIPl1!11 Hi nil
i mi! ii nil nil iiiiii|iiiii Nim iiiiiiiiiiiii i it i iiiiii i in m i iiiiiiiiiiiiii n ill it i in wnini i iii'iii inn111 111 niii in ill
This document is available to the public through the National
Technical Information Service, Springfield, Virginia 22161.
it i, i vi .r *:•' Ki""11!1!.; ,,'1'! I'ltiti v; kit. '• .,,;';,;,'<-. iiasrii n-E'ts.! i-i' iiiiM11
i!l!":' i '!,;l:::i ^ :€'.: ;C; S; P"i:*^^^ iiMIIR^^^^^^^^^^^^^^ !«ii:i!!!ilH^^^^^^^^^^^^^
-------
EPA-600/2-75-062
December 1975
THE HELICAL BEND COMBINED SEWER OVERFLOW REGULATOR
By
Richard H. Sullivan
Ralph R. Boericke
Morris M. Cohn
George Galina
Carl Koch
Fred E. Parkinson
T. M. PrUs-Chacinski
James E. Ure
American Public Works Association
Chicago, Illinois 60637
Contract No. 68-03-0272
Project Officer
Richard Field
Storm and Combined Sewer Section (Edison, N.J.)
Wastewater Research Division
Municipal Environmental Research Laboratory
Cincinnati, Ohio 45268
MUNICIPAL ENVIRONMENTAL RESEARCH LABORATORY
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
CINCINNATI, OHIO 45268
-------
DISCLAIMER
This report has been reviewed by the Municipal Environmental
Research Laboratory, U.S. Environmental Protection Agency, and
approved for publication. Approval does not signify that the
contents necessarily reflect the views and policies of the U.S.
Environmental Protection Agency, nor does mention of trade names
or commercial products constitute endorsement or recommendation
for use.
ii
-------
FOREWORD
Man and his environment must be protected
from the adverse effects of pesticides, radiation,
noise and other forms of pollution, and the unwise
management of solid wastes. Efforts to protect the
environment require a focus that recognizes the
interplay between the components of our physical
environment — air, water, and land. The National
Environmental Research Centers provide this
multidisciplinary focus through programs engaged
in
• studies on the effects of environmental
contaminants on man and the biosphere, and
• a search for ways to prevent contamination
and to recycle valuable resources.
This report on the Helical Bend Combined
Sewer Overflow Regulator presents the basis of
design for an effective device to limit pollution of
receiving waters from combined sewer overflows.
Improvement of the quality of the overflow while
regulating the quantity is an important
performance specification. The use of a
nonmechanical device utilizing fluid dynamics is
particularly valuable in light of the importance of
energy conservation.
111
-------
ABSTRACT
A design procedure and method to calculate
settleable solids retention efficiency has been
prepared for a combined sewer overflow regulator,
using the principle of developing helical flow in an
enlarged, curved extension of a sewer. A curved
side overflow weir is used to draw off the clarified
combined sewer overflow.
Hydraulic and mathematical model studies
were used in developing the helical flow concept
and the design procedure. The reports of each
model study are included.
The helical bend principle was originally
investigated by Dr. Prus-Chacinski of London,
England. A second generation regulator was
constructed in 1971 in Nantwich, England. The
present study modified the design to conform to
American wastewater conditions and practices.
It was found that the helical bend combined
sewer overflow regulator is capable of higher
removal efficiencies with less hydraulic head loss
than the swirl concentrator regulator, although the
overall length of the helical unit appears to make it
more expensive than the swirl concentrator.
The helical bend regulator requires the use of
only a mechanical outlet gate and cleaning
facilities; otherwise it is essentially non-mechanical.
This report is submitted by the American
Public Works Association in partial fulfillment of
Contract No. 68-03-0272, between the U. S.
Environmental Protection Agency and the
American Public Works Association.
IV
-------
TABLE OF CONTENTS
Page
Abstract iv
Section I Conclusions, Recommendations and Overview of the Studies 1
Section II The Study 5
Section III Design Factors for Helical Bend Regulator Separator Facilities 18
Section IV Implementation 46
SectionV Glossary of Pertinent Terms — Helical Bend Report 47
Section VI References : 50
Appendix A Hydraulic Model Study 51
Appendix B Mathematical Model Study 81
Nomenclature 114
Appendix C Background Information on the Helical Bend 116
TABLES
Page
1. Velocities in Transitions 24
Helical Bend Regulator (Design Example) 28-32
2. Helical Bend Regulator Dimensions 33
3. Site Dimensions and Areas for Helical Bend
and Swirl Concentrator Regulators 35
4. Typical Head Losses in Helical Bend and Swirl
Concentrator Regulators 37
5. Construction Cost of Helical Bend Regulator 40
6. Swirl Concentrator Dimensions ".....' 41
7. Construction Cost - Swirl Concentrator 42
8. Comparison of Construction Costs — Helical Bend-
and Swirl Concentrator Regulators . 44
9. Recovery Results — Modifications 1-8 57
10. Recovery Results-Modifications 9, 10, 12, 14-18 59
11. Recovery Efficiency at Various Bend Radii 114
-------
FIGURES
Page
1. Isometric View of Helical Bend Regulator 3
2. Helical Bend Regulator, 1967 7
3. Photograph of Nantwich, England, Helical Bend Regulator 8
4. Helical Bend Regulator — Form 3 (Model Layout) 9
5. Discharge Measuring Weirs 10
6. Prototype Gradation for Grit and Settleable Organic Material 12
7. Particle Settling Velocities for Solids in Still Water 13
8. Gradation Curves for Gilsonite® and Petrothene® Used in Model 14
9. Prototype Grit Sizes Simulated by Gilsonite and Petrothene in Model 15
10. Prototype Settleable Organic Material Sizes Simulated
by Gilsonite and Petrothene in Model 16
11. Predicted Prototype Grit Recovery vs Design Flowrate Ratio 18
12. Predicted Prototype Settleable Organic Material Recovery vs Design Flowrate Ratio . 18
13. Recommended Transition Profile 19
14. Effect of Transition Slope 19
15. Design Flowrate vs Inlet Diameter 21
16. Recommended Plan Layout 22
17. Water Level in Regulator Above Weir 23
18. Recommended Cross Sections 25
19. Spillway Channel Profile 27
20. Site Requirements 34
21. Typical Cross Section — Helical Bend Regulator 38
22. Estimated Construction Costs — Helical Bend
and Swirl Concentrator Regulators 43
23. Predicted Separation Efficiency vs Settling Velocity
at Several Flowrates 45
24. Downstream View of Model 52
25. Upstream View of Model 52
26. Solids Injection Vibrator 53
27. Top and Bed Flow Angles 54
28. Floor Flow Angle Measuring Equipment 54
29. Bed Angle Relationships for Straight-Through Flow 55
30. Modification 6 - Tests 23-29 56
31. Modification 4 - Tapered Outlet 58
32. Modification 5 - Narrow Inlet 2D Wide 58
33. Modification 5 - Tapered Bend 58
34. Minimum Width Inlet, l.SDWide 59
35. Narrowest Bend, l.SDWide 59
36. Modification 9 - Tests 38-39 60
37. Modification 10-Tests 40-41 60
38. Modification 12-Tests 46-47 60
39. Modifications 14 & 15-Tests 53-55 60
40. Petrothene Recovery in Model as Function of Bend Channel Width 61
41. Modifications 17 & 18-Tests 60-61 61
42. Inlet Configurations Tested 62
43. Transition lODLong 62
44. Petrothene Recovery in 60° Bend for Various Inlet Configurations 63
45. Transition 15D Long 64
VI
-------
FIGURES (Continued)
Page
46. Scum Board Locations Tested 64
47. Petrothene Recovery in 50° Bend for Various Scum Board Locations 65
48. Petrothene Recovery in 60° Bend for Various Scum Board Locations ....... 65
49. Floor Flow Angles with Scum Board D/2 Deep 65
50. Floor Flow Angles with Scum Board D/6 Deep . 66
51. Cross Section of Offset Weir in Bend 66
52. 60° Bend with Offset Weir 67
53. Upstream View of Scum Board on Offset Weir 67
54. Downstream View of Bend and Offset Weir 67
55. Scum Board on Offset Weir 67
56. Petrothene Recovery in 60° Bend with Offset Weir
and Various Inlet Configurations 68
57. Proving Tests Layout ; 69
58. Solids Recovery in Model 70
59. Velocity Contours on 0° Cross Section 70
60. Velocity Contours on 10° & 20° Cross Sections . 71
61. Velocity Contours on 30° & 40° Cross Sections 72
62. Solids Concentration Sampling Jig 73
63. Simultaneous Filling of Solids Sampling Bottles 73
64. Suspended Solids Concentrations on 0° Cross Section 74
65. Suspended Solids Concentrations on 10° & 20° Cross Sections 75
66.. Suspended Solids Concentrations on 30° & 40° Cross Sections 76
67. Distribution of Solids Which Escape Over Weir 77
68. Floor Flow Angles 78
69. Comparative Average Floor Angles for Recommended Layout 80
70. Cylindrical Coordinate System 82
71. Definition of Surface Shape in Channel Cross Section 83
72. Coordinate Transformation 85
73. Mathematical Representation of Channel Cross Section . 87
74. Channel Cross Section in Cartesian Coordinates 88
75. Transformation of Channel Cross Section to Rectangular Region . . , 89
76. Rotation of Velocity Components '. 90
77. Computational Grid in Transformed Space 91
78. Fluid Flow Solution Summary 92
79. Location of Computational Grid Points for
Concentration Field Calculation 95,
80. Particle Flow Solution Summary 97
81. Comparison of Measured and Predicted Velocity Profiles at 0° 103
82. Comparison of Measured and Predicted Velocity Profiles at 20° 104
83. Comparison of Measured and Predicted Velocity Profiles at 40° 105
84. Comparison of Measured and Predicted Flow Angles on Bottom
Along Channel Centerline 106
85. Measured Values of Relative Petrothene Concentration at 0° 107
86. Comparison of Measured and Predicted Relative
Petrothene Concentration at 20° 108
VII
-------
FIGURES (Continued)
Page
87. Comparison of Measured Values of Relative
Petrothene Concentration at 40° 109
88. Mass Flux Over Weir as a Function of Angular Position Ill
89. Cumulative Mass Flux Over Weir as a Function of Angular Position Ill
90. Comparison of Predicted and Measured Recovery Efficiencies at Discharges
Greater than the Design Value 0.85m3 /sec (QD =30 cfs) 11.2
91. Predicted Separation Efficiency versus Settling Velocity
at Several Flow Rates 113
92. Alternative Design for Future Investigation 121
via
-------
ACKNOWLEDGEMENTS
The American Public Works Association is
deeply indebted to the following persons and their
organizations for the services they have rendered to
the APWA Research Foundation in carrying out
this study for the U. S. Environmental Protection
Agency.
PROJECT DIRECTOR
Richard H. Sullivan
CONSULTANTS
Dr. Morris M. Cohn, P.E., Consulting Engineer
Paul B. Zielinski, P.E.
C. H. DOBBIE & PARTNERS
Dr. T. M. Prus-Chacinski
ALEXANDER POTTER ASSOCIATES, CONSULTING ENGINEERS
Morris H. Klegerman, P.E.
James E. Ure, P.E.
LA SALLE HYDRAULIC LABORATORY, LTD.
F. E. Parkinson
George Galina
GENERAL ELECTRIC COMPANY
Dr. Ralph R. Boericke
Dr. Carl Koch
U. S. ENVIRONMENTAL PROTECTION AGENCY
Richard Field, P.E., Project Officer, Chief, Storm & Combined
Sewer Section (Edison, New Jersey)
IX
-------
AMERICAN PUBLIC WORKS ASSOCIATION
Board of Directors
Herbert A. Goetsch, President
Ray W. Burgess, Vice President
Erwin F. Hensch, Immediate Past President
Robert D. Bugher, Executive Director
Jean V. Arpin
John T. Carroll
Donald S. Frady
Lambert C. Mims
James J. McDonough
Robert D. Obering
John J. Roark
James E. McCarty
Kenneth A. Meng
Frank R. Bowerman
Wesley E. Gilbertson
Rear Admiral A. R. Marschall
APWA RESEARCH FOUNDATION
Board of Trustees
Samuel S. Baxter, Chairman
Ross L. Clark, Vice Chairman
Robert D. Bugher, Secrtary-Treasurer
Richard H. Sullivan, General Manager
Fred J. Benson
John F. Collins
Jean L. DeSpain
Richard Fenton
W. C. Gribble
John A. Lambie
James E. McCarty
Marc C. Stragier
-------
SECTION I
CONCLUSIONS, RECOMMENDATIONS AND OVERVIEW OF THE STUDIES
Conclusions
1. A helical flow device has been
developed and demonstrated, as an expanded
section of a combined sewer line, for the dual
purpose of removing settleable suspended
solids from combined sewer overflows, and to
serve as an overflow spill weir regulator for
the discharge of clarified effluent to receiving
waters or to retention and/or treatment
facilities.
2. The helical bend regulator can operate
efficiently as a flow-through, nonmechanical
channel bend device in removing settleable
solids and, thereby, reducing the pollutional
impact of periodic overflow incidents.
3. The solids removed from the combined
sewer flows may tend to shoal in the deep
invert section of the helical bend chamber
during storm flow conditions.
4. Model tests showed that these deposits
remain dormant during such high-flow periods
and are then moved or self-scoured out of the
invert by dry-weather flows and discharged,
via a foul sewer outlet, to the interceptor
sewer and downstream wastewater treatment
works. Similar conditions should prevail in
prototype helical bend installations but the
scale representation of combined sewer solids
is not confirmed and the self-scouring action
experienced in the model studies must be
demonstrated under prototype operations.
(Note: A prototype unit in Nantwich,
England has not experienced the buildup of
solids during operation.)
5. Washing or flushing facilities, and
regular maintenance procedures should be
provided in prototype helical bend
installations until the removal of solids
deposits by dry-weather scouring action has
been confirmed.
6. Uniformity of flow as it enters the helical
bend section of the regulator-concentrator
facility is of utmost importance to the
successful operation of the secondary flow
actions induced in the chamber. .
7. A straight section of conduit, ahead of
the helical bend section, is required to assure
effective solids recovery rates. This straight
section length should be five times the
diameter of the inlet sewer (5D).
8. If no straight section is provided, the
recovery rates for organic solids in the helical
bend will be reduced approximately 10
percent below the efficiency attained with the
5D straight approach line.
9. Beyond this 5D straight approach
length, greater lengths will have minimal
effect on solids removal efficiencies. A 10D
straight section will further improve removals
by less than five percent.
10. Grit recoveries will not be greatly
affected by changes in the length of the
straight approach section.
11. The straight horizontal transition
length from the inlet pipe diameter, out to
the full helical bend channel width of three
times the inlet diameter, should be at least 15
times the inlet, diameter.
12. The transition section must be roofed
or covered to provide the proper flow
conditions for secondary motion and helical
action in the bend section. The cover or roof
should rise from the crown of the inlet pipe
to the hydraulic gradient at the end of the
transition section.
13. The drop through the structure can be
confined to the slope of the incoming sewer;
although greater slopes to encourage
self-scouring should not present any
problems.
14. The length of the bend is determined
by the inlet diameter. A radius of 16D should
be used.
15. A 60 degree segment is adequate to
develop the full effects of the helical flow
principle.
16. Most of the floatable solids entering
the helical bend regulator will be retained in
the chamber by the scum board or dip plate,
located ahead of the overflow weir, during
storm flow conditions. As the discharge rate
decreases and the surface of the combined
-------
flow in the helical chamber recedes, the
floatables will drop with the flow level and be
removed with the channel bottom deposits
during dry-weather flow.
17. A design procedure has been
developed for use in dimensioning the
geometrical configuration of full-scale
prototype applications having specific flow
rates, solids characteristics, and percentages of
dry-weather flows diverted to interceptor
sewers and downstream treatment facilities.
The design criteria are based on a maximum
flow of three times dry-weather flow through
the foul sewer outlet at the end of the bend
section.
18. Mathematical techniques have been
developed for computing the particle and
liquid flowfields in open channel bends.
Recommendations
1. The hydraulic model studies and the
computer-mathematical simulation of the
helical bend combined sewer overflow
regulator indicate that this flash method of
solids removal, without use of mechanical
appurtenances, can produce excellent
efficiencies with reasonable size units in
combined sewer systems. The results are so
promising that they should be confirmed in a
full-scale prototype installation at the earliest
possible opportunity.
Such an installation in the field will
provide confirmation of the pattern of solids
deposition in the deeper channel portion of
the helical bend, located along the inner
circumference of the bend section, and of the
ability of dry-weather flows in this restricted
deep channel to self-scour the deposited solids
into the foul sewer outlet.
The full-scale field prototype installation
should be fully instrumented to record
performance under actual combined sewer
conditions. Facilities should be provided for
sampling the stratified bed load of the inlet
flow and the overflow by means of vertical
arrays of sampling probes. The flow over the
weir should be sampled at several locations to
determine the distribution of solids lost with
the overflow, as a function of position along
the side weir length. A full regimen of
laboratory control studies could be planned
and instituted to explore the helical flow
patterns created in the prototype unit under
normal service conditions. Out of such a pilot
prototype investigation will come operating
experience that can assist in the application of
this promising and economical device to meet
the 2 Q challenge of effective
regulator-overflow control — control of both
quantity and quality of combined sewer
discharges to receiving waters.
Overview
A study1 conducted by the American
Public Works Association (APWA) Research
Foundation in 1969-70, Problems of
Combined Sewer Facilities and Overflows,
disclosed that combined sewer
regulator-overflow devices on the American
continent did little to improve the quality of
liquid spills from these facilities into receiving
waters. The environmental impact of these
discharges into the nation's waters was
characterized as a pollutional paradox in view
of current concern over the secondary
treatment and even higher degrees of
purification of sewage and industrial
processing wastewaters.
The findings of the 1969-70 investigation
have catalyzed subsequent studies on how
economical and efficient devices can be
utilized to eliminate part of the suspended
solids contained in admixed sewage storm
flows in combined sewers prior to the
wasting of major volumes into surface water
sources. In short, the endeavors have been
aimed at making such regulator-overflow
facilities perform two functions: Control the
flow rate of combined sewer overflows and
improve the quality of such spill volumes.
A study2 conducted by APWA, The Swirl
Concentrator as a Combined Sewer Overflow
Facility, described the development and
demonstration of a nonmechanical swirl
concentrator as an economical means of
removing suspended solids from combined
sewer overflow wastewaters.
To supplement the search for workable
means of accomplishing the clarification of
combined sewer overflows, and thereby
-------
INLET
CHANNEL FOR
OVERFLOW
WEIR
OUTLET TO
STREAM
TRANSITION SECTION
15D
STRAIGHT
\ SECTION
5D
NOTES:
1. Scum baffle is not shown.
2. Dry-weather flow shown in channel
OUTLET TO PLANT
FIGURE 1 ISOMETRIC VIEW OF HELICAL BEND REGULATOR
reducing the pollution load imposed by such
periodic discharges, the current study has
investigated the feasibility of achieving the
removal of solids from combined sev/er flows
by means of a helical bend incorporated into
an expanded section of such a sewer system.
This device would utilize the secondary flow
patterns produced by helical hydraulic
phenomena to deposit solids in the invert of a
helical channel section, adjacent to the inner
curve of the conduit, and to allow a clarified
effluent to overflow a side spill weir for
discharge to receiving waters or to retention
and/or treatment facilities.
Figure 1. Isometric View of Helical Bend
Regular, shows the important elements of
the facility. The device may require only one
mechanical part, a control on the outlet to
allow a predetermined rate of flow to the
interceptor sewer if this advanced type of
positive control is desired.
Hydraulic model studies carried out at
the LaSalle Hydraulic Laboratory at
Montreal, Canada, and mathematical
model-computer studies performed by the
Re-entry and Environmental Systems Division
of the General Electric Company at
Philadelphia, Pennsylvania, have resulted in
the development of a helical
regulator-separator device. Searching
investigations of the liquid and solids flow
patterns through a helical bend section have
been carried out, preceded by a roofed- or
covered-transition section length 15 times the
inlet sewer diameter, expanding the
flow-through conduit to a diameter three
times that of the inlet sewer followed by a
straight-flow section ahead of the helical bend
channel having a length of ten times the inlet
pipe diameter.
The studies have demonstrated the
efficacy of the helical bend principle in
achieving a flash separation of solids from the
combined sewer flow and the diversion of a
foul flow of up to three times the dry-weather
flow into the interceptor sewer and thence to
downstream treatment facilities.
Investigations have shown the ability of the
-------
helical separator to act as a nonmechanical
overflow regulator by means of a side weir
overflow arrangement along the outer radius
of the treatment chamber.
The results of helical flow patterns for
these purposes have been gratifying. They
confirm and augment previous findings
described by Dr. T. M. Prus-Chacinski in a
thesis at the Imperial College of Science,
London, England;3 demonstrated in a
prototype installation at Nantwich, England;
and outlined in a consulting report to APWA
by Dr. Prus-Chacinski as part of his advisory
services in connection with the current study.
As a result of hydraulic and mathematical
studies, basic design data have been evolved
for use by designers in applying the helical
regulator-separator to prototype installations
in combined sewer systems.
-------
SECTION II
THE STUDY
The study was designed to investigate the
efficiency of pipe and channel bend
configurations in separating suspended solids
in combined sewer flows and diverting such
concentrated solids to the fraction of the
volume which is carried by the interceptor
line to wastewater treatment facilities. The
studies investigated the solids separation
patterns induced by the helical motion
established in the flow-through bend sections
of pipe or channels, and the consequent
clarification of the overflows from the helical
device into receiving waters or into holding or
treatment facilities for the overflow volumes.
The investigative procedures covered both
hydraulic model configurations and results,
and mathematical model evaluations.
Dr. T. M. Prus-Chacinski reported4 at a
1967 Symposium on Storm Sewage Overflows
sponsored by the Institution of Civil
Engineers on the results of a study he had
conducted utilizing secondary motions in a
bend to separate solids from the moving flow
of sewage. In recent years, Dr. Prus-Chasincki
and the firm with whom he is associated, C.
H. Dobbie and Partners, designed a combined
sewer overflow regulator for the City of
Nantwich, England. A hydraulic model study
was conducted by the Mersey and Weaver
River Authority. Financing of the
construction was, in part, by the Construction
Industry Research and Information
Association.
At the initiation of this study, C. H.
Dobbie and Partners were engaged to develop
an initial design for the hydraulic model
which would incorporate the tentative
findings from the operation of the Nantwich
unit and adapt the design to United States
practices, i.e., to provide a maximum of three
times dry-weather flow to downstream
treatment facilities. Appendix C. Background
and Preliminary Design, is the C. H. Dobbie
and Partners report which established the
essentially third generation configuration for
the hydraulic and mathematical model tests.
The flow pattern phenomenon utilized in
the helical bend . principle involves
liquid-solids separation which occurs when
flow moves around a bend and creates
secondary currents which tend to move the
water surface toward the outside of the bend
and to induce bottom current angles which
move to the inside of the curve. This
separation principle is well known in river
sections where the deepest part of the channel
is usually scoured on the outside of any bend
and finer material is lifted by the bottom
currents and deposited along the inside
shoreline.
The hydraulic model developed for the
present studies was intended to impose the
opposite geometry from that experienced in a
natural stream; the deepest section of the
helical bend line was located at the inside
circumference. The intent was to provide the
deepest part of the bend cross section where
the solids removed by the secondary flow
would have opportunity to settle out and be
swept by the bottom flow to the outlet end
of the bend and thence to downstream
wastewater treatment facilities.
Appendix C reviews the established facts
about helical motion in single bends, the areas
where uncertainty exists, and the basis for the
initial design of the hydraulic model.
This study confirmed the findings outlined
by Dr. Prus-Chacinski in a preliminary report
to the American Public Works Research
Foundation on the model of helical combined
sewer overflow regulators. Selected points
from Dr. Prus-Chacinski's report follow. They
are of significance in the translation of the
hydraulic and mathematical model findings
into practical design guidelines. They are used
to preface the design factors evolved by both
the hydraulic and mathematical study personnel.
1. All bends investigated to date (1973)
have not been longer than 180°.
2. It has been established that in such bends
in rectangular channels prominent helical
motion exists in the form of one helix, in
which the direction of flow, about
one-third of the depth near the bottom, is
toward the inner wall; and in another,
about two-thirds of the depth near the
top is toward the outer wall.
-------
3. The angle between the circumferential
direction and the direction of the
transverse flow is very much greater near
the bed than near the surface; therefore,
the resulting helix is assymetric and very
much flattened toward the free surface.
4. In such a helix the strength of the helical
motion in a rectangular channel is
proportional to the angle of the direction
of flow next to the bed.
5. Such a single helix changes the velocity
distribution, both in vertical and in
horizontal along the bend; around 60° to
80° of the bend the slow bed water is
upturned to the top. There the velocity
distribution is different than in a normal
straight open channel: the maximum
velocity occurs in the lower half of the
depth.
6. Around the same region the pattern of
helical flow changes (at least in a
rectangular channel) and a second helix in
the opposite direction starts to form at
the outer wall. This second helix
gradually grows towards the end of the
bend 180° long.
7. Around the same region the area of
maximum velocity across the channel (at
least in a rectangular channel), which
almost from the beginning of the bend is
at the inner wall, shifts to the outer wall
and remains there up to the end of a
180°longbend.
8. The maximum strength of the primary
helix occurs at about 40° of the bend and
decreases very slowly to the end of the
180° long bend.
9. The described pattern of the helical
motion in a rectangular channel can exist
only if the flow at the entry to the bend
is completely free of any transverse
currents. This is a fundamental
requirement.
10. The strength of the helix is relatively
insensitive to the velocity distribution
across the channel at the entry to the
channel (at least in a rectangular
channel).
11. The strength of the initial helix is much
more sensitive to the velocity distribution
in the vertical at the entrance to the
bend. The initial strength can be
increased by artificially distorting the
vertical velocity distribution at the entry
to the bend in such a way that the
velocity near the surface is increased.
12. The initial strength can be decreased by
artificially depressing the maximum
velocity in the vertical at the entry.
13. Such changes of the strength of the
helical flow disappear around 60° to 80°
of the bend; from there on, the motion
returns to normal helical motion.
14. Introduction of an artificial helical flow
at the entry to the bend can produce an
infinite variety of patterns of helical
motion. The strength of one helix can be
slightly increased or decreased.
Introduction of more than one helix at
the entry to the bend can result in
patterns of two or three helices. Such
pattern persists for at least 45° along the
length of the bend and probably for a
greater length.
15. Helical flow similar in shape has been
experimentally documented in the bends
of rectangular channels with
depth-to-width ratio from about 0.12 to
about 1.0; triangular channel with the
symmetrical apex and with the apex near
to the outer wall; trapezoidal channels;
parabolic channels.; eliptical channels;
semicircular channels; and rectangular
channels with semicircular bottom. The
strength of a primary single helix is
probably the greatest in semicircular
channels.
Figure 2, Helical Bend Regulator - Form
1, as investigated by Dr. Prus-Chacinski,5 was
presented at the 1967 Symposium. Figure 3,
Helical Bend Regulator - Form 2, as
investigated by the Mersey & Weaver River
Authority is a photograph of the facility
designed for Nantwich, Figure 4, Helical Bend
Regulator - Form 3, as investigated by
American Public Works Association, is the
configuration chosen for the initial model
tests. A review of these three figures reveals a
change to simplicity of shape and an
improved method of takeoff to the
interceptor sewer to prevent the outlet from
being plugged by rags and debris.
Starting with the configuration in Figure
4, concurrent hydraulic and mathematical
model studies were conducted.
-------
The basic structural features of
significance in the helical bend model are: the
inlet from the entrance sewer section to the
device; the transition section from the inlet to
the expanded cross section of the straight-run
section ahead of the bend; the overflow side
weir and scum board, or so-called dip plate;
and the foul outlet for the concentrated solids
removed in the secondary flow pattern,
together with the means for controlling the
amount of this underflow going to the
treatment works.
The hydraulic model was constructed at
the LaSalle Hydraulic Laboratory.
Configuration changes and internal auxiliary
flow control and effluent improvement
modifications, as well as detailed studies of
the helical flow patterns developed and the
suspended solids removal efficiencies
achieved, were -studied. The report of the
hydraulic model studies is presented as
Appendix A to this report.
The mathematical model and computer
simulation of the device was conducted by
the General Electric Company. The report of
the mathematical model study is presented as
Appendix B to this report.
The ultimate purpose of both studies was
to correlate and confirm findings by both
hydraulic and mathematical means and to
develop design criteria which will enable
engineers to utilize the helical flow principle
for solids removal from combined sewer flows
and to properly regulate the overflow of
clarified wastewater to receiving waters or
points of retention and/or treatment. The
design basis is presented in Section III.
Hydraulic Model Studies
The transition from the sewer line to the
overall helical system was made of Plexiglas,®
20 diameters (20D) long to assure proper flow
deceleration. This transition section was
followed by a straight-flow Plexiglas section
of 1OD length to carry the flow to the start of
the bend. It was built of Plexiglas, using a
2.44 m (8 ft) radius from the extension of the
inlet pipe centerline. Since the model was
developed to demonstrate the basic principles
of the secondary motions created in the
helical bend, the curvature was originally
carried around through an angle of 120°. Side
channels were provided on each side to
accommodate the clear overflow discharge
from either the inside or outside
circumference of the bend. Control of either
Channel for normal
flow and heavy solids
Overflow weirs
with dip plates
Siphonic slots
to remove
heavy solids
Pipe for solids to
foul sewer or tank
for first flush
SECTION A-A
Flume
Control
pipe -
SECTION B-B
Overflow wey-s
with dip plates
Control
pipe
SECTION B-B
inlet
Flow.carried
forward to works
Courtesy of Dr. Prus-Chacinski
FIGURE! HELICAL BEND
REGULATOR, 1967
as presented in 1967 at the Symposium in
London, England sponsored by the
Institution of Civil Engineers
-------
Courtesy of C.H. Dobbie and Partners
FIGURE 3 PHOTOGRAPH OF NANTWICH, ENGLAND,
HELICAL BEND REGULATOR
-------
100
120
TYPiCAL SECTION IN BEND
CROSS-SECTION
IN STRAIGHT
Locoted 5m(l5.2')[A;rt
Upstreom of
Transition
Note: location of perforated slide
gate varied from 60° to 120
INLET PIPE
I5.acm(6")
FLOW
FIGURE 4 HELICAL BEND REGULATOR - FORM 3 (Model Layout)
-------
the entire flow through the bend, or just the
small foul outflow was provided by a
perforated slide gate at the 120° position.
The basic floor cross section inside the
transition section, the straight section and the
bend was constructed of polished cement
mortar, as shown in Figure 4. The original
layout could accommodate simple
modifications of the Plexiglas overflow
weirs for the purpose of changing their
positions horizontally or their crest levels.
Similar variations in scum board location and
depth could be provided.
Unmeasured water flows entered the
helical bend complex through the supply pipe
which was fitted with a discharge control
valve and equipment with which to inject
predetermined amounts and types of solids
whose removal characteristics and
performance were to be studied. The overflow
water was diverted to a calibrated Rehbock
weir basin where water levels could be read
and the rate of discharge determined.
Similarly, the foul outflow from the bottom
of the bend channel, containing the deposited
solids, was diverted to a second basin for
measurement over a small 90° V-notch weir.
These control devices are shown in Figure 5,
Discharge Measuring Weirs.
Injection of Process Black ink on the
floor of the bend, applied with hand syringe,
produced clearly defined traces of the bottom
or floor flow angles which could be measured
by means of a protractor. These trace tests
demonstrated the presence of the anticipated
helical secondary flow patterns.
Tests for the efficiency of solids
separation in the helical flow pattern were
performed in accordance with recognized
standard methods. One liter of the selected
model solids was introduced into the water
supply via the supply line. This solids material
entered the bend flow and was exposed to the
separation process. Any solids that were
swept over the overflow weir were captured
on a screen in the measuring basin. The
balance normally was retained as a deposit in
the bend and recovered in the diversion basin
when it was flushed from the bend section.
Any solids contained in the foul outflow
FIGURES DISCHARGE MEASURING WEIRS
during the model run were also captured and
retained in the collection chamber.
The volumes of the two fractions — the
solids passing over the overflow weir and the
solids contained in the foul outflow — were
measured to determine if any material had
been lost. The recovery rate, or efficiency of
solids separation, was computed as the ratio
of, the solids retained in the bend compared
with the original 1-liter input, expressed as a
percentage.
The model scale selected for comparison
with a rational prototype size was 1:6. The
pertinent scale relations were:
Scale
Time
Discharge
1:6
1A/6, = 1/2.45
1/65/2 = 1/88
The prototype combined sewer overflow
in the corresponding 91.5 cm (3 ft) sewer was
computed at 0.85 m3/sec (30 cfs) flowing
full, which would correspond to 9.1 I/sec
(0.339 cfs) in the model. To simplify
operations in the model, a discharge of 10
I/sec (0.354 cfs) was used, giving a storm
discharge of 0.88 m3/sec (31.1 cfs). Selected
discharges in the model were 10, 7.5, 5, and
10
-------
2.5 I/sec (0.354, 0.265, 0.177 and 0.088 cfs),
corresponding to 0.88, 0.66, 0.44 and 0.22
m3/sec (31.1, 23.3, 15.55 and 7.78 cfs) in the
prototype. Dry-weather flow (dwf) was
defined as one percent of the storm flow, or
on the model 0.1 l/sec(0.004 cfs), and in the
prototype 8.81 I/sec (0.31 cfs). The basic
decision was made that during storm flows, 3
dwf would be withdrawn continuously
through the foul outlet at the end of the
helical bend section — 26.4 I/sec (0.93 cfs).
As stated, the preliminary studies of
secondary flow patterns and solids removal
results were carried out with the original
model, as shown in Figure 4. Subsequently,
modifications were made in the bend
configurations and internal ancillary
structural details. In all, the LaSalle Hydraulic
Laboratory report, Appendix A, covers 16
modifications and 59 test runs, involving
various changes in the body of the bend and
other variations in length of transition and
straight-flow sections ahead of the bend
structure.
Over and above the basic structural
configurations of the helical bend composite
provided for the studies, the LaSalle
investigations placed great stress on the
introduction of the proper solids materials
into the inflow water to simulate in a
prototype installation actual combined
sewage.
SOLIDS SIMULATION
In defining the pollutant materials to be
removed from the regulated overflows, two
categories were considered: grit, with specific
gravity of 2.65 and settleable organics with-
specific gravity of 1.2. Details of the
simulation procedures are given below:
Grit
The prototype gradation of the grit
material in sewage which is to be removed in
the helical structure was chosen as shown in
Figure 6, Prototype Gradation for Grit and
Settleable Organic Material. The outside
grain-size limits of 0.2 and 2.0 mm (No. 70
and No. 10 sieve) represent medium and fine
sand as defined by the Unified Soil
Classification System. The specific gravity of
the grit was assumed as 2.65, and the
straight-line grain-size distribution was
selected as a representative average of grit-size
data reported for existing raw sewage influent
at treatment plants. Concentration was
considered as being from 20 to 360 mg/1.
Particle sizes smaller than about 1 mm
(No. 18 sieve) are known to remain
suspended and transported in flowing water
according to equations of the type reported
by Meyer-Peter and Muller,5 or H. A.
Einstein.6
Between 1 mm and 0.2 mm (No. 18 and
No. 70 sieve) the particles are in the transition
zone between the above equations and the
Stokes relation. Since the particles involved in
both prototype and model extended into
both ranges, across the transition zone, the
above equations could not adequately
describe the scale relations.
It was necessary, therefore, to use curves
of particle settling velocities as .shown in
Figure 7, Settling Velocities for Solids in Still
Water. For a given grit size with specific
gravity 2.65 in prototype, the settling velocity
was determined from Figure 7. Based on
Froude's law of similitude, this was divided
by the square root of the scale being
considered to find the required model settling
velocity. By referring to Figure 7 with this
model settling velocity, the model particle
sizes were found for the simulating materials
— Gilsonite or Petrothene.
The physical relations used here can be
expressed as follows:
Model scale = \ = Lp/Lm
where Lp andXm are corresponding lengths
in the prototype and the model respectively.
From Froude's Law, the velocity
simulation is expressed by the equation:
Ys. -
v
m
and
For example, for the scale ratio of
prototype to model of 6:1, the settling
velocity in the prototype should be divided
by the square root of 6, or 2.45. From Figure
7, the settling velocity of prototype grit of
11
-------
U.S. STANDARD SIEVE NUMBERS
34 6 8 IO 16 20 3O 40 50 70 100 140
€
1
\
> A
\
\
OF
\
GANI
V
\
\
\
\
/
cs-^
SG 1.20
\*
\
\
- X
•**"
\
\
\
\
\
\
f»
3.
s
\
\
k
G 2
\
\
\
65
V
\
^\
2 1 O.6 O.4 O2 0.1-
Grain Size in mm
U.S. SIEVE SIZE
4
10
20
40
50
70
417 C *MH»
ol^c mm
5.0
2.0
0.84
0.42
0.3O
0.20
% FINER. BY WEIGHT
GRIT
100
100
63
31
18
0
ORGANICS
100
53
31
17
14
10
IOO
A/N
9O
o^\
8O
-fr\
f\i
•H
&f\ "S*
DO .51
0)
>.
J2
ert %:
il
S?
jtf\
4O
^IA
3vJ
or«
iCVJ
1 A
IO
FIGURE 6 PROTOTYPE GRADATION FOR GRIT
AND SETTLEABLE ORGANIC MATERIAL
12
-------
Reference: Hydraulique et Granulats - J. Larras.
10
0
0
o
o>
1
CO
0.1
/
L/
3
i
y/
_$L
u
XI
•»•
£ &
0.01
0.01
0.1 1
Particle Diameter, mm
1O
FIGURE 7 PARTICLE SETTLING VELOCITIES FOR SOLIDS IN STILL WATER
13
-------
U.S. Standard
8 10
numbers
50 70 100 140
100
FINE
GRAVEL
COARSE
VJ 1 U 1 II O lit III II
MEDIUM
i in
FINE
SAND
FIGURE 8 GRADATION CURVES FOR GILSONITE
AND PETROTHENE USED IN MODEL
0.2 mm (70 mesh) size is 2.2 cm/sec (0.86
in./sec). The model settling velocity is then
0.90 cm/sec (0.35 in./sec). Thus, in the
model, the grit of 0.2 mm size (70 mesh) can
be simulated by 0.70 mm (0.27 in.) Gilsonite,
or 2.1 mm (0.83 in.) Petrothene.
The Gilsonite and Petrothene available
for test work in the laboratory had the grain
size distributions shown in Figure 8,
Gradation Curves for Gilsonite and
Petrothene Used in Model. Practical limits
represented on these curves were chosen
between 0.5 and 3.0 mm (0.02 in. and 0.12
in.) (No. 35 and 6 sieves), and the
corresponding prototype grit sizes simulated
were calculated. The results are shown in
14
-------
1/4 1/6 1/8 1/12 I/IS 1/20 1/24
Model Scale X
Petrothene
Gilsonite
FIGURE 9 PROTOTYPE GRIT SIZES SIMULATED BY
GILSONITE AND PETROTHENE IN MODEL
Figure 9, Prototype Grit Sizes Simulated by
Gilsonite and Petrothene in Model.
This figure shows that for smaller scales,
up to about 1/16, the Gilsonite did;not
cover the larger prototype particle sizes.
However, it was reasoned that if the structure
under study showed a particular recovery rate
for these scales, the larger particles not
simulated would have settled equally as well.
Conversely, the range covered by the
Petrothene includes all of the Lower
Prototype Limit, and drops down to much
finer particles. In fact, as was demonstrated
later, the Gilsonite settled out completely
for all cases tried, so the Petrothene was
used exclusively for the development tests,
and all comparisons are for Petrothene
recovery. This means that the figures quoted
15
-------
are a severe evaluation of the structure; all the
grit in the defined range was removed, and the
comparisons are for even finer particles.
Organics
The settleable organic material desired to
be removed from the flow in the prototype
was assumed to have a specific gravity of 1.2,
and a grain-size distribution as shown in
Figure 7. The same Froude model scale-up
procedure was followed, and the limits
simulated by the Gilsonite and Petrothene
were calculated as shown in Figure 10,
Prototype Settleable Organic Material Sizes
Simulated by Gilsonite and Petrothene in
Model.
The upper limits offered no difficulty; all
the larger sizes were very easily covered by
the Gilsonite. However, there remained a
small section of the lower limit which was not
covered by the Petrothene for scales smaller
than about 1/10. This portion might be
ro.o
1/4 1/6 1/8 1/12 1/16 1/20
Model Scale X
1/24
Petrothene — —
Gilsonite
FIGURE 10 PROTOTYPE SETTLEABLE ORGANIC MATERIAL SIZES
SIMULATED BY GILSONITE AND PETROTHENE IN MODEL
16
-------
considered as lost but, at most, it actually
represents only 10 percent of the prototype
settleable organics at the smallest scales.
Floatables
Only brief tests were carried out to
simulate floating pollutants which were
assumed to have a specific gravity between
0.9 and 0.998, and a size range between 5 and
25 mm. Concentrations of 10 to 80 mg/1 were
assumed. In the model studies, uniformly
sized polythene particles 4 mm in diameter
and a specific gravity of 0.92 were used.
Mathematical Model Studies
The mathematical studies were carried
out on a theoretical model, as shown in
Figure 4. The mathematical model provided
for the same general features as the hydraulic
model, namely, a fixed inlet diameter, a
transition section, a straight section, and a
helical bend that comprised a 60° section.
The configuration of the bend cross section
provided a deeper area at the inner curve, a
set overflow-weir and scum board for
overflow purposes, and an overflow collection
trough. Underflow from the separator was
provided at the 60° location at the end of the
bend.
The mathematical computer-simulation
studies involved detailed numerical
calculation of the liquid flowfield based on
the equations of motion, together with
appropriate boundary and initial conditions.
The velocity field for the solid particles was
approximated by superimposing the settling
velocity on the calculated liquid velocity
flowfield. The spatial variations m the particle
concentration throughout the bend were then
calculated numerically from the mass
conservation equation with turbulent
diffusion. Mathematical relations governing
the scaling of the liquid flow and particle flow
were evolved from the governing equations.
Of final significance was the
mathematical development of comparisons of
calculated mathematical model results with
test data. This phase of the mathematical
model-computer evaluation study compared
the calculated flowfield with test data and
correlated calculated concentration field
findings with the results developed with the
hydraulic model. The findings of the
mathematical model studies were translatable
into design techniques for the development of
prototype installations of the helical bend as
an overflow regulator and as a solids separator
to improve the quality of overflow
wastewaters from combined sewers.
17
-------
SECTION III
DESIGN FACTORS FOR HELICAL BEND REGULATOR SEPARATOR FACILITIES
The hydraulic studies yielded data which
can be utilized to develop full-scale design
values for reasonable prototype dimension.
The LaSalle Hydraulic Laboratory
investigations were on prototype inlet sewer
0.915 m (3 ft) in diameter. The design
discharge QD, was found by taking a nominal
flood flow velocity in the range of 1.20 - 1.50
m/sec (4 tp 5 ft/sec). For the design discharge
actually simulated in the test procedures, 0.88
m3 /sec (31.1 cfs), the flow velocity would be
1.3 5 m/sec (4.4 ft/sec).
Design Guidelines
The designer should use the following
figures for designing the helical separator.
These figures are:
Figure 11 Predicted Prototype Grit
Recovery vs Design
Flowrate Ratio
Figure 12 Predicted Prototype
Settleable Organic
Material Recovery vs
Design Flowrate Ratio
Figure 13 Recommended
Transition Profile
Figure 14 Effect of Transition
Slope
Figure 15 Design Flowrate vs Inlet
Diameter
(Recovery % \
a a o 1
0 0 g |
j-
-
"""""•
\
•5
-__
O.SQD 0D I.5Q0 2QD
Design Flowrate Ratio
Figure 16 Recommended Plan
Layout
Figure 17 Water Level in Regulator
Above Weir
Figure 18 Recommended Cross
Sections
The first step is to determine the design
flowrate. This might be based on a storm
flow expected to occur with a frequency of 5
or 10 years, or some other period, depending
on the local design criteria. Usually larger
flows than the selected design discharge may
occur, depending on the capacity of the sewer
system when surcharged. Figure 11 indicates
that the grit recovery will decrease to 97,
percent with a flow equal to 1.5 times the
design flowrate, and to 93 percent with a flow
equal to 2 times the design flowrate. Hence,
the efficiency of grit removal is not greatly
affected by flows up to twice the design
flowrate.
Similarly, Figure 12 indicates the decrease
in recovery of organic matter with increase in
flow. Thus, for 1.5 times design flowrate, the
efficiency decreases to about 87 percent; and
for 2 times design flowrate, the efficiency
decreases to about 75 percent.
0.5QD QD 1,5 00
Design Flowrate Ratio
Note:
Organic material as defined on Figure 6
FIGURE 11 PREDICTED PROTOTYPE
GRIT RECOVERY VS
DESIGN FLOWRATE RATIO
FIGURE 12 PREDICTED PROTOTYPE
SETTLEABLE ORGANIC
MATERIAL RECOVERY VS
DESIGN FLOWRATE RATIO
18
-------
Inlet Diameter *1D'
Note: For Profile C-C Location
See Figure 16
FIGURE 13 RECOMMENDED TRANSITION
PROFILE
From the foregoing, it would appear that
a considerable increase in flow above the
design flowrate can occur without greatly
affecting the operating efficiency of the
helical separator.
Transition Slope
As shown in Figure 13, the recommended
transition has a length of 15 D and a height ot
D at the inlet and 2 D at the outlet. In the
laboratory model the invert of the transition
was level. In practice this is not desirable
because this would cause deposition of solids
on the invert during periods of minimum
flow. Therefore, the invert should have some
slope. To prevent any surcharge at the inlet,
the top of the transition should be kept level
and the invert given a slope. The slope should
be either the slope of the inlet, or the slope that
will satisfy the hydraulic slope S in the
Manning Equation, whichever is greater.
The resultant hydraulic conditions either
with the invert level or the top level is shown
in Figure 14, Effect of Transition Slope. The
transition with the level top has the following
Hydraulic grade line
Water surface
• Velocity head
2D
T
Velocity Head
Transition loss
D
Transition 15D
Profile with Level Invert
Hydraulic Grade line
• Water surface
• Velocity head
Velocity head
Transition loss
2D
Profile with Level Top
FIGURE 14 EFFECT OF TRANSITION SLOPE
19
-------
advantages: (1) The sewer is not surcharged
upstream of the transition except for loss of
head in the transition, which may be minor;
and (2) the slope will increase the velocity
through the helical separator as the storm
flow subsides which may aid in flushing
deposits out of the helical section. The chief
disadvantage of providing too great a slope in
the transition is that the outlet pipe to the
stream from the helical separator may be
lowered so much that the extension of the
existing sewer cannot be utilized for this
purpose. Therefore, each situation must be
evaluated before selecting the slope. Again as
a minimum, the transition should have the
same slope as the incoming sewer.
Transition Length
The transition length as given in Figure
13 is 15 D. The value of D is selected from
Figure 15, Design Flowrate vs Inlet Diameter.
Assume the designer selects a D of 1.83 m (6
ft) from Figure 15 as appropriate for the
design flowrate. Then the recommended
transition length is 27.4 m (90 ft). However,
assume the existing sewer has a diameter of
1.52 m (5 ft) rather than 1.83 m (6 ft). The
problem is how to effect the connection from
the existing sewer to the transition. The most
logical way would be to extend the transition
to meet the existing sewer while reducing the
area at the same rate as occurs in the
transition. The area of the transition at the
entrance would be 0.785 D2 and at the exit
4.70 D2. These areas are equivalent to squares
with a side of 0.885 D at the entrance and of
2.16 D at the exit. Accordingly, the slope of
the side of the transition would be equal to
1.28 D divided by 30 D, or 0.0426. This slope
has an angle of 2 degrees 26 minutes. Thus, if
the diameter of the existing sewer is 0.30 m
(1 ft) smaller than the selected D, the
transition should be extended by 3.5 m
(11.6 ft).
It would also seem logical to reduce the
length of the transition by a similar process if
the area of the existing sewer is larger than
the area of the transition inlet selected from
the design charts.
Transition Inlet Size
All dimensions of the helical separator are
related to D, the diameter of the transition
inlet. After determining the design flowrate,
the designer should select the inlet diameter,
D, from Figure 15, which shows the simple
scaled-up values according to the Froude Law,
for the design discharge, QD, as well as 1.5
QD and 2 QD . These have been computed
covering the likely range of applicable flood
discharges and pipe sizes that will be
encountered in any prototype installations.
Seldom will the value of D be that of a
standard pipe size. Hence, the designer should
select the nearest D corresponding to a
standard pipe size. If the indicated D falls
between two pipe sizes, the larger D will give
a separator with greater efficiency than the
smaller D. For instance, if the design flowrate
is 2.83 m3/sec (100 cfs), the indicated D will
be 1.45 m (4.75 ft). The designer can select a
D of 1.52 m (5 ft) which is equivalent to a
design flowrate of 3.11 m3/sec (110 cfs), or a
D of 1.37 m (4.5 ft), equivalent to a design
flowrate of 2.40 m3/sec (85 cfs). If the latter
capacity is chosen, the design flowrate will be
18 percent larger than the separator capacity.
From Figure 11, the grit removal efficiency
will be reduced to 99 percent of the total
gtitload. From Figure 12, the settleable
organic removal will be reduced to about 96
percent.
If, in the example given above, the
existing sewer should have a diameter equal to
one of the possible D selections, then it would
be logical to select the D which matches the
existing sewer size. Otherwise, the transition
should be extended as discussed previously.
The overall length of the helical separator
is approximately 37 D including the transition
and straight sections, as shown in Figure 16,
Recommended Plan Layout. If a D of 1.37 m
(4.5 ft), is selected, the length will be 50.9 m
(167 ft), whereas if a D of 1.52 m (5 ft) is
selected, the length will be 56.4 m (185 ft). A
third possibility, assuming the existing sewer
is 1.37 m (4.5 ft) and the D indicated by the
chart is 1.4 m (4.75 ft), would be to base all
dimensions on the indicated D and to extend
the transition according to the method
explained previously. In this case, the
transition would be extended an amount
equal to one-half the difference in diameters,
divided by 0.0426 or 0.91 m (3 ft). The
overall length in this case would be 37 times
D plus 0.91 m (3 ft), or 5.45 m (179 ft).
20
-------
. 12
10
Q
1
a>
.S
Q 4
Jj
"c
10
20 30 40 50
3.5
3.0
2.5
2.0
1.5
2
5
o>
I
Q
;-,
-2
I
.3
Q
i.o
0.5
IOO 20O 300 400500 100?)-cfs
i i I i i i i I I I t
0.5
2345
Design Flowrate
10 15 2O 25 m3/s
FIGURE 15 DESIGN FLOWRATE VS INLET DIAMETER
21
-------
Scum Board
Weir
Note:
See Figure 18 for
sections A-A & B-B
See Figure 13 for
Profile C-C
Inlet Diameter -
-D-
FIGURE 16
RECOMMENDED
PLAN LAYOUT
Thus, the designer is faced with the
choice of three lengths - either 50.9 m (167
ft), 56.4 m (185 ft) or 54.5 m (179 ft).
Obviously the largest helical separator will
provide the most efficient operation.
Velocities in Transition
As a matter of interest the velocities at
the inlet and outlet ends of the transition
were computed for five values of D, from
Figure 17. The D values selected were 0.91,
1.22, 1.52, 1.83 and 2.13 m (3,4,5,6, and 7
ft). The results are shown in Table 1,
Velocities in Transitions. These data indicate
an outlet velocity from the transition ranging
from 0.22 m/sec (0.71 fps) to 0.36 m/sec
(1.17 fps). This compares with the usual
criteria of velocities between 0.23 m/sec (0.75
fps) and 0.38 m/sec (0.75 fps) and 0.38 m/sec
(1.25 fps) in a rectangular grit channel with
velocity control. In general, the outlet
velocities are about one-sixth the inlet
velocities. All velocities are based on the
sections flowing full.
Channel. Slope
In the laboratory model the channel of
the helical separator was level. In practice the
channel should be given enough slope to
maintain a self-cleansing velocity of 0.61
m/sec (2.0 fps) with DWF (average
dry-weather flow).
Assume the following:
D = 0.91 m (3.0 ft)
Design flowrate = 0.85 m3/sec (30 cfs)
DWF = 1 percent of design
flowrate
= 0.008m3/sec (0.3 cfs)
Peak DWF = 0.025 m3/sec (0.9 cfs)
From a chart showing hydraulic
properties of circular sections when the
flowrate is one percent of the full section,
the depth is seven percent of the full depth
and the velocity is 31 percent of the velocity
of the full section. The desired slope is that of
a section when flowing full at a velocity of
0.61 m/sec (2.0 fps), divided by 0.31, or 1.98
m/sec (6.5 fps). From a nomograph of flow for
Manning n = 0.013, the required slope of the
channel for a diameter of 0.91 m (3.0 ft) is
0.48 percent. If the peak dry-weather flow is
3 times DWF the following data prevail:
When slope is 0.48%
Q = 0.025m3/sec (0.90 cfs)
d =0.11 m (0.36ft)
v =0.5 8 m/sec (2.0 fps)
The foregoing assumes a circular section
in the channel when flow is 1 percent and 3
percent of design flowrate. From a visual
comparison of a large-scale section of channel
with a, circular section it is evident that flow
conditions in a circular section will prevail for
the depths of flow considered above. The
foregoing indicates that peak dry-weather
flows should cause no deposition in the
channel.
22
-------
1.2
1.0
0.9
0.8
£ 0.7
1
^ .
'3
g 0.6
< 0.5
(U
0>
3:
0.4
0.3
0.2
O.I
20 30 40 50
1 1 i
0.5
100 200 300 400 500 1000 cf s
-i I I I I 1 I I I I II _
1
2345
FSowrate
10 15 20 25 m/s
Note: Divide Flow/rates
for D = 0.61 m (2 ft)
Curve by 100
FIGURE 17 WATER LEVEL IN REGULATOR ABOVE WEIR
23
-------
TABLE 1
VELOCITIES IN TRANSITIONS
Inlet Design
Diameter Discharge
0.91m
(3ft)
Area
Velocity
Inlet
1.22m
(4ft)
1.52m
(5ft)
1.83m
(6ft)
2.13m
(7ft)
0.85 m3 Is
(30 cfs)
1.84m3/s
(65 cfs)
3.11m3/s
(110 cfs)
4.96 m3/s
(175 cfs)
7.65 m3 Is
(270 cfs)
0.65 m2
(7.0 sf)
1.17m2
(12.6sf)
1.82m2
(19.6 sf)
2.63 m2
(28.3 sf)
3.58 m2
(38.5 sf)
Note: Inlet Area = 0.785 D
Outlet Area = 4.70 D when depth is 2D.
Weir Discharge
Previous research on side overflow weirs
indicates that with a relatively high weir, as
proposed in the helical separator, the usual
weir discharge equations provide a reasonable
basis of design. The usual equation is as
follows:
Q = CLH3/2
where
Q
c
L
H
The
whether
flowrate in m3/sec (cfs)
coefficient
length of weir in m (ft)
head on weir in m (ft)
coefficient C varies depending on
the weir is sharp crested or broad
crested and depending on the head and width
of the weir.
In the laboratory model the weir was
made of Plexiglas and maintained in level
position. The heads on the weir were
measured at the 35 degree position. The heads
on the weir for various values of D and
flowrates are shown in Figure 17. Heads were
obtained from this figure for various flowrates
and the value of C was computed in the weir
equation. The results indicated a range of C
(based on use of U.S. customary units) from
3.3 to 3.7, with most values being 3.4 or 3.5.
Such values are usually associated with
sharp-crested weirs.
Outlet
3.93 m2
(42.3 sf)
6.98 m2
(75.2 sf)
10.9 m2
(117 sf)
15.7m2
(169sf)
21.4m2
(230 sf)
Inlet
1.31 m/s
(4.3 fps)
1.58 m/s
(5.2 fps)
1.71 m/s
(5.6 fps)
1.89 m/s
(6.2 fps)
2.13 m/s
(7.0 fps)
Outlet
0.22 m/s
(0.71 fps)
0.26 m/s
(0.86 fps)
0.29 m/s
(0.94 fps)
0.30 m/s
(1.0 fps)
0.36 m/s
(1.17 fps)
In practice a sharp-crested weir would be
made of a steel or fiber glass plate. The only
advantage of such a weir over a broad-crested
weir is that it is easier to make adjustable.
However, this is considered of minor
importance.
Experience in Great Britain, where side
overflow weirs have been used more widely
than in the United States, favors the use of a
weir with a semicircular shape. This shape
seems preferable for the helical separator.
The coefficient of a broad-crested weir
varies with the width of crest and head on the
weir. For the widths and heads likely to occur
in the helical separator the value of C (for
U.S. customary units) may range from 2.8 to
3.3. The use of a C value of 3.0 (for U.S.
customary units) for design purposes is
suggested. An example follows:
Design flowrate = 0.85 m3 /sec (30 cfs)
D = 0.91m (3.0 ft)
L (weir length) = 18.83 D
= 17.2m (56.5 ft)
Assume flow to plant = 0.028 m3/sec(1.0 cfs)
Q (over weir) = 0.82 m3 /sec (29.0 cfs)
In English units Q = CLH3/2
If C = 3.0 then H = 0.095m (0.32 ft)
The weir height is 1-5/6 D from Figure
18, Recommended Cross Sections. Therefore:
24
-------
Section A-A
m
Note: For Section Locations
See Figure 16
Typical Section B-B
FIGURE 18 RECOMMENDED CROSS SECTIONS
Weir height = 1.68 m (5.50 ft)
Head on weir = 0.095 m (0.32 ft)
Water depth = 1.77m (5.81 ft)
However, to meet the laboratory
demonstrated requirements that the transition
be flowing full, the.water depth should be 2D
or 1.83 m (6.0 ft). Therefore, in this case the
weir height should be a minimum of 1.73 m
(5.68 ft) so that the transition outlet is
flowing full when design flowrate occurs.
Outlet Control
Various methods of controlling the flow
from combined sewer overflow regulators are
discussed in an EPA report.7 This report
indicates that close control of the outlet flow
requires the use of a sluice gate controlled by
a float and actuated by either water power or
an electric motor. On smaller structures where
the use of such devices is not justified, one
method of control is by use of a manually-
operated gate. The intent is to- only operate
such gates to clear them 'of debris or to
change the opening size.
The use of such gates may result in
considerable variation in the flow diverted to
the treatment plant. This may not be serious
when this flow is only a small percentage of
25
-------
the total flow tributary to the plant. To
indicate the possible range in flow, the
following example is based on the use of a
manually-operated gate on the outlet to the
treatment plant. The minimum size gate used
should be 0.20 m (0.67 ft) square but a gate
with a minimum size of 0.30 m (1.0 ft) square
is preferable.
Legend
A = Cross-sectional Area
D = Diameter
V = Velocity
d = Depth of flow
Q = Quantity of flowrate
b = Width of opening
C = Coefficient
DWF = Average dry-weather flow
g = 9.81 m/sec2 (32.2 ft/sec2)
Pertinent Data
DWF = 0.008 m3 /sec (0.30 cfs)
Peak DWF = 0.025 m3 /sec (0.90 cfs)
Try sluice gate 0.30 m (1.0 ft) square
Assume opening 0.10 m (0.33 ft) high
Then A = 0.03 m2 (0.33 sf)
Determine depth upstream of gate when:
Q = 0.025 m3/sec (0.90 cfs)
Q = C A V 2gH
0.025m3 /sec (0.9 cfs)
= 0.7 x 0.031 m2 (0.33 sf) x
4.43 (8.03) x V H
H = 0.21 m (0.69 ft) on center line of orifice
Depth of flow is H, plus one-half height
of orifice, or 0.26 m (0.86 ft). This is much
greater than the normal depth of flow at the
peak dry-weather flow of 0.11 m (0.36 ft)
computed previously. Therefore, the velocity
will be much less than the 0.58 m/sec (2.9
fps) computed previously and may cause
deposition of grit at peak dry-weather flow.
Determine flow to the treatment plant
when the water level is at weir crest.
Depth of flow in chamber is 1.83 m (6.0 ft)
Head on center of orifice is 1.77 m (5.83 ft)
Q = C AV2gh
= 0.7 x 0.03m2 (0.33sf) x 4.43 (8.03)
1.33m (2.41 ft)
= 0.127 m3 /sec (4.46 cfs) = 15 DWF
Hence, the flow to the plant will exceed
15 DWF during periods of design discharge if
there is no further restriction to flow
downstream of the sluice gate. One way to
restrict the flow is to design a sewer between
the sluice gate manhole and the interceptor in
such a way that it will convey the peak
dry-weather flow without surcharge but will
become surcharged when the flow exceeds the
peak dry-weather flow. This procedure is
described and illustrated by an example in
an EPA report.7
Spillway Channel
The side channel in the helical section
which conveys the overflow from the weir to
the outlet sewer leading to the stream should
be designed for the maximum flow expected
to pass through the separator. The maximum
flow will depend on the storm frequency for
which the combined sewer is designed, as well
as on the extent to which the combined sewer
can be surcharged by storm flows greater than
the design flow. It should also be assumed
that the pipe outlet to the treatment plant is
not in use either by design or by accident. On
this basis it is possible for the maximum flow
to exceed the design flowrate by 50 to 100
percent.
As an example, assume that the design
flowrate is 0.85 m3/sec (30 cfs) and the
maximum flow is 1.27 m3/sec (45 cfs). The
side channel can be designed as a lateral spillway
channel with the weir discharge spilling into it
throughout its length. To aid in self-cleaning,
it is desirable to set the downstream end of
the channel above the invert of the outlet
pipe and to provide a slope in the channel so
that at low depths of flow the velocity will
.exceed 0.31 m/sec (1.0 fps).
The channel should be designed large
enough so that the upstream water surface
will not cause submergence of the weir.
The general equation for determining the
depth of flow in a lateral spillway channel is
the following:8
h0 =
(hl -1/3U)2 -2/3H
hi
where
26
-------
h,
z
/
and
where
= upstream water depth
= critical depth
downstream water depth
when flow is submerged
= slope of channel
= length of channel
gb
- width of channel
Entrance loss =
2g
Elevation invert
outlet pipe
Elevation weir
Elevation water
at entrance
Length of channel =
length of weir, less
outlet diameter
= 0.29m (0.96ft)
= 0.00m (0.00ft)
= 1.83m (6.00ft)
1.21 m (3.96 ft)
= 16.3m (53.5ft)
The factors in the foregoing equation are
depicted in Figure 19, Spillway Channel
Profile. Actually, only 17.83/18.83, or 95
percent of the maximum flow discharges
directly into the spillway channel. In the
following example, however, it is assumed
that all the maximum flow is conveyed by the
channel.
Assume the following data:
Design fiowrate = 0.85 m3 /sec (30 cfs)
Q(side channel) = 1.27m3 /sec (45 cfs)
Outlet pipe diameter = 0.91m (3.0ft)
Weir height = 1.83 m (6.0 ft)
Then
Outlet velocity = 1.95 m/sec (6.4 fps)
Initial computation indicated that a
channel 1.83 m (6.0 ft) deep and 0.31 m (1.0
ft) wide would cause submergence of the
weir. For maintenance purposes a minimum
width of channel of 0.61 m (2.0 ft) is
considered desirable.
Preliminary computation with zero slope
and the downstream end of the channel at
elevation 0.0 indicated a water depth at the
upstream end of the channel of 1.49 m (4.9
ft).
The effect of submergence on
broad-crested weirs is surprisingly small. If
necessary the fall in the water surface over the
weir can be limited to 50 percent of the head on
weir without affecting the discharge over the
18.83D
Freeboard •
D = Transition inlet diameter
Figure assumes .that outlet
sewer has same diameter
FIGURE 19 SPILLWAY CHANNEL PROFILE
27
-------
HELICAL BEND REGULATOR
Design Example
Straight Pipe
Assume pipe designed so that at DWF the velocity is about 2 fps
D = 3.0 ft
Q full = 45 cfs
S = 0.44 %
v = 6.5 fps
n = 0.013
For Q design = 30 cfs
.Qi = 30 = o.67
Q2 45
—l = 0.6 from chart
d2
di =0.6x3.0= 1.8ft
-1 = 1.07 from chart1
V2
v, =1.07x6.5 = 7.0 fps
For Q DWF = 0.3 cfs
OIL = Qi3= 0.007
Qa 45
^=0.06
d2
£>-= 0.29
d, =0.06x3 = 0.18 ft
vi =0.29x6.5= 1.9 fps
OK almost 2.0 fps
For Q3(Peak DWF) = 0.9 cfs
QL, QJ = Q.02
Q2 45
^^ 0.1
dt =0.1 x3 = 0.3 ft
v2
= 0.4
= 0.4 x 6.5 = 2.6 fps
28
-------
Helical Bend Regulator, Design Example (continued)
Straight Pipe D = 3 ft
i = 0.3 ft
j =0.18 ft
Q = 0.9 cfs
Q = 0.3cfs
Helical Bend Regulator
0.3ft
Q = 0.9 cfs
At 3 DWF the area will only be slightly larger in separator than in 3 ft diameter
pipes. Assume velocity the same
v = 2.6 fps
Exit Pipe — Assume 1.0 ft diameter
S = 0.44 %
Q = 2.4 cfs
Qi = PJU
Q2 = 2.4
dj- = 0.42
d2
XL = 0.92
v = 3.0 fps
di = 0.42 ft
v, = 0.92 x 3 = 2.8 fps
i = 0.42 ft Q = 0.9 cfs
I
Lower invert of outlet pipe 0.12 ft below invert of separator so as not to raise water surface
Determine outlet design
when Q = 30 cfs
so that Q outlet - = 0.9 cfs
Weir Length
Angle = 60° . D = 3 ft
29
-------
Helical Bend Regulator, Design Example (continued)
Weir Radius =
Head on Weir
Weir Length =
Q per ft
Use Rehbock K
Q
JJ3/2.
H
Depth of Water
1-5/6 D
Head on weir
Total water depth
Assume short tube exit
oo
16D + 2.5D + D/3
(16) (3)+ (2.5) (3)+1
56.5 ft
60
360 V~"7 6
59ft
-3-0- = 0.51 cfs
M-
KLH3/2 = 3.41 (1)H3/2 =0.51
0.15
0.28 ft
ii- (3) = 5.50ft
/eir
>th
D
A
=
=
1.0ft
- TrD2
0.28
5.78 ft
- n 7«
4
CAV 2gH
(0.7) (0.785) (8.03)vr5T28~
10.1 cfs
Too large
0.12
Outlet Design
Determine orifice area for
Q
0.9
A
A
D
CA V 2gH
(0.7) (A) 8.03
6.9 * 12.9 =0.07
30
-------
Helical Bend Regulator, Design Example (continued)
Orifice should be greater than 0.67 ft. If orifice is made this size so that only 0.9 cfs or
3DWF would pass when tank is full, then it is apparent that the tank would fill up whenever
3DWF occurs which might be the peak daily flow.
To prevent deposition of solids on the separator floor and to prevent cleaning the separator
in dry weather periods, the separator and outlet should pass up to 3 DWF without raising
levels in separator to weir levels.
Try throttle pipe on outlet
00
Assume free fall
Outlet Design
Use 8-inch pipe as minimum
Use s
For n
For
0.4% as minimum
0.013
QDWF
QL
Q= 0.75 cfs
v = 2.2 fps'
2g
= 0.08ft
0.3 cfs
0.3
0.75
= 0.4
If discharge ratio is 0.4, then depth is 44% and velocity 94%
Assume
For
d = .44 D
v = .94 V
L
Q3DWF
v
Entrance & exit loss
Slope hydraulic gradient
Actual pipe slope
Water surface above top pipe
Water surface in tank
Depth water in tank
0.44 x 0.67
0.94 x 2.2
100ft
0.9 cfs
2.5 fps
1.5x0.1
100x0.5%
= 0.3ft
= 2.1 fps
H= 5.78-0.67
si = o.i
2g
0.15
= 0.50
0.65
0.40
= 5.11ft
s = 0.5%
0.25
JL61
0.92 ft
31
-------
Helical Bend Regulator, Design Example (continued)
Determine Q max when tank is full
Assume entrance and exit loss =
slopeH.G _-
Q max
1.6cfs
1.2ft
100+0.4 = 4.3%
v = 7.2 fps
2g
= 0.
Thus with 8 in. throttle pipe 100 ft long the maximum Q will be 1.6 cfs or about
5 times DWF.
Try other lengths of pipe
Lft 5%HG Qcfs vfps -^- Q DWF
200
300
400
600
2.45
1.9
1.45
1.2
1.8
1.6
1.4
1.3
5.4
4.7
4.0
3.8
2g
0.5
0.4
0.3
0.2
6
5.3
4.7
4.3
A length of 400 ft should be the maximum for an 8 in. sewer. Therefore it is obvious
that if the discharge is to be limited to 3 DWF some type mechanical device should be
used to close the outlet opening as the water level rises in the separator.
Determine depth of water in tank with throttle pipe 400 ft long, when Q = 0.9 cfs
Required 5 = 0.5% (See previous page)
400 x 0.5% = 2.0 ft
Slope sewer 400 x 0.4% 1.6
0.4ft
Entrance & exit loss 1.5— 0.15
2g
Head on top of pipe 0.55 ft
Diameter pipe 0.67
Depth of water in tank 1.22 ft
weir. As computed previously the head on the
weir is 0.095 m (0.32 ft) and little elevation
can be gained by assuming a submerged weir.
Therefore design can be based on no
submergence of weir. It is also desirable to
locate the downstream end of the discharge
channel above the outlet pipe invert to
prevent deposition in the channel. Therefore
the downstream end of the channel was set at
elevation 0.30 m (1.0 ft) and the channel
slope set at 0.005. The resultant freeboard of
0.18 m (0.6 ft) indicated that the channel
could have been raised an additional 0.18 m
(0.6 ft). The final data are as follows:
Elevation downstream
invert of channel 0.3 m (1.0 ft)
Rise in channel 0.09 m (0.3 ft)
Elevation upstream
invert of channel 0.40 m (1.3 ft)
h0 1.25m (4.1ft)
Elevation upstream
water surface 1.65 m (5.4 ft)
Elevation weir 1.83 m (6.0 ft)
Freeboard 0.18m (0.6 ft)
32
-------
In the foregoing example, it was assumed
that the outlet pipe was located at the end of
the channel. If the outlet pipe were located at
the center of the side channel this would
result in a lesser rise of the water surface in
the spillway channel.
Typical Dimensions
For the purpose of comparing the helical
separator with the swirl concentrator, three
design flows were selected as follows: 1.42
m3/sec (50 cfs), 2.83 m3./sec (100 cfs), and
5.66 m3/sec (200 cfs). The dimensions of the
helical separator are based on Figures 15, 16,
18, and on Profile with Level Top, in Figure
14, Effect of Transition Slope. The resultant
dimensions are tabulated in Table 2, Helical
Bend Regulator Dimensions. The dimensions
of the swirl concentrator are based on Figures
7, and 11 through 18, of Report
EPA-670/S-74-039, Relationship Between
Diameter and Height for the Design of a Swirl
Concentrator as a Combined .Sewer Overflow
Regulator. (July 1974) The resultant
dimensions - are shown in Table 6, Swirl
Concentrator Dimensions.
It should be noted that the helical
separator sized in Table 2 will remove 100
percent of grit, whereas the swirl concentrator
sized for the same discharge in Table 6 will
remove only 90 percent of grit. This is
because the published design charts for,the
swirl concentrator provide for 70, 80, and 90
percent recovery (or removal) of grit. Later in
this section, under the heading "Comparison
of Costs," a method is developed to indicate
the possible cost of a swirl concentrator to
remove 100 percent of grit.
TABLE 2
HELICAL BEND REGULATOR DIMENSIONS
Design Discharge
Inlet Diameter
Transition Length
Straight Section - Length
Radius
Width
Minimum Wall Height
Channel to Top Weir
Height End of Transition
Scum Baffle Height
Distance from Weir to
Bottom of Baffle
Weir Height
Distance Wall to Weir:
Max
Min
m3/s
(cfs)
D-m
(ft)
15D-m
(ft)
5D-m
(ft)
16D-m
(ft)
3D-m
($)
2.5D - m
(ft)
1-D - m
6 (ft)
2D-m
(ft)
D m
3 (ft)
D m
12 (ft)
D m
3 (ft)
1.42
(50)
1.07
(3.5)
16.0
(52.5)
5.33
(17.5)
17.1
(56.0)
3.20
(10.5)
2.67
(8.75)
1.95
(6.4)
2.13
(7.0)
0.36
(1.2)
0.09
(0.3)
0.36
(1.2)
2.83
(100)
1.52
(5.0)
22.9
(75.0)
7.62
(25.0)
24.4
(80.0)
4.57
(15.0)
3.81
(12.5)
2.77
(9.1)
3.05
(10.0)
0.52
(1.7)
0.13
(0.4)
0.52
(1.7)
4.67
(165)
1.83"-
(6.0)
27.4
(90.0)
9.14
(30.0)
29.3
(96.0)
5.49
(18.0)
4.57
(15.0)
3.35
(11.0)
3.66
(12.0)
0.61
(2.0)
0.15
(0.5)
0.61
(2.0)
D m
3 (ft)
D m
6 (ft)
0.36
(1.2)
0.18
(0.6)
0.52
(1.7)
0.26
(0.85)
0.61
(2.0)
0.30
(1.0)
33
-------
11.8D +30.5m (100 ft)
D2 +31.1m (102 ft)
Notes:
1.
2.
3.
4.
Swirl Concentrator
If D< 0.76 (2.5ft)
the length is 16.9 D +
30.5m (100 ft)
Or 15.2m (50 ft)
whichever is greater
D = Diameter of
Transition inlet
D2= Inside diameter
of Swirl Concentrator
FIGURE 20 SITE REQUIREMENTS
Site Requirements
The size of the buffer or protective zone
required around the helical separator will
depend to a large extent on the environment
of the neighborhood. In any locality, a buffer
zone at least 15.2 m (50 ft) wide would be
desirable. Therefore, the site requirements
given herein are based on a 15.2 m (50 ft)
buffer zone around all open or above-ground
parts of the facility. Because the transition is
below the surface, no buffer zone is required
for that part of the structure; however, it is
assumed that all of the transition section is
located on the site.
On the foregoing basis, the site
requirements for both the helical separator
and the swirl concentrator are shown in
Figure 20, Site Requirements. The required
lot dimensions and areas for three sizes of
each facility are shown in Table 3, Site
Dimensions and Areas for Helical Bend and
Swirl Concentrator Regulators. As explained
previously, the site dimensions are based on a
helical separator to remove 100 percent of
grit .and a swirl concentrator to remove 90
percent of grit.
It is evident from Table 3 that the site
requirements for the helical bend are greater
than for the swirl concentrator and that the
larger the design flow the greater the
difference. For the design flows of 1.42
m3/sec (50 cfs), 2.83 m3/sec (100 cfs) and
5.66 m3/sec (200 cfs), the helical bend
requires a site 49 percent, 106 percent, and
137 percent greater, respectively, than the
swirl concentrator.
34
-------
TABLE 3
SITE DIMENSIONS AND AREAS FOR HELICAL BEND
AND SWIRL CONCENTRATOR REGULATORS
Capacity
Relative Area
Capacity
Site Size
Site Area
Relative Area
Swirl
Concentrator
1.42m3/s(50cfs)
Site Size
Site Area
Relative Area
Capacity
Site Size
Site Area
38.0 mx 38.0m
(124.5 ft x 124.5 ft)
1,440m?
(15,500sf)
1.00
2.83m3/s(100cfs)
40 m x 40 m
(131.5 ft x 131.5 ft)
1,600m2
(17,300 sf)
1.00
4.67m3/s(165cfs)
42 m x 42 m
(138 ft x 138 ft)
1,770m2
(19,000 sf)
1.00
Helical
Separator
43.0m x 54.6 m
(141 ft x 179 ft)
2,340 m2
(25,200 sf)
1.63
48.5 mx 71.5m
(159 ft x 234 ft)
3,460 m2
(37,200 sf)
2.15
52.0 m x 82.8 m
(171 ft x 272 ft)
4,300 m2
(46,500 sf)
2.45
Comparison of Hydraulic Head
Losses in Helical Bend and
Swirl Concentrator Regulator
The available head at a specific site may
be a critical factor in the choice of the
specific type of combined sewer regulator to
be used. The head loss must be considered for
two - conditions: (1) For periods of
dry-weather flow, and (2) for periods of
wet-weather flow. The available head during
dry-weather flow will depend on the
difference in elevation between the combined
sewer and the interceptor that will convey the
flow to the wastewater treatment plant. The
available head during wet-weather flow will
depend on the difference in elevation between
the combined sewer and the water surface of
the receiving stream. A further consideration
in the latter case is whether the existing
combined sewer is to be used to convey the
overflow from the regulator to the receiving
stream or to any holding or treatment
facilities involved.
In periods of dry weather there will be
open channel flow through the regulator, and
hence the drop in invert and hydraulic
gradient will be similar from the regulator
inlet to the foul outlet. In wet-weather
periods there may be a considerable
difference between the drop in invert and the
drop in hydraulic gradient from the regulator
inlet to the clear outlet.
The usual practice is to design combined
sewers to run full at design discharge, i.e., the
hydraulic gradient is at the top of the pipe. In
special -cases, due to topography or for other
reasons, it may be possible to design the sewer
to be surcharged, i.e., the hydraulic gradient is
above the top of the sewer. Both cases are
considered below. To simplify the discussion,
the slope of the channel in the regulator to
maintain a dry-weather flow velocity that will
prevent solids deposition has not been
included in the required drop in the invert.
In the following comparisons, the head
losses are given as a multiple of the inlet
35
-------
dimension: D, the inlet diameter of the
helical separator; and Dt the side of the
square inlet of the swirl concentrator. To
show that D and DI are approximately the
same for the same discharge, their values for
three discharges are given, as follows:
Design
Discharge m3/sec 1.42 2.83 5.66
(cfs) (50) (100) (200)
D
m
(ft)
m
(ft)
1.07
(3.5)
0.90
(3.0)
1.52 1.98
(5.0) (6.5)
1.52 1.83
(5.0) (6.0)
First, consider the case where there is to
be no surcharge on the inlet during design
discharge.
In the helical separator the transition will
have a level top as shown in Figure 14. The
drop in the invert of the transition will be 1.0
D. Therefore, the drop in the invert from the
inlet to the foul outlet will be 1.0 D
(neglecting the slope of the channel through
the regulator). The loss in the hydraulic
gradient will also be 1. O D. The invert of the
clear outlet will be at approximately the same
elevation as the invert of the separator, as
explained previously in the discussion of the
weir overflow spillway channel. Therefore,
the drop in the invert between the inlet and
the clear outlet will also be 1.0 D. The loss in
the hydraulic gradient may be the same as the
drop in the invert or it may be slightly
different depending on outlet design.
The dimensions of the swirl concentrator
as a Combined Sewer Overflow Regulator are
given in Figure 11, General Design Details, of
Report EPA-670/2-74-039, Previously cited.
The general construction details are shown in
Figures 12, 13, and 14 of EPA Report
R2-72-008, titled "The Swirl Concentrator as
a Combined Sewer Overflow Regulator
Facility." (Sept. 1972)
If there is to be no surcharge on the inlet
sewer, the crown must be at a distance above
the invert of the chamber equal to H! (the
height of weir above the chamber invert), plus
the head on the weir. The drop in the invert
of the sewer will be this distance less D! , the
dimension of the inlet. The foul outlet is
located below the chamber bottom. Assuming
a foul outlet diameter of 0.31 m (1.0 ft) and
concrete cover over the outlet to the same
amount, the distance from the chamber invert
to the outlet invert is 0.61 m (2 ft).
Excluding the channel slope through the
concentrator, the drop in the invert from the
inlet to the foul outlet is, therefore, 0.8 Dt to
1.5 DI, plus 0.6 m (2.0 ft). The foul outlet
pipe diameter may exceed 0.31 m (1.0 ft)
diameter for larger flows, thus increasing the
total drop somewhat. The hydraulic gradient
will have a similar drop.
The clear outlet is also located below the
chamber floor and, if a 0.31 m (1.0 ft)
concrete cover is provided over the outlet, the
vertical distance from the chamber invert to
the invert of the clear outlet will be 1.0 Dt,
plus 0.31 m (1.0 ft). Combining this with the
entrance drop of 0.8 D! to 1.5 D! , will result
in a total drop in the invert from the inlet to
the clear outlet of 1.8 Dt to 2.5 Dt, plus 0.31
m (1.0 ft). The drop in the hydraulic gradient
in this case will be different. The circular weir
is set a distance equal to the head on the weir
below the top of the inlet sewer. If there is no
submergence of the weir then the loss in the
hydraulic gradient will be equal to this head.
Trial computations indicate the head on the
weir is about 0.2 D! . Allowing for friction
losses in the outlet pipe and some freeboard
downstream of the weir, the drop in hydraulic
gradient is about 0.4 D! .
When the inlet sewer is surcharged
different hydraulic conditions will exist in the
helical separator. If the inlet is surcharged an
amount equal to D, the transition invert will
be level, as shown in Figure 14. The drop in
the invert from the inlet to the foul outlet
will be zero (neglecting the channel slope
through the separator). The drop in hydraulic
gradient will also be zero. Likewise, the drop
in the invert from the inlet to the clear outlet
will be zero. However, the loss in hydraulic
gradient for this case will be 1.0 D.
In the case of the swirl concentrator, if a
surcharge of Dt is permitted then the crown
of the sewer can be set a distance of D, below
the water surface of the chamber. The drop
from the chamber invert to the foul outlet
will be 0.61 m (2.0 ft), as previously
computed. Therefore, the drop in the invert
from the inlet to the foul outlet will be 0 to
0.5 D1; plus 0.61 m (2.0 ft). The drop in
36
-------
TABLE 4
TYPICAL HEAD LOSSES IN HELICAL BEND
AND SWIRL CONCENTRATOR REGULATORS
Helical Swirl
Separator Concentrator
Dry Weather Flow — Drop in Invert
Helical Separator
Transition invert level
Transition roof level
Swirl Concentrator
Wet Weather Flow
Helical Separator
Transition invert level
Hydraulic grade
Drop in invert
Transition roof level
Hydraulic grade
Drop in invert
Swirl Concentrator
Hydraulic grade
Drop in invert
Note: Friction losses not included in above Table.
hydraulic gradient will be the same. The drop
from the chamber invert to the clear outlet
invert will be 1.0 Dt, plus 0.31 m (1.0 ft), as
before. Therefore, the total drop from the
inlet invert to the clear outlet invert will be
1.0 D, to 1.5 D, plus 0.31 m (1.0 ft). The
drop in hydraulic gradient will be about 0-4
DI as computed previously.
The data relative to the foregoing
discussion are shown in Table 4, Typical Head
Losses in Helical Bend and Swirl Concentrator
Regulators.
From the Table, it is apparent that the
drop in the invert is always greater in the swirl
concentrator than in the helical separator.
When the inlet sewer is not surcharged, the
drop to the foul outlet is only slightly greater
but the drop to the clear outlet is about twice
as great. When the inlet sewer is surcharged an
amount equal to D or D! , the drop in the
invert is zero in the helical separator,
compared to the minimum drop of 0.61 m
(2.0 ft) to the foul outlet and of 1.0 Dt, plus
0.31 m (1.0 ft) to the clear outlet in the swirl
concentrator. For dry-weather flows, the drop
None
1.0 D
1.08 D to 1.33 D
1.0 D
None
1.0 D
1.0 D
0.58 D
2.1 D
in hydraulic gradient is similar to the drop in
the invert. For wet-weather flows, the drop in
hydraulic gradient from the inlet to the clear
outlet in the-swirl concentrator is about
one-half that for the helical separator.
Construction Details
Means of access must be provided to the
curved section of the separator for
maintenance purposes, including possible
washing down after each storm event. The
provision of a superstructure over this section
is desirable for safety and esthetic reasons,
and for confining possible odors. The type of
superstructure used will depend on the
character of the locality. As a minimum and
for purposes of this report, the walls are
assumed to be of concrete block and the roof
of precast concrete units. For roof spans
exceeding about 8.5 m (28 ft), it will be
necessary to provide structural steel framing.
A cross section of the helical separator
with a superstructure is shown in Figure 21,
Typical Cross Section, Helical Bend
Regulator. In this Figure, for cost estimating
37
-------
s?
£
r*
E
(0
CO
CO
Q
CM
s/Stf/
1 Precast Concrete Roof
• J
!v~V
—
r
^_
*'•'
• \
Vr
H
*.<•
.#
•'''i.
•'.77
-'4
«'
**t
;-iS.-4:<' '.'A:'.\':-'i'-A.'
^ Concrete block walls -^_^_^
Stainless steel rail — — »
>^ \
f 4
4' Concrete walk jl 1
| \
Concrete beam __ || \
7 II \
|WiSf»»jsSjMS>»{*Sk:{
Fiber glass baffle ^
r5fl
/^°iJ r
Flushing pipe ^^ ;
(i
*- 'i'?'f
•f t ^t
"<5 ^^r'N.a . i/;
CM ^^sa--"' " '••
• ^f^: «
V 5 ^^' .4:?
:*fi;\ 0 ^<"->i' "'<'--
V«t\ ^f-;- • Y • i-..
,'v-r ;K^ ,___^^ ' ^
_j,
*. —
'•'.V
K15
't.(
f;'A'
"V* '
;4f
M
1 ..«
;-'f^
y/
J " •
\t
**> t * . " <4'« »• *^*« • »rf^ i •• •*•&*/
' ' • — r^ '— '
'\
1.3D
wrT . „ te.
-^ Approximately 4.1 D
>,»
i ,
* *
oo c;
E E
•* «*?
•* o
CM c-
I
^\--4^\V. ^
<_Min.0.3m. (1 ft)
Spillway
Channel
0.67 D
_L
A
f
u
T
t:
CM
£
r—
CO
FIGURE 21 TYPICAL CROSS SECTION - HELICAL BEND REGULATOR
purposes the width of structure is indicated as
4.1 D and the width of spillway channel as
0.67 D. For any specific case the width and
elevation of the spillway channel may vary
from that shown in Figure 21, Typical
Cross-Sectional Helical Separator, as
explained previously.
The hydraulic conditions require that the
transition section be provided with a roof.
These conditions do not apply to the straight
section, having a length of 5 D, preceding the
curved section. It is not apparent from the
model studies that this section will require the
same maintenance as the curved section.
Accordingly, there appears to be no need to
make this section accessible.
For purposes of this report it has been
assumed that the straight section will have
walls 2.5 D high and will be provided with a
concrete roof at that elevation.
Other construction details considered
necessary or desirable are as follows:
38
-------
a. Provide concrete walls with a minimum
thickness of 0.3 m (1 ft), extending above
grade a minimum height of 0.3 m (1 ft).
b. Provide a concrete walk 1.2 m (4.0 ft)
wide.
c. Provide a stainless steel railing on each
side of the walk.
d. Provide a fiberglass scum baffle hung
from the beams or supported from the
weir.
e. Provide a flushing water pipe on the
channel side of the scum baffle and hung
from the beams. Connect this line to the
public supply with a backflow device, if
this is permitted by local code. If this is
not permissible, provide a storage tank to
store overflow from the weir and -a
submersible pump to use for washing
down. The usual criteria of 3.1 1/s (50
gpm) at 28,120 N/m2 (40 psi) for
flushing purposes at treatment plants
should be applicable to the helical
separator facility. Hose connections
should also be provided in case the stream
from the wash water pipe is not effective.
f. Provide concrete block walls with a
height of 2.4 m (8.0 ft).
g. Provide a precast concrete roof.
h. Provide adequate electric lights.
i. Provide roof ventilators,
j. Provide doors at both ends of the
structure for ventilation and access.
Cost estimates of the helical separator
were made for two purposes: (1) to indicate
the probable construction cost of the facility;
and (2) to compare its costs with that of the
swirl concentrator used as a combined sewer
regulator.
The cost estimates are considered to be
reasonable engineers' estimates. However,
during periods of economic inflation, it is not
unusual for contractors' bids to materially
exceed engineers' estimates.
In making a choice between the helical
separator and the swirl concentrator, it is
possible that other factors related to the
specific site of the facility may be of greater
importance than the difference in
construction costs.
Quantity Basis
The estimated quantities are based on the
following:
a. The transition will be constructed with a
drop in the invert equal to D as shown in
Figure 14, so that the sewer upstream of
the transition will not be surcharged.
b. The straight section preceding the curved
section will have walls 2.5 D high and
concrete roof.
c. The superstructure over the curved
section will be as shown in Figure 21.
d. The width of the curved section is
assumed to be 4.1 D; the width of the
spillway channel is assumed to be 0.67 D.
e. The .cover on the sewer at the transition
inlet will be 2.44 m (8.0 ft).
f. The ground is level and the subsurface is
earth with no groundwater problems. -
g. All concrete walls will have a minimum
thickness of 0.3 m (1.0 ft), except the
weir.
h. Sheet piling will be required about 0.6 m
(2.0 ft) outside the structure.
i, Transverse concrete beams will be
required at 4.5 m (15 ft) intervals with
the cross section 0.45 m (1.5 ft) square.
j. The continuous concrete walk will be
1.22 m (4.0 ft) wide and 0.20 m (0.67 ft)
thick.
Cost of Helical Separator
The costs are based on the following:
a. The Engineering News-Record
Construction Cost Index average for the
United States is 2,100
b. Unit prices are as follows:
Steel Sheer Piling $ 86/m2 $ 8/sf'
Excavation $ 16/m3 $ 12/cy
Reinforced Concrete $326/m3 $250/cy
Concrete Block Walls $107/m2 $ 10/sf
Roof $150/m2 $ 14/sf
c. Miscellaneous costs are assumed to be 25
percent of the foregoing items and to
include a manual sluice gate and manhole,
handrail, flushing water facilities, scum
baffle, electrical work, roof ventilators
and doors.
d. The estimated cost of the bypass sewer
during construction is based on providing
a sewer of the same diameter as a
transition inlet around the proposed
separator, plus an allowance for
temporary connections at each end.
e. Contingent and engineering costs will be
25 percent of the foregoing items.
39
-------
The resultant costs are shown in Table 5,
Construction Cost of Helical Bend Regulators.
Cost of Swirl Concentrator
Typical dimensions for three sizes of the
swirl concentrator are given in Table 6, Swirl
Concentrator Dimensions, based on Figures 7
and 11 of Report EPA-670/2-74-039. The
cost estimates are based on these dimensions
and on the construction details shown in
Figures 12, 13, and 14 of Report
EPA-R2-7 2-008.
TABLE 5
CONSTRUCTION COST OF HELICAL BEND REGULATOR
Capacity 1.42 m3 /s (50 cfs)
ITEM
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Miscellaneous Costs
Bypass Sewer
Contingent & Engineering Costs
Capacity 2.83 m3/s dOO cfs)
ITEM
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Miscellaneous Costs
Bypass Sewer
Contingent & Engineering Costs
Capacity 4.67 m3 /s ( 1 65 cfs)
ITEM
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Miscellaneous Costs
Bypass Sewer
Contingent & Engineering Costs
420m2
950m3
250m3
114m2
85 m2
QUANTITY
(4,550 sf)
(l,240cy)
( 330 cy)
(1,230 sf)
( 910 sf)
Sub Total
25%±
Sub Total
Total
QUANTITY
710 m2 (7,700 sf)
2,200m3 (2,800 cy)
475m3 ( 620 cy)
160m2 (l,740sf)
170m2 (l,800sf)
Sub Total
25%±
Sub Total
25%±
Total
QUANTITY
950m2 (10,200 sf)
3,200 m3 (4,180 cy)
679 m3 ( 888 cy)
200m2 (2,130sf)
250 m2 (2,700 sf)
Sub Total
25%±
Sub Total
25%±
Total
AMOUNT
$ 36,150
15,200
81,500
12,200
12.750
$157,800
39,450
20.000
$217,250
54.000
$271,250
AMOUNT
$ 61,050
35,200
154,850
17,150
25.500
$293,700
73,300
39.000
$406,000
102.000
$508,000
AMOUNT
$ 81,700
51,100
221,350
21,400
37.500
$413,150
103,850
56.000
$573,000
143.000
$716,000
40
-------
TABLE 6
SWIRL CONCENTRATOR DIMENSIONS
Design Discharge
Diameter of Chamber
Diameter at Inlet
Height No. 1
1-5/6 D + 2.74 m
(1-5/6 D 2 2.74m
Height No. 2
Floor to Weir dr
Head on Weir for 150%
Design Discharge
Clearance to Walk
Headroom
Total
Use for Height
Note: Height is from invert to
' • m3/3
(cfs)
D2 m
(ft)
D! m
(ft)
m
(ft)
m
(ft)
m
(ft)
m
(ft)
m
(ft)
m
(ft)
m
(ft)
underside of roof.
1.42
(50)
6.85
(22.5)
1.15
(3.75)
4.55
(14.9)
1.71
(5.6)
0.27
(0.9)
0.30
(1.0)
2.44
(8.0)
4.72
(15.5)
4.72
(15.5)
2.83
(100)
9.0
(29.5)
1.50
(4-92)
5.18
(17.0)
2.26
(7.4)
0.40
(1.3)
0.30
(1.0)
2.44
(8.0)
5.40
(17.7)
5.40
(17.7)
4.67
(165)
11.0
(36.0)
1.83
(6.00)
5.79'
,(19.0)
2'. 74
(9.0)
0.55
(1.8)
0.30
(1.0)
2.44
(8.0)
6.03
(19.8)
6.03
(19.8)
The height of the structure from the
invert of the chamber to the underside of the
roof is based on the following criteria:
Criteria 1: The clearance between the
top of the walk and the water surface is 0.31
m (1.0 ft) when the discharge is 150 percent
of design discharge and the foul outlet is not
functioning. The head on the weir is
determined from Figure 11, Head Discharge
Curve for Circular Weir, in Report
EP A-R 2-72-008. It is assumed that the
emergency weir starts overflowing when
inflow to the chamber reaches the maximum
design discharge for the facility.
Criteria 2: The headroom: above the walk
is 2.44 m (,8.0 ft). The depth of the structure
below the ground surface was based on the
following criteria:
Criteria 3: The cover on the crown of the ..
sewer is 2.44 m (8,0 ft).
Criteria. 4: To prevent surcharge on the
inlet sewer, the crown of the inlet sewer shall
be the same elevation as the water surface at
design discharge. :
As shown in Table 6, the depth below the
ground surface, based on Criteria 3 and 4, is
less than the interior height, based on Criteria
1 and 2; therefore, the height is adequate to
provide a structure with a roof above the
ground surface. If the roof slab is assumed to
be 0.25 m (0.83 ft) thick, then the top of
roof is about 0.61 m (2 ft) above the ground
surface, as shown in Table 6.
Additional assumptions for estimation
purposes aie as follows:
a. The walls are 0.30 m (1.0 ft) thick.
b. The roof is of poured concrete, about
0.25 m (0.83 ft) thick, with two beams
0.92 m (3.0 ft) by 0.46 m (1.5 ft).
c. The bottom concrete slab is 0.61 m (2.0
ft) thick.
d. The concrete walk is 1.22 m-(4.0. ft) wide
and supported on concrete beams.
e. The superstructure is 3.96 m (13.0 ft)
long by 1.52 m (5.0 ft) wide, by 2.44 m
(8.0 ft) high.
The unit prices used were the same as
those used for the helical separator.
The miscellaneous cost is taken as 25
percent and is intended to include stairs,
handrail, scum baffle, circular weir, flushing
water system and pipes, a manual sluice gate
and manhole, electrical work, ventilating
work and doors.
41
-------
The cost of the bypass pipe is assumed to
be 50 percent of the bypass pipe for the
helical separator.
Contingent and engineering costs are
taken as 25 percent of the foregoing.
The costs of the three selected sizes of
swirl concentrators, designed for 90 percent
removal of grit, are shown in Table 7,
Construction Cost — Swirl Concentrator.
Comparison of Costs
The estimated costs are shown graphically
in Figure 22, Estimated Construction Costs —
Helical Bend and Swirl Concentrator
Regulator. As explained previously, although
the two structures are sized for the same
TABLE 7
CONSTRUCTION COST - SWIRL CONCENTRATOR
Capacity
1.42m3/s fSOcfs)
Item
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Outlet Pipes
Downshaft and Plate
Miscellaneous Costs
Bypass Sewer
Contingent & Engineering Costs 25%±
Capacity 2.83 m3 Is (IQOcfs)
Item
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Outlet Pipes.
Downshaft and Plate
Miscellaneous Costs
Bypass Sewer
Contingent and Engineering Costs 25%±
Capacity
4.67 m3/s <165cfs>
Item
Sheet Piling
Excavation
Reinforced Concrete
Concrete Block Walls
Roof
Outlet Pipes
Downshaft and Plate
Miscellaneous Costs
Bypass Sewer
Contingent & Engineering Costs
Quantity
200m2
(2,160 sf)
460 m3
(600 cy)
98m3
(128cy)
27m2
290 sf)
6m2
(65 sf)
Sub Total
25%±
Sub Total
Total
Quantity
290 m2
(3,120 sf)
900m3
(l,180cy)
156m3
(204 cy)
27m2
(290 sf)
6m2
(65 sf)
Sub Total
25%±
Sub Total
Total
Quantity
375m2
(4,030 sf)
1,360 m3
(1.780cy)
216m3
(282 cy)
27m2
(290 sf)
6m2
(65 sf)
Sub Total
2S7r±
Sub Total
25#±
Total
Amount
S 17,200
7,360
31,948
2,889
900
1,300
2.000
S 63,597
15,403
10.000
S 89,000
22.000
5111,000
Amount
S 24,940
14,400
50,856
2,889
900
3,000
3.000
S 99,985
24,515
19.500
$144,000
36.000
5180,000
Amount
S 32,250
21,760
70,416
2,889
900
4,000
5.000
5137,215
34.785
28.000
5200,000
50.000
5250,000
42
-------
10-
<=>.
§
Notes:
A
d
o
Helical Separator
100% Grit Removal
Swirl Concentrator
100% Grit Removal
Swirl Concentrator
90% Grit Removal
Helical Separator
100% Grit Removal
Swirl Concentrator
100% Grit Removal
Swirl Concentrator
90% Grit Removal
50 100
150
200
Discharge — CFS
FIGURE 22 ESTIMATED CONSTRUCTION COSTS - HELICAL BEND
AND SWIRL CONCENTRATOR REGULATOR
discharge, the helical separator will remove
100 percent of the grit compared to 90
percent for the swirl concentrator. The use of
90 percent removal is based on Report
EPA-670/2-74-039, dated July, 1974, in
which design curves are presented for 90, 80,
and 70 percent recovery (or removal) of grit.
An earlier report on the swirl concentrator,
Report EPA-R2-72-008 of September, 1972,
includes Figure 22, Separation Efficiency
Curve, in Appendix 1 thereof. This curve
indicates that a swirl concentrator designed
for 90 percent removal of grit will remove
100 percent of grit when the discharge is 60
percent of the design discharge. Therefore, to
remove 100 percent of the grit the structure
should be sized for a discharge equal to 167
percent of the design discharge. Using this
method, the costs of swirl concentrators to
remove 100 percent of grit were estimated
and the results shown in Figure 22.
The summary of all the estimated costs is
shown in Table 8, Comparison of
Construction Costs — Helical Bend and Swirl
43
-------
TABLE 8
COMPARISON OF CONSTRUCTION COSTS - HELICAL
BEND AND SWIRL CONCENTRATOR REGULATORS
Capacity
1.42m3/s(50cfs)
2.83m3/s(100cfs)
4.67m3/s(165cfs)
Note: Land costs not included.
Swirl
Concentrator
$111,000
180,000
250,000
Concentrator Regulators. The results indicate
that, for the sizes used in the Table, the
helical separator will cost 1.7 to 2.5 times as
much as a swirl concentrator designed for 90
percent grit recovery and 1.5 to 1.6 times as
much as one designed for 100 percent grit
recovery.
Alternate Design — Solids Removal Efficiency
Further design techniques and the
rationale based on mathematical data were
evolved as a result of the model studies
carried out by the General Electric Company.
The procedure for designing a helical flow
separator for a given discharge flowrate and
specific suspended material is as follows:
Step 1: Characterize the settling
properties of the suspended solids. A study9
(see reference section) has emphasized the
variability of the settling properties at
different geographical locations. Once the
cumulative distribution of settling velocities
has been determined from column settling
tests, the distribution should be broken down
into 10 or more fractions, and the
performance calculations outlined below
should be repeated for each settling velocity
fraction. Thus, assume that 10 percent of the
material settles at a velocity lower than w2
etc. The performance for each settling
velocity is then determined and the overall
recovery efficiency is computed as the
weighted average of the recovery for the
individual fractions.
Helical
Separator
$271,000
508,000
716,000
Step 2: Determine a scale factor, S, as
the ratio
in which Q^ is the design discharge rate and
Q, is a selected flowrate within the range of
the laboratory experiments 0 to 0.85m3 /sec
(0-31.1 cfs). Selecting Q! near the lower end
of the range will give better recovery
efficiencies, but will also require a larger unit.
The selection of Q! can be iterated, as
described in the subsequent step.
Step 3: Determine a rough estimate of
the prototype performance by calculating
first an equivalent settling velocity we from
= w
.5
where ws is the median (50%) settling
velocity from the cumulative distribution
determined in Step 1. The initial estimate of
efficiency can then be found by entering
Figure 23, Predicted Separation Efficiency
Versus Settling Velocity at Several Flowrates,
with settling velocity we, and interpolating
for the recovery efficiency at the selected
flowrate QA. If this estimated efficiency is
lower than about 80 percent, it is probably
desirable to select a lower value of Q{ (to
provide a larger unit), and repeat steps 2 and
3. Conversely, if the estimated efficiency is
greater than 95 percent, a smaller unit will
probably suffice, and a larger Q, can be
selected.
44
-------
Step 4: A refined estimate of the
recovery efficiency can be obtained by
calculating an equivalent settling velocity for
each settling velocity fraction w,- of step 1 .
w
The recovery efficiency for each of these
fractions is then determined from Figure 23
as in Step 3. The overall recovery efficiency is
the weighted average of the results for each
individual fraction.
Step 5: The inlet diameter, D, should ,
then be
D = 0.95 s (ft)
or
D = 3 S (meters)
The inlet dimensions are given in terms of
D as described in the preceding paragraphs.
100
95
s?
8
0)
cc
85
80
Qi
= 0.1**
I III
0.1
1.0
Settling Velocity, w (cm/sec)
FIGURE 23 PREDICTED SEPARATION EFFICIENCY VERSUS
SETTLING VELOCITY AT SEVERAL FLOWRATES
45
-------
SECTION IV
IMPLEMENTATION
Use of the helical bend combined sewer
overflow regulator should be considered
where there is a requirement to improve the
quality of the combined sewer overflow, or
where there is insufficient hydraulic head to
allow construction of a swirl concentrator.
Table 4 indicated the typical hydraulic head
losses in both the helical bend and swirl
concentrator regulators. The helical bend
regulator, because of its length, can be expected
to have a higher construction cost than the
swirl concentrator. Removal of floatable solids
should be of the same order.
An important difference in the operation
of the two units, however, is that the solids
removed from the overflow by the helical
bend regulator may be expected to be
released to the interceptor sewer at the end of
the storm — the first flush will be at the end
of the storm. This characteristic may result in
less shock loading and more efficient
operation at the treatment facility.
The helical bend regulator is essentially
linear and, thus, construction should be
possible within existing street rights-of-way.
The essential features of the device are as
follows:
Inlet: The inlet is the collector sewer of
any shape.
Transition Section: The helical bend can
curve to the right or left, depending on the
direction of flow in the interceptor sewer or
the available right-of-way. The transition
section flares from the inlet to a width of 3D
in a distance of 15D. The section also deepens
to a depth of 2.0D. The section is covered to
facilitate the development of uniform flow
conditions.
Straight Sectibn: A 5D straight section
with the same configuration as the remainder
of the structure is required to fully control
the flow and eliminate cross currents.
Bend: A 60 degree bend of a radius is
selected on the basis of the design used. The
floor of the bend should be smooth.
Weir: A side-spill, curved weir along the
outside of the bend is used to draw off the
clarified overflow.
Scumboard: A scumboard is set 1 /6D to
1/3D in from the weir, with a minimum
submergence depth, to trap floatables.
Outlet: The outlet should be sized to
allow 4.5 to 5 times the dry-weather flow to
pass without activating the regulation
characteristics. This will allow passage of peak
dry-weather flows. However, maximum flow
to the wastewater treatment plant during
storm events is limited to 1 to 3 times its
dry-weather flow. Thus, a mechanical gate
controlled by the flow in the interceptor is
desirable. The outlet should be either 30.4 cm
(12 in.) or of a size capable of being cleaned
very easily in order that rags and other debris
will not clog the opening.
The helical bend regulator can be
constructed of concrete or any other
appropriate material. Ordinarily, it will be
desirable to construct the regulator within a
chamber in order that maintenance may be
readily performed.
Normal requirements for such chambers,
such as ventilation, lighting, safe access and
other features must be provided.
46
-------
SECTION V
Glossary of Pertinent Terms
Helical Bend Report
Combined Sewer — A pipe or conduit which
collects and transports sanitary sewage, with
its component commercial and industrial
wastes and inflow and infiltration during
dry-weather conditions, and which, in
addition, serves as the collector and conveyor
of stormwater runoff flows from streets and
other sources during precipitation and thaw
periods; thus, a pipe or conduit which handles
all of these waste waters in a single or
"combined" facility
Dip Plate — A vertical plate or baffle which is
partially immersed in flowing liquid in a
manner that will prevent the discharge of
surface or floating materials over an outlet
weir; in the helical bend studies, a baffle plate
placed near the overflow weir, at a
predetermined distance therefrom and with a
predetermined depth of "dip" or immersion,
to prevent the overflow of floating solids with
the clarified effluent
DWF (Dry-Weather Flow) - The flow
through the helical bend during periods when
no stormwater runoff is collected and
transported by the combined sewer which
discharges through the helical bend section
Equation of Motion - A mathematical
equation which expresses the hydraulic
patterns of flow occurring in any vessel,
chamber, conduit or other appurtenant
structure in which fluid flow is occurring
Foul Sewer - The sewer connected to the
foul outlet of the helical bend, through which
the solids slurry deposited in the bottom of
the helical section is discharged to a
downstream wastewater treatment facility or
to any other predetermined point of disposal
Floatables - The lighter-than-water solids,
including congealed floating materials which
rise to the surface of the wastewater flowing
through the helical bend regulator, and which
must be intercepted and prevented from
discharging with the clarified effluent flowing
over the outlet weir; the material retained by
a dip plate or scum board
Floor Flow Angle - The angular pattern of
the flow of wastewater and solids as it
traverses the bottom or floor of the helical
bend separator as part of the helix form of
liquid movement through the device; in this
study, the floor flow angles were disclosed by
the pattern of movement along the bend floor
of Process Black ink injected at appropriate
points along the floor.
Gilsonite® - Synthetic solid material utilized
in the helical bend studies to simulate the grit
particles which will be contained in combined
sewer flows to be handled in prototype units
in actual field practice; material having a
specific gravity of 1.06 and a size range of 0.5
to 3 mm
Grit — Heavier and larger-sized solids
contained in stormwater flows in combined
sewers, or in dry-weather flows, which
because of their size and specific gravity settle
in a helical bend separator more readily than
organic materials contained in the flow;
material having a specific gravity of 2.65,
more or less, and an effective size of 0.2 mm.
Helical Bend - A physical configuration of a
pipe or open channel which results in a bend
or radius through which a liquid flow occurs
in a manner that produces helical, or
secondary flow phenomena, inducing the
rapid separation of solids from the liquid and
the deposition of the solids along the inner
diameter of the radius; in the study, the total
helical bend assembly consisted of a transition
section, a straight section, and the bend
section
Helical Flow — The pattern of liquid flow
induced by the helical bend, characterized by
a helical configuration, or secondary motion,
created in the liquid flow
Hydraulic Head Loss — The lowering of the
hydraulic flow line through a pipeline, device,
chamber or other facility, due to dynamic
conditions which produce friction, turbulence
or other conditions that are translated into
loss of pressure, or head, or free water
gradient surface level
47
-------
Mass Flux — The rate of flow of the total
body of liquid and solids flowing through any
device or facility, such as a helical bend; the
flow energy in the total fluid body in the
device or facility; a measure of the strength of
a field of force in a given area, such as in any
part of a helical bend section, in the case of
the study reported herein
Organic Solids — Solids of a non-grit or
light-weight nature, contained in a sewer flow,
which are subject to decomposition and
consequent oxygen-consumption in receiving
waters into which stormwater overflows from
combined sewers discharge, and which impose
pollutional oxygen-demand loadings on such
water resources; the function of a helical flow
regulator is to control the volume of overflow
wastewaters and to remove grit and organic
solids from the effluent.
Petrothene® — A synthetic plastic material
which has a specific gravity of 1.01 and was
used in the helical separator studies to
simulate the lighter organic materials
contained in combined wastes handled by
prototype helical bend regulators in actual
field practice
Prototype — A full-scale replica of the
laboratory test model (in this case, of the
helical bend model), designed in scaled-up
proportion to the test model, for use in actual
field installations; the study provided design
criteria for the dimensioning of various sizes
of full-scale installations to handle
predetermined volumes of combined sewer
flows
Recovery — The percentage of solids
introduced into the helical bend
separator-regulator with the combined
wastewater flow that will settle in the bend
chamber and be drawn off from the bottom
via the foul outlet and the foul sewer; the
measure of the efficiency of the helical bend
device hi clarifying the overflow effluent and
thereby reducing the pollutional effect of the
overflow incidents on receiving waters
Regulator — Any device in a combined sewer
system which regulates or controls the
amount of concentrated wastewater diverted
to interceptor sewers and downstream
treatment facilities and, consequently, the
amounts of effluent discharged to receiving
waters, or holding or treatment facilities
Scum Board — A vertical plate or baffle which
is partially immersed in the flow liquid and
partially above the flow line, in order to
capture and retain floating scum solids and
prevent their discharge over an overflow weir
with the clarified effluent; in this study, a
scum board and a dip plate have similar
functions
Solids Simulation — The use of test solids in
the model, of a weight and size which
simulate, in scale-down, the composition and
character of solids materials in actual
combined sewer flows to be handled in
prototype devices — in this case, in prototype
helical bend systems in combined sewer lines
Straight Section — The part of the helical
bend structure which precedes the bend
section and delivers the flow uniformly and
without velocity interferences into the helical
section; in the studies of the helical bend
principle, it was determined that the straight
section having a length of five times the
diameter of the sewer pipe will be required
for effective solids recovery in the helical
bend system
Spillway Channel — The channel or conduit
which receives the overflow effluent from the
helical bend weir section and delivers it to a
pipe or conduit leading to receiving waters, or
facilities for the retention and/or treatment of
the clarified wastewater discharge
Swirl Concentrator — A cylindrical chamber
in which hydraulic flows experience
swirl-type liquid flow patterns that induce
relatively rapid separation of solids from the
flow, and, thereby, produce a clarified
effluent which is discharged from the
chamber, freed from solids which are
collected in the bottom of the swirl chamber
and discharged via a foul sewer outlet; in this
study, the helical separator was intended to
perform the same general function as a swirl
concentrator used as a combined sewer
regulator
Transition Section — That portion of the
helical bend composite which carries the
48
-------
combined sewer flow from the entering sewer
pipe section and delivers it to the straight
section and thence to the bend section; the
transition section in the studies had a length
of at. least fifteen times the inlet sewer
diameter and expanded the flow cross section
to three times the inlet diameter
WWF (Wet-Weather Flow) - The flow in the
combined sewer during periods of
precipitation or thaw runoff, composed of the
dry-weather flow plus the storm runoff
volume
The following are registered trademarks:
Gilsonite
Petrothene
Perspex
Plexiglas
49
-------
SECTION VI
REFERENCES
1. Problems of Combined Sewer Overflow
FWPCA - WP-20-11, December, 1967.
2. Swirl Concentrator as a Combined
Sewer Overflow Regulator Facility, EPA —
R2-72-008, September 1972. PB-214 687.
3. The Secondary Flow in a Meandering
Channel, by T.M. Prus-Chacinski,
Unpublished, PhD Thesis, University of
London, 1955.
4. Secondary Motion Applied to Storm
Sewage Overflows, by T.M. Prus-Chacinski
and J. W. Wielogorski, Symposium on Storm
Sewage Overflows, Institute of Civil
Engineers, London, 1967.
5. Formulas for Bed-Load Transport by E.
Meyer-Peter and R. Muller, Proceedings IAHR
Second Meeting, Stockholm, 1948.
6. The Bed-Load Function for Sediment
Transportation in Open Channel Flows by H.
A. Einstein, U.S. Department of Agriculture,
Technical Bulletin No. 1026, September
1950.
7. Combined Sewer Regulation and Man-
agement — Manual of Practice, FWQA —
11022 DMUO, July 1970.
8. Water Supply and Waste Water Disposal,
by G. M. Fair and J. C. Geyer - John Wiley
& Sons, Inc., New York, 1954.
9. Physical and Settling Characteristics of
Particulates in Storm and Sanitary Waste--
waters, by Robert J. Dalrymple, Stephen L.
Hodd, and David C. Morin, EPA-670/2-75-011,
April 1975.
50
-------
APPENDIX A
HYDRAULIC MODEL STUDY
This report covers the hydraulic studies
of the helical bend application to the
separation of solids from combined sewer
flows as carried out at the LaSalle Hydraulic
Laboratory, Montreal, Quebec, Canada.
The most complete investigation of the
helical bend principle, as applied to the
current investigation, was carried out by Dr.
T. M. Prus-Chacinski, as reported in his 1955
thesis at Imperial College of Science, London,
titled The Secondary Flow in a Meandering
Channel.2 He augmented his 1955 thesis with
further studies as a Partner of C. H. Dobbie
and Partners, Consulting Engineers, London,
at the University of Surrey and the Mersey
River Authority Hydraulic Laboratory in
Warrington, near Liverpool. These studies led
to the construction of a prototype at
Nantwich, Cheshire, England, for research on
behalf of the Construction Industry Research
and Information Association.
As flow proceeds around a bend,
secondary currents develop, with the upper
volume tending to move toward the outside,
while the bottom current angle is directed
toward the inside of the bend. Actual trace
studies made of this secondary motion show
that it describes a spiral or helical pattern.
This phenomenon is well known in rivers
where the deepest section of the channel is
usually scoured on the outside of any bends,
and the fine sand has been buoyed up by the
bottom currents and deposited at the inside.
In the present study of the bend
principles, the object was to impose just the
opposite geometry from that of a natural
stream — i.e., to have the deepest point at the
inside of the bend. In this manner, the
pollutants would be carried by the helical
flow down into the deepest section where
they would have optimum opportunity to
settle out. Furthermore, the flow velocities in
the study were to be controlled in the
transition and straight sections so the solids
would begin to concentrate in the lower
stratum of the flow. The combined effect of
these two phenomena would be to settle the
solids out of the liquid flow, and to
concentrate them in the deepest section of
the chamber, where they could be delivered,
via a foul sewer, to the treatment plant. The
major part of the liquid flow would then be
clarified and could be discharged over a weir
for discharge to the receiving stream, or to
subsequent holding or treatment processes.
Starting with a basic structure layout
established by Dr. Prus-Chacinski, the object
of the study was to evaluate the effects of
changes in width, depth, length of approach
channels, and discharge and conduit
configurations on solids separating efficiency.
The inlet pipe to the model was selected
as a 15.2 cm (6 in.) diameter P.V.C. pipe.
Based on the agreed geometric proportions to
start the study, this first dimension
established the model layout as shown on
Figure 4 of the foregoing report, Helical Bend
Regulator—Separator, Form 3 (Model Layout),
and as depicted in Figure 24, Downstream
View of Model, and Figure 25, Upstream
View of Model.
Unmeasured water flows were delivered
through the supply pipe, which was fitted
with a discharge control valve and solids
injection equipment, shown in Figure 26,
Solids Injection Vibrator. The original
structure began at the entrance to the
transition section; this was made of Plexiglas®
and was 20 diameters (20D) long to ensure
adequate flow deceleration. Following this, a
straight length of 10 D, also constructed of
Plexiglas, carried the flow to the start of the
bend.
The bend itself was built of Plexiglas,
using a 2.44-m (8 ft) radius from the
extension of the inlet pipe centerline. Since
the model was intended to demonstrate the
basic principles of the secondary motions
being .studied, the bend was carried around
through an angle of 120°. Side channels were
provided on each side to accommodate the
clear overflow spill from either the inside or
the outside of the bend. Control of either the
whole flow or just the small foul outflow was
provided by a perforated slide gate at the
120° position.
51
-------
FIGURE 24 DOWNSTREAM VIEW OF MODEL
FIGURE 25 UPSTREAM VIEW OF MODEL
52
-------
FIGURE 26 SOLIDS INJECTION VIBRATOR
The basic floor cross section inside the
transition, straight sections, and the bend was
made of polished cement mortar, as shown in
Figure 4 of the foregoing report. Simple
modifications to the Plexiglas overflow
weirs could be made to either change their
positions horizontally or their levels.
Flow spilling over the side weir was
directed to a calibrated Rehbock weir basin,
where the water levels could be read and the
discharge determined. Similarly, the foul
outflow was taken to a second basin for
measurement over a small 90° V-notch weir,
as was shown in Figure 5, Discharge
Measuring Weirs.
Injections of Process Black ink were
introduced on the floor of the bend. They
created clearly defined trace patterns which
could be read with a protractor to measure
the angle of the flow along the floor.
Tests ^involving solids recovery were
performed, using a standard procedure. One
liter of the selected model solids was injected
at a predetermined rate into the supply pipe.
The material entered the structure in the
flow, and was subjected to the helical
separation process. Any of the solids that
went over the weir were caught on a screen in
the measuring basin. The rest normally
remained as a deposit in the bend, and were
recaptured after the test by means of a screen.
The volumes of the two fractions were
measured to determine that no material had
been lost. The recovery rate, or efficiency,
was then computed as the ratio if the solids
retained in the be'nd with respect to the
original 1-liter input, expressed as a
percentage.
The model scale selected for comparisons
to a reasonable prototype was 1/6. The
pertinent scale relations are:
Scale =1/6
Time = 1/^/6 = 1/2.45
Discharge = 1/6 5/2 = 1/88
The prototype storm discharge with the
pipe flowing full in the corresponding
prototype 91.5-cm (3-ft) sewer was computed
as 0.85 m3/sec (30 cfs) which would
correspond to 9.6 I/sec (0.339 cfs). In fact, to
simplify operations on the model, a discharge
of 10 I/sec (0.354 cfs) was used, giving a
storm discharge of 0.88 m3/sec (31.1 cfs).
Dry-weather flow (dwf) was defined as
being one percent of the storm flow, or on
the model 0.1 I/sec (0.004 cfs), and in the
prototype 8.8 I/sec (0.31 cfs) through the
foul outlet, to the interceptor sewer and the
treatment works.
The methods used for solids simulation
were presented in Section II.
The first tests were carried out with the
model as shown on Figure 4, and overflow
spill over the full 120° bend, for discharges of
5, 10, 15 and 20 I/sec (0.177, 0.354, 0.531
and 0.708 cfs, respectively). The helical
motion was found to be very well defined.
Measurements were made of the floor flow
angles at the 35° position for straight-through
flow, i.e., no spill over the weir, as shown on
Figure 27, Top and Bed Flow Angles. The
method for determining the angular flow
along the floor of the helical bend is shown in
Figure 28, Floor Flow Angles Measuring
Equipment.
The resulting angles were treated,
following Dr. Prus-Chacinski's equation, in
the form: _ . ' ,
, Tana0= Cl Wr.Rll*
where: ' c "
&o = floor angle
Ci= coefficient
W = wetted perimeter
rc = central radius of
bend
Rn = Reynolds number
and plotted on a partial copy of Figure 28
from his thesis, and included here as Figure
53
-------
- Surface flow direction
- Bottom flow direction
FIGURE 27 TOP AND BED FLOW ANGLES
29, Bed Angle Relationships for
Straight-Through Flow. The preliminary test
points fell just below Dr. Prus-Chacinski's
work for rectangular and semi-circular
channels, but they were fairly close to those
for triangular channels with the deepest point
at the outside of the bend.
The fact that the points were lower than
the Prus-Chacinski findings would normally
indicate that the flow angle on the floor was
less, from which it would be deduced that the
helical action was weaker. However,
observation in the model indicated just the
opposite. Although the angle was perhaps
lower, the helical flow was very strong, and
the angular flow was acting over the full
width of the floor. Brief tests with shredded
Petrothene ® showed the model to be very
efficient in concentrating the settleables
toward the inner wall of the bend.
Further tests were performed, using
various discharges up to 20 I/sec (0.708 cfs),
first with spill over the full 120° bend. At this
maximum discharge, deposits containing
about 70 percent of the injected material
appeared in the 60° to 70° range. A
temporary wall was placed at the 60° position
and Petrothene® was introduced again. About
60 percent of the material was deposited
between the 30° and 50 ° points. From these
first results, it appeared that careful attention
should be paid to moving the solids through
the structure to the foul outlet during a storm
FIGURE 28 FLOOR FLOW ANGLE
MEASURING EQUIPMENT
so as to prevent a build-up of so large a
deposit that the dry-weather flow might be
unable to move it out of the helical bend
section.
Scale Selection and Discharges
As a result of the preliminary tests
described above, the most likely acceptable
discharge range had been defined in the
model, so a firm scale was selected for further
testing. This was taken as 1/6, meaning that
the prototype inflow pipe would be 91.5-cm
(3 ft) diameter. The selected discharges in the
model were 10, 7.5, 5, and 2.5 I/sec, (0.354,
0.265, 0.177 and 0.088 cfs), corresponding to
0.88, 0.66, 0.44, and 0.22 m3/sec (31.1, 23.3,
15.55 and 7.78 cfs) in the prototype.
The dry-weather flow (dwf) was specified
as being one percent of the storm'flow, or 8.8
I/sec (0.31 cfs) in the prototype. This would
normally flow through the structure at a
shallow depth along the invert, and discharge
through the foul outlet to the treatment
plant. It was further decided that during any
storm event, three times the dry-weather flow
(3 dwf) would be discharged through the foul
outlet - i.e., 26.4 I/sec (0.93 cfs).
A series of rapid modifications were made
in the model and tested in an attempt to
determine the importance of the geometric
54
-------
5.0
4.0
3.O
£
^
$ 2-°
c
£
l.O
<#>
V
-%*
Or-1
°C
^?
OQC
<
«L_J
Vi
6
o
•> o
p —
£
30C
O
'A'
e
1
o
V*
0
?!
c
tJ
0
A
1 °
o n o
^&1° o<
^§3^1^
ooo~^
A 0 <^
e*9«
-^^H
«••
p^fec
£?c§
0 ^
3 a
-^
345 JO 20 30 40 50 607080
Rnx I0~3
o Rectangular Section (Prus)
j&. Triangular Section (Prus)
^p Semi-circular Section (Prus)
*• LHL Tests
FIGURE 29 BED ANGLE RELATIONSHIPS FOR STRAIGHT-THROUGH FLOW
variables. All four discharges were run for a
few tests, but it quickly became evident that
only the 0.88 m3/sec (31.1 cfs) should' be
used for comparisons, or at times, the 0.66
m3/sec (23.3 cfs). Petrothene was used as the
settleable material for most tests, with only
one or two check runs using Gilsonite.
Data given in Table 9, Recovery Results
— Modifications 1-8, show the recovery
results for the first series of modifications
tested. All cases shown are for the 0.88
m3/sec (31.1 cfs) overflow discharge and 26.4
I/sec (0.93 cfs) foul outflow, with
Petrothene used as the test material. All
values are for prototype dimensions.
From the outset it was obvious that the
long transition and straight sections formed a
very efficient settling chamber. With
Modifications 1 and 2, it was noted that all
the Petrothene was deposited on the floor in
the bend, so changes were made in an effort
to cause the material to flow through to the
foul outlet. This pattern seemed difficult to
achieve with a short bend, so the full 120°
length was tested. (Test 9) with the bleeder at
55°.
The recovery efficiency was reduced so
much that this approach was given no further
consideration. The bleeder was a
small-diameter discharge pipe located as
shown in Figure 30, Modification 6 — Tests
23-29.
The taper in Modification 4 induced
much stronger helical motion, to such a
degree as to cause an upflow and to carry
more material up over the weir. The change is
55
-------
Foul
Outflow
Underflow
Bleeder
SECTION A-A
FIGURE 30 MODIFICATION 6 - TESTS 23-29
56
-------
TABLE 9
RECOVERY RESULTS - MODIFICATIONS 1-8
Modification
Number
1
2
3
Orig.
4
5
6
7
8
Ref.
Fig.
9
9
9
1
10
11
12
13
14
Test
No..
1
4
7
9
13
16
17
20
22
23
28 ,
30
32
34
Bend
Length
60°
50°
40°
120°
60°
60°
60°
60°
60°
60°
60°
50°
60°
60°
Weir
Length
60°
50°
40°
120°
10° -50°
10° -50°
60°
15° -60°
60°
60°
60°
50°
60°
60°
Weir Height
m. ft.
1.64(5.38)
1.64(5.38)
1.64(5.38)
1.64(5.38)
1.26(4.13)
1.64(5.38)
1.64(5.38)
1.64(5.38)
1.33(4.38)
1.33 (4.38)
1.64(5.38)
1.64(5.38)
1.64(5.38)
1.39(4.55)
Recovery
%
100
85
82
50
( through bleeder at 55°)
90
90
75
60
60
65
70
(through bleeder at 55°)
85
85
83
shown in Figure 31, Modification 4 -
Tapered Outlet. Li spite of this increased
secondary flow, the Petrothene recovered
was still all on the floor of the bend, most of
it ahead of .the 35° position. Weir height and
length changes did not have much apparent
effect.
Modification 5 was the first attempt to
enter the bend at a higher velocity so the
Petrothene would be carried through to the
foul outlet. The layout is shown in Figure 32,
Modification 5 - Narrow Inlet, 2D Wide; and
Figure 33, Modification 5 — Tapered Bend.
The recovery rates dropped slightly, but the
bottom deposits had moved along to the 35°
to 55° position. The outstanding feature of
this layout was the surface eddy which
formed at the end of the helical bend near the
60° position. At times it entrained material
from the bottom, causing it to upsurge to the
surface and to be carried over the overflow
weir.
Narrowing the end of the bend in
Modification 6 did not result in carrying the-
material any further toward the foul outlet.
Test 28, with the bleeder at 55° operating,
produced a slight improvement.
Shortening the bend to 50° in
Modification 7 gave encouraging results, with
an increased recovery. Modification 8 was
tested as a means of having a simpler,
parallel-sided bend: Recoveries were good, but
no improvement over the tapered format was
evident.
Although slightly less Petrothene was
recovered in a 1.83 m-wide (6 ft) channel (2
D) than the 2.745 m (9 ft) channel (3 D), it
was investigated to see if there was any merit
in achieving a narrower structure. (Gravity
and dynamics of motion provide removal, not
the structure). The channel was finally
reduced to a 1.372 m (4 ft - 6 in.) width (1.5
D), to check on the minimum reasonable
dimensions. This configuration is shown in
Figure 34, Minimum Width Inlet, 1.5 D Wide;
and Figure 35, Narrowest Bend. 1.5 D Wide.
Tests were carried out for the 0.88
m3/sec (31J cfs) and 0.66 m3/sec (23,3 cfs),
cases in each instance For both total
discharges, the foul discharge was maintained
constant at 26.4 I/sec (0.93 cfs), or 3 dwf.
Results of these tests are given in Table 10,
57
-------
•:*•••-•*
•; h
•":: '» • ,„:!' » •— t,: attaLfa :s safl a • "jwitntt ta»!ta 4 sti rtiii '.iiiutfMte
: ii'j"1!" iijTB, !!Ji;\jiii!i:ii;iijiiii!!:i;:^liji;!M:i7i!ii;iJ|ii|!!!i;7 r'i'J!!!1'^ "'Bfiis uf^jjiii:*!1^!;;):!
FIGURE 31 MODIFICATION 4 - TAPERED OUTLET
FIGURE 32 MODIFICATION 5 -
NARROW INLET 2D WIDE
FIGURE 33 MODIFICATION 5 -
TAPERED BEND
58
-------
FIGURE 34 MINIMUM WIDTH INLET,
1.5DWIDE
FIGURE 35 NARROWEST BEND,
1.5DWIDE
Recovery Results, Modifications 9, 10, 12,
14-18.
Figure36, Modification 9 - Tests 38-39,
shows that the bend was shortened to 50°.
The object was to attempt to eliminate the
deposition of material on the sloping floor
between 45° and the end of the bend. In Test
32, this deposit was considerable. Test 38 was
only partially successful in eliminating the
material.
A truncation was made at the end of the
channel in Modification 10, as shown in
Figure 37, Modificiation 10 - Tests 40-41.
This, too, was insufficient to prevent the
deposit. It did, however, slightly reduce the
back-eddy in the end zone.
Figure 38, Modification 12 - Tests
46-47, shows that a Z-shaped partition was
provided to allow flow over 60° of the weir,
with the foul outlet placed at 50°. This offset
arrangement aggravated the back-eddy to the
extent that all material reaching beyond 45°
was uplifted and carried over the weir
immediately without any chance of settling.
At this stage, the first tests of shortening
the inlet transition were tried, as shown in
TABLE 10
RECOVERY RESULTS - MODIFICATIONS 9, 10, 12, 14-18
Modification
Number
9
10
12
14
15
16
17
18
Ref.
Fig.
15
16
17
18
18
19
21
21
Test
No.
38
40
46
53
55
57
59
60
Bend
Length
50°
50°
60°
60°
60°
60°
60°
60°
Channel
Width
m
183
183
183
183
183
1.372
1.372
1.372
ft.
(6)
\6)
(6)
(6)
(6)
(4.5)
(4.5)
(4.5)
Recovery
Q=0.85 I/sec
(30 cfs)
87
88
82
65
63
65
Q=0.64 I/sec
(22.5 cfs)
100
100
96
77
89
85
85
70
59
-------
Foul
Outflow
SECTION A-A
FIGURE 36 MODIFICATION 9 -
TESTS 38-39
Foul
Outflow
SECTION A-A
Sloping Triangle In Corn<
Foul
Outflow
SECTION A-A
FIGURE 37 MODIFICATION 10 -
TESTS 40-41
FIGURE 38 MODIFICATION 12 -
TESTS 46-47
FIGURE 39 MODIFICATIONS 14 & 15 -
TESTS 53-55
60
-------
Figure 39, Modification 14 and 15 — Tests
53-55. This reduced the transition to 10D, or
9.15m (30 ft) prototype. There was no roof
in Modification 14 to guide the flow, as
shown on Figure 39. Test 53 showed a sharp
drop in recovery — i.e., 65 percent as
compared to 85 percent in Test 32. In the
model, it was visually evident that the flow
entering the transition did not expand
efficiently, so that at the 0° position entry to
the bend, considerable turbulence and lack of
flow uniformity across the section were
experienced.
A sloping roof was placed in the
transition for Modification 15. This change
was still insufficient to guide the flow, and
too much irregularity remained at the 0°
position. It can be concluded that with the
asymmetrical entrance, 10 D is an insufficient
length for expansion from the 0.915
m-diameter (3 ft) inlet pipe, D, to the 1.83 m
(6 ft) channel, 2 D wide.
At this stage, it was decided that there
would be insufficient economical justification
for providing a narrower channel. Once the
excavation is made for a particular structure,
the cost difference between a channel 2 D and
3 D wide would be minimal. This reasoning
led to resumption of studies of a wider
section.
However, before reverting to a wider
section, a few brief tests were carried out with
the section narrowed down even further to
1.5 D, or 1.372 m (4.5 ft) prototype. The
same shortened inlet was retained, and flow at
the 0° position appeared to be fairly evenly
distributed. Of course, the flow rate was
much faster than in the preceding tests. The
recovery rates of 65 percent for 0.88 m3/sec
(31.1 cfs) and 85 percent for 0.66 m3/sec
(23.3 cfs) seemed to conform well when
compared with the wider channels, as shown
on Figure 40, Petrothene Recovery in Model
as Function of Bend Channel Width.
Even with this narrowed channel, most of
the material recovered was in the form of
deposits on the bend floor. Two minor
changes were investigated in an effort to move
this material along to the foul outlet. The first
was a series of corner blocks 15.2 cm (6 in.)
high and 22.95 cm (9 in.), projecting out
Inlet Pipe
Diameters
Q=O.66n
3/s (23.3
*'•'
0.88 m% (31.1 cfs)
D ISO i.0
BEND CHANNEL. WIDTH
FIGURE 40 PETROTHENE RECOVERY IN
MODEL AS FUNCTION OF
BEND CHANNEL WIDTH
from the side wall, as shown in Figure 41,
Modifications 17 and 18 - Tests 60-61. Only
the 0.66 m3/sec (23.3 cfs) case was tested; it
Floor Deflectors Tried
in Modification 18
(with Blocks Removed)
'liu0.1" High Corner Blocks
\o I
FIGURE 41 MODIFICATIONS 17 & 18
TESTS 60-61
61
-------
gave an 85-percent recovery. This represented
practically no change from the configuration
without the corner blocks. Next, the blocks
were removed, and a series of floor baffles
7.62 cm (3 in.) high were installed, as shown
in Figure 41.
These baffles generated so strong a
secondary current, that the Petrothene
deposits, were re-suspended and carried over
the weir. The recoveries dropped to 50
percent and 70 percent.
The preceding studies resulted in
resumption of interest in the bend with the
full 3 D or 2.745-m (9 ft) width originally
built in the model. The studies showed that
this bend, 60° long, when supplied through
the original inlet configuration shown on
Figure 42, Inlet Configurations Tested,
recovered 100 percent of the Petrothene for
feat
Outflow
D =
FIGURE 42 - INLET CONFIGURATIONS
TESTED
FIGURE 43 TRANSITION 10D LONG
a 0.88-m 3/sec (31.1 cfs) discharge.
Therefore, tests were undertaken to see how
much the inlet length could be reduced and
still retain an acceptable rate of recovery. No
scum boards were used in these tests.
As shown on Figure 42. Alteration A
consisted of advancing the 0.915 m-diameter
(3 ft) inlet pipe 9.15 m (30 ft), providing the
transition out to the 2.745-m (9 ft) width in
9.15 m (30 ft), or 10 D, and retaining a 9.15
m (30 ft) straight section as depicted in
Figure 43, Transition 10 D Long. It was
evident that the flow jet did not expand well
in this short transition and significant
turbulence entered the bend, disturbing the
helical secondary flow. Recoveries for 0.66
m3/sec (23.3 cfs), and 0.88 m3/sec (31.1 cfs),
respectively, were 95 percent and 89 percent,
as shown on Figure 44, Petrothene
Recovery in 60° Bend for Various Inlet
Configurations.
62
-------
100
£
>l
o
I 90
«
Q:
0
w
JC
% 80
a.
j>
V,
\
-x
i ^
r^
\ /-
X ^/^
"x
V ^
V
\
, Q _ O.S6m3/3
P (23.3 cfs )
I
Q 0.88 m3/s
""" (31.1 cfs)
I
Original ABC
Inlet Configurations- Ref. See Fig 22
Jote: Bend Width- 3D= 2.74m (9-0")
FIGURE 44 PETROTHENE RECOVERY IN 60° BEND
FOR VARIOUS INLET CONFIGURATIONS
Alteration B, Figure 42, provided the
start of the transition from the same point,
but it expanded over 13.73 m (45 ft) or 15 D
to reach the 2.745-m (9 ft) width. This is
shown in Figure 45, Transition 15 D Long.
Just 4.58 m (15 ft) of straight line remained
ahead of the bend. Flow through this
configuration was much improved; it
expanded smoothly to fill the cross section
effectively before entering the bend. As
shown in Figure 44, the recoveries for the
0.66 ma/sec (23.3 cfs) and 0.88 m3 /sec (31.1
cfs) discharges rose to 98 percent and 93
percent, respectively. Compared with the-
original long transition and straight sections,
the flow ebtering the bend now was slightly
more turbulent, but helical motion was still
generated.
Since the transition in Alteration B
appeared to be operating well, it was moved
forward so it reached its full width just at the
bend entrance. This was Alteration C, shown
in Figure 42. The turbulence entering the
bend was too strong in this case; Figure 44
shows the drop in recoveries.
The obvious recommendation was to
retain Alteration B, and this was used for the
following tests with scum boards. The
reasoning for this choice was based on the
fact that it resulted in only a seven percent
drop in recovery from the original case, with a
structure that is 9.15 m (30 ft) shorter.
Implicit in this configuration is the
requirement for a transition 15 D long to
provide satisfactory deceleration of the flow
while expanding to the 3 D channel width in
the straight section. A flat sloping roof covers
the transition, starting at the height D and
rising to 2 D over the 15 D length.
As shown on the cross section in Figure
46, Scum Board Locations Tested, the scum
board was tested at four different depths of
immersion and at three distances out from the
weir lip. For each of these locations, the
Petrothene recovery was determined for the
0.88-m3/sec (31.3 cfs) discharge, with both
the 50° and 60° bends. Figure 47.
Petrothene Recovery in 50° Bend for
Various Scum Board Locations, and Figure
48, Petrothene Recovery in 60° Bend for
Various Scum Board Locations, show the test
results.
63
-------
FIGURE 45 TRANSITION 15D LONG
It was immediately evident that the
deepest scum board locations, 45.8 cm (18
in.), or D/2, below the weir lip, seriously
disturbed the helical flow. Figure 49, Floor
Flow Angles with Scum Board D/2 Deep,
shows the bottom flow pattern. The
recoveries fell off sharply, as well, and brief
tests with Process Black ink showed that the
floor angles were severely modified, even
rising toward the weir at some points.
There was some improvement when the
scum board immersion was reduced to 30.4
cm (1 ft) below the lip, particularly for the
60° bend length.
When the Scum board was just 15.2 cm
(6 in.) below the lip, the recoveries rose
significantly. In the 50° bend length, it
reached about 84 percent, whereas, for the
free overflow case — i.e., with no scum board
— the recovery was 87 percent, the optimum
in these tests. Figure 50, Floor Flow Angles
With Scum Board D/6 Deep, shows the effect;
of deep immersion.
In the 60° bend length, recoveries with
the scum board reached 85 percent, but the
free overflow performance was 93 percent,
leaving a wider gap in performance. On this
basis, the scum board was raised to just 7.6
cm (3 in.) below the lip. However, instead of
I
;a
{J£
•in
o
oo
«
e -
_ o
5
o
o
ro
O
in
-------
Hot.: Bw« WMr» - 30. 1.74,, illtt
0* O.tSm'/t (SI. let.)
... *««•> »oor.] Dipt*
x No Sn. I««M /Strew tNk Up- endn)
•<
£•0
i
t
J70
•O
\
/ f
-i
>^ — ^^
— ^^^^
r-t-
^--wa (•>
X4S«(«)
^».8(a)
!" * *'
—
^
«
-i — . i , i , |
Sew •.«< Oi»l«iM Out Fro. IM, Lip
HOH: B.nJ width - 30= 2.74m .1 ft I
Q» O.86 msA.(31.1 eft)
ScumBoord Otplh
No Scum Board S Btlow Wlir L[p-cm(in)
2
§
* i
|
I,
.f KO scum Board /
es=
Scum Board Diitonci Out From W*[r Lip
FIGURE 47 PETROTHENE RECOVERY
IN 50° BEND FOR VARIOUS
SCUM BOARD LOCATIONS
FIGURE 48 PETROTHENE RECOVERY IN
60° BEND FOR VARIOUS
SCUM BOARD LOCATIONS
FIGURE 49 FLOOR
FLOW ANGLES WITH
SCUM BOARD D/2 DEEP
a recovery increase, a slight decrease was
experienced.
The margin which remained between the
recoveries for the no-scum board
configuration and the shallowest scum board
cases shown for the 60° bend in Figure 48
was from 93 percent down to 85 percent.
This seemed too high a price to pay for
adding the scum board. Observations on the
model indicated that placing the scum board
out in front of the weir, right in the main
flow through the bend, caused significant
disturbance. Also, previous work had shown
that the 50° bend with scum board could not
exceed 85 percent recovery, whereas the 60°
bend showed the possibility of reaching better
than 90 percent efficiency. Therefore, the 60°
length was retained for further development.
Brief intermediate tests with the weir set
back from the bend channel limit, and the
scum board on the limit, led to the offset weir
arrangement shown in Figure 51, Cross
Section of Offset Weir in Bend. Weir details
are shown in Figure 52, 60° Bend With Offset
65
-------
r
FIGURE 50 FLOOR FLOW ANGLES WITH SCUM BOARD D/6 DEEP
Water Level for 0.88 rtiYs (3l.lcfs)
FIGURE 51 CROSS SECTION OF OFFSET WEIR IN BEND
Weir; Figure 53, Upstream View of Scum
Board on Offset Weir; in Figure 54,
Downstream View of Bend and Offset Weir;
and in Figure 55, Scum Board on Offset Weir.
Tests run with this cross section and a
60° bend showed that the scum board did not
adversely affect the recovery rates. At this
stage, therefore, a series of tests was carried
out using various inlet configurations.
Earlier studies had shown that the
transition length had to be 15 D to give
adequate flow deceleration. This form was
retained, and it was tested for a discharge of
0.88 m3/sec (31.1 cfs), first with no straight
section, as shown in Alteration C in Figure
42, then with straight sections 5D and 10D
long (Alterations B and D on Figure 42). The
corresponding Petrothene recovery rates are
66
-------
FIGURE 52 60° BEND WITH OFFSET WEIR
FIGURE 53 UPSTREAM VIEW OF SCUM
BOARD ON OFFSET WEIR
FIGURE 54 DOWNSTREAM VIEW OF BEND
AND OFFSET WEIR
FIGURE 55 SCUM BOARD ON OFFSET WEIR
67
-------
100
8
S" so
>
o
o
0)
DC
-------
FIGURE 57 PROVING TESTS LAYOUT
on Figure 52, the plan dimensions in Figure
57, Proving Tests Layout, were related for
detailed testing over the range of discharges in
the model. Both Gilsonite and Petrothene
were used in these tests and the results are
shown on Figure 58, Solids Recovery in
Model.
In predicting the prototype recoveries
expected with grit, reference was made to
Figure 9 showing the sizes simulated by the
Gilsonite and Petrothene materials used in
the model. It appeared that over practically
the whole range of scales, the Gilsonite
represented the grit satisfactorily. Figure 11,
Predicted Grit Recovery for All Scales, was
then prepared, using the Gilsonite recovery
rate appearing in Figure 58.
Using Gilsonite to simulate prototype
organic material in the model as shown on
Figure 10 covered the major portion of larger
particle sizes. The upper portion appearing in
Figure 10 for each given scale was multiplied
by the Gilsonite recovery rates from Figure
58 for QD , 1.5 QD and 2 QD , and taken as
the prototype recovery for the particle sizes
in that zone.
The particle sizes lying below the dm =
0.5 mm line for Gilsonite in Figure 10 were
assumed to be entirely represented by the
Petrothene. Similarly, this portion for each
scale was multiplied by the Petrothene
recovery rates from Figure 61 for QD, 1.5
QD, and 2 QD, and taken as the prototype
recovery for the lower zone. '
These two prototype recoveries simulated
by Gilsonite and Petrothene in the model
were added together and presented on Figure
12, Predicted Prototype Settleable Organic
Material Recovery.
It is, therefore, possible to find the
predicted recovery rates for grit and organic.
materials as defined in Figure 6, for any scale
of prototype structure by using Figure 11,
and Figure 59, Velocity Contours on 0° Cross
Section.
It was found during the development
testing that the solids materials recovered
were normally in the form of deposits on the
floor of the bend. Only traces, normally less
than two percent, were carried through the
foul outlet during the storm flow incidents.
For most tests, the Petrothene and Gilsonite
were deposited or shoaled in the bend
between the 15° and 60° locations. A fair
proportion of the Gilsonite settled to the
floor in the straight section, but it was
subsequently moved into the bend.
As the discharge was reduced in the
model, as it returned to dry-weather flow,
most of the depostis were scoured out and
discharged through the foul outlet. However,
it is possible that this self-scouring
phenomenon during dry-weather flow in the
model simualtion may not take place in any
prototype installation, so provision should be
made in design for some form of wash-down
facility or for some other form of regular
maintenance.
69
-------
100
90
80
70
55 60
|50
I 40
cc
30
20
10
0'
Gils
Pe
, .^_
mite on Model —
rothene on Mod
\ y
y"x
»l ^ ^
*— — — .
Iv
X
X
\^
^
1
5 10 V 15 20 1/s
Model Discharge —^CT
, > ,
0-1 0.2 0.3 0.4 0.5 0.6 0.7
111!
cfs
0 0.5QD Q0 1.5 QD 2.0Q0
FIGURE 58 SOLIDS RECOVERY IN MODEL
NOTE: Velocities in ft/sec (to get cm/sec , multiply by 30.5)
Discharge = 0.88 mVsec (31.1 cfs )
FIGURE 59 VELOCITY CONTOURS ON 0° CROSS SECTION
70
-------
10°
NOTE: Velocities in ft/sec (to get cm/sec , multiply by 30.5)
Discharge = 0.88 m3/sec (31.1 cfs )
FIGURE 60 VELOCITY CONTOURS ON 10° & 20° CROSS SECTIONS
71
-------
30
NOTE: Velocities in ft/sec (to gel cm/sec,multiply by 30.5)
Discharge = 0.88 mVsec (31.1 cfs)
FIGURE 61 VELOCITY CONTOURS ON 30° & 40° CROSS SECTIONS
72
-------
FIGURE 62 SOLIDS CONCENTRATION SAMPLING JIG
fj£r*~
FIGURE 63 SIMULTANEOUS FILLING OF SOLIDS SAMPLING BOTTLES
73
-------
NOTE: Concentrations given in relative terms only
FIGURE 64 SUSPENDED SOLIDS CONCENTRATIONS ON
0° CROSS SECTION
74
-------
NOTE: Concentratrations given in relative terms only
FIGURE 65 SUSPENDED SOLIDS CONCENTRATIONS ON
10° & 20° CROSS SECTIONS
75
-------
30'
NOTE: Concentrations given in relative terms only
FIGURE 66 SUSPENDED SOLIDS CONCENTRATIONS ON
30° & 40° CROSS SECTIONS
76
-------
Q)
Q.
O
U
U>
Ul
O
CO
0>
(T
10 20 30 40
Position Along Weir
50
20 Sections
J
60 Degrees
FIGURE 67 DISTRIBUTION OF SOLIDS WHICH
ESCAPE OVER THE WEIR
77
-------
Floor flow angles for
120° bend with no
overflow
Floor flow angles for
60° bend, free overflow
by offset weir
Floor flow angles for
60° bend with scum
board on offset weir
FIGURE 68 FLOOR FLOW ANGLES
78
-------
Another approach would be to provide
sufficient slope from the inlet to the foul
outlet of the structure to induce all deposited
solids to be scoured out and flushed away
with the dry-weather flow. The model invert
had no slope in any of the tests.
Testing the Performance of the
Recommended Layout
Utilizing the model forms shown in
Figures 51 and 57, final proving tests were
carried out to define the structure's flow
characteristics, to serve as guidance for the
work being carried out on the mathematical
model.
First, the .tangential flow velocities in the
bend were measured with a midget propeller
current meter for the discharge of 0.88
m3/sec (31.1 cfs). Readings were taken at 14
points on each of the cross sections, at the 0°,
10°, 20°, and 40° positions. The prototype
velocities found were plotted at the
corresponding points, and contours of even
velocity values were drawn as shown on
Figures 59, 60, and 61, Velocity Contours on
0°; on 10° and 20°; and on 30° and 40° Cross
Sections.
The same cross section locations were
used for installing a special suspended solids
sampling device, as depicted in Figure 62,
Solids Concentration Sampling Jig. Twelve
6.4-mm (1/4 in.) I.D. Plexiglas sampling
tubes were fixed on a frame, directed
upstream to intercept the flow. A flexible
tygon tube ran from each sampling tube to a
wooden stand, under which 1-liter sampling
bottles could be introduced, as shown in
Figure 63, Simultaneous Filling of Solids
Sampling Bottles.
With the model subjected to clear water
only, each of the tubes was primed and they
were all adjusted to provide the same
discharge. Petrothene injection was begun by
means of the vibrator, and as soon as the
cloud of suspended material reached an
equilibrium state at the cross section being
tested, the bottles were slipped under the 12
tygon tube discharges. This system provided
simultaneous samples of the suspended solids
in the flow at all 12 points.
Results of these tests are shown in
Figures 64, 65, and 66, Suspended Solids
Concentration Contours on 0°; on 10° and
20°; and on 30° and 40° Cross Sections.
The length of the overflow weir in its
offset position was 302 cm (118.8 in.). A
catch screen was fitted outside the weir and
divided into 19 sections, each 15.2 cm (6 in.)
long, with the twentieth being just 12.2 cm
(4.8 in.) long. The standard Petrothene
injection procedure was followed, and the
particles escaping over the weir were caught in
the corresponding sections on the screen. The
corresponding fractions in each section were
measured and plotted as shown on Figure 67,
Distribution of Solids Which Escape Over
Weir. '
Using the special Process Black ink, the
floor flow angles (a0 on Figure 68) were
measured at six points across the cross
sections at 10°, 20° and 35°, for each of three
separate layouts:
1.
2.
120 bend with straight through flow
60 bend with offset weir, but no scum
board
3. 60° bend with offset weir and scum
board.
The averages of these readings were
computed and plotted in Figure 69,
Comparative Average Floor Angles for
Recommended Layout. The floor flow
phenomena are shown in Figure 68, Floor
Flow Angles for 120° Bend With No
Overflow; for 60° Bend. Free Overflow by
Offset Weir; and for 60° Bend With Scum
Board on Offset Weir. The average of floor
angles is plotted in Figure 69, Comparative
Average Floor Angles for Recommended
Layout.
79
-------
30°
120
Straight
Through
No
Scum Board
60°
Scum
Board
FIGURE 69 COMPARATIVE AVERAGE FLOOR ANGLES FOR
RECOMMENDED LAYOUT
80
-------
APPENDIX B
MATHEMATICAL MODEL STUDY
This report describes the development of
a mathematical model and computer
simulation of an open channel bend with a
side overflow weir for the separation of grit
and organic solids from combined sewer
overflows. A companion, hydraulic laboratory
investigation was conducted by LaSalle
Hydraulic Laboratory of Montreal.
The general features of the helical bend
device are illusrated in Figure 4, Helical Bend
Regulator - Form 3 (Model Layout) in
Section It The dimensions shown are for a
prototype designed for a nominal peak flow
of 0.85 m3/sec, (30 cfs). The flow from the
circular inlet sewer passes through a straight,
rectangular transition section designed to
spread the flow uniformly across,the larger
cross section of the device.
The flow then proceeds through a 60°
bend, with cross section as shown in Figure 4.
The clean overflow is withdrawn via a side
weir located along the outside of the bend,
while the concentrated solids are withdrawn
via a foul sewer orifice in the end plate at the
60° position. Downstream control is used to
limit the amount of flow discharged to the
wastewater treatment plant to three times the
dry-weather flow, or 0.025 m3 /sec, (0.9 cfs),
for the nominal prototype design. The
remainder of the inlet flow is discharged to
receiving waters over the side weir, or to any
holding or treatment facility which is
intended to handle the overflow liquid.
The usefulness of this device for
separating grit and organic solids from
combined sewer overflows is enhanced by the
secondary flow patterns which develop in the
bend. These secondary flows form a clockwise
spiral or helical flow pattern, in the direction
of flow. The helical flow tends to sweep the
bedload to the inside of the bend, so that
relatively clearer liquid is discharged over the
side weir. These secondary flow patterns have
been studied extensively by many
investigators, most notably by Dr. T. M.
Prus-Chacinski,1 who served as consultant on
the present study. His report, Appendix C,
contains an extensive bibliography of the
previous work in this field. •
The application of open channel bends to
separation of solids from storm and combined
sewer overflows has also been studied
experimentally by Dr. Prus Chacinski,2.3 and
a prototype device shown in Figure 3 in the
foregoing report, Photograph of Nantwich,
England, Helical Bend Regulator, has been
constructed by C. H. Dobbie and Partners,
Consulting Engineers, under a research
contract from the Construction Industry
Research and Information Association,
London. The objective of the present study
has been to refine and extend this previous
work in accordance with American conditions
and practices.
LIQUID FLOW CALCULATION
The secondary flow which develops in a
curved channel is of definite significance in
determining the separation performance of
the helicaLbend separator. The helical flow
pattern tends to sweep the bedload to the
inside of the bend, thereby reducing the
amount of solids which pass over the weir,
with the clarified combined sewer flow.
The helical secondary flow is induced by
the vertical velocity profile at the entrance to
the bend. An element of fluid near the center
of the channel experiences an outward
centrifugal force proportional to v2 /r, where v
is the longitudinal velocity and r is the radius.
This centrifugal force tends to raise the
surface level at the outer radius of the bend,
and .lower it at the inside of the bend. Near
the center of the channel, the centrifugal
force is, therefore, balanced by internal
pressure gradients caused by the slope of the
water surface. However, near the bottom
where the longitudinal velocity is lower, the
pressure gradient overbalances the smaller
centrifugal force, and the fluid near the
bottom is accelerated toward the inside face
of the bend. Conversely, at the surface where
the longitudinal velocity is greatest, the
centrifugal forces overbalance the pressure
gradient so that water at the surface is
accelerated toward.the outside of the bend.
8.1
-------
This results in a basic clockwise spiral flow
pattern. For large angle bends, the
redistribution of fluid by the secondary
motions can lead to situations where the
longitudinal velocity no longer increases
monotonically from the bottom - to the
surface. Under these conditions, multiple
helices may develop. However, even under
these conditions, the acceleration of the fluid
is controlled by the local balance of the
centrifugal force (v2 //•) and the pressure
gradient caused by the slope of the surface.
The important point in the preceding
discusssion is that the secondary motion is
essentially an inviscid flow phenomenon. The
action of viscosity is important in establishing
the vertical velocity profile in the flowing
liquid at the entrance to the bend. However,
given the initial vertical variations in velocity,
the development of the helical secondary flow
depends upon the balance of the centrifugal
and pressure gradient terms.in the equations
of motion. The complex flow in the channel
bend can, therefore, be adequately described
without including viscous terms in the
equations of motion.
This approach was adopted by Fox and
Ball4 in their computer model for the
secondary flow in bends with no side weir.
The present formulation represents an
extension of this work to include discharge
over a side weir, a more general channel cross
section, and calculation of the resulting
particulate concentration field and separation
efficiency.
Equations of Motion
The basic equations for steady state flow
of an incompressible fluid are the momentum
and continuity equations:
Momentum V • (VV) + I Vp = 0 (1)
Continuity V • V = 0,
where
(2)
V = velocity vector
p = fluid density
p = pressure
Consider first a curved rectangular
channel. Any point in the channel can be
located by the cylindrical coordinates (r, Q, z)
,as shown in Figure 70, Cylindrical Coordinate
System.
Corresponding respectively to the coordinates
(r,Q ,f) are the velocity components (u, v, w),
also shown in Figure 70. The general
momentum and continuity equations (1) and
(2) may be written in cylindrical coordinates:
Momentum "Ui- +
(3)
u 9v
dr
__
r 30
uv
r
-
Continuity f- + f + -J- f + f - 0 (5)
/' XN^
s' 4
p
\
\
\
\
?
/ CHA
CRO,
SEC1:
FIGURE 70 CYLINDRICAL COORDINATE SYSTEM
82
-------
Here, the original pressure gradient terms
have been replaced with derivatives off, which
represents slight variations in the fluid surface
elevations. 3f/3/- is a positive quantity
resulting from the fact that centrifugal force
causes the surface elevation at the outside of
the channel to be slightly higher than the
elevation at the inside of the channel. The
pressure gradient term of gd$/dr balances the
radial acceleration v2 jr. The factor 9f ,/3 e is
also positive due to the slight rise in surface
elevation caused by the channel blockage.
This factor has the effect of slowing the fluid
(decreasing average v) as the channel blockage
is approached.
For computational purposes, equations
(3), (4), and (5) are rewritten as follows:
|| = . _/• L 3« . _£ + w|«
oO v \ dr r az
3r
30 "y ' "~ +— +
^ +-S.3II
3z r 3flj
w = -
JF+f + ff]^
(6)
(7)
(8)
A simplified summary of the computational
procedure is as follows: At an angle 0,
suppose that at discrete points in r and f'
values are known for z4>,0,?), v(r,8 £), and
w(>,0,f)- The quantities on the right hand
sides of (6) and (7) can be calculated, and
values for 3«(r,0,f)/0 and dv(R,6,{/de
follow. Then, u and v can be projected
through an incremental angle A0 to the 0 +
A0 plane:
ufr, 0 +&B,zJ=u(r, 0, z) +(du/d6JAe (9)
vfr.e +&e,zj-v(re,z)+(dv/de]Ae (10)
These new values in the 0 + A0 plane are
then applied to Equation (8) from which
w(r, 0 + Af) follows. With this the cycle is
completed.
Radial Surface Slope
The terms on the right hand sides of
equations (6), (7), and (8) can all be evaluated
knowing u, v, and w in a given 6 plane, but
the terms df /dr and 3f/80 are not as obvious
as the rest. Figure 71, Definition of Surface
Shape in Channel Cross Section, shows a
channel cross section, still represented by a
rectangle, with an enlarged surface slope
component.
f
f.
J
ti
h ^
±^
y
''
_^<"~§
•— ^ ,
f^^\
\f^— j
V»j
^VX'VVVVVxV/C'/'X////
r
o
z
C (r, 6)
§
^
•4 »-r
f2
/
/
y
W/xxxxx
IN CHANNEL CROSS SECTION
Define q(r, 0) to be the lateral discharge
across a line of constant radius r:
h
qfr,8J = | u(r,e,z) dz
-f
(11)
Velocity components u(r0, 0, z) are
illustrated in Figure 71. It is realistic to
assume that the rate of lateral discharge does
not vary with angle 0, so 90 = 0, or
r h
90
The.Slope 3J/90 is very small, so
h
« dz =
dz
= 0
(12)
(13)
83
-------
However, (3zi/90) is known from Equation
(6), so,
h
**
\~_ru _du + „ _ ny 1" _I2 1£~] cfz = 0,
from which
_
iu u
v dr
v ,
(15)
All of the terms on the right hand side of this
expression follow directly from the
knowledge of u, v, and vv, so the resulting
g9f/9r may be substituted into Equation (6).
This technique for obtaining the radial surface
slope, 9f /9r, was developed by Fox and Ball.4
Longitudinal Surface Slope
The fluid flow solution is treated
computationally as an initial value problem,
so only initial conditions at the upstream end
of the base are needed. However, in making
this simplification, there is no way that the
fluid can feel the presence of the dam at the
downstream end of the channel, so a
mechanism is needed to slow the average
longitudinal velocity as it proceeds around
the channel arc. The first possibility in-
vestigated involved artificially increasing the
surface ?(r, 0) with increasing angle 0 . The
integral of the surface slope, ?t (r, 0), can be
found from (3f/3r) in Equation (15):
(16)
where r! is the inner radius of the channel
bend. The average surface height is then:
(17)
where r2 is the outer radius of the bend.
The resulting surface elevation, ffr-,0), has
a zero average value:
To this zero mean surface, an average value,
increasing with angle 6 , was added. The result
was the desired slowing of the average
longitudinal flow, but at angles greater than
40° local backflows developed. These
backflows are not unrealistic, but are
incompatible with the formulation of the
initial value problem.
An alternative means for slowing the
longitudinal flow is to slightly decrease the
average longitudinal velocity at each
incremental angle 0 at which the velocities
are computed. ?(r,0) in Equation (18) is still
computed to evaluate 9f /90 which is needed
in Equation (8). If the inflow rate is Q6, and
the foul flow withdrawal rate is Qs , then a
flow rate equal to (Q0 - Qs) must flow over
the weir. Assuming a uniform discharge rate
over the weir over the arc length 0 m , the
discharge over the weir per radian is (Q0 —
Qs)/6m. So, at any angle 0<0<0m,the flow
rate in the channel must be:
0
(19)
However, if at any angle 0, vt (r,0(?) the
velocity is calculated from Equation (10), the
channel flow rate would be computed to be:
dz dr
(20)
The desired slowing effect will take place if
Vi(r,6 ,z) at each point is multiplied by (Qc
(9 )/Gc • (0) )• This will always be a positive
constant, very nearly equal to unity.
Calculated in this way, negative values, of v
(backflow) will not develop:
h(r)
(r,B,z) dz dr
(21)
84
-------
This approximation, therefore, allows
treatment of the flow as an initial value
problem, while still partially accounting for
the variation in longitudinal velocity, v, which
arises from the secondary fluid motion (u, iv).
Coordinate Transformation
For computational purposes, it is
convenient to deal with a rectangular array of
points in each 6 plane. If the channel cross
section were rectangular, Equations (6), (7),
(8), and (15) could be used directly. The
channel cross section is, however,
approximately a trapezoid so a
transformation is used to map the trapezoidal
cross section into a rectangular computational
space.
In Figure 72(a), Coordinate
Transformation, if at both r = r^, and r = rz,
the vertical sides are divided into the same
number of equally spaced intervals, the
intervals at r = r^ will obviously be greater
than at r = r2. But in Figure 73(b), the
intervals at R = R1 and R = R2 will be equal.
The transformation (r,d,$) -» (&©,£)
transforms the trapezoid into a rectangle. To
accomplish this transformation, define
Z = zlh(r)
0
-------
From the definition of the transformation:
ZfJl£ and Zz=^-.
h(r) h\r>
(29)
To summarize:
R =r
Z =
9 =
"57
9
3z
-2-
3_
3Z
0=0 S
9 =
30
Simplified Transformed Equations
When the transformed partial derivatives
in equations (26), (27), and (28) are applied
to Equations (6), (7), (8), and (15), the
results are as follows:
38
w
3u1
_
9Z
(30)
39 >'
h(r) 9Z
ae
(31)
w —
h(r) 9Z
.M.I 9v
"
(32)
---
Mr) 3Z
h(r)dZ
(33)
By referring to Figure 72(a) and noting that
in Equation (23) hb is the bottom slope, the
velocities u and w are related approximately
as follows:
^ = hbZ (34)
This shows that in any 0 plane, near the
bottom (Z = 1), the flow is parallel to the
bottom, and near the top (Z = 0), the vertical
velocity component approaches zero. When
this simplification is applied to Equations
(30), (31), and (33), the following equations
result in the form which is used for
computation:
30
dv =^R
30 v
w —
9v 4. uv 4. _£ 3?
9K "R R 9©
__
az R R a©
(35)
(36)
h(r)dZ
(37)
-™™,\ h(R)dZ
(38)
The channel cross section is described by
the radii TJ and r2 , height of sides HI and H2,
and overflow height and position B2 and Bl
as shown in Figure 73, Mathematical
Representation of a Channel Cross Section. If
Q° is the inflow rate, and Qs the foul flow
withdrawal flow rate, the fluid will flow over
the weir with a height ZO, given by the
standard weir formula:
-I 2/3
r Qo -Q,
zo =
3.33(r2)(0
n\
(39)
where ®m is the arc length of the flow
chamber. The coordinate system shown in
Figure 73 has its origin at the inner radius
and at the surface of the fluid. Along the
vertical surfaces, there is no radial flow
except at the overflow outlet, hence:
86
-------
r=r.
Bl + ZO
__£__r = r2
T
T
T
B2
z =H2 + ZO
z = Hl+ZO
FIGURE 73 MATHEMATICAL REPRESENTATION OF CHANNEL CROSS SECTION
u(r,e,z) - 0
r=r2, except oottet
(40)
The outlet has a height B2 and length r2@m ,
so the radial flow rate at the overflow outlet
is as follows:
and
fr
\_fBl +ZOJ < z
(43)
where:
p
k
y0
= bottom shear stress
- density
= von Karman constant
= length related to the
roughness height.
The term vrjp is the friction velocity, Vf,
which is given by:
(44)
where / is the Darcy friction factor from the
Manning equation
(45)
87
-------
In Equation (45), g is the acceleration of
gravity, n is Manning coefficient and H is the
hydraulic radius.
°lfj> Equation (43)
is equivalently written:
v = Cln (yly0) = C (In fy/h)+ln(h/y0)] ,
where h is the channel depth. For this
velocity distribution the average channel
velocity is:
Y = C j" [ln(y/h) + ln(h/y0) ] d(yjh)
)*l.} (46)
Then an equivalent expression for equation
(44) is:
v = V + C ( 1 + In (y/h) )
(47)
In this expression, V is the average channel
velocity, a function of flow rate and channel
cross-sectional area only. The second term is a
perturbation term which predicts slower
moving fluid near the bottom and faster flow
near the surface.
Now suppose that the channel is not
wide, so that the effects of the vertical walls
are felt. A functional form analogous to
Equation (43), taking into account bottom
and both walls is as follows:
v = Cln
in which a is the width of the channel.
Again, expressing the velocity in average
and perturbation terms, the result is
analogous to Equation (47):
= V + C
In ^y/h) fx/a) fij*- Tl
(49)
Finally, consider the channel cxos s
section shown on Figure 74, Channel Cross
Section in Cartesian Coordinates, in a
Cartesian reference frame.
1
*,.
, Y)
FIGURE 74 CHANNEL CROSS SECTION IN
CARTESIAN COORDINATES
Define Y as the height above the sloping
bottom:
Y = y - r^-f-k)* (50)
The appropriate form of the Prandtl-von
Karman law is as follows:
v = C In
[,
(
' '" J (5D
Again, in terms of average velocity, it can be
shown that this velocity distribution can be
written:
= V f Cln .
where the average velocity, F, is obtained by
dividing the total flow by the cross-sectional
area:
— 6
V ~ K(ni+hz )a (53)
This result, in an (x,y) plane, is easily
transformed into the r, z plane defined in
Figure 72(a) by using the following:
r = ri + x
z = hi - y
a = ra - /-!
(54)
88
-------
3. —
h(r)
f
*
(a)
FIGURE 75 TRANSFORMATION OF CHANNEL CROSS SECTION
TO RECTANGULAR REGION
Initial Values for the Cross-Flow Velocities
u and w
The initial values for the cross-flow
velocities describe the flow from the interior
of the channel over the weir at the entrance
to the bend, before the customary helical
flow pattern develops.
The u and w velocity components lie in
the © = 0 plane, and are, respectively, in the
radial (>) and vertical (z) directions. Figure
75(a), Transformation of Channel Cross
Section to Rectangular Region, shows the
channel cross section, while Figure 75(b)
shows a transformed cross section.
The ultimate goal is a velocity field in
Figure 75(a), where flow is parallel to both
the top (z = 0) and the bottom (z = h(r) ),
and no flow crosses the sides (r = r^ and r =
r2) except for a point sink at r = r2, and z = 0,
which represents the flow over the weir.
Several steps are required to reach this point:
(a) In Figure 75(b) a potential flow field
solution is derived which crosses and is
perpendicular to x =' 0 at a uniform rate on
0//=0 (57)
along with the boundary conditions:
y = 0 i// = 0
x = 0 \l> =
y - hi ty — Uoh-z
x = a \jj = U0/z2
where
(58)
(60)
(61)
(discharge/
unit length). (62)
89
-------
These boundary conditions are such that a
uniform flow crosses between y = 0 and
y = /z2 at x = 0. A point sink is located at
x - a, y = 0, while no other flow crosses y = 0,
y = h'-i, and x = a.
The solution of Equation (57), with the
above boundary conditions at points interior
to the boundary, is given by:
sink
where
2U0/z,
. . ,n7Ta
OT S/rtfc (-—
(63)
(64)
In practice, an infinite series cannot be
completely evaluated, so Equation (63)
includes all terms up to the point where:
2U0 h-> sink (rntxlh-,)
I nit sink (
<0.01
(65)
Now, suppose that the velocity
components u, and w, are rotated by a
variable angle so that u\ is parallel to the
channel top and bottom, as shown in Figure
76, Rotation of Velocity Components.
U
FIGURE 77 ROTATION OF VELOCITY
COMPONENTS
The angle 4> is chosen so that there is enough
rotation to make MI parallel to the bottom at
2 = h(r), while there is no rotation at z = 0:
0=74-, tan'1
h(r)
(66)
Then the corresponding velocities (w2,
) in the cylindrical coordinate system are:
= MI cos 0 + Wi sin 0
(67)
W2 = -u\sin 0 + M>I cos 0 (68)
These velocities must next be scaled by a
factor F, so that the same lateral discharge
exists for any value of /•, despite the variable
bottom depth h(r):
u3=F-u2 w3 = F- w2, (69)
where the scale factor, F, is chosen so that
F- [it 2 cos 0-w2 sin 0] dz
But
so
sin
U=~h(7j e
w =
Kr [-MI sin 0 +
'(76)
COS 0]
(77)
90
-------
At r = r2, the scale factor is unity. An
experimentally determined acceptable value
for the constant is Kr = 1.5.
These approximate initial values for u and
w describe the flow in the cross-flow (r, z)
plane at 0 = 0, which arises because the flow
upstream of 0 = 0 has already adjusted to the
presence of the weir. This flow generally
proceeds from the inside toward the outside of
the bend, and upward from the bottom. The
exact values of the velocities are not critical,
because the flow will re-adjust itself as the
computation proceeds downstream.
Numerical Method for Liquid Flow
The transformed space shown in Figure
72(b) is divided into a rectangular grid as
shown in Figure 77, Computational Grid in
Transformed Space. Grid points are spaced by
a distance AR radially, and AZ vertically.
Figure 77 is a plane in R and Z. A complete
flow field solution includes a number of these
planes separated by angular increments of A0-
The cross-flow velocities, u(R,®,Z) and
v(R,@,Z), are calculated from Equations (35)
and (36) by use of a predictor-corrector
scheme. In both stages of this two-stage
process, the right hand sides of these
equations must be evaluated.
R-
i-1, k
i, k-KL
1+1, k
•*>R
FIGURE 77 COMPUTATIONAL GRID IN
TRANSFORMED SPACE
The longitudinal surface slope 9f/9© is
evaluated in both stages as a backward
difference:
The. evaluation of d$/dR is the same in both
stages and is discussed in a subsequent
paragraph. Both u(R,®,Z)) and v(R,@, Z) are
calculated by the predictor-corrector process,
so it is sufficient to discuss only u(R,@,Z).
During the predictor stage, a predicted
value u(R,®,Z) is found from:
^(R, 0 + A0,Z) = u(R,Q,Z) + A0 --|^
d° (79)
Where 9u/9o is evaluated from Equation
(35) in which the terms du/dR is found as a
backward difference for^< 0.5(R1+R2):
du= u(R, 0, Z) -u(R -A.R. 0. Z)
_^R (8Q)
and a forward difference for R>0.5(Rl +
R2):
dji - u(R + A.R.0.Z) - u(R,6,Z)
dR A/? (81)
The new u(R,®,Z) is then the average of the
predicted u(R,® + A0, Z) and a correction
term:
u(R,6 + A0,Z) = I ^(^,0 + A0,Z)
(82)
Here, 3^/3 © is^ evaluated using the
predicted u and v", and the backward
and forward differences are reversed, so that
dH/dr is a forward difference for R< $ (R 1
+ RT), and a backward difference for
R>l/2 (Rl + RT). Both w(R,Z) and 3^/3^ in
Equations (37) and (38) must be calculated
after each predictor-corrector cycle, and both
contain the form:
f(Z)dZ
(83)
zi
90
(78)
This integral is approximated by the
trapezoidal rule:
Z2
/(Z)dZa
(84)
1,AZ)
91
-------
where
Zl =/l -AZ Z2 = -AZ
Within the integrand of Equation (37), 9v/90
is evaluated as a backward difference:
-v(R,0 -A0.Z)
A0
(85)
and the terms du/3R and 9u/9Zare evaluated
as centered differences.
A.R.0..Z) -
2A/?
(86)
JZ
2AZ
Liquid Flow Summary
The equations basic to fluid flow, the
momentum and continuity equations, are
presented as Equations (1) and (2). These
equations are transformed to an equivalent
set, Equations (3), (4) and (5), in cylindrical
coordinates to represent the curvature of the
channel. They are then expressed in
Equations (6), (7), and (8) in a form from
which the velocities (u,v,w), corresponding to
coordinates (r,6,z) can be calculated for a
channel of rectangular cross section.
A surface height, $(r, 6), is introduced to
represent the fact that the fluid surface is not
horizontal, but rises radially to counteract
centrifugal force. Equation (15) provides the
surface slope, 9f/9r, needed to calculate the
radial velocity u. Then fO'.fl) follows in
Equation (8) as the integral of (3f/3r) over r.
Since the problem is formulated as an
initial value problem, the channel blockage
cannot be felt. Longitudinal velocity scaling
in Equation (21) provides the necessary
average velocity reduction. A coordinate
transformation is made to account for the
non-rectangular channel cross section. The
simplified transformed equations needed to
solve for (u, v, w) appear in Equations (35) to
(38).
Boundary conditions define the channel
walls, bottom, fluid surface, and overflow. An
initial longitudinal velocity profile, v(R,Z), is
derived, based on the Prandtl-von Karman
universal-velocity-distribution law. An initial
velocity field for u and w in the plane of the
channel cross section is developed to represent
flow toward the overflow as a point sink.
Figure 78, Fluid Flow Solution
Summary, summarizes the computational
INITIAL VELOCITY v(R, 6= 0, Z)
EQ (52)
'
'
INITIAL VELOCITIES u(R, 6 = 0, Z),
w(R, 6 = 0, Z) EQS (76), (77)
3£
g T7
'
C (r, 6 )
•
EQ (38)
'
EQ (18)
^
VELOCITY COMPONENTS u AND v
BY PREDICTOR -CORRECTOR, EQS (35), (36)
'
r
u BOUNDARY CONDITIONS
EQS (40), (41)
1
w BOUNDARY CONDITION
EQ (42)
\
w (r, f), Z) EQ (37)
\
INCREMENT 6 BY A 6
NO SQ* ,
ST(
r
)m J>
YES
DP
FIGURE 78 FLUID FLOW SOLUTION
SUMMARY
procedure, and shows the equations used in
each step.
92
-------
Particle Flow Calculation
The calculation of the particle flow field
is carried out, based on the assumption that
the particle number density is low so that
particle interaction can be neglected. The
individual particle velocities at each grid point
can be accurately represented by
superimposing the particle settling velocity on
the calculated liquid velocities. The local
concentration of the particle cloud is then
calculated from the continuum continuity
equation. These calculations are performed
for a discrete particle size class, and are
repeated for various settling velocities in order
to represent the mixture of particle sizes to be
encountered in practice.
Particle Flow Equation
The solution of the particle flow problem
is carried, out as an initial value problem,
made possible by the assumption that there is
no diffusion in the longitudinal (0) direction,
and that no back-flow (negative v velocities),
occurs. Written in conservation form, the
following equation describes the steady-state
particle concentration (Q distribution:
=0 (88)
Vp is the particle velocity vector
comprised of the components (up,vpwp).
While u and v are those calculated in the
fluid flow problem, wp is the sum of the
component w from that solution, plus the
particle settling velocity, ws. Written in
cylindrical coordinates the equivalent
equation is:
9r
= e
90
N'
9C
9r
9z
72
902
9z2J
(89)
The eddy diffusivity has been assumed to
be constant. The assumption of no
longitudinal diffusion is equivalent to
assuming that (r2) 92C/902 is small with
respect to the other terms, and may be set to
zero. Having done this, and re-arranging the
factors, the equation becomes:
+
J (90)
The numerical solution basically follows the
same pattern as the solution for the liquid
flow. The right hand side of Equation (90) is
evaluated at grid points in a given 9 plane;
then, .
(vQr,6+&e,z=(vCr,8,z+ ^ \~W~ } r>d>z •
V '
and
(91)
A0,z) (92)
Transformed Equation
As is the case for the fluid flow solution,
the numerical solution of the concentration
equation is made easier by a coordinate
transformation. The transformation is
repeated here:
9 = JL + 7 _JL
R =r
br
0 = & ^- =
9
90
_ 9
90
Z =
where
9 -7
9F ~zz
(93)
(94)
(95)
(96)
(97)
Applying this transformation to Equation
(90):
d(vQ = _ I 9 (RUG) + z
-Zhb
h(r)
and
Z = 1
h(r)
-\-*-
dR
r dZ
9Z9.R
93
-------
+ Z,JK.
__
9Z
M
+ R
9Z
(98)
Boundary Conditions For Particle Flow
The boundary conditions for the particle
concentration field, (Equation 88), require
specification of the particle flux vector F
defined as:
F =
eVC
(99)
If ~n denotes the unit inward normal to one of
the bounding surfaces, then the particle flux
normal to this boundary is:
= n • F =
n - e
dC
dn
(100)
Along the vertical walls (except at the
overflow), the particle flux, normal to the
wall must be zero. Since the velocity normal
to the wall is zero, Vp.n = o, for these surfaces,
so that Equation (100) reduces to:
-|^- = 0 (for vertical surfaces). (101)
Through the open surface of the channel
there is also no particle flux. However, along
this free surface, the vertical velocity of the
liquid, w, is zero, but the particle velocity is
the settling velocity, ws . The rate at which
particles settle from the surface must,
therefore, balance the flux of particles toward
the surface by diffusion, or, from Equation
(100):
= 0 (open surface)
(102)
Along the bottom, it is possible to specify
either of two boundary conditions. The
particle flux through the bottom can be
assumed to be zero (F • "n = 0) in which case
all of the bed load remains in suspension
Alternatively, particles reaching the bottom
can be assumed to settle into a thin, highly
concentrated sediment layer external to the
computational mesh. The first boundary
condition — all bed load remains in
suspension — is probably somewhat closer to
reality and has been used in the present
model. However, some deposition of sediment
is known to occur and, by neglecting this, the
model will tend to under-predict the removal
efficiency. To correct for this tendency, the
eddy viscosity has been empirically calibrated
against laboratory data, as will be described to
reproduce the correct removal efficiency at
the design flow rate.
Initial Concentration Distribution
An analytical derivation of the
concentration in the 0=0 plane is not
as straightforward as the initial velocity
calculation. The actual particle concentrations
which exist at the entrance to the bend are
influenced by the design of the entrance
section. It was originally thought that the
vertical variations in concentration could be
obtained from the analytical solution for
sediment transport in a straight open channel.
However, the laboratory measurements
revealed significant lateral concentration
variations which could not be modeled with
this approach. For this reason, the initial
concentration distribution is based on
measured data at 12 points in the channel
cross section. Initial concentrations at each
computational grid point are calculated as a
linear two-dimensional interpolation of the
four measured values surrounding each grid
point. For grid points not surrounded by four
data points, the interpolation becomes an
extrapolation.
Numerical Method for Particle Flow
The numerical solution for C(.R,®,Z)
proceeds ahead in angle 0, according to
Equations (91) and (92). These require that
the right hand side of Equation (98) be
evaluated at grid points in R and Z at each ©
step. The points (o) at which C is calculated
lie midway between the points (x) at which u,
v, and w are calculated in the (R,Z) plane.
Figure 79, Location of Computational Grid
Points for Concentration Field Calculations,
shows the arrangement.
The boundary conditions for the particle flow
are built into the numerical procedure by
specifying the appropriate particle fluxes to
be zero at the solid surfaces and along the
open top.
94
-------
\
r
R
k
Z i
(i-1, k-1) (i-l, k) (i-1, k+1)
(i-1, k-1)
o
(i-1, k)
o
k-1)
(i, k+1)
(i, k)
O
AZ i o
1(1+1, k-1) (i+1, k) (i+1, k+1)
AR
FIGURE 19 LOCATION OF COMPUTATIONAL GRID POINTS FOR
FOR CONCENTRATION FIELD CALCULATION
The numerical calculation of each of the
eight terms in Equation (98) is described
below. In the first partial derivatives in the
first three terms, care is taken so that mass
flux is conserved.
Term 1
b(RuO = (Ru)n C0 - (Ru)m Cm
9R AK (103)
where
(Ru)p =
(Ru)m = >
i,k
and
= Citkfor(Ru)p>0
Cm = Ci,k-l f°r (Ru1
Cm = Ci,k f°r (^«)m<°
R(k) and Z(i), which will appear later, are
evaluated at the points (x).
Term 2
3( RuC) -
(ZrRu) =
where
(Z^?u)m =
and
Cp =
(104)
1, k
>0
Zr (i,k) and Zz(k) in the following
expression are evaluated at the points (o).
95
-------
TermS
Term 6
= 2
2A z
(105)
where
_ 1
2.
_ 1
(108)
wi
it k
and
Term 7
By substitution forZ,.,
~ ct. k for wp > o
Cm =
, k for
for wm
(109)
0
The diffusion terms, 4 through 8, all span
two increments in A/? or AZ as written in the
expressions to follow. However, when the
points (0) are adjacent to the channel
boundaries, only single spans in &.R or AZ
may be possible, in which case one-sided
ditrerences are used.
Term 4
- q fc) /AZ - z, (qfc - q.1;fc) / Az /A z
J
Term 8
Again, by substituting for Zz,
3C\ _ R 32C
(110)
RZ,
3Z
3Z2
(h(Rk)
(106)
+ R 1 /
^fc+i /. _j
+ h(R ))2
fc-t-i |—
AZ
TermS
9
Az
/AZ..
(in)
*-. i
/2AZ 1/2A.R
Finally, when integrating C(R,@,Z) ahead
in 0 by use of equations (91) and (92), the
velocity term, v, is evaluated as the average of
the four points surrounding each
concentration mesh point.
Particle Flux Calculations
The particle flux crossing the 0 = 0 plane
(107) is the integral of the product of concentration
96
-------
and longitudinal velocity over the cross
section:
v(r,e=o,z)dzdr
(112)
C(r,6 =o,$) dzdr
The flux over the weir is a function of the
radial velocity and concentration at that
point. The flux between angles ©i and ©2 is:
ft^ Bl +.B2 + Zo
(C(HJ J»M^C^
®n Bl + Zo
dzdd.
Then, the total flux over the weir is:
(CQW
W otal
= 6
(1 14)
The removal efficiency is, then, the ratio of
the material which did not pass over the weir
to the input flux.
-(.CQw)total
CQ0
(115)
Particle Flow Summary
The particle flow solution is treated as an
initial value problem, starting with
interpolated experimental data in the 0=0
plane. After each A0 increment where
velocity components u, v, and w are
calculated, the concentration is calculated at
points on a (r,z) plane grid. Equation (88) is
the basic relation describing the concentration
field. It is solved numerically in conservation
form for the product (vC) at each angle 6;
then, by dividing by v, the concentration C
results. Equations (90) to (92) describe this
process. In order to account for the
non-rectangular channel cross section, a
transformed equation for 9 (vC)/d6 is given in
Equation (98).
The numerical solution method is a
conservative upwind differencing scheme.
The numerical calculation of the
individual terms of Equation (98) is
described in Equations (103) to (111).
Finally, the particle flux calculations
needed to determine efficiency are presented
in Equations (112) to (115) Figure 80,
Particle Flow Solution Summary, shows the
overall flow and concentration calculation
procedure. Figure 78 portrays details of the
flow calculation.
INITIAL VELOCITIES u, v, w
\
r
INITIAL CONCENTRATION,
INTERPOLATED MEASURED DATA
\
r
INPUT FLUX CALCULATION
EQ. (112)
\
i
VELOCITY COMPONENTS u, v, w
CALCULATED (FIGURE (11) )
'
3 (vC)
se
t
<
EQ (48)
r
c (R, e, z)
EQS (91, 92)
i
>
INCREMENTAL OUTPUT
FLUX, EQ (113)
r
INCREMENT 6
BY Aft
N0 ^n -
flm ^>
|YES
EFFICIENCY
EQ (115)
t
STOP
FIGURE 80 PARTICLE FLOW SOLUTION
SUMMARY
97
-------
Approximate Solution for
One-Dimensional Flow
In a uniform, wide rectangular channel,
the suspended sediment tends to stratify, with
a higher concentration occurring near the
bottom, and lower concentration at the
surface. The degree of stratification is
influenced by the particle settling velocity
and the turbulence in the channel. The classic
analysis of bed-load stratification was
published by Dobbins6 in 1944. Dobbins
studied a simplified version of Equation (88)
in Cartesian coordinates in which the
longitudinal velocity, v, was assumed to be
constant, and u and w for the liquid were
taken as zero. As in the present study, the
longitudinal diffusion was neglected in
comparison to vertical diffusion, and the
particle velocity was approximated by
superimposing the settling velocity, ws, on the
liquid flow. In the rectangular coordinate
system shown in Figure 74 the simplified
equation becomes:
= e_9^ +
v dc
s by
(116)
Because the longitudinal velocity, V, is
assumed constant, Equation (116) can be
considered as a time-dependent problem in
which the time, t, is the liquid flow time
t =-4-
V
(117)
Dobbins obtained analytic solutions to
Equation (116) for the boundary condition of
no particle flux through the open surface at y
= h. Expressed mathematically, this boundary
condition is
'bC_
y=h.
(118)
in analogy with Equation (102). At the
channel bottom, Dobbins considered the
boundary condition:
•w.
A,
(119)
in which A is an arbitrary constant
representing the asymptotic bottom
concentration at t ->°°. For A=0, the
asymptotic concentration is zero, indicating
that all of the material has settled out of
suspension. The initial concentration
distribution was assumed to be of the
exponential form:
Cty.o) =C0e -by .
(120)
For the special case 5 = 0, Equation
(120) corresponds to an initially uniform
concentration. For this problem, Dobbins
obtained the analytical solution:
CG/,0 = Ae
Cne
n sin«B;T| ,
" J
(121)
where the constants, Cn, are evaluated from
2
in which
(122)
H = cosah +
W.,
sin oLnh .
(123)
The values of «„ are obtained as the solution
to the transcendental equation:
>wsh \
2 cot OL./Z =
2e
WL\ anh
ze
(124)
98
-------
For the boundary conditions studied by
Dobbins, it is necessary to specify the
asymptotic bottom concentration, A. In
practice, unless complete clearing of the flow
is anticipated (so that A = 0), it is difficult to
know what an appropriate value for A might
be. Consequently, Dobbins' solution was
rederived for the bottom boundary condition
of zero mass flux through the, bottom, or:
AC = r
3y y=Q L
•=o (125)
in analogy with Equation (102). This
boundary condition corresponds to that used
in the numerical model, and assumes that the
entire bedload remains in suspension. Under
this condition, it can be shown that the
solution to Equation (116) is:
-ws y -ws y
C(y,t) = Ae
2
46 2
_
n=\
where
A =
\\-e -Bh
' ~e~
' (126)
(127)
(128)
,and
Cn =-
(-ire
2e
2e
(129)
This new solution is closely analogous,
but somewhat simpler than that of Dobbins.
Furthermore, the value of the previously
unknown constant A, is now determined from
the requirement that the longitudinal mass
flux remain unchanged — i.e., all bed load
remains in suspension. Although these
solutions correspond to much simpler
geometry than the helical flow separator,
several useful insights into the more
complicated problem can be obtained from a
study of their analytical form. For example,
in both Equations (122) and (126), the
asymptotic vertical concentration distribution
is shown to be a negative exponential:
C(y, t-»»>)=A exp (-wsy/e), (J3Q)
which depends on the scale height "a"
defined by
a = e/ws. (131)
Thus, sediments with large settling
velocities, ws, will have small scale heights,
corresponding to a very rapid decay, in
concentration with distance above the
bottom. Furthermore, for a given settling
velocity, the scale height is directly
proportional to the eddy viscosity due to
turbulence, e. Thus, as the turbulence level is
increased, the vertical concentration
distribution tends to become more uniform,
with less stratification. In the discussion by
Camp''7'' of Dobbins' paper, the eddy
viscosity was evaluated as the vertical
averaged value corresponding to the universal
logarathmic velocity distribution as discussed
in Equation (43). The resulting expression for
e is:
€ = Thvv{- (132)
where h is the channel depth, Kis the average
velocity, and / is the Darcy friction factor.
This useful result has been retained in the
present analysis. The eddy viscosity increases
directly with the average velocity. Since the
recovery efficiency of the helical flow
separator depends on the degree of bedload
stratification, Equations (130) and (132)
indicate the desirability of reducing -the
average velocity as much as possible by using
a large channel cross section.
It was originally 'thought that the
asymptotic vertical concentration distribution
given by Equation (130) could be used
directly to model the initial distribution at
the entrance to the bend, by appropriate
correlation of the scale height, "a". However,
as noted, the laboratory results showed
99
-------
significant lateral variations which could not
be accounted for with this approach.
Additional useful insights can be obtained by
studying the relationship between time and
distance given by Equation (117). For a
channel with a side weir, the average
longitudinal velocity decreases in the
downstream direction, due to loss of fluid
over the weir. For uniform lateral discharge,
the average velocity is given by:
V = dx
dt
in which Q0
_ Qo -Qw
A "~ (133)
is the flow rate at the start of the
side weir, Qw is the total flow over a weir of
length L, and A is the cross-sectional area of
the channel. Equation (133) can be integrated
to give the following relation between
distance and time:
Q
?0-Q,
(134)
Since the degree of stratification (and hence
separator performance) which can develop
within the separator depends on time (as
shown in Equation 121 or 126), Equation
(134) indicates that better performance can
be achieved with longer retention times, T,
defined in the conventional manner:
. = vol „ LA
Q Qw ' (135)
By combining this result with the
analytical form appearing in Equations (121)
and (126), it can be concluded that a
significant non-dimensional parameter, 0,
affecting the separator performance is:
9 i~~ (136)
where ws is the settling velocity, T is the
retention time defined by Equation (135),
and e is the eddy viscosity. To increase
separator efficiency, 0 should be made as
large as possible, which implies large retention
times, T, and small longitudinal velocities, to
reduce e . Furthermore, since the
concentration decays exponentially with
distance from the bottom, and the clear flow
is withdrawn from the surface, it appears
desirable to achieve a given retention time
with a narrow, deep section, as opposed to a
shallow, wide one. This conclusion strictly
applies only to a straight channel. In the
proposed curved channel, the induced
secondary flows act to reduce the
concentration of particulates in the overflow.
The development of these secondary flows
would be suppressed with a deep, narrow
channel section.
Scaling of the Liquid Flow
The development of the secondary flow
in the bend depends primarily on the
distribution of the approach velocity at the
entrance of the bend. If it can be shown that
the non-dimensionalized inlet velocity profile
is the same in both the model and the
prototype, then maintaining the same Froude
number,
F=
-
Lg
(137)
in the model and the prototype will ensure
that the flow patterns will be the same
throughout the bend.
It was shown that the entrance velocities
for the non-rectangular cross section can be
represented by:
+ 1/2
h, +/z2
In
(138)
where "a" is the width of the separator, hi
and hi are the depths at the inner and outer
walls, and Y is the distance from the bottom.
This formula applies to both the prototype
and the model and, since all terms inside the
brackets are non-dimensional, each of these
terms will be identical if the model and
prototype are geometrically similar.
Consequently, the .same non-dimensional
values of velocity (v/V) will result for both
model and prototype if the term C/V is the
same for model and prototype.
The constant Cis given by:
C=l/k Vf, (139)
100
-------
where k is von Karman's constant, and Vf is
the friction velocity
v, (140)
and where / is the Darcv friction factor.
Therefore, the ratio ofC/Fwill. be the same in
model and prototype if:
(vTT)m = (V/V (Hi)
The friction factor, / can be determined
from the Manning Equation:
f= 0.453
#1/3
(142)
in which g is the acceleration of gravity
(ft2 /sec) n is the Manning coefficient (sec
ft-7/5) and H is the hydraulic radius. Typical
values of n are 0.009 for the model and 0.012
for the concrete prototype. Thus the ratio of
friction factors is:
Jttt f **w\ \ / •"• n \ ^/*^
~fp~ ~\^j (Hr») (143)
For a length scale of 1:6:
fm = (0.009\ 2 <-6W3 = i 025
-Jp— lo.012J (6) 1"025'
and so the ratio of CfV in model and
prototype will be:
=1.013,
(em.
which is very close to unity.
It can, thus, be seen that the lower model
velocities are offset by its smoother surface,
so that to a very close approximation the
non-dimensional entrance velocities, vf'y will
be the same as for the prototype.
The entire flow in the bend will then be
reproduced by maintaining the Froude
number constant between model and
prototype. Thus, if the prototype dimensions
are S times larger than the model, the
following relationship follow from Equation
(137):
Lp=SLm (144)
Vm =
v
(145)
(146)
Scaling of the Particle Flow
Equation (89) for the particle
concentration field applies to both the model
and the prototype. If the equation is
non-dimensionalized by '-'dividing-, the
concentration by the inlet concentration, C0,
dividing all lengths by the width of the
channel, "a", and dividing all velocities by the
average entrance velocity, V, there results:
:i_ _!_ (r 3£ ) + -L
- V
3ac
1
_ lev
,2
(147)
in which the tilde (-) denotes a non-dimen-
sional variable. This equation will be identical
for the model and. the prototype if the
non-dimensional eddy diffusivities are the
same:
ef
-
(148)
The equation used to represent the eddy
diffusivity is Equation (132):
e = -1 h V ^/Ji* ,
6
where Ji is the mean depth. Equation (148),
then, requires that:
h V2
VP)
p
(149)
As was shown previously, fm& fp, and since
the velocities and dimensions are related by
Froude number scaling (Equations 144-146),
p
.2
= s2
= s V
m '
so that the equality in Equation (149) is
preserved.
101
-------
As a consequence, Froude number scale
relationships between the model and
prototype will ensure duplication of the
non-dimensional concentration fields and,
therefore, the recovery efficiency. The
Froude scale relationship, Vp = vT~ Vm,
applies to all velocities, including the particle
settling velocity. Thus, to represent a given
class of prototype particles in the laboratory
model, the laboratory material must have a
settling velocity of S— 'A times the prototype
value, where S is the length scale.
Comparison of Mathematical Model
Results with Test Data
Figures 81-83, Comparison of Measured
and Predicted Velocity Profiles at 0°, 20°,
and 40° respectively, show comparisons of
calculated and measured longitudinal, (v),
velocities in cross sections at the indicated
angular positions. The test data were obtained
by LaSalle Hydraulic Laboratory in a
companion study, as part of this project. The
velocities shown are in ft/sec, prototype scale.
In Figure 81, the calculated velocity
contours are in good agreement with the
laboratory data. The predicted velocities at
this position were calculated from Equation
(52) applied to the prototype, for 0.85
m3/sec (30 cfs) inlet flow, with a Manning
coefficient n = 0.012 in Equation (45). As
was shown, the non-dimensional velocity
profile, (y/V), will be very nearly identical for
both the model and the prototype, hence all
calculation has been performed at prototype
scale.
Figure 82 shows significant shifts in the
longitudinal velocity at the 20° position. Due
to the secondary fluid motion, low energy
boundary layer fluid from the bottom is
swept up to the surface at the inner radius of
the bend, resulting in a low velocity region
there. This is clearly seen in the laboratory
data, but is not reproduced in the computer
simulation. This' deficiency in the
computational technique was also noted in
the results of Fox and Ball,4 and is a result
of incorrectly modeling the convective
acceleration term w3v/3z in Equation (34). In
the work of Fox and Ball this term was
neglected altogether. In the present model,
the ratio of w and u was assumed to be of the
form given Equation (34), which is correct
along the bottom and the surface, but does
not properly account for the upflow which
occurs at the inner radius of the bend. An
attempt was made to. resolve the discrepancy
by retaining all of the terms shown in
Equation (31). However, with the additional
terms, the computational technique became
unstable.
In Figure 83, the laboratory data shows a
separated (reverse) flow region near the inner
bend. This is evidenced by the zero velocities
in Figure 83 (a). Dye studies reveal very small
negative velocities in this region. The reverse
flow region is the consequence of the low
energy fluid being more rapidly decelerated
than the faster moving fluid. The original
technique used for calculating the
longitudinal surface slope did result in the
appearance of a reverse flow region in the
computer simulation. However, the predicted
reverse flow occurred further downstream
because the calculated longitudinal velocities
near the inner bend were greater than in the
model, due to the w9v/9f term discussed
above, and thus, took longer to decelerate.
In any event, the reverse flow region is
incompatible with the initial value problem
formulation and had to be eliminated in order
to proceed. The resulting discrepancies in the
longitudinal velocity field are not particularly
important in computing the concentraion
field. The differences in the longitudinal
velocity do result in further errors in the
cross-flow velocities u and w, which are
potentially more critical. However, the
principal effect of the reverse flow region is to
raise the longitudinal velocity near the outer
radius of the bend and to cause an eddy to
appear just upstream of the baffle at 60°. This
eddy causes a sudden uplift in the material
discharged over the weir near the end of the
bend, but the effect is very local.
Figure 84, Comparisons of Measured and
Predicted Flow Angles on Bottom Along
Channel Centerline, shows a comparison of
the measured and predicted flow angles along
the bottom centerline of the bed, with and
without flow over the weir. To obtain the
102
-------
(a) Measured Velocities
(b) Predicted Velocities
NOTE: Velocities are relative
FIGURE 81 COMPARISON OF MEASURED AND PREDICTED
VELOCITY PROFILES AT 0°
103
-------
I
(a) Measured Velocities
NOTE: Velocities are relative
(b) Predicted Velocities
FIGURE 82 COMPARISON OF MEASURED AND PREDICTED
VELOCITY PROFILES AT 20°
104
-------
0.0 0.1 0.2 0.30.4 0.45. 0.50
41
(a) Measured Velocities
NOTE: Velocities are relative
(b) Predicted Velocities
FIGURE 83 COMPARISON OF MEASURED AND PREDICTED
VELOCITY PROFILES AT 40°
105
-------
Bed Angle (degrees)
a LAB DATA, NO OVERFLOW
O LAB DATA, OVERFLOW
WITH SCUMBOARD
^CALCULATED, NO OVERFLOW
.^-CALCULATED, OVERFLOW
10
35 0 (degrees)
FIGURE 84 COMPARISON OF MEASURED
AND PREDICTED FLOW
ANGLES ON BOTTOM ALONG
CHANNEL CENTERLINE
laboratory results with no overflow, the
partition at 60° was removed so the flow
proceeded uninterrrupted for a full 120°.
With no overflow, the laboratory data
indicate a bed flow angle near 22° along most
of the bed, with 18° bed flow angle at the
10°position. This implies a very rapid initial
growth in the bed flow angle. In contrast, the
computer simulation indicates a continuous
growth in bed flow angle, reaching 22° only
at the 35° position. In this respect the present
results are in agreement with those of Fox
and Ball.4 They showed a similar steady
growth in bed angle in substantial agreement
with their test data, although their work was
for a flat bed.
The measurement of the bed flow angles
is performed by using Process Black ink,
placed along a radius at a given angular
position. The fluid shear along the bottom
causes threads of ink to extend downstream
from the inked line, and the angle of these
threads with respect to the Q direction are
measured with a protractor. The technique is
accurate only to (±3°), but the differences in
in the no overflow results shown in Figure 84
are too large to be explained as measurement
inaccuracies. The very rapid initial growth in
bed angle in the laboratory model is puzzling.
There is no obvious reason why the
laboratory and mathematical models should
not be in better agreement in this respect.
The results for the bed angle in the
presence of overflow are in much better
agreement. For this case, the laboratory
results indicate a slowly growing bed angle.
The computer simulation matches the
laboratory results at the 10° position and
then somewhat overpredicts the growth in
bed angle. However, the differences are less
than 5°. Both the laboratory model and the
-computer simulation indicate substantially
smaller bed flow angles in the presence of the
discharge over the weir. This indicates that
the helical secondary flow is largely destroyed
by the weir and scum board. The reason for
this is that the side weir at the outer bend
effectively controls the water level there and
prevents the development of a large radial
surface slope. The smaller radial surface slope
produces a smaller helical motion with smaller
bed flow angles, as observed.
Comparison of Calculated Concentration Field
with Laboratory Data
It was originally anticipated that the
simple straight channel theory could be
applied to the inlet transition section to
obtain the initial concentration field at the 0°
position. However, the laboratory data
revealed significant lateral variations in
concentration which could not be adequately
accounted for with the simplified theory.
Consequently, the initial conditions for the
mathematical model were obtained by
interpolating the measured concentration data
at the 0° position, shown in Figure 85,
Measured Values of Relative Petrothene
Concentration at 0°. The measured
concentrations are only applicable for
shredded Petrothene at the nominal design
discharge of 0.85 m3/sec, (30 cfs). At higher
flow rates or lower settling velocities, the
concentration distribution will be more
uniform in the vertical direction. Conversely,
at lower flows and higher settling rates, the
bed load will be more stratified. To account
for these differences, a correction was applied
to the laboratory data, based on the theory.
The steady state vertical concentration
distribution is of the form:
C_ =e
-(wsy)
e
(150)
106
-------
FIGURE 85 MEASURED VALUES OF RELATIVE
PETROTHENE CONCENTRATION AT 0°
where y is the distance above the channel
bottom, ws is the settling velocity, and e is
the eddy diffusivity. Consequently, at flow
rates or eddy viscosities, and settling velocities
other than the nominal ones applying to the
laboratory measurements (ws0, and eo), the
initial concentration data were taken as:
C(r,z) = C0 (r,z)
(151)
where Q, (r,z) is the measured concentration
field and h is the average depth. It is readily
apparent that if ws/e = s0/e0\i no alteration
results, whereas for a larger settling velocity
or smaller flow rate, the concentration near
the surface, z = 0, will be reduced, giving a
more stratified distribution. Conversely, at
lower settling velocities or higher flow rates,
the, concentrations near the surface are
increased, giving a more uniform distribution.
Figures 86 and 87, Comparison of
Measured Values of Relative Petrothene
Concentration at 20° and at 40°, show the
measured and predicted concentration fields
for the nominal case of shredded Petrothene
(prototype settling velocity of 1.22 cm/sec
(0.04 ft/sec) at the design discharge of
0.85m3/sec, (30 cfs). Figure 86 shows the
comparison at the 20° position. The measured
and predicted results are similar, having
essentially horizontal isopleths, but the
calculated results indicate a more highly
stratified bed load. The very high predicted
concentrations along the floor of the unit are
the result of assuming that the entire bed load
remains in suspension, where the laboratory
measurements do not reflect the sediment on
the : bottom. These differences can be
neglected. However, the calculated results also
indicate very low concentration near the
surface, whereas the laboratory results show
significant amounts of material there. This
discrepancy is caused by the large spread in
107
-------
(a) Measured Concentration
(b) Predicted Concentration
FIGURE 86 COMPARISON OF MEASURED VALUES OF
RELATIVE PETROTHENE CONCENTRATION
AT 20°
108
-------
(a) Measured Concentration
10.0
(b) Predicted Concentration
FIGURE 87 COMPARISON OF MEASURED VALUES OF
RELATIVE PETROTHENE CONCENTRATION
AT 40°
109
-------
the settling velocities of the Petrothene. The
simulation results were obtained for a single
settling velocity applicable to the mid-range
of the Petrothene size distribution. The
laboratory results indicate that a significant
fraction of the slower settling particles
remains near the surface.
Figure 87 shows a similar comparison at
the 40° position. Both the laboratory and the
simulation results indicate very little material
at the surface, although the laboratory results
still indicate a less stratified bed load, due to
the distribution of settling velocities. The
most important comparison between the
laboratory results and the computer
simulation is the relative rate at which
material is discharged over the weir. Figure
88, Mass Flux Over Weir as a Function of
Angular Position, shows the percent of the
original particles lost over the weir in each of
twenty 3° angular segments. Figure 89,
Cumulative Mass Flux over Weir as a Function
of Angular Position, shows the cumulative
percentage lost over the weir as a function of
angular position. The mathematical model
was calibrated to give the correct total
percentage material lost — 9 percent — for
this nominal case by adjusting the eddy
diffusivity constant, k, in Equation (132).
The selected value was 1.5, compared to a
theoretical value of 0.4.
Figure 88 shows that in the laboratory
model, the rate of material passing over the
weir is initially low, but increases rapidly to a
maximum at the 18° angular position.
Thereafter, the rate of material loss decreases
slowly. The initial transient rise in the
laboratory results is not duplicated in the
computer simulation, which indicates,
instead, a rapid continuous decrease from a
maximum at the 0° position. This significant
difference in behavior is a consequence of
incorrectly modeling the lateral flow over the
weir. In the computer simulation, the lateral
weir flow is assumed to be fully developed by
the 0° position, and is accounted for by the
initial distribution of the cross-flow velocities,
u and w. As the flow enters the bend, the net
lateral discharge, q, at each angular station is
unchanged from the initial distribution, as a
consequence of assuming 9^/90 = 0 in
Equation (12). In fact, the laboratory results
in Figure 88 indicate that the lateral discharge
is not fully developed until near the 19°
position.
The development of the lateral discharges
is controlled by the weir which reduces the
surface level at the outer bend radius, and
thereby accelerates the fluid laterally toward
the weir. Because of its inertia, a finite time is
required for the lateral discharge to fully
develop across the cross section, in contrast to
the instantaneous flow adjustment assumed in
the computer simulation. The reduction in
surface elevation at the outer bend also
inhibits the initial development of the helical
motion.
As a consequence, the bed flow angle
initially grows at a lower rate than indicated
by the computer simulation, as shown
previously in Figure 84. Several attempts were
made in the computer studies to resolve these
differences in flow behavior by properly
simulating the development of the lateral
discharge, q(r,@),along the channel. To do
this, the lateral surface slope must be
computed directly, instead of relying on
Equation (38) which is derived on the
incorrect assumption that 9/90 =0. In
principle, the lateral surface slope, 9f/9r is
given by the ratio of the u and w velocity
components at the surface:
3L= -fJL
dr
(152)
However, attempts to use this formulation
resulted in an unstable calculation. It was
concluded that proper calculation of the
surface slope is incompatible with the initial
value problem formulation. This is to be
expected, because in subcritical flow,
disturbances can propagate upstream. Such
flows are, therefore, not properly posed initial
value problems. To calculate the development
of the lateral discharge would, therefore,
require changing the entire formulation of the
problem. In principle, this could be
accomplished by using a relaxation technique.
However, the fully three-dimensional
character of the flow would entail use of a
two-component stream function, and
impractical increases in computer storage and
run time. Instead, the present model
formulation was retained, and the eddy
110
-------
0 6 12 18 24 30 36 42 48 54 60
ANGULAR POSITION (DEC.)
FIGURE 88 MASS FLUX OVER WEIR AS A FUNCTION OF ANGULAR POSITION
10
I
Hi
PS
£
SIMULATION
I
12 18 24 30 36 42
ANGULAR POSITION (DEC)
48
54 60
FIGURE 89 CUMULATIVE MASS FLUX OVER WEIR AS A FUNCTION
OF ANGULAR POSITION
111
-------
% RECOVERY
100
80
60
40
20
0
w = 0.04 FT/SEC
s
O MODEL
DATA
Q-,
1.5Q
D
2.0Q
D
FIGURE 90 COMPARISON OF PREDICTED AND MEASURED RECOVERY
EFFICIENCIES AT DISCHARGES GREATER THAN THE
DESIGN VALUE 0.85 M3 /SEC (QD = 30 cfs)
diffusivity was adjusted to give the correct
total amount of material discharge over the
weir for the nominal case, as indicated in
Figure 89.
Figure 90, Comparison of Predicted and
Measured Recovery Efficiencies at Discharges
Greater than the Design Value QD = 0.85
m3/sec, (30 cfs) shows the predicted and
measured recovery efficiencies at flow rates of
1.5 and two times the design flow rate of 0.85
m3/sec (30 cfs). It is apparent that the
computer simulation considerably
underestimates the drop in recovery at the
higher flows.
This is evidently the result of the poor
representation of the lateral discharge, which
becomes more critical as the flow rate is
increased. At these high flows, large values of
the radial velocity, u, are induced in the
central region of the cross section, thereby
sweeping a substantial amount of material
over the weir. The present method for
assigning the weir-induced cross-flow
velocities, u and w, does not adequately
account for these large cross-flows. Instead,
most of the weir flow is withdrawn from the
upper right-hand area of the cross section,
where the concentration is lower, resulting in
an under-estimate of the lost material.
At flow rates lower than the design
discharge, the computer simulation should be
substantially more accurate because the
weir-induced cross-flows are then relatively
less significant than the helical motion.
Calculated Results for Various Conditions
The computer model was exercised over a
range of settling velocities and flow rates, to
produce the parametric curves shown in
Figure 91, Predicted Separation Efficiency
Versus Settling Velocity at Several Flow
Rates. These results show improvements in
112
-------
100
95
90
85
80
14 M'
J I I I I I
0.1
1.0
SETTLING VELOCITY , w (CM/SEC)
s
FIGURE 91 PREDICTED SEPARATION EFFICIENCY VERSUS
SETTLING VELOCITY AT SEVERAL FLOW RATES
the recovery efficiency as the settling velocity
is increased and the flow rate is reduced.
Results for flow rates larger than the design
discharge of 0.85 m3/sec (30 cfs) are not
shown because the model does not give
reliable results at large flow rates, as explained
in the preceding section.
An attempt was made to reduce the
parametric curves in Figure 91 to a
single-performance curve, through the use of
an appropriate correlation. As has been noted,
a significant parameter derived from the
analysis for straight channels is the
non-dimensional function, , defined in
Equation (136):
in which r is the retention time.
LAS
T ~ -— , •
where
L is the arc length, As is the surface area, and
Qw is the weir flow. Substituting for e from
Equation (132), and noting that Qw = V.AC
where Ac is the cross-sectional area:
. w L As
(T)
1
hv^f/B (153)
For a fixed geometry, Equation (153)
suggests that the performance should vary
with the ratio of ws/V, which is the same
behavior predicted from the non-dimensional
settling parameter:
(154;
Equation (154) was used successfully to
correlate the results of the swirl concentrator
tests.8 This correlation also works for the
present case, but only if the non-dimensional
113
-------
inlet concentration distribution is used. Under
these conditions, the downstream
development of the vertical concentration
distribution — and hence the recovery
efficiency - is basically a balance between the
settling velocity and eddy diffusivity, as
implied in the use of Equation (153) or (154).
However, because of the long inlet transition
section, significant changes in the initial
concentration profile are expected as the inlet
velocity or the settling velocity is varied.
When these changes are incorporated,
correlation of the results with either Equation
(153) or (154) no longer results in a single
universal curve.
A further important point is that in the
swirl concentrator study,8 Equation (154)
was shown to represent the upper limit on
recovery efficiency, e, so that:
e<-
Q
(155)
In the present case, this result is not
applicable, and efficiencies greater than
indicated by Equation (155) can result,
especially at low flow rates and high settling
velocities. The reason for this is that the inlet
transition section also functions as a settling
chamber under these conditions.
Additional numerical experiments were
conducted to determine the importance of
the bend radius on the performance of the
helical separator. Using the same cross
section, the outer radius of 17.4 m (57 ft) was
increased by factors of two and ten. For these
runs, the weir length was maintained constant
by reducing the total angle of the bend, ®m,
to 30° and 6° for these two cases,
respectively. The results of these calculations
for the design flow rate of 0.85m3/sec (31.1
cfs) are summarized in Table 11, Recovery
Efficiency at Various Bend Radii.
TABLE 11
RECOVERY EFFICIENCY AT
VARIOUS BEND RADII
Radius Recovery
(Initial Radius) Efficiency
1 91%
2 89%
3 82%
These results show significant decreases in
performance efficiency as the bend radius
becomes larger. The decrease in recovery is a
direct result of the reduced helical motion
present in these configurations, since the
other parameters were maintained constant.
Design techniques have been summarized
in Section III.
NOMENCLATURE
A
a
B
C
e
F
F
/
g
H
Area; also constant used in
one-dimensional flow solution
Channel width
Constant used in one-dimensional flow
solution
Depth of scumboard below weir (Fig.
5)
Overflow height (Fig. 5)
Particle concentration
Efficiency
Froude number
Particle flux vector
Darcy friction factor
Acceleration due" to gravity
Head loss through orifice
h(r) intercept (Eq. 22)
h(r) slope (Eq. 23)
Channel maximum depth
/z2 Channel .minimum depth
h(r) Channel depth as a function of radius
k von Karman constant
L Weir length
n Manning coefficient
Irt Vector normal to channel wall
P Pressure
q Lateral discharge
Q0 Inflow rate
Qs Sludge withdrawal rate
R Transformed coordinate r
r Cylindrical radius coordinate
TI Channel inner radius
r2 Channel outer radius
S Scale factor
t Time
U<> Flow velocity at overflow
u Radial velocity component
114
-------
*V Velocity vector
v Longitudinal velocity component
P/ Friction velocity
Vp Particle velocity vector
w Vertical velocity component
we Equivalent settling velocity
ws Settling velocity
Y Vertical height above channel bottom
u0 Scale length related to the roughness
height
z Cylindrical vertical coordinate
ZO Fluid depth over the weir
Zr Zjr
Zz Z/z
AJR Incrementing
AZ . Increment in Z
A© Increment in 0
f Fluid surface elevation
© Transformed coordinate ©
e Cylindrical angular coordinate
@ra Channel angle maximum
p Density
T Retention time
P0 Bottom shear stress
REFERENCES
1. T. M. Prus-Chacinski, "Patterns of Motion
in Open-Channel Bends," Union of
Geodesy and Geophysics, Vol. 10, 9A,
1954, pp 311-318; also Proceedings, 6.
Institute of Water Engineering, London,
Vol. 10, Aug. 1956, pp 420-426.
2. T. M. Prus-Chacinski, and J. W. 7.
Wielogorski, "Secondary Motion Applied
to Storm Sewage Overflows," Symposium
on Storm Sewage Overflows, Institute of
Civil Engineers, London, 1967. 8.
3. J. G. Lloyd, "Report on the Model Tests
of the Nantwich Storm Sewage
Spiral-Flow Separator 5/14B,"
Unpublished report by Mersey and
Weaver River Authority, Warrington,
CIRIA, London.
4. J. A. Fox, and D. J. Ball, "The Analysis 9.
of Secondary Flow in Bends in Open
Channels," Proceedings, Institute of Civil
Engineering, London, Vol. 39, March
1968, pp 467-477.
5. V. T. Chow, Open Channel Hydraulics,
McGraw-Hill Book Company, New York
(1959), p 200.
W. E. Dobbins, "Effect of Turbulence on
Sedimentation," ASCE Transcript, Vol.
109, 1944, pp 629-656.
T. R. Camp, "Discussion of 'Effect of
Turbulence on Sedimentation' by W. E.
Dobbins," ASCE Transcript, Vol. 109,
1944, pp 660-667
R. R. Boericke, et al, "Mathematical
Model of the Swirl Concentrator as
Applied to Primary Separation of Sewage
and Combined Sewer Discharges," Draft
report submitted to American Public
Works Association, by General Electric
Co., June 15, 1974.
R. J. Dalrymple, et al, "Physical and
Settling Characteristics of Particulates in
Storm and Sanitary Wastewater," U.S.
EPA Report No. EPA-670/2-75-011.
115
-------
APPENDIX C
BACKGROUND INFORMATION ON THE HELICAL BEND
Report of C. H. Dobbie and Partners
by Dr. Prus-Chacinski
Experiments on the bends of open
channels during the last 20 years have
narrowed very considerably "the band of
uncertainty." Prior to the time the series of
experiments carried out at the Imperial
College in London between 1950 and
19561'16-18 and the Russian experiments
carried out about the same time6'7'13 were
completed, almost no information was
available on this hydraulic phenomenon.
Interest in the research on the bends of open
channels has been revived in recent years and
the experiments carried out in Karlsruhe are
of particular interest.2 •3 •'4
"The band of uncertainty" is however
still considerable. The Russian studies and the
most recent German studies were oriented
towards the erosion in the bends of natural
rivers and, therefore, the ratio of depth to
width not much less than unity was
experimentally investigated only in the
Imperial College studies, and later by the
studies oriented towards storm-sewage
overflows.1 9>2° Of interest have been recent
experiments on stratified flow5 where the
ratio of depth-to-width was unity. What to
date has not been achieved is a complete
investigaion of the series of bends with
different ratios of width-to-radius and with
varying ratio of depth-to-width. The Russian
studies are perhaps more complete in this
respect than any others, but for deep,.small
channels their results are not well
documented.
Considering the form of cross section and
the nature of previous hydraulic bend studies,
it would seem that virtually all forms of
simple cross sections were investigated with
the exception of the cross section with the
deepest part at the inside bend — the basic
form investigated by this (APWA) project.
Therefore, the present study, apart from its
value with regard to the design of combined
sewer overflows, will also fill at least one gap
in existing information on this subject.
Experimental or mathematical formulae
which would describe the pattern of helical
motion and its strength is open to
speculation. In a small, smooth rectangular
channel the expression which describes
accurately the value of the bed angle for the
fully-developed helical flow is one expressed
by Wadekar1 as:
tana0 =
(1)
where "C" depends on the value of the
secondary depth, that is the depth where the
transverse velocity is nil; i.e., where a. = 0. The
writer's equation:
tana0 =
R
(2)
is a "blanket" expression which seems to
provide a reasonable approximation for
smooth channels in a wide range of
depth-to-width ratio and of different cross
sections. It is very likely that the same
formula written in the form:
tan a0 =
wf/r
(3)
where "f' is the friction factor, may be the
best approximation which would also cover
rough channels. This has not yet been
established by any experiments. The Russian
approximation6:
tana0 =
(4)
may be valid for the natural channels; again, it
is a "blanket" expression covering many cross
sections but mostly with small depth-to-width
ratio. The Russian formula certainly
overestimates the value of "a0" in smooth
laboratory channels. The writer has observed,
however, during re-examination of his own
experimental data1 8 that in certain cases of
small rectangular channel bends, where the
bed was rough and "Re" reasonably large, the
116
-------
Russian formula is a good approximation. It
should be added that the Russian formula (4)
can be easily derived from either formula (2)
or (3): W = B + 2d. In shallow, wide channels,
the influence of "Re" or "f' will be small,
and as "B" (width) will tend to be constant, it
is obvious that formula (4) is a "blanket"
variation of formula (2) or (3).
Because the subject of helical flow is still
not very well known, care must be practiced
in accepting the published results of
experiments and, particularly, of their
interpretation. For example, it is difficult to
accept without further study the helix pattern
in the stratified flow at very low "Re," as
reported in reference 5, Mode A. This applies,
as well, to the statement in reference 2 and
reference 10 that the strength of helical
motion depends only on the ratio of
depth-to-width. Obviously it must also
depend on the ratio of width-to-central radius
and on the resistance to flow. In certain cases
control of the entry to the bend may be
suspected, resulting in unacceptable
experimental results. This is difficult to
establish in some published references. A
mathematical analysis has been
attempted1'6'7'9'12'13 but none of the
results so far achieved seem to be of much
help to the practical design. Recent studies on
turbulence21 by Yalin would seem to
demonstrate the origin of the existence of
helical motion in straight channels. The
existence of helical multicellular flow in
straight channels has been reported in several
references, not quoted here, for a long time.
Therefore, there has been some progress in
the mathematical analysis. On the other hand,
the general information on the behaviour of
liquids, and particularly of water, at the
molecular level, obviously most important in
the boundary layers phenomena, is still, to
say the least, incomplete. The recent opinions
of physicists with regard to the present
progress in this direction are rather
pessimistic.2 2
The Main Overflow Chamber
a) The shape of the main chamber, with
the deepest part at the inside wall is a novelty,
and to the writer's knowledge this shape has
never been investigated. Some Russian studies
included trapezoidal and triangular channels
but, because the research was oriented
towards the meanders of natural rivers, the
deepest part was always at the outside bend.
The triangular channel investigated by the
writer and described in his PhD thesis18 was
also deepest at the outside wall.
b) The reasons for adopting the novel
shape are as follows: It has been found by
Wadekar1 and confirmed by recent studies
carried out under Professor E. Mosonyi at the
University of Karlsruhe in Germany2'3 that in
a single bend the line of the maximum
velocity which at the entry to a curved
channel is at the inside wall, remains there
only up to about 40 degrees to 60 degrees of
the bend. Further, it crosses the channel and
remains at the outside wall up to 180 degrees
of the bend.
To the writer's knowledge, a bend longer
than 180 degrees has not yet been
investigated. Moreover, the region between 40
degrees and 60 degrees of the bend is where
the bottom slow water is already at the top.
In the writer's opinion, this region, in fact,
occurs a little further downstream, perhaps at
about 80 degrees of the bend. At the same
region, a second smaller helix forms at the
outer wall and grows bigger along the bend
downstream. This is now firmly established,
although there is still controversy with regard
to the ratio of width-to-depth which may
have an influence on the formation of the
second helix.2 In any case, at a region around
60 degrees of the bend the motion changes
from that of a free vortex to a composite
vortex. In the proposed novel,shape of cross
section, the line of maximum velocity will be
kept at the inside wall for a greater length of
the bend and so the upturning of the slow
bottom water will occur further downstream.
Therefore, the general condition of helical
flow should be stable for a longer distance
along the bend.
c) The relative strength of helical flow
was approximated by the writer as
S= 10tana0 %
where S = qs/ q or the ratio of the transverse
unit discharge to the forward unit discharge
and «0 is the bottom angle of flow at the
117
-------
center line of the rectangular curved channel.
The unit discharge is the discharge per unit
width at the center line of the channel.
The above-mentioned recent German
work carried out under Professor Mosonyi
(Ref. 3) has confirmed the writer's
approximation, except that S was integrated
across the channel and therefore Mosonyi's S
is only about 50 percent of the writer's S. Of
course,
-------
The Inlet
a) The inlet as now adopted is an open
channel transition 3.04 m (10 ft) long based
on an idea of a long diffuser of five degrees.
In the suggested design the open channel
diffuser is followed by 1.52 m (5 ft) length of
open channel of the same cross section as the
main overflow chamber.
b) To the writer's knowledge the
problem of stratification in storm sewage has
never been properly investigated. The solids in
storm water are seldom greater than 2,000
mg/1; the normal content is around 300 to
600 mg/1, or even less. BOD and ammonia are
linked to the solids in a very complex manner.
If the velocity in the incoming sewer is
around 1.82 m/sec (6 ft/sec) it is probably
not possible to say with any certainty
whether the flow is stratified or not; this will
depend on the content of solids and also on
their size. In addition, there is a problem of
flocculation; in fact, as far as the writer
knows, the subject is still obscure. However,
with a high content of solids it would be
surprising if there were no stratification.
c) The problem of waves of discharge is
further complicated by the well-known fact
that storm discharge is seldom uniform. It
occurs in the waves of high discharge lasting
five to ten minutes, which may be followed
by reasonably steady conditions of flow.
Often one incident of overflow demonstrates
two or even three waves. This means that the
hydrograph of storm runoff may be quite
complex. In long sewers the waves of the
runoff are attenuated, but never completely.
The result of these conditions is that the
content of solids is extremely variable because
each wave may contain more solids than in
steady flow.
d) The intended purpose of the inlet is to
decelerate the flow in a uniform manner and
to direct the region of the maximum velocity
to the inner wall of the main overflow
chamber. The straight section downstream
should result in reasonably steady conditions
with regard to any stratification.
e) The helical flow in open channel
bends with stratified liquids adds a further
complexity which has been described in
Reference 5. It has been found that the
helical flow in the bends of open channels
with two stratified liquids occurs in at least
three modes:
1. If stratification is complete and
velocity is low the heavy bottom
liquid is almost free of helical motion
and behaves as a false bottom. There
is then no helical motion in the
denser liquid but the helix in the
lighter liquid was observed to be of
the opposite nature than in the
normal flow, that is, the top
direction was towards the inner wall.
Such conditions may not be bad for
an overflow but probably they
seldom occur, if ever, in sewers. The
model observations carried out by
the Mersey River Authority
Hydraulic Laboratory have
demonstrated good results with salty
water at the bottom of the overflow
chamber but the Reynolds number in
the above experiments was always
greater than Re = 3500, as quoted in
Reference 5 for mode 1.
2. If the stratification is prominent and
the velocity is greater, the heavier
bottom liquid forms its own helix
which mixes only to a small degree
with the top helix of the lighter
liquid. The heavy liquid tends to
remain at the inner wall for at least
90 degrees of the bend. It is probable
that such conditions occur in the
storm water curved overflows if the
solids content is sufficiently large.
The lighter liquid may demonstrate
the presence of two helices after
about 60 degrees of the bend. It is
hoped that the model may provide
some information on this subject.
3. If the stratification is less prominent
or if the velocity is still higher, both
liquids will mix to a large degree and
only one helix occurs as in an
ordinary non-stratified flow. But
even then the heavier liquid is at the
inner wall for at least 45 degrees of
the bend. This case is probably
common in the curved stormwater
overflows if the content of the solids
is average.
119
-------
f) It is now an open question whether a
curved overflow should be designed for
conditions (2) or (3). The design for
conditions (3) would require a short side weir
starting at about five degrees of the bend and
finishing at about 55 degrees of the bend.
Design for conditions (2) would result in a
longer weir, up to 70 degrees or 80 degrees of
the bend. Again, experiments with the model
should result in a reasonable compromise
recommendation for the design. The model
investigated in Reference 5 was small; the
highest Reynolds number was only 9200.
Therefore, the question of stratification is still
far from being clarified.
g) The behavior of a model with steady
flow will not reproduce the effect of the
discharge waves. Such waves will invariably
disturb the helical motion and they will
decrease the efficiency of the overflow.
However, this problem is common to all the
types of stormwater overflow.
h) In view of the above discussion it is
possible that the proposed novel cross section
of the main overflow chamber may produce
better results than any other cross section;
this was an additional reason for the proposed
novel cross section. But it is evident that the
symmetrical cross section and the cross
section deepest at the inner wall should also
be investigated.
The Outlet
In British practice, a stormwater overflow
passes to the treatment works 6 to 7 DWF, or
sometimes even 8 DWF. It is, then, a
relatively simple matter to design an outflow
chamber in which the floating debris goes
over the end weir and the outflow is
controlled by the outflow sewer of such
dimensions as are dictated by the downstream
conditions. However, even with 7 DWF going
to treatment, variations in the outflow
discharge are observed, which are probably
caused by floating material accumulating in
the outlet chamber and temporarily
decreasing the capacity of the outlet sewer.
This, in turn, results in temporarily larger
overflows. This effect is not very large but it
does exist. The requirement of the APWA
Research Foundation is to pass 3 DWF to the
treatment works, and in such conditions the
effect may be more serious. The design of the
outflow chamber depends much on the
dimensions of the incoming sewer; for large
sewers it will be easier than for small sewers.
Alternative Design
The alternative design represented in
Figure 92, Alternative Design for Future
Investigation, may be an improvement as
compared with a simple overflow.-It would
seem to be possible, practical, effective, and
not particularly expensive.
The general concept is that, in addition to
the main chamber, the side weir and the end
outflow similar to those described already,
three side pipes should be installed at the
inner wall of the main chamber. These pipes
are connected into a side sewer of small
diameter which discharges into an end
manhole and which is closed by a flapgate.
The smallest flapgates obtainable in the
United Kingdom are 7.62 cm (3 in.) in
diameter.
A rough examination of discharges
demonstrates that such a design will work
safely if the inlet pipe is not smaller than
about 60.5 cm (2.0 ft) diameter. Then the
draw-off small sewer would be about 10.18
cm (4 in.) diameter and the off-take pipes
could be either 7.6 or 10.8 cm (3 or 4 in.)
diameter. The flapgate should then be
balanced in such a way that it would open
only if the inflow was almost 3 DWF, say 2.5
DWF. This would prevent the side pipes from
silting.
The actual design may be different than
that shown. For example, it would help if the
invert of the end or foul sewer passing to the
treatment works were approximately 15,24
cm (6 in.) lower than the invert of the main
overflow chamber. Also, the small diameter
side sewer could be straight, discharging
directly without a flapgate.
The main factor in the above alternative
is to avoid the silting up of the side pipes and
of the small diameter collecting sewer. The
flapgate should be balanced in such a way
that the velocity in the side pipes is. not
smaller than 60.5 cm/sec (2 ft/sec). The
alternative suggested in Figure 00 should be
easier to design for large than for small inlet
sewers. For example, if the inlet sewer were
120
-------
To treatment works .
FIGURE 92 ALTERNATIVE DESIGN FOR FUTURE INVESTIGATION
121
-------
90.5 cm (3.0 ft) diameter the draw-off sewer
would be about 20.32 cm (8 in.) diameter and
the side pipes could be 12.7 cm (5 in.)
diameter. It is possible that with the above
alternative, a simple symmetrical cross section
similar to that in the Nantwich overflow
might produce satisfactory results, but the
sharp edges should be eliminated. It is also
possible that the velocity in the main chamber
could be increased.
It is recommended that the above
possibilities should be investigated. It might
be sufficient to have only two side pipes. The
position of the side pipes and the position and
the length of the side weir should also be
investigated. The side pipes would discharge
directly to the side spill channel but some
means for collecting samples should be
provided.
Review of Laboratory Techniques
Known to the Writer
The list below is given in approximate
order of time when the reference was
published.
1. Shukry1 s has used the Pilot-sphere. Some
of his results are open to doubt because it
is evident that there was no proper
control of the entry to the channel bend.
2. Malouf16 used only the black dye —
Process Black — and the local velocity
was not measured.
3. Prus-Chacinski17'18 used Process Black
paint, small lighted floats, and the nylon
threads for direction measurements, using
an accurate protractor. Velocity was
measured by a pitometer.
4. Wadekar1 improved greatly on the
methods of (2) and (3). In addition to
Process Black paint he used droplets of
Meridian-unity oil (S G 1.002) as a tracer
for small velocities, and polystytene sand
(1 mm dia. SGI .06) as a tracer for larger
velocities. For the direction inside of the
channel, nylon threads and an accurate
protractor were used.
The surface direction and velocities were
measured by a sophisticated method of
photographing the surface wave patterns.
He was able to observe accurately the
direction of flow on the walls of his
rectangular channel because it was all
constructed of Perspex.®
5. The Russean researches known to
date6-7-13 have used methods similar to
those described under (2), (3) and (4). In
addition, the density flow effects were
investigated,13 demonstrating patterns, in
general, similar to those described in
Reference 5 case (3).
6. Fox and Bell8 used nylon thread as
direction indicators for velocity
measurement. They used a heat-transfer
recorder (not described) for velocities
smaller than 9.15 cm/sec 0.3 ft/sec), and
a multi-propeler or current meter with an
ultraviolet light recorder (not described)
for larger velocities.
7. Prus-Chacinski and Wielogorski19 used
Process Black paint and an accurate
protractor, plastic sand, orange colored (S
G 1.00), coal dust (particle size 0.1 to 1.0
mm), and polystyrene beads (S G 1.05)
for the bed patterns, and nylon thread for
direction inside of the channel. The
plastic sand also demonstrated the general
pattern of the helix in the body of the
flow. Perspex® shavings were used to
observe the surface flow direction. Local
velocities were not measured.
8. The Mersey River Authority, with the
cooperation of C. H. Dobbie and
Partners2 °, used the Process Black paint
and potassium permanganate. Potassium
permanganate was not satisfactory. The
threads were used, as already mentioned.
Local velocities were measured by a
mini-current meter. For surface flow
tracer, plastic beads were used. In addition
the density flow was investigated by
injecting the salty water colored by po-
tassium permanganate upstream of the
inflow. The results were most satisfactory
but only if the injection point was very
near to the bed. The analysis of overflow
was carried out, using the conductivity
meter.
9. Macagno and Alonso5 measured velocities
and densities with a constant temperature
hot-wire anemometer and electrical
conductometer, respectively. The signals
were processed through an IBM DACS.
122
-------
10. Francis and Safari9-10 used tracers and
sophisticated photographic and lighting
equipment. For laminar flow they used
TF.LCO powder (S G 0.9) with grain size
0.1 to 0.2 mm; turbulent flow droplets
(quite large) were formed by the mixture
of nitrobenzene, olive oil,,and water (S G
1.0). The rectangular channel had glass
windows in the walls. Local velocities
were not measured. It is possible that the
entry to their experimental channel was
not free from helical motions.
11. Mosonyi and- Gotz and Siebart and
Musers2'! 4used the hot-film anemometer
and a fibre-wire probe for -the local
velocity and local direction
measurements. This procedure seemed to
produce excellent results.
SYMBOLS USED
rc central radius of the bend in a
rectangular channel
Re Reynolds Number
C constant dependent upon depth
K constant to describe strength of helical
flow
W wetted perimeter
f friction factor
<*0 bed angle of
s strength of helical flow
qs/q forward unit discharge
B width of rectangular channel
d depth of rectangular channel
REFERENCES
Only a short list of references which introduce some new observations will be quoted below.
Otherwise, a long list of references is contained in the writer's PhD Thesis and in his paper
(together with Wielogorski):
Secondary motions applied to Storm Sewage Overflows
Symposium at the ICE, London, 1967
1. G. T. Wadekar, Secondary Flow in
Curved Channels, Unpublished PhD
Thesis, Imperial College, London, 1956.
2. Emil Mosonyi and Werner Gotz,
Secondary Currents in Subsequent Model
Bends, International Association for
Hydraulic Research, International
Symposium on River Mechanics,
Bangkok. January 1973.
3. E. Mosonyi and G. Meder, Effects of
Scale Distortion on Flow Properties in
Model Meanders, The Rehbock
Laboratory for River Improvement,
University of Karlsruhe, Germany, 1973.
4. P. Ackers, A Theoretical Consideration of
Side Weirs as Storm-Water Overflows,
Proceedings, Institute of Civil
Engineering, London, Febraury 1957.
5. E. O. Macagno and C. V. Alonso,
Two-layer Density-Stratified Flow in an
Open Channel Bend, International
Association for Hydraulic Research,
Proceedings of 14th Congress IAHR,
Paris, 1971
6. I. L. Rozovskii, Flow of Water in Bends
of Open Channel, (Translated from
Russian), The National Science
Foundation, Washington, D. C., 1967
(The Russian original, Kiev, 1957).
7. A. K. Ananyan, Fluid Flow in Bends of
Conduits, (Translated from Russian),
Israel Programme for Scientific
Translations, Jerusalem, 1965. (Russian
original, Erevan, 1957).
8. J. A. Fox and D. J. Bell, The Analysis of
Secondary Flow in Bends in Open
Channels, Proceedings, ICE, London,
March 1968.
9. J. R. D. Francis and A. F. Asfari,
Velocity Distributions in Wide, Curved
Open Channel Flows, Journal of Hyrauh'c
Research No. 1, 1971.
10. J. R. D. Francis and A. F. Asfari,
Visualisation of Spiral Motion in Curved
Open Channels of Large Width, Nature,
London, February, 1970.
11. Yoshio Muramoto, Secondary Flows in
Curved Open Channels, International
Association for Hydraulic Research,
Proceedings, 12th Congress IAHR, Paris,
1967.
123
-------
-------
tfc"..
------- |