# Performance measurement of a 50-square-meter kite set-up on a 13-meter trawler

Morgan Behrel, Kostia Roncin, Damien Grelon, Frédéric Montel, Alain Nême, Jean-Baptiste Leroux, Christian Jochum, Yves Parlier

► To cite this version:

Morgan Behrel, Kostia Roncin, Damien Grelon, Frédéric Montel, Alain Nême, et al. Performance measurement of a 50-square-meter kite set-up on a 13-meter trawler. La Houille Blanche - Revue internationale de l'eau, 2019, 2019 (5-6), pp.67-69. 10.1051/lhb/2019036. hal-02431023

HAL Id: hal-02431023

https://hal.science/hal-02431023

Submitted on 27 Sep 2021

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 Attribution 4.0 International License
# Performance measurement of a 50-square-meter kite set-up on a 13-meter trawler

Morgan Behrel<fnref n="1" />*, Kostia Roncin<fnref n="1" />, Damien Grelon<fnref n="2" />, Frédéric Montel<fnref n="1" />, Alain Nême<fnref n="1" />, Jean-Baptiste Leroux<fnref n="1" />, Christian Jochum<fnref n="1" /> and Yves Parlier<fnref n="3" />

**Abstract** – This paper describes an on board measurement campaign held in Grande-Rivière, Gaspésie, Province of Québec, Canada, in October 2015, involving a 13-meter trawler equipped with a 50-square-meter kite. The aim of the campaign was to access the boat performance when the kite is used for auxiliary propulsion. To achieve this purpose, in addition to the kite control system, a set of sensors was installed. During the trials, runs with kite in static flight were done, with around 12 knots of true wind speed. The data post processing is presented in this paper, and allows an estimate of the lift coefficient and the lift to drag ratio of the kite and the tethers. The collected data are consistent with other experimental data published.

**Keywords:** on board measurements / auxiliary propulsion / kite / sea trials

**Résumé** – Mesure des performances d'un kite de 50 mètres carrés installé sur un chalutier de 13 mètres. Cet article décrit une campagne de mesure en mer réalisée sur un chalutier de 13 mètres équipé d'un kite de 50 mètres carrés. Cette campagne s'est déroulée à Grande Rivière, Gaspésie, Province de Québec, Canada, en octobre 2015. Le but était d'évaluer les performances du navire lorsque le kite est déployé. Pour cela, en plus des dispositifs de contrôle du kite, un ensemble de capteurs a été installé. Au cours de la campagne, des phases rectilignes avec le kite fonctionnant en vol statique ont été effectuées, avec environ 12 nœuds de vent réel. Le post-traitement des données qui est présenté permet d'estimer la finesse du kite en vol statique ainsi que son coefficient de portance. Les données obtenues sont cohérentes avec d'autres résultats expérimentaux publiés.

**Mots clés** : mesures embarquées / propulsion auxiliaire / kite / essais en mer

## 1 Introduction

The current research project, undertaken by the company BEYOND THE SEA®, and managed in partnership with ENSTA Bretagne, aims to develop kites as auxiliary propulsion system for ships. In this context, a sea trial and a measurement campaign have been set up in partnership with a Canadian fishing vessel and the Merinov institute, located in Grande-Rivière, Gaspésie, Province of Québec, Canada in October 2015.

The trials were conducted on a 13-meter fishing vessel. These types of vessel (trawler) are particularly suitable for kite-assisted propulsion, because they are operated at low speed, leading to more appropriate relative wind angles.

Various structure modifications were carried out to allow the use of a towing kite, especially a reinforcement of the foredeck on which kite tethers were attached. Moreover the boat was already fitted with engine sensors (fuel flow, propeller shaft rotation speed and torque) whereby the measurement of the impact of the kite on the propulsion chain could be carried out.

The kite used during the trials had an area of 50 m<fnref n="2" /> (34 m<fnref n="2" /> of projected area), with an inflatable leading edge, and 9 inflatable battens. With this architecture, the kite has its own shape without any aerodynamic load which makes launching and recovering procedures easier. The mass of the deflated kite including bridles is 21 kg.

The measurement system was based on 3 main sensors: a 3D load cell providing orientation and intensity of the forces generated by the kite on the boat, an Inertial Measurement Unit coupled with a GPS device providing boat motion and velocity,

<footnote n="1">¹ ENSTA Bretagne, FRE CNRS 3744, IRDL, Brest, France</footnote>

<footnote n="2">² MERINOV, Grande-Rivière, Quebec, Canada</footnote>

<footnote n="3">³ BEYOND THE SEA®, La Teste-de-Buch, France</footnote>

<footnote>*Correspondence: morgan.behrel@ensta-bretagne.org</footnote>
and an anemometer. The details of the experimental set up can be found in Behrel et al. (2016).

Runs with kite in static flight were carried out. Logged data were used to compute the lift to drag ratio and the lift coefficient of the kite. The presented method is specifically adapted to the collected data. Indeed, due to a lack of information on the kite orientation during flight and on wind along the altitude, strong assumptions and estimations have to be done, leading to results which must be interpreted with caution. However, the obtained lift to drag ratio and lift coefficient make sense with other experimental works like Dadd (2012).

## 2 Post-processing of kite flight data

The aim of the following part is to process data to retrieve the lift to drag ratio of the system {kite + tethers}. However, to compute the lift to drag ratio, the wind direction at the altitude of the kite needs to be known as it directly affects the orientation of the wind window. In the same way, to compute the lift coefficient, the wind intensity at the kite altitude also needs to be known. Unfortunately, there were no reliable means to measure it. Actually, only the relative wind $V_{WR}$ over the boat, at 9 m over the sea, was measured. The true wind vector $V_{WT}$ at the same location can be easily retrieved using the velocity vector of the measurement point $V_M$ (see Eq. (1)), but the true wind at other altitudes can only be estimated. Indeed, the altitude impacts true wind both in terms of strength and direction. With no possibility of getting accurate modelling of the twist of the flow along the altitude, it has been decided to neglect the effect of twist, and only take into account the 2D planar wind distribution profile calculated using the equation (2) (ITTC, 2014). Equations (1)–(3) present the process to obtain the relative wind velocity and boat velocity measurement from relative wind velocity for any altitude z.

$$ \underline{V}_{WT}(Z_0) = \underline{V}_{WR}(Z_0) + \underline{V}_{M}, \quad (1) $$

$$ \underline{V}_{WT}(z) = \left(\frac{z}{Z_0}\right)^{1/7} \underline{V}_{WT}(Z_0), \quad (2) $$

$$ \underline{V}_{WR}(z) = \underline{V}_{WT}(z) - \underline{V}_{A}, \quad (3) $$

Where $Z_0$ is the altitude measurement and z is the altitude above sea level. $V_A$ is the velocity of kite attachment point. However, during trials, the velocities induced by ship motions at any points of the ship were observed to stay under 2% of the ship velocity magnitude. Hence, they were neglected, and $V_A$ and $V_M$ are taken equal to $V_S$, and $V_S$ is assumed to be contained in the horizontal plane, as $V_{WR}$ consequently.

With the relative wind vector defined for any altitude, we can create a new axis system based on it. The Relative Wind reference frame $R_{WR}$ is defined with the X-axis along the relative wind vector at the kite altitude, the Z-axis vertical pointing down and the Y-axis completing the coordinate system to create a direct one. To process the lift to drag ratio and the lift coefficient, the apparent wind on the kite needs to be known. For the present study, the kite is considered as to be in perfect static flight, with no velocity in the ship reference frame. Therefore, the apparent wind vector is equal to the relative wind vector at the kite altitude.

## 2.1 Kite force vector

The kite force vector $F_k$ was measured with the 3D load cell rigidly fixed to the ship, and consequently, the vector is at first expressed in the ship reference frame, denoted with superscript (s). However the following computations require a force vector expressed in the relative wind reference frame $R_{WR}$, which can be done using equation (4):

$$ \underline{F}_{k}^{(WR)} = \underline{M}(\beta_{WR}, \theta_{s}, \phi_{s}) \underline{F}_{k}^{(s)}. \quad (4) $$

With the assumptions presented previously, the transfer matrix $\underline{M}$ depends on the relative wind angle $\beta_{WR}$, and the attitude angles of the boat $\theta_s$ and $\phi_s$, respectively the pitch and roll angles. The relative wind angle comes from the construction of the relative wind presented in equation (3). The pitch and roll values are directly measured by the Inertial Measurement Unit.

## 2.2 Lift to drag ratio estimation

From there, it becomes straightforward to compute the lift to drag ratio. Indeed, the component of the force along $\underline{X}_{WR}$ is the total drag and the projection into the $(\underline{Y}_{WR}, \underline{Z}_{WR})$ plane is the lift. The component of force along $\underline{Z}_{WR}$ is the sum of the vertical aerodynamic force generated by the kite and the weight P of the kite and tethers. The latter is then subtracted to the vertical component of the measured force to get only the aerodynamic force. Finally the lift to drag ratio is achieved by processing equation (5):

$$ \frac{L}{D} = \frac{\sqrt{Fy_{WR}^2 + (Fz_{WR} - P)^2}}{F_{X_{WR}}}. \quad (5) $$

## 2.3 Lift coefficient estimation

The previous part led to the identification of lift and drag components of kite force, from the measured force expressed in the relative wind axis system. From there, and with the same assumptions, lift coefficient can be processed using equation (6).

$$ C_L = \frac{L}{\frac{1}{2} \rho_a A_k V_a^2} = \frac{\sqrt{Fy_{WR}^2 + (Fz_{WR} - P)^2}}{\frac{1}{2} \rho_a A_k V_{WR}^2}. \quad (6) $$

The projected kite area $A_k$ was used to enable comparison with Dadd (2012). The density of the air $\rho_a$ was estimated to $1.22 \text{ kg} \cdot \text{m}^{-3}$ during the measurement day (air temperature $15^\circ\text{C}$, atmospheric pressure 1012 hPa, relative humidity 70%). $V_a$ is the kite apparent wind speed, and is equal to the relative wind speed $V_{WR}$ in the specific case of static flight.

## 3 Results

Equations (5) and (6) have been processed for 3 runs with kite in static flight. Results are presented in Table 1. No kite settings were modified during or between periods. The kite position during the second period (P.2) is presented in Figure 1, expressed in relative wind reference frame.
<table><caption>**Table 1.** Compiled results of three periods with kite in static flight.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">RD</th><th rowspan="2">TWS</th><th rowspan="2">TWA</th><th colspan="2">L/D</th><th colspan="2">C<sub>L</sub></th></tr><tr><th>M</th><th>SD</th><th>M</th><th>SD</th></tr></thead><tbody><tr><td>P. 1</td><td>399</td><td>5.7</td><td>88</td><td>5.9</td><td>2.3</td><td>0.68</td><td>0.21</td></tr><tr><td>P. 2</td><td>569</td><td>6.1</td><td>95</td><td>6.1</td><td>2.5</td><td>0.76</td><td>0.26</td></tr><tr><td>P. 3</td><td>209</td><td>5.7</td><td>304</td><td>5.3</td><td>1.5</td><td>0.59</td><td>0.17</td></tr></tbody></table>

RD = Run duration (s); TWS = True wind speed (m/s); TWA = True wind angle relative to boat axis (°); M = Mean; SD = Standard deviation.

The average lift to drag ratio during the longest period is equal to 6.1, and the average lift coefficient is equal to 0.76. This seems consistent with other experimental data published, like Dadd's ones (2012). Indeed, for a 3-square-meter kite with an aspect ratio AR of 4.9, Dadd got a lift coefficient of 0.78 and a lift to drag ratio of 6.07. To estimate the drag coefficient and thus the lift to drag ratio of other kites with other aspect ratios, Dadd uses the Prandtl's formula as presented in Abbott and Von Doenhoff (1959) and given in equation (7).

$$C'_{D} = C_{D} + \frac{C_{L}^{2}}{\pi} \left( \frac{1}{AR'} - \frac{1}{AR} \right). \quad (7)$$

Assuming both kites are trimmed to produce the same lift coefficient, and applying this method to the kite used for the present study, whose aspect ratio AR is 5.5, the expected lift to drag ratio should be 6.28.

## 4 Discussion

The averaged results of the 3 periods are close, with a good agreement with published data, as it has been shown previously. However, a point-to-point data analysis shows extreme values of lift to drag ratios and lift coefficients that are not realistic. This demonstrates the limits of the various assumptions which have been made, particularly the perfect static flight assumption. Indeed, Figure 1 shows that the kite position during a run is a cloud of points, whereas it should be, theoretically, a single point. The straight-line assumption of the tether could also be a source of uncertainty. The decision to disregard the twist wind flow along the altitude due to a lack of data and models could also be detrimental. Indeed, a fourth period with kite in static flight has been logged, but it had to be discarded due to inconsistent data. One possible explanation is a significant shift between the wind orientation at the measurement point and the wind orientation at the kite position. This eventuality leads us to consider, for future

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Kite position and generated force during Period 2, the red triangle denotes the boat position.</figcaption></figure>

experimental campaign, a duplication of wind measurement systems at different altitude. Moreover, the installation of multiple anemometers at various altitudes could also be a good way to improve wind estimation with altitude. This approach has already been undertaken during another experimental campaign, taking place on shore, with the kite control system firmly fixed to the ground. A wind profiler (SODAR type) was used to get the wind intensity and direction along altitude. The results will be presented and discussed in a future paper.

A brief analysis of power can be carried out: the average value of the propulsive force generated by the kite during period 2 is 505 N, for a speed of the boat equals to 2.4 m/s. The power of the kite is then 1212 W. On the engine side the propeller shaft measurements (rotation speed and torque) led to a transmitted power equals to 38 000 W. Assuming a propulsive efficiency of 50%, the kite therefore provided more than 6% of the propulsive power.

*Acknowledgments.* The authors wish to express their deep appreciation for the full support of the Steven Paul's master-owner, Guy Leblanc, and his great motivation to innovate. The authors are also grateful to the French agency for energy development and control (ADEME) for the funding of this study.

## References

Abbott IH, Von Doenhoff AE. 1959. Theory of wing sections. Mineola: Dover Publications Inc.

Behrel M, Roncin K, Grelon D, et al. 2016. Presentation of an experimental campaign for measuring performances and effects of a 50-square-meter kite put on a 13-meter trawler. In: *15e Journées de l'Hydrodynamique, Brest*.

Dadd GM. 2012. Kite dynamics for ship propulsion. Southampton: University of Southampton.

ITTC. 2014. Symbols and terminology list. In: *International Towing Tank Conference, Copenhagen*.