Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quadratic Optimization”

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.

At least 757 records · Page 42

Leveraging Multitime Hamilton–Jacobi PDEs for Certain Scientific Machine Learning Problems

Hamilton-Jacobi partial differential equations (HJ PDEs) have deep connections with a wide range of fields, including optimal control, differential games, and imaging sciences. By considering the time variable to be a higher dimensional quantity, HJ PDEs can be extended to the multi-time case. In this paper, we establish a novel theoretical connection between specific optimization problems arising in machine learning and the multi-time Hopf formula, which corresponds to a representation of the solution to certain multi-time HJ PDEs. Through this connection, we increase the interpretability of the training process of certain machine learning applications by showing that when we solve these learning problems, we also solve a multi-time HJ PDE and, by extension, its corresponding optimal control problem. As a first exploration of this connection, we develop the relation between the regularized linear regression problem and the Linear Quadratic Regulator (LQR). We then leverage our theoretical connection to adapt standard LQR solvers (namely, those based on the Riccati ordinary differential equations) to design new training approaches for machine learning. Lastly, we provide some numerical examples that demonstrate the versatility and possible computational advantages of our Riccati-based approach in the context of continual learning, post-training calibration, transfer learning, and sparse dynamics identification.

97 MATHEMATICS AND COMPUTING↗

Design of minimax output feedback controller for system with parameter uncertainty

The problem of controlling a time-invariant system with parameter uncertainty is considered with incomplete state feedback. The controller is designed by minimaximizing a quadratic performance criterion and a sensitivity (or loss) criterion, involving the state of the system, the control, and the uncertainty vector. The resulting optimal controller is linear and optimal feedback gain matrix must satisfy a set of nonlinear algebraic equations. some algorithms for algebraic minimax problems are presented.

Basuthakur, S.↗

Active flutter control for flexible vehicles, volume 1

An active flutter control methodology based on linear quadratic gaussian theory and its application to the control of a super critical wing is presented. Results of control surface and sensor position optimization are discussed. Both frequency response matching and residualization used to obtain practical flutter controllers are examined. The development of algorithms and computer programs for flutter modeling and active control design procedures is reported.

Mahesh, J. K.↗

Fabrication of Thermoplastic Composite Laminates Having Film Interleaves By Automated Fiber Placement

Experiments were carried out at the NASA Langley Research Center automated Fiber placement facility to determine an optimal process for the fabrication of composite materials having polymer film interleaves. A series of experiments was conducted to determine an optimal process for the composite prior to investigation of a process to fabricate laminates with polymer films. The results of the composite tests indicated that a well-consolidated, void-free laminate could be attained. Preliminary interleaf processing trials were then conducted to establish some broad guidelines for film processing. The primary finding of these initial studies was that a two-stage process was necessary in order to process these materials adequately. A screening experiment was then performed to determine the relative influence of the process variables on the quality of the film interface as determined by the wedge peel test method. Parameters that were found to be of minor influence on specimen quality were subsequently held at fixed values enabling a more rapid determination of an optimal process. Optimization studies were then performed by varying the remaining parameters at three film melt processing rates. The resulting peel data were fitted with quadratic response surfaces. Additional specimens were fabricated at levels of high peel strength as predicted by the regression models in an attempt to gage the accuracy of the predicted response and to assess the repeatability of the process. The overall results indicate that quality laminates having film interleaves can be successfully and repeatably fabricated by automated fiber placement.

Hulcher, A. B.↗

Representation Learning via Quantum Neural Tangent Kernels

Variational quantum circuits are used in quantum machine learning and variational quantum simulation tasks. Designing good variational circuits or predicting how well they perform for given learning or optimization tasks is still unclear. Here we discuss these problems, analyzing variational quantum circuits using the theory of neural tangent kernels. We define quantum neural tangent kernels, and derive dynamical equations for their associated loss function in optimization and learning tasks. We analytically solve the dynamics in the frozen limit, or lazy training regime, where variational angles change slowly and a linear perturbation is good enough. We extend the analysis to a dynamical setting, including quadratic corrections in the variational angles. We then consider a hybrid quantum classical architecture and define a large-width limit for hybrid kernels, showing that a hybrid quantum classical neural network can be approximately Gaussian. The results presented here show limits for which analytical understandings of the training dynamics for variational quantum circuits, used for quantum machine learning and optimization problems, are possible. These analytical results are supported by numerical simulations of quantum machine-learning experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-J-181: Effect of Prostate Motion On Combined Brachytherapy and External Beam Dose Based On Daily Motion of the Prostate

Purpose: In this study, the adequacy of target expansions on the combined external beam and implant dose was examined based on the measured daily motion of the prostate. Methods: Thirty patients received an I–125 prostate implant prescribed to dose of 90Gy. This was followed by external beam to deliver a dose of 90Gyeq (external beam equivalent) to the prostate over 25 to 30 fractions. An ideal IMRT plan was developed by optimizing the external beam dose based on the delivered implant dose. The implant dose was converted to an equivalent external beam dose using the linear quadratic model. Patients were set up on the treatment table by daily orthogonal imaging and aligning the marker seeds in the prostate. Orthogonal images were obtained at the end of treatment to assess prostate intrafraction motion. Based on the observed motion of the markers between the initial and final images, 5 individual plans showing the actual dose delivered to the patient were calculated. A final true dose distribution was established based on summing the implant dose and the 5 external beam plans. Dose to the prostate, seminal vesicles, lymphnodes and normal tissues, rectal wall, urethra and lower sphincter were calculated and compared to ideal. On 18 patients who were sexually active, dose to the corpus cavernosum and internal pudendal artery was also calculated. Results: The average prostate motion in 3 orthogonal directions was less than 1 mm with a standard deviation of less than +2 mm. Dose and volume parameters showed that there was no decrease in dose to the targets and a marginal decrease in dose to in normal tissues. Conclusion: Dose delivered by seed implant moves with the prostate, decreasing the impact of intrafractions dose movement on actual dose delivered. Combined brachytherapy and external beam dose delivered to the prostate was not sensitive to prostate motion.

60 APPLIED LIFE SCIENCES↗

Singular perturbations and the sounding rocket problem

In this paper, Goddard's problem of maximizing the final altitude of a sounding rocket (a singular problem of optimal control) is analyzed using singular perturbation methods. The problem is first cast in singular perturbation form and then solved to zero order by adding boundary-layer corrections to the reduced solution. For a quadratic drag law, a closed-form solution is obtained, although consideration of a numerical example indicates that this solution is not useful for practical sounding rockets. However, use of state variable transformations allows a very accurate numerical approximation to be constructed. It is concluded that application of singular perturbation methods to the well-known sounding rocket problem indicates that these methods may have utility in dealing with singular problems of optimal control.

Ardema, M. D.↗

Modified optimal control pilot model for computer-aided design and analysis

This paper presents the theoretical development of a modified optimal control pilot model based upon the optimal control model (OCM) of the human operator developed by Kleinman, Baron, and Levison. This model is input compatible with the OCM and retains other key aspects of the OCM, such as a linear quadratic solution for the pilot gains with inclusion of control rate in the cost function, a Kalman estimator, and the ability to account for attention allocation and perception threshold effects. An algorithm designed for each implementation in current dynamic systems analysis and design software is presented. Example results based upon the analysis of a tracking task using three basic dynamic systems are compared with measured results and with similar analyses performed with the OCM and two previously proposed simplified optimal pilot models. The pilot frequency responses and error statistics obtained with this modified optimal control model are shown to compare more favorably to the measured experimental results than the other previously proposed simplified models evaluated.

Davidson, John B.↗

Optimization of cascade blade mistuning. I - Equations of motion and basic inherent properties

Attention is given to the derivation of the equations of motion of mistuned compressor blades, interpolating aerodynamic coefficients by means of quadratic expressions in the reduced frequency. If the coefficients of the quadratic expressions are permitted to assume complex values, excellent accuracy is obtained and Pade rational expressions are obviated. On the basis of the resulting equations, it is shown analytically that the sum of all the real parts of the eigenvalues is independent of the mistuning introduced into the system. Blade mistuning is further treated through the aerodynamic energy approach, and the limiting vibration modes associated with alternative mistunings are identified.

Nissim, E.↗

Balancing Power Absorption Against Structural Loads With Viscous Drag and Power-Takeoff Efficiency Considerations

The focus of this paper is to balance power absorption against structural loading for a novel fixed-bottom oscillating surge wave energy converter in both regular and irregular wave environments. The power-to-load ratio will be evaluated using pseudospectral control (PSC) to determine the optimum power-takeoff (PTO) torque based on a multiterm objective function. This paper extends the pseudospectral optimal control problem to not just maximize the time-averaged absorbed power but also include measures for the surge-foundation force and PTO torque in the optimization. The objective function may now potentially include three competing terms that the optimizer must balance. Separate weighting factors are attached to the surge-foundation force and PTO control torque that can be used to tune the optimizer performance to emphasize either power absorption or load shedding. To correct the pitch equation of motion, derived from linear hydrodynamic theory, a quadratic-viscous-drag torque has been included in the system dynamics; however, to continue the use of quadratic programming solvers, an iteratively obtained linearized drag coefficient was utilized that provided good accuracy in the predicted pitch motion. Furthermore, the analysis considers the use of a nonideal PTO unit to more accurately evaluate controller performance. The PTO efficiency is not directly included in the objective function but rather the weighting factors are utilized to limit the PTO torque amplitudes, thereby reducing the losses resulting from the bidirectional energy flow through a nonideal PTO. Furthermore, results from PSC show that shedding a portion of the available wave energy can lead to greater reductions in structural loads, peak-to-average power ratio, and reactive power requirement.

16 TIDAL AND WAVE POWER↗

Application of multivariable search techniques to the optimization of airfoils in a low speed nonlinear inviscid flow field

Multivariable search techniques are applied to a particular class of airfoil optimization problems. These are the maximization of lift and the minimization of disturbance pressure magnitude in an inviscid nonlinear flow field. A variety of multivariable search techniques contained in an existing nonlinear optimization code, AESOP, are applied to this design problem. These techniques include elementary single parameter perturbation methods, organized search such as steepest-descent, quadratic, and Davidon methods, randomized procedures, and a generalized search acceleration technique. Airfoil design variables are seven in number and define perturbations to the profile of an existing NACA airfoil. The relative efficiency of the techniques are compared. It is shown that elementary one parameter at a time and random techniques compare favorably with organized searches in the class of problems considered. It is also shown that significant reductions in disturbance pressure magnitude can be made while retaining reasonable lift coefficient values at low free stream Mach numbers.

Hague, D. S.↗

On reliable control system designs

A mathematical model for use in the design of reliable multivariable control systems is discussed with special emphasis on actuator failures and necessary actuator redundancy levels. The model consists of a linear time invariant discrete time dynamical system. Configuration changes in the system dynamics are governed by a Markov chain that includes transition probabilities from one configuration state to another. The performance index is a standard quadratic cost functional, over an infinite time interval. The actual system configuration can be deduced with a one step delay. The calculation of the optimal control law requires the solution of a set of highly coupled Riccati-like matrix difference equations. Results can be used for off-line studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗

Closed-loop, pilot/vehicle analysis of the approach and landing task

In the case of approach and landing, it is universally accepted that the pilot uses more than one vehicle response, or output, to close his control loops. Therefore, to model this task, a multi-loop analysis technique is required. The analysis problem has been in obtaining reasonable analytic estimates of the describing functions representing the pilot's loop compensation. Once these pilot describing functions are obtained, appropriate performance and workload metrics must then be developed for the landing task. The optimal control approach provides a powerful technique for obtaining the necessary describing functions, once the appropriate task objective is defined in terms of a quadratic objective function. An approach is presented through the use of a simple, reasonable objective function and model-based metrics to evaluate loop performance and pilot workload. The results of an analysis of the LAHOS (Landing and Approach of Higher Order Systems) study performed by R.E. Smith is also presented.

Anderson, M. R.↗

Scalable Optimization Methods for Distribution Networks With High PV Integration

Here this paper proposes a suite of algorithms to determine the active- and reactive-power setpoints for photovoltaic (PV) inverters in distribution networks. The objective is to optimize the operation of the distribution feeder according to a variety of performance objectives and ensure voltage regulation. In general, these algorithms take a form of the widely studied ac optimal power flow (OPF) problem. For the envisioned application domain, nonlinear power-flow constraints render pertinent OPF problems nonconvex and computationally intensive for large systems. To address these concerns, we formulate a quadratic constrained quadratic program (QCQP) by leveraging a linear approximation of the algebraic power-flow equations. Furthermore, simplification from QCQP to a linearly constrained quadratic program is provided under certain conditions. The merits of the proposed approach are demonstrated with simulation results that utilize realistic PV-generation and load-profile data for illustrative distribution-system test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Isotherm modeling and techno-economic analysis of a TSA moving bed process using a tetraamine-appended MOF for NGCC applications

Tetraamine-appended metal-organic frameworks (MOF) are a new family of amine-functionalized MOF materials that show potential for CO 2 capture from flue gas conditions relevant to natural gas combined cycle (NGCC) power plants. This work presents isotherm modeling of the tetraamine-appended MOF Mg 2 (dobpdc)(3-4-3)(dobpdc 4- =4,4'-dioxidobiphenyl-3,3'dicarboxylate;3-4-3=N,N'bis(3-aminopropyl)-1,4-diaminobutane), process modeling, scale up, and a techno-economic optimization of a moving bed Temperature Swing Adsorption (TSA) process for carbon capture using this sorbent. This MOF exhibits a unique two-step CO 2 adsorption profile in three different pressure ranges. Thus, arctangent-based logistic functions and, quadratic and Langmuir models were employed to represent such isotherm behavior. The results of the isotherm model show good fitting vs the experimental data with an RMSE of 0.41. Here, to model the carbon capture process, the isotherm model was embedded into a moving bed contactor model, and this was used to simulate a TSA CO 2 capture cycle and evaluate cost-optimal designs considering flue gas from a -650 MWe NGCC power plant. The capital cost model consists of CAPEX correlations for reactors, compressors, ducting, etc., while the operating costs include steam, water, chemicals, and electricity (following NETL's quality guidelines for energy systems studies). A techno-economic optimization of the capture system was performed by using NETL's Framework for the Optimization and Quantification of Uncertainty and Surrogates tool (FOQUS). Results suggest that a moving bed carbon capture system with tetraamine-appended MOF can be competitive compared to conventional MEA solvent-based capture processes for NGCC plants, when a heat recovery efficiency of at least 40% is achieved, and MOF materials can be produced at a cost below -$9/kg.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time Domain Identification of an Optimal Control Pilot Model with Emphasis on the Objective Function

A method for the identification of the pilot's control compensation using time domain techniques is proposed. From this information we hope to infer a quadratic cost function, supported by the data, that represents a reasonable expression for the pilot's control objective in the task being performed, or an inferred piloting strategy. The objectives for this method are: (1) obtain a better understanding of the fundamental piloting techniques in complex tasks, such as landing approach; (2) the development of a metric measurable in simulations and flight test that correlate with subjective pilot opinion; and (3) to further validate pilot models and pilot vehicle analysis methods.

Schmidt, D. K.↗

Combined control-structure optimization

An approach for combined control-structure optimization keyed to enhancing early design trade-offs is outlined and illustrated by numerical examples. The approach employs a homotopic strategy and appears to be effective for generating families of designs that can be used in these early trade studies. Analytical results were obtained for classes of structure/control objectives with linear quadratic Gaussian (LQG) and linear quadratic regulator (LQR) costs. For these, researchers demonstrated that global optima can be computed for small values of the homotopy parameter. Conditions for local optima along the homotopy path were also given. Details of two numerical examples employing the LQR control cost were given showing variations of the optimal design variables along the homotopy path. The results of the second example suggest that introducing a second homotopy parameter relating the two parts of the control index in the LQG/LQR formulation might serve to enlarge the family of Pareto optima, but its effect on modifying the optimal structural shapes may be analogous to the original parameter lambda.

Salama, M.↗

Neural Network and Response Surface Methodology for Rocket Engine Component Optimization

The goal of this work is to compare the performance of response surface methodology (RSM) and two types of neural networks (NN) to aid preliminary design of two rocket engine components. A data set of 45 training points and 20 test points obtained from a semi-empirical model based on three design variables is used for a shear coaxial injector element. Data for supersonic turbine design is based on six design variables, 76 training, data and 18 test data obtained from simplified aerodynamic analysis. Several RS and NN are first constructed using the training data. The test data are then employed to select the best RS or NN. Quadratic and cubic response surfaces. radial basis neural network (RBNN) and back-propagation neural network (BPNN) are compared. Two-layered RBNN are generated using two different training algorithms, namely solverbe and solverb. A two layered BPNN is generated with Tan-Sigmoid transfer function. Various issues related to the training of the neural networks are addressed including number of neurons, error goals, spread constants and the accuracy of different models in representing the design space. A search for the optimum design is carried out using a standard gradient-based optimization algorithm over the response surfaces represented by the polynomials and trained neural networks. Usually a cubic polynominal performs better than the quadratic polynomial but exceptions have been noticed. Among the NN choices, the RBNN designed using solverb yields more consistent performance for both engine components considered. The training of RBNN is easier as it requires linear regression. This coupled with the consistency in performance promise the possibility of it being used as an optimization strategy for engineering design problems.

Vaidyanathan, Rajkumar↗