A Robust Modification of a Predictive Adaptive-Optic Control Method for Aero-Optics
|
|
- Barbara Clarke
- 6 years ago
- Views:
Transcription
1 AIAA Aviation June 2016, Washington, D.C. 47th AIAA Plasmadynamics and Lasers Conference AIAA A Robust Modification of a Predictive Adaptive-Optic Control Method for Aero-Optics Robert Burns 1, Eric Jumper 2, Stanislav Gordeyev 3 University of Notre Dame, Notre Dame, IN, Abstract A modification of a previous predictive adaptive-optic controller is presented in this paper. Conventional adaptive-optic controllers suffer from bandwidth limitations caused by latency in their control loops. This latency severely limits their capabilities in aero-optic applications that cannot be overcome with conventional feedback techniques. Our method uses prior knowledge of flow behavior to predict future behavior, and thus overcome deadtime. We have modified our previous neural network controller to use a linearized predictor, which we demonstrate to be more accurate, more robust to noise and flow disturbances, and less computationally expensive. Our previous neural network method showed disturbance rejection in the range of 35-55% in simulation over our test conditions in the most optically-active regions, while the improved method shows disturbance rejection between 45-75% over the same range. Additionally, we demonstrate that the predictive control method is stable, even in the presence of latency uncertainty. Nomenclature α = viewing angle (rad) β = modified elevation angle (rad) or integrator gain (-) ε = vector of wavefront prediction residuals λ = laser wavelength (m) ρ = air density (kg/m 3 ) Σ = diagonal matrix of singular values Φ = matrix of POD modes (-) φ = eigenmode vector (-) A = linear prediction matrix Az = azimuth (rad) a = modal coefficients C = compensator transfer function (-) D = aperture diameter (m) d = single-aperture aero-optic disturbance (m) E = root-mean-square residual error El = elevation (rad) f = disturbance frequency (Hz) 1 Graduate Student, Department of Aerospace and Mechanical Engineering, Hessert Laboratory for Aerospace Research, Notre Dame, IN 46556, AIAA Student Member. 2 Professor, Department of Aerospace and Mechanical Engineering, Fitzpatrick Hall of Engineering, Notre Dame, IN 46556, AIAA Fellow. 3 Research Associate Professor, Department of Aerospace and Mechanical Engineering, Hessert Laboratory for Aerospace Research, Notre Dame, IN 46556, AIAA Associate Fellow. 1 Copyright 2016 by William Robert Burns, Eric Jumper, Stanislav Gordeyev. Published by the, Inc., with permission.
2 f = vector-valued nonlinear prediction function g = function describing flow evolution h = observer function of special flow properties I = on-target beam intensity (W/m 2 ) I 0 = unaberrated beam intensity (W/m 2 ) K GD = Gladstone Dale constant (m 3 /kg) M = Mach number (-) or embedding dimension (-) N d = number of timesteps of delay in the feedback loop (-) n = index of refraction (-) OPD = Optical Path Difference (m) OPL = Optical Path Length (m) S = Strehl ratio (-) St = Strouhal number (-) U = coefficient matrix V = freestream velocity (m/s) v = wavefront vector (m) W = wavefront vector (m) x = vector of coefficients (-) y = wavefront residual output (m) z = vector of flow properties I. Introduction ne significant challenge to the implementation of airborne laser systems is the aero-optics problem 1. Near-field Oturbulence over a turret aperture causes density fluctuations that are imprinted on the laser wavefronts through the relationship between index of refraction and density. These effects can greatly reduce a beam control system s performance in the far-field, especially when pointing through the wake region of a turret. Turret configurations are desirable in beam control systems for their large fields of regard; however, this comes at the expense of significant aero-optic aberrations, especially in the wake region. A schematic of the flow structures that impact aero-optics is shown 2 in Figure 1. Pressure and temperature fluctuations contribute to significant density fluctuations as the boundary layer thickens downstream. Flow separation over the top of the turret results in large separation regions that roll up into horn and secondary vortices. So-called necklace vortices form at the base of the turrets and can be encountered when pointing the laser near the surface of the aircraft. Shocks can form at the top of the turret at Mach numbers as low as All of these flow features cause density fluctuations that adversely impact optical performance. The aero-optics problem in the wake region of a turret (without compensation) is depicted schematically in Figure 2. The beam s unaberrated, planar wavefronts encounter turbulence near the turret s aperture. The turbulence-related density variations are effectively imprinted on the wavefront. As the beam propagates to the far field, these density fluctuations will tend to scatter the beam. The strength of these aberrations scales with the squared Mach number of the emitting aircraft, the freestream density, and the diameter of the turret. The index of refraction is related to the density of the medium through the Gladstone-Dale relationship 3, n( x, y, z) 1 KGD( ) ( x, y, z), where n is the index of refraction, λ is the wavelength, and K GD is a constant. This relationship can be used to compute the exact aberration along a path by first computing the optical path length Figure 1: Optically-aberrating flow over a turret 1. 2
3 z 1 from some point in the aperture, OPL( x, y, t) n( x, y, z, t) dz, and then subtracting the aperture-averaged mean to compute the optical path difference, 0 OPD( x, y, t) OPL( x, y, t) OPL( x, y, t) [1] Figure 2: Schematic of the aero-optics problem for a beam control system 1. The RMS optical path difference (OPD) is a convenient metric for estimating the performance of aero-optic systems since it enables the estimation of the ratio of the Strehl ratio: aberrated on-axis beam intensity to the diffraction-limited on-axis beam intensity. Strehl ratios can be estimated for each aberration source in a beam control system: for example, near-field aero-optics, far-field atmospherics, aero-mechanical jitter, optical heating, etc. These Strehl ratios can be multiplied together to obtain the total loss of intensity at the target. In this paper, we focus on the aero-optic aberrations only. For instantaneous aero-optic aberrations following a Gaussian spatial distribution, the Strehl ratio can be estimated through the Marèchal large aperture approximation 4,5,6, shown in Eq OPDrms ( t) It () S() t e. [2] I () t It is desirable to maximize the Strehl ratio by mitigating aero-optical aberrations. It is conventional to quantify reductions in aero-optical aberration using the following definition, db OPD rms, corrected 10log. [3] 10 OPDrms, uncorrected Since laser wavelength is application-dependent, it is useful to consider the impact of wavefront aberration mitigation in non-dimensional terms by rearranging the equation in terms of a corrected and uncorrected Strehl ratio, as discussed in our previous work 7. A modest -3dB improvement in aperture-averaged disturbance rejection can improve Strehl ratio from 0.2 to 0.67, for example. A large body of work has been performed to characterize aero-optical effects on the Airborne Aero-Optics Laboratory 8 (AAOL), which is shown in Figure 3, and its transonic successor AAOL-T 9. AAOL consists of a pointing and a tracking aircraft. The pointing aircraft emits a diverging 532 nm beam toward the tracking aircraft s turret, so that the beam width at the exit aperture is small relative to the flow structures. This means that the optical aberrations at the exit aperture are primarily tip/tilt and do not significantly impact the wavefronts. Aero-optical aberrations near the aperture of the tracking turret are imprinted on the laser wavefronts and measured by a Shack- Hartmann wavefront sensor in the tracking aircraft. 3
4 Figure 3: Picture of the Airborne Aero-Optics Laboratory (AAOL). The pointing aircraft is depicted on the left and the tracking aircraft on the right. 8 Normalized aero-optical aberrations tend to collapse onto a curve that is a function of two quantities: a lookback angle, α, and a modified elevation angle, β. These two quantities can be expressed 8 in terms of azimuth and elevation as shown in Figure 4, cos cos( )cos( ), and [4] 1 Az El 1 tan( El) tan. [5] sin( Az) OPD It has been demonstrated 8 that aero-optic aberrations can be normalized using the following equation, norm OPDrms, [6] 2 M / D 0 where M is the freestream Mach number, ρ/ρ 0 is the ratio of freestream density to sea-level density, and D is the turret diameter. Figure 5 shows OPD norm as a function of viewing angle and modified elevation angle from AAOL. Notably, the optical aberrations are highest in the separated flow regime, corresponding to viewing angles of greater than 110 degrees. Laser-based systems are often tolerant of small aberrations encountered at forward-looking angles, but not at aft-looking viewing angles. Aero-optical aberrations must be mitigated using some technique for many applications involving either communications or directed energy. There are two primary ways of achieving this: flow control and adaptive optics. Flow control devices are generally very Figure 4: Definition of α and β in terms of azimuth and elevation 7 application-dependent. Active flow control can be used to mitigate aero-optic distortions either by suppressing aberrations or relocating them away from an optical aperture. Some of these methods are open-loop (operating on pre-determined, feedback-free control), while others are closed-loop (sensing and reacting to disturbances to mitigate them). Active flow control has been employed in transonic flows to manipulate shock formation over a cylinder 10 using pulsed jets in open-loop configurations. This type of manipulation could be used to relocate optically-active flow features away from an optical aperture. High-frequency synthetic jets 11 have been used to reduce the optical aberrations present in the wake of a hemisphere. Later, this group 12 used a hybrid of passive and active flow control: by using a forward-protruding plate to decouple the turret wake from the necklace vortices, the effectiveness of the oscillating jets upstream of the optical aperture to control the turret wake was enhanced. Vukasinovic et al 10 have also tested control jets in the transonic regime to control separation both upstream and downstream of the shock to reduce the sharp velocity and density gradients present in the shear layer. This type of flow control was studied in more detail 13 using Particle Image Velocimetry (PIV) and Schlieren imagery. 4
5 Figure 5: Collapse of OPD rms as a function of viewing angle from AAOL flight tests. 8 A modification of proper orthogonal decomposition 14 has been developed to model closed-loop flow control. This method attempts to model both the controlled and natural states of the flow, as well as the transition between these states. The research group demonstrated their temporal POD (TPOD) approach using dielectric discharge barrier (DBD) plasma actuators to control the first TPOD mode of flow over a circular cylinder. In practice it would be necessary to control more than the first mode, but nonetheless there are several ideas from this work that could be incorporated into a model for the physical behavior of various turbulent flows. An explosion of research in the field of machine learning has led to an effort in fluids research 15 to improve flow control using these new methods. One method 16 uses optical sensors for feedback, slotted jets for control actuators, and a genetic algorithm to implement a control law to delay flow separation. Genetic algorithms are used in control systems to search for control actions that minimize a cost function for a given set of observations. The results of their experiment were an 80% reduction in the recirculation zone behind a backward-facing step. Genetic algorithms are often referred to as model-free control. A similar approach based on genetic programming 17 was used to control a mixing layer flow and compared to open-loop forcing. In addition to achieving better performance than the open-loop forcing, it is also adaptable and automatically tunable. Another means of aero-optic mitigation is to use adaptive optics. These systems imprint the conjugate of the currently-present aero-optic aberration onto a deformable mirror, and then send the pre-corrected laser beam out from the turret. When the corrected beam then encounters the aero-optic aberration, it emerges from the flow as nearly planar. These techniques are currently employed on ground-based telescopes to correct for atmospheric optical effects. However, these aberrations occur on the order of 1-10 Hz and thus the computational and system requirements are not extremely demanding for currently-available technology. However, sub-sonic aero-optical phenomena consisting of turbulent boundary layers, shear layers, attached vortices, etc. occur on the order of kilohertz on full-scale turrets at realistic Mach numbers. In particular, the dominant frequency content of separated flow on a turret in the wake region typically occurs near a Strouhal ratio of about 1: fd St, [7] V where f is the frequency of the disturbance in hertz, D is the turret diameter, and V is the freestream velocity. For the Airborne Aero-Optics Laboratory, the disturbance frequencies observed in the wake region were concentrated near 1 khz. Thus, a real-time adaptive optics system appropriate for AAOL must be able to perform corrections at a minimum of 10kHz, assuming 10 discrete corrections per disturbance cycle. This is a demanding computational requirement given that wavefronts must be represented by matrices that are at least 10x10, and sometimes as large as 30x30. A significant obstacle is cumulative latency in the control system. Time spent on digital calculations, physical sensor response, analog-to-digital conversion of sensor signals, electronic amplifiers, and electromechanical response of deformable mirrors will all contribute to a significant degree of delay 18 in a highspeed adaptive optic system. 5
6 This latency problem has previously been addressed using various prediction methods such as dynamic mode decomposition 19, receding-horizon adaptive control 20, and simple frozen-flow advection prediction 21. We recently developed a method 7 for overcoming this latency problem using a predictive neural network controller in conjunction with a conventional feedback loop as shown in Figure 6. The neural network controller uses prior knowledge of flow behavior to predict future behavior over a short temporal horizon equal to the known delay in the feedback loop. We simulated an adaptive-optic control loop operating at 25kHz with varying amounts of latency using aero-optic wavefront disturbance data from the AAOL flight tests. In the most opticallyactive fully separated flow regime, disturbance rejection between 35% and 55% was achieved in our simulations. As discussed previously, even 35% reduction in RMS wavefront aberration can have a large effect on a real system s performance. The presumption for this method was that the nonlinear character of the flow would require a nonlinear predictive element such as a neural network. In this work, we develop a modification for the neural network method that uses a linearized prediction model. We demonstrate that not only is this an acceptable assumption over the amounts of latency considered, but it is actually scalable to a larger number of predictable modes and leads to better prediction accuracy. Additionally, we investigate its favorable robustness to perturbed flow conditions and stability. For further details on dimensionality reduction using POD and the impact of latency on controller bandwidth, refer to our previous work 7. II. Predictive Control System Architecture A block diagram of the predictive controller is shown in Figure 7. The controller is denoted by C, the mirror plant model by G, the aero-optic disturbance by d, the residual wavefront error y, the prediction model P, the lagged wavefront sensor H 1, the lagged mirror sensor by H 2, and the unlagged mirror sensor by H 3. For consistency and to avoid excessive tedium, we assume that all of the latency in the control system occurs in wavefront sensing, reconstruction, and control law computation and is accounted for in the model H 3. The wavefront predictor itself is sensitive only to cumulative latency and its predictions are independent of the individual sources of latencies. Lagged and unlagged mirror sensors are required to reconstruct an estimate of the total wavefront disturbance from the sensed residual wavefront and the deformable mirror position. For the purposes of our analysis, we will assume that the conventional controller is a simple integrator and the deformable mirror has unity transfer function. In the absence of latency this is a trivial control problem, but the bandwidth limitations imposed by pure latency are severe enough to investigate in isolation of mirror mechanics, amplifier electrical dynamics, etc. (see our previous paper for more details 7 ). s H 1 may be decomposed as T t H H e, where ΔT is the cumulative latency in the feedback loop, δt is the 1 1 uncertainty in this latency, and H is the latency-free sensor. Practically, the lagged and unlagged deformable 1 s T mirror sensors may be decomposed as H2 H3e since H 2 is simply a digitally-delayed copy of H 3. An important design goal is to match the latency between the lagged deformable mirror sensor and the wavefront sensor. The latency in H 2 is considered to be a tunable parameter matched to the nominal value of the latency in the wavefront sensor, while the latency in the wavefront sensor, ΔT, is measured a priori with some uncertainty, δt. Due to simple measurement uncertainty as well as digital jitter, this uncertainty will always be non-zero; however, the δt = 0 case is considered first. 6 Figure 6: Conventional AO loop augmented by wavefront predictor 1.
7 Figure 7: Predictive controller block diagram. Consider the single-input-single-output (SISO) version of this system. This system gives the following general transfer function from the output to the reference signal, C sgs y s. [8] r 1 C s G s P s H s C s G s P s H s C s G s H s The rationale for matching the delays between H 2 and H 1 is now clear: in the absence of uncertainty, the complex dynamics of the predictor vanish and the closed loop transfer function becomes simply, C sg s y s, [9] r 1 C s G s H s 3 which can typically be stabilized by a P-I controller in an adaptive-optics control loop. Given the simplicity of our idealized component models, the stability of this loop is trivial to assess in the absence of latency uncertainty. T s Now, consider the multidimensional case. Define the transfer function matrix T(s), y1 y1 y1 r1 r2 rn y2 y2 y2 r1 r2 r N yn yn yn r1 r2 rn s s s s s s s s s, [10] which in the multidimensional case can be solved for as follows if the latency uncertainty is zero, s s s s s s 1 T C G I C G H 3. [11] As in the single-dimensional case, this system can easily be controlled (with varying performance) with classical methods. Recall that adaptive optics system components are in general open-loop stable. Similarly, the disturbance rejection function, 7
8 S s y y y d1 d2 d y y y d1 d2 d y y y d1 d2 d s s s s s s s s s N N N becomes (in the absence of delay uncertainty), N N N, [12] S s I C s G s P s H s C s G s H s I C s G s H s. [13] Hence, the prediction filter must be open-loop stable to ensure closed-loop stability. The interpretation for this necessity is essentially that the predictor acts as a feed-forward element. We will now detail the algorithm of the predictor itself, and then re-assess the stability problem using this framework. III. Prediction Method and Estimation We will first briefly summarize the key assumptions from our nonlinear neural network predictor 7 and then discuss the necessary modifications. We assume that we can approximate the evolution of an aero-optical wavefront v k by some discrete function f that depends on prior wavefronts v k-1 v k-m+1, with some amount of truncation error denoted by ε k+1. k1 k k1 km 1 k1 v f v, v,..., v ε [14] A wavefront predictor should approximate f such that the mean-square of the residual error is minimized, L E ε, k1 k 2 where L is the total number of wavefront observations in a dataset. This is has the structure of a non-linear autoregressive (NAR) problem. For practical applications, it is desirable to reduce the order of the prediction model. Practical applications of adaptive optics in aero-optics will have wavefront measurements on the order of 15x15 up to 30x30 subapertures, ranging from ~200 subapertures to ~700 subapertures for a circular beam inscribed in a rectangular sensor array. Direct wavefront prediction is computationally-intensive and prone to measurement noise. A simple linear prediction matrix of the form v k+1 =Av k alone would contain in the worst case between 10 5 and 10 6 elements, which is impractical for real-time applications. This problem becomes substantially worse if more than a single prior wavefront is needed for a prediction; which, in the general case is true due to the nonlinear nature of the flow. A number of methods can be used for dimensionality reduction. A natural choice of basis functions for many optical applications are the Zernike modes, since they describe physically meaningful optical aberrations such as focus/defocus, coma, astigmatism, etc. and are orthogonal on the unit disk. However, these modes are physically unrelated to the physics of the flow itself. A better modal decomposition for aero-optics is the Proper Orthogonal Decomposition 22, which can be computed via a snapshot method that produces the most rapidly-converging set of orthogonal modes possible. These modes, also known as Karhunen-Loève (K-L) or Principle Component Analysis (PCA) modes, can be calculated for multidimensional discrete-time data using the Singular Value Decomposition (SVD) or by solving the eigenvalue problem on the data matrix s autocorrelation matrix. In this case, let V be a matrix of measurement snapshots of wavefronts organized by column vectors of samples ordered by increasing time as shown, T 1 2 N V v v v. [15] 8 1
9 The data matrix, V, is decomposed using SVD, H V UΣΦ. [16] The spatial modes are extracted from the columns of Φ, 1 2 R Φ φ φ φ, [17] where R rank( V ). The temporal coefficients are then calculated from a projection of the spatial modes onto the original observations. Taking advantage of the fact that the pseudoinverse of an orthonormal matrix is its Hermitian transpose, this can be written as H x Φ V Φ V, [18] where x is a matrix of temporal coefficients. This is an important fact in terms of computational efficiency since it means that the projection of wavefronts onto POD modes is a simple matrix multiplication. The POD modes are ranked by the importance of their contribution to the overall energy of the system. Quite often, the POD modes converge quickly to give a good low-dimensional model. Additionally, in the case of naturally occurring fluid flow, low order modes typically exhibit smooth behavior. If it is assumed that the POD modes do not change, then each wavefront can be decomposed as follows, N v x ( k) Φ, [19] k n n n1 where Φ n are the time-invariant POD modes, x n (k) are the coefficients at timestep k, and N is the desired truncation dimension. Additionally, we assume that M snapshots are sufficient to approximate the next wavefront in the sequence: we refer to M as the embedding dimension 23. In this case, the NAR problem becomes a function of the modal coefficients, as shown, N k1 1 k k1 km 1 k1 v Φ g x, x,..., x ε, [20] N LxN Φ R are the truncated set of POD modes and : where 1 g R Thus, it is necessary to use some method to estimate the function g. NxM N R is the nonlinear prediction function. We recognize at this point that every continuous, differentiable nonlinear function can be linearized at any point along a system s trajectory. For the practical application of latency compensation, we assume that the flow dynamics are approximately linear over the prediction horizon. This assumption allows us to write g as a matrix function, xk,,..., x x g x x x A ε xkm1 k 1 k 1 k k 1 k M 1 k, [21] where A is referred to as the prediction matrix. The prediction matrix can be calculated by setting up the overconstrained linear problem using a set of training input and output matrices X and Y, x1 x2 xlm 2 3 LM1 M 1 M 2... L x x x Y x x x ε AX ε A ε xm xm 1 xl 1, [22] 9
10 and then solving for the A matrix by calculating the pseudoinverse of the right-hand data matrix or similar technique. The resulting matrix will be optimal in the least-squares sense for the POD coefficient prediction; that is, it is desired to solve the convex optimization problem, arg min A LM k1 x k 1 xk A x km1 The pseudoinverse method to solve this problem is,. [23] 1 1 H H H Y AX AUΣV A Y UΣV YU Σ V, [24] or equivalently using ordinary least squares, 1 A YX XX H H. [25] This method is known as Vector Auto-Regression 24 (VAR), and combined with the POD model reduction we refer to it as POD-VAR. It is noteworthy that this solution coincides with the maximum likelihood estimate if the error terms are normally distributed. Additionally, it generates a stable predictor if the modeled process is stationary with time-invariant covariance. Writing the prediction problem such the prediction matrix is square, xk 1 xk A1 A2 AM xk xk xk 1 I 0 0 xk 1 A aug x x 0 I 0 x k M 2 k M 1 km 1 a stability analysis for this autonomous process simply requires that, sup zc : zi A 0 1. [27] aug VAR estimators will always meet this condition for natural stationary processes, such as mean-removed, tip/tiltremoved aero-optic disturbances resulting from an unchanging flow condition. IV. Stability of the Closed-Loop System with Delay Uncertainty The stability problem of Section 2 becomes more difficult to address if there is a mismatch between the prediction horizon and the true latency in the system; that is, if t 0. This mismatch will always be present in a real system while digital systems can be very good at measuring this sort of latency and matching it closely, the effect of digital electronics jitter and some small uncertainty cannot be ignored. The remainder of this section will address the effect of this mismatch. A general framework for investigating the effect of delay uncertainty is described in Chapter 11 of Michiels and Niculescu 25. The key results of their analysis will be applied here. First, make the simplifying assumptions that s s T s t s s s e e H s I. In simple terms, the assumptions are that the G I and H1 H2 H 1, where deformable mirror responds instantly while the wavefront sensor responds with some pure time delay and some uncertainty in that delay. The resulting transfer function becomes, T 1 s C s s s s e e I C C P [28] st s t 1 [26] 10
11 The integrator compensator will have the form, reduces to, C s I, where β is the feedback gain. The transfer function s 1 s 1 s T s t T s I P s e e 1. [29] Hence, the characteristic equation that must be studied for stability is, s st s t det I1 P se e 1 0. [30] Now, choose a proper matrix fractional description 1 P s N s D s such that N P (s) and D P (s) are coprime polynomial matrices and D P (s) is Hurwitz (allowing roots on the imaginary axis). This fractional matrix description will exist if the system is both controllable and observable 26. The stability of this characteristic equation can be treated using the theory of delay differential equations 25. If the above characteristic equation is written in the following form, st s t s s Pse e det D / N 1 0, [31] P then its stability can be analyzed by studying the neutral equation 25, 1 x( t) N 1 0D0 x t T t N0D0 x t T 0, [32] where N 0 and D 0 are the column-reduced coefficient matrices of N P s and D s P s /, respectively. The details of this calculation are not important in this case for reasons that will be shortly addressed, but refer to Michiels and Niculescu 25 for additional information. This neutral equation is strongly stable if 25, 1 i1 1 i i sup r N D e N D e : 0, 2, i 1, 2 1. r The stability condition with respect to infinitesimal delay uncertainty may be equivalently stated as, 1 1 ND 0 0. [33] 2 However, because s s s P N D 1 P P / is a strictly proper transfer function, the spectral radius above is 0 (as discussed in Michiels and Niculescu 25 ) and practical stability is preserved for infinitesimal delay mismatches. In summary, the delay compensation system outlined in this section is practically stable for infinitesimal delay mismatch if: 1. P(s) is proper. 2. P(s) is open-loop stable. 3. P(s) is observable and controllable. These conditions will be validated later in this section. While Eq. 26 is convenient for practical estimation, it is helpful to put the predictor into state-space form for application and further analysis. One such state space model is, x A x B u k 1 P k P k yk CPxk DPu. [34] k P 11
12 In this model, the input vector u k is the most recent measurement of POD coefficients, the output vector y k is the POD coefficient prediction at timestep k, and x k are a set of stored POD coefficient vectors that are used in the N NM prediction. Decompose the prediction matrix, NN A A A where A R. The state space matrices themselves have dimension N N DP R and are chosen to be, I 0 0 A P 0 I 0 0 0NN 0NN 0 0 I 0 NN NN NN NN NN NN NN NN NN NN NN NN NN NN A R as, N ( M 1) N ( M 1) AP R, 1 M N ( M 1) N BP R,, [35] k N N ( M 1) CP R, and B P I 0 0 NN NN NN, [36] 2 3 M C A A A, [37] P and, D A. [38] P 1 This form of the predictor allows a convenient description of the transfer function matrix in discrete form, P(z), as follows: z 1 P() z C I A B D. [39] P P P P As discussed previously in this section, this system must be controllable and observable for stability in the presence of delay uncertainty. Otherwise, internally-cancelling modes will exist in P(z) that appear in both the poles and zeros, and hence the resulting fractional decomposition P(z) = N(z)D(z) -1 will not be coprime. Controllability is tested with the following condition 26, N( M1) 1 rank[ P P P P P] NM ( 1) A A B A B, [40] and observability is tested with, CP rank CA P P NM 1. [41] NM1 1 CA P P For this system, 0 k AP BP I 0 knn NN ( M k 1) NN for k > 2, [42] 12
13 so controllability is guaranteed analytically through the structure of the system matrices (in particular, M 2 rank[ AP APBP AP B P] NM ( 1) ). The observability matrix will also be full-rank via a similar argument if C P is full-rank. Since C P is estimated from a natural process, C p will in general be full-rank. Now it is desired to assess the stability of this discrete predictor as it interacts with a real adaptive optics system. The approach taken here will use conformal mapping to transform discrete space to continuous space in such a manner that conventional s-plane analysis tools can be applied. A better (but far more complex) approach would be to use the theory of sampled data systems to accurately model the interaction of the digital predictor with an analog system. Future research should investigate this. Ts The exact mapping between the z-plane and the s-plane is z e. The difficulty with this mapping is that conventional analysis tools cannot be used with non-polynomial expressions. There are several useful approximations to this mapping that could be chosen to map between the s-plane and the z-plane. One approximation that is closely related to this is the Padé approximation 27, Ts 1 e, [43] 1 2 n k1s k0s kns 2 n k1s k0s kns which is used to approximate time delay in continuous systems. This approximation is derived through a Taylor series expansion of the exponential. A similar approximation is the Tustin transformation (or bilinear transformation), defined as 28, Ts 1 Ts / 2 z e. [44] 1 Ts / 2 An important characteristic of this mapping is that the stability boundary (the imaginary axis in the s-plane and the unit circle in the z-plane) is mapped exactly. Minimum-phase systems in the s-plane are also mapped to minimum phase systems in the z-plane, however the converse is not guaranteed. In fact, time delay systems are in general non-minimum phase. Minimum phase is not a requirement for this analysis, but it is useful to understand the limitations of this approximation. Finally, some frequencies will be significantly distorted as a consequence of this transformation particularly those near the Nyquist frequency, 2/T. The dominant disturbance frequencies addressed in this research are between 800 and 2000 Hz, while the control loop samping frequency is 25 khz. Since the primary disturbances occur at an order of magnitude lower than the sampling frequency, high-frequency distortion will amplify noise but should not significantly impact overall predictor response. Using the bilinear transformation, P(z) can be approximated in continuous space as P(s). Because the stability boundary is preserved through the transformation and P(z) is stable, P(s) is also stable (but not minimum phase). Because P(s) is controllable and observable, there exists a proper and stable matrix fractional description P(s) = N P (s)d P (s) -1 such that N P (s) and D P (s) are co-prime and D P (s) is Hurwitz. Therefore, the stability conditions outlined in this section are valid for this predictor and the system is stable in both the nominal case and in the presence of infinitesimal delay uncertainty. V. Simulation Results Our simulation software is written in C++ using LAPACK for linear algebra routines. The wavefront disturbance data was taken from Notre Dame s Airborne Aero-Optics Laboratory (AAOL) flight test program. Each test condition consists of 15,000 total wavefront snapshots. The first 3,000 wavefront snapshots are used to calculate POD modes and train the coefficient prediction matrix. This value was selected as this is sufficient many snapshots for well-converged POD modes. Then, the next 12,000 wavefronts are used to evaluate controller performance independently. Thus, a separate predictive controller is trained for each flow condition. 13
14 Rejection (db) An important selection during the training phase of the controller is to determine a priori the optimal number of POD modes to include in the controller. The lowest-frequency (and thus, more predictable) content in the disturbances tends to occur in the low-order POD modes. As the prediction window increases, the higher frequency content becomes harder to predict and thus diminishing returns are observed. An optimization study was performed to determine the optimal number of POD modes to include for varying amounts of latency. The results are shown in Figure 8. There is not much practical benefit to including more than about 20 modes for if larger amounts of latency are anticipated in the controller. For small amounts of latency, improvements are observed up to 64 modes. Since there is no significant performance penalty for using a larger number of modes, we will select N=16 (to match our prior work) and N=64. A practical application may require fewer modes to reduce computational cost Num Modes (N) Figure 8: Ensembled mean disturbance rejection of all viewing angles for latencies between 1 to 5 timesteps (ascending order). We may now expand the results of the previous optimization study to further investigate how the controller behaves. The full aperture-averaged disturbance rejection results for N=16, N=64 and M=4 are shown in Figure 9. Similar to the results obtained in our previous work 7 (labeled as POD-ANN ), we observe better disturbance rejection in the highly-separated flow regime around a viewing angle of 120 deg. However, we now observe disturbance rejection near -6.5 db for one timestep of latency with a viewing angle beyond 120 deg., whereas our previous controller achieved roughly -4dB. For a large amount of latency, we observe approximately -3dB of disturbance rejection while the neural network controller achieved roughly -2dB Figure 9: Performance comparison of the linear controller with the prior neural network controller for several amounts of latency. 14
15 That a linear matrix predictor performs better than a nonlinear predictor in this highly-separated flow regime is consistent with expectations, since a separated shear flow is dominated by linear dynamics. We may also investigate the structure of the prediction matrices as shown in Figure 10 to gain a better understanding of the dynamic dependencies between POD modes. Naturally, the strongest dependence of a prediction (timestep k+1) is on the immediately preceding coefficient (timestep k) as indicated by the diagonals. The next-strongest dependence is on the diagonals of the previous coefficient (timestep k-1). The reader will observe that this is approximately a discrete solution of the continuous differential equation, ax x 0. Higher-order derivatives (up to the 3 rd -order) are also automatically approximated in the matrix as evidenced by the diagonal terms. The influence matrices also indicate that strong dynamic coupling is present in the fluid flow. This is an important feature of the predictor. Capturing the coupling between each mode improves prediction quality. This is part of the reason that simple single-input-singleoutput conventional PID controllers do not perform well in aero-optic applications. An additional observation is that the prediction matrix is primarily lower-triangular, meaning that higher-order modes tend to depend on lower-order modes. A possible modification to the predictive controller would be to force the prediction matrix to be purely lower-triangular, which would cut the computational cost of the prediction roughly by half. Figure 10: Influence matrices at various viewing angles. VI. Robustness and Sensitivity We may investigate the robustness of the controller by training the predictor on a data set and establishing baseline performance as per the previous section, and then testing the controller once again using a data set from a perturbed flight condition. In order to make this possible, we make the assumption that the flow direction is consistent with respect to the plane of the optical aperture. We use the spectral four-beam Malley probe technique 29 to estimate the mean flow direction and then rotate the wavefronts with bilinear interpolation such that the flow direction is in a consistent direction. This step is important because low-order POD modes (i.e., the most energetic and predictable modes) tend to convey information about the convective nature of coherent structures while higher-order modes contain information about their evolution. If the flow angle changes substantially, the modes would attempt to predict an incorrect direction for the convective portions of the disturbances. While the four-beam Malley technique is computationally expensive, the aperture-plane flow angle could also be estimated geometrically from azimuth and elevation. For initial robustness test, we trained the predictive controller on 3,000 wavefront snapshots using flight test data at α = deg, β = 72.4 deg, M = 0.51, and ρ =.78 kg/m3 (14kft altitude). The controller was then evaluated using the following 12,000 wavefront snapshots to establish baseline performance. We then used exactly the same controller on a different set of wavefronts from a different flight test. The robustness test was performed on a data set with α = deg, β = 65.3 deg, M = 0.50, and ρ =.86 kg/m3 (11.5kft altitude). This is a substantially perturbed 15
16 Disturbance Rejection (%) flow condition from the training condition; however, it is still in the same highly-separated flow regime. The results of the robustness test are shown in Figure 11 for varying amounts of feedback latency. There is essentially no difference between the baseline and robust performance for small amounts of latency. For larger amounts of latency, errors in the training condition build up and disturbance rejection quality drops by up to 7%. 90% 80% 70% 60% 50% 40% 30% 20% 10% 0% Number Timesteps of Feedback Delay Baseline Rejection Robust Rejection Figure 11: Effect of perturbing the flow condition on controller performance for varying amounts of latency. We may further investigate what is really happening by comparing the POD modes of each flow condition. The POD modes are reordered such that they are arranged by similarity rather than energetic contribution, as shown in Figure 12. While some differences are certainly present, most of the modes from the original predictive controller are still present in the perturbed flow condition. This observation gives insight into why the predictive controller tends to perform well. Model reduction via POD not only reduces computational cost but also increases robustness by reducing sensitivity to perturbed flow conditions. Figure 12: POD mode comparison using mode matching. (Mode 1 comparison at the top-left, Mode 4 comparison at the top-right, Mode 32 comparison at the bottom-right.) 16
17 We now evaluate controller sensitivity further for several combinations of baseline cases and perturbed flow cases. The method for testing this is to train the predictive controller on a subset of wavefront data in a baseline condition, evaluate the controller on the subsequent wavefronts from the baseline condition to establish baseline disturbance rejection, and then evaluate the controller again using the wavefront disturbances from a perturbed flow condition. The sensitivity to each parameter can be determined by calculating the change in controller performance as a function of the change in each parameter. Ideally, the sensitivity to each parameter could be measured by perturbing one parameter per test case, and calculating the sensitivity directly. Since this data is not available, an indirect method will be required. The fully separated flow regime will be examined here since it is of primary interest for real applications. It is assumed that the sensitivity is constant in this regime; i.e., that performance degradation will follow linearly with parameter changes. For convenience, define the function f as the mean aperture-averaged RMS disturbance rejection of the controller over L test wavefronts, 1 L yn, rms f L d. [45] n 1 n, rms Then, define s as a vector describing the change in viewing angle and modified elevation angle, s. [46] The directional derivative of f in the direction of s will define the sensitivity of the controller to a perturbation in that direction, s s f f. [47] s Recognizing that the change in RMS disturbance rejection from the baseline to the perturbed condition is approximately equal to the LHS times the length of the perturbation vector, the following equation can be written, u f s f perturbed f. [48] baseline The above quantity may be measured directly from experiment and simulation. Each test condition then gives a new set of equations, f f fperturbed fbaseline, [49] which in turn leads to the following system of equations for N test sensitivity tests, fperturbed f baseline f fperturbed f baseline [50] f f Ntest Ntest perturbed f baseline N test The full set of test conditions used for this sensitivity analysis is listed in Table 1. 17
18 Table 1: Sensitivity test conditions (flat window). Condition # α (deg) β (deg) ρ (kg/m3) M (-) The resulting estimated sensitivity parameters are shown in Table 2. An important result of this analysis is that the POD-VAR control algorithm is nearly 4 times more sensitive to changes in modified elevation than to changes in viewing angle. This is a consequence of the fact that POD modes do not significantly change with α since this direction is closely aligned with the flow direction. Changes in β can result in encountering different parts of flow structures such as horn vortices, necklace vortices, secondary vortices, etc. Table 2: Sensitivity of mean disturbance rejection to perturbations in viewing and modified elevation angles. f/ α f/ β 0.10 (%/deg) 0.38 (%/deg) This has practical importance for system design. If the POD-VAR predictors were to be trained and stored as a function of α and β rather than dynamically updated in real-time, then it would be necessary to have roughly 4 times finer resolution in the β lookup than in the α lookup. On the other hand, if the POD-VAR predictor is updated in real-time in an outer loop, then these values can establish the drift in controller performance as a function of the predictor update loop and turret slew rates. VII. Conclusions Latency in adaptive-optic control systems significantly limits controller performance in aero-optic applications due to the high-frequency nature of disturbances. We have presented a modification to our previous neural network controller that focuses on mitigating this limitation of adaptive-optic systems using flow prediction. The new linear POD-VAR controller improves disturbance rejection from 35%-55% in the case of the neural network controller to 45%-75% over the same range of test conditions in simulation. While a nonlinear predictor may in general be more accurate under ideal conditions, cumulative error in multistep prediction problems tends to build up more rapidly in a nonlinear predictor while error buildup is not as rapid in the linear model. We have demonstrated good robustness to perturbed flow conditions for the most optically-active case and discussed some of the physical reasons for this characteristic. We have also shown that our controller is stable in both the nominal condition and in the presence of delay uncertainty. Future research should include realistic mirror and amplifier models, and re-evaluate the performance of the controller under these conditions. The problem of varying amounts of latencies for each component in the feedback loop does not affect the performance of the predictor in the nominal condition, but the stability analysis of this paper should be expanded to address this more realistic case. References 1 Wang, M., Mani, A., and Gordeyev, S., Physics and Computation of Aero-Optics, Annual Review of Fluid Mechanics, vol. 44, Jan. 2012, pp Gordeyev, S., and Jumper, E., Fluid dynamics and aero-optics of turrets, Progress in Aerospace Sciences, vol. 46, Nov. 2010, pp Gladstone, J. H., and Dale, T. P., Researches on the Refraction, Dispersion, and Sensitiveness of Liquids, Philosophical Transactions of the Royal Society of London, vol. 153, 1863, pp Mahajan, V. N., Strehl ratio for primary aberrations: some analytical results for circular and annular pupils, JOSA, vol. 72, 1982, pp
19 5 Ross, T. S., Limitations and applicability of the Maréchal approximation, Applied optics, vol. 48, 2009, pp Porter, C., Gordeyev, S., and Jumper, E., Large-aperture approximation for not-so-large apertures, Optical Engineering, vol. 52, 2013, pp Burns, R., Jumper, E., and Gordeyev, S., A Latency-Tolerant Architecture for Airborne Adaptive Optic Systems, 53rd Aerospace Sciences Meeting, Jumper, E. J., Zenk, M. A., Gordeyev, S., Cavalieri, D., and Whiteley, M. R., Airborne aero-optics laboratory, Optical Engineering, vol. 52, 2013, pp Jumper, E. J., Gordeyev, S., Cavalieri, D., Rollins, P., Whiteley, M. R., and Krizo, M. J., Airborne Aero-Optics Laboratory-Transonic (AAOL-T), AIAA Paper, vol. 675, 2015, p Vukasinovic, B., Gissen, A., Glezer, A., and Gogineni, S., Fluidic Control of Transonic Shock-Induced Separation,, Vukasinovic, B., Glezer, A., Gordeyev, S., Jumper, E., and Kibens, V., Active control and optical diagnostics of the flow over a hemispherical turret, AIAA Paper, vol. 598, 2008, p Vukasinovic, B., Glezer, A., Gordeyev, S., Jumper, E., and Bower, W. W., Flow control for aero-optics application, Experiments in Fluids, vol. 54, Mar Gissen, A. N., Vukasinovic, B., Glezer, A., Gogineni, S., Paul, M. C., and Wittich, D. J., Active Transonic Shock Control,, Gordeyev, S., and Thomas, F. O., A temporal proper orthogonal decomposition (tpod) method for closed-loop flow control, AIAA Paper, vol. 359, Bunton, Steven, and Noack, Bernd, Closed-Loop Turbulence Control: Progress and Challenges, Applied Mechanics Reviews, vol. 67, Sep Gautier, N., Aider, J.-L., Duriez, T., Noack, B. R., Segond, M., and Abel, M., Closed-loop separation control using machine learning, Journal of Fluid Mechanics, vol. 770, May 2015, pp Parezanović, V., Duriez, T., Cordier, L., Noack, B. R., Delville, J., Bonnet, J.-P., Segond, M., Abel, M., and Brunton, S. L., Closed-loop control of an experimental mixing layer using machine learning control, arxiv preprint arxiv: , Whiteley, M. R., and Gibson, J. S., Adaptive laser compensation for aero optics and atmospheric disturbances, 38th AIAA Plasmadynamics and Lasers Conf, Goorskey, D. J., Schmidt, J., and Whiteley, M. R., Efficacy of predictive wavefront control for compensating aero-optical aberrations, Optical Engineering, vol. 52, 2013, pp Tesch, J., Gibson, S., and Verhaegen, M., Receding-horizon adaptive control of aero-optical wavefronts, Optical Engineering, vol. 52, 2013, pp Page, K. A., Applications of Linear Predictors in Adaptive Optics, The University of Wyoming, Berkooz, G., Holmes, P., and Lumley, J. L., The Proper Orthogonal Decomposition in the Analysis of Turbulent Flows, Annual Review of Fluid Mechanics, vol. 25, 1993, pp Kennel, M. B., Brown, R., and Abarbanel, H. D. I., Determining embedding dimension for phase-space reconstruction using a geometrical construction, Physical Review A, vol. 45, 1992, pp Lütkepohl, H., New Introduction to Multiple Time Series Analysis, Berlin: Springer, Michiels, W., and Niculescu, S., Stability, Control, and Computation for Time-Delay Systems: An Eigenvalue- Based Approach, Society for Industrial and Applied Mathematics, Antsaklis, P. J., and Michel, A. N., A Linear Systems Primer, Boston, MA: Birkhäuser,
Improvement in Optical Environment over Turrets with Flat Window Using Passive Flow Control
41th AIAA Plasmadynamics and Lasers Conference 8 Jun 1 Jul 010, Chicago, IL AIAA-010-449 Improvement in Optical Environment over Turrets with Flat Window Using Passive Flow Control Stanislav Gordeyev *,
More informationApproach for two-dimensional velocity mapping
Approach for two-dimensional velocity mapping Shaddy Abado Stanislav Gordeyev Eric Jumper Optical Engineering 52(7), 071402 (July 2013) Approach for two-dimensional velocity mapping Shaddy Abado Stanislav
More informationAero-Optical Mitigation of Shocks Around Turrets at Transonic Speeds Using Passive Flow Control
51st AIAA Aerospace Sciences Meeting 7-10 Jan 2013, Grapevine, TX 9-12 Jan 2012, Nashville, TN Aero-Optical Mitigation of Shocks Around Turrets at Transonic Speeds Using Passive Flow Control F Stanislav
More informationFluidic Control of a Turret Wake, Part II: Aero-Optical Effects
47 th AIAA Aerospace Sciences Meeting and Exhibit AIAA-009-0817 5-8 January 009, Orlando, FL Fluidic Control of a Turret Wake, Part II: Aero-Optical Effects Stanislav Gordeyev, Eric Jumper, Department
More informationPassive Mitigation of Aero-Induced Mechanical Jitter of Flat-Windowed Turrets
42th AIAA Plasmadynamics and Lasers Conference 27 30 Jun 2011, Honolulu, HI Passive Mitigation of Aero-Induced Mechanical Jitter of Flat-Windowed Turrets Zach B. Ponder 1, Stanislav Gordeyev 2, Eric J.
More informationAero-Optical Investigation of Shock-Related Effects on Hemisphere- On-Cylinder Turrets at Transonic Speeds
Aero-Optical Investigation of Shock-Related Effects on Hemisphere- On-Cylinder Turrets at Transonic Speeds Jacob Morrida, a Stanislav Gordeyev, b Nicholas De Lucca, a Eric Jumper c a University of Notre
More informationEffects of Engine Acoustic Waves on Optical Environment in. Flight on AAOL-T
Effects of Engine Acoustic Waves on Optical Environment in Flight on AAOL-T Nicholas De Lucca 1, Stanislav Gordeyev. 2 And Eric Jumper 3 University of Notre Dame, Notre Dame, IN, 46556 and Donald J. Wittich
More informationShack-Hartmann Wavefront Measurements of Supersonic Turbulent Boundary Layers in the TGF
45th AIAA Plasmadynamics and Lasers Conference 16-0 June 014, Atlanta, GA AIAA-014-493 Shack-Hartmann Wavefront Measurements of Supersonic Turbulent Boundary Layers in the TGF Adam E. Smith 1, Stanislav
More informationDirected Energy Beam Improvement Using Binary Control for the Advanced Tactical Laser (DEBI-BATL)
Directed Energy Beam Improvement Using Binary Control for the Advanced Tactical Laser (DEBI-BATL) William W. Bower, Valdis Kibens, The Boeing Company, MC S306-430, P.O. Box 516, Saint Louis, MO 63166 David
More informationFar-Field Laser Intensity Drop-Outs Caused by Turbulent Boundary Layers
Beam Control Conference 3-6 May 011, Orlando, FL Far-Field Laser Intensity Drop-Outs Caused by Turbulent Boundary Layers Stanislav Gordeyev, Jacob Cress and Eric Jumper, Department of Aerospace and Mechanical
More informationWavefront Measurements of a Supersonic Boundary Layer Using a Laser-Induced Breakdown Spark
AIAA Aviation 22-26 June 2015, Dallas, TX 46th AIAA Plasmadynamics and Lasers Conference AIAA 2015-2804 Wavefront Measurements of a Supersonic Boundary Layer Using a Laser-Induced Breakdown Spark Minh
More informationAero-Optical Measurements of High-Mach Supersonic Boundary Layers
46th AIAA Plasmadynamics and Lasers Conference 6 June 015, Dallas, TX Aero-Optical Measurements of High-Mach Supersonic Boundary Layers Stanislav Gordeyev 1, R. Mark Rennie, Department of Aerospace and
More informationSubsonic Boundary-Layer Wavefront Spectra for a Range of Reynolds Numbers
45th AIAA Plasmadynamics and Lasers Conference 16-0 June 014, Atlanta, GA AIAA-014-491 Subsonic Boundary-Layer Wavefront Spectra for a Range of Reynolds Numbers Adam E. Smith 1, Stanislav Gordeyev, University
More informationWavefront reconstruction for adaptive optics. Marcos van Dam and Richard Clare W.M. Keck Observatory
Wavefront reconstruction for adaptive optics Marcos van Dam and Richard Clare W.M. Keck Observatory Friendly people We borrowed slides from the following people: Lisa Poyneer Luc Gilles Curt Vogel Corinne
More informationAOL Spring Wavefront Sensing. Figure 1: Principle of operation of the Shack-Hartmann wavefront sensor
AOL Spring Wavefront Sensing The Shack Hartmann Wavefront Sensor system provides accurate, high-speed measurements of the wavefront shape and intensity distribution of beams by analyzing the location and
More informationStudies of the large-scale structure in adiabatic and moderately-wall-heated subsonic boundary layers
June 30 - July 3, 2015 Melbourne, Australia 9 7A-3 Studies of the large-scale structure in adiabatic and moderately-wall-heated subsonic boundary layers Stanislav Gordeyev Department of Aerospace and Mechanical
More informationInvestigation of Transonic Flow Behavior around a Three- Dimensional Turret Using Particle Image Velocimetry
Investigation of Transonic Flow Behavior around a Three- Dimensional Turret Using Particle Image Velocimetry Carlos Caballero College of Engineering, University of Florida Light distortions produced by
More informationControl of the Keck and CELT Telescopes. Douglas G. MacMartin Control & Dynamical Systems California Institute of Technology
Control of the Keck and CELT Telescopes Douglas G. MacMartin Control & Dynamical Systems California Institute of Technology Telescope Control Problems Light from star Primary mirror active control system
More informationSky demonstration of potential for ground layer adaptive optics correction
Sky demonstration of potential for ground layer adaptive optics correction Christoph J. Baranec, Michael Lloyd-Hart, Johanan L. Codona, N. Mark Milton Center for Astronomical Adaptive Optics, Steward Observatory,
More informationOptical Characterization of Nozzle-Wall Mach-6 Boundary Layers
54th AIAA Aerospace Sciences Meeting 4-8 Jan 016, San Diego, CA Optical Characterization of Nozzle-Wall Mach-6 Boundary Layers Stanislav Gordeyev 1 and Thomas J. Juliano, Department of Aerospace and Mechanical
More informationResponse of DIMM turbulence sensor
Response of DIMM turbulence sensor A. Tokovinin Version 1. December 20, 2006 [tdimm/doc/dimmsensor.tex] 1 Introduction Differential Image Motion Monitor (DIMM) is an instrument destined to measure optical
More informationActive Control and Optical Diagnostics of the Flow over a Hemispherical Turret
46th AIAA Aerospace Sciences Meeting and Exhibit 7-1 January 28, Reno, Nevada AIAA 28-598 Active Control and Optical Diagnostics of the Flow over a Hemispherical Turret Bojan Vukasinovic, Ari Glezer, Woodruff
More informationNumerical atmospheric turbulence models and LQG control for adaptive optics system
Numerical atmospheric turbulence models and LQG control for adaptive optics system Jean-Pierre FOLCHER, Marcel CARBILLET UMR6525 H. Fizeau, Université de Nice Sophia-Antipolis/CNRS/Observatoire de la Côte
More informationPolarization Shearing Interferometer (PSI) Based Wavefront Sensor for Adaptive Optics Application. A.K.Saxena and J.P.Lancelot
Polarization Shearing Interferometer (PSI) Based Wavefront Sensor for Adaptive Optics Application A.K.Saxena and J.P.Lancelot Adaptive Optics A Closed loop Optical system to compensate atmospheric turbulence
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More informationUsing a Membrane DM to Generate Zernike Modes
Using a Membrane DM to Generate Zernike Modes Author: Justin D. Mansell, Ph.D. Active Optical Systems, LLC Revision: 12/23/08 Membrane DMs have been used quite extensively to impose a known phase onto
More informationAerodynamic Design of an Aircraft-Mounted Pod for Improved Aero-Optic Performance
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 010, Orlando, Florida AIAA 010-437 Aerodynamic Design of an Aircraft-Mounted Pod for Improved
More informationWavefront Sensing using Polarization Shearing Interferometer. A report on the work done for my Ph.D. J.P.Lancelot
Wavefront Sensing using Polarization Shearing Interferometer A report on the work done for my Ph.D J.P.Lancelot CONTENTS 1. Introduction 2. Imaging Through Atmospheric turbulence 2.1 The statistics of
More information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationA NEW METHOD FOR ADAPTIVE OPTICS POINT SPREAD FUNCTION RECONSTRUCTION
Florence, Italy. May 2013 ISBN: 978-88-908876-0-4 DOI: 10.12839/AO4ELT3.13328 A NEW METHOD FOR ADAPTIVE OPTICS POINT SPREAD FUNCTION RECONSTRUCTION J. Exposito 1a, D. Gratadour 1, G. Rousset 1, Y. Clénet
More informationHypersonic Flight Effects on Optical Sensors
A Tutorial Of: Hypersonic Flight Effects on Optical Sensors Matt Salem The University of Arizona: OPTI 521 12/4/2016 Background: In recent years hypersonic vehicles have received a lot of attention from
More informationSPATIO-TEMPORAL PREDICTION FOR ADAPTIVE OPTICS WAVEFRONT RECONSTRUCTORS
SPATIO-TEMPORAL PREDICTION FOR ADAPTIVE OPTICS WAVEFRONT RECONSTRUCTORS Michael Lloyd-Hart and Patrick McGuire Center for Astronomical Adaptive Optics, Steward Observatory, University of Arizona, Tucson,
More informationUsing Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver
Using Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708 0300 A novel nonlinear
More informationPredictions of Aero-Optical Distortions Using LES with Wall Modeling
Predictions of Aero-Optical Distortions Using LES with Wall Modeling Mohammed Kamel, Kan Wang and Meng Wang University of Notre Dame, Notre Dame, IN 46556 Large-eddy simulation (LES) with wall-modeling
More informationPreliminary Study of Single Particle Lidar for Wing Wake Survey. M.Valla, B. Augère, D. Bailly, A. Dolfi-Bouteyre, E. Garnier, M.
Preliminary Study of Single Particle Lidar for Wing Wake Survey M.Valla, B. Augère, D. Bailly, A. Dolfi-Bouteyre, E. Garnier, M. Méheut Context of research Clean Sky Joint Technology Initiative : Aims
More informationExperimental investigation of flow control devices for the reduction of transonic buffeting on rocket afterbodies
Experimental investigation of flow control devices for the reduction of transonic buffeting on rocket afterbodies F.F.J. Schrijer 1, A. Sciacchitano 1, F. Scarano 1 1: Faculty of Aerospace Engineering,
More informationAPPLICATION OF ADAPTIVE CONTROLLER TO WATER HYDRAULIC SERVO CYLINDER
APPLICAION OF ADAPIVE CONROLLER O WAER HYDRAULIC SERVO CYLINDER Hidekazu AKAHASHI*, Kazuhisa IO** and Shigeru IKEO** * Division of Science and echnology, Graduate school of SOPHIA University 7- Kioicho,
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationEXPERIMENTS OF CLOSED-LOOP FLOW CONTROL FOR LAMINAR BOUNDARY LAYERS
Fourth International Symposium on Physics of Fluids (ISPF4) International Journal of Modern Physics: Conference Series Vol. 19 (212) 242 249 World Scientific Publishing Company DOI: 1.1142/S211945128811
More informationThe IPIE Adaptive Optical System Application For LEO Observations
The IPIE Adaptive Optical System Application For LEO Observations Eu. Grishin (1), V. Shargorodsky (1), P. Inshin (2), V. Vygon (1) and M. Sadovnikov (1) (1) Open Joint Stock Company Research-and-Production
More informationA Comparative Study on Automatic Flight Control for small UAV
Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic
More informationIMPROVEMENTS IN ACTIVE NOISE CONTROL OF HELICOPTER NOISE IN A MOCK CABIN ABSTRACT
IMPROVEMENTS IN ACTIVE NOISE CONTROL OF HELICOPTER NOISE IN A MOCK CABIN Jared K. Thomas Brigham Young University Department of Mechanical Engineering ABSTRACT The application of active noise control (ANC)
More informationAnalysis of Communication Systems Using Iterative Methods Based on Banach s Contraction Principle
Analysis of Communication Systems Using Iterative Methods Based on Banach s Contraction Principle H. Azari Soufiani, M. J. Saberian, M. A. Akhaee, R. Nasiri Mahallati, F. Marvasti Multimedia Signal, Sound
More informationSky Projected Shack-Hartmann Laser Guide Star
Sky Projected Shack-Hartmann Laser Guide Star T. Butterley a, D.F. Buscher b, G. D. Love a, T.J. Morris a, R. M. Myers a and R. W. Wilson a a University of Durham, Dept. of Physics, Rochester Building,
More informationSimultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors
Simultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors Sean Borman and Robert L. Stevenson Department of Electrical Engineering, University of Notre Dame Notre Dame,
More informationApplications of the Abbe Sine Condition in Multi-Channel Imaging Systems
Applications of the Abbe Sine Condition in Multi-Channel Imaging Systems Barbara Kruse, University of Arizona, College of Optical Sciences May 6, 2016 Abstract Background In multi-channel imaging systems,
More informationWavefront errors due to atmospheric turbulence Claire Max
Wavefront errors due to atmospheric turbulence Claire Max Page 1 Kolmogorov turbulence, cartoon solar Outer scale L 0 Inner scale l 0 h Wind shear convection h ground Page Atmospheric Turbulence generally
More informationActive optics challenges of a thirty meter segmented mirror telescope
Active optics challenges of a thirty meter segmented mirror telescope George Z. Angeli 1, Robert Upton 1, Anna Segurson 1, Brent Ellerbroek 1 1 New Initiatives Office, AURA Inc. ABSTRACT Ground-based telescopes
More informationError Budgets, and Introduction to Class Projects. Lecture 6, ASTR 289
Error Budgets, and Introduction to Class Projects Lecture 6, ASTR 89 Claire Max UC Santa Cruz January 8, 016 Page 1 What is residual wavefront error? Telescope AO System Science Instrument Very distorted
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationFURTHER RESULTS ON ROBUST CONTROL OF MICROVIBRATIONS ON MASS LOADED PANELS
FURTHER RESULTS ON ROBUST CONTROL OF MICROVIBRATIONS ON MASS LOADED PANELS G. S. Aglietti, J. Stoustrup, E. Rogers, R. S. Langley, S. B. Gabriel, Depts. of Aero & Astro/Electronics and Computer Science,
More informationMAE 142 Homework #5 Due Friday, March 13, 2009
MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.
More informationOptical Interface for MSC.Nastran
Optical Interface for MSC.Nastran Victor Genberg, Keith Doyle, Gregory Michels Sigmadyne, Inc., 803 West Ave, Rochester, NY 14611 genberg@sigmadyne.com Abstract Thermal and structural output from MSC.Nastran
More informationShallow Water Fluctuations and Communications
Shallow Water Fluctuations and Communications H.C. Song Marine Physical Laboratory Scripps Institution of oceanography La Jolla, CA 92093-0238 phone: (858) 534-0954 fax: (858) 534-7641 email: hcsong@mpl.ucsd.edu
More informationAFRL-SR-AR-TR b. GRANT NUMBER F c. PROGRAM ELEMENT NUMBER
REPORT DOCUMENTATION PAGE Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instrut maintaining the data needed, and
More informationOpto-Mechanical I/F for ANSYS
Opto-Mechanical I/F for ANSYS Victor Genberg, Gregory Michels, Keith Doyle Sigmadyne, Inc. Abstract Thermal and structural output from ANSYS is not in a form useful for optical analysis software. Temperatures,
More informationAdaptive Optics for the Giant Magellan Telescope. Marcos van Dam Flat Wavefronts, Christchurch, New Zealand
Adaptive Optics for the Giant Magellan Telescope Marcos van Dam Flat Wavefronts, Christchurch, New Zealand How big is your telescope? 15-cm refractor at Townsend Observatory. Talk outline Introduction
More informationWavefront Correction of Model-based Sensorless Adaptive Optics System
Wavefront Correction of Model-based Sensorless Adaptive Optics System Huizhen Yang 1*, Jian Wu 2 1. School of Electronic Engineering, Huaihai Institute of Technology, Lianyungang, China 222005; 2. School
More informationWhite Paper FINAL REPORT AN EVALUATION OF THE HYDRODYNAMICS MECHANISMS WHICH DRIVE THE PERFORMANCE OF THE WESTFALL STATIC MIXER.
White Paper FINAL REPORT AN EVALUATION OF THE HYDRODYNAMICS MECHANISMS WHICH DRIVE THE PERFORMANCE OF THE WESTFALL STATIC MIXER Prepared by: Dr. Thomas J. Gieseke NUWCDIVNPT - Code 8233 March 29, 1999
More informationStat 159/259: Linear Algebra Notes
Stat 159/259: Linear Algebra Notes Jarrod Millman November 16, 2015 Abstract These notes assume you ve taken a semester of undergraduate linear algebra. In particular, I assume you are familiar with the
More informationDeformable mirror fitting error by correcting the segmented wavefronts
1st AO4ELT conference, 06008 (2010) DOI:10.1051/ao4elt/201006008 Owned by the authors, published by EDP Sciences, 2010 Deformable mirror fitting error by correcting the segmented wavefronts Natalia Yaitskova
More informationExpected Performance From WIYN Tip-Tilt Imaging
Expected Performance From WIYN Tip-Tilt Imaging C. F. Claver 3 September 1997 Overview Image motion studies done at WIYN show that a significant improvement to delivered image quality can be obtained from
More informationOrthonormal vector polynomials in a unit circle, Part I: basis set derived from gradients of Zernike polynomials
Orthonormal vector polynomials in a unit circle, Part I: basis set derived from gradients of ernike polynomials Chunyu hao and James H. Burge College of Optical ciences, the University of Arizona 630 E.
More informationLecture 12. AO Control Theory
Lecture 12 AO Control Theory Claire Max with many thanks to Don Gavel and Don Wiberg UC Santa Cruz February 18, 2016 Page 1 What are control systems? Control is the process of making a system variable
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationInternational Conference on Methods of Aerophysical Research, ICMAR 2008
International Conference on Methods of Aerophysical Research, ICMAR 8 EXPERIMENTAL STUDY OF UNSTEADY EFFECTS IN SHOCK WAVE / TURBULENT BOUNDARY LAYER INTERACTION P.A. Polivanov, А.А. Sidorenko, A.A. Maslov
More informationShallow Water Fluctuations and Communications
Shallow Water Fluctuations and Communications H.C. Song Marine Physical Laboratory Scripps Institution of oceanography La Jolla, CA 92093-0238 phone: (858) 534-0954 fax: (858) 534-7641 email: hcsong@mpl.ucsd.edu
More informationIntegrating MD Nastran with Optical Performance Analysis
Integrating MD Nastran with Optical Performance Analysis Victor Genberg, Gregory Michels Sigmadyne, Inc., 803 West Ave, Rochester, NY 14611 genberg@sigmadyne.com Abstract The development of products in
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationWave-front generation of Zernike polynomial modes with a micromachined membrane deformable mirror
Wave-front generation of Zernike polynomial modes with a micromachined membrane deformable mirror Lijun Zhu, Pang-Chen Sun, Dirk-Uwe Bartsch, William R. Freeman, and Yeshaiahu Fainman We investigate the
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationComparison of the Performance of Modal Control Schemes for an Adaptive Optics System and Analysis of the Effect of Actuator Limitations
Comparison of the Performance of Modal Control Schemes for an Adaptive Optics System and Analysis of the Effect of Actuator Limitations Masaki Nagashima, Jae Jun Kim, Brij Agrawal Abstract In this study,
More informationApplied Mathematics 205. Unit II: Numerical Linear Algebra. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit II: Numerical Linear Algebra Lecturer: Dr. David Knezevic Unit II: Numerical Linear Algebra Chapter II.3: QR Factorization, SVD 2 / 66 QR Factorization 3 / 66 QR Factorization
More informationOverview of sparse system identification
Overview of sparse system identification J.-Ch. Loiseau 1 & Others 2, 3 1 Laboratoire DynFluid, Arts et Métiers ParisTech, France 2 LIMSI, Université d Orsay CNRS, France 3 University of Washington, Seattle,
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More informationMachine Learning. B. Unsupervised Learning B.2 Dimensionality Reduction. Lars Schmidt-Thieme, Nicolas Schilling
Machine Learning B. Unsupervised Learning B.2 Dimensionality Reduction Lars Schmidt-Thieme, Nicolas Schilling Information Systems and Machine Learning Lab (ISMLL) Institute for Computer Science University
More informationCHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER
114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationB553 Lecture 5: Matrix Algebra Review
B553 Lecture 5: Matrix Algebra Review Kris Hauser January 19, 2012 We have seen in prior lectures how vectors represent points in R n and gradients of functions. Matrices represent linear transformations
More informationUNIVERSITY OF CALGARY. Base Region Topology of Turbulent Wake around Finite Wall-Mounted Cylinder with. Application of Low Order Flow Representation
UNIVERSITY OF CALGARY Base Region Topology of Turbulent Wake around Finite Wall-Mounted Cylinder with Application of Low Order Flow Representation by Golriz Boorboor A THESIS SUBMITTED TO THE FACULTY OF
More informationControl-Structure Interaction in Extremely Large Telescopes
Control-Structure Interaction in Extremely Large Telescopes A. Preumont, B. Mokrani & R. Bastaits Université Libre de Bruxelles (ULB)-Active Structures Laboratory, Brussels, Belgium Abstract: The next
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More information1. INTRODUCTION ABSTRACT
Simulations of E-ELT telescope effects on AO system performance Miska Le Louarn* a, Pierre-Yves Madec a, Enrico Marchetti a, Henri Bonnet a, Michael Esselborn a a ESO, Karl Schwarzschild strasse 2, 85748,
More informationChapter 9 Robust Stability in SISO Systems 9. Introduction There are many reasons to use feedback control. As we have seen earlier, with the help of a
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Robust
More informationJournal of Fluid Science and Technology
Bulletin of the JSME Vol.9, No.3, 24 Journal of Fluid Science and Technology Re-evaluating wake width in turbulent shear flow behind an axisymmetric cylinder by means of higher order turbulence statistics
More information7 Principal Component Analysis
7 Principal Component Analysis This topic will build a series of techniques to deal with high-dimensional data. Unlike regression problems, our goal is not to predict a value (the y-coordinate), it is
More informationLECTURE NOTE #11 PROF. ALAN YUILLE
LECTURE NOTE #11 PROF. ALAN YUILLE 1. NonLinear Dimension Reduction Spectral Methods. The basic idea is to assume that the data lies on a manifold/surface in D-dimensional space, see figure (1) Perform
More informationMaking FEA Results Useful in Optical Analysis Victor Genberg, Gregory Michels Sigmadyne, Inc. Rochester, NY
Making FEA Results Useful in Optical Analysis Victor Genberg, Gregory Michels Sigmadyne, Inc. Rochester, NY Keith Doyle Optical Research Associates,Westborough, MA ABSTRACT Thermal and structural output
More informationCHARA Meeting 2017 Pasadena, California
MORE AUTOMATION Laszlo Sturmann M7 ACTUATORS LAB. LASER ALIGNMENT TELESCOPE OPTICAL ALIGNMENT NEW ACTUATORS REMOTELY ACTUATED M7 MOUNT MOTIVATION THE PRECISION OF THE COUDE ALIGNMENT WAS NOT SUFFICIENT
More informationPupil matching of Zernike aberrations
Pupil matching of Zernike aberrations C. E. Leroux, A. Tzschachmann, and J. C. Dainty Applied Optics Group, School of Physics, National University of Ireland, Galway charleleroux@yahoo.fr Abstract: The
More informationAstronomy 203 practice final examination
Astronomy 203 practice final examination Fall 1999 If this were a real, in-class examination, you would be reminded here of the exam rules, which are as follows: You may consult only one page of formulas
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationNA LASER GUIDE STAR AO WITH DYNAMICAL REFOCUS
Florence, Italy. Adaptive May 2013 Optics for Extremely Large Telescopes III ISBN: 978-88-908876-0-4 DOI: 10.12839/AO4ELT3.13893 NA LASER GUIDE STAR AO WITH DYNAMICAL REFOCUS Sebastian Rabien 1,a, Fernando
More informationRobust Internal Model Control for Impulse Elimination of Singular Systems
International Journal of Control Science and Engineering ; (): -7 DOI:.59/j.control.. Robust Internal Model Control for Impulse Elimination of Singular Systems M. M. Share Pasandand *, H. D. Taghirad Department
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationChapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 7 Interconnected
More informationCDS 101/110a: Lecture 8-1 Frequency Domain Design
CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve
More information