This paper presents the experimental validation and dynamic similarity analysis for a lab-scale version of an airborne wind energy (AWE) system executing closed-loop motion control. Execution of crosswind flight patterns, achieved in this work through the asymmetric motion of three tethers, enables dramatic increases in energy generation compared with stationary operation. Achievement of crosswind flight in the lab-scale experimental framework described herein allows for rapid, inexpensive, and dynamically scalable characterization of new control algorithms without recourse to expensive full-scale prototyping. We first present the experimental setup, then derive dynamic scaling relationships necessary for the lab-scale behavior to match the full-scale behavior. We then validate dynamic equivalence of crosswind flight over a range of different scale models of the Altaeros Buoyant airborne turbine (BAT). This work is the first example of successful lab-scale control and measurement of crosswind motion for an AWE system across a range of flow speeds and system scales. The results demonstrate that crosswind flight can achieve significantly more power production than stationary operation, while also validating dynamic scaling laws under closed-loop control.

## Introduction

Wind energy is one of several renewable energy sources that continue to play an increasing role in satisfying the world's energy demand. However, tower costs and installation logistics typically limit towered wind turbine installations to nonremote and shallow-water offshore locations. Furthermore, towers limit hub heights to approximately 100 m, which does not provide access to stronger high-altitude winds. Airborne wind energy (AWE) systems, which replace towers with tethers and a lifting body (a kite, wing, or aerostat), enable operation at altitudes up to 600 m, where winds are stronger and more consistent, all while eliminating up to 90% of the material required for towered systems (see Ref. [1]). Simultaneously, these systems produce power at a lower price to the consumer than diesel generators (see Refs. [2] and [3]). This makes AWE systems a good fit for remote and off-grid communities. Inspired by these possibilities, numerous companies and research organizations (see Refs. [1] and [4–9]) have pioneered unique designs for AWE systems over the past decade. Several AWE systems are depicted in Fig. 1. These systems either feature airborne power generation (as is the case with Refs. [1] and [5], where power is generated by an on-board turbine(s) and transmitted to the ground) or ground-based power generation (as is the case with Ref. [4], where power is generated on the ground as a result of tension in the tethers as they are spooled out).

In addition to enabling high-altitude operation, the replacement of a tower with tethers enables motion through the air. While this may initially be viewed as a drawback, appropriately designed periodic motions have been shown in Ref. [10] to significantly enhance the power output of an AWE system. Specifically, by flying an AWE system perpendicular to the wind direction, the apparent wind speed presented to the airborne system is increased. This motion, termed *crosswind* flight, translates into *significantly* increased power production, as the power production from a wind turbine is proportional to the apparent wind speed *cubed*. The promise of significantly increased energy generation has led to a growing body of literature that focuses on crosswind flight control algorithms, including [11–16].

The vast majority of existing literature on AWE system flight control has been simulation based. It is well established that while existing models capture the key variables and phenomena involved with AWE flight dynamics, they fall short of predicting significant dynamic characteristics [17,18]. The limitations of existing models arise from necessary model order reduction, use of lumped, steady aerodynamic models, and unknown model parameters.

Given the nascent state of AWE control systems, inaccuracies in existing simulation models, and significant full-scale experimental prototyping costs, an inexpensive, small-scale experimental framework for AWE system flight characterization is highly desirable. However, achieving such a framework in a way that is approximately dynamically scalable to full-scale systems is a challenging engineering feat. This has limited the number of efforts in the AWE community to develop a lab-scale framework.

Within the limited body of literature that addresses small-scale experimental characterization of crosswind systems, the authors of Ref. [19] built a *small*-scale (approximately 30 m maximum line length) system for characterizing and optimizing figure-eight crosswind motions. Meanwhile, the authors of Ref. [20] have developed a wind tunnel-based *lab*-scale experimental framework for a rotating cylinder that uses the Magnus effect to generate lift. However, the use of a heavier-than-air model in a wind tunnel at this scale does not allow for dynamic similarity between lab-scale and full-scale flight, due to the tremendous resulting differences in characteristic oscillation frequencies.

The authors of Refs. [17], [18], and [21–24], on the other hand, have developed a *water channel*-based framework for lab-scale experimental characterization of AWE systems, with the goal of establishing dynamic similarity results between this lab-scale setup and its full-scale counterpart. Work in Refs. [17] and [18] proposed an initial framework necessary for measuring position and orientation of a lab-scale system during *passive* flight. This work was then extended in Refs. [21–23] to include active control to *constant* set points, resulting in nominally stationary flight. In order to implement crosswind flight at lab scale, a significantly enhanced motion capture system was introduced in Ref. [24] for the purpose of measuring model position and orientation over a much larger range of dynamic motions. This led to the first demonstration of lab-scale crosswind flight in Ref. [24], with an evaluation of crosswind flight over multiple flow conditions (and a single model) presented in Ref. [16].

While our previous conference publications [16,24] demonstrate successful crosswind flight under a limited number of flow conditions, these references focus on a single model, do not include a controller scaling analysis that demonstrates equivalence between lab-scale and full-scale, and do not account for the energy *consumed* in crosswind flight. In fact, all of the previous scaling analysis for the water channel system, which is summarized in Ref. [23], is restricted to the passive system, without consideration of control parameters. The present paper is the first to establish dynamic scaling results for control parameters, for both stationary and crosswind closed-loop flight. Thus, it serves as a greatly extended version of our previous conference results.

The contributions of this work are as follows:

The first lab-scale demonstration of closed-loop crosswind flight at lab scale, for a variety of flow conditions and model scales.

An enhanced motion capture system capable of operating under rapidly varying trajectories and lighting conditions.

Dynamic scaling relationships that indicate how controller parameters need to be scaled to achieve similarity between lab-scale and full-scale flight.

Experimental validation and error analysis of dynamic scaling laws under closed-loop crosswind flight control.

Energy generation performance analysis.

After describing the full-scale system, dynamic model, and crosswind flight control strategy in Sec. 2, we detail the lab-scale system and present the augmented motion capture system in Sec. 3. Section 4 details the dynamic similarity between the lab-scale and full-scale systems, presenting scaling laws necessary to achieve dynamic similarity under closed-loop control. Results in Sec. 5 show that:

Crosswind flight is effectively induced and controlled via this strategy;

this crosswind flight can produce significant increases in energy generation performance; and

dynamic behavior across three model scales obeys scaling laws derived via Buckingham-Pi analysis to within experimental tolerances.

We conclude with a discussion of possible sources of error inherent in the water channel setup.

## Full-Scale System: Modeling and Flight Control

The various designs of AWE systems can be controlled though either active aerodynamic surfaces (as is the case with the Google-Based Makani Power [1] and Ampyx [6] systems, among others) or tether articulation (used by Altaeros Energies [5] and KITEnrg [4], among others). The control strategy utilized in this work induces crosswind motion through the asymmetric articulation of tethers. For this purpose, we focus specifically on the Altaeros Buoyant air turbine (BAT). The BAT, shown in Fig. 2, is a lighter-than-air design wherein a horizontal-axis turbine is suspended within an annular shell. The system is intended for flight at altitudes up to 600 m, where it can harness energy from strong and consistent winds. The BAT consists of four main components, namely, a helium-filled annular shroud that holds the turbine, the turbine itself, three tethers (one of which is a conductive cable), and a ground station. One tether is attached at the center of the rear (aft) of the lifting body. The other two are attached on either side (port and starboard) of the front of the lifting body. While the BAT is nominally designed for stationary flight, the cross section of the annular shell is composed of lifting airfoils, thereby producing the aerodynamic forces that make it possible to realize increased power generation under crosswind motion.

### Dynamic Model.

The dynamic scaling analysis described in this paper is based on a relatively simple model, discussed in Ref. [17], that captures the key variables involved in the AWE system but does not possess the fidelity of real flight experiments. Hence, the model contains the necessary information for conducting a dynamic scaling analysis; nevertheless, it is by no means a replacement for experiments.

*L*and a spherical bridle joint as shown in Fig. 3. This approximation results in a model that is fully described by ordinary differential equations, without algebraic constraints. The model derived therein has the form

_{t}*v*

_{wind}is the wind speed,

*ψ*

_{wind}is the wind heading, and

*L*, are shown in Fig. 3. The vector

_{t}*X*includes the ground frame positions, denoted by the subscript “

*g*,” and the body frame translational velocities, denoted by the subscript “

*b*.” The vector Ω includes the rotational degrees-of-freedom (Euler angles), along with the body frame rotational velocities. The quantities $\theta \u2032$ and $\varphi \u2032$ are termed the

*induced*pitch and the induced roll, respectively. Along with yaw, they describe the orientation of the lifting body relative to its position vector. These angles are approximately related to the tether lengths by

For the full-scale system, the tether speed commands are realized through a variable frequency drive system. For the lab-scale system, the tether release speed commands are realized through direct current (DC) motors whose input voltages are proportional to the tether release speeds.

### Basic Flight Control.

The goal of the BAT's control system is to use three tethers to impart desirable flight characteristics to the system. In the case of stationary flight, this involves using the three tethers to keep the altitude, pitch angle, and roll angle fixed at desirable values. For crosswind flight, detailed in Sec. 2.3, the side tethers will be used to periodically alter the roll angle setpoint for the purpose of power augmentation.

*θ*, is controlled by spooling out on the aft tether while spooling in on the front tethers, or vice versa. Roll, $\varphi $, is controlled by spooling out on the port tether and spooling in on the starboard tether, or vice versa. Altitude,

*z*, is controlled by spooling all tethers in or out simultaneously. Each of three controllers, namely, the altitude controller, pitch angle controller, and roll angle controller, is composed of a lead filter (i.e., a filtered proportional plus derivative controller). A block diagram of the control strategy is shown in Fig. 4. The individual controller transfer functions are

_{g}*v*

_{1},

*v*

_{2}, and

*v*

_{3}through a linear combination (see Fig. (4)) where

Because there exists a pure integrator between the motor voltage commands and the tether lengths (*l*_{1}, *l*_{2}, and *l*_{3}, which dictate *r _{z}*,

*θ*, and $\varphi $), an integrator is not included in the controller itself, as it is unnecessary for tracking sufficiently slowly varying commands and would lead to reduced stability margins if included.

### Crosswind Augmentation.

The Altaeros BAT, which is used as a case study in this work, is designed for nominally stationary operation; thus, it does not possess the high lift/drag ratios and nimble dynamics of some other lifting bodies, which are required to realize several times as much power under crosswind flight than in stationary operation. Nevertheless, its annular design does allow for significant lift and side force, thereby allowing for moderately increased power production under crosswind motions.

*v*, is a combination of the wind velocity vector $vf$ and the velocity vector of the turbine $v$. It is computed by a simple vector difference. This is shown in Fig. 5 (left) and given by

_{a}*power ratio*,

*P*. The power ratio characterizes the amount of additional power that can be created by the crosswind motion as compared with the power that could be generated by the same turbine under stationary flight, assuming that the turbine is operating below its rated wind speed at all times. This ratio is given by

*ρ*represents the fluid density,

*C*represents the turbine's power coefficient,

_{p}*A*

_{ref}represents the rotor area, and

*τ*represents individual tether tensions. The term $kw(ui(t))$ represents a lumped winch efficiency term (accounting for regenerative capabilities) given by

_{i}where *η _{m}* and

*η*represent the motor efficiency and regenerative efficiency, respectively.

_{g}*z*and

_{sp}*θ*constant, while $\varphi sp$ is varied periodically, according to a square wave profile given by

_{sp}*a*

_{0}, and period,

*T*, represent tunable parameters that have a significant impact on the energy that can be produced through crosswind flight trajectories.

## Lab-Scale Experimental Setup

This work is based on a small-scale setup for characterizing the dynamics and control of an AWE system in a 1 m × 1 m *water channel*, as shown in Fig. 6. Three-dimensional (3D) printed models, whose cross-sectional areas comprise less than 1.5% of the total water channel cross-sectional area, are tethered and “flown” in the water channel. For the full-scale system, resolving position and orientation simply involves a global positioning system and inertial measurement unit. Since that equipment is impractical for a small-scale model in water, it was necessary to replace the motion capture system with a wholly different system. The method of motion capture presented here produces the same outputs as the full-scale system, namely, the position and orientation of the BAT. This makes it possible to accurately model the structure of the controller used at full scale. Figure 7 shows a detailed comparison between the control strategies of the full-scale and lab-scale systems.

During an experiment, an individual model is tethered and “flown” in the water channel experimental setup shown in Fig. 6. The experimental equipment consists of three DC motors for tether actuation, three high-speed cameras for image acquisition, and a high-performance host/target computer pair for real-time motion capture and closed-loop control. The motor time constants for the micro-DC motors used in this work are more than an order of magnitude faster than the desired closed-loop system time constants. The image processing and real-time controller are executed on a dedicated personal computer, using Mathworks (Natick, MA) Simulink Real-Time^{®} software. The system uses Basler ACE 340 km grayscale cameras to obtain the images of the model.

The use of three cameras at different orientations (two perpendicular to the flow and one at a 45 deg “slant” underneath the water channel) provides sufficient imaging data for resolving both the position and orientation of the model. In this work, the model position vector, $r=[xg\u2009yg\u2009zg]$, points from the origin of the ground-fixed coordinate system to the origin of the body-fixed coordinate system at the models center of mass. The model orientation (attitude) is characterized through roll ($\varphi $), pitch (*θ*), and yaw (*ψ*) Euler angles.

In order to design a controller capable of producing crosswind flight, it is necessary to resolve the position and orientation of the model in real time. A schematic representing the system is shown in Fig. 8. It shows the ground-fixed coordinate system, denoted by the subscript “g,” as well as the body-fixed coordinate system, denoted by the subscript “b.” All of the key variables used in this work are summarized in Table 1.

Variable | Description | Units |
---|---|---|

$r$ | Position of model | cm |

$v$ | Velocity of model | $cm\xb7s\u22121$ |

$vf$ | Free stream flow velocity | $cm\xb7s\u22121$ |

$va$ | Velocity of apparent wind | $cm\xb7s\u22121$ |

$\Phi $ | Azimuth angle | — |

Λ | Zenith angle | — |

$\varphi ,\theta ,\psi $ | Roll, pitch, and yaw angles | — |

$x\u0302b,y\u0302b,z\u0302b$ | Unit vectors of body-fixed coordinate system | cm |

$l1,l2,l3$ | Center, starboard, and port tether lengths | cm |

$r1,r2,r3$ | Camera #1, #2, and #3 positions in ground-fixed coordinates | cm |

$u1,u2,u3$ | Commanded tether release speeds | $cm\xb7s\u22121$ |

Variable | Description | Units |
---|---|---|

$r$ | Position of model | cm |

$v$ | Velocity of model | $cm\xb7s\u22121$ |

$vf$ | Free stream flow velocity | $cm\xb7s\u22121$ |

$va$ | Velocity of apparent wind | $cm\xb7s\u22121$ |

$\Phi $ | Azimuth angle | — |

Λ | Zenith angle | — |

$\varphi ,\theta ,\psi $ | Roll, pitch, and yaw angles | — |

$x\u0302b,y\u0302b,z\u0302b$ | Unit vectors of body-fixed coordinate system | cm |

$l1,l2,l3$ | Center, starboard, and port tether lengths | cm |

$r1,r2,r3$ | Camera #1, #2, and #3 positions in ground-fixed coordinates | cm |

$u1,u2,u3$ | Commanded tether release speeds | $cm\xb7s\u22121$ |

The unit vectors representing the body fixed *x*, *y*, and *z* axes are represented as three-component vectors in the ground-fixed system. Figure 8 also shows the lengths of the center tether, *l*_{1}, starboard tether, *l*_{2}, and port tether, *l*_{3}.

The azimuth angle characterizes the extent by which the lifting body deviates from a direct downwind configuration. Under successful crosswind flight, $\Phi $ will exhibit periodic oscillations about Φ = 0.

### Motion Capture System.

Successful crosswind flight in the lab-scale setup requires continuous, accurate measurements of the lifting body's position and orientation (attitude). This means that the three cameras depicted in Fig. 8 must continually track the white dots on the model, as shown in Fig. 9. A basic threshold-based image processing algorithm was implemented in the previous work [22], for accomplishing this on a stationary AWE system; however, successful motion capture for crosswind flight places increasing demands on the underlying image processing algorithms, as they must track a fast-moving system under varying lighting conditions. The enhanced system has augmented the effective field of view and increased the degree of dynamic behavior that can be accurately tracked. This subsection details the enhanced motion capture system that has been implemented for these purposes, which consists of two key steps:

- (1)
Dot centroid detection in a single still frame.

- (2)
Derivation of AWE position and orientation from dot centroid locations.

The location of this improved motion capture algorithm within the overall control scheme is shown in the feedback loop of the lab-scale system in Fig. 7.

#### Dot Centroid Detection.

**G**is a two-dimensional matrix of gray-scale values ranging from 0 to 255. To classify each pixel in the ROI as belonging to a dot or the background, a threshold-based technique known as

*Otsu's method*[25] is used. This method works by generating an intensity histogram from the ROI matrix and then calculating the intra-class variance, $\sigma w2$ for each intensity value

*ξ*. The intra-class variance is defined in the following equation:

*ω*is the weighted probability of the class to the left of the point being considered on the histogram and

_{L}*ω*is the weighted probability of the class on the right of the point being considered on the histogram. Variances appear in Eq. (17) as $\sigma L2$ and $\sigma R2$. These represent the variance of the left and right classes, respectively. It is shown in Ref. [25] that the value

_{R}*ξ*

_{0}that minimizes $\sigma w2$

*ξ*

_{0}, all pixels in the ROI are classified as 1 or 0 in a binary matrix,

**B**, whose entries are specified as follows:

**B**(horizontal if the matrix is divided horizontally; vertical if the matrix is divided vertically). Once $B1$ and $B2$ have been obtained, coordinates of the dot centers, denoted by $xc1,yc1,xc2,$ and $yc2$, are located by calculating individual components of a discrete two-dimensional geometric centroid according to the following equation:

Here, $b1ij$ is the component located in the *i*th row and the *j*th column of the first binary matrix, $B1$, and $b2ij$ is the component located in the *i*th row and the *j*th column of the second binary matrix, $B2$. The individual position components are then transformed back into coordinates in the original image by simple addition with the position coordinates of the ROI within the overall image.

where the index $i\u2208{1,2,3,4}$ specifies the combination of camera and dots being tracked.

#### Derivation of Orientation and Position.

Each set of dots is aligned with either the body-fixed *x* or *y* axis; specifically, two of the ROIs track dots oriented along the *x _{b}* axis, whereas the other two track dots that are oriented along the

*y*axis. Table 2 delineates the precise geometric information that is provided by each camera and each set of dots. This information is sufficient for the calculation of the unit vectors $x\u0302b$ and $y\u0302b$ in ground-fixed coordinates. Once $x\u0302b$ and $y\u0302b$ have been determined, $z\u0302b$ is determined from a simple cross product; specifically, $z\u0302b=x\u0302b\xd7y\u0302b$. The complete set of body-fixed unit vectors in the ground-fixed coordinate system forms the columns of a rotation matrix that encodes the Euler angles, namely, roll ($\varphi $), pitch (

_{b}*θ*), and yaw (

*ψ*).

Unit vector | Component | Camera | Dot set |
---|---|---|---|

$x\u0302b$ | $xg$ | 1 | 1 |

$x\u0302b$ | $yg$ | 2 | 3 |

$x\u0302b$ | $zg$ | 1 | 1 |

$y\u0302b$ | $xg$ | 2 | 2 |

$y\u0302b$ | $yg$ | 2 | 2 |

$y\u0302b$ | $zg$ | 3 | 2 |

Unit vector | Component | Camera | Dot set |
---|---|---|---|

$x\u0302b$ | $xg$ | 1 | 1 |

$x\u0302b$ | $yg$ | 2 | 3 |

$x\u0302b$ | $zg$ | 1 | 1 |

$y\u0302b$ | $xg$ | 2 | 2 |

$y\u0302b$ | $yg$ | 2 | 2 |

$y\u0302b$ | $zg$ | 3 | 2 |

where $r2,x$ is the *x* component of the bottom, upstream camera position vector, $r2,y$ is the *y* component of the bottom upstream camera position vector, and $r1,z$ is the *z* component of the side camera position vector. Furthermore, *N _{H}* and

*N*are the number of horizontal and vertical pixels in an image, and

_{V}*γ*is the view angle of the camera relative to the centerline. The position of the dots within each frame is denoted as $vdots,1$ for the vertical position of the dots in the frame from the bottom upstream camera, $hdots,2$ for the horizontal position of the dots in the frame from the bottom upstream camera, and $vdots,1$ for the vertical position of the dots in the frame from the front camera. Note that the vertical position in the frame is measured from the top of the image down.

Because this system relies on digital images, it produces quantization error. In general, the magnitude of this error will be a nonlinear function of the true position and orientation. However, it is possible to provide an estimate of the typical (average) value of this error, which will not vary by more than approximately 10% over the course of flight. Typically, the total distance between dot centers in an image is image approximately 80 pixels. Therefore, the quantization error in orientation, calculated through a simple arc tangent, is approximately 0.7 deg. The quantization error in position can be approximated by using the derivative of Eq. (22) with respect to $r$. For a set of dots in the middle of the camera field of view, and the center of the water channel, it is found to be approximately 1 mm.

## Dynamic Scaling Analysis for Lab-Scale/Full-Scale Equivalence

In order to compare the closed-loop dynamic behavior across model scales, we utilize the Buckingham-Pi theorem of Ref. [26] to derive relationships that describe how to appropriately scale plant and controller parameters. Although the dynamic model detailed in Eqs. (1)–(6) does not provide the level of detail that actual flight tests provide (hence the desire for an experimental setup), it does provide enough information to identify key variables in the system dynamics which is essential for applying the Buckingham-Pi theorem. This section first reviews the methods and results of Ref. [23], which are used to scale the *plant*. Then, we derive a new set of scaling laws that describe how the controller of the lab-scale system should be scaled to achieve dynamic equivalence. Together, these two sets of scaling laws provide a set of guidelines for achieving dynamic similarity between the lab-scale and full-scale systems under closed-loop control.

### Dynamic Scaling Analysis of Plant Parameters.

*t*), through a scalar parameter,

*α*

Parameter | Variable | Units |
---|---|---|

Mean tether length | $Lt$ | $m$ |

Aerodynamic chord length | $Lc$ | $m$ |

Tether diameter | $dt$ | $m$ |

Mass of aerostat | $Ma$ | $kg$ |

Fluid density | $\rho $ | $kg\xb7m\u22123$ |

Wind or flow speed | $vf$ | $m\xb7s\u22121$ |

Gravitational acceleration | $g$ | $m\xb7s\u22122$ |

Azimuth rate of change | $\Phi \u02d9$ | $s\u22121$ |

Azimuth acceleration | $\Phi \xa8$ | $s\u22122$ |

Zenith rate of change | $\Lambda \u02d9$ | $s\u22121$ |

Zenith angle acceleration | $\Lambda \xa8$ | $s\u22122$ |

Tether twist angle rate of change | $\Psi \u02d9$ | $s\u22121$ |

Tether twist acceleration | $\Psi \xa8$ | $s\u22122$ |

Linear tether density | $\lambda $ | $kg\xb7m\u22121$ |

Parameter | Variable | Units |
---|---|---|

Mean tether length | $Lt$ | $m$ |

Aerodynamic chord length | $Lc$ | $m$ |

Tether diameter | $dt$ | $m$ |

Mass of aerostat | $Ma$ | $kg$ |

Fluid density | $\rho $ | $kg\xb7m\u22123$ |

Wind or flow speed | $vf$ | $m\xb7s\u22121$ |

Gravitational acceleration | $g$ | $m\xb7s\u22122$ |

Azimuth rate of change | $\Phi \u02d9$ | $s\u22121$ |

Azimuth acceleration | $\Phi \xa8$ | $s\u22122$ |

Zenith rate of change | $\Lambda \u02d9$ | $s\u22121$ |

Zenith angle acceleration | $\Lambda \xa8$ | $s\u22122$ |

Tether twist angle rate of change | $\Psi \u02d9$ | $s\u22121$ |

Tether twist acceleration | $\Psi \xa8$ | $s\u22122$ |

Linear tether density | $\lambda $ | $kg\xb7m\u22121$ |

*α*), at small scale versus full-scale. The new dimensionless groups, where derivatives are calculated relative to the normalized time, $t\xaf$, are

These dimensionless groups lead to the following guidelines for achieving dynamic similarity:

- (1)
Individual tether lengths

*l*_{1},*l*_{2}, and*l*_{3}must be scaled by the same factor as the model length, thereby maintaining geometric similarity. - (2)
Net buoyancy $(Fbouyant/mg)$ must be preserved at lab-scale and full scale.

- (3)
Taking

*L*as the length-scale factor (i.e., a 1:100-scale model would have $L=(1/100)$), flow speed should be scaled according to the*square root*of*L*, i.e., $vf\u221dLvfull$, where*v*_{full}is the target velocity of the full-scale system for which we wish to replicate flight characteristics. This scaling rule ensures preservation of the fifth dimensionless group (the inverse of the Froude number squared) at scale.

These scaling rules are valid under two assumptions:

- (1)
Mass distribution is identical at lab-scale and full-scale, which means that the moments of inertia are completely determined by the mass and scale factor,

*L*. This assumption is satisfied through proper design of 3D printed models. - (2)
Aerodynamic coefficients are independent of the Reynolds number within the Reynolds number ranges seen in the water channel (75,000 at a 1 m/s flow speed, using chord length as the reference length) and full-scale system (5,000,000 at a 10 m/s wind speed). Hence, the viscosity of the medium is neglected as a physical dependency. The approximate validity of this assumption has been verified experimentally in Ref. [23] for the Altaeros BAT.

### Dynamic Scaling Analysis of Controller Parameters.

In order to compare dynamic behavior under closed-loop control across model scales, it is necessary to extend the dimensional analysis of Ref. [23] to include controller parameters. From the controller detailed in Eq. (10) and Fig. 4, we identify the gains, time constants, and setpoints of the controllers as the relevant parameters. A detailed explanation of each parameter, along with associated units, is provided in Table 4.

Parameter | Variable | Units |
---|---|---|

Roll controller proportional gain | $kp\varphi $ | $m\xb7s\u22121$ |

Pitch controller proportional gain | $kp\theta $ | $m\xb7s\u22121$ |

Altitude controller proportional gain | $kpz$ | $s\u22121$ |

Roll controller derivative gain | $kd\varphi $ | m |

Pitch controller derivative gain | $kd\theta $ | m |

Altitude controller derivative gain | $kdz$ | — |

Roll controller time constant | $\tau \varphi $ | s |

Pitch controller filter time constant | $\tau \theta $ | s |

Altitude controller filter time constant | $\tau z$ | s |

Pitch setpoint | $\theta SP$ | — |

Altitude setpoint | $ZSP$ | m |

Roll setpoint amplitude | $a0$ | — |

Roll setpoint period | $T$ | s |

Tether release/retract speed limit | $vlsat$ | $m\xb7s\u22121$ |

Parameter | Variable | Units |
---|---|---|

Roll controller proportional gain | $kp\varphi $ | $m\xb7s\u22121$ |

Pitch controller proportional gain | $kp\theta $ | $m\xb7s\u22121$ |

Altitude controller proportional gain | $kpz$ | $s\u22121$ |

Roll controller derivative gain | $kd\varphi $ | m |

Pitch controller derivative gain | $kd\theta $ | m |

Altitude controller derivative gain | $kdz$ | — |

Roll controller time constant | $\tau \varphi $ | s |

Pitch controller filter time constant | $\tau \theta $ | s |

Altitude controller filter time constant | $\tau z$ | s |

Pitch setpoint | $\theta SP$ | — |

Altitude setpoint | $ZSP$ | m |

Roll setpoint amplitude | $a0$ | — |

Roll setpoint period | $T$ | s |

Tether release/retract speed limit | $vlsat$ | $m\xb7s\u22121$ |

*k*,

_{dz}*θ*, and

_{SP}*a*

_{0}are already dimensionless, they can be removed from the analysis, and the resulting analysis produces 11 dimensionless groups

Given the analysis results of Eqs. (24) and (25) in Sec. 4.1, which demonstrate that time (and all time constants) scale according to $t\u221dLc$, along with the results of Eq. (26), we arrive at the following scaling laws:

- (1)
Scale the derivative gains of the roll and pitch controllers, $kd\theta $ and $kd\varphi $, according to the length scale,

*L*._{c} - (2)
Scale the proportional gains of the roll and pitch controllers, $kp\theta $ and $kp\varphi $, along with the tether speed saturation limit, $vlsat$ according to $Lc$.

- (3)
Scale the proportional gain of the altitude controller,

*k*, according to $(1/LC)$._{pz} - (4)
Scale the filter time constants, $\tau z,\theta ,\varphi $ as well as the roll setpoint period,

*T*, according to $Lc$.

Under the aforementioned scaling laws for both the plant and the controller, the full-scale dynamic behavior will be replicated at lab scale with the only exception being uniformly accelerated time constants. In particular, the lab-scale time constants will scale uniformly according to the square root of the length scale, i.e., $\tau \u221dLc$, where *τ* represents any system time constant (Table 5).

Model size | Flow speed, $vf\u2009(cm/s)$ | Normalized flow speed, $v\xaff$$(cm/s)$ | Flow speed index | Roll setpoint period range $(s)$ | Normalized roll setpoint period range $(s)$ |
---|---|---|---|---|---|

Small | 22.0 | 25.4 | A | 4.330–7.578 | 5.000–8.750 |

Medium | 25.0 | 25.0 | A | 6.000–8.750 | 6.000–8.750 |

Large | 28.0 | 25.0 | A | 6.708–10.90 | 6.000–9.750 |

Small | 16.5 | 19.1 | B | 4.330–6.495 | 5.000–7.500 |

Medium | 18.8 | 18.8 | B | 4.500–8.500 | 4.500–8.500 |

Large | 21.0 | 18.8 | B | 6.708–9.783 | 6.000–8.750 |

Model size | Flow speed, $vf\u2009(cm/s)$ | Normalized flow speed, $v\xaff$$(cm/s)$ | Flow speed index | Roll setpoint period range $(s)$ | Normalized roll setpoint period range $(s)$ |
---|---|---|---|---|---|

Small | 22.0 | 25.4 | A | 4.330–7.578 | 5.000–8.750 |

Medium | 25.0 | 25.0 | A | 6.000–8.750 | 6.000–8.750 |

Large | 28.0 | 25.0 | A | 6.708–10.90 | 6.000–9.750 |

Small | 16.5 | 19.1 | B | 4.330–6.495 | 5.000–7.500 |

Medium | 18.8 | 18.8 | B | 4.500–8.500 | 4.500–8.500 |

Large | 21.0 | 18.8 | B | 6.708–9.783 | 6.000–8.750 |

## Experimental Validation Results and Discussion

To experimentally evaluate crosswind flight and dynamic similarity at different scales, we consider three different scale models, as shown in Fig. 10, rather than comparing a single lab-scale model to a full-scale model (which we have very limited access to and no means of controlling the experiment, since wind is random). These models are termed the small, medium, and large models. The small-scale model has a chord length, $Lcs$, equal to 75% of the chord length of the medium-scale model, $Lcm$, and the large-scale model has a chord length, $Lcl$, equal to 125% of the medium-scale model. The medium model is a 1:100 scale model of the 2013 Altaeros BAT prototype. The largest of these models presents a blockage ratio of approximately 1.5% in the UNC Charlotte water channel.

Each of the three scale models was subjected to the crosswind flight control algorithm of Sec. 2.3 at two different flow speeds, for a duration of at least 20 periods corresponding to between 90 s and 210 s depending on the roll setpoint period being tested. This was found to be sufficient for the mean power ratio to converge. At both normalized flow speeds, each model was subjected to between 11 and 18 experiments wherein some experiments were run multiple times to ensure consistency. Periods were chosen to provide a normalized period resolution of 0.25 s. The organization of individual experiments is shown in Table 6. Note that the exact matching of normalized flow speeds is limited by quantization of the flow speed measurement, which is discussed in detail in Sec. 5.4.

Model size | Normalized | Roll angle variance (rad^{2}) | Yaw angle variance (rad^{2}) | Azimuth angle variance (rad^{2}) |
---|---|---|---|---|

Small | No | 0.0209 | 0.0309 | 0.0756 |

Small | Yes | 0.00949 | 0.00463 | 0.0141 |

Medium | No | 0.0221 | 0.0365 | 0.0971 |

Medium | Yes | 0.00533 | 0.00809 | 0.0207 |

Large | No | 0.0228 | 0.0287 | 0.108 |

Large | Yes | 0.00913 | 0.00903 | 0.0462 |

Model size | Normalized | Roll angle variance (rad^{2}) | Yaw angle variance (rad^{2}) | Azimuth angle variance (rad^{2}) |
---|---|---|---|---|

Small | No | 0.0209 | 0.0309 | 0.0756 |

Small | Yes | 0.00949 | 0.00463 | 0.0141 |

Medium | No | 0.0221 | 0.0365 | 0.0971 |

Medium | Yes | 0.00533 | 0.00809 | 0.0207 |

Large | No | 0.0228 | 0.0287 | 0.108 |

Large | Yes | 0.00913 | 0.00903 | 0.0462 |

The results presented in this section have three purposes:

- (1)
To demonstrate successful crosswind flight at many different operating conditions and model scales.

- (2)
To validate previously derived dynamic scaling laws for the controller parameters.

- (3)
To demonstrate an increase in power production potential from crosswind flight, quantified by the power ratio of Eq. (13).

### Dynamic Scaling of Controller for Experiments.

During the course of a single experiment, the controller for each model was held constant. It was then varied across model scales according to the scaling laws derived in Sec. 4.2, using the scale factors of Eq. (28). The exact values used for controller parameters are given in Table 7.

Parameter | Value–small model | Value–medium model | Value–large model |
---|---|---|---|

$kp\varphi $ | 8.66 | 10.0 | 11.18 |

$kp\theta $ | 4.33 | 5.00 | 5.59 |

$kpz$ | 0.115 | 0.10 | 0.089 |

$kd\varphi $ | 2.25 | 3.00 | 3.75 |

$kd\theta $ | 0.75 | 1.00 | 1.25 |

$kdz$ | 0.00 | 0.00 | 0.00 |

$\tau \varphi $ | 0.13 | 0.15 | 0.168 |

$\tau \theta $ | 0.173 | 0.20 | 0.224 |

$\tau z$ | 0.217 | 0.250 | 0.280 |

$\theta sp$ | 5.00 | 5.00 | 5.00 |

$zsp$ | 45.0 | 60.0 | 75.0 |

$a0$ | 15.0 | 15.0 | 15.0 |

$vlsat$ | 0.0677 | 0.0606 | 0.0524 |

Parameter | Value–small model | Value–medium model | Value–large model |
---|---|---|---|

$kp\varphi $ | 8.66 | 10.0 | 11.18 |

$kp\theta $ | 4.33 | 5.00 | 5.59 |

$kpz$ | 0.115 | 0.10 | 0.089 |

$kd\varphi $ | 2.25 | 3.00 | 3.75 |

$kd\theta $ | 0.75 | 1.00 | 1.25 |

$kdz$ | 0.00 | 0.00 | 0.00 |

$\tau \varphi $ | 0.13 | 0.15 | 0.168 |

$\tau \theta $ | 0.173 | 0.20 | 0.224 |

$\tau z$ | 0.217 | 0.250 | 0.280 |

$\theta sp$ | 5.00 | 5.00 | 5.00 |

$zsp$ | 45.0 | 60.0 | 75.0 |

$a0$ | 15.0 | 15.0 | 15.0 |

$vlsat$ | 0.0677 | 0.0606 | 0.0524 |

### Comparison of Time-Dependent Quantities.

In order to compare individual quantities such as roll, yaw, or azimuth under crosswind flight, it is also necessary for the normalized period of the roll setpoint to be the same across all three model scales. By characterizing each model over the same range of *normalized* crosswind periods, we are able to compare the dynamic behavior during crosswind flight across all three model scales. An example of this comparison for a single normalized roll setpoint period is shown in Figs. 11–13. For each of the curves in Figs. 11–13, the average of the three curves was calculated, and the mean squared error from this averaged was subsequently calculated, in order to provide a signal variance that quantified the closeness of the curve to the mean. The signal variances for each of the experiments are reported in Table 6. These results demonstrate that the time-dependent, nondimensional quantities of roll angle, heading (yaw angle), and azimuth angle exhibit very similar amplitudes, phase shifts, frequencies, and nondimensionalized transient response characteristics across all scales.

### Comparison of Aggregate Quantities.

where *F _{b}* represents the net buoyant force. The net lift and drag forces,

*F*and

_{L}*F*, were estimated based on known characterizations of the lifting body. It is important to note that this method contains many simplifying assumptions, and therefore, should be considered as a coarse estimate for the purpose of illustration. Nonetheless, for the results shown here, the magnitude of the control power calculated via Eq. (30) rarely exceeded 5% of the power generated or about 25% of the excess power created through crosswind flight.

_{D}*T*, using an exponentially weighted regression function. An example form of this function for the small model at flow speed A is given by

*A*, $PAs={P\xaf1,As,P\xaf2,As,\u2026,P\xaf\eta ,As}$ is the set of corresponding measured mean power ratios given by Eq. (13), and $\eta \u2208\mathbb{N}$ is the number of periods tested for the small model at flow speed

*A*. The individual weight for a particular roll setpoint period,

*T*, $w(T|TAs)$, is given by

where *c* is a tuning parameter, chosen to be 5 $s\u22122$ in this work. Equivalent equations for the best curve for the medium and large models at specific flow speeds are given by changing the superscript *s* to *m* or *l* and the subscript *A* to *B* as appropriate. The raw data, along with the fit curves, are shown for both the mean power factor and 85 percentile power factor in Fig. 14.

*A*reads

where $T\u0303As={\alpha \u0303sT1,As,\alpha \u0303sT2,As,\u2026,\alpha \u0303sT\eta ,As}$ denotes the set of all normalized roll setpoint periods for the small model at flow speed *A* calculated using the variable scale factor $\alpha \u0303s$, and $T\xafAm={\alpha mT1,As,\alpha mT2,As,\u2026,\alpha mT\eta ,As}$ denotes the set of all normalized roll setpoint periods for the medium model at flow speed A calculated using the nominal scale factor *α ^{m}* calculated directly from the length scale, given in Eqs. (27) and (28). Note that since $\alpha m=1,\u2009T\xafjm=Tjm\u2009\u2009\u2009\u2009\u2009\u2200\u2009j\u2208{A,B}$.

Plots of the best fit curve for the mean power factor of the small, medium, and large models are shown by the darker curves in Fig. 14. The average power ratio is limited by frequent tacking, which arises due to tether length limitations in the water channel. In fact, the bottom plots of Fig. 14 show the 85th percentile power ratio, demonstrating that the system indeed achieves power ratios close to two over substantial portions of the crosswind motion. Maximizing the fraction of crosswind flight over which these larger power ratios are realized represents an important topic of ongoing control design research.

We compare the measured scale factors, given by Eq. (33), $\alpha \u0302s$ and $\alpha \u0302l$ to the predicted values, *α ^{s}* and

*α*, calculated directly from the characteristic lengths in Table 8. Figure 15 shows plots of the power ratio as a function of the normalized roll setpoint period, as calculated using both the predicted and measured scale factors.

^{l}Model size | Flow speed | Predicted scale factor, $\alpha $ | Experimentally identified scale factor, $\alpha \u0302$ | Error (%) |
---|---|---|---|---|

Small | B | 1.15 | 1.30 | 13 |

Large | B | 0.894 | 0.846 | 5.4 |

Small | A | 1.15 | 1.22 | 6.7 |

Large | A | 0.894 | 0.806 | 10 |

Model size | Flow speed | Predicted scale factor, $\alpha $ | Experimentally identified scale factor, $\alpha \u0302$ | Error (%) |
---|---|---|---|---|

Small | B | 1.15 | 1.30 | 13 |

Large | B | 0.894 | 0.846 | 5.4 |

Small | A | 1.15 | 1.22 | 6.7 |

Large | A | 0.894 | 0.806 | 10 |

### Discussion and Error Analysis.

Because the roll angle is controlled to a setpoint that is scaled according to the known characteristic lengths of the models, it is expected that the roll signals will overlap well when the horizontal axis of Fig. 11 is normalized using nominal scaling. Similarly, it is unsurprising that the yaw and azimuth angles exhibit the same normalized frequencies across different scales. The most significant result is the similarity in *amplitude* and general transient response characteristics of the uncontrolled signals, yaw, and azimuth. Because these quantities are nondimensional, their magnitudes should be independent of the model scale. This in fact conforms to the result shown in Figs. 12 and 13.

The dynamic scaling laws derived in this work and Ref. [23] are valid under the condition that *all* variables appearing the dimensionless groups are scaled precisely. If any of these variables is not scaled precisely, it can undermine the accuracy of the results. In this work, there are three possible sources of significant error, all of which relate to inherent limitations in the water channel setup.

First, imprecision in the flow speed measurement may have resulted in variation of the dimensionless groups across experiments. In this work, flow speed is measured by a point probe with a resolution of ±1.5 cm/s. Depending on which experiment is being conducted, this corresponds to a possible error in the flow speed of between 5.4% and 9.1%. The impact of this quantization-type error is particularly visible in the normalized flow speeds of Table 6.

Second, camera lens distortion effects can have a nontrivial effect on results. In particular, the approximations of Eq. (22) are accurate when the tracking targets appear near the center of the frame. Whereas lens distortion will cause them to become less accurate as the model deviates from the center of the frame, as is common with crosswind flight.

Finally, it is understood (and demonstrated in Refs. [21–23]) that the lateral stability of the model depends heavily on the achieved pitch angle. Generally, higher pitch results in greater lateral movement. Small inaccuracies in placement of the image tracking targets on the models can result in slight differences between the measured and achieved pitch angle across model scales. This in turn can lead to inconsistencies between model behavior at different scales.

## Conclusions

This paper presented a framework for lab-scale control and measurement of crosswind airborne wind energy systems. The framework is capable of simultaneous control of altitude, pitch, and roll, while measuring important variables upon which energy generation depends. The paper followed with a detailed discussion of crosswind flight performance across a range of model scales and flow speeds, including experimental validation of the dynamic scaling laws derived in this and previous work. Specifically, a dynamic model-based dimensional analysis from Ref. [23] suggested that properly scaled plant and controller parameters will lead to dynamic equivalence across all model scales during closed-loop control of crosswind flight. It was shown that the dynamics scaled properly across a range of available model scales, within experimental uncertainty.

## Funding Data

Entitled Collaborative Research: Self-Adjusting Periodic Optimal Control with Application to Energy-Harvesting Flight, National Science Foundation (Grant No. 1538369).