0
Research Papers

# Experimentally Infused Plant and Controller Optimization Using Iterative Design of Experiments—Theoretical Framework and Airborne Wind Energy Case StudyOPEN ACCESS

[+] Author and Article Information
Nihar Deodhar

Department of Mechanical Engineering,
University of North Carolina at Charlotte,
Charlotte, NC 28223
e-mail: ndeodhar@uncc.edu

Joseph Deese

Department of Mechanical Engineering,
University of North Carolina at Charlotte,
Charlotte, NC 28223
e-mail: jdeese23@uncc.edu

Christopher Vermillion

Mem. ASME
Department of Mechanical Engineering,
University of North Carolina at Charlotte,
Charlotte, NC 28223
e-mail: cvermill@uncc.edu

1Corresponding author.

Contributed by the Dynamic Systems Division of ASME for publication in the JOURNAL OF DYNAMIC SYSTEMS, MEASUREMENT, AND CONTROL. Manuscript received December 20, 2016; final manuscript received June 1, 2017; published online August 29, 2017. Assoc. Editor: Dumitru I. Caruntu.

J. Dyn. Sys., Meas., Control 140(1), 011004 (Aug 29, 2017) (10 pages) Paper No: DS-16-1603; doi: 10.1115/1.4037014 History: Received December 20, 2016; Revised June 01, 2017

## Abstract

This research presents an iterative framework for optimizing the plant and controller for complex systems by fusing expensive but valuable experiments with cheap yet less accurate simulations. At each iteration, G-optimal design is used to generate experiments and simulations within a prescribed design space that is shrunken in size after each successful iteration. The shrinking of the design space is determined through statistical characterization of a response surface model, and further shrinking is achieved at successive iterations through a numerical model correction factor that is driven by the results of experiments. An initial validation of this iterative design optimization framework was performed on an airborne wind energy (AWE) system, where tethers and an aerostat are used in place of a tower to elevate the turbine to high altitudes. Using a unique lab-scale setup for the experiments, the aforementioned iterative methodology was used to optimize the center of mass location and pitch angle set point for the airborne wind energy system. The optimum configuration yielded a substantial improvement in system responses as compared to a numerically optimized configuration. The framework was recently extended to include four variables (horizontal and vertical stabilizer areas, center of mass location, and pitch angle set point).

<>

## Introduction

In numerous controlled systems, the optimal controller depends on the physical system design (i.e., the plant) and vice versa. In response to this reality, a significant body of literature has emerged in the past two decades on the combined plant and controller design, sometimes termed the codesign, of such systems. For example, a nested combined plant and controller optimization has been presented for an active suspension system in Refs. [1] and [2], whereas Ref. [3] uses the same techniques for optimization of an elevator. The combined plant and controller optimization is achieved in Refs. [4] and [5] using a sequential design approach, along with a control proxy function, to achieve optimal or near-optimal designs. References [6] and [7] use a decomposition-based codesign approach, whereas Ref. [8] uses direct transcription to tackle the combined plant and controller design problem.

The aforementioned works differ in their approaches to codesign but are similar in the sense that they all rely on a dynamic model of the system at hand to perform their optimization, rather than integrating experiments into the design approach. A more limited set of work, most notably [9], has demonstrated the use of experiments for plant and controller optimization. The work of Ref. [9] is entirely experimentally driven, using experiments from one iteration to identify elements of a Jacobian matrix that is used to determine the direction that the design variables need to be varied at the following iteration.

Airborne wind energy (AWE) systems, which are the focus of the experiments and simulations in this work, feature complex dynamics and coupling between the plant and the controller. These systems replace conventional towers with tethers and a lifting body (an aerostat, kite, or rigid wing). Replacement of traditional tower-mounted turbines with flexible tethers and a lifting body has enabled wind energy generation at altitudes up to 600 m [10]. Over the past decade, AWE technology has moved from paper to practice, with several companies and research organizations pioneering a variety of AWE systems. Google-owned Makani Power [11] and Altaeros Energies [12] have developed AWE systems with on-board wind energy generation. Ampyx Power [13], Kitegen [14], and Windlift [15] represent just a few examples of organizations that rely on cyclic motions (high-tension crosswind spool-out and low-tension spool-in) and ground-based energy generation.

A substantial body of research focuses on control design for AWE systems. For example, the authors in Refs. [16] and [17] have aimed to optimize energy production for kite-based systems using crosswind flight patterns. Control design for the Altaeros Buoyant Airborne Turbine (BAT) is introduced in Ref. [18]. Net power generation for the same system is maximized in Ref. [19] by determining instantaneous optimal altitude using a switched extremum seeking algorithm. A much smaller body of literature addresses the design of both the AWE system plant and the controller, with Refs. [18] and [20] being some recent examples of such work.

Most of the literature on AWE systems focuses on numerical simulations, with some validation results presented on full-scale systems. However, AWE systems fall under the broad category of complex nonlinear systems with imperfect numerical models. Although AWE numerical models are helpful in informing control and plant design decisions, different phenomena such as wake effects, vortex shading, added mass effects, and tether/fluid interaction exist in reality, yet are not captured by relatively low-order numerical models. Consequently, it is very desirable for AWE systems to fuse numerical modeling tools with inexpensive experimental tools for optimizing system performance. This fusion of plant and controller optimization has been absent in the codesign literature until recently.

In response to the desire to fuse experiments with numerical models in optimizing AWE system designs, the authors of Refs. [2124] have developed and exploited a lab-scale setup to characterize different AWE plant and control system designs. The system has evolved from a passive test setup at the University of Michigan to a controlled test setup at UNC Charlotte. In the most recent version of this lab-scale setup, AWE lifting body models are three-dimensional (3D) printed and “flown” in a water channel, cameras are used for real-time motion capture, direct current (DC) motors are used to regulate tether lengths, and control algorithms are implemented on a high-performance target computer. The lab-scale AWE prototyping setup at UNC Charlotte enables the fusion of experiments with simulations in the optimization of AWE system controllers and plant designs. Deodhar et al. [25] present an initial codesign result using the lab-scale setup, taking the Altaeros BAT as a case study. The approach followed in Ref. [25] is local in nature. It uses experiments in the vicinity of the estimated optimum design to determine the gradient of performance with respect to design variables, then moves along the direction of that gradient in subsequent iterations. While the initial results of a two-variable optimization in Ref. [25] were positive, this gradient-based approach can easily result in convergence to local optima in the presence of highly nonlinear systems with complex dynamics, especially when more design variables are considered. To address this, our recent conference paper, Ref. [26], proposed an alternative framework that uses optimal design of experiments (DOE) to choose control and plant design variables that span the design space, then uses a statistical quality of fit characterization to shrink the design space at the next iteration, as well as a correction term to improve the numerical model prediction at the next iteration. By ensuring coverage of the design space at each iteration, the optimization approach is much less likely to get stuck in local optima.

While the contribution of Ref. [26] presented a theoretical framework, it lacked actual experimental validation. In this paper, we build upon our initial work done in Ref. [26] by validating the framework with actual experiments and modifying the framework to handle larger numbers of optimization variables. The results presented herein focus on the Altaeros BAT [12] as a case study.

The contributions of this paper are threefold:

1. (1)To validate the framework presented in Ref. [26] with physical experiments performed in the UNC-Charlotte lab-scale AWE test setup.
2. (2)To extend the numerical model and optimization framework to include additional design variables (specifically, stabilizer areas, which now exist as explicit design variables in the numerical model).
3. (3)To validate the effectiveness of the design optimization approach for larger numbers of design variables, using an increased-fidelity model as a surrogate for experiments.

## Framework for Experimentally Infused Optimization With Iterative Design of Experiment

The work presented herein represents a departure from the vast majority of codesign literature, which assumes that an accurate closed-form representation for the system's dynamic equations is available. The availability of a closed-form model allows for a variety of numerical tools to be used to optimize the plant and control system design variables, denoted by $pp$ and $pc$, respectively. Our work, on the other hand, assumes that while a model of the system dynamics exists, it is imperfect. Thus, converging upon an optimal set of design variables requires fusion of experiments with numerical modeling. The design optimization problem that we are pursuing in this work can be formally summarized as follows: Display Formula

(1)$Minimize J(pp,pc)=∫0Tfl(x(t),u(t);pp,pc)dt$
Display Formula
(2)$subject to:x˙=ftrue(x,u)$
Display Formula
(3)$pp∈P$
Display Formula
(4)$pc∈C$

Here, x and u represent the state and control input vectors, respectively, and Eq. (2) represents the true system dynamics, as represented by experimental results. It is assumed that the mathematical structure of ftrue (x, u) is not known, but an approximate numerical model for the dynamics exists. The constraint sets P and C represent hard constraints on the plant and controller design variables, respectively.

Because the true system dynamics are not fully known but instead can only be revealed through experiments, the infusion of numerical modeling tools with experimental results is necessary to address the optimization problem of Eq. (1). In particular, two broad questions may be asked regarding any approach that aims to tackle the optimization problem of Eq. (1) through the fusion of experiments:

1. (1)Which experiments should be run in order to reveal the most valuable information for minimizing $J(pp,pc)$?
2. (2)How can the numerical model of the system be improved after each run of experiments in order to facilitate convergence of $pp$ and $pc$ to their optimal values?

We address the aforementioned questions through an iterative optimization process. This process, termed experimentally infused optimization, involves the generation of experiments that span a design space, followed by the reduction of that design space and improvement of the numerical model. The entire process is divided into four main parts as follows:

1. (1)Perform an optimal DOE that generates a set of design variables to be characterized through numerical simulations and another set of design variables to be characterized through experiments.
2. (2)Characterize a response surface from numerical simulations and experiments.
3. (3)Learn from experiments to improve the numerical model for the next iteration.
4. (4)Reduce the size of the design space based on the quality of the response surface approximation, then repeat the process.

A flowchart of this process is presented in Fig. 1. Each of the aforementioned processes is explained in detail in this section.

###### Optimal Design of Experiment.

The overarching objective of DOE is to design experiments over a prescribed design space that maximizes the amount of information gathered by a fixed number of experiments. In this work, G-optimality was used in designing two sets of designs, namely, the numerical simulations and experiments, where NS and NE represent the number of simulations and experiments, respectively. In the proposed framework, each of the NS and NE simulations and experiments is comprised of a different set of controller and plant design variables; hence, the G-optimal design explores the control and plant design spaces simultaneously. G-optimal design chooses design points in an attempt to minimize the maximum prediction variance across all experiments. It is known to do an effective job of filling the design space without repeating experiments. G-optimal design relies on the specification of an m element column vector known as the regressor vector, $r(pp,pc)$, which captures the anticipated dependencies of the system performance on the design variables. For example, for a system with a single scalar plant design variable and control design variable to be optimized, a regressor vector of $r(pp,pc)=[pppp2pcpc2]T$ would be chosen if the system designer anticipated that the performance was quadratic with respect to both the plant and controller design variable.

The matrix of n regressor vectors can be expressed as Display Formula

(5)$R=[r1 r2… rn]$

Here, n corresponds to the number of design points (either numerical simulations or experiments). G-optimality criteria minimize the maximum prediction variance across all the possible design points in the design space, M. The resulting design points are given by Display Formula

(6)$[pp,pc]=argmin(pp,pc)∈MriT(RRT)−1ri$

Here, ri denotes the regressor vector corresponding to the ith design point.

At every iteration, the NE experimental design points are split into an equal number of training and validation points. In this work, the experimental points were sorted based on their distance from each other, and alternate points were used for training and validation to get better coverage of the design space. The set of training, validation, and simulation points resulting from the G-optimal DOE are denoted by $ptrain, pval$, and $psim$, respectively.

###### Numerical Simulations and Experiments.

Experiments and numerical simulations are conducted for every $ptrain, pval$, and $psim$ using identical operating conditions (i.e., if a controlled perturbation is applied in numerical simulations, then that same controlled perturbation is applied in experiments). An identical performance index, denoted by J(p), is used to characterize the resulting performance under all experiments and simulations.

###### Response Surface Characterization and Numerical Model Correction.

Three different response surfaces are characterized at every iteration of the optimization process, specifically:

1. (1)Pure experimental response using both training and validation experiments, $Ĵexp$;
2. (2)Pure numerical simulations response, $Ĵsim$; and
3. (3)Combined experimental training and numerical simulations response, $Ĵ$.

Each response surface characterization takes the data from experiments or simulations and then uses this data to fit a function over the entire design space. Two techniques are proposed in this work for that response surface characterization.

###### Parametric Response Surface Characterization Based on a Regressor Vector.

In cases where the general structure of the dependencies of performance on the design variables is known, and in cases where the number of design variables is small, it is often possible to relate the performance to the design variables through a regressor vector. Coefficients of the regressor vector can then be determined through weighted least squares regression.

Using this approach, the response surface is estimated using the expression Display Formula

(7)$Ĵ(p)=r(p)Tb$

where $Ĵ(p)$ is the estimated performance at any point p, and b is the vector of regression coefficients obtained by minimizing the weighted sum of squared errors given by Display Formula

(8)$b=argminb∈IRn∑i=1nw(pi)(J(pi)−r(pi)Tb)2$

Here, $w(pi)$ is the weight on ith point. These weights can be used to value simulations and experiments differently in arriving at the combined response surface, $Ĵ$.

###### Distance-Weighted Response Characterization.

With a large number of design variables, and in instances when it is difficult to speculate a regressor vector based on known properties of the system, the previous approach may result in poor characterization of the response surface. A more general distance-weighted response surface characterization can yield a more accurate representation of performance over the design space in these cases. In this approach, the performance associated with any point p is modeled as the weighted sum of all measured points, placing more weight on the points that are closest to p. The response surface is approximated as follows: Display Formula

(9)$Ĵ(p)=∑i=1nW(pi)J(pi)∑i=1nwi$
Display Formula
(10)$W(pi)=1d(p,pi)2$

where $d(p,pi)$ is the Euclidean distance between p and $pi$.

###### Numerical Model Improvement.

In order to continue to reduce the size of the design space at each iteration, the numerical model must be improved at each iteration. This improvement is accomplished by comparing the response surface generated with experiments to the one generated with the numerical model, then adjusting the numerical model so that it conforms more closely with experimental results at the next iteration. Once the experimental ($Ĵexp$) and numerical ($Ĵsim$) response surface characterizations are obtained, a correction term ($Ĵcor$) is applied to numerical simulations. The form for this correction, at iteration i, is given as follows: Display Formula

(11)$Ĵnum(p)=Ĵsim(p)+Ĵcor,i(p)$
Display Formula
(12)$Ĵcor,i(p)=k(Ĵ exp (p)−Ĵsim(p))+∑j=1i−1Ĵcor,i(p)$

Here, $Ĵcor,i(p)$ represents the numerical correction at step i, which is obtained by adding a correction from the current iteration to the cumulative correction from all previous iterations (which is reflected in $Ĵsim$ for the current iteration). The variable k, which is taken to be between 0 and 1, limits the correction to a certain percentage of the gap between experiments and the numerical model at each iteration.

###### Design Space Reduction.

The next goal is to converge close to the optimum design point p* by rejecting all candidate configurations that are decidedly suboptimal. Determination of what is “decidedly suboptimal” is achieved using hypothesis testing.

The job of the design space reduction is to determine whether one can conclude, for each candidate design within the design space, that the candidate design cannot possibly (within some confidence interval) be the true optimum. If the candidate design cannot possibly be optimal, it is eliminated from the design space; otherwise, it is retained. For these purposes, the null (H0) and alternate (H1) hypotheses are defined as Display Formula

(13)$H0:J(p)=Ĵ(p*)$
Display Formula
(14)$H1:J(p)>Ĵ(p*)$

The z test is used to reject a candidate configuration by comparing the corresponding z-score to a threshold z0 for a given confidence interval Display Formula

(15)$z(p)=Ĵ(p)−Ĵ(p*)S1n1+1n2$
Display Formula
(16)$S=∑i=1NE/2(Ĵ(pval)−J(pval))2NE2−1$

The term S is referred to as the quality of fit between combined response at validation points and the true experimental response at validation points. n1 and n2 denote the number of times each experiment/simulation is repeated.

## Case Study on an Airborne Wind Energy System

For the case study in this paper, the Altaeros BAT is considered, for which the system model and controller are detailed in this section. This system elevates a horizontal axis turbine to high altitudes through the use of a lighter-than-air annular aerostat. The system is intended to operate in a stationary fashion, and it is essential that the system maintains its altitude and attitude (orientation), while keeping control inputs (tether motions) to a minimum, in order to maximize its energy production.

###### Dynamic Model of Airborne Wind Energy System.

The dynamic model of Altaeros system is described in detail in Ref. [22] but is also summarized here for purposes of self-containment. Figure 2 shows the relevant variables used in the dynamic model.

The three tethers are approximated as a single tether with length Lt plus a bridle joint at the top that splits in three attachment points on the BAT. The single tether model circumvents the need to solve differential algebraic equations and instead, results in a set of ordinary differential equations that describe the equations of motion. This effect multiple tethers is approximated by a two degrees-of-freedom bridle joint, whose angular degrees-of-freedom are given by $ϕ′$ and $θ′$. We refer to $ϕ′$ and $θ′$ as tether-induced angles. These angles are approximated by following expressions, which relate the actual individual tether lengths to $ϕ′$ and $θ′$: Display Formula

(17)$ϕ′=tan−1(l3−l2lseplat)$
Display Formula
(18)$θ′=tan−1(l1−0.5(l2+l3)lseplong)$

where $lseplong$ and $lseplat$ are longitudinal and lateral tether attachment separation distances, respectively. l1, l2, and l3 are the lengths of the individual tethers. The center of mass location is a function of Φ (zenith angle), Θ (azimuth angle), and Lt (average tether length). The control inputs are the tether release speeds, $u¯i,$ which are related to the released tether lengths by Display Formula

(19)$ddtli=u¯i$

The governing equations of motion below are derived using the Euler–Lagrange formulation Display Formula

(20)$D(Q)Q¨+C(Q,Q˙)Q˙+g(Q)=τ(Q,Q˙,Vflow,ψflow)$
Display Formula
(21)$X=f(Q,Q˙)$
Display Formula
(22)$Ω=g(Q,Q˙)$

where Display Formula

(23)$Q=[Φ Θ Ψ Lt θ′ ϕ′]$
Display Formula
(24)$X=[xg yg zg u v w]$
Display Formula
(25)$Ω=[ϕ θ ψ p q r]$

Here, Vflow is the flow (wind in the case, of the full-scale system, water in the case of the lab-scale water channel-based system) speed, and ψflow and τ represent flow direction and vector of generalized external forces. Aerodynamic forces and moments are functions of α (angle of attack) and β (side slip angle). Each of the variables in Eq. (23) is summarized in Nomenclature section.

###### Closed-Loop Controller Structure.

The flight controller for the BAT, first described in Ref. [18], uses the multiple tethers to track three different set points, namely, altitude (zsp), pitch (θsp), and roll ($ϕsp$). Figure 3 presents a block diagram of this controller structure. The controller first computes “virtual” control inputs, denoted by $v¯z, v¯θ$, and $v¯ϕ$, which reflect the need to adjust all three tethers synchronously to change altitude, adjust the forward and aft tether asynchronously to change the pitch angle, and adjust the port and starboard tethers asynchronously to adjust the roll angle. These virtual control inputs are the output of filtered proportional plus derivative controllers (lead filters). Because the plant itself contains pure integrators (i.e., a “type 1” system) in the dynamics between the control inputs ($u¯center, u¯stbd$, and $u¯port$) and the variables to be controlled (z, θ, and $ϕ$), integral control is not necessary for tracking in this case, since the plant includes three pure integrators between the motor speeds and tether lengths. The three control inputs represent tether release speeds, which are controlled in the full-scale system through alternating current induction models and in the lab-scale, water channel-based system by micro DC motors.

The average tether release speed ($v¯z$), difference between central/peripheral tether release speeds ($v¯θ$), and port/starboard tether release speeds ($v¯ϕ$) are related to the tracking errors in altitude (ze), pitch (θe), and roll ($ϕe$) by the transfer functions Gz(s), Gθ(s), and $Gϕ(s)$, respectively. These transfer functions are given by the following expressions: Display Formula

(26)$Gz(s)=kdzs+kpzτzs+1$
Display Formula
(27)$Gθ(s)=kdθs+kpθτθs+1$
Display Formula
(28)$Gϕ(s)=kdϕs+kpϕτϕs+1$

Legacy work by Vermillion et al. (see Ref. [18]) discusses this control structure in more detail, along with mechanisms for calibrating the feedback control parameters. Scaling to the lab-scale environment is discussed in Ref. [22]. Given the symmetrical design of the BAT and the desire for stationary operation, $ϕsp$ is always equal to zero. The tether release speeds ($u¯center, u¯stbd$, and $u¯port$) are related to the virtual control inputs ($v¯z, v¯θ$, and $v¯ϕ$) as follows: Display Formula

(29)$u=Av$
Display Formula
(30)$A=[1 −1 01 1 11 1 −1]$

Here, u represents the vector of control inputs ($u¯center, u¯stbd$, and $u¯port$), whereas v denotes the vector of virtual control inputs ($v¯z, v¯θ$, and $v¯ϕ$).

## Airborne Wind Energy Case Study—Two Design Variables

In this section, we show how the numerical model of Sec. 3 has been fused with lab-scale experiments in order to optimize design variables for the BAT. In this section, we show the results of the experimentally infused optimization framework for two design variables, namely, the longitudinal location of center of mass (xcm) and the trim pitch angle (θsp). In Sec. 5, we extend the framework to four design variables.

###### Experimental Setup of the Lab-Scale System.

Although it is ultimately necessary to validate the plant and control design on a full-scale AWE system, extensive full-scale experiments are expensive and time-consuming. Furthermore, the consequences of a failure are significant, making it inadvisable to explore design configurations that are on the verge of instability for the full-scale system. For these reasons, a lab-scale setup consisting of water channel of 1 m × 1 m cross section and 1/100th model of the BAT was developed. This setup evolved from an initial passive system at the University of Michigan (described in Ref. [21]) to the present actively controlled system at UNC Charlotte (initially described in Ref. [23]). The present setup is shown in Fig. 4 and has four major components:

1. (1)Three high-speed motion capture cameras for tracking the position and orientation of the BAT,
2. (2)DC motors to independently control the lengths of three tethers,
3. (3)A high-performance target computer for intense real-time computations and feedback control, and
4. (4)A host computer to communicate data between the user and the target computer.

The BAT model is 3D printed out of acrylonitrile butadiene styrene plastic and shown in Fig. 5. A drag screen slit accommodates a mesh screen that replicates the drag properties of a spinning rotor. The center of mass location and net buoyancy can be adjusted by placing lead wires in selected ballast holes. Full position and orientation information is resolved from the location of three sets of white dots on the model, projected onto the field of view of three cameras. Specifically, the cameras capture instantaneous pixel locations of the center of every pair of dots, and, based on the body-frame locations of the dots, the Euler angles and center of mass locations are computed. The technical details of this process are highlighted in Ref. [24].

Finally, a separate dimensional analysis was conducted by the authors in Refs. [21] and [22], demonstrating the dynamic equivalence between the lab-scale and full-scale system, with the exception of uniformly accelerated time constants in the water channel. Specifically, this analysis reveals that, under the assumption of consistent aerodynamic (hydrodynamic) coefficients for the lab-scale and full-scale model (these are not guaranteed due to a Reynolds number mismatch, but have been shown through wind tunnel testing at different Reynolds numbers to be similar), dynamic similarity can be achieved under the following conditions:

1. (1)Maintaining geometric similarity between the lab-scale and full-scale system (requires all lengths to be scaled according to Lt),
2. (2)Maintaining a consistent net buoyancy between the lab-scale and full-scale system (requires the quantity $m/ρLt3$ to be maintained between the lab-scale and full-scale system), and
3. (3)Maintaining a consistent Froude number between the lab-scale and full-scale system (requires the flow velocity, Vflow to be proportional to $Lt$).

Following these scaling rules, with consistent aerodynamic (hydrodynamic) coefficients between scales, the dynamics of the lab-scale system will be dynamically similar to those of the full-scale systems, with time constants that are accelerated according to the square root of the length scale ($τ∝Lt$).

###### Performance Index.

Optimizing the performance of the BAT via the techniques of Sec. 2 requires us to quantify performance. For the BAT, this was done through a cost function, J, which we desire to minimize. The cost function used for this case included six terms that penalize tracking errors (on roll, pitch, and altitude), ground footprint (using zenith angle) and energy consumption (through the control inputs). The cost function structure is given as follows: Display Formula

(31)$J=∫TiTf(k1(ψ−ψflow)2+k2(θe)2+k3(ϕe)2+k4(ze)2+k5(uTu)+k6(Φ)2) dt$

Here, k1,…, k6 are the weights on each term in the cost function. These weights were chosen to emphasize some components of the cost over other and account for the different engineering units in each of the terms.

###### Experimental Profile and Perturbations.

Every training and validation experiment was conducted using identical flow speeds and lateral perturbations to the ones used in numerical simulations. Specifically, every experiment/simulation was conducted for a period of 300 s, wherein three different flow speeds of 0.205, 0.245, and 0.285 m/s were used over each 100 s interval. Twenty seconds after each speed change, the BAT was subjected to a sequence of lateral perturbations, which were induced by pulling down on the port tether at full speed for 1 s, followed by a 1 s delay before resuming closed-loop control. This perturbation profile can be represented as follows: Display Formula

(32)$u¯port={−umax,if t∈Tu¯z−u¯θ−u¯ϕ,otherwise$
Display Formula
(33)$T=[20 s 22 s]∪[120 s 122 s]∪[220 s 222 s]$

###### Case Study Results.

The framework in Sec. 2 was used to optimize center of mass location and the trim pitch angle. The weights w(p) given by Eq. (8) were chosen as follows: training experiments $w(ptrain)$ had a weight of 80% as compared to only 20% for numerical simulations $w(psim)$ in the response surface characterization. For response surface characterization, we used a quadratic response surface for which the regressor vector is given by Display Formula

(34)$r=[1 pp2 pp pc2 pc]T$

A 90% confidence interval was used for determining whether or not to reject points from the design space based on the z-score of Eq. (15).

The results of three iterations of the experimentally infused optimization process, along with the optimal configuration given purely by the numerical model, are presented in Table 1. The left half of Fig. 6 shows the G-optimal experiments and numerical simulations generated at every iteration, and the right half shows the corresponding $Ĵ$ at for every combination of design variables in the new design space. It is evident from the table and figures that the design space is shrinking between every iteration. Thereafter, the variation in experiments is dictated by the noise in the experimental data, which limits further shrinking of the design space. The shrinking at very first iteration is not as significant as that of the second iteration since the numerical model remains uncorrected during the first iteration. At the onset of second iteration, the numerical model is corrected using experiments, which leads to a better quality of fit between combined numerical plus training data and the validation data. Figure 7, which compares the quality of flight under purely numerical optimization to the quality of flight after three iterations of the experimentally infused process, shows a significant improvement as a result of the incorporation of experiments.

## Extended Optimization Framework for Four Design Variables

Figure 7 shows that a notable improvement in the AWE system performance results from the optimization of just two design variables. However, to further improve the system response, this section shows how two design variables can be added to the previous ones to improve the performance even more. In particular, this section focuses on the optimization of four design parameters, namely, the longitudinal center of mass location, trim pitch angle, horizontal stabilizer area, and vertical stabilizer area. Thus, the additional design variables chosen were:

1. (1)Horizontal stabilizer surface areas (AH) and
2. (2)Vertical stabilizer surface areas (AV)

The pairs of horizontal and vertical stabilizers are critical in reducing heading oscillations. In order to characterize the impact of varying stabilizer surface areas, the model presented in Sec. 3 was modified such that the aerodynamic force and moment coefficients for the overall body were partitioned between the main body (the “fuselage”) and the stabilizers. These equations in the modified numerical model are given by Display Formula

(35)$CD,L,S(α,β)=CD,L,SF(α,β)+CD,L,SH(α,β)AHAref+CD,L,SV(α,β)AVAref$
Display Formula
(36)$CMx,My,Mz(α,β)=CMx,My,MzF(α,β)+CMx,My,MzH(α,β)AHlHAreflref+CMx,My,MzV(α,β)AVlVAreflref$

where CD, CL, and CS are the drag, lift, and side force coefficients, whereas $CMx, CMy$, and $CMz$ are the roll, pitch, and yaw moment coefficients. α and β are angles of attack and the side slip angle, and Aref and lref are the reference area and reference length for the full system. The variables AH and AV represent the horizontal and vertical stabilizer surface areas, respectively, which serve as design variables to be optimized. The parameters lH and lV represent reference lengths for the horizontal and vertical stabilizers, respectively. These parameters were adjusted such that a constant aspect ratio was maintained for any combination of stabilizer surface areas (AH and AV, respectively).

###### Numerical Alternative to Experiments.

For our initial optimization of four design variables, we used an increased-fidelity model as a surrogate for expensive experiments. The alternate model incorporated stall effects for the fuselage and the horizontal stabilizer that were omitted from the numerical model. Previous wind tunnel testing on a $1/20th$ scaled BAT model in Ref. [21] showed that the body and horizontal stabilizers would stall beyond a particular angle of attack. To replicate these stall effects, the lift coefficients on body and stabilizer were modeled as follows: Display Formula

(37)$CLF={CL0F+CLαFαF,if |αF|≤αstallFCL0F+2CLαFαstallF−CLαFαF,otherwise$
Display Formula
(38)$CLH={CL0H+CLαHαH,if |αH|≤αstallHCL0H+2CLαHαstallH−CLαHαH,otherwise$

Here, $CLF(αF)$ and $CLH(αH)$ denote the lift coefficients for main body and the horizontal stabilizer, respectively; $CL0F$ and $CL0H$ are the main body and horizontal stabilizer lift coefficients at zero angle of attack; $CLαF$ and $CLαH$ represent the linear lift coefficients; αF and αH are the angles of attack; and $αstallF and αstallH$ are the angles of attack at which the main body and horizontal stabilizers stall.

###### Results Using Alternative Numerical Model as Experiments.

The experimentally infused optimization framework was applied to modified numerical model and surrogate experiments using identical test profile as in the case of a two-variable optimization. Unlike previous weights of 80/20 on experiments versus numerical simulations, the weight on the experiment was reduced to 60% since the actual experiments were replaced by a surrogate model. The results are shown in Table 2, where performance is improved at each iteration and the size of the design space is reduced. The zeroth iteration represents the results of numerical optimization alone. The comparison of system responses for optimal design given at zeroth iteration (i.e., the purely numerically optimized system) and third iteration is shown in Fig. 8. The figure shows that the purely numerically optimized design yields a limit-cycle behavior after the first perturbation. It can be seen that the experimentally infused optimization process produces a design that performs significantly better than the design that is optimized based on the simple numerical model alone. A tradeoff between quality of flight given by tracking errors and ground footprint given by the zenith angle is observed in the case of four-variable optimization (also seen in Fig. 9). The quality of flight is improved greatly at the cost of a slight increase in the ground footprint.

## Conclusions and Future Work

A formal iterative method for combined plant and controller optimization, fusing both simulations and experiments, was proposed in this research. The framework was validated using an AWE application, initially using the center of mass location and pitch angle set point as optimization variables. By applying the experimentally infused optimization process, using a unique lab-scale setup for the experiments themselves, we demonstrated that we can obtain much better system performance through the infusion of experiments than with the numerical model alone. We also extended the results of our framework to include additional design variables, adding horizontal and vertical stabilizer surface areas to the list of design variables to be optimized. We performed an initial demonstration of the effectiveness of our approach for larger numbers of design variables (four, in this case), using a higher-fidelity “surrogate experimental” model in place of experiments. Future work will replace the surrogate experimental model with actual experiments using rapidly reconfigurable 3D printed models that enable us to “plug-and-play” different stabilizer geometries. In addition, the controller set points will be optimized in real time based on flow conditions and power production. The stability and robustness of the experimentally infused framework will be assessed through different convergence criteria.

## Acknowledgements

The authors would like to thank the team members at UNC Charlotte who participated in this research.

## Funding Data

• National Science Foundation (1453912).

## Nomenclature

• Lt =

average tether length

• p =

angular velocity about body fixed x axis

• q =

angular velocity about body fixed y axis

• r =

angular velocity about body fixed z axis

• u =

body frame x velocity

• v =

body frame y velocity

• w =

body frame z velocity

• xg =

ground frame x position

• yg =

ground frame y position

• zg =

ground frame z position

• θ =

Euler pitch angle

• Θ =

azimuth angle

• $θ′$ =

tether-induced pitch angle

• $ϕ$ =

Euler roll angle

• $ϕ′$ =

tether-induced roll angle

• ψ =

Euler yaw angle

• Ψ =

zenith angle

• Ψ =

“twist” angle (around tether)

## References

Fathy, H. K. , Papalambros, P. Y. , Ulsoy, A. G. , and Hrovat, D. , 2003, “ Nested Plant/Controller Optimization With Application to Combined Passive/Active Automotive Suspensions,” American Control Conference (ACC), Denver, CO, June 4–6, pp. 3375–3380.
Fathy, H. K. , Papalambros, P. Y. , and Ulsoy, A. G. , 2003, “ Integrated Plant, Observer, and Controller Optimization With Application to Combined Passive/Active Automotive Suspensions,” ASME Paper No. IMECE2003-42014.
Fathy, H. , Bortoff, S. , Copeland, S. , Papalambros, P. , and Ulsoy, A. , 2002, “ Nested Optimization of an Elevator and Its Gain-Scheduled LQG Controller,” ASME Paper No. IMECE2002-39273.
Peters, D. L. , Papalambros, P. Y. , and Ulsoy, A. G. , 2013, “ Sequential Co-Design of an Artifact and Its Controller Via Control Proxy Functions,” J. Mechatronics, 23(4), pp. 409–418.
Peters, D. L. , Papalambros, P. Y. , and Ulsoy, A. G. , 2011, “ Control Proxy Functions for Sequential Design and Control Optimization,” ASME J. Mech. Des., 133(9), p. 091007.
Alexander, M. J. , Allison, J. T. , and Papalambros, P. Y. , 2012, “ Decomposition-Based Design Optimization of Electric Vehicle Powertrains Using Proper Orthogonal Decomposition,” Int. J. Powertrains, 1(1), pp. 72–92.
Allison, J. T. , and Nazari, S. , 2010, “ Combined Plant and Controller Design Using Decomposition-Based Design Optimization and the Minimum Principle,” ASME Paper No. DETC2010-28887.
Allison, J. T. , Guo, T. , and Han, Z. , 2014, “ Co-Design of an Active Suspension Using Simultaneous Dynamic Optimization,” ASME J. Mech. Des., 136(8), p. 081003.
Pil, A. C. , and Asada, H. H. , 1996, “ Integrated Structure/Control Design of Mechatronic Systems Using a Recursive Experimental Optimization Method,” IEEE/ASME Trans. Mechatronics, 1(3), pp. 191–203.
Archer, C. L. , Monache, L. D. , and Rife, D. L. , 2014, “ Airborne Wind Energy: Optimal Locations and Variability,” Int. J. Renewable Energy, 64, pp. 180–186.
MAKANI, 2006, “ Energy Kites,” Makani Power, Alameda, CA, accessed June 15, 2017,
Altaeros Energies, 2010, “ The Next Transformation in Rural Infrastructure Efficiency,” Altaeros Energies, Somerville, MA, accessed June 15, 2017,
Ampyx Power, 2008, “ The Energy of Tomorrow,” Ampyx Power, The Hague, The Netherlands, accessed June 15, 2017,
Kitegen, 2007, “ Energy is in the Air,” Kitegen Research, Torino, Italy, accessed June 15, 2017,
Windlift, 2017, “ Airborne Wind Energy Systems,” Windlift, Raleigh, NC, accessed June 15, 2017,
Fagiano, L. , Milanese, M. , and Piga, D. , 2012, “ Optimization of Airborne Wind Energy Generators,” Int. J. Robust Nonlinear Control, 22(18), pp. 2055–2083.
Williams, P. , Lansdorp, B. , and Ockels, W. , 2008, “ Optimal Crosswind Towing and Power Generation With Tethered Kites,” J. Guid., Control, Dyn., 31(1), pp. 81–93.
Vermillion, C. , Grunnagle, T. , Lim, R. , and Kolmanovsky, I. , 2014, “ Model-Based Plant Design and Hierarchical Control of a Prototype Lighter-Than-Air Wind Energy System, With Experimental Flight Test Results,” IEEE Trans. Control Syst. Technol., 22(2), pp. 531–542.
Bafandeh, A. , and Vermillion, C. , 2016, “ Real-Time Altitude Optimization of Airborne Wind Energy Systems Using Lyapunov-Based Switched Extremum Seeking Control,” American Control Conference (ACC), Boston, MA, July 6–8, pp. 4990–4995.
Nikpoorparizi, P. , Deodhar, N. , and Vermillion, C. , 2016, “ Combined Plant and Controller Performance Analysis and Optimization for an Energy-Harvesting Tethered Wing,” American Control Conference (ACC), Boston, MA, July 6–8, pp. 4089–4094.
Vermillion, C. , Glass, B. , and Greenwood, S. , 2014, “ Evaluation of a Water Channel-Based Platform for Characterizing Aerostat Flight Dynamics: A Case Study on a Lighter-Than-Air Wind Energy System,” AIAA Paper No. 2014-2711.
Vermillion, C. , Glass, B. , and Szalai, B. , 2014, “ Development and Full-Scale Experimental Validation of a Rapid Prototyping Environment for Plant and Control Design of Airborne Wind Energy Systems,” ASME Paper No. DSCC2014-5907.
Deese, J. T. , Muyimbwa, T. , Deodhar, N. A. , Vermillion, C. R. , and Tkacik, P. , 2015, “ Lab-Scale Characterization of a Lighter-Than-Air Wind Energy System—Closing the Loop,” AIAA Paper No. 2015-3350.
Cobb, M. , Vermillion, C. , and Fathy, H. , 2016, “ Lab-Scale Experimental Crosswind Flight Control System Prototyping for an Airborne Wind Energy System,” ASME Paper No. DSCC2016-9737.
Deodhar, N. , Vermillion, C. , and Tkacik, P. , 2015, “ A Case Study in Experimentally-Infused Plant and Controller Optimization for Airborne Wind Energy Systems,” American Control Conference (ACC), Chicago, IL, July 1–3, pp. 2371–2376.
Deodhar, N. , and Vermillion, C. , 2016, “ A Framework for Fused Experimental/Numerical Plant and Control System Optimization Using Iterative G-Optimal Design of Experiments,” ASME Paper No. DETC2016-60488.
View article in PDF format.

## References

Fathy, H. K. , Papalambros, P. Y. , Ulsoy, A. G. , and Hrovat, D. , 2003, “ Nested Plant/Controller Optimization With Application to Combined Passive/Active Automotive Suspensions,” American Control Conference (ACC), Denver, CO, June 4–6, pp. 3375–3380.
Fathy, H. K. , Papalambros, P. Y. , and Ulsoy, A. G. , 2003, “ Integrated Plant, Observer, and Controller Optimization With Application to Combined Passive/Active Automotive Suspensions,” ASME Paper No. IMECE2003-42014.
Fathy, H. , Bortoff, S. , Copeland, S. , Papalambros, P. , and Ulsoy, A. , 2002, “ Nested Optimization of an Elevator and Its Gain-Scheduled LQG Controller,” ASME Paper No. IMECE2002-39273.
Peters, D. L. , Papalambros, P. Y. , and Ulsoy, A. G. , 2013, “ Sequential Co-Design of an Artifact and Its Controller Via Control Proxy Functions,” J. Mechatronics, 23(4), pp. 409–418.
Peters, D. L. , Papalambros, P. Y. , and Ulsoy, A. G. , 2011, “ Control Proxy Functions for Sequential Design and Control Optimization,” ASME J. Mech. Des., 133(9), p. 091007.
Alexander, M. J. , Allison, J. T. , and Papalambros, P. Y. , 2012, “ Decomposition-Based Design Optimization of Electric Vehicle Powertrains Using Proper Orthogonal Decomposition,” Int. J. Powertrains, 1(1), pp. 72–92.
Allison, J. T. , and Nazari, S. , 2010, “ Combined Plant and Controller Design Using Decomposition-Based Design Optimization and the Minimum Principle,” ASME Paper No. DETC2010-28887.
Allison, J. T. , Guo, T. , and Han, Z. , 2014, “ Co-Design of an Active Suspension Using Simultaneous Dynamic Optimization,” ASME J. Mech. Des., 136(8), p. 081003.
Pil, A. C. , and Asada, H. H. , 1996, “ Integrated Structure/Control Design of Mechatronic Systems Using a Recursive Experimental Optimization Method,” IEEE/ASME Trans. Mechatronics, 1(3), pp. 191–203.
Archer, C. L. , Monache, L. D. , and Rife, D. L. , 2014, “ Airborne Wind Energy: Optimal Locations and Variability,” Int. J. Renewable Energy, 64, pp. 180–186.
MAKANI, 2006, “ Energy Kites,” Makani Power, Alameda, CA, accessed June 15, 2017,
Altaeros Energies, 2010, “ The Next Transformation in Rural Infrastructure Efficiency,” Altaeros Energies, Somerville, MA, accessed June 15, 2017,
Ampyx Power, 2008, “ The Energy of Tomorrow,” Ampyx Power, The Hague, The Netherlands, accessed June 15, 2017,
Kitegen, 2007, “ Energy is in the Air,” Kitegen Research, Torino, Italy, accessed June 15, 2017,
Windlift, 2017, “ Airborne Wind Energy Systems,” Windlift, Raleigh, NC, accessed June 15, 2017,
Fagiano, L. , Milanese, M. , and Piga, D. , 2012, “ Optimization of Airborne Wind Energy Generators,” Int. J. Robust Nonlinear Control, 22(18), pp. 2055–2083.
Williams, P. , Lansdorp, B. , and Ockels, W. , 2008, “ Optimal Crosswind Towing and Power Generation With Tethered Kites,” J. Guid., Control, Dyn., 31(1), pp. 81–93.
Vermillion, C. , Grunnagle, T. , Lim, R. , and Kolmanovsky, I. , 2014, “ Model-Based Plant Design and Hierarchical Control of a Prototype Lighter-Than-Air Wind Energy System, With Experimental Flight Test Results,” IEEE Trans. Control Syst. Technol., 22(2), pp. 531–542.
Bafandeh, A. , and Vermillion, C. , 2016, “ Real-Time Altitude Optimization of Airborne Wind Energy Systems Using Lyapunov-Based Switched Extremum Seeking Control,” American Control Conference (ACC), Boston, MA, July 6–8, pp. 4990–4995.
Nikpoorparizi, P. , Deodhar, N. , and Vermillion, C. , 2016, “ Combined Plant and Controller Performance Analysis and Optimization for an Energy-Harvesting Tethered Wing,” American Control Conference (ACC), Boston, MA, July 6–8, pp. 4089–4094.
Vermillion, C. , Glass, B. , and Greenwood, S. , 2014, “ Evaluation of a Water Channel-Based Platform for Characterizing Aerostat Flight Dynamics: A Case Study on a Lighter-Than-Air Wind Energy System,” AIAA Paper No. 2014-2711.
Vermillion, C. , Glass, B. , and Szalai, B. , 2014, “ Development and Full-Scale Experimental Validation of a Rapid Prototyping Environment for Plant and Control Design of Airborne Wind Energy Systems,” ASME Paper No. DSCC2014-5907.
Deese, J. T. , Muyimbwa, T. , Deodhar, N. A. , Vermillion, C. R. , and Tkacik, P. , 2015, “ Lab-Scale Characterization of a Lighter-Than-Air Wind Energy System—Closing the Loop,” AIAA Paper No. 2015-3350.
Cobb, M. , Vermillion, C. , and Fathy, H. , 2016, “ Lab-Scale Experimental Crosswind Flight Control System Prototyping for an Airborne Wind Energy System,” ASME Paper No. DSCC2016-9737.
Deodhar, N. , Vermillion, C. , and Tkacik, P. , 2015, “ A Case Study in Experimentally-Infused Plant and Controller Optimization for Airborne Wind Energy Systems,” American Control Conference (ACC), Chicago, IL, July 1–3, pp. 2371–2376.
Deodhar, N. , and Vermillion, C. , 2016, “ A Framework for Fused Experimental/Numerical Plant and Control System Optimization Using Iterative G-Optimal Design of Experiments,” ASME Paper No. DETC2016-60488.

## Figures

Fig. 1

Flowchart demonstrating the sequence of events involved with an iteration of the experimentally infusedoptimization process. In short, the process involves the design of experiments, running the simulations/experiments and generating a response surface, learning from the experiments, and reducing the size of the design space for the next iteration.

Fig. 3

Block diagram of the flight controller used to track altitude and attitude set points [24]. The present research uses fixed roll and altitude setpoint values.

Fig. 4

Water channel-based test platform at UNC Charlotte

Fig. 5

Isometric and bottom views of the 1/100-scale BAT model used for testing. The ballast holes provide the ability to modify the center of mass location between iterations.

Fig. 2

Left full-scale prototype of an Altaeros BAT used during 2013 flight testing [12]. Right: ground-frame and body-frame coordinates, along with key variables used in deriving equations of motion for the BAT.

Fig. 7

The left and right halves of this figure show a comparison between the numerically optimal configuration and experimentally infused optimum, respectively. A pronounced improvement in pitch angle and heading angle error with experimental infusion is evident from the plots.

Fig. 9

The experimentally infused optimum results in a slightly higher zenith angle than the numerically optimized design. Thus, the improvement in flight quality seen in Fig. 8 comes at the (relatively minor) cost of an increase in ground footprint.

Fig. 6

The experiments and numerical simulation points generated by G-optimal design at every iteration are shown in the left half. The right half shows surface plots of combined numerical and experimental response surface (Ĵ) for every (pp, pc) combination in the shrunken design space.

Fig. 8

This figure shows a comparison between pitch angle and heading angle responses for the numerically optimized design (left) versus the design resulting from the experimentally infused optimization approach (right)

## Tables

Table 2 Table showing estimated optimal design variables (values of p*) at each iteration
Note: Iteration zero indicates optimum design variables obtained from pure numerical simulations.
Table 1 Table showing estimated optimal design variables (values of p*) at each iteration, with the zeroth iteration corresponding to the purely numerically optimized design

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections