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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Sensitivity Analysis of Linear Programming and Quadratic Programming Algorithms for Control Allocation

The Next Generation (NextGen) transport aircraft configurations being investigated as part of the NASA Aeronautics Subsonic Fixed Wing Project have more control surfaces, or control effectors, than existing transport aircraft configurations. Conventional flight control is achieved through two symmetric elevators, two antisymmetric ailerons, and a rudder. The five effectors, reduced to three command variables, produce moments along the three main axes of the aircraft and enable the pilot to control the attitude and flight path of the aircraft. The NextGen aircraft will have additional redundant control effectors to control the three moments, creating a situation where the aircraft is over-actuated and where a simple relationship does not exist anymore between the required effector deflections and the desired moments. NextGen flight controllers will incorporate control allocation algorithms to determine the optimal effector commands and attain the desired moments, taking into account the effector limits. Approaches to solving the problem using linear programming and quadratic programming algorithms have been proposed and tested. It is of great interest to understand their relative advantages and disadvantages and how design parameters may affect their properties. In this paper, we investigate the sensitivity of the effector commands with respect to the desired moments and show on some examples that the solutions provided using the l2 norm of quadratic programming are less sensitive than those using the l1 norm of linear programming.

Frost, Susan A.↗
DOE OSTI2017

Causal control design for wave energy converters with finite stroke

Control design is of central importance in the optimization of power generation from ocean wave energy converters (WECs). The problem is made challenging by the fact that the resource is vibratory and stochastic, with a low quality factor. It has recently been shown that when constrained to causality, the optimal feedback controller for a WEC is solvable in closed-form as a sign-indefinite LQG problem, for the case in which the WEC exhibits linear dynamics and quadratic power transmission dissipation. However, most WECs exhibit stroke limits, which must be explicitly accounted for in control design to avoid end-stop collisions. This is important because such collisions lead to damage, and also because the introduction of unmodeled stroke limits can actually destabilize the otherwise-optimal closed loop. We illustrate a control design technique to address this issue, which consists of two design steps. In the first step, LMI techniques are used to optimize the mean power harvested from waves with a known power spectral density, subject to several competing objectives. In the second step, the controller is augmented via a nonlinear passivity-based outer-loop design which guarantees to prevent stroke saturation while also guaranteeing to preserve closed-loop stability. Furthermore, the technique is illustrated using a model of a simple heaving point absorber, and verified in simulation.

16 TIDAL AND WAVE POWER↗
NASA NTRSDate not supplied

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗
DOE OSTI2017

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Optical neural network system for pose determination of spinning satellites

An optical neural network architecture and algorithm based on a Hopfield optimization network are presented for multitarget tracking. This tracker utilizes a neuron for every possible target track, and a quadratic energy function of neural activities which is minimized using gradient descent neural evolution. The neural net tracker is demonstrated as part of a system for determining position and orientation (pose) of spinning satellites with respect to a robotic spacecraft. The input to the system is time sequence video from a single camera. Novelty detection and filtering are utilized to locate and segment novel regions from the input images. The neural net multitarget tracker determines the correspondences (or tracks) of the novel regions as a function of time, and hence the paths of object (satellite) parts. The path traced out by a given part or region is approximately elliptical in image space, and the position, shape and orientation of the ellipse are functions of the satellite geometry and its pose. Having a geometric model of the satellite, and the elliptical path of a part in image space, the three-dimensional pose of the satellite is determined. Digital simulation results using this algorithm are presented for various satellite poses and lighting conditions.

Lee, Andrew↗
NASA NTRS1990

Second-order design sensitivities for linear elastic problems

The subject of this paper is the efficient and accurate determination of second-order design sensitivities in elastic bodies. The approach being carried out here is the direct differentiation of the governing derivative boundary element method (DBEM) formulation of the problem. Second-order sensitivities of boundary stresses are obtained in an elegant manner. A numerical implementation of the method is carried out with isoparametric quadratic boundary elements and numerical results are presented for several sample problems. Considerable savings in computational effort for an optimization procedure is possible through the use of efficiently determined accurate second-order design sensitivities.

Zhang, Qing↗
NASA NTRS1989

Integrated control-system design via generalized LQG (GLQG) theory

Thirty years of control systems research has produced an enormous body of theoretical results in feedback synthesis. Yet such results see relatively little practical application, and there remains an unsettling gap between classical single-loop techniques (Nyquist, Bode, root locus, pole placement) and modern multivariable approaches (LQG and H infinity theory). Large scale, complex systems, such as high performance aircraft and flexible space structures, now demand efficient, reliable design of multivariable feedback controllers which optimally tradeoff performance against modeling accuracy, bandwidth, sensor noise, actuator power, and control law complexity. A methodology is described which encompasses numerous practical design constraints within a single unified formulation. The approach, which is based upon coupled systems or modified Riccati and Lyapunov equations, encompasses time-domain linear-quadratic-Gaussian theory and frequency-domain H theory, as well as classical objectives such as gain and phase margin via the Nyquist circle criterion. In addition, this approach encompasses the optimal projection approach to reduced-order controller design. The current status of the overall theory will be reviewed including both continuous-time and discrete-time (sampled-data) formulations.

Bernstein, Dennis S.↗
DOE OSTI2019

Optimizing multi-modal, non-axisymmetric plasma response metrics with additional coil rows on DIII-D

The expansion of the DIII-D tokamak operational space via the addition of new non-axisymmetric coil rows is investigated using the linear MHD plasma response code MARS-F. The plasma response to magnetic perturbations imposed by five distinct low- and high-field side coil geometries is modeled for toroidal mode numbers up to n = 6. A set of linear (~δB^2: resonant, kink mode amplitude) and quadratic (~δB: core and edge integrated NTV torque) plasma metrics are defined and compared directly for the various coil geometries. Coil size is also scanned as a secondary optimization variable. An in-vessel coil located at the low-field side midplane is identified as the most efficient plasma response actuator across all metrics and toroidal mode numbers. The multi-row combination of this new midplane (M) coil with the existing DIII-D non-axisymmetric coil set is considered in numerical optimization studies to quantify the impact of this new coil on modeled plasma response metrics. The addition of the M coil is shown to significantly increase the actuator efficiency in maximizing individual plasma response metrics as well as maximizing multi-modal ratios of response metrics, e.g. maximizing resonant modes or edge NTV torque while minimizing kink modes or core NTV torque. The ex-vessel midplane (C) coil is also observed to have an outsize effect on multi-modal optimization due to its uniquely narrow poloidal spectra, and the expansion of optimization space due to the combined use of the two existing in-vessel coil rows with both the M and C coils is notably larger than either coil's individual contribution even though both share similar poloidal locations. This four-row system is expected to greatly extend abilities to understand and optimize the control of transients and rotation on the DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Optimal control laws for microburst encounter

Simplified structures for longitudinal control laws that reduce an aircraft's response to the strong head-tailwind and downdraft variations associated with microbursts are presented. They are based on non-zero-set-point linear-quadratic regulators that command throttle setting and angle of attack as functions of velocity and flight path angle, and they can incorporate direct measurements of the wind profile if available. Selection of cost functions to be minimized by feedback control has been aided by a prior study of classical control laws and exact nonlinear-optimal flight paths through realistic microburst wind profiles. The resulting optimal control laws have an adaptive, dual-mode structure that can be implemented either in flight-director logic or in an autopilot.

Stengel, R. F.↗
NASA NTRS1972

Davidon-Broyden rank-one minimization methods in Hilbert space with application to optimal control problems

The Davidon-Broyden class of rank one, quasi-Newton minimization methods is extended from Euclidean spaces to infinite-dimensional, real Hilbert spaces. For several techniques of choosing the step size, conditions are found which assure convergence of the associated iterates to the location of the minimum of a positive definite quadratic functional. For those techniques, convergence is achieved without the problem of the computation of a one-dimensional minimum at each iteration. The application of this class of minimization methods for the direct computation of the solution of an optimal control problem is outlined. The performance of various members of the class are compared by solving a sample optimal control problem. Finally, the sample problem is solved by other known gradient methods, and the results are compared with those obtained with the rank one quasi-Newton methods.

Straeter, T. A.↗
NASA NTRS1984

Application of optimization techniques to the design of a flutter suppression control law for the DAST ARW-2

The design of a candidate flutter suppression (FS) control law for the symmetric degrees of freedom for the DAST ARW-2 aircraft is discussed. The results illustrate the application of several currently employed control law design techniques. Subsequent designs, obtained as the mathematical model of the ARW-2 is updated, are expected to employ similar methods and to provide a control law whose performance will be flight tested. This study represents one of the steps necessary to provide an assessment of the validity of applying current control law synthesis and analysis techniques in the design of actively controlled aircraft. Mathematical models employed in the control law design and evaluation phases are described. The control problem is specified by presenting the flutter boundary predicted for the uncontrolled aircraft and by defining objectives and constraints that the controller should satisfy. A full-order controller is obtained by using Linear Quadratic Gaussian (LQG) techniques. The process of obtaining an implementable reduced-order controller is described. One example is also shown in which constrained optimization techniques are utilized to explicitly include robustness criteria within the design algorithm.

Adams, W. M., Jr.↗
DOE OSTI2017

Effect of nano-SiO{sub 2} particles and curing time on development of fiber-matrix bond properties and microstructure of ultra-high strength concrete

Bond properties between fibers and cementitious matrix have significant effect on the mechanical behavior of composite materials. In this study, the development of steel fiber-matrix interfacial bond properties in ultra-high strength concrete (UHSC) proportioned with nano-SiO{sub 2} varying between 0 and 2%, by mass of cementitious materials, was investigated. A statistical model relating either bond strength or pullout energy to curing time and nano-SiO{sub 2} content was proposed by using the response surface methodology. Mercury intrusion porosimetry (MIP) and backscatter scanning electron microscopy (BSEM) were used to characterize the microstructure of the matrix and the fiber-matrix interface, respectively. Micro-hardness around the embedded fiber and hydration products of the matrix were evaluated as well. Test results indicated that the optimal nano-SiO{sub 2} dosage was 1% in terms of the bond properties and the microstructure. The proposed quadratic model efficiently predicted the bond strength and pullout energy with consideration of curing time and nano-SiO{sub 2} content. The improvement in bond properties associated with nano-silica was correlated with denser matrix and/or interface and stronger bond and greater strength of hydration products based on microstructural analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Self‐Potential Tomography Preconditioned by Particle Swarm Optimization—Application to Monitoring Hyporheic Exchange in a Bedrock River

Abstract A self‐potential (SP) data‐inversion algorithm was developed and tested on an analytical model of electrical‐potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP‐monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock‐lined riverbed. The algorithm combined stochastic source‐localization by particle‐swarm‐optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least‐squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical‐model data to evaluate the algorithm performance. Results indicated that true parameters of single‐source models were inverted and approximated with small residual error, whereas inversion of analytical‐model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non‐uniqueness of the inverse‐model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point‐scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP‐monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Evidence for a universal Fermi-liquid scattering rate throughout the phase diagram of the copper-oxide superconductors

Abstract The phase diagram of the cuprate superconductors continues to pose formidable scientific challenges. While these materials are typically viewed as doped Mott insulators, it is well known that they are Fermi liquids at high hole-dopant concentrations. It was recently demonstrated that at moderate doping, in the pseudogap (PG) region of the phase diagram, the charge carriers are also best described as a Fermi liquid. Nevertheless, the relationship between the two Fermi-liquid (FL) regions and the nature of the strange-metal (SM) state at intermediate doping have remained unsolved. Here we show for the case of the model cuprate superconductor HgBa 2 CuO 4+ δ that the normal-state transport scattering rate determined from the cotangent of the Hall angle remains quadratic in temperature across the PG temperature, upon entering the SM state, and that it is doping-independent below optimal doping. Analysis of prior transport results for other cuprates reveals that this behavior is universal throughout the entire phase diagram and points to a pervasive FL transport scattering rate. These observations can be reconciled with a variety of other experimental results for the cuprates upon considering the possibility that the PG phenomenon is associated with the gradual, non-uniform localization of one hole per planar CuO 2 unit.

Barišić, N.↗
NASA NTRS1991

Optimal digital redesign of continuous-time controllers

This paper proposes a new optimal digital redesign technique for finding a dynamic digital control law from the available analog counterpart and simultaneously minimizing a quadratic performance index. The proposed technique can be applied to a system with a more general class of reference inputs, and the developed digital regulator can be implemented using low cost microcomputers.

Shieh, Leang S.↗
NASA NTRS1989

Optimization of structure and control system

The objective of this study is the simultaneous design of the structural and control system for space structures. This study is focused on considering the effect of the number and the location of the actuators on the minimum weight of the structure, and the total work done by the actuators for specified constraints and disturbance. The controls approach used is the linear quadratic regulator theory with constant feedback. At the beginning collocated actuators and sensors are provided in all the elements. The actuator doing the least work is removed one at a time, and the structure is optimized for the specified constraints on the closed-loop eigenvalues and the damping parameters. The procedure of eliminating an actuator is continued until an acceptable design satisfying the constraints is obtained. The study draws some conclusions on the trade between the total work done by the actuators, and the optimum weight and the number of actuators.

Khot, N. S.↗
NASA NTRS2019

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗
DOE OSTI2019

Improving small-scale CMB lensing reconstruction

Over the past decade, the gravitational lensing of the Cosmic Microwave Background (CMB) has become a powerful tool for probing the matter distribution in the Universe. The standard technique used to reconstruct the CMB lensing signal employs the quadratic estimator (QE) method, which has recently been shown to be suboptimal for lensing measurements on very small scales in temperature and polarization data. We implement a simple, more optimal method for the small-scale regime, which involves taking the direct inverse of the background gradient. We derive new techniques to make continuous maps of lensing using this "Gradient-Inversion" (GI) method and validate our method with simulated data, finding good agreement with predictions. For idealized simulations of lensing cross- and autospectra that neglect foregrounds, we demonstrate that our method performs significantly better than previous quadratic estimator methods in temperature; at $L=5000-9000$, it reduces errors on the lensing auto-power spectrum by a factor of $\sim 4$ for both idealized CMB-S4 and Simons Observatory-like experiments and by a factor of $\sim 2.6$ for cross-correlations of CMB-S4-like lensing reconstruction and the true lensing field. Finally, we caution that the level of the neglected small-scale foreground power, while low in polarization, is very high in temperature; though we briefly outline foreground mitigation methods, further work on this topic is required. Nevertheless, our results show the future potential for improved small-scale CMB lensing measurements, which could provide stronger constraints on cosmological parameters and astrophysics at high redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗