Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

(Swanson 1961) Swanson Magnus Effect Review Paper

Download as pdf or txt
Download as pdf or txt
You are on page 1of 10
At a glance
Powered by AI
The paper summarizes experimental results on the Magnus effect and mathematical attempts to model it. The Magnus effect causes spinning objects to experience a lateral force when traveling through a fluid and is responsible for deviations seen in balls, missiles, and other projectiles.

The Magnus effect is a lateral aerodynamic force experienced by a spinning object moving through a fluid. It is caused by the interaction between the object's rotation and its velocity through the fluid, creating an asymmetrical pressure distribution via the Bernoulli effect.

For cases where the axis of rotation is not aligned with the flight path, or where there is a large yaw angle between them, the flow separates at moderate to high Reynolds numbers, making mathematical modeling very difficult. The problem becomes even more complex for three-dimensional objects with different geometries.

W. M.

SWANSON
Associote Professor o f M e c h a n i c a l Engineering,
W a s h i n g t o n U n i v e r s i t y , St. Louis, M o .
Assoc. M e m . A S M E

The Magnus Effect: A Summary of


Investigations to Date
The Magnus force on a rotating body traveling through a fluid is partly responsible for
ballistic missile and rifle shell inaccuracies and dispersion and for the strange deviational behavior of such spherical missiles as golfballs and baseballs. A great deal of
effort has been expended in attempts to predict the lift and drag forces as functions of the
primary parameters, Reynolds number, ratio of peripheral to free-stream velocity, and
geometry. The formulation and solution of the mathematical problem is of sufficient
difficulty that experimental results give the only reliable information on the phenomenon.
This paper summarizes some of the experimental results to date and the mathematical
attacks that have been made on the problem.

Introduction

cover so many different geometr'es that they are only briefly con
sidered here.

A SPINNING missile or body traveling through the air


in such a way that the body axis of rotation is at an angle with the
flight path will experience a Magnus force component in a direction perpendicular to the plane in which the flight path and rotational axis lie. The magnitude of this force is a function of the
spin rate, the flight velocity, and the shape of the missile.
This Magnus force has been the subject of much investigation
in order to determine its effect on the flight path of spinning shells
and as a possible high lift device. Many ballistic missiles are given
some spin for stabilization even though some may have fins for
guidance in the propelled and guided stages. For such missiles
the angle of yaw between the flight path and spin axis is usually
small so that a small perturbation analysis is possible. For large
yaw, or when the axis is at right angles to the flight direction, the
flow separates (in the case of moderate or large Reynolds numbers) and any method of analysis, even approximate, is virtually
impossible at the present time.
A review of some of the solutions is briefly presented along with
the experimental results for the particular case of a two-dimensional (infinitely long) cylinder. The three-dimensional results
Contributed by the Fluid Mechanics Subcommittee of the Hydraulic Division and presented at the Winter Annual Meeting,
N e w York, N . Y . , November 27-December 2, 1960, of THE AMERICAN
SOCIETY OF MECHANICAL ENGINEERS.

Headquarters, August 1, 1960.

M a n u s c r i p t received at A S M E

Paper No. 60WA-150.

Historical Notes
The first record of the drift deviation of a spinning body was
described b y G. T . Walker in 1671. The body was a "sliced"
tennis ball. G. Magnus [l], 1 in faying to account for the drift of
spinning projectiles, performed some crude experiments with
musket balls and with a cylinder in an air jet. He correctly
ascribed the drift to an aerodynamic force produced by the
interaction between the rotation and flight velocity which gave
rise to an unsymmetrical pressure distribution produced by the
Bernoulli effect. His projectile drift experiments were performed
b y laying musket balls with their center of mass not coincident
with the geometric center in a musket with their center of mass
in a known orientation. If the center of mass were to the
right, an impulsive clockwise rotation about a vertical axis would
be produced b y the center of pressure acting through the geometric center and the inertial reaction acting through the center
of mass. When fired, the ball should drift to the rightit did.
Lord Rayleigh was the first to set up the ideal flow representation. In 1877 he published a paper "On the Irregular Flight of a
Tennis Ball" [2] in which he presented the mathematical model
of the flow as the classical potential flow around a cylinder with
circulation. A t the time he also stated that it was not possible
to give a complete mathematical formulation of the actual physi1

Numbers in brackets designate References at end of paper.

.Nomenclature,
a

cylinder radius for case of unbounded flow around a


cylinder
c = radial distance to external vortex, z = eeiy
CD = drag coefficient
CDo
= static drag coefficient
CL = lift coefficient
Cp

L = lift
p =

r = radius, radial co-ordinate


R = Reynolds number
u = radial velocity component, u(r, 6), general notation
U

- pressure coefficient =

drag

= free stream, uniform velocity of approach at x =

v = tangential velocity component, v(r, 6); general notation

i = imaginary quantity
1
K = a factor relating velocity ratio and circulation, Ka
1
a

= radial velocity component in potential, or nearly irrotational flow region

[/
jU =

pressure

= tangential velocity component, V(r, 6); in potential or

= complex velocity, W = U + iV, also velocity magnitude

nearly irrotational flow region

aUw

(Continued on next page)

Journal of Basic Engineering

SEPTEMBER

Copyright 1961 by ASME


Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

19 6 1 /

461

cal process since no mathematical methods were available to express the manner in which friction between the fluid and the
rotating cylinder would produce circulation.
The first quantitative data on record were determined by Lafay
in Paris between 1910 and 1912 [3, 4]. He also obtained negative
lift forces at low rotational speeds.
Many attempts have been made to use the high lift forces o b tainable on a spinning cylinder in the wind, or in an air stream,
but none has ever led to a financially successful venture. The
most notable attempt was made by Anton Flettner in Germany
during the period after the First World War. Flettner consulted
with Prandtl and the Gottingen group on the idea of replacing a
ship's sails with rotors. In a cross wind, the Magnus effect would
produce a thrust many times that for equivalent sail area. It
was, of course, necessary to drive the rotors, but the power was
only a small fraction of the power required for screw propulsion.
Ackeret and Busemann [5, 6, 7] conducted a series of tests on
cylinders with and without end plates which indicated that the
method was feasible. T w o Flettner Rotor ships, the Buckau and
the Barbara [8, 9, 10, 11) were built and one sailed across the
Atlantic. Although this method of propulsion was quite inexpensive, the speed and reliability of screw propulsion was more
than competitive.
The Madaras Rotor Power Plant, patented in 1926, was proposed as a rather complicated method of generating power on a
large scale using large vertical spinning rotors propelling themselves around a round or oval track by means of the Magnus
thrust. Generators attached to the wheels generated the power
output.
Attempts have been made to employ spinning cylinders as
portions of high lift wings, but the difficulties of complicated
drives and excessive weight have eliminated this application.
The first Magnus investigation in this country was made by
E. G. Reid in 1924, at the Langley Field N A C A Laboratory [12],
using a single cylinder projecting through both sides of the fiveoot diameter tunnel. The most interesting of these results were
never published.
The most complete experimental work was done by A. Thom
at the University of Glasgow and reported in his doctor's dissertation and in five Reports and Memoranda of the British Aircraft Research Council covering a period of nine years from 1925
to 1935 [13-18]. The effects of Reynolds number, surface condition, aspect ratio, end conditions, etc., were investigated.
Pressure, velocity, and circulation data were also obtained.
A comparison of the Magnus lift force coefficient as a function
of peripheral to free-stream velocity ratio as obtained by several
investigators under a number of conditions is shown in Fig. 1.
One general trend indicated by these data is that due to aspect
ratio: A finite and smaller aspect ratio produces a smaller lift
at a given velocity ratio and also gives a peak lift which has not
been obtained with an infinite aspect-ratio cylinder.

7T

/'d

Fig. 1 S u m m a r y of previous lift-coefficient


a data

Reynolds
no.

Aspect
ratio

Curve investigator a n d
ref.

CL v e r s u s

velocity-ratio

Remarks

( I d e a l fluid)

oo

CO

Thom [ 1 8 ]

12.5
and 2 6

5.3 103

8.8

Reid [ 1 2 ]

13.3

3.9 104

1 1.6

(Gottingen) [21]

4.7

5.2 X

104

Thom [ 1 4 ]

104

Swanson

OO

1.6 X

3.5 X 104 X 105

Thom [ 1 7 ]

5.7

3 -

9 X

104

Rough ( " s a n d e d " )


surface

Thom [ 1 7 ]

5.7

3 -

9 X

104

Smooth s u r f a c e

(Gottengen) [37]

4.7

5.2 X

Schwartzenberg

4.5

5.4,18.6 X

Swanson

5 X

X cyl-dia
disks

end

1.7 X cyl-dia end


disks
3

104

104

104

Unpublished (Case
Inst.)
Continuous
sections

end

-NomenclatureX

X
y
Y

horizontal Cartesian co-ordinate

velocity ratio, a =

= argument of location of external vortex, z =

r =

cely

circulation

= general angular co-ordinate in z = re'e measured from


positive x axis which is parallel to the free stream flow
direction at x = <*>
vortex strength, circulation

462 /

Subscripts

v0/U

p = density
<p = stream function
co = angular velocity

= force component in ^-direction


= vertical Cartesian co-ordinate
= force component in indirection
complex variable, z = x + iy

SEPTEMBER

196 1

0 (zero) denotes conditions at surface of cylinder


r - r0)
co (infinity) denotes conditions far upstream

(i.e., at

Superscripts

' (prime) denotes and implies Cartesian co-ordinates, e.g., u' is


the velocity component in the x-direction
denotes complex conjugate

Transactions of the AS ME

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

One of the first attempts at obtaining an analytical representation was made by W. G. Bickley [19] in 1928. This is probably
one of the most useful analyses with respect to providing an explanation for some of the observed flow phenomena. Bickley
considered the potential flow resulting from a vortex in the neighborhood of a cylinder with circulation in a free stream.
The most complicated mathematical model was formulated
and worked out by Torsten Gustafson [20] and published in 1933.
In a certain sense it may be considered as an extension of Bickley's method. Instead of a single vortex used to represent the
vorticity shed in the wake, " . . . a flow where vortices of constant
intensity are distributed on the logarithmically deflected streamlines in the wake" is considered. This method is based on
Oseen's methods of approximate analyses of viscous flows as extended by Zeilon. Gustafson's method is quite tedious and requires numerical solution as well as determination of several
parameters from experimental data.
A recent work by E. Krahn [21] considers the very low
Reynolds number region in which no wake is formed. A laminar
boundary-layer analysis produces the circulation as a function of
rotational speed. The method is carried through only for the
case where the t wo stagnation points coincide at the bottom of
the cylinder.
The foregoing investigations have been carried out primarily on
the case of a cylinder normal to the flow direction. The Magnus
force and tumbling moment produced on spinning projectiles
and missiles at small yaw angles are also of interest and have been
the subject of much investigation [22, 23, 24, 25, 26, 27, 28, 29].
The Magnus effect is one of the primary factors contributing to
the drift and dispersion of spinning missiles.
More recently Glauert [30, 31, 32] and Moore [33] have made
some analyses for the small perturbation case of low peripheral to
free-stream velocity ratio.

Experimental Measurements of the Two-Dimensional


Magnus Forces

cylinders running on shafts concentric with the main shaft. All


three sections were spun simultaneously using couplings that
transmitted torque, but negligible transverse thrust. A very
close clearance (0.010 to 0.015 in.) was maintained between the
6-in. diameter cylinder sections. The dummy cylinders were also
extended through the wind-tunnel walls with a close clearance to
obtain minimum end effects. Data were obtained for as broad a
range of Reynolds numbers (based on free-stream velocity and
cylinder diameter) and velocity ratios as was possible with the
apparatus. The force coefficient data are presented in Figs. 2, 3, 4,
and 5.
One of the primary objectives of this investigation was to determine whether or not a maximum (peak) lift coefficient indicated by Prandtl [35] would be obtained for an infinite aspectratio cylinder. None was obtained and it can be seen that the
Magnus lift was still increasing uniformly at a velocity ratio of 17.
Contrary to this behavior, Prandtl had predicted that a maximum
in lift coefficient would exist at a velocity ratio of about 4. His
prediction was based on the belief that the stagnation points
would coincide at the bottom of the cylinder at this velocity
ratio, beyond which no more vorticity could be shed; the lift
would then remain constant at a value near that corresponding
to an ideal velocity ratio of 2, i.e., C L u m x = 4ir. His prediction
was based on his well-known observation and photographs [36,
37] of flow patterns which indicated that the real fluid flow pattern for a velocity ratio of 4 corresponded to the ideal fluid flow
at 2. The apparatus with which these observations were made
had two restrictions which have subsequently proved them to
produce incorrect deductions regarding the actual flow patterns:
First, the patterns are of the flow on a free surface which is at a
uniform pressure and, second, the Reynolds number (free-stream
velocity, cylinder diameter) was in the neighborhood of 4000.
A comparison of Prandtl's flow patterns with those obtained by
Prof. F. N. M. Brown [33] at Notre Dame University shows a
striking difference. In the actual flow with a wake behind the
cylinder, the flow patterns will be much different from those
obtained by Prandtl on which he made his hypothesis. One interesting fact shown by the smoke tunnel patterns obtained by
Brown is the lack of rotation of the forward stagnation point predicted by the ideal flow representation of the flow by a doublet
and single lifting vortex in a free stream. This phenomenon will
be discussed further.

As previously mentioned, Lafay was the first to make quantitative measurements of Magnus forces. He also made measurements of the pressure distribution around the rotating cylinder.
The extensive data gathered by Thom for several end conditions
including various end shapes and combinations of end disks give
a variety of data for many basic shapes for the case of 90-deg
yaw. These results along with those of other investigators are of
primary interest in indicating the effect of finite aspcct ratio.
The smaller the aspect ratio, the smaller the maximum lift obtained and the smaller is the velocity ratio at which this maximum
is reached. Leakage flow and consequent pressure equalization
around the ends of the cylinder is responsible for this aspect-ratio
effect. The best approximation to two-dimensional flow (i.e.,
independent of end effects) is obtained by extending the cylinder
through the tunnel walls with a very small clearance. The only
investigation using such an apparatus seems to have been made by
Reid [12] in 1924. Many attempts have been made to approach
two-dimensional or infinite cylinder end conditions by adding end
disks. End disks, however, give entirely different flow conditions
from those of an infinite aspect-ratio cylinder so that the flow
pattern will be one due to a finite cylinder plus that due to the end
disks resulting in secondary axial flows and end flows around
the disks.
Under such circumstances, no combination of disks
on a finite cylinder woidd be expected to produce conditions
similar to those for ail infinite cylinder.

The variation of the lift and drag at low Reynolds numbers and
velocity ratios is also interesting. The variation in behavior with
Reynolds number is attributed to the nonsymmetric boundarylayer separation from the top and bottom surfaces of the cylinder
as explained by Krahn [21] and Swanson [34]. As one example,
consider curve g of Fig. 3 for which the Reynolds number is subcritical (with respect to laminar-turbulent boundary-layer transition) at 1.52 X 105.
At low cylinder to free-stream velocity
ratios, a, between 0 and 0.1, both the upper and lower boundary
layers remain fully laminar up to their separation points. The
upper separation point moves rearward (in the direction of the free
stream) and the lower separation point moves forward. The additional length of boundary attachment on the top surface gives a
greater region of negative pressure coefficient here than over the
bottom where the length of attached boundary layer has decreased. As a consequence of the pressure distribution, a positive
lift is produced.

The closest approach to infinite cylinder conditions is believed


to have been obtained with a three-section apparatus used by the
author [34]. A "live" cylinder section mounted on a long shaft
supported by cantilever strain-gage beams was flanked by dummy

2 This CD versus R plot for the static cylinder is of interest in that


it indicated the degree of turbulence in the tunnel. Also, the transition region of what is usually referred to as that where the drag coefficient suddenly decreases actually covers about three quarters of
the Reynolds number range up to the point of minimum CD-

Journal of Basic Engineering

The boundary layer on the stationary cylinder at this Reynolds


number is already slightly transitional as is indicated on Fig. 6.2
As a is considered to further increase beyond 0.1, the relative

SEPTEMBER

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

19 6 1 /

463

/
/

F
a

\p

' A

<

"N.

4
.5
Fig. 2

Group 1

Group 2

Group

SEPTEMBER

1 9 6 1

.8

.6

.9

1.0

Lift c o e f f i c i e n t Cl v e r s u s v e l o c i t y r a t i o a

Curve designation

464

"""

Reynolds
3.58
4.9
6.07
7.91
9.9
12.8
15.2
18.15
22.5
26
29.5
32.5
36.5

no.
104

42
45
50.1

Transactions of the AS ME

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

a
F i g . 3 D r a g c o e f f i c i e n t Co
R e y n o l d s n u m b e r s as Fig. 2 .

v e r s u s v e l o c i t y r a t i o a.

Letter c u r v e d e s i g n a t i o n s

imply

same

Experimental

V,

(Eouation ( 7 )
\ x a = 2.89, c/a = 3.51, T = 242.1
(Eouation ( 7 )
( K a = 2 . 8 9 , c/a = 4 , y 224

>

/
/
0

7
/

\
j

i
\

Fit;. 4

\
9

10

11

Lift c o e f f i c i e n t C i v e r s u s v e l o c i t y r a t i o

12

15

14

15

16

17

a.

I
1/

//

V
N

f ,
ii / I
/

i
i t
/J
w
/

Experimental

(Equation ( 7 )
.
l K a = 2.89, c / a - 3 . 5 1 , Y = 242.1
(Eouation ( 7 )

U 0 - 2 . 8 9 , c/a = 4, r = 224

./
0

3
Fig. 5

10

11

12

13

14

D r a g coefficient Co versus v e l o c i t y ratio a .

Journal of Basic Engineering

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

SEPTEMBER

I 96I

/ 465

5
B x

10"

F i g . 6 D r a g c o e f f i c i e n t C o v e r s u s R e y n o l d s n u m b e r R p l o t t e d o n a l i n e a r s c a l e . T h e letters
a, c, e, e t c . , c o r r e s p o n d t o t h e f r e e - s t r e a m R e y n o l d s n u m b e r s f o r t h e c u r v e s s h o w n i n F i g . 2 .

velocity over the top of the cylinder is decreased, but it is increased over the lower half. It is useful to introduce here the concept of a "relative Reynolds number." Where the surface is
traveling with the free-stream velocity, the transition can be expected to be delayed, and conversely. The relative velocity R r o i
= R ( l a) ( + sign for bottom and sign for top). The
boundary-layer behavior can then be correlated with the relative
Reynolds number and a qualitative estimation of the relative portions of laminar and turbulent boundary layer and total length
of boundary layer can be obtained. Referring back to curve g
of Fig. 3, for a increasing from 0.1, the upper relative Reynolds
number decreases so that conditions in the upper boundary layer
can be considered to be similar to those for Reynolds numbers
less than 1.5 X 10E on Fig. 6. The lower boundary-layer conditions correspond to Reynolds numbers greater than 1.5 X 105.
For example, for a = 0.2, the upper relative Reynolds number
would be 1.5 X 10s (0.8) = 1.2 X 106 and the lower would be 1.8 X
106. Fig. 6 then indicates a longer attached boundary layer on the
bottom than on the top. This effect counteracts that previously
mentioned and the lift begins to decrease. A t higher rotational
speeds, a greater portion of the lower boundary layer becomes
turbulent and therefore more reattached. The greater reattachment on the lower surface produces the negative lift effect shown
and. also results in a lower drag as shown in Fig. 3. The lower
boundary layer will finally reach a fully developed turbulent state
near the point of maximum negative C L .
At this point it is helpful to introduce the concept of a "boundary-layer origin.'' The positions of the boundary-layer instability
and separation points are a function of a boundary-layer length
Reynolds number. For a stationary surface such as a flat plate,
airfoil, or nonrotating cylinder, the boundary-layer length is
measured from the front stagnation point. The shear (or rotation) in the boundary layer has opposite direction (or is of opposite sign) on opposite sides of the body from this point (Fig. 7).
It seems logical to define the beginning of the boundary layer, or
a boundary-layer origin, with respect to such a condition if we
accept the definition and idea of a boundary layer as a shear layer.
For a stationary surface, this boundary-layer origin is coincident
with the stagnation point; however, for a moving surface this
boundary-layer origin will, in general, not coincide with the stagnation point. On the rotating cylinder, the forward stagnation
point is translated in a direction opposite to the direction of rotation. (In addition, it is interesting to note that the stagnation
point as a point of zero velocity can no longer lie on the surface
466 /

SEPTEMBER

Fig. 7 Representation of b o u n d a r y layers a r o u n d rotating cylinder for


a = 0.2 and R
X I04

of the body.) If the boundary or shear layers are now investigated, we find that the upper and lower layers conforming to the
required conditions for a boundary layer begin at a point where
the cylinder surface velocity is equal to the fluid velocity (i.e.,
where the relative velocity is zero). The boundary layers are
then as shown in Fig. 7. It is seen that the boundary-layer
origin is translated in the same direction as the direction of motion of the surface (direction of rotation) but opposite to that of
the stagnation point. Boundary-layer profiles obtained for the
cylinder rotating at two different velocity ratios are presented in
Figs. 8 and 9. The outlines of the boundary layers are shown.
It is interesting to note that, for a = 1, both boundary layers are
about the same length. In order to calculate characteristics of the
boundary layer, it is necessary to have a set of initial conditions.
The previously described boundary-layer origin is the most logical
position at which to attempt to start a boundary-layer solution
for a geometry where the wall is in a nonuniform motion with
respect to the external flow.
As the rotational speed is further increased (beyond a value
corresponding to a ~ 0.5) the boundary-layer origin moves
further back on the upper surface with a consequent increased
length of flow attachment producing more Magnus lift. It is
believed that both boundary layers are in the fully developed
turbulent state at a velocity ratio of 1.0. An idea of the flow
pattern may be obtained by referring to Fig. 8. Although the
relative velocity between the surface and free stream is much

196 1

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

Transactions of the AS ME

Fig. 8 B o u n d a r y - l a y e r p r o f i l e s o b t a i n e d a t p o s i t i o n s e v e r y 3 0 d e g a r o u n d
the cylinder for a = 1. T h e dotted line represents t h e a p p r o x i m a t e limit
f o r t h e t h r e e b o u n d a r y l a y e r s . T h e v e r t i c a l s c a l e o n t h e p r o f i l e s is t e n
t i m e s t h a t f o r t h e c y l i n d e r . S c a l e s a r e v/toro v e r s u s y / r o .

Fig. 9

S a m e as Fig. 8 ; a =

smaller in the upper boundary layer than in the lower boundary


layer, the velocity gradients appear to be about the same at equal
distances from the boundary-layer origin. The main indication
that both the layers are fully developed turbulent at and beyond
this point lies in the fact that the drag in this region of a is a
minimum corresponding to the fully developed turbulent boundary layer over the nonrotating cylinder at R ~ 3.5 X 10s; in
addition, no further dependence of either CL or CD on R is indicated beyond a = 1. The foregoing is only a very qualitative
description of the assumed flow behavior in the light of the experimental results available.
Since no hot-wire data have been
taken within the boundary layer, the foregoing assumptions have
not been verified.

Journal of Basic Engineering

For velocity ratios greater than one, the force behavior is indicated to be essentially independent of Reynolds number for 2 X
10' < R < 5 X 106. In this region, both boundary layers are
apparently fully developed. Since the surface and fluid velocities
over the top are in the same direction producing a low relative
velocity, it is possible that the upper boundary layer in its fully
developed state at a > 1 is laminar. However, boundary-layer
instability and transition, and the transition from laminar to
turbulent flow, in general, are functions of Reynolds number based
on absolute velocity (considering the Reynolds number and its
consequences to be representative of the ratio of inertial to viscous
forces). Over the upper surface the absolute velocity is quite
high so that, for free-stream Reynolds numbers corresponding to
laminar separation for the static cylinder, it is possible that the
upper boundary layer is turbulent for the case of the rotating
cylinder. One of Brown's photographs for 1 < a < 2 (Fig. 14,
ref. [36]) favors the second argument: The smoke filament on
the upper surface indicates a turbulent boundary layer. Similar
arguments can be presented for the behavior of the lower boundary layer.
The knee of the lift curve near a 3 is produced when the
boundary-layer origin reaches the top of the cylinder. At any a
greater than this value, the cylinder is everywhere traveling at a
velocity greater than that attainable by the external flow; consequently the boundary-layer origin as defined previously no longer
exists. Prandtl considered that the flow pattern produced by the
stagnation point coincidence at the bottom of the cylinder would
produce a flow pattern and force system that would not change
by any further increase in rotational speed. He reasoned that no
further vorticity could be shed; therefore the CL would level off
at a value of 4ir. The foregoing condition where the boundarylayer origin reaches the top may be considered to pose a similar
limitation, except that the CL is seen to continuously increase, but
at a slower rate than before. As the rotational speed is further increased beyond a value corresponding to a = 3, the vorticity shed
in the upper boundary layer goes from positive (clockwise) to
negative (counterclockwise) while that shed in the lower boundary layer is always negative. In a real fluid, the lift and shed
vorticity will always have opposite signs (for the co-ordinate
system orientation shown) and there will be a correspondence between the magnitude of the shed vorticity and the lift. As a
is further increased, the vorticity shed in both boundary layers
is negative and increases in magnitude. The field induced by this
shed vorticity has the effect of rotating the flow pattern around
the cylinder in a counterclockwise direction. (An approximate
calculation of this effect will be given later.) The pressure distribution produced by this pattern is seen to increase the lift. If
this is the primary effect producing the lift rise beyond the knee,
it appears as if some limit should exist as the wake formed by the
separating boundary layers rotates around to a position near the
front of the cylinder. N o such limiting condition was reached
within the limits of data taken.
The drag behavior is also indicated as a consequence of the
foregoing flow behavior. As a increases beyond 1, the drag, surprisingly, increases to a value much greater than the drag on the
nonrotating cylinder, even though the wake area is decreasing.
This large drag increase with increasing a is produced by flow reattachment over the rear of the cylinder accompanied by a movement of the rear stagnation point and the wake in a counterclockwise direction into the region near the bottom of the cylinder.
The drag peaks in the region where the lift knee occurs. The
boundary-layer origin is at the top of the cylinder and the separation points and wake are near the bottom of the cylinder. An increase of a as described produces a further rotation of the
wake toward the front of the cylinder. The resultant flow pattern
and pressure distribution produce a decrease in CD.

SEPTEMBER

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

19 6 1 /

467

A Proposed Ideal Flow Analysis

dF
W > - U > - i V > - - - - U

The potential flow pattern of a free stream, doublet, and bound


vortex combined to yield the flow on a cylinder with lift is quite
familiar. This simple potential function and its solution yield a
lift, CL = 2ira (a = v0/U), but no drag is predicted.
As one step beyond the simple Kutta-Joukowski theorem of
lift, consider the force on the cylinder produced by the circulation
from a singularity (potential vortex T) within the boundary as in
the Kutta-Joukowski theory and an additional vortex downstream
from the cylinder representing the shed vorticity in the wake that
is continuously being carried downstream. This moving vortex
sets up an induced velocity field which acts to produce a force in
addition to that arising from the interaction of the free-stream
velocity and the circulation within the body contour.
In a very interesting paper, W. G. Bickley [19] has set up an
approximate mathematical model for this problem in which the
net vorticity moving in the wake is considered to be concentrated
in a single vortex of opposite direction to the lift-circulation
around the cylinder (Fig. 10). Expressions for the transportation
velocity components of the "shed" vortex are given and from
these and the velocity potential for this flow model, the lift and
drag are determined for an arbitrarily prescribed location of the
shed vortex.

_ t(r +

K)

t
(l

jk /

a2 \
-)

i_\

2ir \ 2 zA

2TZ

z zB )

For ideal motion of an incompressible flow the Blasius theorem


yields the force on the cylinder boundary:

d<f> -

~dT dz

The first contour integral on the right-hand side of (3) gives the
convective or W2 contribution to the force and the second term
gives the unsteady or time-dependent contribution. The contour
of integration is the cylinder surface, z = a.
After performing the necessary calculations, the lift and drag
coefficients can be determined as:

S T . (' "

- -

idir.

0 ^ 0

( ' - ? )

" ^ f " H

Ctt7- -

i r ) " "

<4

The three parameters, o / c , K (or K/T), and y cannot directly be


prescribed. By associating the vorticity shed in the wake with
the lift circulation, one may arbitrarily take K = I\ This
assumption gives

>

( >

vortex"

Consider the flow around a circular cylinder of radius a with


positive circulation and superimposed on this, a vortex of strength
k at 2 = aA = ce' 7 , and the necessary vortexes + k at the inverse
a2
point (image point in the circle) z = zB = e17, and K at 2 = 0
c
to maintain the circle z = a as a streamline (Fig. 10). The vortex
at A is a free vortex; those at z = 0 are bound, and the vortex at B
is restricted to move such that it is always at the inverse point of
A. With this arrangement z = a is preserved as a streamline of
the flow as Ka moves according to the velocity induced at A by the
rest of the field. This potential flow representation is an approximation to the actual flow where the net vorticity continuously shed by the cylinder and washed downstream is considered
to be instantaneously concentrated with strength K at A, and
with an instantaneous velocity as determined from the rest of the
field. The net vorticity within the cylinder contour is the " l i f t "
circulation.
The complex potential function for the flow of Fig. 10 is
,

Tr

a2 \

i(T + k)

in

2 - ZA

(1)
a2

where z., .= ce,y and zB = ely.


c2

468 /

SEPTEMBER

196 1

The complex velocity is given

(f)

a y

The angle 7 is in the third quadrant so that CD is a positive


quantity and the second term in CL is negative.
The most interesting characteristic of this solution is the prediction of a drag that increases with the square of the circulation.
The lift is now seen to be parabolic with T, instead of linear as predicted b y the Kutta-Joukowski theory.
A polar CL versus CD relation is obtained by eliminating T:

= [ ' " (7)']

<6'

which is seen to be a family of parabolas with parameters cja


and 7 . Fig. 11 shows a plot of CL versus CD for the experimental
data and for this relationship with c / o = 4 and 7 = 210 deg.
(These values were arbitrarily chosen to give the best fit of the
foregoing equation with the experimental results.) Curves for two
other sets of parameters are also shown. The lift and drag curves
for the experimental data and for these sets of parameters in the
foregoing expressions for CL and CD are shown on Figs. 4 and 5.
The circulation T (and k) is a function of a . On the KuttaJoukowski case, T/(aUa)
= 27ra; however, since this model
replaces the Kutta-Joukowski, T/(AU) is not determined, but
can only be said to be a function of A: T/(AUa)
is therefore
arbitrarily taken as equal to Ka, i.e., T/(aU m) = Kaa where Ka
is obviously always less than 27T. This substitution gives:

Transactions of the AS ME

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

(7)

cos y / a \ ,

- - ^ T (7)

Even considering that this model presents a good approximation to the actual flow condition, the parameters c/a, 7 , and Ka
would be expected to be functions of or; consequently, if these
relationships could be determined, a better " f i t " to the experimental results would result. Unfortunately, there is no way by
which c/a(a), 7 ( a ) , and Ka(a) can be determined.
It is interesting to note that the field induced by the trailing
vorticity has the effect of rotating the entire flow pattern in a
counterclockwise direction about the center of the cylinder (in
the same direction as the lift circulation T). The forward stagnation point is then located closer to the front (6 = 0) of the cylinder and the rear stagnation point is closer to the bottom (6 = 90
deg).
For example, for the values of c/a = 4, 7 = 224 deg, and
Kaa = 2.89 cited in the foregoing example, the stagnation points
for a = 1 are located at 9 deg and at 180 + 21 deg; for a = 2,
at 19 deg and 180 + 50 deg. This predicted effect is in agreement with visual observations by Brown [36] and Gustafson
[20] and measurements by Swanson [34] (Figs. 8, 9). In reality,
the rear stagnation point does not exist because of the wake region
formed. It is of significance that for the range of Reynolds numbers investigated (3 X 104 to 5 X 106) the wake is always present
at all values of a. Therefore, for a real fluid flow at high Reynolds
numbers, the condition where the two stagnation points can
coincide will not exist and a separated wake will always be present.
Indeed, some smoke studies by Gustafson [20] at very large a (5

\
\

/>

/
A

//

Experimental
(Y = 210
/ c/a - '<0
jy = 22'*
1 c/a = 3 . M
\ y = 2'l2.)

to 7) indicate that the wake has moved around to a position such


that both separation points are within the fourth quadrant.
The mathematical model described in the foregoing is of interest primarily in that it predicts some of the observation
flow pattern characteristics. It cannot be used for quantitative
predictions of the force behavior since, first, it is an ideal
(potential) flow representation of a flow that is principally determined by viscous effects.
The circulation around the actual rotating cylinder, which can
be calculated and plotted from velocity field data, is a consequence of the unsymmetrical flow pattern produced by the upper
and lower boundary layers separating at different positions.
The circulation is then a consequence of the flow pattern as determined by the boundary-layer behavior. The circulation at the
cylinder surface is, of course, T = 27roor02 which would yield a lift
coefficient corresponding to that for ideal flow (e.g., C l = 2xa for
the Kutta-Joukowski model). This value of F falls off through
the viscous (boundary-layer) region to a value that decreases very
slowly outside this region. The value of F determined experimentally at the outer edge of the boundary layer agrees closely
with that calculated from the ideal flow model.

The Fundamental Problem


The most simple problem in the range of Reynolds number
where separation exists would be that where the Reynolds number indicates fully laminar boundary-layer flow up to the point of
separation. Since the problem of separating flow on a stationary
cylinder has not yet been solved without first knowing the answer
(i.e., the actual pressure distribution), the possibility of a solution
to this problem does not appear imminent. On the other hand,
with the actual pressure distribution given, solution of the
boundary-layer equations for the case of the nonrotating cylinder
gives very good agreement with respect to velocity profile and
separation point determination. The analogous problem for the
rotating cjdinder is currently under attack.
Experimental observations of the actual flow patterns often
provide valuable or interesting information with respect to
analytical formulations. In particular, it appears misleading to
present only the mean values of data such as force coefficient
(C L and CD) and pressure distribution data. Even though the
range of Reynolds numbers discussed is much higher than that in
which a discrete vortex street is formed, the flow pattern is far
from steady or stationary. The random fluctuating wake pressure (and velocity) signals are propagated throughout the entire
flow field. The varying pressure distribution in the neighborhood
of the separation points produces relatively large variations in the
position where separation occurs. Consequently, the net force on
the body experiences large random fluctuations in both magnitude
and direction. In the case of the cylinder (rotating or not) the
vertical (lift) force was observed to vary in magnitude by an
amount on the same order as the average value of the net force.
The variations in drag force were not as great as the vertical force
variations. The fundamental problem is then one of an unsteady flow; however, solutions neglecting the random flow variations do give satisfactory results that agree with the physical
average value results obtained from experiments. A solution to
the complete problem of the flow around a bluff body still requires solution of nonlinear partial differential equations. In
the case of large Reynolds numbers, the combined or simultaneous
solutions to the boundary-value problem in the external flow and
the boundary-layer flow are required.

References

Fig. 11

L i f t - d r a g p o l a r p l o t w i t h c o m p a r i s o n s f r o m Eq. ( 6 )

Journal of Basic Engineering

1 G. Magnus, " O n the Deflection of a Projectile," (a) "Abhandlung der Akademie der Wissenschaftern," Berlin, Germany, 1852,
English translation in Taylor's Scientific Memoirs, 1853. (b) " P o g gendorffs Annalen der Physik und Chemie," vol. 88, no. 1, 1853.
SEPTEMBER

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

1 96 1 /

469

Lord Rayleigh, "On the Irregu1ar Flight. of a Tennis Ball,"


" Messenger oC

(a) "Scientific Papers," vol. 1, 1B57, pp. 344-6.


Mathematics," vol. 7, 1B77, p. 14.

3 LaCay, "Sur l'Inversion du Ph6nomene de Magnus," Comptes


Rend1l8, vol. 161, HHO.
4 Lafay, "Contribution Experiment.ale n l'Acrodynamique du
Cylindrc," Rermea Mechaniqlle, vol. 30, 1912.
5 A. Bet.z, "Del' Magnuscffckt., die Grundlage del' F lettnerWalze," VOl, 1925, p. 9. English translation, NACA TM 310,
1925.
6 L. PrandU, "Mngnuseffekt uud Windkrnftsehiff," Die Naturwissenschaften, 1925, p. 93. English translation, NACA TM 367,
1926.
7 A. Buseman, "Mcssungen an Rolicrcnden Zylindern," Ergibnisse del' Aerodynamik Versuehsanstalt zu Gottingen, IV Lieferung,
1932, p. 101.
B R. Rizzo, "The Flettner Rotor Ship in the Light of the KuttaJoukowski Theory and oC Experimental Results," NACA TN 22B,
1925.
9 H. Ziekelldraht, "Magnuseffekt, Flettner-Rot.or und DllsenFIUgel," Schweiz Tech7i. Zcitschr. , vol. I, 1926, pp. B37-843, B55-857.
10 F. Ahlborn, "Del' MagnuseO'ekt in Theorie und Wirklichkeit,"
Zeilachr. /. Flllotechm'k laid Motorlll/lschiJJahrt, Dec. 28, 1929, pp.
642- 653. English translation, NACA TM 567.
11 A. Flettner. "The Flettncr Rotor Ship," l!-'1I0illce,.i,10, vol. 19,
January 23, 1925, pp. 117- 120.
12 E. G. Reid, "Tests of Rotating Cylinders," NACA TN 209,
1924.
13 A. Thom, "The Aerodynamics of a Rotating Cylinder," thesis,
University of Glasgow, 1926.
14 A. Thorn, "Experiments ou the Air Forces on Rotating
CyHllders," ARC Rand M lOIS, 1925.
15 A. Thorn, "The Pressures Round a Cylinder Rotating in au
Ail' Current," ARC Rand M lOS2, 1926.
16 A. Thorn, "Experiments on the Flow Pnst a Rotating Cylinder." ARC Rand M 14 10,1931.
17 A. Thorn aud S. R. Sengupta, "AirTorquc on a Cylinder Rotating in an Air Stream," ARC Rand M 1520, 1932.
18 A. Thorn, "Effects of Discs 011 the Air Forces on a Rotatillg
Cylinder," ARC Rand M 1623. 1934.
19 W. G. Bickley, "The InBuence ofVorticcs Upon thc Hesistallce
Experienced by Solids Moving Through a Liquid," Proceccii1lgs, Royal
Society (London), vol. 119, scries A, 102S, pp. 146-156.
20 T. Gustafson, "On the Magnus Effect According to the Asymptotic Hydrodynamic Theory," Hakan Ollissons Buchdruckel'ei, Lund
(Sweden), 1933. NACA trans!. N-2M121, 1954.
21 E. Krahn, "Tho Laminar Boundary Layer on a Rotat.ing
Cylinder in CrossBow," NAVORD Rep. 4022, Aerobal. Res. Rep. 28B,
U. S. NOL, Maryland, June, 1955.
22 J. C . Martin, "On Magnus Effects Caused by the Boundary
Layer Displacement Thickness on Bodics of Revolution at Small
Angles of Attack," BRL Rep. 870, D ept. of the Army Prou. 5B03-03~Ol, Ord. Res. and Dev. Proj. TB3-0108, Aberdeen Proving Ground,
Maryland, June, 1955.
23 Karl H. Stefan , "Magnus Effect on a Body of Revolution at.
Various Angles of Att.ack," Master's t.hesis, Dept. of Aero. Engr ..
University of Minn., 1949.
2'1 William E. Buford, "Magnus Effect in the Case of Rotating
Cyli nders and Shell," BRL Memo Report No. 821, July, 1954.
25 Wind Tunnel Group, COil vail' Wind Tunnel Handbook, vol. L,
COil vail' Repl. 2T-043, April, 1955.
26 John W. MaeeoD, "Aerodynamics of a Spin ning Sphere,"
Jour1lal a/the Royal Aeronautical Society, vol. 32, 1928, p. 777.
27 J. M. Davies, "The Aerodynamics o f Golf Balls," Journal 0/
Applied Physica, vol. 20, 1949, p. 821.
28 R. H. Heald, J. G. Logan, Jr., H. Spivak. and W. Squire,
" Aerodynamic Characteristics of a Hot.ating Model of t he 5.0-Inch

470 I

SEPTEMBER

1961

Spin-Stabilized Rocket Model 32," U . S. Navy BuOrd Tech. Note


No. 37, June, 1957.
29 M. J. Queijo Rnd H. S. Fletcher, "Low Speed Experimental
Investigation of the Magnus Effect on VarioW! Sections of 11 Body of
R evoluLion With and Without a Propeller," NACA TN 4013, August,
1957.
30 M. B. Glauert, "A Boundary Layer Theory With Application
to Rotating Cylinders," Joumal 0/ Flltid .M echa1l'lca, yo l. II, Part I,
1957, p. 89.
31 M. B. Glauert, "The Laminar Boundary Layer Oil Oscillating
Plates and Cylinders," JOllnlal 0/ Fluid M echallic8, \'01. I, Part I,
May, 1956, p. 97.
32 M. B. Glauert, "The l~low Past a Rapidly Rotating Cylinder,"
Proceedi1l01J 0/ ths Royal Society, vol. 242, series A, 1957, pp. 108-115.
33 D. W. Moore, "The Flow Past a Rapidly Rotating Circular
Cylinder in a Uniform Stream," Joltr1lal 0/ Fflt'id J[echa l~ics, vol. II,
Part 6, 1957, p. 541.
34 W. M. Swanson, "An Experimental I nvestigation of tbe Magnus Effect," Final Report, OOR Proj. No. 1082, Case Inst. of Tech.,
Deecmber, 1956.
35 S. Goldstein , "Modern Development.s in Fluid Dynamics,"
vols. 1 and 2, Clarendon Press, Oxford, England, 1938.
36 R. W. Van Aken a nd H. R. Kelly, "The Magnus Effect on
Spinning Cylinders, " Reprint No. 712, paper presented at the 25th
Annual Meeting, LA.S., 1957.
37 L. Prandtl and O. Tietjens, "Applied Hydro and Aerody.
namics," English translation by J. P. den Hartog, McGraw-Hili Book
Company, Inc., New York, N. Y ., 1934.

Transactions 01 the AS ME

Downloaded From: http://fluidsengineering.asmedigitalcollection.asme.org/ on 06/08/2013 Terms of Use: http://asme.org/terms

You might also like