First published online May 2, 2008
Journal of Experimental Biology 211, 1541-1558 (2008)
Published by The Company of Biologists 2008
doi: 10.1242/jeb.015644
Numerical investigation of the hydrodynamics of carangiform swimming in the transitional and inertial flow regimes
Iman Borazjani and
Fotis Sotiropoulos*
St Anthony Falls Laboratory, University of Minnesota, Minneapolis, MN
55402, USA

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 1. (A,B) Different views of the virtual carangiform swimmer closely modeled
after a mackerel and meshed with triangular elements as needed for the
sharp-interface immersed boundary numerical method.
|
|

View larger version (9K):
[in this window]
[in a new window]
|
Fig. 2. Midlines of the fish for different instants during one tail beat cycle (A)
and the amplitude envelope profile (B).
|
|

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 3. Time history of the force coefficient normalized by the rigid body drag for
different St at Re=4000. Positive and negative values
indicate that the net force is of thrust- and drag-type, respectively.
|
|

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 4. Effect of Re and St on the mean force coefficient
produced by the tethered fish. The force coefficient is time-averaged and
normalized by the rigid body drag coefficient. The lower broken line shows the
rigid body drag coefficient and the upper broken line shows the zero mean
force coefficient, i.e. self-propulsion limit.
|
|

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 5. Variation of the skin (friction) drag, form drag and total drag with
Strouhal number at two Reynolds numbers: Re=300 (A) and
Re=4000 (B). The drag forces are calculated using
Eqn 17.
|
|

View larger version (70K):
[in this window]
[in a new window]
|
Fig. 6. Pressure contours and streamlines on the midplane of the fish relative to a
frame moving with the body wave phase speed V (Re=300). (A)
Rigid body (St=0). (B) Flow separates for St=0.1,
U/V=2.11. (C) Flow does not separate for St=0.3,
U/V=0.7.
|
|

View larger version (121K):
[in this window]
[in a new window]
|
Fig. 7. Calculated out-of-plane vorticity contours with velocity vectors for the
Re= , St=0.26 case (A) on the horizontal
(x1–x3) mid-plane and (B) the
vertical (x2–x3) plane.
|
|

View larger version (66K):
[in this window]
[in a new window]
|
Fig. 8. Instantaneous streamlines with vorticity contours showing (A) a single row
regular Karman street (Re=4000, St=0.2); (B) singe row
reverse Karman street (Re= , St=0.26); and (C) double
row reverse Karman street (Re=4000, St=0.7). The red arrows
show the general direction of the wake flow.
|
|

View larger version (74K):
[in this window]
[in a new window]
|
Fig. 9. Instantaneous streamlines and vorticity contours at the horizontal
mid-plane for St=0.3 highlighting the effect of Reynolds number on
the wake structure. (A) Re=300; (B) Re=4000 (middle); (C)
Re= .
|
|

View larger version (33K):
[in this window]
[in a new window]
|
Fig. 10. Three-dimensional (3D) vortical structures visualized using the
q-criterion showing 3D wake structures simulated for the
Re=300 case. (A) Double row wake at St=1.2; (B) single row
wake at St=0.3.
|
|

View larger version (42K):
[in this window]
[in a new window]
|
Fig. 11. Three-dimensional (3D) vortical structures visualized using the
q-criterion showing 3D wake structures simulated for the
Re=4000 case. (A) Double row wake at St=0.7; (B) single row
wake at St=0.2.
|
|

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 12. Three-dimensional (3D) vortical structures visualized using the
q-criterion showing 3D wake structures simulated for the inviscid
case. (A) Double row wake at St=0.7; (B) single row wake at
St=0.26.
|
|

View larger version (58K):
[in this window]
[in a new window]
|
Fig. 13. The iso-surfaces of q cut with (A)
x1–x3 and (B)
x2–x3 mid-planes showing the
footprints of the wake structure using the out-of-plane vorticity contours.
The inset (C) gives a closer look at the structures cut with both
x1–x3 and
x2–x3 mid-planes
(Re=4000, St=0.7).
|
|

View larger version (28K):
[in this window]
[in a new window]
|
Fig. A1. Time history of the inline force coefficient (solid lines) and its pressure
(dotted lines) and viscous (broken lines) components from the present
computations (red lines) compared with the computations of Dutsch et al.
(Dutsch et al., 1998 ) (black
lines).
|
|

View larger version (21K):
[in this window]
[in a new window]
|
Fig. A6. Time history of force coefficient normalized by the rigid body drag
coefficient for three different grid sizes at Re=4000 and
St=0.5. Positive and negative values indicate that the net force is
of thrust- and drag-type, respectively.
|
|

CiteULike
Complore
Connotea
Del.icio.us
Digg
Reddit
Technorati
Twitter What's this?
© The Company of Biologists Ltd 2008