SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The development of a digital adaptive flight control algorithm is described that combines a weighted least squares estimator with optimal linear regulator control logic designed to minimize a quadratic model following performance index. Very promising results for a typical fighter aircraft's trajectory over six distinct flight conditions are reported.
This paper considers the problem of designing certain structurally constrained optimal regulators for linear systems subjected to additive white process noise and measurement noise. Three types of controller structures are considered, using direct output feedback, prespecified time constant filters, and optimal dynamic compensators. Necessary conditions are obtained for minimizing quadratic performance criteria. The techniques are demonstrated by application to a helicopter/slung load system, and a flexible space station.
A preliminary investigation into the use of modern control theory for the design of controls for a supersonic inlet is described. In particular, the task of controlling a mixed-compression supersonic inlet is formulated as a linear optimal stochastic control and estimation problem. An inlet can exhibit an undesirable instability due to excessive inlet normal shock motion. For the optimal control formulation of the inlet problem, a non quadratic performance index, which is equal to the expected frequency of inlet unstarts, is used. This physically meaningful performance index is minimized for a range of inlet disturbance and measurement noise covariances.
The problem of closed-loop control of maneuvers between two states for linear dynamical systems, subject to an arbitrarily specified terminal state, is considered. The feedback controller design is based on finite-time quadratic regulator theory. Closed-form expressions for the optimal control law are developed. Solutions are presented for both conventional and smoothed control profiles with fixed and/or free end condition problems. In the maneuvers using control-rate penalties, smooth profiles are generated throughout the maneuvers, in the sense that the initial condition jump discontinuities have been eliminated. Several examples involving large-angle maneuvers of a spacecraft are demonstrated. Results include control maneuvers from one state to another such as rest to rest and spin to rest, which effectively justify the solutions developed in this paper.
This paper presents a direct digital control law synthesis procedure for a large order, sampled data, linear feedback system using constrained optimization techniques to meet multiple design requirements. A linear quadratic Gaussian type cost function is minimized while satisfying a set of constraints on the design loads and responses. General expressions for gradients of the cost function and constraints, with respect to the digital control law design variables are derived analytically and computed by solving a set of discrete Liapunov equations. The designer can choose the structure of the control law and the design variables, hence a stable classical control law as well as an estimator-based full or reduced order control law can be used as an initial starting point. Selected design responses can be treated as constraints instead of lumping them into the cost function. This feature can be used to modify a control law, to meet individual root mean square response limitations as well as minimum single value restrictions. Low order, robust digital control laws were synthesized for gust load alleviation of a flexible remotely piloted drone aircraft.
Methods for multirate digital control system design are discussed. A simple method for sampling rate selection based on control bandwidths is proposed. Methods for generating a discrete-time state model of a sampled-data plant and a discrete-time equivalent to an analog cost function for a sampled-data plant are described. The succesive loop closures and linear quadratic Gaussian synthesis methods are reviewed, and a constrained optimization synthesis method is introduced. The proposed sampling rate selection, discretization, and synthesis methods are applied to two example design problems. Multirate and single-rate compensators synthesized by the different methods are compared, based on closed-loop responses, with compensators having the same real-time computation load.
The problem of controller design for flexible spacecraft is addressed. Model-based compensators, which rely on the knowledge of the system parameters to tune the state estimator, are considered. The instability mechanisms resulting from high sensitivity to parameter uncertainties are investigated. Dissipative controllers, which use collocated actuators and sensors, are also considered, and the robustness properties of constant-gain dissipative controllers in the presence of unmodeled elastic-mode dynamics, sensor/actuator nonlinearities, and actuator dynamics are summarized. In order to improve the performance without sacrificing robustness, a class of dissipative dynamic compensators is proposed and is shown to retain robust stability in the presence of second-order actuator dynamics if acceleration feedback is employed. A class of dissipative dynamic controllers is proposed which consists of a low-authority, constant-gain controller and a high-authority dynamic compensator. A procedure for designing an optimal dissipative dynamic compensator is given which minimizes a quadratic performance criterion. Such compensators offer the promise of better performance while still retaining robust stability.
Purpose: The current practice of considering the relative biological effectiveness (RBE) of protons in intensity modulated proton therapy (IMPT) planning is to use a generic RBE value of 1.1. However, RBE is indeed a variable depending on the dose per fraction, the linear energy transfer, tissue parameters, etc. In this study, we investigate the impact of using variable RBE based optimization (vRBE-OPT) on IMPT dose distributions compared by conventional fixed RBE based optimization (fRBE-OPT). Methods: Proton plans of three head and neck cancer patients were included for our study. In order to calculate variable RBE, tissue specific parameters were obtained from the literature and dose averaged LET values were calculated by Monte Carlo simulations. Biological effects were calculated using the linear quadratic model and they were utilized in the variable RBE based optimization. We used a Polak-Ribiere conjugate gradient algorithm to solve the model. In fixed RBE based optimization, we used conventional physical dose optimization to optimize doses weighted by 1.1. IMPT plans for each patient were optimized by both methods (vRBE-OPT and fRBE-OPT). Both variable and fixed RBE weighted dose distributions were calculated for both methods and compared by dosimetric measures. Results: The variable RBE weighted dose distributions were more homogenous within the targets, compared with the fixed RBE weighted dose distributions for the plans created by vRBE-OPT. We observed that there were noticeable deviations between variable and fixed RBE weighted dose distributions if the plan were optimized by fRBE-OPT. For organs at risk sparing, dose distributions from both methods were comparable. Conclusion: Biological dose based optimization rather than conventional physical dose based optimization in IMPT planning may bring benefit in improved tumor control when evaluating biologically equivalent dose, without sacrificing OAR sparing, for head and neck cancer patients. The research is supported in part by National Institutes of Health Grant No. 2U19CA021239-35.
The solutions to the optimal l2 to l2 disturbance rejection problem H(infinity) as well as to the LQG (linear quadratic Gaussian) H2 problem in periodic systems using the lifting technique are presented. Both problems involve a causality condition on the optimal LTI (linear time invariant) compensator when viewed in the lifted domain. The H(infinity) problem is solved using the Nehari's theorem, whereas in the H2 problem the solution is obtained using the projection theorem in Hilbert spaces. The authors demonstrate that exactly the same method of solution to the H(infinity) and H2 problems in periodic systems applies when considering the same problems in multirate sampled data systems.
The "direct control method" is a novel concept that is an attractive alternative and competitor to the differential-equation-based methods. The direct method is equally well applicable to nonlinear, linear, time-varying, and time-invariant systems. For all such systems, the method yields explicit closed-form control laws based on minimization of a quadratic control performance measure. We present an application of the direct method to the dynamics and optimal control of the Duffing system where the control performance measure is not restricted to a quadratic form and hence may include a quartic energy term. The results we present in this report also constitute further generalizations of our earlier work in "direct optimal control methodology." The approach is demonstrated for the optimal control of the Duffing equation with a softening nonlinear stiffness.
With increasing penetration of distributed generation in the distribution networks (DN), the secure and optimal operation of DN has become an important concern. In this paper, an iterative quadratic constrained quadratic programming model to minimize voltage deviations and maximize distributed energy resource (DER) active power output in a three phase unbalanced distribution system is developed. The optimization model is based on the linearized sensitivity coefficients between controlled variables (e.g., node voltages) and control variables (e.g., real and reactive power injections of DERs). To avoid the oscillation of solution when it is close to the optimum, a golden search method is introduced to control the step size. Numerical simulations on modified IEEE 13 nodes test feeders show the efficiency of the proposed model. Compared to the results solved by heuristic search (harmony algorithm), the proposed model converges quickly to the global optimum.
We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using a quadratic flux term in the objective function. Four use cases are presented to demonstrate the effectiveness of the approach, including simple and complex stellarator geometries. As shown here, this method outperforms previous two-stage approaches, achieving smaller plasma objective function values when coils are taken into account.
This paper presents an adaptive control approach that modifies a reference model in order to satisfy a multi-objective optimization of two performance signals: a linear signal and a quadratic signal. The time varying reference model modification is accomplished by the real-time solutions of the time-varying Riccati and Sylvester equations coupled with the least-squares parameter estimation of the sensitivities of the performance signals to be minimized. The effectiveness of the proposed method is to be demonstrated by a flight control application of maneuver load alleviation and drag minimization.
The Modified Barrier Functions (MBF) have elements of both Classical Lagrangians (CL) and Classical Barrier Functions (CBF). The MBF methods find an unconstrained minimizer of some smooth barrier function in primal space and then update the Lagrange multipliers, while the barrier parameter either remains fixed or can be updated at each step. The numerical realization of the MBF method leads to the Newton MBF method, where the primal minimizer is found by using Newton's method. This minimizer is then used to update the Lagrange multipliers. In this paper, we examine the Newton MBF method for the Quadratic Programming (QP) problem. It will be shown that under standard second-order optimality conditions, there is a ball around the primal solution and a cut cone in the dual space such that for a set of Lagrange multipliers in this cut cone, the method converges quadratically to the primal minimizer from any point in the aforementioned ball, and continues to do so after each Lagrange multiplier update. The Lagrange multipliers remain within the cut cone and converge linearly to their optimal values. Any point in this ball will be called a "hot start". Starting at such a "hot start", at most Omicron(1n 1n epsilon(exp -1)) Newton steps are sufficient to perform the primal minimization which is necessary for the Lagrange multiplier update. Here, epsilon > 0 is the desired accuracy. Because of the linear convergence of the Lagrange multipliers, this means that only Omicron(1n epsilon(exp -1))omicron(ln 1n epsilon(exp-1)) Newton steps are required to reach an epsilon-approximation to the solution from any "hot start". In order to reach the "hot start", one has to perform Omicron(square root(m) 1n C) Newton steps, where m characterizes the size of the problem and C > 0 is the condition number of the QP problem. This condition number will be characterized explicitly in terms of key parameters of the QP problem, which in turn depend on the input data and the size of the problem.
Finding the optimal position for the individual cells (also called functional modules) on the chip surface is an important and difficult step in the design of integrated circuits. This paper deals with the problem of relative placement, that is the minimization of a quadratic functional with a large, sparse, positive definite system matrix. The basic optimization problem must be augmented by constraints to inhibit solutions where cells overlap. Besides classical iterative methods, based on conjugate gradients (CG), we show that algebraic multigrid methods (AMG) provide an interesting alternative. For moderately sized examples with about 10000 cells, AMG is already competitive with CG and is expected to be superior for larger problems. Besides the classical 'multiplicative' AMG algorithm where the levels are visited sequentially, we propose an 'additive' variant of AMG where levels may be treated in parallel and that is suitable as a preconditioner in the CG algorithm.
The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.
New control synthesis procedures for digital flight control systems were developed. The theoretical developments are the solution to the problem of optimal disturbance suppression in the presence of windshear. Control synthesis is accomplished using a linear quadratic cost function, the command generator tracker for trajectory following and the proportional-integral-filter control structure for practical implementation. Extensions are made to the optimal output feedback algorithm for computing feedback gains so that the multirate and optimal disturbance control designs are computed and compared for the advanced transport operating system (ATOPS). The performance of the designs is demonstrated by closed-loop poles, frequency domain multiinput sigma and eigenvalue plots and detailed nonlinear 6-DOF aircraft simulations in the terminal area in the presence of windshear.
We report that extragalactic foregrounds in cosmic microwave background (CMB) temperature maps lead to significant biases in CMB lensing reconstruction if not properly accounted for. Combinations of multifrequency data have been used to minimize the overall map variance (internal linear combination, or ILC), or specifically null a given foreground, but these are not tailored to CMB lensing. In this paper, we derive an optimal multifrequency combination to jointly minimize CMB lensing noise and bias. We focus on the standard lensing quadratic estimator, as well as the "shear-only"and source-hardened estimators, whose responses to foregrounds differ. We show that an optimal multifrequency combination is a compromise between the ILC and joint deprojection, which nulls the thermal Sunyaev-Zel'dovich (tSZ) and cosmic infrared background (CIB) contributions. In particular, for a Simons Observatory-like experiment with ℓ max,T =3000, we find that profile hardening alone (with the standard ILC) reduces the bias to the lensing power amplitude by 40%, at a 20% cost in noise, while the bias to the cross-correlation with a LSST-like sample is reduced by nearly an order of magnitude at a 10% noise cost, relative to the standard quadratic estimator. With a small amount of joint deprojection the bias to the profile hardened estimator can be further reduced to less than half the statistical uncertainty on the respective amplitudes, at a 20% and 5% noise cost for the auto- and cross-correlation respectively, relative to the profile hardened estimator with the standard ILC weights. Finally, we explore possible improvements with more aggressive masking and varying ℓ max,T .