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

First published online January 31, 2007
Journal of Experimental Biology 210, 685-698 (2007)
Published by The Company of Biologists 2007
doi: 10.1242/jeb.02692

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 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 Peng, J. Articles by Lauder, G. V. Search for Related Content PubMed PubMed Citation Articles by Peng, J. Articles by Lauder, G. V.

Non-invasive measurement of instantaneous forces during aquatic locomotion: a case study of the bluegill sunfish pectoral fin

Jifeng Peng¹, John O. Dabiri^1,2,*, Peter G. Madden³ and George V. Lauder³

¹ Bioengineering, California Institute of Technology, Pasadena, CA 91125, USA
² Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena, CA 91125, USA
³ Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA

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

Accepted 12 December 2006

	Summary

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 
Swimming and flying animals generate unsteady locomotive forces by delivering net momentum into the fluid wake. Hence, swimming and flying forces can be quantified by measuring the momentum of animal wakes. A recently developed model provides an approach to empirically deduce swimming and flying forces based on the measurement of velocity and vortex added-mass in the animal wake. The model is contingent on the identification of the vortex boundary in the wake. This paper demonstrates the application of that method to a case study quantifying the instantaneous locomotive forces generated by the pectoral fins of the bluegill sunfish (Lepomis macrochirus Rafinesque), measured using digital particle image velocimetry (DPIV). The finite-time Lyapunov exponent (FTLE) field calculated from the DPIV data was used to determine the wake vortex boundary, according to recently developed fluid dynamics theory. Momentum of the vortex wake and its added-mass were determined and the corresponding instantaneous locomotive forces were quantified at discrete time points during the fin stroke. The instantaneous forces estimated in this study agree in magnitude with the time-averaged forces quantified for the pectoral fin of the same species swimming in similar conditions and are consistent with the observed global motion of the animals. A key result of this study is its suggestion that the dynamical effect of the vortex wake on locomotion is to replace the real animal fin with an `effective appendage', whose geometry is dictated by the FTLE field and whose interaction with the surrounding fluid is wholly dictated by inviscid concepts from potential flow theory. Benefits and limitations of this new framework for non-invasive instantaneous force measurement are discussed, and its application to comparative biomechanics and engineering studies is suggested.

Key words: swimming, force, locomotion, fish, particle image velocimetry, pectoral fin, vortex, fluid dynamics

	Introduction

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 
The challenge of quantifying the locomotive forces exerted by fluids on swimming and flying animals has attracted the interest of both biologists and engineers for many years. Early investigations of swimming and flying animals estimated locomotive forces from kinematic data, due to the lack of experimental techniques to quantitatively measure properties of the fluid wake. Quantification and validation of locomotive forces using these approaches relies heavily on semi-empirical models (for reviews, see Webb and Blake, 1985; Vogel, 1994).

Recent studies have taken an alternative approach. Newton's second and third laws dictate that flying and swimming animals generate net locomotive forces by transferring momentum into the wake. The hallmark of the momentum transfer is the creation of vortices. Hence, studies have aimed to quantify the momentum of fluid vortices surrounding animals in order to estimate the associated locomotive forces. Over the past decade, the development of the flow visualization technique digital particle image velocimetry (DPIV) has enabled researchers to realize the direct visualization of flow in the wakes of animals and to measure the corresponding velocity and vorticity fields (e.g. Drucker and Lauder, 1999; Lauder and Drucker, 2002; Nauwelaerts et al., 2005; Spedding et al., 2003; Warrick et al., 2005). These fluid dynamic data have provided information necessary for the estimation of locomotive forces based on kinematics and dynamics of the animal wake.

Several models have been proposed to evaluate locomotive forces using the data from wake measurement data. These studies typically estimate the time-averaged force rather than instantaneous forces. For example, studies have estimated the fluid forces based on measurements of the near-appendage circulation of vortices created by the animals (Dickinson, 1996; Dickinson and Götz, 1996). Locomotive forces experienced by the animal are calculated as the reaction to the momentum of vortex loops shed into the wake (e.g. Drucker and Lauder, 1999; Drucker and Lauder, 2001; Johansson and Lauder, 2004; Stamhuis and Nauwelaerts, 2005). In these cases, the momentum of the vortex is usually measured at the time instant when the vortex ring has just detached from the animal fin/wing. The time-averaged locomotive force over the stroke cycle is then determined by dividing the momentum of the shed vortex by the time duration of the stroke cycle. In other studies, the locomotive forces have been evaluated by examining the wake far downstream, which is equivalent to taking the time-average of what occurred at the site of force generation (e.g. Spedding et al., 2003) (cf. Walker, 2004).

As discussed elsewhere (Peng and Dabiri, 2007a), viscous dissipation and vorticity cancellation will reduce the efficiency of the momentum transfer process from 100%, resulting in an `information loss' in the record of locomotive dynamics contained in by wake. However, a straightforward viscous scaling argument shows that these effects are usually negligible on the time scale of individual stroke cycles. In particular, the distance over which viscosity will act during a single stroke of duration T_S goes as (T_S)^1/2, where is the kinematic viscosity of the fluid (Rosenhead, 1963). Regions of opposite-signed vorticity (e.g. shed from the dorsal and ventral edges of a pectoral or caudal fin) must be within this distance (T_S)^1/2 from each other in order to undergo vorticity cancellation and the associated information loss in the wake. For repeated swimming or flying motions at frequency f_S, the scaling is equivalently (/f_S)^1/2. Hence, information loss in the wake becomes important if the ratio /L(/f_S)^1/2/L is of order one or larger, where L is the characteristic length scale of the appendage. A 1 Hz swimming motion in water (10^–2 cm² s^–1) corresponds to a characteristic viscous length scale 1 mm, which is substantially smaller than the length scales of most fish appendages (although not necessarily small for swimming microorganisms). In air (10^–1 cm² s^–1) at 1 Hz, 3 mm, which is also smaller than the length scales of most bird appendages. Insects may have appendage length scales on this order, but will also operate at much higher frequencies. For example, in Drosophila (/f_S)^1/2/L10^–1, indicating a limited role for vorticity cancellation during the transfer of momentum from the animal to the wake in a single stroke. Therefore, for the near-wake (vis-à-vis far downstream) studies of concern here, we will assume no information loss between the dynamics of the animal and the wake it generates.

Previous studies analyzing wake vorticity have found the measured time-averaged forces to be comparable with the necessary [but not sufficient (Dabiri, 2005)] lift and thrust required to sustain flying and swimming. However, time-averaged forces provide little information about the dynamics of swimming and flying. It is the instantaneous forces that dictate important dynamics of locomotion such as the trajectory, speed and efficiency of swimming and flying. In addition, these previous methods implicitly assume that the flow is steady so that the vortex momentum can be determined from the distribution of vorticity alone. For these reasons, an approach toward the task of estimating instantaneous, unsteady locomotive forces of freely moving animals was recently developed (see Dabiri, 2005; Dabiri, 2006; Dabiri et al., 2006). This model provides a method to empirically deduce unsteady swimming and flying forces based on the measurement of velocity and vortex added-mass in the animal wake. In this method, the momentum of the vortex in animal wake is evaluated as the sum of the linear momentum of fluid inside the vortex and the linear momentum of fluid surrounding the wake vortex, i.e. the added-mass of the wake vortex. Given velocity field measurements in the wake, the vortex boundary in the wake can be determined and the momentum of the wake vortex and its added-mass can be calculated, leading to a quantitative evaluation of instantaneous locomotive forces.

In the present study, this method is applied to quantify the instantaneous locomotive forces generated by the pectoral fins of the bluegill sunfish Lepomis macrochirus Rafinesque during labriform locomotion. The overall goal of this study is to develop a case study in which the aforementioned method can be used to analyze velocity field measurements in order to estimate instantaneous locomotive forces in freely moving animals. We use a new data set from pectoral fin locomotion in a well-studied fish, the bluegill sunfish, to explore the application of this method to fish locomotion. Due to previous work on this species using the time-averaged vortex approach (e.g. Drucker and Lauder, 1999; Drucker and Lauder, 2000; Drucker and Lauder, 2001), good estimates of stroke averaged forces are available for comparison to the instantaneous forces calculated here. In addition, the mean vertical force is known, as the weight of sunfish underwater is a previously measured quantity (Drucker and Lauder, 1999), which allows a further check on the instantaneous force calculations. Furthermore, we aim to explore the utility of this method with a data set of the kind that is typically available to investigators studying animal locomotion: a time series of two-dimensional DPIV vector fields.

The outline of our approach is as follows: (1) the time-dependent velocity field of the wake in freely swimming sunfish is measured at high temporal and spatial resolution using DPIV; (2) a dynamical systems theory developed recently (Haller, 2000; Haller, 2001; Shadden et al., 2005; Shadden et al., 2006) is used to identify the boundary of the vortex wake on the two-dimensional (2-D) plane from the velocity field data; (3) a three-dimensional (3-D) approximation based on the 2-D vortex boundary is used determine the momentum of the vortex itself and its added-mass; (4) the corresponding locomotive forces acting on the pectoral fins are calculated. Unlike previous studies that estimated the time-averaged forces over the stroke cycle, this study provides detailed information on how locomotive forces evolve within the fin stroke cycle.

Like most studies involving time-averaged locomotive force estimation of swimming and flying animals, analysis in the present study is based on the flow velocity data from 2-D DPIV measurements. 3-D velocity field data with a large control volume including both the appendage and fish body would give more accurate force estimation [e.g. using the exact equations of motion derived by Noca et al. (Noca et al., 1997; Noca et al., 1999)], but such measurements are not yet currently possible. Limitations associated with the use of flow velocity data from a 2-D DPIV measurement to evaluate instantaneous locomotive forces are investigated in this study, and directions for future studies are suggested to enable more accurate force measurements.

	Materials and methods

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 
Fishes and experimental conditions
Bluegill sunfish Lepomis macrochirus Rafinesque (N=7, total body length L=17.3±1.0 cm, mean ± s.d.) were collected by seine from ponds in Concord, MA, USA and maintained at an average temperature of 20°C in 40 l aquaria.

In experiments, sunfish swam individually in the center of the working area (28 cmx28 cmx80 cm) of a variable-speed freshwater flow tank under conditions similar to those described previously (Drucker and Lauder, 1999; Drucker and Lauder, 2000; Drucker and Lauder, 2001; Drucker and Lauder, 2003). The sunfishes were trained to hold station in a current with a velocity of 0.5 L s^–1. At this relatively low speed, sunfish swimming usually involves use of the pectoral fins to generate locomotive forces (labriform locomotion). Only steady rectilinear swimming, during which the fish maintained a speed within 5% of the flow tank's current speed, was considered for analysis. To minimize wall effects, the fish were required to swim near the center of the volume of the working area. Thus the flow structures are assumed to result directly from movements of the fish pectoral fins.

Wake visualization and measurement
DPIV was used to visualize and measure the wake of the sunfish pectoral fin. General details of the method are provided elsewhere (Drucker and Lauder, 1999), although the data generated for this paper resulted from a modified approach to produce a temporally and spatially much higher resolution data set. The DPIV technique provides empirical velocity field data for flow in two-dimensional sections of the swimming fish wake (for details, see Willert and Gharib, 1991; Drucker and Lauder, 1999; Lauder, 2000). A 10 W continuous-wave argon-ion laser (Coherent Inc., Santa Clara, CA, USA) was focused into a thin light sheet 1–2 mm thick and 10 cm wide, which illuminated reflective silver-coated glass spheres (mean diameter 12 µm, density 1.3 g cm^–3) suspended in the water (concentration 14 mg l^–1). Particle motion induced by pectoral fin activity was recorded by imaging the laser sheet with a high-speed video camera (Photron Fastcam-ultima APX, 1024x1024 pixels at 500 frames s^–1); a second camera (Photron Fastcam-X 1280PCI, 1280x1024 pixels, 500 frames s^–1) synchronously recorded a perpendicular reference view showing the position of the fin relative to the visualized transection of the wake. In the majority of experiments described herein, the laser was oriented to reveal the flow in the transverse plane, which crosses the left pectoral fin (Fig. 1) and intersects the sunfish body perpendicularly. This provided a full field view of the pectoral fin wake moving toward the camera, and the high temporal sampling rate provided a detailed image of pectoral fin wake flow patterns as the wake moved through the laser plane. Camera images of the transverse light sheet plane were obtained through a mirror located in the flow tank, well downstream (5–10 fin chord lengths) of the sunfish pectoral fin. A total of 55 sequences were obtained for analysis from the seven sunfish studied. Considerable effort was made to obtain DPIV images in which the fish was swimming as steadily as possible, not maneuvering, and in which the fin was located in a variety of locations relative to the horizontal light sheet. In some sequences, we maneuvered the fish into a location where the laser sheet sliced through the fin itself so that the boundaries of vortex structures relative to the fin as it moved during the fin beat cycle could be determined. The sunfish pectoral fin is translucent and the small supporting fin rays do not significantly obstruct the laser light, so no shadows were formed that might inhibit vector calculation between the fin and the body. In other sequences, sunfish were induced to swim upstream of the light sheet so that the pectoral fin wake alone was imaged.

View larger version (9K): [in this window] [in a new window]   Fig. 1. The side and top views of a bluegill sunfish and the laser plane (shown in blue) for the DPIV experiments that generated the data used for the analysis in this paper. Note that the camera viewed flow from behind the fish at 500 frames s–1, and pectoral fin wakes thus move toward the field of view allowing a complete view of the wake at high temporal resolution.

We focused the analysis in this paper on the vertical (lift or dorsal–ventral) and side (lateral) forces, for several reasons. (1) Time averaged calculations of sunfish pectoral fin forces in both of these directions are available from previous work (Drucker and Lauder, 1999); (2) mean side forces should be near zero if the sunfish is swimming steadily, and the mean lift force should balance the weight of sunfish in water [sunfish are slightly negatively buoyant (Drucker and Lauder, 1999)], providing validation of the calculated instantaneous forces; (3) the DPIV data in the transverse plane were most accurately calculated in the side and vertical directions, while thrust data estimates were not available as flow was moving through the imaging plane.

2-D velocity fields in the wake of the swimming sunfishes were calculated from consecutive digital video images (1024x1024 pixels, 8-bit grayscale) using DaVis software (LaVision GmbH, Göttingen, Germany). Image spatial cross-correlation was performed with convolution filtering iterating from a 64x64 pixel to a 16x16 pixel interrogation area with 50% window overlap, giving about 10 000 vectors per time sample. Erroneous vectors and outliers were removed automatically by removing vectors more than 2 standard deviations away from the mean of the neighbors. Vectors were smoothed with a 3x3 vector averaging filter.

For the purpose of demonstration in this paper, only one representative set of velocity data was chosen for the following analyses. Accordingly, the results shown correspond to the specific set of data analyzed, not a composite analysis of all of the measurement samples.

Force estimation
Locomotive forces experienced by the fin were calculated as the reaction to the total momentum imparted to the vortex wake. The momentum of the vortex wake consists of two components: the linear momentum of the fluid inside the vortex and the linear momentum of the fluid surrounding the vortex that moves in the same direction as the vortex. The latter component is the added-mass of the vortex in the wake and is identical to the added-mass traditionally associated with fluid surrounding solid bodies in potential flow (Dabiri, 2006). The expression for the wake momentum and its derivation were given previously (Dabiri, 2005) [see also Dabiri et al. (Dabiri et al., 2006) for correction]. Quantification of the momentum of the wake requires not only experimental data with sufficient spatial and temporal resolution, but also an involved mathematical analysis. Thus, approximations and simplifications were made where possible.

If the wake vortex does not deform rapidly, the impulse I of the fluid circulating inside the vortex can be simplified as:

(1)

where is water density (1000 kg m^–3 at 20°C), V_V the volume of the vortex and U_V the velocity of the wake vortex center of mass. Furthermore, the impulse of the wake vortex added-mass can be approximated based on the added-mass tensor C of the vortex, as well as the volume and velocity of each vortex as it is formed in the wake:

(2)

The contribution of the added-mass of the vortex depends on the vortex shape and the type of motion the vortex is experiencing (i.e. linear or angular). Depending on whether or not bulk rotational motion of the vortex (i.e. rotation of the principal axes of the vortex volume) is accounted for, the velocity U_V will be either a 3x1 or 6x1 vector, and the vortex added-mass tensor C is correspondingly a 3x3 or 6x6 matrix, respectively, with elements c_ij, which are the dimensionless added-mass coefficients that relate acceleration in the ith direction to the resultant forces in the jth direction (where i and j can assume translation in x-, y- and z-axis directions, or rotation in the xy-, xz- and yz-planes. Repeated subscripts c_ii do not indicate summation). For the purpose of demonstration in this paper, the added-mass effect from bulk rotational motion of the vortex is not considered here. Also, the added-mass effect is assumed to be dominated by the c_ii added-mass coefficients [suggesting the existence of two planes of vortex spatial symmetry (Dabiri et al., 2006)]; therefore, the off-diagonal elements of the added-mass tensor C are neglected. The total impulse I of the moving vortex can then be simplified as:

(3)

where the scalar C represents the added-mass coefficient for linear motion along the direction of U_V. By Newton's second and third laws, the locomotive force exerted by the fluid on the fin is equal and opposite to the rate at which the momentum of the vortex I changes due to the interaction between the fluid and the fin:

(4)

To make the expression in Eqn 4 compatible with the format of typical wake measurements, the time derivative can be written in terms of the finite difference of data collected at discrete time points t_j, with t_j+1=t_j+t and j=0, 1, 2, etc:

(5)

As t decreases to zero, Eqn 5 becomes an estimate of the instantaneous force generated by the swimming or flying animal. As t increases to T, the duration of the propulsive stroke, Eqn 5 becomes an estimate of the time-averaged locomotive force. In this study, t=10 ms was used. For comparison, the duration of an entire stroke cycle was approximately 600 ms. In practice, the duration t is only limited by the temporal resolution of the DPIV data.

Vortex boundary identification
It can be seen from Eqn 4 that the estimate of locomotive forces requires the determination of the physical boundary separating the vortex from the surrounding flow. The volume, velocity and the added-mass coefficient of the vortex structure all depend on identification of the vortex boundary. For steady flows, the flow structure can be generally identified by examining streamlines derived from velocity fields measured by DPIV, because these streamlines are also fluid particle trajectories. For the unsteady flows of most animal wakes, however, defining the boundary between a vortex and the surrounding flow is not an obvious task. In some simple unsteady flows, streamlines may still reveal the boundary of the vortex if plotted in a reference frame moving with the vortex (Dabiri and Gharib, 2004). While this method can be effective for studying radially symmetric vortex ring wakes such as those generated by some jellyfish [cf. fig. 9 in Dabiri (Dabiri, 2005)], squids and salps, it cannot be used to elucidate the structure of more complex wakes of other swimming and flying animals (Dabiri, 2005).

A more objective, frame-independent technique was recently developed to determine vortex boundaries in unsteady vortex flows measured empirically (Shadden et al., 2006). In this approach, the boundary of the vortex is determined by tracking the relative Lagrangian trajectories of individual fluid particles in the flow, rather than by analyzing the velocity or vorticity fields or streamlines at each time instant as in the Eulerian perspective. The method is briefly summarized here; however, the reader is directed to the papers of Shadden et al. for greater technical detail (Shadden et al., 2005; Shadden et al., 2006).

Given a time-dependent velocity field u(x,t), the trajectory of a fluid particle x(t) can be determined by the ordinary differential equation:

(6)

with given initial conditions. The flow map, which maps fluid particles from their initial location at time t to their location at time t+T, can be expressed as:

(7)

where describes the current location of a fluid particle advected from the location x(t) at time t after a time interval T. Consider two adjacent fluid particles x(t) and y(t)=x(t)+x(0) in the flow at time t, where x(0) is infinitesimal. Their locations after a time interval T are and . The distance between the two fluid particles at time t+T is therefore:

(8)

where the second equality comes from taking the Taylor series expansion of the flow about point x, and the order O of the last term is indicated in parentheses. By dropping the small term O[||x(0)||²] in Eqn 8, a parameter is introduced that represents the rate of change of the distance between two initially adjacent fluid particles:

(9)

The parameter , which is called the finite-time Lyapunov exponent (FTLE), measures the linearized growth rate of the small perturbation x over the interval T of trajectories starting near x(t). In other words, it characterizes the amount of fluid particle separation, or stretching, about the trajectory of point x over the time interval [t,t+T].

View larger version (8K): [in this window] [in a new window]   Fig. 2. A schematic drawing of a vortex boundary in a flow. Circles with inscribed arrows indicate vortex cores and their rotational sense. A pair of adjacent fluid particles close to but on different sides of the vortex boundary separate from each other faster than other arbitrary pairs of fluid particles, giving a larger value of the FTLE at the boundary. Trajectories can be followed in backward-time (A) to reveal attracting LCS, and in forward-time (B) to reveal repelling LCS boundaries.

It is important to note that though the FTLE is a function of position variable x and time t, it is thought of as a Lagrangian quantity since it is derived from fluid particle trajectories over the time interval [t,t+T]. The absolute value |T| is used instead of T in Eqn 10 because FTLE can be computed for T>0 and T<0 (Fig. 2). Forward-time integration (T>0) reveals repelling LCS because particles straddling this type of LCS separate faster than other arbitrary point pairs. Backward-time integration locates attracting LCS because particles straddling this type of LCS converge faster than other arbitrary point pairs when advected forward in time, i.e. separate faster when advected backward in time (Haller, 2001). The integration time |T| is chosen according to the particular flow being analyzed. If a smaller integration time is used, then less of the boundary is revealed, whereas if a longer integration time is used, more of the boundary is revealed. Generally, if the integration time |T| is sufficiently long, the repelling and the attracting LCS usually intersect to give the boundary of the vortex [in cases where a vortex is known to be present; cf. fig. 6 in Shadden et al. (Shadden et al., 2006)]. A larger integration time |T| also gives LCS with higher spatial resolution. However, the choice of |T| is sometimes limited in practice by the availability of data.

The FTLE fields in the transverse plane were calculated by using the LCS MATLAB Kit Version 1.0 (freeware download at http://dabiri.caltech.edu/software.html) to analyze the 2-D experimental DPIV data. A Cartesian grid was used for the FTLE computations in the study, with uniform spacing of 0.5 mm. The flow map at each node was calculated with a 4th-order Runge–Kutta integration algorithm. Since the velocity data are discrete, a 3rd-order spatial interpolation was used to provide necessary spatial resolution. An integration length |T|=200 ms was used due to the nature of the very short stroke cycle (about 600 ms); this limited availability of velocity data prevented a calculation with a larger integration time. Longer integration times result in more clearly defined maxima in the FTLE field, i.e. more well-defined vortex boundaries. Positive and negative integration time intervals were used to determine forward- and backward-time FTLE fields, respectively. Repelling and attracting LCS were determined by locating ridges (i.e. contours of local maxima) of the forward- and backward-time FTLE fields, respectively.

Model approximations
The boundary of the forming wake vortex in the transverse plane was determined by using the method described above. This 2-D boundary of the vortex was used to approximate the following quantities.

First, the 2-D vortex boundary on the transverse plane was used to estimate the vortex volume V_V. It has been shown in a previous study (Drucker and Lauder, 1999) that an isolated vortex ring is generated by the pectoral fin of the bluegill sunfish during labriform fin kinematics. However, the vortex boundary may not possess a simple geometry. To approximate the 3-D shape of the vortex, the calculated vortex boundary on the transverse plane was mapped into an ellipse, which represents the cross-section of the vortex in the transverse plane. The calculated vortex boundary and the model ellipse have the same long axis length and also the same width (Fig. 3). Assuming the elliptical vortex cross-section was perpendicular to and symmetrical about the transverse plane, the volume of the vortex can be approximated by:

(10)

where w is the width of the wake vortex on the laser plane view and A is the projected area of the vortex ring on the plane perpendicular to the laser plane. The added-mass coefficient C of the vortex body can be determined from an equivalent solid body calculation of the approximated 3-D ellipsoidal vortex shape (Lamb, 1932; Dabiri, 2006).

View larger version (9K): [in this window] [in a new window]   Fig. 3. Illustration of the approximation of the vortex volume. The calculated vortex boundary on the transverse plane was modeled by an ellipse, which represents the cross-section of the vortex in the transverse plane. The calculated vortex boundary and the model ellipse have the same long axis length and also the same width. The long axes of the calculated boundary and the approximated shape are parallel. The volume of the vortex is approximated by the product of the center-width w of the vortex in the laser plane view and the projected area A of the vortex ring onto a plane perpendicular to the laser sheet.

The motion of the LCS was approximated by the motion of the model ellipse. The long axes of the LCS and the model ellipse are always parallel; hence, the model is consistent with the known motion of a vortex ring in the direction normal to its long axis (Lamb, 1932). The projection of the vortex body velocity U_V on the transverse plane was approximated by the velocity of the centroid of the fluid inside the 2-D vortex boundary. Since the velocity component perpendicular to the transverse plane was not available, the locomotive forces in that direction were not evaluated.

Uncertainty analysis
The uncertainty in this study comes from velocity measurements using DPIV, interpolation of discretized velocity data when calculating FTLE field, identification of the vortex boundary from the FTLE field, and the 3-D approximation of the vortex structure based on the boundary information on a 2-D cross-sectional plane. The velocity measurement error using DPIV is on the order of 1–3% for flows at these velocities and accelerations (Willert and Gharib, 1991). This error includes uncertainty due to the fact that the seed particles in the flow are slightly negatively buoyant. Adrian (Adrian, 1995) shows that the characteristic particle lag time is a function of the fluid viscosity , particle diameter d_p, and particle-fluid density ratio . For the present experiments, _lag=41 µs, or 2% of the temporal frame spacing (i.e. 1 frame/500 frames s^–1=0.002 s) used to process the DPIV data.

Given that the FTLE is calculated from particle trajectories (integration of velocity), a reasonable concern is that local error in velocity measurement may accumulate upon integration in time. However, it has been rigorously shown (Haller, 2002) that even large local (in space and time) velocity errors do not prevent reliable calculations of the position of LCS, as long as the errors remain small in a special time-weighted norm.

Hence, the major source of uncertainty arises from identification of LCS from the FTLE calculations and from 3-D approximation of the wake geometry based on the 2-D measurements. An uncertainty of 1 mm is approximated for measurements of position and dimensions of the LCS. The uncertainty in parameters derived from these fundamental measurements is determined by the rules of error propagation, i.e.

(11)

and

(12)

where is the uncertainty (Taylor, 1997).

View larger version (154K): [in this window] [in a new window]   Fig. 4. Velocity and vorticity fields of the pectoral fin wake in the transverse plane. (A) Early downstroke; (B) late downstroke and stroke reverse; (C) early upstroke; (D) late upstroke. Red colors represent negative or clockwise fluid rotation, while blue colors indicate positive vorticity or counterclockwise fluid rotation. The camera recorded a posterior view of the left pectoral fin and the fish body. Every other vector is shown.

	Results

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

View larger version (66K): [in this window] [in a new window]   Fig. 5. A snapshot of color contour plots of the FTLE fields computed from DPIV. (A) Backward FTLE; (B) forward FTLE. Position coordinates are specified in mm.

View larger version (107K): [in this window] [in a new window]   Fig. 6. The boundary of the vortex derived from LCS. The left solid line shows the attracting LCS from backward FTLE calculation while the right-hand-side solid line shows the repelling LCS from forward FTLE calculation. Broken lines are spline lines connecting the LCS. The fin (the curve with high brightness inside the lines) can be seen embedded inside the vortex. The attracting and repelling LCS do not intersect to give the entire vortex boundary because of the limitation in integration time T. Every other vector is shown.

View larger version (18K): [in this window] [in a new window]   Fig. 7. Time evolution of the vortex boundary. Vortex boundaries at 11 different time instances are plotted from red to blue with a time interval of 30 ms.

The FTLE fields shown in Fig. 5 do not give a vortex boundary that is as sharply defined as in previous studies of isolated vortex rings [cf. fig. 6 in Shadden et al. (Shadden et al., 2006)]. This is due to the aforementioned limitation in integration time |T|. Nonetheless, the location of the LCS boundary can still be approximated to lie at the centerline of the FTLE contours in each frame. The backward- and forward-time LCS derived from Fig. 5A,B are plotted together to give the vortex boundary in Fig. 6. A best-fit spline connection was used if the backward- and forward-time LCS did not intersect in a given frame. Ideally the attracting and repelling LCS should always completely intersect to give a well-defined wake vortex boundary when a wake vortex is present. However, due to the aforementioned short integration time, in some frames less of the boundary is revealed by the LCS. In Fig. 6, the boundary of the vortex at the transverse plane is superimposed on the DPIV velocity field data at t=250 ms. Notice that it is impossible to define a vortex boundary from inspection of the velocity field alone, whereas the theory governing the LCS ensures that the boundary is captured by the present FTLE measurements (Shadden et al., 2005; Shadden et al., 2006).

The boundary of the vortex on the transverse plane does not have a regular elliptic shape, a possibility anticipated in the previous section. The fin (the portion with high brightness alongside the sunfish body in Fig. 6) can be seen embedded inside the vortex, i.e. the shape of the vortex corresponds to the shape and location of the fin. This result that the fin is enclosed in the vortex is in accordance with the fact that the pair vortices of the vortex ring are generated at the dorsal and the ventral edges of the fin. Given the fact that the fin density is only 10% higher than water and its thickness is on the order of 100 µm (Alben et al., 2007; Lauder et al., 2007), its contribution to the total mass of the vortex can be neglected in the subsequent force calculations.

View larger version (8K): [in this window] [in a new window]   Fig. 8. Trajectory of the projection of the vortex centroid on the transverse plane. Squares: calculated data at each time instance. Solid line: spline fitting of the data using a centered moving average method with a span of five data points. Error bars indicate measurement uncertainty. The designations are the same for Figs 9, 10, 11, 12.

View larger version (11K): [in this window] [in a new window]   Fig. 9. The evolution of (A) x and (B) y positions of the projection of the vortex centroid on the transverse plane. Error bars indicate measurement uncertainty. Note that due to limitations on FTLE integration time, these plots are for the first 400 ms of a 600 ms fin beat cycle; see also Figs 10, 11, 12. Squares, calculated position of the vortex centroid; solid line, spline fitting of the data.

View larger version (11K): [in this window] [in a new window]   Fig. 10. The velocity of the projection of the vortex centroid on the transverse plane. (A) Horizontal component Ux; (B) vertical component Uy.

View larger version (23K): [in this window] [in a new window]   Fig. 11. (A) Volume of the vortex; (B) width of the vortex; (C) cross-sectional area of the vortex; and (D) added-mass coefficient of the vortex.

View larger version (10K): [in this window] [in a new window]   Fig. 12. The locomotive force (mN) in (A) the horizontal and (B) the vertical directions. Squares, calculated locomotive forces; broken line, time-averaged forces calculated using the vorticity method (Drucker and Lauder, 1999).

The position of the projection of the vortex centroid on the transverse plane was determined by calculating the centroid of the area enclosed in the vortex boundary. The trajectory of the projected vortex centroid on the transverse plane is plotted in Fig. 8 and the change of the position with time is plotted in Fig. 9. The movement of the vortex is consistent with the fin kinematics. The velocity of the projected vortex centroid on the transverse plane was calculated from the trajectory by the change of position in unit time (Fig. 10). It can be seen from the velocity profile that during the downstroke, the vortex accelerates vertically most significantly during t=100 to 150 ms and accelerates horizontally most significantly during t=160–210 ms. It will be shown later that this pattern of movement plays a major role in the force generation.

View larger version (22K): [in this window] [in a new window]   Fig. 13. Analysis of the contributing factors to the locomotive force. (A) Horizontal direction; (B) vertical direction. The changes in the logarithm of each parameter to the vortex momentum are plotted as well as the total change in the vortex momentum, in order to determine the most dominant contributing parameters to the change in momentum.

Vortex added-mass
The added-mass coefficient of the vortex was determined by the approximated 3-D shape of the vortex ring (Lamb, 1932). The relationship between added-mass coefficient and time is plotted in Fig. 11D. In the early part of the downstroke, the vortex resembles a thin disk, giving a larger added-mass coefficient, while in the latter part of the downstroke, the vortex is more sphere-like, giving an added-mass coefficient with a value approaching 0.5, which is the added-mass coefficient of a sphere. Since the added-mass coefficient is determined by the shape of the vortex, it follows the change in vortex shape closely.

Locomotive force
The forces exerted by the fluid on the pectoral fin over the time span of interest are plotted in Fig. 12, with horizontal component (lateral force) and vertical component (lift force) separated. Also plotted are the time-averaged forces calculated using the vorticity method (Drucker and Lauder, 1999). The positive horizontal force is directed toward the body while the positive vertical force is directed upward. The data indicate that at the early downstroke there is a relatively large lateral force in the direction opposite to the fin movement. The fin also generates a significant lift at the early downstroke. In the late part of the downstroke, the horizontal force is close to zero and the magnitude of the vertical force also becomes much smaller than it was in the early part of the downstroke.

There is a phase difference between when lift force and lateral force are generated during the downstroke. Comparing Fig. 12A with Fig. 12B, it is apparent that the fin generates significant lift (at t=100–160 ms) before it generates lateral force (at t=160–210 ms). This effect arises because the acceleration of the vortex in the vertical direction occurs prior to its acceleration in the horizontal direction (Fig. 10).

	Discussion

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 
In this paper we present a case study to illustrate a framework for combining traditional DPIV measurements with a new class of Lagrangian analysis tools that analyze the wake as a distinct fluid structure with its own added-mass dynamics. The fin of the bluegill sunfish was observed to be embedded within this vortex wake structure. It is expected that a similar phenomenon will observed for other animal appendages moving in air or water. This result suggests that the dynamical effect of the vortex wake on locomotion is to replace the real animal fin with an `effective appendage,' whose geometry is dictated by the forward- and backward-time Lagrangian Coherent Structures (LCS) and whose interaction with the surrounding fluid is wholly governed by inviscid concepts from potential flow theory (Peng and Dabiri, 2007a; Peng and Dabiri, 2007b). A similar concept, the `displacement thickness', is well known in the theory of steady aero- and hydrodynamic flows (e.g. Rosenhead, 1963).

The benefit of this perspective based on an `effective appendage' is that it facilitates a direct correlation between the morphology and kinematics of the real appendage and the dynamics of locomotion. The expression for the locomotive force (Eqn 4) indicates that there are three parameters contributing to the force generation: the vortex added-mass coefficient C, volume V_V, and velocity U_V of the wake vortex body. Of these three parameters, C and V_V are more directly related to the geometry of the vortex while U_V is more closely related to the fin motion. Since locomotive forces are generated proportionally to the change of momentum [(1+C)V_VU_V] in unit time, the contribution of each parameter can be evaluated by comparing the three logarithmic terms on the right-hand side of the following equation:

(9)

Fig. 13 shows the change of the logarithm of each individual parameter and the total change in wake momentum. It can be seen that early in the stroke cycle the changes in the shape of the vortex are small and the total change in momentum follows closely the change in vortex velocity. As mentioned previously, the vortex velocity closely follows the fin motion during the stage of the stroke cycle. Therefore, force production during the early fin motion is primarily dependent on the fin kinematics. After the initial fin motion, force generation is primarily dictated by changes in the vortex shape; hence, fin morphology (which dictates vortex shape) governs force production later in the fin stroke cycle. Analyses of the vertical and horizontal forces are consistent, showing that the shape of the vortex remains relatively constant early in the stroke cycle but changes in the latter stages. Hence we conclude that locomotive forces can be generated not only by accelerating or decelerating the fin, but also by changing the shape of the vortex via fin morphology. The clarity of the parameter relationships in Fig. 13 is a benefit of the `effective appendage' approach taken presently. Since the timing of forces on the `effective appendage' is coincident with those on the real appendage (i.e. instantaneous forces are being evaluated), this perspective is compatible with investigation of muscle dynamics and activity patterns (Lauder et al., 2007; Mittal et al., 2007).

The ability of Lagrangian methods as implemented here to identify the boundary of this `effective appendage' in a fully unsteady flow without direct appeal to the vorticity field makes it possible to develop new models of fluid dynamic locomotion that are simultaneously more accurate and less complex than existing ones. Furthermore, the current heavy reliance on vorticity and other Eulerian field concepts, which are known to be inadequate for vortex identification in unsteady flows (Haller, 2005), is made unnecessary. The present results clearly illustrate that, as suggested by Haller (Haller, 2005), the correlation between the spatial distribution of vorticity and the actual vortex boundary can be quite poor in an unsteady flow, especially where global rotational motion is present. The downstroke of the bluegill sunfish pectoral fin is dominated by such rotational kinematics.

In this study, locomotive force is quantified only on the downstroke and early portions of the stroke reversal, when the vortex structure is primarily attached to the fin. The forces estimated are consistent with the magnitude of the time-averaged forces calculated using the same vorticity method as in Drucker and Lauder (Drucker and Lauder, 1999). The time-averaged vertical and lateral forces are 2.23 mN and 4.55 mN for this study (Fig. 12), which are consistent with the time-averaged forces quantified in Drucker and Lauder (Drucker and Lauder, 1999) for the same species (3.24 mN vertical and 6.96 mN lateral), noting that the data used in that study were taken from a different animal sample and using a different experimental arrangement. The evolution of forces is consistent with the motion of the animal. The lateral forces from the paired pectoral fins are expected to approximately cancel, resulting in the absence of substantial lateral motion as was observed empirically. The resultant vertical force, i.e. lift minus weight, is positive during the early phase but negative during the later phase, consistent with the kinematic measurements of Gibb et al. (Gibb et al., 1994). Since Gibb et al. only present kinematic descriptions of the locomotion (i.e. no force measurements) (Gibb et al., 1994), the comparison that can be made with the present work is only qualitative. Nonetheless, the force peaks in Fig. 12 are suggestive of a rise during the early phase of the fin beat and sinking during the late phase, as seen by Gibb et al. (Gibb et al., 1994). Perhaps coincidentally, the evolution of the locomotive forces resembles the pattern of forces generated by an insect wing, which has an early peak when the wing starts from rest followed by decay to a stable level (Birch and Dickinson, 2003).

The greatest challenge in this study was to identify the 3-D boundary of the vortex and its motion given 2-D planar DPIV data. As seen from Eqn 4, only by identifying the 3-D boundary of the vortex can the volume, added-mass and velocity of the vortex body be determined to evaluate locomotive forces. Since only 2-D velocity data are available in the present study, approximations were required.

A concern in the calculations is the assumption that the vortex ring is symmetrical about the transverse laser plane. This may not be true during the entire stroke cycle, and an asymmetrical distribution of the vortex ring on either side of the laser plane would cause underestimation of the vortex volume. Since the vortex cannot move far from the laser plane as long as it is attached to the fin, which is oscillating near the laser plane, the approximation used in this study of the attached vortex wake is reasonable. However, after the vortex is shed from the fin, it is advected out of the laser plane by the ambient current, making the approximation invalid. This is likely the effect observed toward the end of the present measurements (e.g. Fig. 11A,C).

Attempts were made to increase the accuracy of the 3-D vortex boundary approximation. FTLE fields were calculated on horizontal planes from DPIV data, with the aim of identifying the vortex boundary corresponding to the transverse plane measurement and thereby constructing the 3-D vortex boundary. However, boundaries were not clearly revealed in the corresponding FTLE fields due to the effect of net ambient flow on these two planes, which dominated the vortical motions. Even if it is possible to identify the boundary of the vortex on three perpendicular planes, approximation is still required to construct the 3-D boundary of the vortex (though the result may exhibit higher accuracy than the present methods when measuring more complex wake geometries). The identification of 3-D vortex boundaries would also enable evaluation of the effect of rotational added-mass, which is neglected in the present study. The determination of the 3-D vortex boundary requires the development of flow visualization techniques or numerical methods that can provide 3-D velocity field information, which can then be analyzed using the LCS method implemented in this study.

Given volumetric DPIV data for swimming or flying animals (or corresponding data computed in numerical simulations), a true validation of non-invasive force measurements demands a comparison of the animal body trajectory predicted by the force measurements with the body trajectory measured empirically. To our knowledge, such a comparison of measured and predicted kinematic data has not yet been performed by any study. Ideally, wake vortex dynamics and body kinematics should be measured simultaneously and in three dimensions. A demonstration that the estimated forces agree with the time-averaged force required to sustain lift of neutral buoyancy is necessary, but not sufficient by itself to validate instantaneous force estimates (cf. Dabiri, 2005). An additional avenue for validation would be comparison of calculated force profiles with those resulting from computational fluid dynamic analysis of the same fin beats using measured 3-D kinematics from that beat (with coupled fluid-structure interactions included in the computation for the swimming case, due to the comparable density of the fluid medium and the appendage). Once the force model has been validated, other aspects of animal behavior such as energetics can be examined according to established models (Schultz and Webb, 2002).

Finally, the relative importance of unsteady effects in this study was deduced a posteriori via direct examination of the contribution from wake vortex added-mass. Alternatively, the wake vortex ratio Wa introduced in Dabiri (Dabiri, 2005) and refined in Dabiri et al. (Dabiri et al., 2006) provides an a priori indication of the need to examine a particular flow with the level of scrutiny applied here.

	List of symbols

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 

A

area

C

added-mass coefficient for linear motion along the direction of U_V

C

added-mass tensor of the vortex

d_p

particle diameter

DPIV

digital particle image velocimetry

F_L

locomotive force

f_S

stroke frequency

FTLE

finite-time Lyapunov exponent

I

impulse

L

length

LCS

Lagrangian Coherent Structures

t

time

T_S

stroke duration

U_V

velocity of the wake vortex center of mass

V_V

volume of the vortex

w

width

distance

Lagrangian flow map

water density

uncertainty

fluid viscosity

	Acknowledgments

 
The authors thank Shawn Shadden, Jerry Marsden and Parviz Moin for discussions of the LCS theory. This work was supported by NSF grant OCE0623475 to J.O.D.; an ONR-MURI Grant N00014-03-1-0897 on fish pectoral fin function, monitored by Dr Thomas McKenna and initiated by Dr Promode Bandyopadhyay; and NSF grant IBN0316675 to G.V.L.

	References

TOP Summary Introduction Materials and methods Results Discussion List of symbols References

 

Adrian, R. J. (1995). Limiting resolution of particle image velocimetry for turbulent flow. Proc. 2nd TRA Conference (Pohang, S. Korea). Adv. Turbulence Res.1 -19.

Alben, S., Madden, P. G. A. and Lauder, G. V. (2007). The mechanics of active fin-shape control in ray-finned fishes. J. R. Soc. Interface doi:10.1098/rsif.2006.0181 . In press.

Birch, J. M. and Dickinson, M. H. (2003). The influence of wing–wake interactions on the production of aerodynamic forces in flapping flight. J. Exp. Biol. 206,2257 -2272.[Abstract/Free Full Text]

Dabiri, J. O. (2005). On the estimation of swimming and flying forces from wake measurements. J. Exp. Biol. 208,3519 -3532.[Abstract/Free Full Text]

Dabiri, J. O. (2006). Note on the induced Lagrangian drift and added-mass of a vortex. J. Fluid Mech. 547,105 -113.[CrossRef]

Dabiri, J. O. and Gharib, M. (2004). Fluid entrainment by isolated vortex rings. J. Fluid Mech. 511,311 -331.[CrossRef]

Dabiri, J. O., Colin, S. P., Costello, J. H. and Gharib, M. (2005). Flow patterns generated by oblate medusan jellyfish: field measurements and laboratory analyses. J. Exp. Biol. 208,1257 -1265.[Abstract/Free Full Text]

Dabiri, J. O., Colin, S. P. and Costello, J. H. (2006). Fast-swimming jellyfish exploit velar kinematics to form an optimal vortex wake. J. Exp. Biol. 209,2025 -2033.[Abstract/Free Full Text]

Dickinson, M. H. (1996). Unsteady mechanisms of force generation in aquatic and aerial locomotion. Am. Zool. 36,537 -554.

Dickinson, M. H. and Götz, K. G. (1996). The wake dynamics and flight forces of the fruit fly Drosophila melanogaster. J. Exp. Biol. 199,2085 -2104.[Abstract]

Drucker, E. G. and Lauder, G. V. (1999). Locomotor forces on a swimming fish: three-dimensional vortex wake dynamics quantified using digital particle image velocimetry. J. Exp. Biol. 202,2393 -2412.[Abstract]

Drucker, E. G. and Lauder, G. V. (2000). A hydrodynamic analysis of fish swimming speed: wake structure and locomotor force in slow and fast labriform swimmers. J. Exp. Biol. 203,2379 -2393.[Abstract]

Drucker, E. G. and Lauder, G. V. (2001). Wake dynamics and fluid forces of turning maneuvers in sunfish. J. Exp. Biol. 204,431 -442.[Abstract]

Drucker, E. G. and Lauder, G. V. (2003). Function of pectoral fins in rainbow trout: behavioral repertoire and hydrodynamic forces. J. Exp. Biol. 206,813 -826.[Abstract/Free Full Text]

Gibb, A., Jayne, B. C. and Lauder, G. V. (1994). Kinematics of pectoral fin locomotion in the bluegill sunfish Lepomis macrochirus. J. Exp. Biol. 189,133 -161.[Abstract]

Haller, G. (2000). Finding finite-time invariant manifolds in two-dimens