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 739 records · Page 41

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discrete Time McKean–Vlasov Control Problem: A Dynamic Programming Approach

We consider the stochastic optimal control problem of nonlinear mean-field systems in discrete time. We reformulate the problem into a deterministic control problem with marginal distribution as controlled state variable, and prove that dynamic programming principle holds in its general form. We apply our method for solving explicitly the mean-variance portfolio selection and the multivariate linear-quadratic McKean–Vlasov control problem.

97 MATHEMATICS AND COMPUTING↗

Optimization of the fractionated irradiation scheme considering physical doses to tumor and organ at risk based on dose–volume histograms

Purpose: Radiotherapy of solid tumors has been performed with various fractionation regimens such as multi- and hypofractionations. However, the ability to optimize the fractionation regimen considering the physical dose distribution remains insufficient. This study aims to optimize the fractionation regimen, in which the authors propose a graphical method for selecting the optimal number of fractions (n) and dose per fraction (d) based on dose–volume histograms for tumor and normal tissues of organs around the tumor. Methods: Modified linear-quadratic models were employed to estimate the radiation effects on the tumor and an organ at risk (OAR), where the repopulation of the tumor cells and the linearity of the dose-response curve in the high dose range of the surviving fraction were considered. The minimization problem for the damage effect on the OAR was solved under the constraint that the radiation effect on the tumor is fixed by a graphical method. Here, the damage effect on the OAR was estimated based on the dose–volume histogram. Results: It was found that the optimization of fractionation scheme incorporating the dose–volume histogram is possible by employing appropriate cell surviving models. The graphical method considering the repopulation of tumor cells and a rectilinear response in the high dose range enables them to derive the optimal number of fractions and dose per fraction. For example, in the treatment of prostate cancer, the optimal fractionation was suggested to lie in the range of 8–32 fractions with a daily dose of 2.2–6.3 Gy. Conclusions: It is possible to optimize the number of fractions and dose per fraction based on the physical dose distribution (i.e., dose–volume histogram) by the graphical method considering the effects on tumor and OARs around the tumor. This method may stipulate a new guideline to optimize the fractionation regimen for physics-guided fractionation.

60 APPLIED LIFE SCIENCES↗

In-flight adaptive performance optimization (APO) control using redundant control effectors of an aircraft

Practical application of real-time (or near real-time) Adaptive Performance Optimization (APO) is provided for a transport aircraft in steady climb, cruise, turn descent or other flight conditions based on measurements and calculations of incremental drag from a forced response maneuver of one or more redundant control effectors defined as those in excess of the minimum set of control effectors required to maintain the steady flight condition in progress. The method comprises the steps of applying excitation in a raised-cosine form over an interval of from 100 to 500 sec. at the rate of 1 to 10 sets/sec of excitation, and data for analysis is gathered in sets of measurements made during the excitation to calculate lift and drag coefficients C.sub.L and C.sub.D from two equations, one for each coefficient. A third equation is an expansion of C.sub.D as a function of parasitic drag, induced drag, Mach and altitude drag effects, and control effector drag, and assumes a quadratic variation of drag with positions .delta..sub.i of redundant control effectors i=1 to n. The third equation is then solved for .delta..sub.iopt the optimal position of redundant control effector i, which is then used to set the control effector i for optimum performance during the remainder of said steady flight or until monitored flight conditions change by some predetermined amount as determined automatically or a predetermined minimum flight time has elapsed.

Gilyard, Glenn B.↗

Global optimal consensus for discrete-time multi-agent systems with bounded controls

In this paper, we consider the global optimal consensus problem for discrete-time multi-agent systems with bounded control protocols over xed undirected communication networks. Each agent is described by a discrete-time single integrator and endowed with a quadratic objective function which is private to itself. For each agent, we develop two bounded distributed protocols: bounded proportional-integral (PI) protocol and bounded integral (I) protocol, based on information receiving from its neighboring agents through the communication network and the gradient of its own objective function. We show that the proposed bounded distributed protocols with properly chosen parameters solve the global optimal consensus problem if the communication network is connected, i.e., the multi-agent system achieves consensus at a state that minimizes the sum of the objective functions of all agents. The theoretical results are illustrated by numerical simulations.

Consensus, multi-agent system, optimal control↗

A Comfort Model Simplification for Tight Integration with Grid Service Optimizations

Localized generation, storage, and intelligent power electronics are increasingly being integrated into building. The excess storage, thermal, and generation capacity can be used to develop new markets for ancillary services such as peak reduction, voltage and frequency support, or load flattening. Widespread participation in these markets will make the power grid more resilient. However, widespread participation in these markets will require a simple and intuitive control over the excess capacity by the building occupants. Most importantly for adopting these technologies, participation in these markets should not require sacrificing occupant comfort. Many of the proposed control and optimization technologies for supplying the ancillary services do not explicitly incorporate occupant comfort into their models and when they do, generally utilize a rigid constraint on the indoor air temperature set a priori. This reduces flexibility in using excess thermal capacity. It also cannot take into account uncertainties in local factors that affect thermal comfort such as clothing insulation or relative humidity. In this paper we present and justify some assumptions that simplify the standard predictors of thermal comfort. We also develop a regression model of the thermal comfort that is linear with respect to the indoor air temperature. Next we develop the concept of thermal comfort variation that eliminates the dependence of the thermal comfort on uncertain thermal comfort factors. Then, we present a quadratic program utilizing the thermal comfort variation in both the objective function and as a constraint. Finally, we show the results of a time-of-use cost optimization that utilizes the simplified thermal comfort model.

Melin, Alexander↗

A reliable algorithm for optimal control synthesis

In recent years, powerful design tools for linear time-invariant multivariable control systems have been developed based on direct parameter optimization. In this report, an algorithm for reliable optimal control synthesis using parameter optimization is presented. Specifically, a robust numerical algorithm is developed for the evaluation of the H(sup 2)-like cost functional and its gradients with respect to the controller design parameters. The method is specifically designed to handle defective degenerate systems and is based on the well-known Pade series approximation of the matrix exponential. Numerical test problems in control synthesis for simple mechanical systems and for a flexible structure with densely packed modes illustrate positively the reliability of this method when compared to a method based on diagonalization. Several types of cost functions have been considered: a cost function for robust control consisting of a linear combination of quadratic objectives for deterministic and random disturbances, and one representing an upper bound on the quadratic objective for worst case initial conditions. Finally, a framework for multivariable control synthesis has been developed combining the concept of closed-loop transfer recovery with numerical parameter optimization. The procedure enables designers to synthesize not only observer-based controllers but also controllers of arbitrary order and structure. Numerical design solutions rely heavily on the robust algorithm due to the high order of the synthesis model and the presence of near-overlapping modes. The design approach is successfully applied to the design of a high-bandwidth control system for a rotorcraft.

Vansteenwyk, Brett↗

Optimizing supercritical carbon dioxide in the inactivation of bacteria in clinical solid waste by using response surface methodology

Highlights: • Supercritical carbon dioxide sterilization of clinical solid waste. • Inactivation of bacteria in clinical solid waste using supercritical carbon dioxide. • Reduction of the hazardous exposure of clinical solid waste. • Optimization of the supercritical carbon dioxide experimental conditions. - Abstract: Clinical solid waste (CSW) poses a challenge to health care facilities because of the presence of pathogenic microorganisms, leading to concerns in the effective sterilization of the CSW for safe handling and elimination of infectious disease transmission. In the present study, supercritical carbon dioxide (SC-CO{sub 2}) was applied to inactivate gram-positive Staphylococcus aureus, Enterococcus faecalis, Bacillus subtilis, and gram-negative Escherichia coli in CSW. The effects of SC-CO{sub 2} sterilization parameters such as pressure, temperature, and time were investigated and optimized by response surface methodology (RSM). Results showed that the data were adequately fitted into the second-order polynomial model. The linear quadratic terms and interaction between pressure and temperature had significant effects on the inactivation of S. aureus, E. coli, E. faecalis, and B. subtilis in CSW. Optimum conditions for the complete inactivation of bacteria within the experimental range of the studied variables were 20 MPa, 60 °C, and 60 min. The SC-CO{sub 2}-treated bacterial cells, observed under a scanning electron microscope, showed morphological changes, including cell breakage and dislodged cell walls, which could have caused the inactivation. This espouses the inference that SC-CO{sub 2} exerts strong inactivating effects on the bacteria present in CSW, and has the potential to be used in CSW management for the safe handling and recycling-reuse of CSW materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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.↗