# Aerodynamic Ground Effect in Fruitfly Sized Insect Takeoff

Dmitry Kolomenskiy, Masateru Maeda, Thomas Engels, Hao Liu, Kai Schneider, Jean-Christophe Nave

► To cite this version:

Dmitry Kolomenskiy, Masateru Maeda, Thomas Engels, Hao Liu, Kai Schneider, et al. Aerodynamic Ground Effect in Fruitfly Sized Insect Takeoff. PLoS ONE, 2016, 11 (3), pp.e0152072. 10.1371/journal.pone.0152072. hal-01299261

HAL Id: hal-01299261

https://hal.science/hal-01299261v1

Submitted on 20 Sep 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Aerodynamic Ground Effect in Fruitfly Sized Insect Takeoff

Dmitry Kolomenskiy¹,⁶*, Masateru Maeda², Thomas Engels³,⁴, Hao Liu¹,⁵*,
Kai Schneider³, Jean-Christophe Nave⁶

¹ Graduate School of Engineering, Chiba University, Chiba, Japan, ² Department of Biology, Lund University, Lund, Sweden, ³ M2P2-CNRS, Université d'Aix-Marseille, Marseille, France, ⁴ Institut für Strömungsmechanik und Technische Akustik (ISTA), TU Berlin, Berlin, Germany, ⁵ Shanghai-Jiao Tong University and Chiba University International Cooperative Research Center, Shanghai, China, ⁶ Department of Mathematics and Statistics, McGill University, Montreal, QC, Canada

* dkolom@chiba-u.jp (DK); hliu@faculty.chiba-u.jp (HL)

## Abstract

Aerodynamic ground effect in flapping-wing insect flight is of importance to comparative morphologies and of interest to the micro-air-vehicle (MAV) community. Recent studies, however, show apparently contradictory results of either some significant extra lift or power savings, or zero ground effect. Here we present a numerical study of fruitfly sized insect takeoff with a specific focus on the significance of leg thrust and wing kinematics. Flapping-wing takeoff is studied using numerical modelling and high performance computing. The aerodynamic forces are calculated using a three-dimensional Navier-Stokes solver based on a pseudo-spectral method with volume penalization. It is coupled with a flight dynamics solver that accounts for the body weight, inertia and the leg thrust, while only having two degrees of freedom: the vertical and the longitudinal horizontal displacement. The natural voluntary takeoff of a fruitfly is considered as reference. The parameters of the model are then varied to explore possible effects of interaction between the flapping-wing model and the ground plane. These modified takeoffs include cases with decreased leg thrust parameter, and/or with periodic wing kinematics, constant body pitch angle. The results show that the ground effect during natural voluntary takeoff is negligible. In the modified takeoffs, when the rate of climb is slow, the difference in the aerodynamic forces due to the interaction with the ground is up to 6%. Surprisingly, depending on the kinematics, the difference is either positive or negative, in contrast to the intuition based on the helicopter theory, which suggests positive excess lift. This effect is attributed to unsteady wing-wake interactions. A similar effect is found during hovering.

## Introduction

The aerodynamic forces of an air vehicle or an animal may be affected by the ground proximity. This phenomenon, known as the ground effect, has been extensively studied for aircraft [1] and rotorcraft [2]. Although the effect varies depending on many design parameters, the

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Kolomenskiy D, Maeda M, Engels T, Liu H, Schneider K, Nave J-C (2016) Aerodynamic Ground Effect in Fruitfly Sized Insect Takeoff. PLoS ONE 11 (3): e0152072. doi:10.1371/journal.pone.0152072

**Editor:** Roi Gurka, Coastal Carolina University, UNITED STATES

**Received:** November 18, 2015

**Accepted:** March 7, 2016

**Published:** March 28, 2016

**Copyright:** © 2016 Kolomenskiy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Datasets from this study can be downloaded from http://osf.io/a4thw (Open Science Framework is listed as a recommended repository, DOI will be generated after acceptance of the manuscript for publication to make sure that no further changes are required). The Navier-Stokes solver employed is an open source software available on https://github.com/pseudospectators/FLUSI.git The DOI/Accession numbers associated with the data are as follows: DOI 10.17605/OSF.IO/WURBC ARK c7605/osf.io/wurbc.

**Funding:** DK gratefully acknowledges the financial support from the CRM-ISM (Centre de Recherches Mathematiques - Institut des Sciences
Mathematiques) Postdoctoral Fellowship and from the JSPS (Japan Society for the Promotion of Science) Postdoctoral Fellowship P15061. TE and KS thank the Deutsch-Französische Hochschule / Université Franco-Allemande (DFH-UFA) for financial support. HL was partly supported by the Grant-in-Aid for Scientific Research on Innovative Areas of No. 24120007, JSPS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

general trend is an increase in lift and pitching moment, and a decrease in drag. The effect decays as the distance from the ground increases, and vanishes at a distance slightly larger than the characteristic length of the vehicle. For example, for a hovering helicopter, the excess thrust vanishes if the distance to the ground exceeds 1.25 times the diameter of the main rotor [2].

Rayner [3] proposed a fixed wing lifting line theory for forward flight of birds, bats and insects. His analysis suggested that flight in ground effect provides performance improvements, if the flight speed is not too low. However, this theory could not be applied to hovering or slow forward flight at very low height, since it neglected flapping motion. Normal hovering in ground effect was considered by Gao and Lu [4]. They carried out two-dimensional numerical simulations of hovering and identified three regimes: force enhancement, force reduction, and force recovery, depending on the distance from the ground. Liu et al. [5] considered clap-and-fling near the ground and found force enhancement at all distances. A three-dimensional numerical simulation of fruitfly hovering was carried out by Maeda and Liu [6]. An increase in lift and a reduction in power was found. A significant vertical force was generated on the insect's body due to the 'fountain effect'. Energetic savings have also been reported for a hummingbird hovering in ground effect [7].

Several studies considered pitching-plunging foils near a solid wall or a free surface [8–10]. This configuration is relevant to fish swimming as well as forward flapping flight. The ground effect mainly consists in enhanced propulsive force. However, it also generates a non-zero vertical force due to asymmetry.

The main motivation for this study comes from the fact that the ground proximity is natural for takeoff and landing. These manoeuvres, unlike hovering or forward flight, are characterized by gradual change of distance to the ground. The 'dynamic' ground effect in these circumstances may be different from the 'static' effect at a constant distance [11]. This difference may be even larger for flapping wings than for fixed wings, because animals vary their wing kinematics during takeoff.

So far, the ground effect during takeoff has been assessed for very few insects only. It was found negligible for butterflies (*Pieris rapae* [12], *Papilio xuthus* [13]), a dronefly (*Eristalis tenax*) [14], and a fruitfly (*Drosophila virilis*) [15], but significant for a beetle (*Trypoxylus dichotomus*) [16]. The disparity can be attributed to significant differences in the size, morphology and kinematics of these insects. Thus, our work is motivated by the apparently contradictory conclusions on the significance of the ground effect that could be found in the animal flight literature. It is important to identify the parameters that make the ground effect strong or negligible.

In the present study, we consider a numerical model having the morphology of a fruitfly, with variable wing kinematics and leg parameters. Our objective is to determine if the ground effect can be significant for this model, and which conditions can lead to it. We thus explore the parameter space of the model and perform numerous numerical simulations using FLUSI [17], which is an open source software available on https://github.com/pseudospectators/FLUSI/tree/plos_one__ground__effect. First, for completeness, we revisit the voluntary takeoff of a fruitfly analyzed in [15]. The main difference with respect to [15] is the use of a flight dynamics solver. We then compare takeoffs with modified parameters of the leg thrust model and wing kinematics. Finally, we consider hovering as a limiting case of very slow takeoff.

The paper is organized as follows. In section, we describe our computational approach and the takeoff parameters used in this study. The results are presented in section Leg thrust, first for a natural voluntary takeoff, then for modified takeoffs and for hovering flight. The main conclusions are summarized in section Conclusions.
# Methods

## Morphology and kinematics

In this work, we consider a fruitfly having mass *m* = 1.2 mg and wing length *R* = 2.83 mm, which are the values reported by Chen and Sun [15]. The body is modelled as a rigid solid, and the wings are modelled as rigid flat plates. This approximation is accurate for *Drosophila* during voluntary takeoff [15, 18], though it occasionally fails during fast escape manoeuvres [18]. The wing contour used in this study is shown in Fig 1a. It is adapted from [15]. Its mean chord length is equal to *c* = 0.85 mm. The body is generated by sweeping a circular section of variable radius along a curvilinear centreline (an arc). The body has approximately the same dimensions as in [15]. The side view of the body is shown in Fig 1b. Even though the yaw and roll angles can eventually become large during takeoff, there is no significant trend for all takeoffs. Hence, to simplify the problem, we assume bilateral symmetry. Therefore, the body orientation is fully defined by the pitch angle *β* between the body and the horizontal axis, see Fig 1c. The wing kinematics is described by three angles: *φ*, *α* and *θ*, measured with respect to the stroke plane, as shown in Fig 1d. The positional angle *φ* defines the motion of the wing tip projection on the stroke plane. The deviation (elevation) angle *θ* defines the deviation of the wing tip from the stroke plane. The feathering angle *α* defines the rotation about the longitudinal axis of the wing, and it is related to the geometrical angle of attack (*AoA*) as *α* = 90° – *AoA* during down-stroke and as *α* = 90° + *AoA* during upstroke. It is convenient to refer to an ‘anatomical’ stroke plane angle *η*, i.e., to assume that the inclination of the stroke plane against the body axis is held at a constant angle for any motion of the body.

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Schematic drawing of the morphological model. (a) Wing contour. Coordinates are normalized to the wing length *R*. (b) The body is generated by circular sections of variable radius, which changes depending on the position along the centre line (dash-dotted arc). The body axis *O<sub>b</sub>*x<sub>b</sub>*y<sub>b</sub>*z<sub>b</sub> is the thorax-abdomen principal axis, approximately. The coordinates of the wing pivot points in the body frame of reference *O<sub>b</sub>*x<sub>b</sub>*y<sub>b</sub>*z<sub>b</sub> are (–0.07*R*, ±0.18*R*, 0.115*R). Note that, since *β*(*t*) is prescribed in the simulations and only (*x*<sub>c</sub>, *z*<sub>c</sub>) are dynamically calculated, the position of the body point of reference with respect to the body contour is chosen arbitrarily. (c) The insect’s position with respect to the ground is described by the body point of reference coordinates (*x*<sub>c</sub>, *z*<sub>c</sub>) and the position angle *β*. (d) Definition of the wing’s angles with respect to the stroke plane frame of reference *O<sub>sp</sub>*x<sub>sp</sub>*y<sub>sp</sub>*z<sub>sp</sub>. The origin *O<sub>sp</sub>* is the wing pivot point.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g001
Since the main focus of this study is the ground effect, it is important to ensure that the time evolution of the distance to the ground is consistent with the forces acting on the insect. For this reason, in our computations, unlike in [15], the position of the insect is dynamically computed as opposed to be prescribed. We compute the position of the body point of reference ($x_c, z_c$), see Fig 1b, from Newton's 2nd Law,

$$m \frac{d^2 x_c}{dt^2} = F_{ax} + F_{\ell x}, \quad m \frac{d^2 z_c}{dt^2} = F_{az} + F_{\ell z} - mg, \quad (1)$$

where ($F_{ax}, F_{az}$) is the aerodynamic force, ($F_{\ell x}, F_{\ell z}$) is the leg thrust, subscripts x and z correspond, respectively, to the horizontal and vertical components, m is the insect's mass and g is the gravitational acceleration. Eq (1) are integrated using the adaptive second order Adams-Bashforth scheme [19], simultaneously with the incompressible Navier-Stokes equations. We defined the positive z direction to be upwards and the positive x direction to be forwards (see Fig 1c).

## Aerodynamics

The aerodynamic forces $F_{ax}$ and $F_{az}$ are obtained by solving the three-dimensional incompressible Navier-Stokes equations. The no-slip boundary condition at the body and wings surfaces is imposed using the volume penalization method [20], and the penalized equations are solved using a classical Fourier pseudo-spectral method. More details about the solver and the generic insect model, including a numerical validation case of fruitfly hovering, can be found in [17]. Numerical validation of the ground plane modelling using the volume penalization method is described in S1 Appendix.

The computational domain in the present study is a rectangular box with sides $L_x, L_y$ and $L_z$. Suitable values of $L_x, L_y$ and $L_z$, in terms of accuracy and computational efficiency, depend on the motion of the insect within the domain. Therefore, different values are used in different simulations, as described later in the text. The domain is discretized using a uniform Cartesian grid. Periodic boundary conditions are applied on all sides of the domain, as required by the Fourier discretization. Vorticity sponge boundary conditions are imposed at the left, right, rear and front sides of the domain, as explained in [17, 21], in order to minimize the effect of the finite domain size. The ground surface is modelled as a solid layer at the bottom of the domain, which by periodicity also imposes the no-slip on the top of the fluid domain. We have carried out numerical experiments to ensure that, in the numerical simulations presented in this paper, the domain size is sufficiently large, i.e., its further increase does not change significantly the forces. The dimensions that we chose are also comparable with the size of the mineral oil tank used in the experiments with a mechanical model [22].

## Leg thrust

The model of the leg thrust employed in the present study is a slight modification of the compression spring model proposed in [12]. We assume that takeoff begins from rest and starts at time $t = t_\ell$, which can be estimated from the initiation of the legs motion in the video sequences shown in [15]. The two components of the force are given by

$$F_{\ell x} = Z \cot \phi_{\ell}, \quad F_{\ell z} = Z. \quad (2)$$

The magnitude of the leg force is assumed to depend on the vertical component of the leg extension $\zeta = z_c(t) - z_c(t_\ell)$ only. The force is supposed to be distributed between the three pairs
of legs such that its change with horizontal displacement can be neglected,

$$Z = \begin{cases} (L_{\ell} - \zeta)K_{\ell} & \text{for } \zeta < L_{\ell}, \\ 0 & \text{for } \zeta \ge L_{\ell}, \end{cases} \qquad (3)$$

where $L_\ell$ is the maximum leg extension length, i.e., the difference between the values of $z_c$ when the legs are fully extended at takeoff and when the insect is at rest. When the legs are fully extended, $\zeta = L_\ell$, the legs lose contact with the ground and the force drops to zero. This length is estimated using video sequences in [15] to be equal to $L_\ell = 1.24$ mm. The spring stiffness $K_\ell$ varies in time: it increases from $K_\ell^-$ before takeoff to $K_\ell^+$ after takeoff. The initial value $K_\ell^- = mg/L_\ell$ ensures that the insect is in equilibrium before takeoff, when the aerodynamic force is zero. The final value $K_\ell^+$ is a parameter of the model that controls the maximum leg thrust. Its value can be estimated from the climb velocity at the beginning of takeoff, shown in, e.g., [15]. It may also be estimated from jumps of wingless flies [23, 24] for a slightly different fruitfly, *D. melanogaster*. We assume the time evolution of $K_\ell$ of the form

$$K_{\ell} = \begin{cases} K_{\ell}^{-} & \text{for } t < t_{\ell}, \\ K_{\ell}^{-} + \frac{K_{\ell}^{+} - K_{\ell}^{-}}{\tau_{\ell}}(t - t_{\ell}) & \text{for } t_{\ell} \le t < t_{\ell} + \tau_{\ell}, \\ K_{\ell}^{+} & \text{for } t \ge t_{\ell} + \tau_{\ell}. \end{cases} \qquad (4)$$

The transition time $\tau_\ell$ can be equal to zero, in which case the leg force increases impulsively at the beginning of takeoff. However, measurements of the leg force [23] suggest a gradual increase which can be accounted for by setting $\tau_\ell$ to a value larger than zero. The value $\tau_\ell = 1.3$ms results in the gradient $dF_\ell/dt$ consistent with the experimental data shown in [23]. The direction $\phi_\ell$ also changes in time. Before takeoff, when the insect is at rest, the force is applied only in the vertical direction, i.e., $\phi_\ell^- = 90^\circ$. During takeoff, the horizontal component is non-zero, in general. We assume a time evolution of the form

$$\phi_{\ell} = \begin{cases} \phi_{\ell}^{-} & \text{for } t < t_{\ell}, \\ \phi_{\ell}^{-} + \frac{\phi_{\ell}^{+} - \phi_{\ell}^{-}}{\tau_{\ell}}(t - t_{\ell}) & \text{for } t_{\ell} \le t < t_{\ell} + \tau_{\ell}, \\ \phi_{\ell}^{+} & \text{for } t \ge t_{\ell} + \tau_{\ell}. \end{cases} \qquad (5)$$

The values of the leg thrust model parameters used in our numerical simulations are given in Table 1.

<table><caption>Table 1. Kinematics and leg model parameters of the takeoffs considered in the present study.</caption><thead><tr><th rowspan="2">Case<br>Name</th><th colspan="4">Kinematics</th><th colspan="5">Leg model</th></tr><tr><th>φ, α, θ</th><th>β</th><th>η, °</th><th>z<sub>c</sub>(0), mm</th><th>L<sub>ℓ</sub>, mm</th><th>K<sub>ℓ</sub><sup>+</sup>, N/m</th><th>φ<sub>ℓ</sub><sup>+</sup>, °</th><th>t<sub>ℓ</sub>, ms</th><th>τ<sub>ℓ</sub>, ms</th></tr></thead><tbody><tr><td>Voluntary</td><td>Fig 2</td><td>Fig 2</td><td>62</td><td>1.08</td><td>1.24</td><td>0.165</td><td>84</td><td>4.2</td><td>1.3</td></tr><tr><td>Slow</td><td>Fig 2</td><td>Fig 2</td><td>62</td><td>1.08</td><td>1.24</td><td>0.041</td><td>84</td><td>4.2</td><td>1.3</td></tr><tr><td>Simplified</td><td>Fig 7a</td><td>46.3°</td><td>32</td><td>3.11</td><td>1.24</td><td>0.0095...0.043</td><td>84</td><td>0</td><td>1.3</td></tr><tr><td>Hovering</td><td>Fig 7a and b</td><td>55°</td><td>55</td><td>2.07</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td></tr></tbody></table>

doi:10.1371/journal.pone.0152072.t001
# Results and Discussion

## Summary of the numerical simulations

The starting point for our study is the voluntary takeoff, as it is shown in section Voluntary takeoff (in agreement with [15]) that the ground effect is very small in that case. It is much smaller than during hovering (cf. [6]). We conjecture that this difference is due to the large takeoff vertical velocity, which is mainly the result of the leg thrust. To test this hypothesis, in section Slow takeoff, we discuss a situation in which the legs produce less force and the insect takes off slower. The ground effect becomes significant. The vertical force increases during the first two wingbeats due to the ground effect, but slightly decreases later on. We then carry out a parametric study using periodic wing kinematics in section Takeoffs with simplified kinematics, and find an even stronger adverse ground effect. Finally, in section Ground effect in hovering flight, we find similar trends during the first wingbeats in hovering flight, which can is considered as a limiting case of takeoff with zero rate of climb. Tables 1 and 2 summarize the parameters of the different cases considered in the present study. Datasets for the ‘voluntary’ and ‘simplified’ cases can be downloaded from [25].

## Voluntary takeoff

In this section, we consider voluntary takeoff of a fruitfly with the parameters as in the first lines in Tables 1 and 2. This case shows some important general features of fruitfly takeoff such as the first wingbeat cycles beginning while the legs extend. Therefore it is likely that, despite some variability in voluntary takeoffs, the ground effect in general remains of the same order of magnitude in natural circumstances.

The values of the body and wing angles are taken from one of the cases documented in [15]. However, the wing motion in [15] is not exactly symmetric. Therefore, the time series of $\phi$, $\alpha$ and $\theta$ that we use for both wings correspond to the left wing data shown in [15]. Fig 2 presents the time evolution of the wing positional angle $\phi(t)$, the feathering angle $\alpha(t)$, the elevation angle $\theta(t)$ and the body pitch angle $\beta(t)$, which are prescribed in our numerical simulations. The angle between the horizontal plane and the stroke plane, $\eta - \beta$, is also shown for reference.

Even though the wing motion is not exactly periodic, it is useful to introduce the wing beat frequency. When calculated using the average wing beat cycle period over the five cycles shown in Fig 2, it is equal to $f = 169$ Hz. Similarly, the average wing beat amplitude is equal to $\Phi = 134^{\circ}$, and the characteristic wing tip velocity is $U = 2\Phi Rf = 2.23$ m/s. The kinematic viscosity of air, equal to $\nu = 1.45 \cdot 10^{-5}$ m$^2$/s yields the Reynolds number $Re = Uc/\nu = 131$. Note that $U$ and $Re$ do not account for the forward speed of the body.

The computational domain size is equal to $L_x = L_y = 5R$, $L_z = 8R$, where $R$ is the wing length. The influence of the domain size in the vertical direction is discussed in S2 Appendix. The number of grid points in each direction, respectively, is $N_x = N_y = 640$ and $N_z = 1280$. The penalization parameter is $\epsilon = 2.5 \cdot 10^{-4}$ (for details see, e.g., [21]).

<table><caption>Table 2. Numerical parameters of the takeoffs considered in the present study.</caption><thead><tr><th rowspan="2">Case</th><th colspan="6">Numerical parameters</th></tr><tr><th>L<sub>x</sub></th><th>L<sub>y</sub></th><th>L<sub>z</sub></th><th>N<sub>x</sub></th><th>N<sub>y</sub></th><th>N<sub>z</sub></th><th>ε</th></tr></thead><tbody><tr><td>Voluntary</td><td>5R</td><td>5R</td><td>8R</td><td>640</td><td>640</td><td>1280</td><td>2.5 · 10<sup>-4</sup></td></tr><tr><td>Slow</td><td>5R</td><td>5R</td><td>6R</td><td>640</td><td>640</td><td>768</td><td>2.5 · 10<sup>-4</sup></td></tr><tr><td>Simplified</td><td>4R</td><td>4R</td><td>6R</td><td>512</td><td>512</td><td>768</td><td>2.5 · 10<sup>-4</sup></td></tr><tr><td>Hovering</td><td>8R</td><td>8R</td><td>4R</td><td>864</td><td>864</td><td>432</td><td>2.5 · 10<sup>-4</sup></td></tr></tbody></table>

doi:10.1371/journal.pone.0152072.t002
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Time evolution of the angular position of the body and of the wings during the voluntary takeoff. Gray shaded regions correspond to downstrokes. $\eta - \beta$ is the angle between the horizontal plane and the stroke plane, i.e., the global stroke plane angle [26].</figcaption></figure>

doi:10.1371/journal.pone.0152072.g002

The aerodynamic ground effect is evaluated by comparing two numerical simulations with two different values of the initial distance from the body point of reference to the ground: $z_c(0) = 0.38R$ and $2R$, which we denote ‘in ground effect’ (IGE) and ‘out of ground effect’ (OGE), respectively. The first case corresponds to a takeoff from a flat ground surface, with $z_c(0)$ being consistent with the data in [15]. In the second case, the leg model behaves as during takeoff from the ground, but the aerodynamic interaction between the insect and the ground is weak because of the large distance. This case may be interpreted as takeoff from a perch that provides enough support for the legs but has a small surface, such that the aerodynamic interactions are negligible. With the distance equal to $2R$ or greater, the ground effect is negligible during hovering [13]. The circulation of the wake vortices is mainly determined by the integral aerodynamic force, therefore it is not larger during takeoff than during hovering, and the spatial rate of decay of the induced velocity is the same. Hence, the ground effect with the distance equal to $2R$ is likely to be negligible during takeoff. The influence of the ground on the shape of the vortices is only visible during the 2nd wingbeat and later on. This difference is localized to the vicinity of the ground plane. Since the insect is relatively far from the ground by that time, this difference is unlikely to have any influence on the aerodynamic forces.

Fig 3a shows the fruitfly model and the wake, IGE and OGE, at 4 subsequent time instants. The vortices created by the wings and the body are identified as the volume of fluid enclosed by the iso-surfaces of the Q-criterion. At $t = 0$, the air is at rest. The insect body is almost horizontal. The wings are in a pre-takeoff position from which they begin the first downstroke after $t = 4.1$ ms. The time $t = 9.2$ ms corresponds to the first reversal from downstroke to upstroke. Because of the small body pitch angle $\beta$, the stroke plane is effectively vertical. In addition, the wing tip speed during the first downstroke is smaller than during all subsequent strokes. Therefore, the vertical aerodynamic force is small, but the body lifts noticeably because of the leg thrust. The time $t = 12.8$ ms corresponds to the second upstroke. At this point, the distance from the body point of reference to the ground $z_c$ is already larger than the wing length $R$. Therefore, the aerodynamic interference with the ground is expected to be very small. Note that the kinematics during the first two wingbeat cycles are a transient. After that, the time evolution of the wing angles approaches a periodic regime and the stroke plane becomes less inclined with respect to the ground, see Fig 2.
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig 3. Voluntary takeoff.** (a) Visualization of the wings, body and ground surface, and the wake at 4 subsequent time instants. Blue semi-transparent iso-surfaces show the Q-criterion, $Q/f^2 = 15$. (b) Vertical and horizontal displacement. To obtain distance $z_c$ from the ground for the IGE case, add 1.08 mm. The black dash-dotted line indicates $z_c = R$. (c) Components of the leg force. (d) horizontal and (e) vertical components of the aerodynamic force and (f) the aerodynamic power. The black dash-dotted line in panel e indicates the weight. Solid circles connected by dotted lines show wingbeat cycle averages. The results for OGE and IGE are shown, but the curves in panels b to f overlap because the difference is negligible.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g003
The displacement of the body point of reference is shown in Fig 3b. It presents the evolution of the vertical component $\zeta(t) = z_c(t) - z_c(0)$ and the horizontal component $\xi(t) = x_c(t) - x_c(0)$ over time for the cases IGE and OGE. The curves overlap. In both cases, at the end of the 4th wingbeat cycle, $t = 28.4$ ms, the insect gains 8.7 mm of altitude and propels 1.6 mm forward. These numbers are consistent with the trajectories shown in [15]. The displacements IGE and OGE differ by less than 1%. Therefore, the ground effect on $\zeta$ and $\xi$ is indeed negligible.

Fig 3c shows the two components of the leg force. At $t = 0$, the vertical component of the leg force is equal to the weight and the horizontal component is zero. The jump is triggered at $t_e = 4.2$ ms. At time $t_e + \tau_e = 5.5$ ms, both components reach their peaks. After that the force decreases and vanishes at $t = 9.3$ ms, when the legs lose contact with the ground. Note that the leg thrust can, in principle, be different for the takeoffs IGE and OGE, because the leg model depends on the aerodynamic force via $z_c(t)$. However, for the voluntary takeoff considered here, there is no influence of the ground effect.

The horizontal and the vertical components of the aerodynamic force are shown in Fig 3d and 3e, respectively. Over the first four wingbeat cycles, the wingbeat averaged aerodynamic forces are significantly lower than the weight. This can be explained by the large initial rate of climb due to the leg thrust, which cannot be supported by the wings. Even during the fourth wingbeat, the wing force is equal to 29% of the weight. The vertical acceleration is therefore negative after the legs lose contact with the ground, and the rate of climb slowly decreases. The ground effect is, again, negligible. Even during the first wingbeat cycle, when the wings approach the ground surface, the difference in the instantaneous vertical force between IGE and OGE is at most 0.0005 mN, i.e., about 4% of the weight. The wingbeat cycle averaged forces differ by less than 1% of the weight.

Fig 3f displays the time evolution of the aerodynamic power, when operating IGE and OGE. Note that, in this study, we do not consider the inertial power because the wings have the same kinematics in both cases, IGE and OGE. Therefore, the inertial power is the same. The aerodynamic power is the aerodynamic component of the power required to actuate the wings,

$$P = -M_l \cdot (\Omega_l - \Omega_b) - M_r \cdot (\Omega_r - \Omega_b). \qquad (6)$$

In Eq (6), $M_l$ and $M_r$ are the aerodynamic moments of the left and of the right wing, respectively, relative to the corresponding pivot point. $\Omega_l$ and $\Omega_r$ are the angular velocities of the wings and $\Omega_b$ is the angular velocity of the body. All vectors are taken in the laboratory frame of reference. P is positive if power is consumed. We find that it is positive during most part of the takeoff (see Fig 3f). Only at the reversals during the first two cycles, when the body velocity is still small, P is slightly negative. During the 2nd wingbeat, the mean body-mass specific aerodynamic power is equal to $P_b^* = P_{ave}/m = 21$ W/kg. Assuming that the muscles contribute to 30% of the body mass, the mean muscle-mass specific aerodynamic power is equal to $P_m^* = P_{ave}/(0.3m) = 69$ W/kg. The relative difference in the cycle averaged values between IGE and OGE is less than 0.5%.

We conclude that the ground effect is unimportant for the voluntary takeoff, a result which is in agreement with [15]. This is mainly a consequence of rapid acceleration during the first wingbeat cycle, when the legs produce a large vertical force. The main question of the next section is whether this scenario changes if the takeoff is slower and the insect remains near the ground for a longer time. The rate of climb at the beginning of takeoff is controlled by the leg model stiffness coefficient $K_l^+$, and the horizontal velocity is controlled by the leg angle $\phi_l^+$.

## Slow takeoff

This section describes a modified takeoff with the leg thrust coefficient decreased to $K_l^+ = 0.041$ N/m (see the second line in Table 1). Smaller $K_l^+$ results in less leg thrust and slower
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Slow takeoff. (a) Vertical and horizontal displacement. Dash-dotted line indicates $z_c = R$. (b) Components of the leg force. (c) horizontal and (d) vertical components of the aerodynamic force. The black dash-dotted line in panel d indicates the weight. Solid circles connected by dotted lines show wingbeat cycle averages. Note that, in panels b and c, the lines for IGE and OGE almost coincide.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g004

climb, compared to the natural voluntary takeoff. Therefore we refer to this case as a 'slow takeoff'. Any further decrease of $K_l^+$ would result in a very close approach of the insect to the ground surface and, ultimately, to a collision which would require special treatment. Modelling of such collisions could be an interesting topic for future research (see [27] for a review of structural modelling of insect wings, including impact modelling). For a fruitfly, collisions between the wings and the ground may be undesirable because of the large wingbeat frequency and light wing structure. Therefore, $K_l^+ = 0.041$ N/m is an interesting limiting case.

In the present numerical simulations, the computational domain size in x and y directions, the discretization grid step size and the penalization parameter are the same as in the previous section. The domain size in the vertical direction z is reduced to $6R$ because the insect gains much less altitude by the end of the simulation.

Fig 4a shows the displacement of the body point of reference. The rate of climb is about one third of its original value and the insect only gains 3.9 mm by the end of the 4th wingbeat cycle. This is just slightly larger than the wing length $R$ (2.83 mm). The displacement is slightly larger for IGE than for OGE in both directions, horizontal and vertical. The time evolution of leg thrust is given in Fig 4b. There is no visible difference between the two cases. The peak of the vertical force is equal to 0.051 mN, which is about four times less than in the original voluntary takeoff discussed in section Voluntary takeoff.

The time evolution of the instantaneous aerodynamic force, shown in Fig 4c and 4d, is qualitatively similar to the voluntary takeoff case considered previously. The difference between the cases OGE and IGE is negligible for the horizontal force (Fig 4c), but for the vertical force it reaches values as large as 0.0027 mN, i.e., 23% of the weight (Fig 4d). Fig 5 shows the difference between the wingbeat averaged forces in the cases IGE and OGE, normalized by the weight. The vertical force difference is shown in Fig 5b. During the 1st wingbeat, the ground effect makes the total vertical force increase by almost 6% of the weight (red line). However, during
<figure><img src="image_9.png" /><figcaption>**Fig 5. Slow and voluntary takeoffs.** The difference between the cases IGE and OGE, in terms of the wingbeat cycle averaged aerodynamic force normalized by the body weight (a) horizontal and (b) vertical components. The forces acting on the wings and the body are shown separately. The total force, which is their sum, is also shown.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g005

the 2nd wingbeat, the extra force decreases to only 2% of the weight. During the 3rd, the 4th and the 5th wingbeats, the difference between the vertical forces IGE and OGE is very small and negative. The increase of the vertical force during the first wingbeat is mainly due to the wings (blue line). The extra force acting on the body is only about 1% of the weight. However, later on, the contribution of the body becomes important because it remains positive, whereas for the wings it becomes negative. The horizontal force difference, shown in Fig 5a, is positive, i.e., the propulsive force increases due to the ground effect by about 2% of the weight, for all wingbeats. The contribution of the body is up to 1% of the weight. For reference, Fig 5 also shows the force differences during the voluntary takeoff. They are all smaller than 1%.

The aerodynamic power, in the cases IGE and OGE, is compared in Fig 6. The maximum difference is of about 3% in magnitude for the slow takeoff, but less than 1% for the voluntary takeoff. Considering the slow takeoff, the insect consumes more power when operating in ground effect (IGE) during the first two wingbeat cycles, but less power during the subsequent cycles. Overall, we find that the differences in the power are small.

<figure><img src="image_10.png" /><figcaption>**Fig 6. Slow and voluntary takeoffs.** Aerodynamic power ratio IGE/OGE.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g006
<figure><img src="image_11.png" /><figcaption>**Fig 7. Periodic wing kinematics.** (a) Data obtained by periodization of the last wingbeat shown in [15]. It is used in section Takeoffs with simplified kinematics. Also, in section Ground effect in hovering flight, this kinematics is referred to as 'P1' ('P' for 'Periodic'). The cycle begins from the upstroke. (b) Data adapted from [6]. In section Ground effect in hovering flight, it is referred to as 'P2'. The cycle begins from the downstroke.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g007

## Takeoffs with simplified kinematics

In the previous sections we noticed that the ground effect depends on the takeoff kinematics. We are mainly interested in the effects that might be generally applicable to a fruitfly sized insect takeoff. Therefore, in this section we consider parametric studies. They are performed using simplified periodic wing kinematics. The time evolution of the three wing angles over one wing-beat period, obtained by periodization of the last wingbeat in [15], is shown in Fig 7a. The wing-beat frequency is equal to $f = 210$ Hz. This takeoff mode can be relevant to MAVs, for which the wing kinematics and the body angle during takeoff do not change as much as for the fruitfly (see, e.g., [28]). The leg strength parameter $K_l^+$ is varied, resulting in a variation of the takeoff rate of climb $V_{t.o.}$. The body angle is constant and equal to $\beta = 46.3^\circ$, the anatomical stroke plane angle is equal to $\eta = 32^\circ$. In these computations, we use $L_x = L_y = 4R$, $L_z = 6R$, $N_x = N_y = 512$ and $N_z = 768$, corresponding to more than 200 million grid points. The penalization parameter is equal to $\epsilon = 2.5 \cdot 10^{-4}$.

Smaller $K_l^+$ implies smaller rate of climb (Fig 8c) which leads to a more significant ground effect (Fig 8d). A striking feature of Fig 8b is a significant decrease of the vertical force during the 4th, 5th and 6th wingbeats, by up to 6%. The horizontal force varies slightly, by about 1% (see Fig 8a).

The decrease in vertical force on the wings found in the IGE cases compared to OGE (Fig 8b) is an adverse ground effect which may be the result of complex wing-wake-ground interactions that depend on wing and body kinematics. Adverse ground effects have been reported previously for fixed-wing aircraft [11], but they consist in increased drag together with increased lift. However, for a two-dimensional ellipse with normal hovering kinematics [4], the mean vertical force decreases when the height from the ground is between approximately 1.5D and 4D, where D is the chord length of the ellipse. For flapping wings, an adverse ground effect was found by Quinn et al. [9]. They considered an airfoil undergoing pitch oscillations in a closed-loop water channel with prescribed free-stream velocity. Such a configuration represents a section of a bird wing in forward flight or a fish fin. Even if the pitching motion was symmetric, the proximity of the ground broke the symmetry of the flow. Thus, the airfoil produced non-zero lift. The lift was positive if the distance to the ground was less than 40% of the chord length, but it became negative at larger distances, such that the lift force pulled the airfoil towards the ground. Nevertheless, the extra propulsive force due to the ground effect was positive in all cases. Note, however, that the flows considered in [4] and [9] are effectively two-dimensional.
<figure><img src="image_12.png" /><figcaption>**Fig 8. Takeoffs with simplified kinematics.** (a) Horizontal and (b) vertical difference between the wingbeat-averaged force IGE and OGE, normalized to the insect weight. (c) Vertical velocity of the body point of reference (rate of climb) versus time. (d) Maximum normalized force difference versus takeoff rate of climb at the moment when the legs lose contact with the ground.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g008

In the present work, importantly, we find an adverse ground effect in a three-dimensional configuration, which has not been previously recognized. The wingbeat cycle averaged vertical force of the wings in the ground effect is slightly larger during the first two cycles, but, as the insect flies away from the ground, the vertical force in the case IGE becomes less than that in the case OGE.

Fig 9 shows a comparison of the wake at the end of the 5th wingbeat cycle ($t = 23.78$ ms) in the 4 different takeoffs IGE with different values of $K_l^+$. It corresponds to the maximum decrease of the vertical force. There are noticeable differences between the vortices when the takeoff is fast and when it is slow. The part of the wake that approaches the ground deforms when it impinges on the ground. It then rolls up in a pair of vortex rings. Similar "ground vortices" are known in the context of helicopter rotor aerodynamics. In each of the 4 cases shown in Fig 9, they have different strength and position with respect to the wings. Therefore, they induce the downwash of different strength.

A detailed view of the flow near the wings is presented in Fig 10. It shows the pressure and the vorticity magnitude during the 6th wingbeat for two IGE cases with different $K_l^+$. Five time instants are visualized in five rows, respectively. The left column shows the pressure distribution over the surface of the insect, as well as over a hemisphere of radius 0.9R centred the body reference point, for $K_l^+ = 0.0095$ N/m. Fig 10a is at $t = 24.73$ ms, during upstroke. The dark blue area near the leading edge that expands towards the wing tip is the trace of the leading-edge vortex (LEV), similar to the one discussed in [29]. The LEV at the downstroke is evident in Fig 10d at $t = 27.58$ ms. The pressure distributions during the reversals are more complex (Fig 10b, 10c and 10e).
<figure><img src="image_13.png" /><figcaption>**Fig 9.** Flow visualization at the end of the 5th wingbeat (t = 23.78 ms) for the simplified kinematics cases. Iso-surfaces Q/f² = 15 are shown for 4 different takeoffs IGE: (a) Kᵢ⁺ = 0.0095 N/m; (b) Kᵢ⁺ = 0.0184 N/m; (c) Kᵢ⁺ = 0.0307 N/m; (d) Kᵢ⁺ = 0.0430 N/m.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g009

The middle column of Fig 10 shows the pressure iso-contours sampled on the surface of the semisphere. Two cases are compared, $K_i^+ = 0.0095$ N/m and 0.0430 N/m. This choice of $K_i^+$ corresponds to the largest difference in the vertical force (see Fig 8b). For each of the cases, two isolines are drawn, $\rho/\rho R^2 f^2 = -2$ and -5. There is a large difference between the contours for different $K_i^+$ in the far wake (wake far from the wings). However, in the near wake of a wing including the LEV, the difference is much smaller, which is consistent with the force changing by only a few per cent.

The iso-contours of the vorticity magnitude are compared in the right column. Here again, despite significant differences in the far wake, the LEV contours virtually overlap for $K_i^+ = 0.0095$ N/m and 0.0430 N/m. This shows that the ground effect has almost no influence on the dynamics of the LEV.

## Ground effect in hovering flight

In this section, we simplify the kinematics even further. We consider hovering with the insect body being fixed. The flight dynamics solver is not used in this case. The distance from the body centre to the ground is equal to 0.48*R (where R is the wing length) for hovering in ground effect (IGE) and 2.4*R for hovering out of ground effect (OGE). The body pitch angle and the anatomical stroke plane angle are both constant and equal to 55°, such that the stroke plane is horizontal. The wing kinematics is the same as in the previous section, see Fig 7a. We denote it
<figure><img src="image_14.png" /><figcaption>**Fig 10.** Visualization of the leading-edge vortex during the 6th wingbeat for the simplified kinematics cases, IGE. Left column (a-e) shows the dimensionless pressure distribution over the surface of the model and over a hemisphere of radius 0.9R around the body point of reference. Middle column (f-j) shows the pressure iso-contours for two different takeoffs: $K_i^+ = 0.0095$ N/m and $0.0430$ N/m. Right column (k-o) compares iso-contours of the dimensionless vorticity magnitude for the same two takeoffs. Time instants are $t = 24.73, 25.68, 26.63, 27.58$ and $28.53$ ms ($tf = 5.2, 5.4, 5.6, 5.8$ and $6$, where $f$ is the wingbeat frequency).</figcaption></figure>

doi:10.1371/journal.pone.0152072.g010
<figure><img src="image_15.png" /><figcaption>**Fig 11. Hovering flight.** (a) The ratio of wingbeat averaged vertical force IGE to OGE, *F*<sub>z,ave,IGE</sub>/*F*<sub>z,ave,OGE</sub>. (b) The ratio of wingbeat averaged power IGE to OGE, *P*<sub>ave,IGE</sub>/*P*<sub>ave,OGE</sub>.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g011

as ‘P1’ kinematics. The wingbeat frequency is equal to *f* = 218 Hz. The first wingbeat starts from the upstroke, as done in [15].

In these numerical simulations we are interested in the long-time evolution of the aerodynamic forces, which after the initial transient eventually reach a periodic state. Most of the results known from the helicopter rotor theory [2] are obtained in reference to the periodic state, while the takeoffs considered in the previous sections of this paper (the slow takeoffs, in particular) last only for a few wingbeats. Therefore, it is instructive to consider the time evolution of the aerodynamic forces during hovering from *t* = 0 until the time when the periodic state is reached.

Since the time span of the numerical simulations presented in this section is large, it is necessary to increase the domain size in the horizontal directions. We set *L*<sub>x</sub> × *L*<sub>y</sub> × *L*<sub>z</sub> = 8*R* × 8*R* × 4*R*, where *z* is the vertical direction. The number of grid points is *N*<sub>x</sub> × *N*<sub>y</sub> × *N*<sub>z</sub> = 864 × 864 × 432. The penalization parameter is equal to *ε* = 2.5 · 10<sup>-4</sup>.

The quantity of interest is the ratio of the wingbeat averaged forces, IGE to OGE: *F*<sub>z,ave,IGE</sub>/*F*<sub>z,ave,OGE</sub>. This quantity is shown in Fig 11a. The red solid line with “+” symbols corresponds to the total vertical force ratio. As already noticed in [6], the vertical force during hovering in ground effect, *F*<sub>z,IGE</sub>, reaches its periodic state significantly later than during hovering out of ground effect, *F*<sub>z,OGE</sub>. Therefore, the ratio of their wingbeat averages, *F*<sub>z,ave,IGE</sub>/*F*<sub>z,ave,OGE</sub> converges slowly with the number of wingbeats. It oscillates between 102% and 107%. After 27 wingbeats it reaches 106%.

A similar comparison for the force generated by the wings is shown with a red dot-dashed line in Fig 11a. It was calculated by integration of the distributed forces over the wings only, in the same numerical simulations. Therefore, the aerodynamic interaction between the body and the wings is included. This force ratio drops from 100.6% to 92.5% during the first 6 wingbeats, oscillates and then increases to 96.3%.

The time evolution of the wake vortices generated by the insect is shown in Fig 12. The visualized time instants correspond to the end of the 1st, 2nd, ..., 25th downstroke. There are significant differences between the positions of the vortices over the first four time instants. The first wingbeat generates very strong vortex rings, that collide with the ground. Then they rebound during the second wingbeat, and parts of them moving upwards are still visible during the third wingbeat. The downwash produced by these vortices influences the nearer wake dynamics and it is likely to be responsible for the decrease of the vertical force during the first few wingbeats. After the 10th
<figure><img src="image_16.png" /><figcaption>**Fig 12.** Time evolution of the wake during hovering. Iso-surfaces of the Q-criterion, $Q/f^2 = 15$, are shown at the end of the downstroke. 'P1' kinematics with 100% wingbeat amplitude.</figcaption></figure>

doi:10.1371/journal.pone.0152072.g012

wingbeat, the wake approaches its quasi-periodic state. There are almost no visible differences between the visualizations at the end of the 20th wingbeat and at the end of the 25th wingbeat.

The pair of numerical simulations (cases IGE and OGE) that we have discussed in the above paragraphs leads to the following conclusions.

1. Over the first 27 wingbeats, $F_{z,ave,IGE}/F_{z,ave,OGE}$ varies within about 5% for the total force and 9% for the wings force.

2. The wing generates less vertical force in the case IGE than in the case OGE (adverse ground effect).

3. The body makes an important contribution to the total vertical force when operating IGE, which results in the excess total vertical force (positive ground effect).

These conclusions are, of course, only valid for the particular wing shape and kinematics used in the simulation. Periodic flapping is only an approximation to the real insect wing motion, which varies from one wingbeat to another, and depends on many different conditions. To determine the effect of all existing fruitfly wing kinematics is beyond the reach of our numerical simulations. However, it is useful to compare a few different cases.

We carried out numerical simulations with the wing kinematics used in [6] (abbreviated as 'P2' in the figures). Note that, in this case, the first wingbeat begins from the downstroke, as
shown in Fig 7b. The wingbeat frequency is the same, $f = 218$ Hz. The results of these numerical simulations are shown in Fig 11 with green lines. They are qualitatively similar to the previously shown ‘P1’ case, but the values are systematically larger. $F_{z,ave,IGE}/F_{z,ave,OGE}$ reaches 111% for the total force and 100% for the wings, such that there is no adverse ground effect after the periodic state is established.

The adverse ground effect is rarely encountered in the aircraft or rotorcraft aerodynamics literature. However, in the context of flapping wings, it is not unusual. In the two-dimensional numerical simulations [4], a U-shape profile of the force ratio $F_{z,ave,IGE}/F_{z,ave,OGE}$ versus $h/c$ was found, where $h$ is the distance from the wing centre to the ground and $c$ is the wing chord. $F_{z,ave,IGE}/F_{z,ave,OGE}$ was greater than 100% for $h/c < 1.5$, but less than 100% for $h/c > 1.5$, and the minimum ratio was of about 54%.

In the three-dimensional model considered in the present paper, it is possible to partially reduce the three-dimensional effects by decreasing the wing beat amplitude. One can then expect the adverse ground effect to be amplified. Indeed, this is what we find by decreasing the wingbeat amplitude by a factor of 2 (by rescaling the positional angle shown in Fig 7 such that the amplitude of the positional angle is halved but the mean positional angle is unchanged, and the wingbeat frequency remains unchanged). The results are shown in Fig 11 with blue lines. Now we have $F_{z,ave,IGE}/F_{z,ave,OGE} < 100\%$ at intermediate times, for wings force and for the total force. The final periodic state produces a very slight excess of the total vertical force (less than 1%). Note that, in this reduced-amplitude case, the total force ratio decreases more than the wing force ratio. This indicates, not surprisingly, that the feathering effect becomes weaker when the wingbeat amplitude is reduced. Similar computations with the ‘P2’ kinematics show the same trend with an even larger decrease of $F_{z,ave,IGE}/F_{z,ave,OGE}$.

The wingbeat averaged aerodynamic power ratio $P_{ave,IGE}/P_{ave,OGE}$ is shown in Fig 11b. Its variation is smaller than the variation of the force, and the computations suggest that its long-time limit is between 97% and 99%, in all cases that we have considered. The shape of the time evolution profiles of the power ratio is approximately similar to the time profiles of the wings vertical force ratio. This means that a local decrease of the vertical force ratio is accompanied by a decrease of the power ratio. Therefore, if the kinematics of the wings operating in ground effect is adjusted such that $P_{ave,IGE}/P_{ave,OGE} = 100\%$ at any time, the force ratio $F_{z,ave,IGE}/F_{z,ave,OGE}$ is likely to increase. Among other factors, the feathering angle is very likely to change passively, when in ground effect, due to compliance of the wing [30–32]. Such effects would need further investigation.

## Conclusions

The aerodynamic ground effect in fruitfly sized insect takeoff has been studied numerically using high performance computing. The three-dimensional incompressible Navier-Stokes equations were solved using a pseudo-spectral method with volume penalization [33, 34] using the FLUSI open source code [17], in order to obtain the flow field and the aerodynamic forces acting on the insect. The takeoff trajectories were calculated using a simple flight dynamics solver that accounts for the body weight, inertia, and the legs thrust. A series of computations has been carried out to explore the parametric space of the model. A natural voluntary takeoff of a fruitfly, modified takeoffs with different kinematics and leg model parameters, and hovering flights have been compared.

We found that the ground effect during the natural voluntary takeoff is negligible. The wingbeat averaged forces only differ by less than 1% of the weight. The aerodynamic power differs by less than 0.5%.
In the modified takeoffs, we decreased the leg strength. As a consequence, the rate of climb decreased and the ground effect became significant. Surprisingly, the vertical force did not always increase. It even dropped in some of the cases that we considered. This is an unsteady effect related to the vortex rings bouncing off the ground surface.

To better understand the mechanism of the adverse ground effect, we considered hovering near a flat ground surface, being the limiting case of zero rate of climb. In that case, the fountain effect produced a large upward force on the insect's body. The net ground effect was therefore positive. However, the aerodynamic force acting on the wings in ground effect was sometimes less than when the wings operate out of ground effect. The most significant decrease was observed during the first 15 wingbeats. Note that this is a much longer time period than a typical takeoff. At long time hovering, the effect was either positive or negative, depending on the wings kinematics.

The parameter space in the takeoff problem is very large. In the present study, we focused on the legs thrust and wing kinematics. However, the aerodynamic ground effect may also be sensitive to the Reynolds number, because the structure of the wake at high *Re* is significantly different from that at low *Re*. Since high Reynolds number computations are costly, they are beyond the scope of the present study, but it is an important question for future research.

## Supporting Information

**S1 Appendix. Numerical validation of the ground plane modelling using the volume penalization method.** See file S1_Appendix.pdf.
(PDF)

**S2 Appendix. Influence of the domain size in the vertical direction.** See file S2_Appendix.pdf.
(PDF)

## Acknowledgments

Numerical simulations were carried out using HPC resources of IDRIS, Paris (project 81664) and of Aix-Marseille Université (project Equip@Meso). DK gratefully acknowledges the financial support from the CRM-ISM Postdoctoral Fellowship and from the JSPS Postdoctoral Fellowship. TE and KS thank the Deutsch-Französische Hochschule / Université Franco-Allemande (DFH-UFA) for financial support. HL was partly supported by the Grant-in-Aid for Scientific Research on Innovative Areas of No. 24120007, JSPS.

## Author Contributions

Conceived and designed the experiments: DK MM HL. Performed the experiments: DK TE. Analyzed the data: DK TE MM. Contributed reagents/materials/analysis tools: TE KS JCN. Wrote the paper: DK MM TE HL KS JCN.

## References

1. Asselin M. An introduction to aircraft performance. AIAA education series. American Institute of Aeronautics & Astronautics; 1997.

2. Leishman JG. Principles of helicopter aerodynamics. Cambridge University Press; 2006.

3. Rayner JMV. On the Aerodynamics of Animal Flight in Ground Effect. Phil Trans R Soc B. 1991; 334 (1269):119-128. doi: 10.1098/rstb.1991.0101

4. Gao T, Lu X. Insect normal hovering flight in ground effect. Phys Fluids. 2008; 20:087101. doi: 10.1063/1.2958318
5. Liu Y, Liu N, Lu X. Numerical study of two-winged insect hovering flight. Adv Appl Math Mech. 2009; 1 (4):481-509.

6. Maeda M, Liu H. Ground Effect in Fruit Fly Hovering: A Three-Dimensional Computational Study. J Bio-mech Sci Eng. 2013; 8(4):344-355. doi: 10.1299/jbse.8.344

7. Kim EJ, Wolf M, Ortega-Jimenez VM, Cheng SH, Dudley R. Hovering performance of Anna's hummingbirds (*Calypte anna*) in ground effect. J R Soc Interface. 2014; 11(98):20140505. doi: 10.1098/rsif.2014.0505 PMID: 24990291

8. Tanida Y. Ground Effect in Flight. JSME Int J B. 2001; 44(4):481-486. doi: 10.1299/jsmeb.44.481

9. Quinn DB, Moored KW, Dewey PA, Smits AJ. Unsteady propulsion near a solid boundary. J Fluid Mech. 2014; 742:152-170. doi: 10.1017/jfm.2013.659

10. Wu J, Shu C, Zhao N, Yan W. Fluid Dynamics of Flapping Insect Wing in Ground Effect. J Bionic Eng. 2014; 11:52-60. doi: 10.1016/S1672-6529(14)60019-6

11. Jones BL. Experimental investigation into the aerodynamic ground effect of a tailless chevron-shaped UCAV [Ph.D. thesis]. Department of the Air Force, Air University, Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio, USA; 2005.

12. Bimbard G, Kolomenskiy D, Bouteleux O, Casas J, Godoy-Diana R. Force balance in the take-off of a pierid butterfly: relative importance and timing of leg impulsion and aerodynamic forces. J Exp Biol. 2013; 216(18):3551-3563. doi: 10.1242/jeb.084699 PMID: 23788714

13. Maeda M. Aerodynamics of flapping flight interacting with environments [Ph.D. thesis]. Biomechanical Engineering Laboratory, Graduate School of Engineering, Chiba University; 2014.

14. Chen MW, Zhang YL, Sun M. Wing and body motion and aerodynamic and leg forces during take-off in droneflies. J R Soc Interface. 2013; 10(89). doi: 10.1098/rsif.2013.0808

15. Chen MW, Sun M. Wing/body kinematics measurement and force and moment analyses of the takeoff flight of fruitflies. Acta Mech Sinica. 2014; 30(4):495-506. doi: 10.1007/s10409-014-0059-0

16. Truong TV, Byun D, Kim MJ, Yoon KJ, Park HC. Aerodynamic forces and flow structures of the leading edge vortex on a flapping wing considering ground effect. Bioinspir Biomim. 2013; 8(3):036007. doi: 10.1088/1748-3182/8/3/036007 PMID: 23851351

17. Engels T, Kolomenskiy D, Schneider K, Sesterhenn J. FluSI: A novel parallel simulation tool for flapping insect flight using a Fourier method with volume penalization. arXiv:150606513. 2015.

18. Fontaine EI, Zabala F, Dickinson MH, Burdick JW. Wing and body motion during flight initiation in *Drosophila* revealed by automated visual tracking. J Exp Biol. 2009; 212(9):1307-1323. doi: 10.1242/jeb.025379 PMID: 19376952

19. Schneider K. Numerical simulation of the transient flow behaviour in chemical reactors using a penalis-sation method. Comput Fluids. 2005; 34:1223-1238. doi: 10.1016/j.compfluid.2004.09.006

20. Angot P, Bruneau CH, Fabrie P. A penalisation method to take into account obstacles in viscous flows. Numer Math. 1999; 81:497-520. doi: 10.1007/s002110050401

21. Engels T, Kolomenskiy D, Schneider K, Sesterhenn J. Numerical simulation of fluid-structure interac-tion with the volume penalization method. J Comput Phys. 2015; 281:96-115. doi: 10.1016/j.jcp.2014.10.005

22. Dickinson MH, Lehmann FO, Sane SP. Wing Rotation and the Aerodynamic Basis of Insect Flight. Science. 1999; 284:1954-1960. doi: 10.1126/science.284.5422.1954 PMID: 10373107

23. Zumstein N, Forman O, Nongthomba U, Sparrow JC, Elliott CJH. Distance and force production during jumping in wild-type and mutant *Drosophila melanogaster*. J Exp Biol. 2004; 207(20):3515-3522. doi: 10.1242/jeb.01181 PMID: 15339947

24. Card G, Dickinson M. Performance trade-offs in the flight initiation of *Drosophila*. J Exp Biol. 2008; 211 (3):341-353. doi: 10.1242/jeb.012682 PMID: 18203989

25. Kolomenskiy D, Maeda M, Engels T, Liu H, Schneider K, Nave JC. Aerodynamic Ground Effect in Fruit-fly Sized Insect Takeoff: Source Codes and Datasets; 2016. doi: 10.17605/OSF.ITWURBC

26. 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(13):2368-2382. doi: 10.1242/jeb.005686 PMID: 17575042

27. Wootton RJ, Herbert RC, Young PG, Evans KE. Approaches to the structural modelling of insect wings. Phil Trans Royal Soc B. 2003; 358(1437):1577-1587. doi: 10.1098/rstb.2003.1351

28. Chirarattananon P, Ma KY, Wood RJ. Adaptive control for takeoff, hovering, and landing of a robotic fly. In: IEEE International Conference on Intelligent Robots and Systems; 2013. p. 3808-3815.

29. Aono H, Liang F, Liu H. Near- and far-field aerodynamics in insect hovering flight: an integrated computational study. J Exp Biol. 2008; 211(2):239-257. doi: 10.1242/jeb.008649 PMID: 18165252
30. Ishihara D, Horie T, Denda M. A two-dimensional computational study on the fluid x2013; structure inter-action cause of wing pitch changes in dipteran flapping flight. J Exp Biol. 2009; 212(1):1-10. doi: 10.1242/jeb.020404 PMID: 19088205

31. Tanaka H, Whitney JP, Wood RJ. Effect of Flexural and Torsional Wing Flexibility on Lift Generation in Hoverfly Flight. Integr Comp Biol. 2011; 51(1):142-150. doi: 10.1093/icb/icr051 PMID: 21622947

32. Ishihara D, Horie T, Niho T. An experimental and three-dimensional computational study on the aerodynamic contribution to the passive pitching motion of flapping wings in hovering flies. Bioinspir Biomim. 2014; 9(4):046009. doi: 10.1088/1748-3182/9/4/046009 PMID: 25378268

33. Kolomenskiy D, Schneider K. A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles. J Comput Phys. 2009; 228:5687-5709. doi: 10.1016/j.jcp.2009.04.026

34. Kolomenskiy D, Moffatt HK, Farge M, Schneider K. Two- and three-dimensional numerical simulations of the clap-fling-sweep of hovering insects. J Fluids Struct. 2011; 27(5-6):784-791. doi: 10.1016/j.jfluidstructs.2011.05.002