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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Modeling and analysis of mixed convection stagnation point flow of nanofluid towards a stretching surface: OHAM and FEM approach

A theoretical study of mixed convection stagnation point flow towards a stretching surface is presented. The governing boundary layer equations are transformed into a set of highly nonlinear ordinary differential equations using suitable similarity transforms. The semi-analytical solution is obtained using optimal homotopy analysis method (OHAM) and the numerical solution is obtained via finite element method (FEM). Solutions obtained via two different approaches are in excellent agreement, which validates the accuracy of present analysis. In a special case, the present OHAM solution is also validated with the earlier available results. Effect of pertinent flow parameters on the skin friction coefficient and Nusselt number is presented in tabular form, whereas the velocity, temperature and nanoparticle distribution are presented in graphical forms. Further, a quadratic multiple regression analysis on numeric data of skin friction coefficient and Nusselt number is performed. The findings suggest that velocity slip assists the fluid motion in presence of buoyancy forces, whereas it exhibits a retarding nature on fluid motion when no buoyancy forces exist.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A State-Space Model of an Inverter-Based Microgrid for Multivariable Feedback Control Analysis and Design

In this work, a synchronous model for grid-connected and islanded microgrids is presented. The grid-connected model is based on the premise that the reference frame is synchronized with the AC bus. The quadrature component of the AC bus voltage can be cancelled, which allows to express output power as a linear equation for nominal values in the AC bus amplitude voltage. The model for the islanded microgrid is developed by integrating all the inverter dynamics using a state-space model for the load currents. This model is presented in a comprehensive way such that it could be scalable to any number of inverter-based generators using inductor–capacitor–inductor (LCL) output filters. The use of these models allows designers to assess microgrid stability and robustness using modern control methods such as eigenvalue analysis and singular value diagrams. Both models were tested and validated in an experimental setup to demonstrate their accuracy in describing microgrid dynamics. In addition, three scenarios are presented: non-controlled model, Linear-Quadratic Integrator (LQI) power control, and Power-Voltage (PQ/Vdq) droop–boost controller. Experimental results demonstrate the effectiveness of the control strategies and the accuracy of the models to describe microgrid dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1987

Control of Spacecraft Control Laboratory Experiment (SCOLE) slew maneuvers

This is the second report of a set of two reports on the dynamics and control of slewing maneuvers of NASA Spacecraft Control Laboratory Experiment (SCOLE). The control problem of slewing maneuvers of SCOLE is developed in terms of an arbitrary maneuver about any given axis. The control system is developed for the combined problem of rigid-body slew maneuver and vibration suppression of flexible appendage. The control problem is formulated by incorporating the nonlinear equations derived in the previous report and is expressed in terms of a two-point boundary value problem utilizing a quadratic type of performance index. The two-point boundary value problem is solved as a hierarchical control problem with the overall system being split in terms of two subsystems, namely the slewing of the entire assembly and the vibration suppression of the flexible antenna. The coupling variables between the two dynamical subsystems are identified and these two subsystems for control purposes are treated independently in parallel at the first level. Then the state-space trajectory of the combined problem is optimized at the second level.

Kakad, Y. P.↗
NASA NTRS1989

Dynamics and control of flexible spacecraft during and after slewing maneuvers

The dynamics and control of slewing maneuvers of NASA Spacecraft COntrol Laboratory Experiment (SCOLE) are analyzed. The control problem of slewing maneuvers of SCOLE is formulated in terms of an arbitrary maneuver about any given axis. The control system is developed for the combined problem of rigid-body slew maneuver and vibration suppression of the flexible appendage. The control problem formulation incorporates the nonlinear dynamical equations derived previously, and is expressed in terms of a two-point boundary value problem utilizing a quadratic type of performance index. The two-point boundary value problem is solved as a hierarchical control problem with the overall system being split in terms of two subsystems, namely the slewing of the entire assembly and the vibration suppression of the flexible antenna. The coupling variables between the two dynamical subsystems are identified and these two subsystems for control purposes are treated independently in parallel at the first level. Then the state-space trajectory of the combined problem is optimized at the second level.

Kakad, Yogendra P.↗
DOE OSTI2026

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗
DOE OSTI2016

On the constrained classical capacity of infinite-dimensional covariant quantum channels

The additivity of the minimal output entropy and that of the χ-capacity are known to be equivalent for finite-dimensional irreducibly covariant quantum channels. In this paper, we formulate a list of conditions allowing to establish similar equivalence for infinite-dimensional covariant channels with constrained input. This is then applied to bosonic Gaussian channels with quadratic input constraint to extend the classical capacity results of the recent paper [Giovannetti et al., Commun. Math. Phys. 334(3), 1553-1571 (2015)] to the case where the complex structures associated with the channel and with the constraint operator need not commute. In particular, this implies a multimode generalization of the “threshold condition,” obtained for single mode in Schäfer et al. [Phys. Rev. Lett. 111, 030503 (2013)], and the proof of the fact that under this condition the classical “Gaussian capacity” resulting from optimization over only Gaussian inputs is equal to the full classical capacity. Complex structures correspond to different squeezings, each with its own normal modes, vacuum and coherent states, and the gauge. Thus our results apply, e.g., to multimode channels with a squeezed Gaussian noise under the standard input energy constraint, provided the squeezing is not too large as to violate the generalized threshold condition. We also investigate the restrictiveness of the gauge-covariance condition for single- and multimode bosonic Gaussian channels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Electrooptical effect in the plasmon structure glass–In{sub 2}O{sub 3}: Sn–ferroelectric–Al with a subwavelength grating

The spectral features of the electrooptical effect, in which an applied voltage changes the refractive index of a ferroelectric copolymer, are studied. Three nanostructures are investigated: two nanostructures play an auxiliary role, and the basic structure consists of a glass substrate, a transparent ITO (In{sub 2}O{sub 3}: Sn) layer, active copolymer layer, and an Al layer with a nanograting. This grating with a period of 400 nm meets the conditions of excitation of plasmon resonances. The light transmission coefficients of all structures are analyzed in the spectral range 400–900 nm. The transmission spectra have two characteristic plasmon dips, one of which is related to the copolymer–ITO interface and the other, to the copolymer–Al interface. The aluminum and ITO layers play the role of electrodes, which supply voltage pulses (from 0 to 15 V) to the copolymer layer. When studying the electrooptical effect, we detected spectral shifts in plasmon resonance bands when the amplitudes of both positive and negative voltage pulses (quadratic effect) increase. These shifts change the effective refractive indices (n{sub eff}) of the structural elements, reaching the minimum negative increment Δ{sub eff} =–0.06.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

E-PVT: enhanced position-velocity-time scheduler for computer-controlled optical finishing with comprehensive considerations of dynamics constraints, continuity and efficiency

Deterministic computer-controlled optical finishing is an essential approach for achieving high-quality optical surfaces. Its determinism and convergence rely heavily on precise and smooth motion control to guide the machine tool over an optical surface to correct residual errors. One widely supported and smooth motion control model is position-velocity-time (PVT), which employs piecewise cubic polynomials to describe positions. Our prior research introduced a PVT-based velocity scheduling method, demonstrating sub-nanometer level convergence in ion beam figuring (IBF) processes. However, three challenges remained. Firstly, this method relies on quadratic programming, resulting in computational intensiveness for dense tool paths. Secondly, the dynamics constraints and velocity and acceleration continuities are not comprehensively considered, limiting the full potential of PVT-based control. Thirdly, no compensation mechanism existed when dynamics constraints are exceeded. In this study, in response to these challenges, we proposed the Enhanced PVT (E-PVT) method, reducing the time complexity from O ( n 3 ) to O ( n ) while fully addressing dynamics constraints and continuities. A novel compensation method utilizing particle swarm optimization was proposed to address situations where dynamics constraints might be exceeded while maintaining the overall processing efficiency. Validation through simulation and experimentation confirmed the improved performance of E-PVT.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Pilot/vehicle control optimization using averaged operational mode and subsystem relative performance index sensitivities

A method is presented for designing optimal feedback controllers for systems having subsystem sensitivity constraints. Such constraints reflect the presence of subsystem performance indices which are in conflict with the performance index of the overall system. The key to the approach is the use of relative performance index sensitivity (a measure of the deviation of a performance index from its optimum value). The weighted sum of subsystem and/or operational mode relative performance index sensitivies is defined as an overall performance index. A method is developed to handle linear systems with quadratic performance indices and either full or partial state feedback. The usefulness of this method is demonstrated by applying it to the design of a stability augmentation system (SAS) for a VTOL aircraft. A desirable VTOL SAS design is one that produces good VTOL transient response both with and without active pilot control. The system designed using this method is shown to effect a satisfactory compromise solution to this problem.

Leininger, G. G.↗
DOE OSTI2022

Data-Driven Optimal Control of Traffic Signals for Urban Road Networks

This paper studies the issue of data-driven optimal control design for traffic signals of oversaturated urban road networks. The signal control system based on the store and forward model is generally uncontrollable for which the controllable decomposition is needed. Instead of identifying the unknown parameters like saturation rates and turning ratios, a finite number of measured trajectories can be used to parametrize the system and help directly construct a transformation matrix for Kalman controllable decomposition through the fundamental lemma of J. C. Willems. On top of that, an infinite-horizon linear quadratic regulator (LQR) problem is formulated considering the constraints of green times for traffic signals. The problem can be solved through a two-phase data-driven learning process, where one solves an infinite-horizon unconstrained LQR problem and the other solves a finite-horizon constrained LQR problem. The simulation result shows the theoretical analysis is effective and the proposed data-driven controller can yield desired performance for reducing traffic congestion.

Liu, Tong↗
DOE OSTI2024

Simultaneous control of the electron temperature and safety factor profiles in DIII-D using model-based optimal control techniques

Future tokamak power plants will likely operate using a single, well-defined plasma scenario, either in steady state or for very long pulse lengths. In order to enhance the robustness of the scenario, feedback controllers for a variety of plasma properties will be necessary to counteract any disturbances and ensure safe operation. However, only a limited set of actuators will be available to control many different quantities. Because of this, it is necessary to develop controllers that are able to regulate multiple plasma properties using a limited set of actuators. To this end, a controller has been developed for the simultaneous regulation of both the electron temperature and safety factor profiles in DIII-D. This algorithm uses a linear quadratic integral control synthesis approach based on a linearized model of the dynamics of the two profiles. Two neural network surrogate models, NubeamNet and MMMnet, are included to improve the fidelity of the model. Furthermore, the controller has been tested in simulation using COTSIM, and has demonstrated the ability to simultaneously track changes in both the electron temperature and safety factor targets, including changes in both the magnitude and the shape of the profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantum-Inspired Power System Reliability Assessment

To enable an in-depth study of power system operation and planning, the assessment of standard reliability indices is inevitable. The Monte Carlo Simulation (MCS) approach is a broadly used method in replacing the analytical methods in reliability indices assessment. The accuracy of MCS, however, highly depends on the sampling size, and hence, a complicated system with large number of components requires a large sampling size and daunting computational effort. To address this shortcoming, we, in this paper attempt to take advantage of potentials of the quantum computing (QC) for power system reliability assessment by realizing the following contributions: 1) an innovative quantum model designed for reliability assessment; 2) a quantum circuit that achieves the quadratic speed up compared to the classical MCS method; 3) an efficient quantum amplitude estimation (QAE) algorithm to accurately evaluate the reliability indices. The accuracy and efficacy of the quantum reliability method are extensively verified and demonstrated on both radial and mesh distribution systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2013

Repopulation of interacting tumor cells during fractionated radiotherapy: Stochastic modeling of the tumor control probability

Purpose: Optimal treatment planning for fractionated external beam radiation therapy requires inputs from radiobiology based on recent thinking about the “five Rs” (repopulation, radiosensitivity, reoxygenation, redistribution, and repair). The need is especially acute for the newer, often individualized, protocols made feasible by progress in image guided radiation therapy and dose conformity. Current stochastic tumor control probability (TCP) models incorporating tumor repopulation effects consider “stem-like cancer cells” (SLCC) to be independent, but the authors here propose that SLCC-SLCC interactions may be significant. The authors present a new stochastic TCP model for repopulating SLCC interacting within microenvironmental niches. Our approach is meant mainly for comparing similar protocols. It aims at practical generalizations of previous mathematical models. Methods: The authors consider protocols with complete sublethal damage repair between fractions. The authors use customized open-source software and recent mathematical approaches from stochastic process theory for calculating the time-dependent SLCC number and thereby estimating SLCC eradication probabilities. As specific numerical examples, the authors consider predicted TCP results for a 2 Gy per fraction, 60 Gy protocol compared to 64 Gy protocols involving early or late boosts in a limited volume to some fractions. Results: In sample calculations with linear quadratic parameters α = 0.3 per Gy, α/β = 10 Gy, boosting is predicted to raise TCP from a dismal 14.5% observed in some older protocols for advanced NSCLC to above 70%. This prediction is robust as regards: (a) the assumed values of parameters other than α and (b) the choice of models for intraniche SLCC-SLCC interactions. However, α = 0.03 per Gy leads to a prediction of almost no improvement when boosting. Conclusions: The predicted efficacy of moderate boosts depends sensitively on α. Presumably, the larger values of α are the ones appropriate for individualized treatment protocols, with the smaller values relevant only to protocols for a heterogeneous patient population. On that assumption, boosting is predicted to be highly effective. Front boosting, apart from practical advantages and a possible advantage as regards iatrogenic second cancers, also probably gives a slightly higher TCP than back boosting. If the total number of SLCC at the start of treatment can be measured even roughly, it will provide a highly sensitive way of discriminating between various models and parameter choices. Updated mathematical methods for calculating repopulation allow credible generalizations of earlier results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Spatial Study 2021: Sensor-Based Time Series of Surface Water Temperature, Specific Conductance, Total Dissolved Solids, pH, and Dissolved Oxygen from across Multiple Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides two-hour time series hydrological and water chemistry sensor data, manual chamber open channel respiration data, handheld sensor water chemistry data, river substrate grain size photos, general environmental context photos, and field metadata (including qualitative information on instream and river corridor characteristics) collected during the same two-week period at 47 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Grain size photos can be used to improve estimates of channel substrate D50 data. Related sample-based water chemistry data are published separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1898914.This dataset is comprised of four main folders, one containing three sensor-specific subfolders and the others containing photographs. The SFA_SpatialStudy_2021_SensorData main data folder includes file-level metadata (FLMD), data dictionary (dd), installation methods, field metadata, Ultrameter water chemistry data, field data collection protocols, international generic sample number (IGSN) mapping file, and a readme file. The “Sensor_Manual_Specifications” subfolder contains pdf files from the manufacturer of each sensor with details on the sensor specifications. Each sensor subfolder (BarotrollAtm, MantaRiver, and MinidotManualChamber) contains a sensor data subfolder for timeseries data and a subfolder for plots and summary statistics. The BarotrollAtm Data subfolder contains In Situ Rugged BaroTROLL pressure and temperature data. The MantaRiver Data subfolder contains Eureka Manta+ 35B multisonde temperature, specific conductance, and pH data. The MinidotManualChamber Data subfolder contains PME MiniDOT Logger dissolved oxygen (mg/L and percent saturation) and temperature data. The folder SFA_SpatialStudy_2021_EnvironmentalContextPhotos contains environmental context photographs and videos. The folders SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part1 and SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part2 contain sediment quadrat photographs. All files are .csv, .pdf, .R, .jpg, .jpeg, .mp4, or .mov. This data package was originally published September 2022. It was updated January 2023 (modified files) and June 2024 (new and modified files). See the change history in data package readme for more details.We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unrolled Video Super-Resolution Network with Autoregressive Prior for the Case of Known Motion

Real-time detection and classification of distant objects is necessary for many national security applications. However, when objects are far from the sensor, they occupy only a small number of pixels in the captured video, limiting the amount of visual detail available for recognition. State-of-the-art classification methods typically rely on high-resolution (HR) video streams to capture characteristic object features, but obtaining such detail is challenging for distant objects that occupy only a few pixels. This motivates the development of video super-resolution (VSR) methods that enhance object classification by recovering fine details from low-pixel representations. Current VSR methods rely either on model-based optimization, which is interpretable but computationally expensive, or on learning-based approaches, which are efficient and high-performing but often lack flexibility and interpretability. In this report, we propose an end-to-end trainable unrolled VSR network, UVSRNet, which super-resolves each frame in a video by exploiting sub-pixel motion between neighboring low-resolution (LR) frames as well as incorporating high-frequency detail from previously super-resolved frames. In particular, by unrolling a plug-and-play (PnP) half-quadratic splitting (HQS) algorithm, we leverage a model-based data-fitting module alongside a learning-based autoregressive prior module. This combination yields a method that maintains the flexibility and interpretability of model-based methods while achieving the performance advantages of learning-based methods.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2017

Stabilized Finite Elements in FUN3D

A Streamlined Upwind Petrov-Galerkin (SUPG) stabilized finite-element discretization has been implemented as a library into the FUN3D unstructured-grid flow solver. Motivation for the selection of this methodology is given, details of the implementation are provided, and the discretization for the interior scheme is verified for linear and quadratic elements by using the method of manufactured solutions. A methodology is also described for capturing shocks, and simulation results are compared to the finite-volume formulation that is currently the primary method employed for routine engineering applications. The finite-element methodology is demonstrated to be more accurate than the finite-volume technology, particularly on tetrahedral meshes where the solutions obtained using the finite-volume scheme can suffer from adverse effects caused by bias in the grid. Although no effort has been made to date to optimize computational efficiency, the finite-element scheme is competitive with the finite-volume scheme in terms of computer time to reach convergence.

Anderson, W. Kyle↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

The dynamics and control of large flexible space structures - 13

The optimal control of three-dimensional large angle maneuvers and vibrations of a Shuttle-mast-reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam subject to three-dimensional deformations. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem is then solved by using the quasilinearization algorithm and the method of particular solutions. The study of the large angle maneuvering of the Shuttle-beam-reflector spacecraft in the plane of a circular earth orbit is extended to consider the effects of the structural offset connection, the axial shortening, and the gravitational torque on the slewing motion. Finally the effect of additional design parameters (such as related to additional payload requirement) on the linear quadratic regulator based design of an orbiting control/structural system is examined.

Bainum, Peter M.↗