Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Results of an integrated structure-control law design sensitivity analysis

Next generation air and space vehicle designs are driven by increased performance requirements, demanding a high level of design integration between traditionally separate design disciplines. Interdisciplinary analysis capabilities have been developed, for aeroservoelastic aircraft and large flexible spacecraft control for instance, but the requisite integrated design methods are only beginning to be developed. One integrated design method which has received attention is based on hierarchal problem decompositions, optimization, and design sensitivity analyses. This paper highlights a design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changess in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient that finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗
DOE OSTI2014

SU-E-T-444: Quantitatively Comparison of Low Dose Spillage Outside of PTV Edge in Arc Therapy Modalities

Purpose: To quantitatively compare low dose spillage outside of PTV edge in arc therapy modalities Methods: The machines used in the study are Tomotherapy Hi-Arc and Varian 21EX with millennium120 MLC. TPS are TomoPlaning and RayStation for VMAT, respectively. The phantom is a 30cm diameter cylindrical solid water (TOMOTHERAPY, TOMOPHANTOM ASSY). The PTV is 4cm length with ellipsoidal sectional shape with major axis=5cm, minor axis=3cm in the axial plane and reversed in the coronal plane. The PTV volume is created with interpolation. It is located at the center of the phantom. The prescribed dose is 1000x5 cGy to 95% the PTV. The isocenter is set co-centered with the PTV. EBT-3 film was used to measure iso-dose lines at the center plane. Film dosimetry is performed with the RIT, v6.2. Results: the study shows: (1) dose falloff gradient is usually uneven, depending on the PTV shape in the gantry rotation plane. For an elliptical shape, the low dose spillage is wider in the minor axis direction than that in the major axis direction. The more a shape is closer to circular, the more even gradient is all directions; (2)for a circular shape (CAX plane in this study), the maximum dose in % of Rx dose at 2cm from PTV is 55% for Tomo, vs. 70% for VMAT (3) the most rapid dose falloff rang is between 95%–80% IDL for both modalities. Conclusion: Tomo has more rapid dose falloff outside of PTV. In some areas, the gradient is double for Tomo helical than that for LINAC VMAT at same points. Future work will examine the differences between optimization of doses and inherent delivery limitations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗
DOE OSTI2016

The impact of short-term stochastic variability in solar irradiance on optimal microgrid design

This paper proposes a new methodology to capture the impact of fast moving clouds on utility power demand charges observed in microgrids with photovoltaic (PV) arrays, generators, and electrochemical energy storage. It consists of a statistical approach to introduce sub-hourly events in the hourly economic accounting process. The methodology is implemented in the Distributed Energy Resources Customer Adoption Model (DER-CAM), a state of the art mixed integer linear model used to optimally size DER in decentralized energy systems. Results suggest that previous iterations of DER-CAM could undersize battery capacities. The improved model depicts more accurately the economic value of PV as well as the synergistic benefits of pairing PV with storage.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Domain Knowledge Guided Bayesian Optimization For Autonomous Alignment Of Complex Scientific Instruments

Bayesian Optimization (BO) is a powerful tool for optimizing complex non-linear systems. However, its performance degrades in high-dimensional problems with tightly coupled parameters and highly asymmetric objective landscapes, where rewards are sparse. In such needle-in-a-haystack scenarios, even advanced methods like trust-region BO (TurBO) often lead to unsatisfactory results. We propose a domain knowledge guided Bayesian Optimization approach, which leverages physical insight to fundamentally simplify the search problem by transforming coordinates to decouple input features and align the active subspaces with the primary search axes. We demonstrate this approach's efficacy on a challenging 12-dimensional, 6-crystal Split-and-Delay optical system, where conventional approaches, including standard BO, TuRBO and multi-objective BO, consistently led to unsatisfactory results. When combined with an reverse annealing exploration strategy, this approach reliably converges to the global optimum. The coordinate transformation itself is the key to this success, significantly accelerating the search by aligning input co-ordinate axes with the problem's active subspaces. As increasingly complex scientific instruments, from large telescopes to new spectrometers at X-ray Free Electron Lasers are deployed, the demand for robust high-dimensional optimization grows. Our results demonstrate a generalizable paradigm: leveraging physical insight to transform high-dimensional, coupled optimization problems into simpler representations can enable rapid and robust automated tuning for consistent high performance while still retaining current optimization algorithms.

FOS: Computer and information sciences↗
DOE OSTI2024

Control and Optimization of Energy Storage System in Power Distribution System

The widespread adoption of electric vehicles (EVs) and transportation electrification is encumbered by two chief barriers: i) the limited driving range of EVs in the market today and ii) inadequate fast-charging infrastructure for long-distance trips. Extreme fast charging (XFC) technology can recharge EVs in less than 10 minutes for 200 miles range. Firstly, a novel robust optimization-based mixed integer linear programming model is proposed to size a battery energy storage system (BESS) and PV system in an XFCS. In this part, it is assumed that the sizing and location of the XFCS are known. Secondly, the aforesaid assumption is relaxed, and a strategic multi-period coordinated planning model is proposed to optimally site and size BESS-assisted charging stations in a highway transportation network and PV systems in a power distribution network by considering the coupling between both networks. Optimal operation and control of BESS-assisted EV charging stations are vital to alleviate the adverse impact of extreme fast charging of EVs on the host power network. A joint solution is proposed to mitigate the steady state and transient impact of extremefast charging of EVs and ensure grid-friendly integration of XFCSs with the host grid. Lastly, to make the operation of the XFCS cost-effective, a multi-layered energy management framework is proposed for the XFCS by considering forecast uncertainties, monthly demand charges reduction, and BESS degradation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1975

A linear feedback guidance law for an Aeromaneuvering Orbit-to-Orbit Shuttle

The sensitivity of the trajectory of an Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) to unpredictable variations in atmospheric density is shown and a compensating guidance law is developed. The sensitivity was determined using a varying atmosphere model whose statistical properties matched those of actual atmosphere data. This study showed that a compensating guidance law is crucial to the feasibility of AMOOS. Using an optimal regulator approach, a linear feedback law for bank angle was developed and checked out with the atmosphere model. This law uses the altitude profile of the trajectory to control the amount of energy dissipated within the atmosphere.

Rehder, J. J.↗
DOE OSTI2016

SU-F-T-344: Commissioning Constant Dose Rate VMAT in the Raystation Treatment Planning System for a Varian Clinac IX

Purpose: To demonstrate the commissioning of constant dose rate volumetric modulated arc therapy (VMAT) in the Raystation treatment planning system for a Varian Clinac iX with Exact couch. Methods: Constant dose rate (CDR) VMAT is an option in the Raystation treatment planning system, enabling VMAT delivery on Varian linacs without a RapidArc upgrade. Raystation 4.7 was used to commission CDR-VMAT for a Varian Clinac iX. Raystation arc model parameters were selected to match machine deliverability characteristics. A Varian Exact couch model was added to Raystation 4.7 and commissioned for use in VMAT optimization. CDR-VMAT commissioning checks were performed on the linac, including patient-specific QA measurements for 10 test patients using both the ArcCHECK from Sun Nuclear Corporation and COMPASS from IBA Dosimetry. Multi-criteria optimization (MCO) in Raystation was used for CDR-VMAT planning. Results: Raystation 4.7 generated clinically acceptable and deliverable CDR-VMAT plans for the Varian Clinac. VMAT plans were optimized including a model of the Exact couch with both rails in the out positions. CDR-VMAT plans generated with MCO in Raystation were dosimetrically comparable to Raystation MCO-generated IMRT plans. Patient-specific QA measurements with the ArcCHECK on the couch showed good agreement with the treatment planning system prediction. Patient-specific, structure-specific, multi-statistical parameter 3D QA measurements with gantry-mounted COMPASS also showed good agreement. Conclusion: Constant dose rate VMAT was successfully modeled in Raystation 4.7 for a Varian Clinac iX, and Raystation’s multicriteria optimization generated constant dose rate VMAT plans which were deliverable and dosimetrically comparable to IMRT plans.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1978

Specific-optimal control with a dual minimal-order observer-based compensator

The steady-state optimal control of a linear time-invariant stochastic system by means of a minimal-order dual-observer-based compensator is considered in this paper. The structure of the compensator is fixed while the associated gains are to be chosen so as to minimize a quadratic penalty on the plant state. Necessary and sufficient conditions for optimality are given, and an explicit solution is displayed. Salient features pertaining to the optimal system are: a decoupling property, a projection property, and an innovation property. Finally, it is shown that this design corresponds to a singular LQG problem, which is precisely the dual of another singular LQG problem: namely Newmann's problem. A complete picture is then given showing clearly the correspondence between the two designs.

Blanvillain, P. J.↗
DOE OSTI2016

SU-F-T-497: Spatiotemporally Optimal, Personalized Prescription Scheme for Glioblastoma Patients Using the Proliferation and Invasion Glioma Model

Purpose: To investigate a spatiotemporally optimal radiotherapy prescription scheme and its potential benefit for glioblastoma (GBM) patients using the proliferation and invasion (PI) glioma model. Methods: Standard prescription for GBM was assumed to deliver 46Gy in 23 fractions to GTV1+2cm margin and additional 14Gy in 7 fractions to GTV2+2cm margin. We simulated the tumor proliferation and invasion in 2D according to the PI glioma model with a moving velocity of 0.029(slow-move), 0.079(average-move), and 0.13(fast-move) mm/day for GTV2 with a radius of 1 and 2cm. For each tumor, the margin around GTV1 and GTV2 was varied to 0–6 cm and 1–3 cm respectively. Total dose to GTV1 was constrained such that the equivalent uniform dose (EUD) to normal brain equals EUD with the standard prescription. A non-stationary dose policy, where the fractional dose varies, was investigated to estimate the temporal effect of the radiation dose. The efficacy of an optimal prescription scheme was evaluated by tumor cell-surviving fraction (SF), EUD, and the expected survival time. Results: Optimal prescription for the slow-move tumors was to use 3.0(small)-3.5(large) cm margins to GTV1, and 1.5cm margin to GTV2. For the average- and fast-move tumors, it was optimal to use 6.0cm margin for GTV1 suggesting that whole brain therapy is optimal, and then 1.5cm (average-move) and 1.5–3.0cm (fast-move, small-large) margins for GTV2. It was optimal to deliver the boost sequentially using a linearly decreasing fractional dose for all tumors. Optimal prescription led to 0.001–0.465% of the tumor SF resulted from using the standard prescription, and increased tumor EUD by 25.3–49.3% and the estimated survival time by 7.6–22.2 months. Conclusion: It is feasible to optimize a prescription scheme depending on the individual tumor characteristics. A personalized prescription scheme could potentially increase tumor EUD and the expected survival time significantly without increasing EUD to normal brain.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Optimization of the artificial viscosity in Lagrangian staggered discretization codes. Modeling 1D stand-alone shock - case study

We have developed new measures of errors for numerical shock. The new approach is based on analysis of the structure function, and separation of the errors related to oscillations and shock width, which also include error in the position of the ”center” of the numerical shock. We have demonstrated that those measures correctly characterize the numerical solution. We introduced an objective function in, which both types of errors are weighted, and presented optimal values of the coefficients of the linear and quadratic viscosity for different weights and different Mach numbers.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Temporal Gain Correction for X-Ray Calorimeter Spectrometers

Calorimetric X-ray detectors are very sensitive to their environment. The boundary conditions can have a profound effect on the gain including heat sink temperature, the local radiation temperature, bias, and the temperature of the readout electronics. Any variation in the boundary conditions can cause temporal variations in the gain of the detector and compromise both the energy scale and the resolving power of the spectrometer. Most production X-ray calorimeter spectrometers, both on the ground and in space, have some means of tracking the gain as a function of time, often using a calibration spectral line. For small gain changes, a linear stretch correction is often sufficient. However, the detectors are intrinsically non-linear and often the event analysis, i.e., shaping, optimal filters etc., add additional non-linearity. Thus for large gain variations or when the best possible precision is required, a linear stretch correction is not sufficient. Here, we discuss a new correction technique based on non-linear interpolation of the energy-scale functions. Using Astro-HSXS calibration data, we demonstrate that the correction can recover the X-ray energy to better than 1 part in 104 over the entire spectral band to above 12 keV even for large-scale gain variations. This method will be used to correct any temporal drift of the on-orbit per-pixel gain using on-board calibration sources for the SXS instrument on the Astro-H observatory.

X-ray detector X-ray spectrometer X-ray calorimete↗
NASA NTRS1990

Optimal fabrication processes for unidirectional metal-matrix composites: A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with non-linear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗
NASA NTRSDate not supplied

Hardware-Efficient Quantum Optimization Layered Algorithms and Experiments

Quantum optimization algorithms, such as QAOA, that implement parametrized stochastic optimization solvers attempt to identify low-energy solutions of Ising systems by exploiting available quantum effects in noisy-intermediate scale machines. Engineering a well-performing parametrized quantum optimization circuit is indeed an exercise in balancing the trade-off between expressivity and implementation complexity. We show that, for MaxCut QAOA circuits defined on native hardware topology (Rigetti’s Aspen Quantum Processors), error-mitigation techniques recover simulated features of the noiseless theory. Moreover, we explore a design space for QAOA-like ansatze that perform well in theory as well as in hardware for fully-connected problems. We also discuss how efficient coherence and entanglement detection methods that could be coupled with quantum optimization experiments require only linear overhead in benchmarking time.

quantum computing↗
DOE OSTI2015

Inconsistent Investment and Consumption Problems

In a traditional Black–Scholes market we develop a verification theorem for a general class of investment and consumption problems where the standard dynamic programming principle does not hold. The theorem is an extension of the standard Hamilton–Jacobi–Bellman equation in the form of a system of non-linear differential equations. We derive the optimal investment and consumption strategy for a mean-variance investor without pre-commitment endowed with labor income. In the case of constant risk aversion it turns out that the optimal amount of money to invest in stocks is independent of wealth. The optimal consumption strategy is given as a deterministic bang-bang strategy. In order to have a more realistic model we allow the risk aversion to be time and state dependent. Of special interest is the case were the risk aversion is inversely proportional to present wealth plus the financial value of future labor income net of consumption. Using the verification theorem we give a detailed analysis of this problem. It turns out that the optimal amount of money to invest in stocks is given by a linear function of wealth plus the financial value of future labor income net of consumption. The optimal consumption strategy is again given as a deterministic bang-bang strategy. We also calculate, for a general time and state dependent risk aversion function, the optimal investment and consumption strategy for a mean-standard deviation investor without pre-commitment. In that case, it turns out that it is optimal to take no risk at all.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

TH-CD-BRA-09: Towards Absolute Dose Measurement in MRI-Linac and Gamma-Knife: Design and Construction of An MR-Compatible Water Calorimeter

Purpose: The purpose of this work is to design and optimize a portable Water Calorimeter (WC) for use in a commercial MRI-linac and Gamma-knife in addition to conventional radiotherapy linacs. Water calorimeters determine absorbed dose to water at a point by measuring radiation-induced temperature rise of the volume (the two are related by the medium specific heat capacity). In this formalism, one important correction factor is heat transfer correction k-ht. It compensates for heat gain/loss due to conductive and convective effects, and is numerically calculated as ratio of temperature rise in the absence of heat loss to that in the presence of heat loss. Operating at 4°C ensures convection is minimal. Methods: A commercial finite element software was used to evaluate several WC designs with different insulation materials and thicknesses; channels allowing coolant to travel around WC (to sustain WC at 4°C) were modeled, and worst-case scenario variation in the temperature of the coolant was simulated for optimization purposes (2.6 mK/s). Additionally, several calorimeter vessel design parameters (front/back glass thickness/separation, diameter) were also simulated and optimized. Optimization is based on minimizing long term calorimeter drift (24h) as well as variation and magnitude of k-ht. Results: The final selected WC design reached a modest drift of 11µK/s after 15h for the worst-case coolant temperature variation. This design consists of coolant channels being encompassed on both sides by cryogel insulation. For the MRI-linac beam, glass thickness plays the largest effect on k-ht with variation of upto 0.6% in the first run for thicknesses ranging between 0.5–1.7mm. Subsequent runs vary only within 0.1% with glass thickness. Other factors such as vessel radius and top/bottom glass separation have sub 0.1% effects on k-ht. Conclusion: An MR-safe 4°C stagnant WC appropriate for dosimetry in MRI-linac and Gamma-Knife was designed, optimized, and construction is nearly completed. NSERC Discovery Grant RGPIN-435608.

07 ISOTOPE AND RADIATION SOURCES↗