QUICK SEARCH:   [advanced] Author: Keyword(s):
Home     Help     Feedback     Subscriptions     Archive     Search     Table of Contents

First published online June 13, 2008
Journal of Experimental Biology 211, 2087-2100 (2008)
Published by The Company of Biologists 2008
doi: 10.1242/jeb.016279

This Article Summary Figures Only Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Related articles in JEB Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Zhu, Q. Articles by Shoele, K. Search for Related Content PubMed PubMed Citation Articles by Zhu, Q. Articles by Shoele, K. Social Bookmarking                         What's this?

Propulsion performance of a skeleton-strengthened fin

Qiang Zhu^* and Kourosh Shoele

Department of Structural Engineering, University of California, San Diego, La Jolla, CA 92093, USA

^* Author for correspondence (e-mail: qizhu{at}ucsd.edu)

Accepted 8 April 2008

	Summary

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
We examine numerically the performance of a thin foil reinforced by embedded rays resembling the caudal fins of many fishes. In our study, the supporting rays are depicted as nonlinear Euler–Bernoulli beams with three-dimensional deformability. This structural model is then incorporated into a boundary-element hydrodynamic model to achieve coupled fluid–structure interaction simulation. Kinematically, we incorporate both a homocercal mode with dorso-ventral symmetry and a heterocercal mode with dorso-ventral asymmetry. Using the homocercal mode, our results demonstrate that the anisotropic deformability of the ray-reinforced fin significantly increases its capacity of force generation. This performance enhancement manifests as increased propulsion efficiency, reduced transverse force and reduced sensitivity to kinematic parameters. Further reduction in transverse force is observed by using the heterocercal mode. In the heterocercal model, the fin also generates a small lifting force, which may be important in vertical maneuvers. Via three-dimensional flow visualization, a chain of vortex rings is observed in the wake. Detailed features of the wake, e.g. the orientation of the vortex rings in the heterocercal mode, agree with predictions based upon particle image velocimetry (PIV) measurements of flow around live fish.

Key words: fish locomotion, flexible fin, skeleton-reinforced membrane, fluid–structure interaction

	INTRODUCTION

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
Skeleton-reinforced bio-membranes are ubiquitous and play critical roles in many biological functions. For example, the flexibility of the wings of many insects is determined primarily by the architecture of the embedded veins. By experimentally measuring and numerically characterizing the wing structure of 16 different species of insects, Combes and Daniel concluded that an anisotropic flexibility is achieved through the architecture of these composite structures (Combes and Daniel, 2003a; Combes and Daniel, 2003b). In all the species they studied, the spanwise bending stiffness was found to be at least 1–2 orders of magnitude higher than the chordwise bending stiffness (see also Newman and Wootton, 1986; Wootton, 1992). It was attributed mostly to the existence of clustered or thickened veins near the leading edge. This suggests that, through specific arrangement of supporting veins, these insect wings allow certain passive deformations while discouraging others. Structures with similar architecture are also found in the fins of aquatic animals. Many fishes possess caudal, dorsal, pectoral or ventral fins that are stiffened by fin rays embedded in a thin layer of collagenous membrane (Lauder, 1989; Lauder and Drucker, 2004; Fish and Lauder, 2006; Lauder and Madden, 2007). Overall, there are over 25 000 species of fish with these ray-strengthened fins. According to morphological studies (Harder, 1971; Kardong, 1998), a fin ray contains a central bundle of collagen surrounded by small segmented bony elements called hemitrichs. These elements are paired and resemble a bimetallic strip with two elongated bony elements separated by the central collagen core. Along the length of a hemitrich, and at the ends, there exist short ligaments and elastic fibers. In addition, the basal end of each ray attaches to four separate muscles. This architecture enables a fish to control the motion of each individual ray and thus achieve three-dimensional fin deformations. This property is expected to have fundamental effects on the hydrodynamic performance of the fins in generating propulsion force during cruising, bursting and maneuvering.

The exact function of the anisotropic flexibility of skeleton-strengthened membranes in insect flying and fish swimming remains to be fully understood. Despite a limited number of investigations (Cheng et al., 1998; Pedley and Hill, 1999), existing studies in fish swimming and insect flying are predominantly based upon the assumption that the wings and fins are rigid (Liu and Kawachi, 1998; Ramamurti and Sandberg, 2002; Ramamurti and Sandberg, 2007; Wang, 2000; Pullin and Wang, 2004) or with prescribed deformations (Liu and Bose, 1997; Ramamurti et al., 2005; Mittal et al., 2006). It has always been speculated that the structural flexibility of these composite bio-structures, with their stiffness determined by the coupling of the embedded skeletons and the surrounding materials, attributes to the highly efficient force generation through flapping motions. Evidence of the beneficial effects of structural flexibility in the generation of lift and propulsion forces comes from different sources, in particular the analytical, numerical and experimental investigations of flapping foils, as discussed below.

Unlike conventional aero- and hydro-foils, which utilize steady fluid dynamics for force generation, flapping foils are modeled on fish fins and insect wings so that they produce aerodynamic or hydrodynamic lift, propulsion or maneuvering forces using unsteady foil oscillations (Triantafyllou et al., 1991; Tuncer and Platzer, 1996). Through various investigations, it has been shown that structural flexibility, per se, increases the propulsion efficiency in many cases. For example, by experimentally studying the performance of a semi-flexible flapping foil, Yamamoto et al. reported up to a 27% increase in propulsion efficiency compared with a rigid foil (Yamamoto et al., 1995). Similar performance enhancement was observed by Heathcote et al. in their investigation of a thin steel plate undergoing periodic heaving motion in still water (Heathcote et al., 2004). By theoretically examining the oscillation of a two-dimensional flexible plate, Katz and Weihs were able to obtain a 20% increase in propulsion efficiency (Katz and Weihs, 1978). The performance of a flexible foil, by contrast, has been found to depend not only on its elasticity but also on its inertia (Pederzani and Haj-Hariri, 2006). This conclusion was confirmed in our recent work (Zhu, 2007). In that study, we investigated the three-dimensional fluid–structure interaction of a flapping foil. A number of extreme, yet representative, cases were examined, including chordwise vs spanwise flexibility and fully-coupled (the foil deformation is driven mostly by the fluid) vs partially decoupled (the foil deformation is driven mostly by its own inertia) interactions. By using a boundary-element model coupled with a two-dimensional plate model, we were able to show that both the structural flexibility and the inertia of a flapping foil had significant effects on its capability for thrust generation. Among all the numerical tests we have conducted, two cases – one with chordwise flexibility and fluid-driven deformation and the other with spanwise flexibility and inertia-driven deformation – yield higher propulsion performance than rigid foils. The former case shows increased efficiency (by 20%, consistent with previous studies) and the latter shows increased thrust (by almost twofold in some cases) without any significant loss in efficiency.

Structural flexibility has also been found to effectively reduce fluid drags. Studying the two-dimensional soap film flow around an elastic glass fiber as an example, Alben et al. characterized the drag reduction effect of structural flexibility, illustrating the transition from a rigid-body drag scaling law of U² to U^4/3 (U is the speed of incoming flow) when the fiber flexibility exceeds a threshold (Alben et al., 2002; Alben et al., 2004).

Despite the abovementioned efforts to characterize dynamics of flexible foils or fibers, the hydrodynamic significance of three-dimensional anisotropic flexibility is still not illustrated. In particular, comprehensive studies are required to understand the coupled fluid–structure interactions of bio-membranes with morphological relevance, e.g. those with mechanical properties similar to skeleton-reinforced fins or wings.

We herein present a numerical investigation by coupling fluid dynamics and structural mechanics to examine the propulsion performance of a skeleton-reinforced fin that is geometrically and structurally similar to the caudal fin of a fish. We study cases in which the unsteady motions of the rays are individually controlled by imposing sway–yaw motions at their basal ends. The rest of the ray may deform under the combined effect of its own inertia, the hydrodynamic load and constraint from the membrane. Our approach includes an Euler–Bernoulli beam model for the strengthening rays, as well as a boundary-integral equation representation for the surrounding fluid. To account for large three-dimensional deformations, the beam model is fully nonlinear and allows stretching, bending and twisting. Via systematic simulations, dynamics of a flexible fin and that of a rigid one are compared with each other using the performance metrics, including the capacity and efficiency of the fin to generate thrust and the reduction of transverse force. In addition, we investigate both the homocercal mode, in which the fin deformation is symmetric across the span, and the heterocercal mode, in which the deformation of the dorsal side of the fin is different from that of the ventral side.

The rest of the paper will be organized as follows. In the next section we define the geometry and internal structure of the fin, together with the mathematical formulations and the numerical methods to solve the coupled fluid–structure interaction problem. Following that, numerical results, including the thrust force, the transverse force and the propulsion efficiency of the flexible fin in comparison with a rigid fin, will be presented. Also included will be comparisons between the homocercal and the heterocercal modes. Finally, conclusions will be drawn.

	MATERIALS AND METHODS

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
Fin structure and kinematics
As a model of the real ray-strengthened fin (Fig. 1A), we consider a geometrically and structurally simplified fin as shown in Fig. 1C. This trapezoid fin has chord length c and leading edge span. In a Cartesian coordinate system (x, y, z), the center plane of the fin at its mean position locates within the x–z plane and is symmetric with respect to the x–y plane. This fin is reinforced by nine embedded rays, numbered 1 to 9. The basal ends of the rays are distributed evenly on the leading edge, as are the angles they form with respect to the x-axis, ranging from 30° (Ray 1) to –30° (Ray 9). The thickness of the fin, d, is assumed to be uniform except for tapering at the leading and the trailing edges (note that the value of d is only relevant in the hydrodynamic problem). To put it in perspective, this configuration is similar to the homocercal caudal fin of species such as bluegill sunfish (Lepomis macrochirus).

View larger version (22K): [in this window] [in a new window]   Fig. 1. (A) Schematic showing a typical ray-strengthened caudal fin [modified from Alben et al. (Alben et al., 2007)]. We note that the rays are segmented and are tapered and branched near the trailing edge. (B) The dorsal view of the internal structure of a ray [courtesy of S. Alben and the Royal Society (Alben et al., 2007)]. (C) A geometrically and structurally simplified fin employed in the present modeling study. See List of symbols and abbreviations for definitions.

We consider a three-dimensional flapping motion of the fin in combination with a constant forward speed, U, in the –x direction. The flapping motion is actuated by imposing sway–yaw motions on each ray. By assumption, the rays share the same sway motion y(t)=y₀sin(t) following the motion of the caudal peduncle. The yaw motion with respect to the z-axis of Ray i is (t)=_isin(t–_i), i=1,..., 9. The Strouhal number, a non-dimensional parameter characterizing the unsteady fluid-dynamic effects, is defined to be S_t=2y₀f/U, where f=1/T=/2o is the frequency of the motion and T is the period.

In the following, we elaborate the mathematical formulations of the ray motion and the fluid dynamics.

Mathematical formulations and numerical algorithms
We apply a three-dimensional nonlinear formulation and model each ray as an Euler–Bernoulli beam that could be stretched, bent and twisted (Tjavaras et al., 1998). Based upon the Euler–Lagrangian dual-coordinate approach and the robust Euler parameters, this method achieves fully-nonlinear simulation of arbitrary ray deformations. A finite difference algorithm, the modified box method, is then employed to solve these nonlinear equations (see Appendix A for details).

The flow around the fin is assumed to be irrotational except for zero-thickness sheets of vorticity shed from the trailing edge so that it can be described by a flow potential (x,t), where x=(x, y, z) is the position vector. The near-body flow field as well as the hydrodynamic loading on the fin are then predicted using a boundary-element approach based upon the boundary-integral framework by segmenting the fin surface S_b into N_b panels S_bj, j=1,..., N_b (Appendix B).

By solving the potential flow around the body and determining the flow potential , the hydrodynamic pressure p on the fin surface S_b is readily obtained through Bernoulli's equation. We have:

(1)

where is the density of water (10³ kg m^–3). Integrating p over the fin surface, we obtain the hydrodynamic force as:

(2)

The power required to drive the fin is given as:

(3)

By definition, the thrust force, F_t, is the –x component of F. The propulsion efficiency, , is defined as:

(4)

where _t and represent the thrust force and power consumption averaged over one period T, respectively. The fin oscillation also creates a transverse force F_r, which is the y component of F. In the heterocercal mode, there also exists a non-zero lifting force, F_z, leading to a mean lifting coefficient _z/U²c².

We thus characterize the performance of the fin by its mean thrust coefficient, _t/U²c², the propulsion efficiency, , the mean lifting coefficient, and the transverse force coefficient, F_r⁽¹⁾/U²c², where:

(5)

As mentioned before, the structural function of the membrane in which the rays are embedded is simplified as linear spring-dampers spanning between neighboring rays. Specifically, in the following calculations, we chose the spring constant to be 20 N m^–2, and the damping constant to be 2 Ns m^–2 (both are measured per unit length). With such a constraint, in the following simulations the maximum variation of the surface area of the fin is around 5–10%.

We employ an iteration method to solve the fully coupled fluid–structure interaction problem. Below we list the primary steps of this algorithm.

1. At time t, we start from an initial guess, , of the configuration of the fin.
   
2. The hydrodynamic problem is solved based on .
   
3. After the hydrodynamic pressure, p, within panel j is evaluated using Eqn 1, we assume that the hydrodynamic force acting on this panel, pnS_bi (S_bi is the area of the panel), is transferred equally to four grid points that are closest to the panel on the two neighboring rays (two points on each ray). Summing up the forces on all the panels, the hydrodynamic loads, F_h, on the nine rays are thus determined.
   
4. The Bernoulli–Euler beam equations (EqnA9) are numerically integrated to obtain the updated positions of the rays. Based upon these, the updated configuration of the fin, S_b, is determined via linear interpolation.
   
5. If S is not sufficiently close to , we set and repeat Steps 2–4.
   

To eliminate the well-known added mass induced instability, we apply an implicit added mass scheme. In this approach, the added mass effects are subtracted from both sides of Eqn A9 during time integration. As long as convergence of the iteration is achieved, the exact values of the added mass will not affect the results. A detailed description of this method is provided in Connell and Yue (Connell and Yue, 2007).

Convergence tests
The validity and accuracy of both the Euler–Bernoulli beam solver and the boundary integral solver have been extensively tested and documented through numerous convergence studies and comparisons with experiments or other numerical/theoretical predictions (Tjavaras et al., 1998; Zhu et al., 2002; Zhu et al., 2006; Zhu, 2007). These are omitted here for brevity. Instead, we will concentrate primarily on the convergence of the algorithm with respect to the time step, as it reflects the accuracy of our iteration scheme discussed above, which is the core of the interaction solver.

Convergence tests are performed by examining the propulsive force generated by a fin with leading-edge span s₀=0.02 m, chord length c=0.1 m, and thickness d=0.002 m. The rays are structurally identical except for their lengths. Each has a cross-sectional diameter of 1 mm, a Young's modulus, E, of 0.4 GPa, and a hysteretic (material) damping coefficient, D, of 0.4 sGPa. The density of the rays is chosen to be the same as that of the water, i.e. 10³ kg m^–3.

Kinematically, we set the forward speed to be U=0.1 ms^–1. The sway amplitude, y₀=0.5c. The Strouhal number, S_t=0.2. We also assume that the nine rays shear the same yaw motion, with _0i=₀=10⁰ and _i=90° (i=1,..., 9).

Table 1 shows the mean thrust coefficient, _t/U²c², as a function of the number of body elements, N_b, and the time step t. Linear convergence is achieved with respect to both parameters.

View this table: [in this window] [in a new window]   Table 1. Convergence of the mean thrust coefficient with respect to the number of body panels (Nb) and time step (t)

	RESULTS

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
Lauder (Lauder, 2000) studied kinematics of caudal fins of three species: the leopard shark (Triakis semifasciata), the sturgeon (Acipenser transmontanus) and the bluegill sunfish (Lepomis macrochirus). Two of these fishes, the shark and the sturgeon, possess caudal fins with dorso-ventral asymmetry (heterocercal fins). The bluegill sunfish, on the other hand, has caudal fins that are geometrically symmetric across the span (homocercal fins). Surprisingly, his observations illustrated that, despite the geometric difference between homocercal and heterocercal fins, their motions are nevertheless heterocercal in many cases. Specifically, the dorsal side of the fin was found to be ahead of the ventral side. In such a heterocercal mode, a small vertical force is generated, which has to be balanced by forces generated by the anterior part of the fish body. The heterocercal motion was attributed to the embedded muscle through which the kinematics of the rays can be controlled. Similar heterocercal modes in homocercal fins have been found on chub mackerel (Scomber japonicus) (Gibb et al., 1999). On the other hand, in another experiment of a bluegill sunfish, it was observed that the caudal fin performed mostly homocercal motions (Tytell, 2006). The implication is that a fish with a homocercal caudal fin may be able to switch between different swimming modes.

By specifying the motion of each individual reinforcing ray, in our model we are able to achieve controlled anisotropic deformations, both homocercal and heterocercal, of the fin. We assume that the basal ends of the rays share the same sway motion following the undulation of the anterior part of the fish body. The yaw motions, i.e. the rotation of the rays around their basal ends, can be varied individually. This imitates the muscle action that controls the ray kinematics (Lauder and Drucker, 2004). In the following, we will present our numerical results of the force generation by a rigid fin and a flexible fin with identical geometries. The flexible fin has the same geometric and structural properties as the one applied for the convergence tests. In addition, we assume that in all the simulations the forward speed, U, is 0.1 m s^–1, the sway amplitude, y₀, is 0.05 m, and the phase angle between sway and yaw, _i, is 90° (i=1,..., 9). The focus of our numerical simulation will be the comparison between force generation by a rigid fin and that by a flexible fin. We will also study the differences and similarities between homocercal and heterocercal modes.

Numerically, 1000 boundary elements are applied on the fin surface. On each ray, the number of grids is N_r=32. Finally, a time step of t=T/128 is chosen.

Homocercal modes
In the homocercal motion, we assume that the yaw amplitudes of all the rays are the same, i.e. _0i=₀ (i=1,..., 9). A typical sequence of the deformation of a flexible fin over period is displayed in Fig. 2. At t=0, when the fin is at its mean sway position (y=0) and maximum yaw angle (=10°), we observe that there is no pronounced spanwise (dorso-ventral) bending. As t increases towards T/8 and T/4, accompanying the increasing sway, significant dorso-ventral bending is observed, creating a bow-shaped trailing edge so that the central part of the fin bends into the flow. If we plot the trajectories of the trailing ends of Ray 5 (at the center) and Ray 1 (at the upper edge), it is clear that the upper and lower edges of the fin undergo smaller and lagged-behind undulations compared with the center part (Fig. 3). In Fig. 2, we also observe that the dorso-ventral bending reverses its direction near the upper and lower edges. There are consequently two inflection points.

View larger version (15K): [in this window] [in a new window]   Fig. 2. Deformation of the fin in the homocercal mode within period. The arrow illustrates the direction of sway. The forward motion is subtracted. St=0.2, 0=10°. t, time; T, period.

View larger version (5K): [in this window] [in a new window]   Fig. 3. The y-displacements (normalized by c) of the trailing ends of Ray 1 and Ray 5 during a period in the homocercal mode. St0.2, 0=10°. t, time; T, period.

View larger version (25K): [in this window] [in a new window]   Fig. 4. Distribution of the pressure coefficient Cp=p/U2 on (A) the left side and (B) the right side of the surface of a flexible fin in the homocercal mode at t=T/8. St=0.2, 0=10°. t, time; T, period.

View larger version (14K): [in this window] [in a new window]   Fig. 5. The thrust coefficient of (A) a rigid fin and (B) a flexible fin as a function of the Strouhal number, St, and the yaw amplitude, 0, in the homocercal mode.

The bow-shaped trailing edge of the fin, as well as the existence of inflection points near the upper and lower edges, is qualitatively consistent with a recent report about kinematics of the caudal fin of a bluegill sunfish (Tytell, 2006) (see fig. 3 in that paper).

View larger version (17K): [in this window] [in a new window]   Fig. 6. The propulsion efficiency of (A) a rigid fin and (B) a flexible fin as a function of the Strouhal number, St, and the yaw amplitude, 0, in the homocercal mode.

View larger version (9K): [in this window] [in a new window]   Fig. 7. The transverse force coefficient of (A) a rigid fin and (B) a flexible fin as a function of the Strouhal number, St, and the yaw amplitude, 0, in the homocercal mode.

View larger version (11K): [in this window] [in a new window]   Fig. 8. (A) The thrust coefficient, (B) the transverse force coefficient and (C) the propulsion efficiency of rigid and flexible fins in the homocercal mode at St=0.3.

Both the increase in propulsion efficiency and the decrease in transverse force indicate performance enhancement through three-dimensional structural flexibility. These characteristics are more clearly illustrated in Fig. 8, where cross-sections of Figs 5, 6, 7 are plotted.

Heterocercal modes
In our simulations, the heterocercal mode, i.e. a mode with dorso-ventral asymmetry, is achieved by varying the amplitudes (but not the frequency and phase) of yaw motions of the rays. Specifically, we assume that the ray at the upper edge on the dorsal side (Ray 1) performs no yaw motion (₀₁=0°). For the rest of the rays, the yaw amplitude increases linearly with the largest yaw motion occurring at the lower edge (Ray 9), i.e. ₀₁=i/9x₀. Hereafter, we take a closer look at a flexible fin.

As shown in Fig. 9, the ray at the upper edge (Ray 1) leads both the central ray (Ray 5) and the ray at the lower edge (Ray 9). This trend is more clearly shown in Fig. 10, where the trajectories of the trailing ends of these three rays are plotted. Such behavior is consistent with experimental observations of bluegill sunfish as reported by Lauder (Lauder, 2000).

View larger version (12K): [in this window] [in a new window]   Fig. 9. Deformation of the fin in the heterocercal mode within period. The arrow demonstrates the direction of sway. The forward motion is subtracted. St=0.2, 0=35°. t, time; T, period.

View larger version (5K): [in this window] [in a new window]   Fig. 10. The y-displacements (normalized by c) of the trailing ends of Ray 1, Ray 5, and Ray 9 during a period in the heterocercal mode. St=0.2, 0=35°. t, time; T, period.

View larger version (28K): [in this window] [in a new window]   Fig. 11. (A) Thrust coefficient, (B) propulsion efficiency, (C) transverse force coefficient and (D) lifting coefficient of a flexible fin in the heterocercal mode. St, Strouhal number; 0, yaw amplitude.

Fig. 11D shows the existence of a mean lifting force as generated by the heterocercal motion of the fin. Although this force is 1–2 orders of magnitude smaller than the thrust, it may be an important source of maneuvering force in the vertical direction.

Visualization of three-dimensional flow field
To further investigate the fluid–structure interaction problem involved in the sway/yaw motion of the fin, and to compare the homocercal mode with the heterocercal mode, we numerically visualized both the near-body flow field and the vorticity distribution in the wake.

As illustrated in previous studies, an important phenomenon concerning flexible objects in flow is that these objects often bend their surfaces to match the incoming stream. This has been confirmed in various experiments, including, for example, the bending of a glass fiber by a soap film flow (Alben et al., 2004) and the meandering motion of a dead fish between incoming vortices (Liao et al., 2003). To confirm if this is the case in our problem, in Fig. 12 we checked the inplane streamlines within two horizontal planes, z=0 and z=0.3c, around a rigid fin and a flexible fin undergoing sway motions. It was found that at the mid-span (z=0), the fin is bent just a little towards the incoming flow. This results from the fact that the central ray (Ray 5) is the shortest and least deformable. On the plane z=0.3c, however, two different tendencies are observed. The front part of the fin tends to align with the incoming flow, echoing features of the two-dimensional fluid-driven foil deformations (Zhu, 2007). The hind part of the fin, however, undergoes much smaller reorientation. This is caused by the three-dimensional effect, in particular interactions with rays closer to the center via the membrane.

View larger version (38K): [in this window] [in a new window]   Fig. 12. Inplane streamlines on z=0 (plane a) and z=0.3c (plane b) around (A) a rigid fin and (B) a flexible fin in the homocercal mode. St=0.3, 0=0°. The streamlines are plotted in a reference system moving with the forward and the sway motions. The locations of plane a and plane b are illustrated in the inset.

Fig. 13A,B shows the iso-surfaces of vorticity in the wake behind a rigid fin and a flexible fin undergoing a homocercal flapping motion, respectively. It is evident (especially from the insets of Fig. 13A,B) that the wake consists of a sequence of vortex rings symmetric in the dorsal/ventral direction (in two-dimensional experimental visualizations or numerical simulations, these vortex rings appear as pairs of vortices). The wakes of the rigid fin and the flexible fin resemble each other except for a subtle difference at the connection between neighboring rings. Specifically, in the flexible case, the figure shows elongated forklike structures between neighboring rings (Fig. 13B). These structures resemble the hairpin structures proposed by Tytell (Tytell, 2006) based upon analysis of PIV data obtained from the near-body flow around a bluegill sunfish. The hairpins were attributed to the cupping of the fin during flapping, as is also shown in our simulations (see Fig. 2).

View larger version (41K): [in this window] [in a new window]   Fig. 13. Iso-surfaces of vorticity in the wakes behind (A) a rigid fin in the homocercal mode, (B) a flexible fin in the homocercal mode and (C) a flexible fin in the heterocercal mode. The iso-surfaces are plotted at the level of 0.9 (after normalization by U/c). St=0.3. For the homocercal mode, 0=10°; for the heterocercal mode, 0=35°. U, forward speed; c, chord length of the foil.

View larger version (90K): [in this window] [in a new window]   Fig. 14. Flow field within the y=0 plane behind a flexible fin in (A) the homocercal mode and (B) the heterocercal mode. xl is the x-location of the leading edge. The contour displays the y-component of the vorticity (normalized by U/c). St=0.3. For the homocercal mode, 0=10°; for the heterocercal mode, 0=35°.

	DISCUSSION AND CONCLUSIONS

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
By using a fluid–structure interaction model combining a nonlinear Euler–Bernoulli beam model and a boundary-integral hydrodynamics model, we have investigated the force generation by a flexible fin made of a thin membrane reinforced by embedded rays. This structure resembles the ray-strengthened fins of many fishes. Structurally, the distribution of reinforcing rays creates a anisotropic flexibility. Two different locomotion modes, the homocercal (dorso-ventral symmetric) motion and the heterocercal (dorso-ventral asymmetric) motion, are studied.

Our results show that compared to rigid fins with the same geometry, a fin with anisotropic flexibility yields higher propulsion efficiency. The structural flexibility also reduces the transverse force, another beneficial effect in fish swimming. Further reduction of transverse force is observed by applying the heterocercal mode.

An interesting phenomenon is that in terms of propulsion efficiency, the fluid–structure coupling greatly reduces the sensitivity of the fin performance to kinematic parameters such as the Strouhal number and the amplitude of yaw motion. Even in cases when these parameters are not optimized, the flexible fin is able to function well as a propeller. This property may have significant effects on fish locomotion. The biological implication is that it may not be critical for a fish to control the motion of its caudal fin accurately to achieve high propulsion performance. In biomimetic applications, this suggests that in the future development of mechanical devices imitating the swimming mechanism of fishes, the controlling system could be significantly simplified provided that the structural properties of the propeller are designed to mimic the ray-reinforced fins.

In addition to the locomotion modes examined in our work, there may exist other combinations of fin structures and kinematics, leading to different fin deformations in the swimming process. For example, by observing the fin actuation of a dace (Leuciscus leuciscus), Bainbridge concluded that the maximum dorso-ventral curvature occurred when the caudal fin is close to its maximum sway displacements (Bainbridge, 1963). His observations also show that the upper and lower edges of the fin usually precede the center part in their lateral motions. This contradicts what we found in the homocercal mode when the upper/lower edges are behind the central part. Similarly, although the fin kinematics predicted by our model appear to resemble the experimental observations of Tytell (Tytell, 2006), subtle differences, such as the relationship between the center part and the upper/lower edges in terms of amplitudes of deflection and phases, suggest more complicated fin kinematics in his experiment than in our simulations. In species with heterocercal tails, such as spiny dogfish (Squalus acanthias), a sophisticated ring-within-a-ring vortex structure has been observed (Wilga and Lauder, 2004). This was attributed to the asymmetry in vortex formation at the upper and the lower lobes. A more sophisticated ray model, with actively controlled curvatures, is also required. To study cases with closer biological relevance, and to fully understand the benefits of different fin architectures and locomotion modes, interdisciplinary investigations including fish morphology and physiology, fish behavior, control, as well as fluid–structure interactions are required.

The perspective gained from our work and previous studies makes it appropriate to envisage future development of biomimetic propellers imitating flexible flapping fins made of skeleton-reinforced membranes. Compared with the conventional rigid foil design, these novel biomimetic propellers possess the following advantages.

1. Enhanced performance: as illustrated in our study, a flexible foil with three-dimensional deformability yields better performance than a rigid one, manifested in the increased efficiency, reduced sensitivity to kinematic parameters and diminished force in the transverse direction.
   
2. Controllability: with strengthening skeletons, the structural properties of the composite membrane can be adjusted in different operation modes. For example, our results indicate that a more rigid fin may generate larger thrusts in a certain parameter range. Therefore, in burst swimming, it may be desirable to increase the stiffness of the supporting rays.
   
3. Structural strength and lightness: the biomimetic composite membrane provides a light structure with high strength.
   
4. Deployability: the skeleton-reinforced membrane resembles deployable structures such as cable roof buildings. These structures can be easily folded and unfolded, making them an ideal design as a portable device.
   

It is necessary to point out the limitations of the current model. First, our work concerns a geometrically and structurally simplified caudal fin. Effects of the fish body and attachments such as dorsal, pectoral or anal fins are not considered. As illustrated by Wilga and Lauder (Wilga and Lauder, 2002), in leopard sharks (Triakis semifasciata) and bamboo sharks (Chiloscyllium punctatum), the inclination angle of the body plays a pivotal role in balancing the torque created by the tail. In addition, as shown by Tytell (Tytell, 2006), vortices shed from the dorsal and the anal fins of a bluegill sunfish may contribute significantly to the near-body flow. Recent experiments with a mechanical fish have shown that the fish body itself might also generate a sequence of weaker vortex rings in the wake (Brucker and Bleckmann, 2007). The dynamic effects of the vortices generated upstream of the caudal fin depend not only on their strength but also on their phase with respect to the caudal fin (Zhu et al., 2002). This, in turn, is determined by the exact morphology of the fish as well as its kinematics. Another effect not considered in the present study is vorticity generation from the fin surface in locations other than the trailing edge, in particular the upper and lower edges. This vorticity generation will not only change the local dynamics (e.g. the pressure distribution) around the upper/lower edges but will also affect the overall performance of the fin in propulsion. To fully account for this effect, a Navier–Stokes solver capable of studying three-dimensional fluid–structure interactions is required [for a possible method to achieve this, see Connell (Connell, 2006)]. A fully viscous study will not only illustrate effects of vortex generation from places other than the trailing edge but will also provide more accurate predictions of the flow field and the hydrodynamic forces, especially in cases when the Reynolds number is low. These Navier–Stokes simulations, however, are computationally expensive. In that respect, a potential flow method like ours provides an alternative way to conduct systematic tests over a wide range of parameters. To improve the accuracy of our current method, investigations are underway to develop methods to study vortex generation from the leading edge area or the upper/lower edges of the fin. The knowledge gained from these studies will pave the way for future simulations with more comprehensive numerical tools and experiments using mechanical devices. These studies may eventually contribute to the creation of biomimetic apparatus that will be installed on AUVs (autonomous underwater vehicles) or other vehicles (see Tangorra et al., 2007).

Finally, we would like to mention that, in addition to insect wings and fish fins, there exist a large number of bio-structures in nature with similar structural characteristics, which can thus be categorized as skeleton-reinforced membranes. A typical example is found in the membrane of the erythrocyte (red blood cell), which features a composite structure including a lipid bilayer and a protein skeleton consisting primarily of actin and spectrin (Mohandas and Evans, 1994). The combined visco-elasticity of the bilayer and the skeleton, as well as their dynamic inter-connectivity, impart remarkable stability and deformability essential for the cell functionality when it circulates around the body and squeezes through capillaries half its own diameter. Similar structures are found in the biopolymer membranes within the shells of mollusks such as the red abalone (Haliotis rufescens) and the nautilus (Nautilus pompilius). Remarkable mechanical properties of these structures (e.g. structural strength, stability, durability and deployability) suggest biomimetic applications and justify concentrated research efforts to understand the detailed structural mechanics and fluid–structure interactions of these bio-structures.

	APPENDIX A. MODEL OF RAYS VIA NONLINEAR BEAM DYNAMICS

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 
Our ray model is based upon a three-dimensional nonlinear Euler–Bernoulli beam formulation reported by Tjavaras et al. (Tjavaras et al., 1998). This paradigm employs an Euler–Lagrangian dual coordinate system, including the above-defined global coordinate system (x, y, z) and a local reference system with unit vectors , n and b in the tangential, normal and bi-normal directions of the ray, respectively. The two coordinate systems are related by:

(A1)

where i, j and k are unit vectors in the x, y and z directions, respectively. C is the orthogonal tensor representing the change of orientation from the global to the local reference frames. C is constructed by using the singularity-free method of Euler parameters based upon Euler's principal rotation theorem, which depicts an arbitrary reorientation as a single rotation by a principal angle about the principal unit vector of C, i.e. l. Correspondingly, four Euler parameters – β₀, β₁, β₂ and β₃ – are defined as:

(A2)

The rotation matrix, C, thus becomes:

(A3)

Dynamic equations of the ray movement are derived by considering the conservation of translational and angular momenta of an infinitesimal segment of the ray. We have:

(A4)

and

(A5)

where m is the mass per unit length, and s is the distance from this point to the basal end along the unstretched ray. V(s,t)=V+V_nn+V_bb and (s,t)=+_nn+_bb are the translational and angular velocities, respectively. T(s,t)=T+T_nn+T_bb is the internal force inside the ray. M(s,t)=M+M_nn+M_bb is the internal moment. (s,t) is the strain. (s,t)=+_nn+_bb is the Darboux vector measuring the material torsion and curvatures of the ray. F_c is the constraining force from the membrane that controls the distance between neighboring rays. F_h is the hydrodynamic force to be defined later.

The compatibility relations, which state that the ray's configuration must be continuous in both space and time, should also be satisfied as:

(A6)

The internal forces and moments are related to the strain and the Darboux vector through the constitutive relations T=EA, M=GJ, M_n=EI_n, M_b=EI_b. By definition, A is the cross-sectional area. EI and GJ are the bending and torsional stiffnesses, respectively. Our current simulations do not consider cases with ray twisting; therefore, the value of GJ is irrelevant. At this point, a hysteretic (material) damping can be incorporated by replacing the Young's modulus E with E+D/t, where D is the damping coefficient. In addition, the angular velocity is expressed in terms of the Euler parameters by:

(A7)

The governing equations are closed by applying the spatial derivative of the Euler parameters. We have:

(A8)

Finally, the system of equations is summarized into a vector form, i.e.:

(A9)

where Y=[T_nT_bVV_nV_bβ₀β₁β₂β_3nb]^T. Detailed description of the matrices H and P is provided by Tjavaras et al. (Tjavaras et al., 1998).

Overall, 13 boundary conditions are required at the two ends of each ray. At the basal end, we apply: (1) β₀=0; (2) the velocities V, V_n and V_b, as determined by the prescribed sway motion; (3) M=0; and (4) the three-dimensional orientation of the ray as determined by the prescribed yaw motion. In practice, condition 4 is enforced by adding a stiff rotational spring between the basal end of the ray and a reference direction representing the yaw motion. This leads to moments M_n and M_b at the basal end, whose values increase as the ray deviates from the reference direction. Six additional boundary conditions are applied at the trailing end, including: (1) ; (2) the internal force T=0; and (3) the internal moments M_n and M_b disappear so that _n=_b=0.

To solve Eqn A9 and determine the dynamics of a single ray, we numerically segment this ray into N_r–1, each with length s, by distributing N_r points s_k (k=1,..., N_r) along its unstretched length. A modified box method is then applied to integrate the system of equations from time step t_i–1 to t_i=t_i–1+t. We have:

(A10)

An iterative algorithm is employed to solve this system of implicit equations (Tjavaras et al., 1998). The method possesses second-order accuracy in space. The accuracy of time integration, on the other hand, is first order. The purpose of selecting a lower-order time integration scheme rather than the more conventional box method, which possesses second-order accuracy, is to providing better numerical stability.

APPENDIX B. THE BOUNDARY-INTEGRAL FORMULATION AND THE BOUNDARY-ELEMENT METHOD
The flow around the fin is assumed to be irrotational except for zero-thickness sheets of vorticity shed from the trailing edge. The flow field is then easily described by a velocity potential (x,t) (the position vector x(x, y, z) represents an arbitrary point in space). We further express as the linear superposition of two parts, =_b+_w, where _b represents the contribution from the fin body and _w represents the contribution from the wake. The problem is thus decomposed into two interdependent ones: the boundary value problem for _b and the evolution of the wake that determines _w.

The boundary-value problem for _b consists of the Laplace equation ²_b=0 inside the fluid domain, the far-field condition |_b|0 as |x|, and the no-flux condition on the fin surface S_b:

(B1)

where V_b is the instantaneous velocity at the surface of the fin, n is the unit normal vector pointing into the fin. The integral equation for _b is then obtained by invoking Green's theorem. At any point x on the fin surface we have:

(B2)

Note that the term on the right-hand side is known since the value of n·_b on S_b is given through Eqn B1. Solving the integral equation (B2), we obtain the unknown potential _b on S_b. The value of _b at any point x inside the fluid is then evaluated via Green's theorem, i.e.:

(B3)

To determine _w, we model the wake as a distribution of dipoles on a sheet shed from the sharp trailing edge of the fin. With the continuous generation of vorticity from the trailing edge, this wake sheet elongates over time. The strength of newly shed wake is determined by the Kutta condition, i.e. the strength of the newly shed wake _w equals the instantaneous difference of body influence potential _b between the upper and lower surfaces near the edge so that the local singularity is eliminated. The rest of the wake, treated as a material surface, is carried downstream by the collective effect of the upcoming flow and the self-induced velocity, and its strength is kept unchangeable due to the absence of dissipation. In order to stabilize the evolution of the discretely distributed vorticity in the wake, it is necessary to introduce a desingularization algorithm, in which we assign a finite core radius to each (body or wake) panel of dipole distribution (Krasny, 1986). It has been demonstrated that the hydrodynamic forces on the body are not sensitive to the value of (Zhu et al., 2002). The desingularization also eliminates singularity in the wake vorticity distribution, allowing us to visualize the near-body flow field.

A low-order boundary-element algorithm is employed to find the unknown body influence potential _b on S_b by solving Eqn B2 (Katz and Plotkin, 1991). To achieve this, the fin surface S_b is discretized into N_b panels. The low-order panel method refers to the fact that in this algorithm within each panel S_bj (j=1,..., N_b), the body potential _b and its normal derivative are assumed to be constants. By collocating at the centroids of the panels, Eqn B2 is transformed into a system of linear equations and can be numerically solved. Afterwards, a forward Euler algorithm is applied for time integrations. As an illustration, the distribution of boundary elements on the fin surface and the wake is demonstrated in Fig. B1. The flow potential is then determined by summing up the body influence potential _b and the wake influence potential _w.

View larger version (22K): [in this window] [in a new window]   Fig. B1. Distribution of boundary elements on the fin surface and the wake.

LIST OF SYMBOLS AND ABBREVIATIONS

c

chord length of the foil

C

rotation matrix

d

thickness of the fin

D

material damping coefficient of the rays

E

Young's modulus of the rays

EI

bending stiffness of the rays

F

hydrodynamic force on the fin

F_h

hydrodynamic force on rays

F_r^(r)

first-harmonic transverse force

F_t(_t)

(mean) thrust force

GJ

torsional stiffness of the rays

(i, j, k)

unit vectors of the global coordinate system

N_b

number of boundary elements on the fin surface

N_r

number of grids on a ray

p

hydrodynamic pressure

P()

(mean) power spent

s

distance from a point on a ray to the basal end measured along the unstretched ray

s₀

span at the leading edge

S_t

Strouhal number

T

period

T=(T, T_n, T_b)

internal force inside a ray

U

forward speed

V=(V, V_n, V_b)

translational velocity of a ray

(x, y, z)

global coordinate system

y₀

sway amplitude

β0, β1, β2, β3

Euler parameters

strain in a ray

t

time step

propulsion efficiency

0i

yaw amplitude of the i-th ray

(, n, b)

tangential, normal and bi-normal unit vectors in the Lagrangian coordinate system

b

potential due to the fin

w

potential due to the wake

flow potential

i

phase between sway and yaw of the i-th ray

frequency

=(, n, b)

rotational velocity of a ray

=(, n, b)

Darboux vector measuring material torsion and curvatures of a ray

	References

TOP Summary INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION AND CONCLUSIONS APPENDIX A. MODEL OF... References

 

Alben, S., Shelley, M. and Zhang, J. (2002). Drag reduction through self-similar bending of a flexible body. Nature 420,479 -481.[CrossRef][Medline]

Alben, S., Shelley, M. and Zhang, J. (2004). How flexibility induces streamlining in a two dimensional flow. Phys. Fluids 16,1694 -1713.[CrossRef]

Alben, S., Madden, P. G. and Lauder, G. V. (2007). The mechanics of active fin-shape control in ray-finned fishes. J. R. Soc. Interface 4, 243-256.[Abstract/Free Full Text]

Bainbridge, R. (1963). Caudal fin and body movement in the propulsion of some fish. J. Exp. Biol. 40, 23-56.[Abstract]

Birch, J. M., Dickson, W. B. and Dickinson, M. H. (2004). Force production and flow structure of the leading edge vortex on flapping wings at high and low Reynolds numbers. J. Exp. Biol. 207,1063 -1072.[Abstract/Free Full Text]

Brucker, C. and Bleckmann, H. (2007). Vortex dynamics in the wake of a mechanical fish. Exp. Fluids 43,799 -810.[CrossRef]

Cheng, J. Y., Pedley, T. J. and Altringham, J. D. (1998). A continuous dynamic beam model for swimming fish. Philos. Trans. R. Soc. Lond. B Biol. Sci. 353,981 -997.[Abstract/Free Full Text]

Combes, S. A. and Daniel, T. L. (2003a). Flexible stiffness in insect wings. I. Scaling and the influence of wing venation. J. Exp. Biol. 206,2979 -2987.[Abstract/Free Full Text]

Combes, S. A. and Daniel, T. L. (2003b). Flexible stiffness in insect wings. II. Spatial distribution and dynamic wing bending. J. Exp. Biol. 206,2989 -2997.[Abstract/Free Full Text]

Connell, B. S. H. (2006). Numerical investigation of the flow-body interaction of thin flexible foils and ambient flow. PhD thesis, Massachusetts Institute of Technology, USA.

Connell, B. S. H. and Yue, D. K. P. (2007). Flapping dynamics of a flag in a uniform stream. J. Fluid Mech. 581,33 -67.[CrossRef]

Fish, F. E. and Lauder, G. V. (2006). Passive and active flow control by swimming fishes and mammals. Annu. Rev. Fluid Mech. 38,193 -224.[CrossRef]

Furman, A. and Breitsamter, C. (2004). Delta wing steady pressure investigations for sharp and rounded leading edges. In New Results in Numerical and Experimental Fluid Mechanics V: Contributions to the 14th STAB/DGLR Symposium Bremen, Germany (ed. H.-J. Rath, C. Holze, H.-J. Heinemann, R. Henke and H. Hönlinger), pp. 77-84. Berlin, Heidelberg: Springer.

Gibb, A. C., Dickson, K. A. and Lauder, G. V. (1999). Tail kinematics of the chub mackerel Scomber japonicus: testing the homocercal tail model of fish propulsion. J. Exp. Biol. 202,2433 -2447.[Abstract]

Harder, W. (1975). Anatomy of Fishes. Stuttgart: E. Schweizerbart'sche Verlagsbuchhandlung.

Heathcote, S., Martin, D. and Gursul, I. (2004). Flexible flapping airfoil propulsion at zero freestream velocity. AIAA J. 42,2196 -2204.[CrossRef]

Kardong, K. V. (1998). Vertebrates: Comparative Anatomy, Function, Evolution (2nd edn). Dubuque: W. C. Brown.

Katz, J. and Plotkin, A. (1991). Low-speed Aerodynamics: From Wing Theory to Panel Methods. New York: McGraw-Hill.

Katz, J. and Weihs, D. (1978). Hydrodynamic propulsion by large amplitude oscillation of an airfoil with chordwise flexibility. J. Fluid Mech. 88,485 -497.[CrossRef]

Krasny, R. (1986). Desingularization of periodic vortex sheet roll-up. J. Comput. Phys. 65,292 -313.[CrossRef]

Lauder, G. V. (1989). Caudal fin locomotion in ray-finned fishes: historical and functional analyses. Am. Zool. 29,85 -102.

Lauder, G. V. (2000). Function of caudal fin during locomotion in fishes: kinematics, flow visualization, and evalutional patterns. Am. Zool. 40,101 -122.[CrossRef]

Lauder, G. V. and Drucker, E. G. (2004). Morphology and experimental hydrodynamics of fish fin control surfaces. IEEE J. Ocean Eng. 29,556 -571.[CrossRef]

Lauder, G. V. and Madden, P. G. A. (2007). Fish locomotion: kinematics and hydrodynamics of flexible foil-like fins. Exp. Fluids 43,641 -653.[CrossRef]

Liao, J. C., Beal, D. N., Lauder, G. V. and Triantafyllou, M. S. (2003). Fish exploiting vortices decrease muscle activity. Science 302,1566 -1569.[Abstract/Free Full Text]

Liu, H. and Kawachi, K. (1998). A numerical study of insect flight. J. Comput. Phys. 46,124 -156.

Liu, P. and Bose, N. (1997). Propulsive performance from oscillating propulsors with spanwise flexibility. Proc. R. Soc. Lond. A Math. Phys. Sci. 453,1763 -1770.

Mittal, R., Dong, H., Bozkurttas, M., Lauder, G. V. and Madden, P. (2006). Locomotion with flexible propulsors. II. Computational modeling of pectoral fin swimming in sunfish. Bioinspir. Biomim. 1,35 -41.[CrossRef]

Mohandas, N. and Evans, E. (1994). Mechanical properties of the red cell membrane in relation to molecular structure and genetic defects. Annu. Rev. Biophys. Biomol. Struct. 23,787 -818.[Medline]

Nauen, J. C. and Lauder, G. V. (2002). Hydrodynamics of caudal fin locomotion by chub mackerel, Scomber japonicus (Scombridae). J. Exp. Biol. 205,1709 -1724.[Abstract/Free Full Text]

Newman, D. J. S. and Wootton, R. J. (1986). An approach to the mechanics of pleating in dragonfly wings. J. Exp. Biol. 125,361 -372.[Abstract/Free Full Text]

Pederzani, J. and Haj-Hariri, H. (2006). Numerical analysis of heaving flexible airfoils in a viscous flow. AIAA J. 44,2773 -2779.[CrossRef]

Pedley, T. J. and Hill, S. J. (1999). Large-amplitude undulatory fish swimming: fluid mechanics coupled to internal mechanics. J. Exp. Biol. 202,3431 -3438.[Abstract]

Pullin, D. I. and Wang, Z. J. (2004). Unsteady forces on an accelerating plate and application to hovering insect flight. J. Fluid Mech. 509,1 -21.[CrossRef]

Ramamurti, R. and Sandberg, W. C. (2002). A three-dimensional computational study of the aerodynamic mechanisms of insect flight. J. Exp. Biol. 205,1507 -1518.[Abstract/Free Full Text]

Ramamurti, R. and Sandberg, W. C. (2007). A computational investigation of the three-dimensional unsteady aerodynamics of Drosophila hovering and maneuvering. J. Exp. Biol. 210,881 -896.[Abstract/Free Full Text]

Ramamurti, R., Sandberg, W. C., Lohner, R., Walker, J. A. and Westneat, M. W. (2005). Fluid dynamics of flapping aquatic flight in the bird wrasse: three-dimensional unsteady computations with fin deformation. J. Exp. Biol. 205,2997 -3008.

Tangorra, J. L., Davidson, S. N., Hunter, I. W., Madden, P. G. A., Lauder, G. V., Dong, H., Bozkurttas, M. and Mittal, R. (2007). The development of a biologically inspired propulsor for unmanned underwater vehicles. IEEE J. Ocean Eng. 32,533 -550.[CrossRef]

Tjavaras, A. A., Zhu, Q., Liu, Y., Triantafyllou, M. S. and Yue, D. K. P. (1998). The mechanics of highly-extensible cables. J. Sound Vibration 213,709 -737.[CrossRef]

Triantafyllou, M. S., Triantafyllou, G. S. and Gopalkrishnan, R. (1991). Wake mechanics for thrust generation in oscillating foils. Phys. Fluids 3,2835 -2837.

Triantafyllou, M. S., Techet, A. H. and Hover, F. S. (2004). Review of experimental work in biomimetic foils. IEEE J. Ocean Eng. 29,585 -594.[CrossRef]

Tuncer, I. H. and Platzer, M. F. (1996). Thrust generation due to airfoil flapping. AIAA J. 34,324 -331.[CrossRef]

Tytell, E. D. (2006). Median fin function in bluegill sunfish Lepomis macrochirus: streamwise vortex structure during steady swimming. J. Exp. Biol. 209,1516 -1534.[Abstract/Free Full Text]

Wang, Z. J. (2000). Vortex shedding and frequency selection in flapping flight. J. Fluid Mech. 410,323 -341.[CrossRef]

Wilga, C. D. and Lauder, G. V. (2002). Function of the heterocercal tail in sharks: quantitative wake dynamics during steady horizontal swimming and vertical maneuvering. J. Exp. Biol. 205,2365 -2374.[Abstract/Free Full Text]

Wilga, C. D. and Lauder, G. V. (2004). Hydrodynamic function of the shark's tail. Nature 430, 850.[CrossRef][Medline]

Wootton, R. J. (1992). Functional morphology of insect wings. Annu. Rev. Entomol. 37,113 -140.[CrossRef]

Yamamoto, I., Terada, Y., Nagamatu, T. and Imaizumi, Y. (1995). Propulsion system with flexible rigid oscillating fin. IEEE J. Ocean Eng. 20,23 -30.[CrossRef]

Zhu, Q. (2007). Numerical simulation of a flapping foil with chordwise or spanwise flexibility. AIAA J. 45,2448 -2457.[CrossRef]

Zhu, Q., Wolfgang, M., Yue, D. K. P. and Triantafyllou, M. S. (2002). Three-dimensional flow structures and vorticity control in fish-like swimming. J. Fluid Mech. 468, 1-28.[CrossRef]

Zhu, Q., Liu, Y. and Yue, D. K. P. (2006). Dynamics of a three-dimensional oscillating foil near the free surface. AIAA J. 44,2997 -3009.[CrossRef]

CiteULike    Complore    Connotea    Del.icio.us    Digg    Reddit    Technorati    Twitter    What's this?

Related articles in JEB:

FLEXIBLE FINS BEAT RIGID FINS

Kathryn Phillips
JEB 2008 211: i. [Full Text]  

This article has been cited by other articles:

				K. Shoele and Q. Zhu Fluid-structure interactions of skeleton-reinforced fins: performance analysis of a paired fin in lift-based propulsion J. Exp. Biol., August 15, 2009; 212(16): 2679 - 2690. [Abstract] [Full Text] [PDF]

				B. E. Flammang and G. V. Lauder Caudal fin shape modulation and control during acceleration, braking and backing maneuvers in bluegill sunfish, Lepomis macrochirus J. Exp. Biol., January 15, 2009; 212(2): 277 - 286. [Abstract] [Full Text] [PDF]

				K. Phillips FLEXIBLE FINS BEAT RIGID FINS J. Exp. Biol., July 1, 2008; 211(13): i - ii. [Full Text] [PDF]

This Article Summary Figures Only Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Related articles in JEB Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Zhu, Q. Articles by Shoele, K. Search for Related Content PubMed PubMed Citation Articles by Zhu, Q. Articles by Shoele, K. Social Bookmarking                         What's this?