Forward flight of birds revisited. Part 1: aerodynamics and performance

This paper is the first part of the two-part exposition, addressing performance and dynamic stability of birds. The aerodynamic model underlying the entire study is presented in this part. It exploits the simplicity of the lifting line approximation to furnish the forces and moments acting on a single wing in closed analytical forms. The accuracy of the model is corroborated by comparison with numerical simulations based on the vortex lattice method. Performance is studied both in tethered (as on a sting in a wind tunnel) and in free flights. Wing twist is identified as the main parameter affecting the flight performance—at high speeds, it improves efficiency, the rate of climb and the maximal level speed; at low speeds, it allows flying slower. It is demonstrated that, under most circumstances, the difference in performance between tethered and free flights is small.


Introduction
Aerodynamics and flight mechanics of flapping flight have drawn research attention since the beginning of the aviation era. The complexity of aerodynamic models involved progressively increased, until recent advances in computing power have made full Reynolds-averaged Navier-Stokes simulations within reach [1]. It seems, however, that in the race for fidelity, a few fundamental problems became buried under excessive details. Two of these problems, performance and short-term dynamic stability of birds in forward flight, are revisited in this study with the simplest aerodynamic model feasible. The performance is addressed in this part; the short term dynamic stability is addressed in part 2 [2]. The aerodynamic model, serving both parts, is constructed herein. By 'simplest model feasible', we understand a model that can furnish the aerodynamic loads in a closed analytical form, and is accurate enough to capture their behaviour in flapping flight. Taking the cue from [3], we construct this model in the framework of the basic lifting line theory. The present model differs from that of [3]   for the wings' motion, in furnishing the forces acting on a single wing as well as their first moments, and in remarkable simplicity of its final expressions. Using the vortex lattice method to provide the reference, the model is shown in § §4.5 and 5.2 to be accurate enough for all the purposes of this two-part study. It intrinsically limits it, however, to those flyers that generate lift without flow separation-that is, sufficiently large birds and bats (but not insects).
A wing invariably twists during flapping. The effect of twist on flight performance was addressed in many papers (e.g. [1,3,4]), but the twist never received the credit of being the most important parameter affecting the performance. It could have happened because of the historical definition of the propulsion efficiency: the ratio of the average power made good (the product of force produced by the wings in the direction of flight and airspeed) and the average power spent (see for example [4,5]). This definition ignores the double role played by flapping wings in flight. Stopping the wings makes the force generated by them in the direction of flight negative-in fact, this force becomes drag, parasite and induced combined. In comparison, it is commonly accepted that stopping the propeller of an aeroplane makes its thrust vanish-the drag of the aeroplane's wings is considered an inseparable part of the aeroplane's total drag. As propulsion efficiency of a propeller is defined irrespective of the aerodynamic characteristics of the aeroplane it propels, so the propulsion efficiency of flapping wings should be defined irrespective of the aerodynamic characteristics of the bird in non-flapping flight. Propulsion efficiency and its dependence on twist are addressed in §5.1.
A common measure of performance of a fixed wing aeroplane is the specific excess power [6]essentially, it is the maximal sustained rate of climb at a given airspeed. We did not see the use of this measure in relation to flight performance of birds. Specific excess power of an aeroplane is limited by the engine power throughout the flight envelope. Specific excess power of a bird is also limited by the maximal power the bird can generate, but only at high speeds. Increasing the power at low speeds stalls the wings, putting an additional limit on performance. Specific excess power of a bird and its dependence on twist are addressed in §5. 4.
Bird's body invariably pitches and heaves during flapping. Effects of the body motion on performance are addressed in § §6.1-6.3.

Kinematics
Consider a simplified symmetric bird in symmetric flight with constant velocity v. The mass of the bird is m; the density of the air in which it flies is ρ; acceleration of gravity is g; the length of a single wing (the semi-span) is s; its area is S; the aspect ratio of the two wings is A = 2s 2 /S. s, v, vs, s/v, v/s, ρsv 2 , ρs 2 v 2 , ρSv 2 , ρSv 2 s and ρSv 3 will serve as convenient units of length, velocity, circulation, time, frequency, force per unit span, moment per unit span, force, moment and power, respectively. Note that although S is half the quantity commonly used as the wing area, the units of force and power are standard. Use of dimensionless quantities is implicitly understood hereafter. Should a dimensional quantity (other than ρ, g, m, s, S and v) be required, it will be marked by an asterisk. A list of nomenclature can be found in table 1. Each wing is allowed to flap, sweep fore and aft, twist, heave and pitch. It is assumed that the wing twists in such a way that its sections do not deform and remain parallel to each other; moreover, the twist axis crosses all sections at their respective quarter-chord points and remains straight at all times. The sweep angle of the twist axis is λ (positive aft), flapping angle is φ (positive down); the twist angle is α g (positive for leading edge up); pitch angle relative to the average flight path is τ (positive for nose up); vertical translation of the twist axis is h (figure 1). It is assumed that the twist varies linearly along the span, with α g = α g0 + α g1 y, (3.1) where y ∈ (0, 1) is the spanwise coordinate; the description here pertains to the right wing. The twist can be active (through muscle contraction) or passive (through aerodynamic twisting moment); no attempt is made to model its intricate details. Two right-handed Cartesian reference frames will be used. Frame L is rigidly connected to a local chord; its y-axis coincides with the twist axis of the right wing; its x-axis points backwards parallel to the chord; its z-axis points upwards and its origin rests in the symmetry plane of the bird (where the two wings meet). Frame C follows the bird with constant velocity (v) along a straight path; its x-axis points backwards along that path; its y-axis points right and its z-axis points upwards, through the origin of L. Unit vectors along the axes of L are e L x , e L y , e L z ; unit vectors along the axes of C are e x = e C x , e y = e C y , e z = e C z . Heave is manifested in time dependence of the distance h between the origins of L and C. Rotation         Figure 1. Reference frames and wing motion parameters. The wing is allowed to heave (h), flap (φ), sweep (λ), pitch (τ ) and twist (α g ) about the quarter-chord line. The twist shown is highly exaggerated; sweep and pitch are not shown. Frame L is rigidly connected to a local section of the wing; the axes shown correspond to the mid span of the right wing. Frame C is an inertial reference frame that follows the bird along a straight path at the distance h beneath the wing. from C to L is a series of four Euler's rotations: about the y-axis through angle τ ; about the new z-axis through angle −λ; about the new x-axis through angle −φ and about the new y-axis through angle α g . A comparable approach can be found in [5].
Explicit expressions for the components of e L x , e L y and e L z in C are lengthy, and hence are not written here for the general case. In the particular case where all angles are small when compared with unity and e L z = (α g + τ )e x + φe y + e z + · · · , (3.4) where the ellipses stand for the higher-order terms with respect to angles and their time derivatives. Concurrently, the velocity of a point on the wing relative to C is where x and y are (by interpretation) the distances from the twist and flapping axes respectively, and an over-dot stands for derivative with respect to (dimensionless) time.

Assumptions
It is assumed that A −1 , α g , τ ,ḣ, φ,φ, λ andλ are small when compared with unity. The first assumption underlies the lifting line theory; the remaining assumptions underlie linearization. For the sake of simplicity of the following discussion, it is assumed that α g , τ ,ḣ andφ are of comparable magnitudes, say . The magnitude of λ is assumed not to exceed A −1 and φ;λ is assumed to be a second-order quantity with respect to φ, and their products.
It is postulated that a vortical wake exists past the wing, starting at the trailing edge and extending to infinity. It is assumed that the vorticity is constant along that portion of the wake adjacent to the trailing edge that affects the flow over the wing. This assumption implies that the flapping frequency ω is sufficiently small; it is plausible-in fact, it will be shown by example-that ω can be of the order of unity. The aerodynamic model developed herein is coherent only in the leading order with respect to φ, , and their products with ω. It is assumed that the wing has an elliptical plan-form with chord length prescribed by c 0 = 8/π A, because S = s 2 c 0 π/4 and A = 2s 2 /S by definition. Exploiting the symmetry of the problem, the range of y in all subsequent equations is extended to (−1, 1); negative values corresponding to the left wing.

Fundamentals
As already mentioned in §2, the aerodynamic model for this study is based on the classical (quasi-steady) lifting line theory ( [7], p. 586). In brief, this theory associates the lift of a wing section, represented by the right-hand side of the following equation, with the lift of an equivalent vortex, represented by the left-hand side Γ is the circulation of that vortex, a = 2π is the lift-slope coefficient of the wing section, α is the effective angle of attack the wing section relative to unperturbed fluid (that can be considered known), and α i is the angle of attack induced by the wake (figure 2). The closure of (4.2) is obtained by relating α i with Γ by Biot-Savart's law ( [7], p. 94) which transforms (4.2) into an integro-differential equation for Γ . Once solved, the aerodynamic loads follow by quadratures. The non-intuitive elements in this paradigm are the angles α and α i . In the lifting-line limit, the wing reduces to a vortex in the y-z plane whereas the wake reduces to the sheet of vortices starting at the wing vortex and extending to infinity in the positive x-direction [8]. Coherent with this model, α is approximated by its leading-order term (with respect to the aspect ratio), lim For a non-cambered wing, it yields α = τ + α g +φ|y| −ḣ + · · · = α 0 + α 1 |y| + · · · , (4.3) and by (3.1), (3.4) and (3.5); it can be extended ad hoc for a cambered wing by defining α g as the angle between the x-axis of C and the zero-lift line (rather than the chord) of the respective section. α 0 will be recognized as the spanwise uniform constituent of the angle of attack-it is associated with pitch and heave of the body and twist at the shoulder; α 1 is the gradient of the angle of attack along the span-it is associated with the flapping rate and with spanwise-variable twist. α i is approximated by its leading-order term with respect to the aspect ratio as well-that is, the normal-to-the-wing (reduced) velocity component in the y-z plane induced by the wake [8]. Because φ is small when compared with unity, and the wake vorticity is supposedly constant along that portion of the wake that affects the flow about the wing (see §4.1), in which α is given by (4.3), c is given by (4.1), and, in general, a = 2π . Its relevant solution is where θ = cos −1 (−y), (4.9)

Forces and moments
The force per unit span acting on a section of the right wing, f, and its couple about the origin of L, m, (4.13) and m = ye L y × f. (4.14) With (3.3) and (3.5), they take on the explicit forms f = Γ ((ḣ −φy + α i )e x + φe y + e z + · · · ) (4.15) and m = yΓ (e x − (λ + φ(ḣ −φy + α i ))e y − e z (ḣ −φy + α i ) + · · · ), (4. 16) where the ellipses stand for higher-order terms in φ and ; they are omitted hereafter. Consequently, is the power needed to move both wings. Because in the dimensionless representation, the power made good,Ẇ mg , and the thrust, T, are equivalent, the conjunction of (4.21) and (4.18), where is coherent with the interpretation of D i as the induced drag. In (4.19), the first term in the right-hand side, the one involving λ, reflects the contribution of the lift to the pitching moment; the remaining terms reflect the contribution of the thrust. With   (independent of the aspect ratio) in Despite its simple form (it equals −1/3 by (A 7)), I 11 was left unevaluated in many of the subsequent equations to facilitate tracking of the particular terms. Concluding this list is the explicit expression for the power required to move the wings, it follows from (4.21) by (4.26) and (4.24). Explicit expression for the power made good,Ẇ mg = T, immediately follows from (4.32) and (4.25) by (4.22), and hence is missed out.

Local angle of attack
The local angle of attack during flapping, α − α i , can be found directly, from (4.3) and (4.10)-(4.12):  The equivalence between the two approaches is proved in appendix C. Maximal angle of attack along the span, α c , will be needed in §5.3. It occurs at the wing tips, θ = 0 and θ = π -this conjecture is proved in appendix D; consequently, where the first term in the right-hand side follows the respective term in (4.34) by (4.11), whereas is yet another function of the aspect ratio. Its Padé approximation can be found in (B 2).

Corroboration
The accuracy of (4.   Recalling that the present model ignores unsteady and non-planar effects, its adequacy in predicting aerodynamic loads implies that at least in some cases, relevant to the flight of birds, these effects are indeed secondary to the effects of finite span. We will return to this point in §5.2.

Propulsion efficiency
The preceding sections were based on the traditional definitions of thrust and power made good. As already mentioned in §2, these definitions have to be modified for flapping flight, where the wings serve the dual role of providing lift and generating thrust. When the wings stop flapping, traditionally defined thrust becomes negative, it turns into induced drag-(see (4.22)). If it were a fixed wing aircraft, stopping the engine would have made its thrust vanish-the induced drag would have been an inseparable part of the aircraft drag, rather than part of its thrust. This inconsistency is removed here by differentiating between the proper thrust T = T +D i , (5.1) and the proper excess thrust, The propulsion efficiency can now be defined as the ratio of the average proper power made good, Ẇ mg = T , to the average power needed to move the wings, Ẇ : the angular brackets mark the respective period averages. The two forms follow by (5.1) and by the variant of (4.22), respectively. The first one is better suited for numerical computations, where the average thrust and the induced drag in the adjoint flight are readily available. The second one is clearer conceptually, as it explicitly associates the efficiency with the drag added by flapping, D i −D i . The ratio, T ex / Ẇ , commonly used instead of (5.4)-for example, in [4,5]-is a viable figure of merit, but it is not efficiency.
Marking the respective parameters in the adjoint flight by over-bars, the adjoint flight is formally defined as the flight in whichᾱ 0 = α 0 ,ᾱ 1 = α 1 . (5.5) by (4.24) and (4.25), whereas . Three assumptions will greatly simplify the following discussion. The first one fixes the twist at the shoulder withα g0 = 0. (5.9) The second one slaves the wing's twist to the flapping rate with where ε ∈ (0, 1) is a certain proportionality coefficient (as in [1,4]  it is equivalent to placing the bird on a sting in a wind tunnel. This assumption will be released in §6.1. A posteriori, its effect on performance is small. Because of (5.9), (5.13) and (4.4), Note that in the framework of the linear theory, the pitch (τ ) and the spanwise-constant constituent of the twist (α g,0 ) are equivalent. With (5.11)-(5.14), equations (5.8) and (4.32) reduce to expressions in the right-hand side follow those on their left by (A 7) and (4.31). Both parameters are shown in figure 5. The efficiency increases with increasing twist and aspect ratio; the power dramatically diminishes with increasing twist, and slightly increases with the aspect ratio. Increasing the aspect ratio reduces the velocity induced by the wake; smaller induced velocity increases the undulatory constituent of the lift and reduces the respective constituent of the drag. Larger lift increases the power needed to flap the wings; smaller drag improves the efficiency. Increasing the twist reduces the undulatory constituents of the lift and the drag alike. Power required is associated with lift; power losses are associated with drag. Because lift is linear in the angle-of-attack and drag is quadratic, increasing the twist improves the efficiency.  A flapping wing leaves vortical wake behind it. Some of the vortices comprising the wake are associated with finiteness of the wing; others are associated with the wing's motion. Vortices of the first type are roughly parallel to the direction of the flow relative to the wing; they are created as long as the circulation varies along the span. Vortices of the second type are perpendicular to the direction of the flow; they are created any time there is a change in the flow about the wing. Both types of the vortices carry energy with them; energy loss to the wake is manifested as drag in the definition of the propulsion efficiency (5.4). Equation (5.15) accounts only for the energy carried away by vortices of the first type, and because the energy carried away by the vortices of the second type increases with the flapping frequency, the accuracy suffers.

Nonlinear effects
The detrimental effect of increasing frequency on the propulsion efficiency offsets somewhat the beneficial effect of the increasing twist. For the same flapping amplitude and the same average power, increasing the twist parameter from 0.3 to 0.7, necessitates a 50% increase in the flapping frequency. This conjecture follows from (5.16), because in harmonic flapping φ 2 = φ 2 0 ω 2 /2. This increase in the twist parameter increases the propulsion efficiency from the bottom line in figure 5b to the top linesay, from 0.83 to 0.93 (at A = 8). The 50% increase in the flapping frequency decreases the efficiency by the difference between the circles (ω = 0.8) and the squares (ω = 1.2) in figure 5b, approximately by 0.015. Hence, the decrease of the propulsion efficiency with frequency is a second-order effect-in fact, consistent with the assumptions underlying the present model.
The adjoint flight has been formally defined in (5.5) by specifying the angle of attack and the twist. In the framework of the linear theory, this definition yieldsL = L by (4.24) (see (5.6)). However, in the framework of a nonlinear theory-as the one represented by our numerical solution-the lift in the adjoint flight also depends on the dihedral angle of the wings,φ. To obtainL = L , |φ| should be, roughly, two-thirds of the flapping amplitude; stopping the wings atφ = 0 yieldsL > L . We could have defined the adjoint flight differently, by specifyingᾱ 1 = α 1 ,L = L andφ = 0. It would have changed nothing in the linear theory, but it would have impliedᾱ 0 < α 0 (and henceD < lim In turn, deficient drag in the adjoint flight would have yielded negative propulsion efficiency at slow flapping rates. In choosing (5.5), we have avoided this outcome, but there are pros and cons for each one of these two definitions.

Level flight
In straight-and-level-constant-speed flight-it will be referred to as 'SL flight' below-the average lift counterbalances the weight, is shown in figure 6a. In turn, φ 2 is related to flapping amplitude (φ 0 ) and angular frequency (ω) by where k φ is a certain parameter depending on the flapping pattern. It equals 1/2 when φ varies with time as a sine, and 1/(π 2 τ d (1 − τ d )) when it varies as a 'saw-tooth' with τ d being the relative part of the down stroke in the flapping period. Combination of (5.22) and (5.24) furnishes the flapping frequency in SL flight: where C 1 and C 2 are certain constants. Equation (3) in reference [9] is the Taylor series of this equation. Different behaviour of the flapping frequency with airspeed will imply that the twist parameter ε and/or the flapping amplitude φ 0 and/or the flapping pattern k φ change in flight.
The drag-to-lift ratioD/L has a minimum 2 √ D 0 /π A at v = mg/ρS/ 4 √ π AD 0 ; the existence of this (shallow) minimum is manifested in the left-hand side of the flapping frequency curves in figure 7d (at 11 m s −1 ). 3 Assuming that a bird cruises where the cost of locomotion is minimal, the drag-to-lift ratio at cruise should not significantly differ from its minimal value. The flapping frequency at cruise is, 3 By differentiating (5.27) with respect toL, and equating the result to zero, one will find thatD/L has a minimum atL = therefore, by (5.26) and by definition of the aspect ratio. In (5.28), H 1 (A, ε) √ A is only a weak function of aspect ratio and twist (figure 6b), whereas k φ is, practically, a universal constant. In fact, with k ω = 1.82, approximates (5.28) to within 9% for every A in (6,16) and every ε in (0.1, 0.9). Equation (5.29) can be seen as a rational variant of equations (10) in [10] and (9) in [11].  figure 8b. We found no reliable data on the twist parameter. Nonetheless, based on a few available observations [12][13][14][15], φ 0 is of the order of 0.5 rad. It implies that the twist parameter ε for all the species represented in figure 8 exceeds 0.5.

Specific excess power
Common measures of flight performance are the specific excess thrust, F ex , and the specific excess power, P * s [6]. Specific excess thrust is defined as the ratio between the proper excess thrust (the difference between the maximal available thrust and drag) and weight. Specific excess power is defined as the At given flight conditions, F ex can be interpreted as the level acceleration (in 'g' units) and hence is equivalent to the maximal sustained climb angle (if it is sufficiently small); P * s can be interpreted as the rate of change of the energy altitude and hence is equivalent to the maximal sustained climb rate. Formally, where P max is the maximal available power. When F ex = 0, equation (5.31) becomes a variant of (5.18) or a variant of (5.19).
In flapping flight, the maximal available power can be limited either by the maximal sustained power, P * a , or by stall. In the first case, P max = P * a ρSv 3 , (5.32) and, consequently, Because for a constant ε, η is independent of airspeed (see (5.15)), equations (5.33) and (5.30) are the same as for a propeller-driven fixed-wing aeroplane. An example is shown in figure 7a,b. Both F ex and P * s diminish with airspeed; they vanish at the maximal speed that can be obtained in horizontal flight (in this example, 27-29 m s −1 , depending on the twist parameter).
The stall limit has no analogy with fixed-wing aeroplanes. It was shown in §4.4 that the maximal angle of attack along the span is this equation follows (4.35) by (5.11), (5.12), (5.14) and (5.6). If α c is not to exceed the onset of stall at α max = L max /2π ,φ should be limited from above bẏ . (5.35) This limit onφ can be translated into a limit kφφ 2 max on the variance φ 2 of the flapping rate (kφ is akin to k φ in (5.24)-both equal 1/2 if the flapping is harmonic); and, in turn, into a limit on the maximal power that can be supplied to the wings without stalling a part of them (see (5.16)). Consequently, by (5.31) and (5.15). It manifests a rapidly increasing function of airspeed (figure 7a). Reiterating, a bird is limited by stall at low speed, and by power at high speed. The best climb rate and the best climb angle are achieved where both limitations meet. Increasing the twist decreases the minimal flight speed (see (5.37)), and increases both the maximal flight speed and the maximal climb rate (see (5.15) and (5.33))). Increasing the twist above ε = 0.5 has practically no effect on the specific power-and hence on the climb rate-but has a dramatic effect on the reduced frequency, especially at low speed. Increasing the reduced frequency eventually makes the velocity of the wing owing to flapping comparable to the flight velocity (figure 7c), and it is here where the present theory fails. Lowspeed phases of flight require different analysis (e.g. [16][17][18]). Nonetheless, the stall limit does have significance-it marks the transition between forward and hovering flight. Flight beyond the stall limit is possible only by tilting the flapping plane so as to direct most of the thrust upwards [19].

.1. Formulation
We release now the degrees of freedom of the bird's body (that is, take the bird off the sting) and seek the resulting effect on performance. The degrees of freedom involve heave and pitch; the relationship between the two depends on the active control strategy adopted by the bird. The two most obvious strategies, keeping the angle-of-attack (α 0 ) constant, and keeping the pitch angle (τ ) constant, are addressed in § §6.2 and 6.3 below. For the sake of simplicity, we assume that the twist at the shoulder does not change with time-that is, we adopt (5.9); we also assume that the mass of the wings is negligible when compared with the mass of the body.
Under these assumptions, the heave is governed by the equation where is a parameter. It can be interpreted as the ratio of the lift slope coefficient, L ,α = π AA 1 (see (4.24)), and the reduced mass, m/ρSs. For most of the species compiled in appendix D of part 2, s 1 varies in the interval (0.1, 1).
If φ = φ 0 cos ωt, by (6.2), and, concurrently, by (5.11). Combination of (6.9), (6.11), (6.12) and (6.13) furnishes a quadratic equation for (ω/s 1 ) 2 : is a parameter; the right-hand side of (6.15) follows by (6.3). The relevant solution of (6.14) is Substituting it in (6.12), and the result in (6.10) and (6.11), yields K η and KẆ as functions of d. They are shown in figure 9. The limit d → ∞ corresponds to m → ∞ by (6.15), and hence lim d→∞ḣ = 0 by (6.2). Indeed, in this case, both correction factors tend to unity. For the hypothetical bird of appendix G, d ≈ 300 Ẇ . Because Ẇ is comparable to the (reduced) drag, a few hundredths, d turns out to be of the order of 10. It implies negligible effects of free heave on the propulsion efficiency and the input power. Nonetheless, a combination of small twist, low wing loading and large flapping amplitude may reduce d by an order of magnitude. Should this be the case, flying with constant pitch may still have an imperceptible effect on efficiency, but may have a profound effect on the input power.

Concluding remarks
An aerodynamic model based on the lifting line theory seems to provide a good balance between accuracy ( § §4.5 and 5.2) and simplicity ( § §4.3 and 4.4). Its main limitation is the low-speed phases of  flight, where the velocity of the wing owing to flapping is comparable to (or exceeds) the flight velocity. Better theories in this respect (e.g. [16]) are too complicated to be used effectively in analysis of the type that was carried through here, and will be carried through in part 2.
Wing twist was identified as the central parameter affecting performance ( § §5.1 and 5.4). It improves efficiency, and hence improves both the maximal rate of climb and the maximal level speed. It also reduces the angle of attack at the wing's tip, and hence allows the bird to fly slower. For 46 species of birds analysed in this study, the twist parameter at cruise should have been more than 0.5; that is, the twist at the tip should have been larger than half the angle of attack induced by flapping.
The flapping frequency required to sustain a level-constant-speed flight is proportional to the square root of the drag-to-lift ratio in the adjoint flight. It is minimal when the drag-to-lift ratio is minimal-in fact, where the cost of locomotion is minimal.
The following list summarizes the most important of the odd 100 equations of the paper. Aerodynamic loads: Data accessibility. This paper uses no external data. Two files have been uploaded with this article as the supplementary material. One file contains 21 cases comparing the aerodynamic theory developed herein with numerical simulations based on the vortex lattice method. The other file contains a verification of the numerical code that was used to generate these simulations.

Appendix G. Data for performance analysis
A hypothetical mid-sized bird was generated for analysis, resembling Cape petrel (Daption capense) in mass and dimensions (table 3). The maximal lift coefficient chosen for this bird probably exceeds the maximal lift coefficient that can be expected at the relevant Reynolds number (which is about 70 000 at 10 m s −1 ), even when corrected for possible delay of stall because of flapping. Increasing the maximal lift coefficient allowed for a clearer demonstration of the minimum of the flapping frequency in level flight. The parasite drag coefficient was chosen on the lower side of what could have been expected based on parasite drag of a typical profile (see appendix D of part 2); this choice served the same purpose of clearer visualization. Maximal power was assumed after [12]; the flapping was harmonic, φ(t) = φ 0 cos ωt, with (reduced) frequency ω and amplitude φ 0 ; there was no pitch and no heave; the air density was standard, ρ = 1.225 kg m −3 .