See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/228581326

The Turn of the Screw: Optimal Design of an Archimedes Screw

Article  in  Journal of Hydraulic Engineering · January 2000

DOI: 10.1061/(ASCE)0733-9429(2000)126:1(72)

CITATIONS READS

130 4,468

1 author:

Chris Rorres
Drexel University
65 PUBLICATIONS   1,605 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Archimedes in the 21st Century View project

All content following this page was uploaded by Chris Rorres on 24 December 2018.

The user has requested enhancement of the downloaded file.

 THE TURN OF THE SCREW:
OPTIMAL DESIGN OF AN ARCHIMEDES SCREW
By Chris Rorres1

ABSTRACT: The geometry of an Archimedes screw is governed by certain external parameters (its outer radius,
length, and slope) and certain internal parameters (its inner radius, number of blades, and the pitch of the blades).
The external parameters are usually determined by the location of the screw and how much water is to be lifted.
The internal parameters, however, are free to be chosen to optimize the performance of the screw. In this paper
the inner radius and pitch that maximize the volume of water lifted in one turn of the screw are found. The
optimal parameter values found are compared with the values used in a screw described by the Roman architect
and engineer Vitruvius in the first century B.C., and with values used in the design of modern Archimedes
screws.

INTRODUCTION (or ‘‘pitch’’) of these blades was equal to the circumference

of the cylindrical core. Finally, an outer cylindrical covering
One of the oldest machines still in use is the Archimedes
(the ‘‘outer cylinder’’) of wooden planks was nailed to the
screw, a device for lifting water for irrigation and drainage
helical blades. Liquid pitch was smeared over all of the parts
purposes. Its invention has traditionally been credited to Ar-
during their assembly to make the screw watertight. The rigid
chimedes (circa 287–212 B.C.). For example, Diodorus Si-
screw was mounted so that it could be rotated along its length
culus (Greek historian, circa first century B.C.) writes
(Fig. 2) and was tilted in the direction of the hypotenuse of a
3-4-5 triangle (a ‘‘Pythagorean right-angled triangle,’’ as Vi-
men easily irrigate the whole of it [an island in the delta of truvius calls it). Its bottom end was immersed in a body of
the Nile] by means of a certain instrument conceived by
water (the ‘‘lower reservoir’’), and through the rotation of the
Archimedes of Syracuse, and which gets its name [cochlias]
because it has the form of a spiral or screw.
screw water was lifted to the top end of the screw into an
‘‘upper reservoir.’’
When an Archimedes screw is tilted, ‘‘buckets’’ that can
And from Athenaeus of Naucratis (Greek historian, circa A.D. trap water are formed between the blades. These buckets ap-
200): pear to move upward when the screw is rotated, carrying the
water within them. The screw collects water from the lower
The bilge-water [of the ship Syracusia], even when it be- reservoir, where the buckets are formed, and empties it into
came very deep, could easily be pumped out by one man
the upper reservoir, where the buckets are unformed. When
with the aid of the screw, an invention of Archimedes.
operated manually it is rotated by a crank (Fig. 3) or by a man
walking around the circumference of the outer cylinder in a
(See references under Diodorus Siculus and Athenaeus of Nau- treadmill manner.
cratis.) Archimedes, however, made no reference to it in his In modern industrial screws, the outer cylinder is usually
extant works [cf. Heath (1897) and Dijksterhuis (1938)], and fixed and the blades attached to the inner cylinder are rotated
it may be that he simply transmitted its knowledge from Egypt within it (Fig. 4). This allows the top half of the outer cylinder
(where it is believed he studied in Alexandria under the stu-
to be eliminated so that a stationary trough is formed from the
dents of Euclid) to Syracuse (his native city in Sicily). On the
bottom half of the outer cylinder. Such a construction permits
other hand, in defense of Archimedes: no mention of the de-
easy access to the interior of the screw, in order to remove
vice exists before his time, its design involves the type of
debris and for routine maintenance. In addition, the stationary
geometry in which he excelled, and his abilities as an inventor
outer cylinder relieves the moving blades and inner cylinder
of mechanical and military machines are well documented.
of some of the weight of the water. A disadvantage of this
The Roman engineer and architect Vitruvius gave a detailed
design is that water can leak down through the small gap be-
and informative description of the construction of an Archi-
tween the moving blades and the stationary trough. However,
medes screw in his De Architectura, written in the first century
this leakage can be considered an advantage in that it allows
B.C. (See reference under Vitruvius.) Vitruvius’s description
the screw to drain when it stops rotating.
contributed greatly to keeping the device well known through-
out the ages, and the particular screw he described will be The Archimedes screw has had a resurgence in recent years
used throughout this paper as a test case. because of its proven trouble-free design and its ability to lift
Vitruvius’s screw began with a tree trunk shaped into a cy- wastewater and debris-laden water effectively. It has also
lindrical core (the ‘‘inner cylinder’’) whose length is 16 times
its diameter (Fig. 1). On this cylindrical core eight intertwined
helical blades (also called ‘‘flights’’ or ‘‘starts’’) were con-
structed by nailing withes (slender flexible willow twigs) to-
gether up to a height equal to the radius of the core. The period
1
Prof., Dept. of Math. and Comp. Sci., Drexel Univ., Philadelphia, PA
19104. E-mail: crorres@mcs.drexel.edu
Note. Discussion open until June 1, 2000. To extend the closing date
one month, a written request must be filed with the ASCE Manager of
Journals. The manuscript for this paper was submitted for review and
possible publication on February 19, 1998. This paper is part of the Jour-
nal of Hydraulic Engineering, Vol. 126, No. 1, January, 2000. qASCE, FIG. 1. Vitruvius’s Eight-Bladed Screw as Described in his De
ISSN 0733-9429/00/0001-0072–0080/$8.00 1 $.50 per page. Paper No. Architectura. Diagram by Morris Hicky Morgan (Vitruvius Ref-
17708. erence)

72 / JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000

 FIG. 2. Profile of Three-Bladed Archimedes Screw

proved valuable in installations where damage to aquatic life

must be minimized.
The purpose of this paper is to examine the optimal design
of an Archimedes screw. By this is meant that geometry that
will maximize the amount of water delivered to the upper res-
ervoir in one turn of the screw. This optimization will be per-
formed under the assumption that the outer radius of the screw
is specified, for clearly the amount of water lifted in each
revolution can be continually increased by continually increas-
ing this radius.
The amount of water lifted per unit time can also be in-
creased by increasing the rotational velocity of the screw.
However, there is a practical limit to how fast one can rotate
the screw. A handbook on the design and operation of Archi-
medes screws (Nagel 1968, p. 37) states that, based on field
experience, the rotational velocity of a screw in revolutions
per minute should be no larger than 50/D2/3, where D is the
diameter of the outer cylinder in meters. Thus a screw with
FIG. 3. Two-Bladed Archimedes Screw in Use in Modern an outside diameter of 1 m should have a maximum rotational
Times in Egypt to Irrigate Fields velocity of 50 rpm. If the screw is rotated much faster, tur-

FIG. 4. Seven 96-Inch-Diameter, 3-Flight Screw Pumps That Lift Wastewater from Pretreatment Section to Balance of Treatment Plant
in Memphis, Tenn. [Each Pump Has Capacity of 19,900 gal./min (Courtesy Lakeside Equipment Corp., Bartlett, Ill.)]

JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000 / 73

bulence and sloshing prevent the buckets from being filled and volume of one bucket. The volume of one bucket is denoted
the screw simply churns the water in the lower reservoir rather by VB. It is a complicated function of N, K, L, Ro and Ri.
than lifting it. The restriction L # 2pRo /K given above on the pitch L of
the screw needs some discussion. In order for water to be
PROBLEM FORMULATION trapped in the screw, it is necessary that the sinusoidal curve
in Fig. 5 defining the outer edge of a blade tilts downward as
Fig. 5 is a profile view of a segment of an Archimedes screw it crosses the axis of the screw. In terms of the angles u and
showing how the buckets are formed between pairs of adjacent a in the figure, it is necessary that u # a or, equivalently,
blades. The buckets move up the screw as it is rotated clock- tan u # tan a. Now tan u = K by definition, and tan a =
wise when viewed from the upper reservoir. (Throughout this Ro(2p/L) since the sinusoidal curve has amplitude R0 and
paper the screws are oriented for clockwise rotation). The in- period L. The condition tan u # tan a is thus K # Ro(2p/L),
ner and outer edges of each blade determine two sinusoidal which is just L # 2pRo /K.
curves of the same period (or pitch) and phase. One has an Fig. 5 also shows the angle b that the sinusoidal curve de-
amplitude equal to the inner radius of the screw and the other fining the inner edge of a blade makes with axis of the screw.
has an amplitude equal to the outer radius. The angle u that Because the inner sinusoidal has amplitude Ri and period L,
the screw makes with the horizontal determines its slope tan u. it follows that tan b = Ri(2p/L). If u # b, so that the inner
Now defined are the following three ‘‘external’’ parameters: sinusoidal curve tilts downward as it crosses the axis of the
screw, then the horizontal water level in a bucket is tangent
Ro = radius of screw’s outer cylinder (m) to this inner sinusoidal curve. As above, the condition u # b
can be expressed as L # 2pRi /K.
L = total length of screw (m)
If L lies in the interval (2pRi /K, 2pRo /K ), then the hori-
K = slope of screw (dimensionless) zontal water level is right at the point where the inner and
outer sinusoidal curves cross the axis of the screw. At that
These external parameters are usually determined by the site point the inner sinusoidal curve has positive slope (with re-
of the screw and the materials available for its construction. spect to the horizontal) and the outer sinusoidal curve has neg-
In this paper, these three parameters are taken as fixed. In ative slope. The buckets of water that form in this case are
addition to these parameters, the following three ‘‘internal’’ rather small and in, in fact, will not be in contact with the
parameters are needed to completely specify the geometry of inner cylinder. Summarizing, the horizontal water level of a
the screw: bucket is
Ri = radius of screw’s inner cylinder (m) (0 # Ri # Ro) 1. Tangent to the inner sinusoidal curve in a profile view if
L [ (0, 2pRi /K ) [Fig. 6(a)].
L = pitch (or period) of one blade (m) (0 # L # 2pRo /K )
2. At the intersection of the outer and inner sinusoidal
N = number of blades (dimensionless), N = 1, 2, . . . , curves in a profile view if L [ (2pRi /K, 2pRo /K ) [Fig.
6(b)].
By ‘‘one cycle of the screw’’ is meant a segment of the screw
whose length is equal to one pitch of the screw. The volume Returning to our optimization problem, the following quantity
of one cycle of the screw is pR 2oL. is now defined:
By a ‘‘chute’’ is meant a region of the screw bounded by
VT = volume of water in one cycle of the screw (m3 )
two adjacent blades and the inner and outer cylinders. The
region between the inner and outer cylinders of the screw con- This quantity is also the volume of water emptied into the
sists of N disjoint congruent chutes separated by the N blades. upper reservoir with each turn of the screw—precisely the
The volume of each chute is p(R 2o 2 R 2i )L/N. quantity to be maximized. Note also that VT = NVB, expressing
By a ‘‘bucket’’ is meant one of the maximally connected the fact that N buckets of water are emptied into the upper
regions occupied by the trapped water within any one chute. reservoir with each turn of the screw. The basic problem of
Each bucket is filled with water at the lower reservoir and is this paper can now be stated as follows:
emptied into the upper reservoir. Within one cycle of the screw
the volume of the water within all the buckets is N times the Problem: Given N, Ro, and K, find values of Ri and L that
maximize VT.

In words: Given the number of blades, outer radius, and slope

FIG. 6. Views of Single Bucket Looking Down Two Screws

with Transparent Blades and Different Pitches (Both Screws Are
8-Bladed and Their Inner Radius is Half Their Outer Radius;
Dark-Shaded Region Is Horizontal Water Surface and Light-
FIG. 5. Profile View of Segment of Two-Bladed Archimedes Shaded Region Is Portion of Bucket in Contact with Blade
Screw Nearer the Top of Screw)

74 / JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000

of an Archimedes screw, find the inner radius and pitch that
will maximize the volume of water emptied into the upper
reservoir with each turn of the screw.
V*
T = S D
2p 2 R o3
K
l*n(N, r*, l*) (7)

VOLUME OF A BUCKET
In this paper it is assumed that the blades have negligible
thickness. With this assumption the volume of water in one The algorithm for computing the volume ratio n(N, r, l)
cycle monotonically increases with the number N of blades. If for a specific set of screw parameters N, r, and l is described
the blades had some nonnegligible thickness, however, they in Appendix I. As an example of this algorithm, consider the
would occupy an increasing fraction of the volume of the volume ratio for Vitruvius’s 8-bladed screw. Since r = 1/2 and
screw as their number increased and a point of diminishing l = 3/8 for his screw, the desired volume ratio is n(8, 1/2,
returns would be reached. One could then include the deter- 3/8). The algorithm gave 0.1703 for this volume ratio, so that
mination of an optimal N in the problem statement. But in 17% of Vitruvius’s screw was occupied by water when oper-
modern screws the number of blades is usually 1, 2, or 3, ating. With this value of n(8, 1/2, 3/8), Eqs. (1–3) and the
because of manufacturing, weight, and cost constraints. It is fact that K = 3/4 leads to
thus assumed that the number of blades is predetermined, and
their thickness is neglected in the computations. VT = 1.68R 3o (8)
In the next section this formula is compared with the corre-
DIMENSIONLESS PARAMETERS sponding formula that the optimization algorithm provides.
Fig. 7 is plot of n(N, r, l) for N = 8 and all values of r
The analysis of the optimization problem is begun by de-
and l between 0 and 1. Plots of n(N, r, l) for other values
fining a dimensionless parameter n as the ratio of VT (the vol-
of N are similar, with n(N, r, l) monotonically increasing as
ume of water in one cycle of the screw) to pR 2oL (the total
N increases for each fixed r and l.
volume of one cycle of the screw):
Fig. 7 shows that the volume ratio has a global maximum
VT for some value of r when l = 0. The case l = 0 arises when
n= = volume ratio (1) either the slope or the pitch of the screw is zero. Although
pR 2oL
neither situation is practical, they may be considered as the
= fraction of the volume of one cycle of the screw occupied limiting cases when the slope or the pitch approaches zero.
by water. Fig. 8 is a typical cross section of a screw when l = 0, show-
This volume ratio is also the fraction of the entire screw ing how the water level is tangent to the top of the inner
occupied by water if the ends where water is either entering cylinder. Every cross section of the screw has the same water
or leaving the screw are neglected. By its definition n is a profile when l = 0 and this is true regardless of the number
number between 0 and 1. Notice that maximizing VT is not the
same as maximizing n.
Two more dimensionless parameters are defined as follows:
Ri
r= = radius ratio (0 # r # 1) (2)
Ro
and
KL
l= = pitch ratio (0 # l # 1) (3)
2pRo
Because Ro = maximum value of the inner radius Ri, the radius
ratio r is the ratio of the actual inner radius to its maximum
possible value. Similarly, because 2pRo /K = maximum pitch
for which buckets will form, the pitch ratio l is the ratio of
the actual pitch to its maximum possible value. In the case of
Vitruvius’s screw: Ri = (1/2)Ro, L = 2pRi, K = 3/4, and so r
= 1/2 and l = 3/8.
A dimensional analysis shows that n depends only on N, r,
and l. This quantity can be written as n(N, r, l) to emphasize
this fact. From the above three equations it follows that

VT = S D
2p 2 R o3
K
ln(N, r, l) (4)
FIG. 7. Graph of n(N, r, l) (Volume Ratio) versus r (Radius Ra-
tio) and l (Pitch Ratio) for 8-Bladed Screw (N = 8)
Given N, Ro, and K, the problem of maximizing VT with re-
spect to Ri and L can then be reduced to maximizing ln(N,
r, l) with respect to r and l, each restricted to the interval
[0, 1]. Let the values of r and l that maximize ln(N, r, l)
be denoted by r* and l*, respectively. Then the optimal values
of Ri, L, and VT are given by

R*
i = r*Ro (5)
2pRol*
L* = (6)
K
FIG. 8. Cross Section of 8-Bladed Screw with Zero Pitch Ratio
and (l = 0)

JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000 / 75

of blades the screw has. A little geometry then shows that the OPTIMAL SCREW
fraction of the area of the cross section of the screw occupied
by water is given by The determination of the optimal screw is now considered.
This is the problem of finding the values of r and l in the
1 r interval [0, 1] that maximize the function of ln(N, r, l). Fig.
A0(r) = (1 2 r2) 2 cos21(r) 1 Ï1 2 r2 (9) 9 is a plot ln(N, r, l) for an 8-bladed screw. As with n(N, r,
p p
l), plots of ln(N, r, l) for different values of N are similar
Now, because each cross section of the screw is occupied by to each other, with ln(N, r, l) monotonically increasing as N
this same fraction of water, the volume ratio of the screw must increases.
also be given by this expression. Consequently, (9) is the for- It can be seen from Fig. 9 that ln(N, r, l) has a unique
mula for the curve obtained by slicing the surface in Fig. 7 maximum inside the unit square. This maximum was com-
with the plane l = 0; that is, A0(r) = n(N, r, 0). A little calculus puted numerically for various N using the ‘‘fmin’’ function of
shows that A0(r) is maximized when r = 1/Ï1 1 p2 = MatLab (Hanselman and Littlefield 1995), which implements
0.30331. . . . The corresponding value of n is 0.598093. . . . a Nelder-Mead simplex search (Nelder and Mead 1965). For
Thus no Archimedes screw can have more than 60% of its an 8-bladed screw, the peak value of the volume-per-turn ratio
volume occupied by water, and this maximum value is ap- was found to be 0.0771, attained when r* = 0.5354 and l* =
proached when the slope or pitch is very small and the inner 0.2957. Eq. (7) then shows that the volume per turn in the
radius of the screw is about 30% of the outer radius. optimal case is given by
1.52R 3o
V*
T = (10)
K
Table 1 gives the optimal values of the pitch and radius ratios
for screws with 1 to 25 blades, together with the corresponding
optimal values of the volume-per-turn and the volume ratios.
The last row gives the limiting value of these quantities as the
number of blades approaches infinity. It provides an upper
bound for the amount of water that can be lifted in one turn
since the volume-per-turn ratio monotonically increases as N
increases.
Fig. 10 represents the data in Table 1 in graphical form.
Notice how r* varies only in the interval (0.5352, 0.5394),
and so to two decimal places the optimal ratio is 0.54 for any
number of blades. It may be that the optimal radius is the same
for any number of blades and this slight variation of r* may

TABLE 1. Optimal Ratio Parameters of Archimedes Screw for

Various Numbers of Blades

Optimal
Number of Optimal Optimal volume-per- Optimal
blades radius ratio pitch ratio turn ratio volume ratio
N r* l* l*n(N, r*, l*) n(N, r*, l*)
(1) (2) (3) (4) (5)
1 0.5358 0.1285 0.0361 0.2811
2 0.5369 0.1863 0.0512 0.2747
3 0.5357 0.2217 0.0598 0.2697
4 0.5353 0.2456 0.0655 0.2667
5 0.5352 0.2630 0.0696 0.2647
6 0.5353 0.2763 0.0727 0.2631
7 0.5354 0.2869 0.0752 0.2619
8 0.5354 0.2957 0.0771 0.2609
9 0.5356 0.3029 0.0788 0.2601
10 0.5356 0.3092 0.0802 0.2592
11 0.5358 0.3145 0.0813 0.2586
12 0.5360 0.3193 0.0824 0.2580
13 0.5360 0.3234 0.0833 0.2574
14 0.5360 0.3270 0.0841 0.2571
15 0.5364 0.3303 0.0848 0.2567
16 0.5362 0.3333 0.0854 0.2562
17 0.5362 0.3364 0.0860 0.2556
18 0.5368 0.3380 0.0865 0.2559
19 0.5364 0.3404 0.0870 0.2555
20 0.5365 0.3426 0.0874 0.2551
21 0.5370 0.3440 0.0878 0.2553
22 0.5365 0.3465 0.0882 0.2544
23 0.5369 0.3481 0.0885 0.2543
24 0.5367 0.3500 0.0888 0.2538
25 0.5371 0.3507 0.0891 0.2542
? ? ? ? ?
? ? ? ? ?
FIG. 9. Two views of Graph of ln(N, r, l) (Volume-per-Turn Ra- ? ? ? ? ?
tio) versus r (Radius Ratio) and l (Pitch Ratio) for 8-Bladed ` 0.5394 0.3953 0.0977 0.2471
Screw (N = 8)

76 / JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000

be a numerical artifact resulting from the approximate numer-
ical integration needed to compute n(N, r, l) and the approx-
imate numerical optimization of ln(N, r, l).
Table 1 also shows that n(N, r*, l*) decreases from 0.2811
to 0.2471 as N increases from one to infinity. Thus, a decreas-
ing fraction of the screw is occupied by water as the number
of blades increases.
Fig. 11 shows a view of a single bucket looking down the
screw for the optimal screws with 1, 2, 3, or 4 blades. Notice
how the horizontal water surface is in one piece for a 1- or 2-
bladed piece but is in two pieces for a 3- or 4-bladed screw.
When the water surface is in one piece, there is a clear air
passage from the top to the bottom of the screw if the screw
has a complete watertight outer cylinder. However, when the
water surface is broken into two pieces by the inner cylinder,
the buckets of water close off the chutes and so air is trapped
in the spaces between the buckets in any one chute.

FIG. 11. Views Looking down Optimal 1-, 2-, 3-, and 4-Bladed
Screw of Single Bucket Lying Between Two Adjacent Blades
(Dark-Shaded Region Is Horizontal Water Surface and Light-
Shaded Region Is Portion of bucket in Contact with Blade nearer
Top of Screw)

CONCLUSIONS
Vitruvius’s screw configuration can be compared with its
optimal configuration as determined by the calculations above.
As shown in (8), the volume per turn of the Vitruvius screw
is governed by VT = 1.68R 3o. The optimal volume per turn is
3
governed by (10) with K = 3/4, which gives V * T = 2.03R o. The
optimal design thus results in a fractional increase in the water
lifted per turn given by V *T /VT = 2.03/1.68 = 1.21, or a per-
centage increase of 21%. Consequently, the output of Vitru-
vius’s screw is fairly close to that of the optimal 8-bladed
screw. In addition, its radius ratio is within 7% of the optimal
value (0.5 versus 0.5354), its pitch ratio is within 27% of the
optimal value (0.375 versus 0.2957), and the construction lines
associated with its design are much simpler than those that
would be needed to construct the optimal screw. No doubt
many generations of experience went into the design of the
screw that Vitruvius described.
Fig. 12 shows one bucket in a screw of Vitruvius’s design
and one bucket in a screw of the corresponding optimal design.
Ritz-Atro Pumpwerksbau of Nürnberg, Germany, a manu-
facturer of Archimedes screws, prepared an Archimedean
Screw Pump Handbook in 1968 to ‘‘provide the technical in-
formation needed for the calculation, planning, construction,
and operation of water pumping installations using Archime-
dean screw pumps’’ (Nagel 1968). This handbook gives cer-
tain rules of thumb for maximizing the volume of water raised
with each turn of the screw based on heuristic arguments and
field experience. In particular, it states that the ratio of the
inner radius to the outer radius should be between 0.45 and
0.55 for most conditions (Nagel 1968, p. 30). The results of
this paper show that it should be 0.54 under all circumstances
when the corresponding optimal pitch is used.
For the pitch, the handbook gives the following heuristic
rule for a screw tilted at an angle u (Nagel 1968, p. 31):
L = 2.4Ro if u < 307
L = 2.0Ro if u = 307
L = 1.6Ro if u > 307
The result obtained in this paper is given in (6): L = (2pl*
FIG. 10. Graphical Representation of Data in Table 1 cot u)Ro, where the numerical value of l* is given in Table 1
JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000 / 77
and depends on N. Fig. 13 is a graph of L as a function of u APPENDIX I. VOLUME OF A BUCKET (DETAILS)
as given by this formula for the cases N = 1, 2, and 3. Also
in Fig. 13 is a graph of the piecewise constant formula deter- This appendix describes the algorithm used in the paper for
mined by the handbook’s heuristic rule. The heuristic rule can the computation of the volume ratio n(N, r, l). Run an x-axis
be seen to be a pretty good match for the exact relationship, along the axis of the screw, with x increasing going up the
especially for a 2-bladed screw. The handbook further states screw. Select any particular bucket of the screw at any partic-
that to avoid a proliferation of screw designs and the concom- ular time and set
itant manufacturing costs, the choice L = 2.0Ro should be AB(x) = Area (m2) of the water in the cross section of the screw
made for all angles. Fig. 13 justifies this particular simplifying
choice for 1-, 2-, and 3-bladed screws and for angles around at position x (m) in the selected bucket.
307. The volume of VB of the bucket is the integral of AB(x) over
Finally, it seems appropriate to close this paper with the the extent of the bucket along the x-axis. Thus, if (x0, x1) is
following words that Vitruvius used to close the description the interval along the x-axis over which the bucket is defined,
of his screw: then

E
x1
I have now described as clearly as I could, to make them
VB = AB(x) dx (11)
better known, the principles on which wooden engines for x0
raising water are constructed, and how they get their motion
so that they may be of unlimited usefulness through their Next define a new (dimensionless) variable by
revolutions.
2px
f= (12)
L
called the blade angle. The blade angle goes through 2p rad
over one cycle of the screw. With respect to the blade angle,
(11) becomes

E
f1
L ˆ B(f) df
VB = A (13)
2p f0

where ÂB(f) = function AB(x) expressed as a function of f,

and the new limits of integration are f0 = 2px0 /L and f1 =
2px1/L. Next define the dimensionless parameter
ÂB(f)
gB(f) = = area ratio (14)
pR o2
FIG. 12. Two Views Looking down Two 8-Bladed Screws of
Single Bucket Lying Between Two Adjacent Transparent Blades
(Dark-Shaded Region Is Horizontal Water Surface and Light-
Shaded Region Is Portion of Bucket in Contact with Blade
nearer Top of Screw)

FIG. 13. Graph of Ratio of Pitch of Screw to Its Outer Radius

versus Angle of Inclination of Screw When Number of Blades Is
1, 2, or 3 [Also Shown Is Piecewise Constant Approximation FIG. 14. Area Ratio of Typical Bucket of 3-Bladed Screw as
Discussed in Archimedean Screw Pump Handbook (Nagel Function of Blade Angle
1968)]

78 / JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000

which is the ratio of the cross-sectional area of the water in repeats every 1207 or, more generally, every 2p/N radians for
one bucket to the cross-sectional area of the screw. Combining an N-bladed screw. Notice that for angles close to 3307 the
(13) and (14) with (1), then, results in cross section of the bucket completely fills the cross section

E
f1 of the chute so that the value of gB is 1/4, since for this screw
VT NVB N the cross-sectional area of each chute is one-fourth the cross-
n= = = gB(f) df (15)
pR 2o L pR 2o L 2p f0 sectional area of the entire screw.
Fig. 16 illustrates how certain angles associated with the
The above concepts are illustrated in Fig. 14 with a 3-bladed geometry of a bucket are determined. In both the profiles and
screw for which K = 3/8, L = pRo, and Ri = (1/2)Ro (thus, cross-sectional views only a single bucket is shown and only
N = 3, r = 1/2, and l = 3/16). The top diagram shows the the two blades of the chute containing the bucket are shown.
profile of the screw with the bucket shaded. Of course, the Additionally, the surfaces of the two blades and the inner and
screw is actually tilting upward with slope K = 3/8 and the outer cylinders have been made transparent. Looking down
straight lines slanting downward in the diagram are the hori- into the screw in the cross-sectional view, the darker region is
zontal water levels in the buckets. The middle diagram shows the top horizontal water surface of the bucket and the lighter
a single bucket with the angles f0 and f1 indicated. The origin region is the portion of the bucket in contact with the upper
has been chosen along the f-axis so that the equation of the blade of the chute. Letting y denote distance above the axis of
sine curve forming the outer edge of the lower blade of the the screw in the profile view, the following curves in this fig-
bucket is Ro sin f. The bottom diagram is a graph of gB(f) ure have the following equations:
versus f of the particular bucket selected. (The top curve in
the bottom graph is the cumulative area ratio of all of the Outer edge of the lower (left) blade: y = Ro sin f
buckets in the screw. It is a periodic curve with period 2p/N Inner edge of the lower (left) blade: y = Ri sin f
and its average value is n.)
Fig. 15 shows cross sections of this 3-bladed screw for val- Outer edge of the upper (right) blade: y = Ro sin(f 2 2p/N)
ues of the blade angle from 907 to 4207. The chute containing Inner edge of the upper (right) blade: y = Ri sin(f 2 2p/N)
the selected bucket is outlined heavily, and the cross section
of the bucket itself is shaded darker than the buckets in the KL
Water level in the bucket: y = 2 (f 2 f0) 1 Ri sin f0
other two chutes. The cycle of cross sections in Fig. 15 repeats 2p
with each pitch of the screw, that is, for every 3607 increase
of the blade angle. In fact, the pattern of cross sections actually The angle f0, where the bucket begins, is where the water
level is tangent to the inner edge of the lower blade. Thus
KL l
Ri cos f0 = 2 or cos f0 = 2 (16)
2p r
(This is assuming l < r; otherwise f0 = p.) The angle f1
where the bucket ends is where the water level hits the outer
edge of the upper blade for a third time. Thus

R0 sin Sf1 2
2p
N D =2
KL
2p
(f1 2 f0) 1 Ri sin f0

or S
sin f1 2
2p
N D = 2l(f1 2 f0) 1 r sin f0
(17)
The angles denoted fL or fR determine the portion of the outer
cylinder in contact with the water. If the outer cylinder is an
open trough, then they determine the minimum angular bound-
aries of the trough needed to contain the buckets of water.
Both angles are determined by points of intersection of the
water level and the outer edge of the lower blade, so that they
are both solutions of the equation
KL
Ro sin f = 2 (f 2 f0) 1 Ri sin f0
FIG. 15. Cross Section of Screw for Various Values of Blade 2p
Angle (Water in One Chute Is Darkened for Easier Identification)
or sin f = 2l(f 2 f0) 1 r sin f0 (18)

FIG. 16. Profile and Cross-Sectional Views of Screw in Figs. 14 and 15

JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000 / 79

For Vitruvius’s screw (r = 1/2, l = 3/16) these angles turn out by Ivor Thomas in Greek mathematical works, Loeb Classical Library,
to f0 = 112.027, f1 = 441.967, fL = 162.687, and fR = 343.007. Harvard University Press, Cambridge, Mass., 1941.
Hanselman, D., and Littlefield, B. (1995). The student edition of
Depending on the parameters of the screw, there are many MATLAB: Version 4: User’s guide. Prentice-Hall, Englewood Cliffs,
other possibilities as to how the water level can hit the four N.J., 419–421.
sinusoidal curves and the inner and outer cylinders in Fig. 16. Heath, T. L. (1897). The works of Archimedes, originally published by
One must be careful with the bookkeeping to make sure the Cambridge University Press. Republished by Dover, New York, 1953,
correct angles are chose in the correct ranges. pages xx–xxi.
As seen in Fig. 15, each cross section of a bucket consists Nagel, G. (1968). Archimedean screw pump handbook. Prepared for Ritz-
Atro Pumpwerksbau GMBH Roding, Nürnberg, Germany.
of a portion of the cross section of a chute cut off by a straight Nelder, J. A., and Mead, R. (1965). ‘‘A simplex method for function
line representing the horizontal water level of the bucket in minimization.’’ Computer J., 7, 308–313.
that cross section. The cross section of a chute, in turn, consists Vitruvius (circa first century B.C.). De Architectura, Book X, Chapter VI,
of a region bounded by two concentric circles and two rays. ‘‘The Water Screw.’’ Translation and diagram by Morris Hicky Morgan
Basic geometry can be used to find exact formulas for the in Vitruvius: The ten books on architecture, Harvard University Press,
function gB(f), although the formulas are quite messy because Cambridge, 1914. Republished by Dover, New York, 1960, 295–297.
there are many ways that the horizontal water line can hit the
two circular boundaries and two straight-line boundaries of a APPENDIX III. NOTATION
chute cross section. The following symbols are used in this paper:
Even though exact formulas for gB(f) for each f can be
generated, it is not possible to obtain an exact expression for A0 = fraction of area of cross section of horizontal screw occu-
the integral in (15) that determines the volume ratio. This is pied by water (dimensionless);
because the integral in (15) must be written as a sum of sub- AB = area of water in cross section of one bucket (m2);
integrals, each for a different shape of the bucket cross section, K = slope of screw (dimensionless);
and the angles needed as limits for these subintegrals (such as L = total length of screw (m);
the angles f1, fL, and fR) cannot be expressed in closed form. N = number of blades in screw (dimensionless);
Consequently, a MatLab computer program was written to de- Ri = radius of screw’s inner cylinder (m);
Ro = radius of screw’s outer cylinder (m);
termine the volume ratio numerically. This program imple-
VB = volume of one bucket of water (m3);
mented Simpson’s Rule to evaluate the integral in (15) using Vt = volume of water lifted in one turn of screw (m3);
a blade-angle spacing of 27. x = position along axis of screw (m);
y = distance perpendicular to axis of screw (m);
ACKNOWLEDGMENTS a = angle of incline of spiral intersection of blade and outer
The writer would like to thank John Urban of Jupiter, Florida, for cylinder with respect to axis of screw (rad);
bringing the Archimedean Screw Pump Handbook (Nagel 1968) to his b = angle of incline of spiral intersection of blade and inner
attention and for securing him a copy. cylinder with respect to axis of screw (rad);
u = angle of incline of screw (rad);
APPENDIX II. REFERENCES L = pitch (or period) of blade (m);
l = pitch ratio = KL/2pRo (dimensionless);
Athenaeus of Naucratis. (circa A.D. 200). Deipnosophistae, Book V. n = fraction of volume of one cycle of screw occupied by water
Translation by Charles Burton Gulick in Athenaeus, Deipnosophists, (volume ratio) (dimensionless);
Loeb Classical Library, Harvard University Press, Cambridge, Mass., r = radius ratio = Ri /Ro (dimensionless); and
1928.
Dijksterhuis, E. J. (1938). Archimedes, original Dutch edition published
f = blade angle = 2px/L (rad).
by I. P. Noordhoff, Groningen. English translation by C. Dikshoorn,
Princeton University Press, Princeton, N.J., 1987, 21–23. (Note: The optimal value of a parameter is denoted by an
Diodorus Siculus (circa first century B.C.). Bibliotheke, i.34.2. Translation asterisk; e.g., V *
T .)

80 / JOURNAL OF HYDRAULIC ENGINEERING / JANUARY 2000

View publication stats