Viv Mit
Viv Mit
Viv Mit
13.42 Lecture:
Vortex Induced Vibrations
Prof. A. H. Techet
21 April 2005
Offshore Platforms
2
Fixed Rigs Tension Leg Platforms
Spar Platforms
3
Genesis Spar Platform
VIV Catastrophe
If neglected in design, vortex induced vibrations can prove
catastrophic to structures, as they did in the case of the Tacoma
Narrows Bridge in 1940.
4
I n another city, the John Hancock tower wouldn't be anything special -- just another
reflective glass box in the crowd. But because of the way Boston and the rest of New
England has grown up architecturally, this "70's modern" building stands out from the
rest. Instead of being colonial, it breaks new ground. Instead of being quaint, it soars
and imposes itself on the skyline. And Instead of being white like so many buildings in
the region, this one defies the local conventional wisdom and goes for black. For these
reasons and more the people of Boston have fallen in love with the 790-foot monster
looming as the tallest building in New England at the time of its completion. In the
mid-1990's, The Boston Globe polled local architects who rated it the city's third best
architectural structure. Much like Boston's well-loved baseball team, the building has
had a rough past, but still perseveres, coming back stronger to win the hearts of its
fans. The trouble began early on. During construction of the foundation the sides of the
pit collapsed, nearly sucking Trinity Church into the hole. Then in late January, 1973
construction was still underway when a winter storm rolled into town and a 500-pound
window leapt from the tower and smashed itself to bits on the ground below. Another
followed. Then another. Within a few weeks, more than 65 of the building's 10,344
panes of glass committed suicide, their crystalline essence piling up in a roped-off area
surrounding the building. The people of Bean Town have always been willing to kick a
brother when he's down, and started calling the tower the Plywood Palace because of
the black-painted pieces of wood covering more than an acre of its faade. Some
people thought the building was swaying too much in the wind, and causing the
windows to pop out. Some thought the foundation had shifted and it was putting stress
of the structural geometry. It turns out the culprit was nothing more than the lead
solder running along the window frame. It was too stiff to deal with the kind of
vibrations that happen every day in thousands of office buildings around the world. So
when John Hancock Tower swayed with the wind, or sighed with the temperature, the
windows didn't and eventually cracked and plummeted to Earth. It cost $7,000,000.00
to replace all of those panes of glass. The good news is, you can own a genuine piece
of the skyscraper. According to the Globe, the undamaged sheets were sold off for use
as tabletops, so start combing those garage sales. For any other skyscraper, the
hardship would end there. But the Hancock building continued to suffer indignities.
The last, and most ominous, was revealed by Bruno Thurlimann, a Swiss engineer who
determined that the building's natural sway period was dangerously close to the period
of its torsion. The result was that instead of swaying back-and-forth like a in the wind
like a metronome, it bent in the middle, like a cobra. The solution was putting a pair of
300-ton tuned mass dampeners on the 58-th floor. The same engineer also determined
that while the $3,000,000.00 mass dampeners would keep the building from twisting
itself apart, the force of the wind could still knock it over. So 1,500 tons of steel braces
were used to stiffen the tower and the Hancock building's final architectural indignity
was surmounted.
Reprinted from http://www.glasssteelandstone.com/BuildingDetail/399.php
John Hancock Building
Classical Vortex Shedding
Von Karman Vortex Street
l
h
Alternately shed opposite signed vortices
5
Potential Flow
U() = 2U
sin
P() = 1/2 U()2 = P
+ 1/2 U
2
Cp = {P() - P
}/{1/2 U
2
}= 1 - 4sin
2
s
= 2 f
s
DRAG = D(t) = Do cos (2
s
t+ )
Heave Motion z(t)
2
( ) cos
( ) sin
( ) cos
o
o
o
z t z t
z t z t
z t z t
=
=
=
&
&&
11
Lock-in
A cylinder is said to be locked in when the frequency of
oscillation is equal to the frequency of vortex shedding. In this
region the largest amplitude oscillations occur.
v
= 2 f
v
= 2 S
t
(U/d)
n
=
k
m + m
a
Shedding
frequency
Natural frequency
of oscillation
Equation of Cylinder Heave due
to Vortex shedding
Added mass term
Damping If L
v
> b system is
UNSTABLE
k
b
m
z(t)
( ) mz bz kz L t + + = && &
( ) ( ) ( )
a v
L t L z t L z t = + && &
( ) ( ) ( ) ( ) ( )
a v
mz t bz t kz t L z t L z t + + = + && & && &
{
( ) ( ) ( ) ( ) ( ) 0
a v
m L z t b L z t k z t + + + = && &
14243 14243
Restoring force
12
LIFT FORCE:
Lift Force on a Cylinder
( ) cos( )
o o
L t L t = +
v
if <
( ) cos cos sin sin
o o o o
L t L t L t =
2
cos sin
( ) ( ) ( )
o o o o
o o
L L
L t z t z t
z z
= + && &
where
v
is the frequency of vortex shedding
Lift force is sinusoidal component and residual force. Filtering
the recorded lift data will give the sinusoidal term which can
be subtracted from the total force.
Lift Force Components:
Lift in phase with acceleration (added mass):
Lift in-phase with velocity:
2
( , ) cos
o
a o
L
M a
a
=
sin
o
v o
L
L
a
=
Two components of lift can be analyzed:
(a = z
o
is cylinder heave amplitude)
Total lift:
( ) ( ( , ( ) ( ) , ) )
a v
L t z t L a M z a t = + && &
13
Total Force:
If C
Lv
> 0 then the fluid force amplifies the motion
instead of opposing it. This is self-excited
oscillation.
C
ma
, C
Lv
are dependent on and a.
( ) ( ( , ( ) ( ) , ) )
a v
L t z t L a M z a t = + && &
( )
( )
2
4
2
1
2
( , ) ( )
( , ) ( )
ma
Lv
d C a z t
dU C a z t
=
+
&&
&
Coefficient of Lift in Phase with
Velocity
Vortex Induced Vibrations are
SELF LIMITED
In air:
air
~ small, z
max
~ 0.2 diameter
In water:
water
~ large, z
max
~ 1 diameter
14
Lift in phase with velocity
Gopalkrishnan (1993)
Amplitude Estimation
=
b
2 k(m+m
a
*
)
m
a
*
= V C
ma
; where C
ma
= 1.0
Blevins (1990)
a
/
d
=
1.29
/[1+0.43 S
G
]
3.35
~
S
G
=2 f
n
2
2m (2)
d
2
; f
n
= f
n
/f
s
; m = m + m
a
*
^
^
_
_
15
Drag Amplification
Gopalkrishnan (1993)
C
d
= 1.2 + 1.1(
a
/
d
)
VIV tends to increase the effective drag coefficient. This increase
has been investigated experimentally.
Mean drag: Fluctuating Drag:
C
d
occurs at twice the
shedding frequency.
~
3
2
1
C
d
|C
d
|
~
0.1 0.2 0.3
fd
U
a
d
= 0.75
Single Rigid Cylinder Results
a) One-tenth highest
transverse
oscillation amplitude
ratio
b) Mean drag
coefficient
c) Fluctuating drag
coefficient
d) Ratio of transverse
oscillation frequency
to natural frequency
of cylinder
1.0
1.0
16
Flexible Cylinders
Mooring lines and towing
cables act in similar fashion
to rigid cylinders except that
their motion is not spanwise
uniform.
t
Tension in the cable must be considered
when determining equations of motion
Flexible Cylinder Motion Trajectories
Long flexible cylinders can move in two directions and
tend to trace a figure-8 motion. The motion is dictated by
the tension in the cable and the speed of towing.
17
Shedding patterns in the wake of oscillating
cylinders are distinct and exist for a certain range
of heave frequencies and amplitudes.
The different modes have a great impact on
structural loading.
Wake Patterns Behind
Heaving Cylinders
2S 2P
f , A
f , A
U
U
Transition in Shedding Patterns
W
i
l
l
i
a
m
s
o
n
a
n
d
R
o
s
h
k
o
(
1
9
8
8
)
A/d
f* = fd/U
Vr = U/fd
18
Formation of 2P shedding pattern
End Force Correlation
Uniform Cylinder
Tapered Cylinder
H
o
v
e
r
,
T
e
c
h
e
t
,
T
r
i
a
n
t
a
f
y
l
l
o
u
(
J
F
M
1
9
9
8
)
19
VIV in the Ocean
Non-uniform currents
effect the spanwise vortex
shedding on a cable or
riser.
The frequency of shedding
can be different along
length.
This leads to cells of
vortex shedding with
some length, l
c
.
Three Dimensional Effects
Shear layer instabilities as well as longitudinal (braid)
vortices lead to transition from laminar to turbulent flow in
cylinder wakes.
Longitudinal vortices
appear at R
d
= 230.
20
Longitudinal Vortices
C.H.K. Williamson (1992)
The presence of
longitudinal vortices leads
to rapid breakdown of the
wake behind a cylinder.
Longitudinal Vortices
21
Strouhal Number for the tapered
cylinder:
S
t
= fd / U
where d is the average
cylinder diameter.
Oscillating Tapered Cylinder
x
d(x)
U
(
x
)
=
U
o
Spanwise Vortex Shedding from
40:1 Tapered Cylinder
T
e
c
h
e
t
,
e
t
a
l
(
J
F
M
1
9
9
8
)
d
max
R
d
= 400;
St = 0.198; A/d = 0.5
R
d
= 1500;
St = 0.198; A/d = 0.5
R
d
= 1500;
St = 0.198; A/d = 1.0
d
min
No Split: 2P
22
Flow Visualization Reveals:
A Hybrid Shedding Mode
2P pattern results at
the smaller end
2S pattern at the
larger end
This mode is seen to
be repeatable over
multiple cycles
Techet, et al (JFM 1998)
DPIV of Tapered Cylinder Wake
2S
2P
Digital particle image
velocimetry (DPIV)
in the horizontal plane
leads to a clear
picture of two distinct
shedding modes along
the cylinder.
Rd = 1500; St = 0.198; A/d = 0.5
z
/
d
=
2
2
.
9
z
/
d
=
7
.
9
23
Evolution of the
Hybrid Shedding Mode
2P 2S
Rd = 1500; St = 0.198; A/d = 0.5
z/d = 22.9 z/d = 7.9
Evolution of the
Hybrid Shedding Mode
2S 2P
Rd = 1500; St = 0.198; A/d = 0.5
z/d = 22.9 z/d = 7.9
24
Evolution of the
Hybrid Shedding Mode
2S 2P
Rd = 1500; St = 0.198; A/d = 0.5
z/d = 22.9 z/d = 7.9
Evolution of the
Hybrid Shedding Mode
2S 2P
Rd = 1500; St = 0.198; A/d = 0.5
z/d = 22.9 z/d = 7.9
25
NEKTAR-ALE Simulations
Objectives:
Confirm numerically the existence of a stable,
periodic hybrid shedding mode 2S~2P in the
wake of a straight, rigid, oscillating cylinder
Principal Investigator:
Prof. George Em Karniadakis, Division of Applied
Mathematics, Brown University
Approach:
DNS - Similar conditions as the MIT experiment
(Triantafyllou et al.)
Harmonically forced oscillating straight rigid
cylinder in linear shear inflow
Average Reynolds number is 400
Vortex Dislocations, Vortex Splits & Force
Distribution in Flows past Bluff Bodies
D. Lucor & G. E. Karniadakis
Results:
Existence and periodicity of hybrid mode
confirmed by near wake visualizations and spectral
analysis of flow velocity in the cylinder wake and
of hydrodynamic forces
Methodology:
Parallel simulations using spectral/hp methods
implemented in the incompressible Navier- Stokes
solver NEKTAR
VORTEX SPLIT
Techet, Hover and Triantafyllou (JFM 1998)
VIV Suppression
Helical strake
Shroud
Axial slats
Streamlined fairing
Splitter plate
Ribboned cable
Pivoted guiding vane
Spoiler plates
26
VIV Suppression by Helical Strakes
Helical strakes are a
common VIV suppresion
device.
Galloping
Galloping is a result of a wake instability.
m
y(t), y(t)
.
Y(t)
U
-y(t)
.
V
Resultant velocity is a combination of the
heave velocity and horizontal inflow.
If
n
<< 2 f
v
then the wake is quasi-static.
27
Lift Force, Y()
Y(t)
V
C
y
=
Y(t)
1
/
2
U
2
A
p
C
y
Stable
Unstable
Galloping motion
m
z(t), z(t)
.
L(t)
U
-z(t)
.
V
b
k
a
mz + bz + kz = L(t)
.. .
L(t) =
1
/
2
U
2
a C
lv
- m
a
y(t)
..
C
l
() = C
l
(0) +
C
l
(0)
+ ...
Assuming small angles, :
~ tan = -
z
U
.
=
C
l
(0)
V ~ U
28
Instability Criterion
(m+m
a
)z + (b +
1
/
2
U
2
a )z + kz = 0
.. .
U
~
b +
1
/
2
U
2
a
U
< 0 If
Then the motion is unstable!
This is the criterion for galloping.
is shape dependent
U
1
1
1
2
1
2
1
4
Shape
C
l
(0)
-2.7
0
-3.0
-10
-0.66
29
b
1
/
2
a
( )
Instability:
=
C
l
(0)
<
b
1
/
2
U a
Critical speed for galloping:
U >
C
l
(0)
Torsional Galloping
Both torsional and lateral galloping are possible.
FLUTTER occurs when the frequency of the torsional
and lateral vibrations are very close.
30
Galloping vs. VIV
Galloping is low frequency
Galloping is NOT self-limiting
Once U > U
critical
then the instability occurs
irregardless of frequencies.
References
Blevins, (1990) Flow Induced Vibrations,
Krieger Publishing Co., Florida.