Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A computational investigation of lift generation and power expenditure of Pratt’s roundleaf bat (Hipposideros pratti) in forward flight

  • Peter Windes,

    Roles Conceptualization, Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia, United States of America

  • Xiaozhou Fan,

    Roles Data curation

    Affiliation Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia, United States of America

  • Matt Bender,

    Roles Resources, Writing – review & editing

    Affiliation Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia, United States of America

  • Danesh K. Tafti ,

    Roles Conceptualization, Supervision, Writing – review & editing

    dtafti@exchange.vt.edu

    Affiliation Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia, United States of America

  • Rolf Müller

    Roles Resources, Supervision, Writing – review & editing

    Affiliations Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia, United States of America, Shandong University Virginia Tech International Laboratory, Jinan, China

Abstract

The aerodynamic mechanisms of bat flight have been studied using a numerical approach. Kinematic data acquired using a high resolution motion capture system was employed to simulate the unsteady air flow around a bat’s wings. A flapping bat wing contains many degrees of freedom, which make 3D motion tracking challenging. In order to overcome this challenge, an optical motion capture system of 21 cameras was used to reduce wing self-occlusion. Over the course of a meter-long flight, 108 discrete marker points on the bat’s wings (Pratt’s roundleaf bat, Hipposideros pratti) were tracked. The time evolution of the surface of each wing was computationally reconstructed in 3D space. The resulting kinematic model was interfaced with an unsteady incompressible flow solver using the immersed boundary method (IBM) and large eddy simulation (LES). Verification and validation of the flow simulation were conducted to establish accuracy. The aerodynamic forces calculated from the simulation compared well to the forces theoretically needed to sustain the observed flight trajectory. The transient flow field generated by the simulation allowed for the direct calculation of lift, drag, and power output of the bat during flight. The mean lift coefficient was found to be 3.21, and the flap cycle averaged aerodynamic power output was 1.05 W. Throughout the flap cycle, the planform area of the wings varied up to 46% between the largest and smallest values. During the upstroke, wing rotation was found to mitigate negative lift thereby improving overall flight efficiency. The high resolution motion capture and flow simulation framework presented here has the potential to facilitate the understanding of complex bat flight aerodynamics for both straight and maneuvering flight modes.

Introduction

Recently, there has been growing interest in understanding and optimizing low Reynolds number flight in the regime, . For example, micro air vehicles (MAVs) have been proposed as sensing and monitoring platforms for agriculture, security, and other applications [1]. However, the aerodynamics of small flight vehicles is fundamentally different from large aircraft which operate in the high Reynolds number regime, Re ~ (106) to (109). Looking to nature, biological flapping fliers such as insects, birds, and bats can lend insight towards the design of MAVs [2,3].

One hallmark of flapping flight is that it remains effective at low Reynolds numbers, where the efficiency of fixed wing systems begin to decline. Flow separation over an airfoil at high angles of attack leads to the stall condition, resulting in a large drop in the lift to drag ratio [2]. In contrast, flapping flight allows for much more versatility in controlling the behavior of vortices in the flow [1,4,5]. This translates to improved flight performance over a range of Reynolds numbers.

Bats are unique from other biological fliers due to their highly articulated wing structure, and their ability to actively contour their wing membrane during flight [68]. Compared to birds and insects, bats exhibit greater changes in wing surface area during flight, and their wings have more degrees of freedom [9]. These features are evidenced by their ability to weave through dense forests while flying, as well as chase and capture insects [10].

The study of bat flight over the last several decades has evolved through several phases, each limited by the technology available to researchers at the time. Early research focused on morphological measurements derived from specimens, and documenting physiological parameters related to flight performance [6,1113]. Additionally, video or multiple exposure still images were used to study basic wing kinematics during flight [1421]. The measured flight kinematics allowed researchers to model aerodynamic forces using a steady-state flow approximation [14].

Between 2005 and 2010, improvement in high-speed digital videography allowed researchers to employ 3D stereoscopic motion capture systems to measure bat wing kinematics in wind tunnels [22,23]. Around the same time, particle image velocimetry (PIV) was employed to measure air flow in the wake of flying bats [9,2429]. Estimates of lift and power output were made from measurements of the circulation and kinetic energy in the wakes [8,30]. This was a significant improvement over previous steady-state models, and allowed for the investigation of unsteady flow phenomena.

However, recent research has raised questions regarding the accuracy of wake based methods by suggesting that some of the underlying assumptions may not hold [31,32]. The frozen turbulence hypothesis—a prerequisite for application of the Kutta-Joukowski theorem, vortex ring model, and the actuator disk model—posits that vortices remain intact as they advect downstream into the wake, which is rarely the case in animal flight [31]. In their 2015 review on bat flight, Hedenström et al. noted that numerical simulation, when they become computationally feasible, have the potential to advance the state of the art and augment existing bat flight research [8].

Numerical analysis is adding a new dimension to the study of unsteady aerodynamic flows. Recently, significant progress has been made towards the simulation of animal flight aerodynamics—for example, dragonflies [33,34], cicadas [35,36], and hummingbirds [37,38]. One feature all these animals have in common is their relatively stiff wings in comparison to bats. Fewer degrees of freedom greatly simplifies both data collection as well as the flow simulations. In contrast, progress towards the simulation of bat flight is lacking in the literature. In 2014, Viswanath et al. conducted a first-of-its-kind aerodynamic simulation of the left wing of a fruit bat in climbing flight using the immersed boundary method [39]. The time-varying position of 49 discrete points on the bat wing drove the wing kinematics in the model. In 2015, Wang et al. conducted a numerical simulation of bat flight, however the wing kinematics were only defined by 10 discrete points. That is, the wing motion was dictated by the position of the shoulders, legs, and three points per wing. Based on estimates by Riskin et al., the skeleton alone within a bat wing contains over 20 degrees of freedom [23]. When membrane deformation is included, a minimum of 30–40 discrete points are required to fully capture the complexity of bat wing kinematics.

The goal of the present study is to examine the underlying aerodynamic mechanisms, which make bats effective fliers in the low Reynolds number regime. In the process, we aim to establish and validate a framework for capturing bat flights in 3D, and investigate it using fluid simulations. In order to achieve these goals, we demonstrate two novel advances in the study of bat flight—the use of multiple video cameras for high spatial resolution stereoscopic wing tracking, and the use of numerical flow simulations to allow for detailed examination of unsteady flight aerodynamics. These are some of the first high-resolution studies of bat flight, which give an unprecedented window into the physics of bat flight by direct determination of the time-dependent forces generated at the wing surface unlike experimental studies which attempt to deduce the mean lift forces generated by measurements in the wake.

Motion capture of bat flight

There are several inherent challenges to conducting aerodynamic simulations of bat flight. Unlike insect wings which primarily flex, bat wings are comprised of a highly articulated skeleton which is spanned by an anisotropic, active wing membrane. These significant biological differences allow the bat wing to deform substantially during even the simplest of flight regimes. This has two major implications. First, the wings self-occlude during the flap stroke blocking aspects of the kinematics from observation. This can lead to voids in the final 3D kinematic data set if not properly addressed. Second, no simple parametric model can describe the surface of the bat wing, making the mathematical characterization of the wing shape non-trivial. For these two reasons, other methods which have been recently employed to motion capture wing kinematics of other biological fliers such as parrotlets [40], humming birds [33,41], dragonflies [34], and bees [42] are unsuitable for capturing bat flight.

In order to overcome the aforementioned challenges, a highly redundant camera array was constructed which greatly ameliorated the wing self-occlusion problem. Flight kinematic data was captured in a 1.2 m × 1.2 m cross section flight tunnel fitted with 21 cameras arranged in 3 rings about 40cm apart. The large camera array allowed for an expansive motion capture volume in addition to reducing the number of possible wing occlusions. As demonstrated in the present work, this system is capable of collecting wing kinematic data over an extended flight relative to other existing methods. Furthermore, the motion captured bat flights are not restricted in speed or direction as they are in wind tunnels, allowing for the study of maneuvering flights using the same setup.

Aerodynamic simulation of bat flight

Another challenge is the reconstruction of the bat wing surface from the 3D kinematic data. The interplay between passive and active control of the bat wing membrane is currently not known. Thus, kinematic data collection methods which only track sparse features on the bat wing and reconstruct the surface using a priori assumptions about the behavior of the membrane risk mischaracterization of the wing shape. This challenge has been addressed in the present work by increasing the spatial resolution of the wing motion capture relative to other similar studies. Using a high marker density across the bat wings allowed for the direct measurement of all aspects of the wing kinematics—both skeletal joint motion as well as membrane stretching.

The present work uses the immersed boundary method to represent the bat wing in the fluid flow simulation. As other researchers have pointed out, this is an effective strategy for simulating flapping flight as it avoids the need for complex and computationally expensive re-meshing techniques [43]. Many bats, such as the one studied in the present work, fly at a Reynolds number of order 104. This is a challenging intermediate Reynolds number flow to simulate since viscous effects are not negligible, requiring the solution of the full Navier-Stokes equations. Additionally, the wake does not remain strictly laminar. A parallel numerical computation on approximately 200 to 300 CPU cores was required to resolve the flow in one to four days depending on the grid resolution.

Methods

Flight data collection

The motion capture system consisted of a large array of action cameras arranged inside a rectangular tunnel, as shown in Fig 1. The approximate dimensions of the tunnel were 1.2m × 1.2m × 4m. The camera model used was the GoPro Hero 3+ Black, with the resolution set to 1280×720 (720p) and the frame rate set to 120 fps. This setup allowed for the capture of up to 4 meters of bat flight. The cameras were arranged on the walls of the rectangular flight tunnel in order to observe the flying bat from various viewing directions and to minimize self-occlusion of the wings.

thumbnail
Fig 1. Motion capture system, consisting of 21 GoPro cameras inside a 1.2 x 1.2 meter cross section flight tunnel.

The bat was trained to fly through the tunnel and land on the cork perch shown.

https://doi.org/10.1371/journal.pone.0207613.g001

The cameras were synchronized using a central control unit such that all video frames were captured simultaneously (±0.1 ms) on all 21 cameras. The tunnel was illuminated using six lights mounted in the upper corners to mitigate motion blur while filming at high frame rates. The motion tracking system used in the present work is described in more detail in Bender et al., 2015 [44] and Bender et al., 2016 [45].

In the present study, bat wing kinematic data was obtained from an adult male Pratt"s roundleaf bat (Hipposideros pratti), an insectivorous bat indigenous to China and neighboring countries in southeast Asia. The bat was collected from a cave in southern China. Data acquisition was conducted at the Shandong University–Virginia Tech International Laboratory in Jinan, China, with oversight from Virginia Tech’s Institutional Animal Care and Use Committee (IACUC) under protocol number 15–067.

When no experiments were being conducted, the bat was kept in a small group of conspecifics in an indoor flight room (1.3 m wide, 6 m long, and 3 m high). In order to keep the bats active during the daytime experiments, the daytime and nighttime of bats were switched by lighting the bat room at nighttime for 10 hours and darkening the bat room at daytime for 14 hours. The bats were fed a daily diet of mealworms with vitamins and mineral supplements and were provided water ad libitum. Measurements of the particular specimen that participated in the experiments are given in Table 1.

Multiple small white “marker points” were placed on each wing as landmarks for tracking, as shown in Fig 2. The spatial resolution of the reconstructed kinematic data was determined by how many markers were used. Using too few points to measure bat flight inadequately captures all the degrees of freedom of the bat wing. In the present work, 108 points were tracked over the course of several wing beat cycles. The final wing kinematic data set consisted of 108 discrete points at 52 instances in time.

thumbnail
Fig 2. Small white markers were placed on each bat wing to aid 3D reconstruction of the wing kinematics during flight.

https://doi.org/10.1371/journal.pone.0207613.g002

The cameras were calibrated using the Svoboda Multi-Camera Self-Calibration toolbox for Matlab [46]. The calibration solves for the position, orientation, and radial distortion of each of the cameras in the array such that the 3D position of any point in the scene can be reconstructed from two camera views. If more than two cameras observe the same point, multiple reconstructions of that point are generated corresponding to each camera pair. The final 3D representation was taken to be the median of the point cluster. Once the data was digitized for the entire flight, principal component analysis (PCA) was used to filter out erroneous points in a data cleaning process which is described in Fan et al., 2018 [47]. Finally, a time series of coordinates of each of the marker points on the bat wing was obtained. Fig 3 shows five snapshots from the flight video along with the corresponding 3D reconstruction.

thumbnail
Fig 3. Five snapshots of the bat in flight along with the corresponding 3D reconstructed data.

The radial distortion of the camera lens slightly changes the visual perspective of the bat in the video, however that is mathematically accounted for during the camera calibration process. In the digital reconstruction, the magenta trace is the path of the bat body, and the green traces are the path of the wing tips.

https://doi.org/10.1371/journal.pone.0207613.g003

Flight simulation methodology

The aerodynamic simulations were conducted using an in-house incompressible flow solver, GenIDLEST (Generalized Incompressible Direct and Large Eddy Simulation of Turbulence) [48,49]. The flow of air over bat wings is governed by the incompressible Navier-Stokes equations, (1) where is the Cartesian velocity vector, p is pressure, t is time, ρ is air density, and v is kinematic viscosity. In GenIDLEST, solution of the Navier-Stokes equations is achieved using a finite volume framework with a second-order central difference discretization scheme for the convective and viscous terms. The pressure correction method is used with preconditioned BiCGSTAB linear solvers. Message Passing Interface (MPI) is employed for parallelization, and up to 640 CPU cores were used in the present work. Turbulence modeling was done using Large Eddy Simulation (LES). Specifically, the dynamic Smagorinsky subgrid stress model was used, where the model coefficient is calculated locally at each grid point as the simulation progresses [50,51]. The adaptive nature of the subgrid model is important due to the transitional Reynolds number of the flow (Re = 11,680).

A sharp interface immersed boundary method (IBM) was used to impose the wing surface boundary in the flow (Fig 4). The Navier-Stokes equations were solved on a structured Cartesian background grid, while the bat wing is represented by an unstructured triangular surface grid. Details of the particular implementation and validation of this IBM algorithm can be found in Nagendra et al., 2014 [52]. A thin surface approximation was used for the immersed surface, thus the net force on the wing was taken as the vector sum of the force on each side of a given surface element.

thumbnail
Fig 4. Immersed boundary method setup.

A thin interface (red) is embedded into the background fluid grid. The no slip boundary condition is enforced on the interface by applying an appropriate velocity and pressure values to the IB nodes (green) at each time step. The property values at each IB node are determined by creating velocity profiles between the fluid nodes and the no slip surface.

https://doi.org/10.1371/journal.pone.0207613.g004

Numerical implementation of flapping wings

In order to represent the experimentally generated wing kinematics using the immersed boundary method (IBM), a procedure was developed to interpolate the wing surface spatially and temporally. The 3D reconstructions of the observed wing marker points are referred to in this section as “control points.” First, the control point cloud was triangulated using the Delaunay method, resulting in a coarse triangular surface mesh. This surface mesh, however, was comprised of elements with edge lengths of approximately 2 cm—too coarse for IBM simulations. Thus, the planar inner region of each control triangle was sub-meshed resulting in a finer mesh with each element edge length on the order of the background fluid mesh edge length. Finally, the wing mesh consisted of 182 control triangles, and 42,839 fine mesh elements depicted in Fig 5. The coarse mesh was used to impart the proper kinematic motion to the immersed surface, while the interpolated fine mesh was used to specify the wing surface location on the background fluid grid.

thumbnail
Fig 5. Wing surface mesh viewed from above during the most outstretched point of the downstroke.

The coarse and fine wing surface meshes are showed overlaid with the fluid grid in the background. The vertices of the coarse triangular mesh correspond to the white marker points on the bat wings.

https://doi.org/10.1371/journal.pone.0207613.g005

The location of each fine mesh vertex within a given control triangle can be described as a linear combination of two parameters—α and β—ranging between 0 and 1. (2) where is the interior point to be interpolated, and , and are the three surrounding control point vertices, as shown in Fig 6. The fine mesh is generated based on the initial position and shape of the wings. The alpha and beta values are calculated once from the initial wing state and remain constant as the wing moves, causing the interior points to maintain their relative positions throughout the flapping cycles.

thumbnail
Fig 6. A given control triangle consists of three control points, and several interior points.

The area inside the triangle is parameterized into a 2D surface, where every location is described by a unique pair of α and β values. As the wing flaps, the triangles deform and the locations of the interiors points are updated accordingly.

https://doi.org/10.1371/journal.pone.0207613.g006

Movement of the immersed surface wing mesh is achieved using piecewise spline equations which describe the continuous evolution of each control point in time. At each time step, the spline equations are evaluated at the given simulation time. Subsequently, all fine mesh points are updated based on their respective stored α and β values as well as the updated control point locations.

Once the wing surface mesh reaches its new location at the next time instance, the centroid of each surface element is located on the computational grid using a search algorithm—that is, the element’s i,j,k computational coordinates are determined using its Cartesian x,y,z physical coordinates. This allows for proper data exchange between the fluid and the solid surface to ensure that the presence of the solid surface is manifested in the fluid field. The no-slip boundary condition is applied at the surface using the immersed boundary method. The wing surface velocity is calculated for each element using a second order backward difference scheme. Similarly, the wing acceleration is calculated using the second time derivative of the surface element locations.

Aerodynamic simulation setup

The size of the fluid domain was set 32cm × 16cm × 16cm along the x-, y-, and z-directions respectively. The dimensions of the cross-section of the domain approximately match the physical size of the flight tunnel. The fluid domain represents a moving reference frame that follows the bat’s flight in the x-direction at 2.5 m/s. Thus, 2.5 m/s was subtracted from the absolute x-velocity of the measured wing kinematic data. The displacement of the bat body within the fluid domain during the course of the simulation represented a perturbation relative to the moving reference frame.

The inlet plane and four side-wall faces of the domain were set to a constant velocity boundary condition that was equal to the reference frame velocity of 2.5 m/s. The outlet plane (+x) of the domain was set to a zero gradient boundary condition, allowing air to flow out of the domain. The bat was initially located at the origin, and the domain extends 8 chord lengths upstream of the bat and 24 chord lengths downstream of the bat as shown in Fig 7. During the course of the simulation, no fluid disturbances were observed to propagate through either the inlet or outlet boundary planes. Thus, the 32cm length of the domain in the streamwise direction was deemed sufficient to prevent unphysical interference by the boundary.

thumbnail
Fig 7. The full fluid grid viewed from above shows the refined region in the vicinity of the bat.

The y- and z-faces of the domain reflect the physical size of the experimental flight tunnel, while the x-faces were positioned sufficiently far from the bat to avoid interference from the boundary.

https://doi.org/10.1371/journal.pone.0207613.g007

The fluid domain was discretized using a non-uniform, orthogonal, Cartesian mesh for optimal computational performance. The grid contained 88.1 million fluid cells and the refined region had an average grid spacing of, Δfluid = 0.020 chord lengths. Thus, on average there were 50 fluid cells per chord length.

The mesh for the immersed boundary method was an unstructured triangular surface mesh which contained 21,354 vertices, and 54,238 triangular elements. The initial average edge length was Δsurface = 0.022. As the wing flaps, the area changes and the relative position of the vertices is maintained using the movement algorithm described in Section 2.3.

Grid independence verification

In order to validate the computational grid used in the present work, a grid independence study was conducted. Five grid resolutions were compared to determine a sufficient level of refinement. The number of fluid cells along the wing chord was used to define the grid levels. The five grids used, ranging from 12 million to 210 million fluid cells, are outlined in Table 2.

thumbnail
Table 2. Several different grids ranging from 20 fluid cell lengths per wing chord length to 100 were compared to ensure grid independence.

The CPU architecture used was Intel Xeon E5-2670 (Sandy Bridge).

https://doi.org/10.1371/journal.pone.0207613.t003

The fluid grid was refined in the region in which the wing surface passes, and coarsened near the edges of the domain as show in Fig 8. In doing so, the entire swept volume of the wings falls within the refined region, properly resolving the unsteady flow.

thumbnail
Fig 8. The grid spacing distribution for the five grids described in Table 2.

https://doi.org/10.1371/journal.pone.0207613.g008

Simulations of the bat flight were run on each of the five grids. The difference in cycle averaged force magnitude between each case and the finest grid is shown in Table 2. There is a noticeable difference in force of 5.9% between the 20 and 100 cells per chord grids. However, the difference drops to 1.6% when comparing the forces from the 50 and 100 cells per chord grids which indicates close agreement between the grids. Given the uncertainties in the camera calibration and stereo triangulation, the small differences between the two finest grids can be considered negligible.

In order to further ensure grid convergence, the transient surface integrated fluid forces on the wings were compared for the 20, 50, and 100 cells per wing chord cases as shown in Fig 9. The y-dir force is close between all three grids. The x-dir force is close aside from the point of maximum thrust (-x force) between 120 ms to 160 ms. Within this time frame, the 20 cells/chord grid slightly under predicts thrust, however the 50 and 100 cells/chord grids compare closely. The z-dir force compares closely on the upstroke (180 ms to 250 ms), but differences can be observed on the downstroke around the time of peak lift production (120 ms to 180 ms). The 20 cells/chord grid consistently under predicts lift force through most of the downstroke. The 50 and 100 cells/chord grids compare well with only a few small differences, which do not materially impact the results.

thumbnail
Fig 9. Comparison of the unsteady fluid force components for a complete flap cycle.

Some discrepancy can be observed between the coarsest (blue) and finest (black) grids, however the two finest grids—50 and 100 cells per chord—show close agreement throughout the flap cycle for all three force components.

https://doi.org/10.1371/journal.pone.0207613.g009

Results

Kinematic analysis

In the present study, a 1 m long flight segment was chosen for analysis. The trajectory was mostly straight and level, however the bat was slightly descending at a 5° angle and accelerating from 2.2 m/s to 2.7 m/s. This flight spanned 50 frames of video with a total duration of 0.42 s corresponding to about 3 flap cycles. Several snapshots from the flight are shown in Fig 10.

thumbnail
Fig 10. Nine frames sampled from two complete flap cycles are shown.

t/T is the normalized time starting from the downstroke, and normalized by the flap period of 138 ms.

https://doi.org/10.1371/journal.pone.0207613.g010

The bat had a wingspan of 52.0 cm, and an average wing chord of 7.4 cm. Thus the Reynolds number based on the mean wing chord was, Re = 11,680, and is defined as, (3) where ρ is the air density, U is the average flight velocity, cm is the chord length, and μ is the air viscosity. The flow is dominated by inertia, however the Reynolds number is low enough that viscous effects cannot be neglected. The Strouhal number of the flight was, St = 0.61, defined as, (4) where f is the flapping frequency, 7.25 Hz, ha is the stroke or plunge amplitude, 10.9 cm, and U is the average flight velocity, 2.57 m/s. All the flight parameters are summarized in Table 3.

thumbnail
Table 3. Flight parameters from a 1 meter long straight flight by an H. pratti.

https://doi.org/10.1371/journal.pone.0207613.t004

The plunge amplitude is defined by half the maximum minus minimum wing tip z-positions. This value varies between the left and right wing, and between flaps so the mean was taken. The stroke plane angle is defined by the angle of the stroke plane with respect to horizontal, and describes the balance between lift and thrust. Hovering flight requires low stroke plane angles which correspond to nearly horizontal stroke planes. Conversely, accelerating flight requires a high stroke plane angle to achieve greater thrust. The value of 54.1 degrees for the present straight and level flight is typical of bats in this flight regime [53].

Variations in wing surface area during flight is another aspect of the wing kinematics which can only be studied using a sufficiently spatial resolution of the wings. Fig 11 shows the time evolution of the total wing surface area along with a breakdown between different regions of the bat. The mean total surface was 398 cm2, and the variation was +17% and -35% for the minimum and maximum, respectively. The minimum total surface area was observed during the latter third of the upstroke, while the maximum was observed near the middle of each downstroke. The wing area begins to increase during the second half of the upstroke. This is possibly due to either active muscular control of the wing, higher air pressure under the wing as it rotates at the top of the upstroke thereby catching the flow, or a combination of both factors.

thumbnail
Fig 11. The total surface area variation of the bat wings, tail, and body is shown over the course of the flight (solid black).

For context, the wing tip positions are shown (dotted black) to indicate the upstroke and downstroke.

https://doi.org/10.1371/journal.pone.0207613.g011

The trend observed in the wing surface area change is consistent for the left and right wing, as well as the inner and outer portions of each wing. Thus, no section of the wing was observed to experience vastly different percent change in area compared to the wing as a whole.

Validation using lumped mass dynamics

In order to verify the flight simulation, the fluid forces exerted on the bat were examined in relation to the observed flight trajectory. The transient force output from the flow simulation was used to predict the flight trajectory of an equivalent point mass. The bat’s body mass of 55g was approximated as being concentrated between the shoulder points at the same location that the net force was applied. This allowed for a comparison that would highlight the accuracy of the numerical solution.

The transient force curve (Fig 9) was calculated by integrating the pressure and shear force over the surface of the wing at each time step. As expected, the dominate force component was lift countering the weight of the bat. For a perfectly straight and level flight, the streamwise and lateral components should be zero. Absent acceleration, the thrust and drag should exactly balance for a net zero x-force. For the present flight, the net x and y forces were indeed near zero as seen in Fig 9.

The lift, thrust, drag, lateral, and weight forces were applied to the bat. The expected acceleration of the bat due to the fluid forces at each time step was calculated by dividing the net fluid force by the mass of the bat.

(5)(6)

A Runge-Kutta numerical integration was performed on Eq (6) to predict the flight trajectory based on the transient lift and propulsion curves generated by the flow simulation. The trajectory, , predicted from the fluid forces was then compared to the observed flight trajectory of the bat as shown in Fig 12.

thumbnail
Fig 12.

Left: Comparison of the simulated and observed velocity in x, y, and z. Right: The simulated flight trajectory is compared to the observed flight trajectory. The vertical and lateral predictions are very close, however the streamwise position is under predicted by 15 cm.

https://doi.org/10.1371/journal.pone.0207613.g012

The dynamics analysis approximates the mass of the bat to be a single rigid body, and ignores the wing inertia. For the assumptions made, the predicted trajectory is remarkably close to the observed flight path of the bat. The one notable discrepancy is an under prediction of the thrust force. Minimal acceleration is predicted, yet the bat accelerated from 2.2 to 2.7 m/s. The grid independence study indicated good convergence of the numerical simulation, thus the slight under-prediction of the x-location is more likely due to small errors in the wing kinematic data. On the other hand, the predicted y- and z- trajectories are very close to the observed values. Overall the dynamics analysis indicates that the forces calculated by the flow simulation are close to the actual fluid forces generated by the bat.

Force analysis

During flight, forces are manifested as a pressure differential across the wings. At the Reynolds number, Re ~ O(104), shear stress at the wing is relatively minor compared to the pressure force, however both are included in the present analysis. The total aerodynamic force on the wings at each time instant is the surface integral, (7) where p is fluid pressure, is the wing surface normal, and is the viscous stress tensor. The resulting force curve is shown in Fig 13.

thumbnail
Fig 13. The time variation of aerodynamic force is shown along with the wing tip position for context.

The peak lift force was around 1.2 N, and the cycle averaged mean was 0.525 N. The streamwise and lateral forces were both close to zero since the flight was approximately straight and level.

https://doi.org/10.1371/journal.pone.0207613.g013

The cycle averaged force values are , , and . Due to the low accelerations experienced during a mostly straight and level flight, the streamwise and lateral forces are expectedly near zero. The much larger z-component represents the lift necessary for flight. The lift can be normalized by standard methods to obtain the lift coefficient, (8) where Fz is the dimensional vertical force, ρ is the air density, U is the freestream velocity, and S is the maximum observed planform area of the wing during the flap cycle. The resulting values are plotted over two flap cycles in Fig 14.

thumbnail
Fig 14. The lift coefficient is plotted along with the wing tip location for context.

The mean value was CL = 3.21.

https://doi.org/10.1371/journal.pone.0207613.g014

The mean value of the lift coefficient was calculated to be 3.21. For fixed wing aircraft, the lift coefficient is only a function of airfoil shape and angle of attack. For this reason, it provides a convenient method to calculate lift at a given angle of attack, or calculate the proper angle of attack given a velocity. However, for bat flight, the wing shape, wing area, angle of attack, and instantaneous velocity vary continuously throughout the stroke. Others have noted the limitations of quasi-steady theory in the study of bat flight as well [8].

To understand the aerodynamic phenomena in flapping flight the unsteady flow features in relation to force production are examined. This has been the focus of recent research reviewed by Shyy et al. in 2010 [1]. Coherent vorticity is shown via Δ-criterion isosurfaces in Fig 15. Four instances in time are shown along with the corresponding flight velocity, planform area, instantaneous span, aerodynamic force, and aerodynamic power.

thumbnail
Fig 15.

Top: Coherent vorticity (iso-surfaces of Δ-criterion) is shown along with the wing surface pressure at four snapshots throughout the flap cycle. Bottom: Flight velocity, planform area, wing span, aerodynamic force, and aerodynamic power are plotted with each of the four snapshot locations indicated. (a) the top of the upstroke, (b) the point of maximum lift production, (c) the bottom of the downstroke, and (d) the midpoint of the upstroke.

https://doi.org/10.1371/journal.pone.0207613.g015

In Fig 15, snapshot (a) corresponds to the top of the upstroke or also the beginning of the downstroke. The lift force at (a) is about 0.45 N, or 35% of maximum lift, despite the wing having no downward velocity. Wing rotation towards the end of the upstroke as the wing readies itself for the downstroke results in a leading edge vortex (LEV) to begin to form. The aerodynamic power at time (a), as shown in Fig 15, is quite low at about 0.3 W, which is less than 10% of the maximum power expenditure. This is due to relatively low instantaneous wing flap velocity at the top of the upstroke. The planform area at time (a) has already increased to a near maximum indicating that either the freestream flow is impacting the underside of the wing stretching the membrane, or that the bat is actively stretching its wing with muscular activation. At time (a), the tip vortex (TiV) shedding is observed to be minimal. Between snapshots (a) and (b), the downward velocity of the wings accelerates rapidly. The majority of lift generation is in a 34 ms time interval from 258 ms to 292 ms which is highlighted in detail in Fig 16.

thumbnail
Fig 16.

Left: Aerodynamic pressure on the top and bottom of the bat wing is shown at 5 instances during the downstroke. Right: Pressure difference is shown on the wing surface with iso-surfaces of coherent vorticity.

https://doi.org/10.1371/journal.pone.0207613.g016

Fig 16 shows the coefficient of pressure on the bottom and top surfaces of the wing. In addition, the pressure difference—defined as the pressure under the wing minus the pressure on the top of the wing—is shown along with iso-surfaces of coherent vorticity. During both downstrokes, a notable pattern of three local maxima in the lift force can be observed (110–160 ms and 260–290 ms). Since this pattern emerges in different flap cycles, it is more likely a feature of the aerodynamics rather than random fluctuations in force. In order to explore the phenomenon in more detail, each of the five time instances shown in Fig 16 correspond to a local maxima or local minima in the force curve.

The first local maximum depicted in Fig 16, t = 258 ms, shows a coherent low pressure band on the leading edge of the wing. Moving along the chord from the leading to trailing edge, there is a distinct location, about a third of the way down the wing chord, at which the pressure rises, representing the bound of the influence of the LEV. At this instant in time, the stagnation pressure under the wing is maximum causing the large lift force—approximately 81% of the absolute maximum lift. The first local minimum in lift, t = 267 ms (Fig 16), shows the LEV convecting down the wing and breaking up slightly. Comparing to t = 258 ms, the low pressure band on the leading edge is not as complete. The stagnation pressure is also not as high at t = 267 ms compared to the previous local maximum at t = 258 ms. As such, the lift force is only 62% of the absolute maximum. At the second local maximum of lift (t = 275 ms in Fig 16), the low pressure region on the top of the wing grows in area as the LEV continues to convect along the wing. This in combination with the stagnation pressure rising slightly from t = 267 ms to t = 275 ms, causes the lift to rise to 90% of the absolute maximum. At t = 286 ms, the LEV becomes somewhat less coherent on the top surface and the low pressure region decreases in strength causing another local minimum if lift at 79% of the peak lift value.

Finally, the fifth and last instant shown in Fig 16, t = 292 ms, corresponds to the absolute maximum of lift which is also snapshot (b) from Fig 15. At this instant, a strong LEV can be seen which helps keep the flow remain attached despite a relatively high angle of attack. The LEV also maintains a strong low pressure region on the top of the wing directly contributing to lift. Flow stagnation under the wing creates an elevated pressure, which is especially pronounced towards the outer region of each wing. Additionally, prominent TiVs are shed which has also been observed in PIV experiments of bat flight [22,24,25,29,5457]. In contrast to fix wing high Reynolds number flight, in which TiVs are known to be detrimental, flapping fliers can use TiVs to augment lift and thrust due to interactions with the leading and trailing edge vortices [1]. The point of peak lift force also is the point of maximum power expenditure at about 3.6 watts. The vertical velocity component at this time is also rising. The cyclical nature of lift production in flapping flight means that the vertical position of the bat’s body moves slightly up and down during flight.

Spanshot (c) from Fig 15 is taken at the bottom of the downstroke. Lift production drops rapidly towards the end of the downstroke—from (b) to (c) there is a nearly 90% drop in lift. The prominent TiVs can be seen convecting downstream into the wake, yet have little impact on the flight once they are shed and are significantly distanced from the wing surface. The bat begins to rotates its wings aligning them to the flow and drawing the wing tips together. At time (c) the span is about 33 cm compared to the maximum span of 52 cm. Despite the wingtip velocity dropping to zero at the bottom of the downstroke, the power expenditure is still around 1 watt, likely due to wing inertia and elastic potential being converted to aerodynamic power as the wings stop and reverse direction. The vertical velocity of the bat body is near maximum at this point, and then begins to reverse.

Spanshot (d) from Fig 15 corresponds to the mid-upstroke. The upstroke is shorter in duration than the downstroke, which allows the bat to quickly reset the wing position and prepare for the next downstroke where the majority of lift and thrust are generated. During the upstroke, the bat positions its wing perpendicular to the flow to mitigate negative lift. In fact, the lift force never drops below zero during the upstroke. The vortex structures in the flow can be seen to be minimal at this point. The planform area drops to a minimum of around 260 cm2 midway through the upstroke after which it starts increasing again.

Fight power analysis

Power consumption is a particularly important metric in flight because the flight fuel must be held aloft by the flying animal or vehicle. This results in a non-linear feedback between weight and power consumption—more fuel requires more lift power which in turn requires even more fuel. Bats store fuel from food to efficiently power their wing muscles. Pennycuick [58] and others [5963] conducted early research into the power requirements of animal flight using respiratory experiments and fixed wing airfoil theory. Other attempts to understand flapping flight power consumption involved measuring metabolic power during flight [12,64], and measuring kinetic energy in the wake of a flying animal using PIV [64].

Aerodynamic simulations offer a straight forward and direct method to calculate power usage by a flying animal regardless of the complexity of the wing kinematics. At each instant in time, the aerodynamic power equals the instantaneous rate of work done by the wing on the fluid, and can be determined by, (9) where p is fluid pressure, is the wing surface normal, is the viscous stress tensor, and is the velocity of the wing surface. In the context of the immersed boundary method, Eq (9) is evaluated by summing over the discrete surface elements comprising the wings. The resulting power is shown for two complete flap cycles in Fig 17.

thumbnail
Fig 17. Aerodynamic power is plotted for two complete flap cycles along with wing tip position for context.

Power expenditure was maximum during the second half of each downstroke at around 2.5 to 3.5 Watts. The cycle averaged value over both flaps was 1.05 Watts.

https://doi.org/10.1371/journal.pone.0207613.g017

The cycle averaged power expenditure was 1.00 W and 1.10 W for the two flaps respectively. Peak power output reached 3–3.5 W in the second half of each downstroke. In the present analysis, the wing inertia and elastic strain energy are assumed to have minimal impact on power, in line with the findings of other researchers [64]. This is largely because the wing velocity diminishes to zero at the top and bottom of the stroke, thus the inertial energy is returned to the flow. Similarly, the wing stretches during the downstroke thereby storing potential energy, but as the membrane relaxes on the upstroke the energy is returned to the flow. These factors may impact the temporal accuracy of the power values, but should not affect the cycle averaged results.

Conclusion

Aerodynamic simulations of bat flight allow for the direct investigation of the fluid phenomena responsible for generating of lift and thrust by a flexible flapping membrane. The ability to analyze transient flow fields alongside detailed wing kinematics can advance the understanding of flapping flight beyond prior methods. Previous work in the field has leaned on the well-known theory of conventional airfoils which is not always applicable to bat flight. For example, the angle of attack, camber, and chord length vary significantly both along the span of a bat wing as well as temporally throughout the flap cycle in contrast to the well-defined and predictable behavior of a conventional airfoil. Other research has used PIV methods to measure the flow field in the wake of the bat, however this method has restrictions on the size of the bat and the ability to study flight maneuvers. Additionally, pseudo-steady wake models for force calculations rely on the frozen vortex hypothesis which does not always hold in animal flight.

Conducting aerodynamic simulations to numerically obtain the flow field around bat wings during flight is an attractive alternative to other methods of analysis. However, both high spatial resolution motion capture, as well as the numerical simulations remain a technical challenge. In the present work, 108 discrete marker points were tracked in order to fully capture the many degrees of freedom of the bat wings. After motion capture, the kinematic data was interfaced with a Navier-Stokes solver using the immersed boundary method as well as LES. A fluid grid of 88.1 million cells was used to resolve the unsteady flow of Re = 11,680.

Aerodynamic forces and the associated power were directly calculated from the resulting flow field using surface integrations. Coherent vorticity visualization was used in conjunction with wing area, force, and power data to investigate some of the underlying mechanism of bat flight. For example, while still in the upstroke, the bat starts generating lift due to wing rotation which results in the formation of a LEV. However, the majority of lift is generated during the downstroke when the LEV mitigates massive flow separation from the wing. During the upstroke, negative lift is completely avoided by virtue of wing rotation and a drastic reduction in planform area providing a net performance benefit.

There are still many challenges and open questions to be considered in the field of bat flight. For example, what is the interplay between wing elasticity and unsteady flow features? How do different bat wing morphologies manifest as optimizations to different types of flight (speed versus maneuverability versus efficiency)? What is the underlying control scheme to actuate flapping flight? How are bats so adept at maneuvering flight? Computational modeling can be a critical tool in uncovering answers to these and other open questions in the realm of bat flight. The 3D motion capture and computational framework described and validated in this study has the potential to analyze a wide variety of bat flights, including the study of maneuvering bat flight.

Supporting information

S1 CSV File. Time series of force and power data from the bat flight simulation.

https://doi.org/10.1371/journal.pone.0207613.s001

(CSV)

S1 Video Frames. Series of images from the bat flight video.

Images are provided from the view of two separate cameras. Each successive frame is separated by 1/120 seconds. The cameras were synchronized in their recording (i.e. frame 380 represents the same moment for both cameras).

https://doi.org/10.1371/journal.pone.0207613.s002

(RAR)

Acknowledgments

The financial support of the NSF CBET Grant No. 1510797, NSF IRES Grant No. 1658620, support from VT ICTAS/BIST Center, National Natural Science Foundation of China (Grant Nos. 11374192 & 11574183), Chinese Ministry of Education Tese Grant for international faculty exchange is gratefully acknowledged. The authors would also like to acknowledge the computational resources provided by Advanced Research Computing (ARC) at Virginia Tech.

References

  1. 1. Shyy W, Aono H, Chimakurthi SK, Trizila P, Kang CK, Cesnik CES, et al. Recent progress in flapping wing aerodynamics and aeroelasticity. Prog Aerosp Sci. Elsevier; 2010;46: 284–327.
  2. 2. Shyy W, Berg M, Ljungqvist D. Flapping and flexible wings for biological and micro air vehicles. Prog Aerosp Sci. 1999;35: 455–505.
  3. 3. Lentink D, Jongerius SR, Bradshaw NL. The scalable design of flapping micro-air vehicles inspired by insect flight. Flying Insects and Robots. 2010. pp. 185–205.
  4. 4. Ho S, Nassef H, Pornsinsirirak N, Tai YC, Ho CM. Unsteady aerodynamics and flow control for flapping wing flyers. Progress in Aerospace Sciences. 2003. pp. 635–681.
  5. 5. Chin DD, Lentink D. Flapping wing aerodynamics: from insects to vertebrates. J Exp Biol. 2016;219: 920–932. pmid:27030773
  6. 6. Norberg UM, Rayner JM V. Ecological Morphology and Flight in Bats (Mammalia; Chiroptera): Wing Adaptations, Flight Performance, Foraging Strategy and Echolocation. Philos Trans R Soc B Biol Sci. 1987;316: 335–427.
  7. 7. Swartz SM, Konow N. Advances in the study of bat flight: the wing and the wind. Can J Zool. 2015;93: 977–990.
  8. 8. Hedenstrom A, Johansson LC. Bat flight: aerodynamics, kinematics and flight morphology. J Exp Biol. 2015;218: 653–663. pmid:25740899
  9. 9. Hedenström A, Johansson LC, Spedding GR. Bird or bat: comparing airframe design and flight performance. Bioinspir Biomim. 2009;4: 15001. pmid:19258691
  10. 10. Kunz TH, Fenton MB. Bat Ecology. the University of Chicago Press. 2003. pp. 107–112.
  11. 11. Hartman F a. Some Flight Mechanisms of Bats. Ohio J Sci. 1963;63: 59–65.
  12. 12. Thomas SP, Suthers RA. The physiology and energetics of bat flight. J Exp Biol. 1972;57: 317–335.
  13. 13. Bullen R, McKenzie NL. Bat airframe design: Flight performance, stability and control in relation to foraging ecology. Aust J Zool. 2001;49: 235–261.
  14. 14. Norberg UM. Aerodynamics, kinematics, and energetics of horizontal flapping flight in the long-eared bat Plecotus auritus. J Exp Biol. 1976;65: 179–212. pmid:993701
  15. 15. Rayner JM V, Aldridge HD. Three-dimensional reconstruction of animal flight paths and the turning flight of microchiropteran bats. J Exp Biol. 1985;118: 247–265.
  16. 16. Aldridge HD. Kinematics and aerodynamics of the greater horseshoe bat, Rhinolophus ferrumequinum, in horizontal flight at various flight speeds. J Exp Biol. 1986;126: 479–97. pmid:3806000
  17. 17. Aldridge HD. Turning flight of bats. J Exp Biol. 1987;128: 419–425. pmid:3559468
  18. 18. Aldridge HD. Body accelerations during the wingbeat in six bat species: the function of the upstroke in thrust generation. J Exp Biol. 1987;130: 275–293.
  19. 19. Aldridge HD. Flight kinematics and energetics in the little brown bat, Myotis lucifugus (Chiroptera: Vespertilionidae), with reference to the influence of ground effect. J Zool. 1988;216: 507–517.
  20. 20. Aldridge HD. Vertical Flight in the Greater Horseshoe Bat Rhinolophus Ferrumequinum. J Exp Biol. 1991;157: 183–204.
  21. 21. Norberg UM, Kunz TH, Steffensen JF, Winter Y, Von Helversen O. The Cost of Hovering and Forward Flight in a Nectar-Feeding Bat, Glossophaga Soricina, Estimated From Aerodynamic Theory. J Exp Biol. 1993;182: 207–227. pmid:8228780
  22. 22. Tian X, Iriarte-Diaz J, Middleton K, Galvao R, Israeli E, Roemer A, et al. Direct measurements of the kinematics and dynamics of bat flight. Bioinspir Biomim. 2006;1: S10–S18. pmid:17671313
  23. 23. Riskin DK, Willis DJ, Iriarte-Diaz J, Hedrick TL, Kostandov M, Chen J, et al. Quantifying the complexity of bat wing kinematics. J Theor Biol. 2008;254: 604–615. pmid:18621062
  24. 24. Hedenstrom A, Johansson LC, Wolf M, von Busse R, Winter Y, Spedding GR. Bat Flight Generates Complex Aerodynamic Tracks. Science (80-). 2007;316: 894–897. pmid:17495171
  25. 25. Muijres FT, Johansson LC, Barfield R, Wolf M, Spedding GR, Hedenstrom A. Leading-Edge Vortex Improves Lift in Slow-Flying Bats. Science (80-). 2008;319: 1250–1253. pmid:18309085
  26. 26. Johansson LC, Wolf M, von Busse R, Winter Y, Spedding GR, Hedenström A. The near and far wake of Pallas’ long tongued bat (Glossophaga soricina). J Exp Biol. 2008;211: 2909–18. pmid:18775928
  27. 27. Johansson LC, Wolf M, Hedenström A. A quantitative comparison of bird and bat wakes. J R Soc Interface. 2010;7: 61–66. pmid:19324669
  28. 28. Hubel TY, Hristov NI, Swartz SM, Breuer KS. Time-resolved wake structure and kinematics of bat flight. Anim Locomot. 2010; 371–381.
  29. 29. Hedenström A, Muijres FT, Von Busse R, Johansson LC, Winter Y, Spedding GR. High-speed stereo DPIV measurement of wakes of two bat species flying freely in a wind tunnel. Animal Locomotion. 2010. pp. 361–370.
  30. 30. Spedding GR, Hedenström A. PIV-based investigations of animal flight. Animal Locomotion. 2010. pp. 187–201.
  31. 31. Gutierrez E, Quinn DB, Chin DD, Lentink D. Lift calculations based on accepted wake models for animal flight are inconsistent and sensitive to vortex dynamics. Bioinspir Biomim. 2016;12: 16004. pmid:27921999
  32. 32. Dabiri JO. On the estimation of swimming and flying forces from wake measurements. J Exp Biol. 2005;208: 3519–3532. pmid:16155224
  33. 33. Dong H, Koehler C, Liang Z, Wan H, Gaston Z. An Integrated Analysis of a Dragonfly in Free Flight. AIAA J. 2010; 1–10.
  34. 34. Li C, Dong H. Wing kinematics measurement and aerodynamics of a dragonfly in turning flight. Bioinspir Biomim. 2017;12: 26001. pmid:28059781
  35. 35. Wan H, Dong H, Gai K. Computational investigation of cicada aerodynamics in forward flight. J R Soc Interface. 2014;12: 20141116–20141116. pmid:25551136
  36. 36. Liu G, Dong H, Li C. Vortex dynamics and new lift enhancement mechanism of wing–body interaction in insect forward flight. J Fluid Mech. 2016;795: 634–651.
  37. 37. Ren Y, Dong H, Deng X, Tobalske B. Turning on a dime: Asymmetric vortex formation in hummingbird maneuvering flight. Phys Rev Fluids. 2016;1: 50511.
  38. 38. Song J, Luo H, Hedrick TL. Three-dimensional flow and lift characteristics of a hovering ruby-throated hummingbird. J R Soc Interface. 2014;11: 20140541–20140541. pmid:25008082
  39. 39. Viswanath K, Nagendra K, Cotter J, Frauenthal M, Tafti DK. Straight-line climbing flight aerodynamics of a fruit bat. Phys Fluids. 2014;26.
  40. 40. Deetjen ME, Biewener AA, Lentink D. High-speed surface reconstruction of a flying bird using structured light. J Exp Biol. 2017;220: 1956–1961. pmid:28348041
  41. 41. Tobalske BW, Warrick DR, Clark CJ, Powers DR, Hedrick TL, Hyder GA, et al. Three-dimensional kinematics of hummingbird flight. J Exp Biol. 2007;210: 2368–2382. pmid:17575042
  42. 42. Altshuler DL, Dickson WB, Vance JT, Roberts SP, Dickinson MH. Short-amplitude high-frequency wing strokes determine the aerodynamics of honeybee flight. Proc Natl Acad Sci U S A. 2005;102: 18213–8. pmid:16330767
  43. 43. Bozkurttas M, Mittal R, Dong H, Lauder G V., Madden P. Low-dimensional models and performance scaling of a highly deformable fish pectoral fin. J Fluid Mech. 2009;631: 311.
  44. 44. Bender MJ, McClelland HG, Bledt G, Kurdila A, Furukawa T, Mueller R. Trajectory Estimation of Bat Flight Using a Multi-View Camera System. Model Simul Technol Conf SciTech 2015. 2015; 1–13.
  45. 45. Bender MJ, McClelland HM, Kurdila A, Mueller R. Recursive Bayesian Estimation of Bat Flapping Flight Using Kinematic Trees. AIAA Model Simul Technol Conf. 2016; 1–12.
  46. 46. Svoboda T, Martinec D, Pajdla T. A convenient multicamera self-calibration for virtual environments. PRESENCE teleoperators virtual Environ. 2005;14: 407–422.
  47. 47. Fan X, Windes P, Tafti D, Sekhar S, Bender M, Kurdila A, et al. Proper orthogonal decomposition of straight and level flight kinematics in an insectivorous bat. 2018 AIAA Model Simul Technol Conf. 2018; 1–22.
  48. 48. Tafti DK. GenIDLEST: A scalable parallel computational tool for simulating complex turbulent flows. ASME-PUBLICATIONS-FED. 2001. pp. 347–356.
  49. 49. Tafti DK. Time-accurate techniques for turbulent heat transfer analysis in complex geometries. In: Amano R, Sunden B, editors. Computational Fluid Dynamics and Heat Transfer. Southampton, UK: WIT Press; 2011. pp. 217–264.
  50. 50. Germano M, Piomelli U, Moin P, Cabot WH. A dynamic subgrid-scale eddy viscosity model. Phys Fluids A Fluid Dyn. 1991;3: 1760–1765.
  51. 51. Najjar FM, Tafti DK. Study of discrete test filters and finite difference approximations for the dynamic subgrid-scale stress model. Phys Fluids. 1996;8: 1076–1088.
  52. 52. Nagendra K, Tafti DK, Viswanath K. A new approach for conjugate heat transfer problems using immersed boundary method for curvilinear grid based solvers. J Comput Phys. 2014;267: 225–246.
  53. 53. Wolf M, Johansson LC, von Busse R, Winter Y, Hedenstrom A. Kinematics of flight and the relationship to the vortex wake of a Pallas’ long tongued bat (Glossophaga soricina). J Exp Biol. 2010;213: 2142–2153. pmid:20511529
  54. 54. Hubel TY, Riskin DK, Swartz SM, Breuer KS. Wake structure and wing kinematics: the flight of the lesser dog-faced fruit bat, Cynopterus brachyotis. J Exp Biol. 2010;213: 3427–3440. pmid:20889823
  55. 55. Hubel TY, Hristov NI, Swartz SM, Breuer KS. Wake structure and kinematics in two insectivorous bats. Philos Trans R Soc B Biol Sci. 2016;371: 20150385. pmid:27528775
  56. 56. Hubel TY, Hristov NI, Swartz SM, Breuer KS. Changes in kinematics and aerodynamics over a range of speeds in Tadarida brasiliensis, the Brazilian free-tailed bat. J R Soc Interface. 2012;9: 1120–1130. pmid:22258554
  57. 57. Muijres FT, Johansson LC, Winter Y, Hedenström A, Muijres FT, Johansson LC, et al. Comparative aerodynamic performance of flapping flight in two bat species using time-resolved wake visualization. J R Soc Interfac Interfac. 2011; pmid:21367776
  58. 58. Pennycuick C. Bird Flight Perfomance. Oxford University Press; 1989.
  59. 59. Withers BYPC Timko PL. The significance of ground effect to the aerodynamic cost of flight and energetics of the black skimmer (Rhyncops Nigra). J Exp Biol. 1977;70: 13–26.
  60. 60. Rayner JM V., Viscardi PW, WARD S, Speakman JR. Aerodynamics and energetics of intermittent flight in birds. Am Zool. 1999;39: 4A–5A.
  61. 61. Rayner JM V. A new approach to animal flight mechanics. J Exp Biol. 1978;80: 17–54.
  62. 62. Ellington C. Limitations on animal flight performance. J Exp Biol. 1991;91: 71–91.
  63. 63. Ellington CP. The Aerodynamics of Hovering Insect Flight. I. The Quasi-Steady Analysis. Philos Trans R Soc B Biol Sci. 1984;305: 1–15.
  64. 64. von Busse R, Waldman RM, Swartz SM, Voigt CC, Breuer KS. The aerodynamic cost of flight in the short-tailed fruit bat (Carollia perspicillata): comparing theory with measurement. J R Soc Interface. 2014;11: 20140147. pmid:24718450