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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An investigation of Newton-Sketch and subsampled Newton methods

Sketching, a dimensionality reduction technique, has received much attention in the statistics community. In this paper, we study sketching in the context of Newton's method for solving finite-sum optimization problems in which the number of variables and data points are both large. In this work, we study two forms of sketching that perform dimensionality reduction in data space: Hessian subsampling and randomized Hadamard transformations. Each has its own advantages, and their relative tradeoffs have not been investigated in the optimization literature. Additionally, our study focuses on practical versions of the two methods in which the resulting linear systems of equations are solved approximately, at every iteration, using an iterative solver. The advantages of using the conjugate gradient method vs. a stochastic gradient iteration are revealed through a set of numerical experiments, and a complexity analysis of the Hessian subsampling method is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition

In “Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition,” we review the state of the art in practical algorithms for scheduling power-systems operations in the short term and the results of the recent competition organized by the U.S. Advanced Research Projects Agency–Energy. We explain the mixed-integer nonlinear formulation used in the competition for nonspecialists in electrical engineering, the context and organization of the competition, and the performance of competitors. We find that the collective approaches and results of competitors provide support for efforts to move nonlinear optimization techniques into industrial applications, as they have proven to be a robust and efficient alternative to current linear approximation techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Socially optimal replacement of conventional with electric vehicles for the US household fleet

In this study, a framework is proposed for minimizing the societal cost of replacing gas-powered household passenger cars with battery electric ones (BEVs). The societal cost consists of operational costs of heterogeneous driving patterns' cars, the government investments for charging deployment, and monetized environmental externalities. The optimization framework determines the timeframe needed for conventional vehicles to be replaced with BEVs. It also determines the BEVs driving range during the planning timeframe, as well as the density of public chargers deployed on a linear transportation network over time. We leverage datasets that represent U.S. household driving patterns, as well as the automobile and the energy markets, to apply the model. Results indicate that it takes 8 years for 80% of our conventional vehicle sample to be replaced with electric vehicles, under the base case scenario. The socially optimal all-electric driving range is 204 miles, with chargers placed every 172 miles on a linear corridor. All of the public chargers should be deployed at the beginning of the planning horizon to achieve greater savings over the years. Sensitivity analysis reveals that the timeframe for the socially optimal conversion of 80% of the sample varies from 6 to 12 years. The optimal decision variables are sensitive to battery pack and vehicle body cost, gasoline cost, the discount rate, and conventional vehicles' fuel economy. In conclusion, faster conventional vehicle replacement is achieved when the gasoline cost increases, electricity cost decreases, and battery packs become cheaper over the years.

25 ENERGY STORAGE↗
DOE OSTI2016

An online replanning method using warm start optimization and aperture morphing for flattening-filter-free beams

Purpose: In a situation where a couch shift for patient positioning is not preferred or prohibited (e.g., MR-linac), segment aperture morphing (SAM) can address target dislocation and deformation. For IMRT/VMAT with flattening-filter-free (FFF) beams, however, SAM method would lead to an adverse translational dose effect due to the beam unflattening. Here the authors propose a new two-step process to address both the translational effect of FFF beams and the target deformation. Methods: The replanning method consists of an offline and an online step. The offline step is to create a series of preshifted-plans (PSPs) obtained by a so-called “warm start” optimization (starting optimization from the original plan, rather than from scratch) at a series of isocenter shifts. The PSPs all have the same number of segments with very similar shapes, since the warm start optimization only adjusts the MLC positions instead of regenerating them. In the online step, a new plan is obtained by picking the closest PSP or linearly interpolating the MLC positions and the monitor units of the closest PSPs for the shift determined from the image of the day. This two-step process is completely automated and almost instantaneous (no optimization or dose calculation needed). The previously developed SAM algorithm is then applied for daily deformation. The authors tested the method on sample prostate and pancreas cases. Results: The two-step interpolation method can account for the adverse dose effects from FFF beams, while SAM corrects for the target deformation. Plan interpolation method is effective in diminishing the unflat beam effect and may allow reducing the required number of PSPs. The whole process takes the same time as the previously reported SAM process (5–10 min). Conclusions: The new two-step method plus SAM can address both the translation effects of FFF beams and target deformation, and can be executed in full automation except the delineation of target contour required by the SAM process.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Optimal Energy Management for Microgrid with Stationary and Mobile Storage

This paper studies energy management in a Microgrid (MG) with solar generation, Battery Energy Management System (BESS) and gridable (V2G) EVs. A two-stage stochastic optimization method is proposed to capture the intermittent solar generation and random EV user behaviors. It is subsequently formulated as a Mixed Integer Linear Programming (MILP) problem. To evalutate the proposed method, real solar generation, loads, BESS and EV data is used in Sample Average Approximation (SAA). Computational results show the correctness of the proposed method as well as steady and tightly bounded optimality gap. Comparisons demonstrate that the proposed stochastic method outperforms its deterministic counterpart at the expense of higher computational cost. It is also observed that moderate number of EVs helps to reduce the overall operational cost of the MG, which sheds light on future EV integration to the smart grid.

Battery Energy Management System (BESS)↗
DOE OSTI2015

Trajectory Modulated Arc Therapy: A Fully Dynamic Delivery With Synchronized Couch and Gantry Motion Significantly Improves Dosimetric Indices Correlated With Poor Cosmesis in Accelerated Partial Breast Irradiation

Purpose: To develop planning and delivery capabilities for linear accelerator–based nonisocentric trajectory modulated arc therapy (TMAT) and to evaluate the benefit of TMAT for accelerated partial breast irradiation (APBI) with the patient in prone position. Methods and Materials: An optimization algorithm for volumetrically modulated arc therapy (VMAT) was generalized to allow for user-defined nonisocentric TMAT trajectories combining couch rotations and translations. After optimization, XML scripts were automatically generated to program and subsequently deliver the TMAT plans. For 10 breast patients in the prone position, TMAT and 6-field noncoplanar intensity modulated radiation therapy (IMRT) plans were generated under equivalent objectives and constraints. These plans were compared with regard to whole breast tissue volume receiving more than 100%, 80%, 50%, and 20% of the prescription dose. Results: For TMAT APBI, nonisocentric collision-free horizontal arcs with large angular span (251.5 ± 7.9°) were optimized and delivered with delivery time of ∼4.5 minutes. Percentage changes of whole breast tissue volume receiving more than 100%, 80%, 50%, and 20% of the prescription dose for TMAT relative to IMRT were −10.81% ± 6.91%, −27.81% ± 7.39%, −14.82% ± 9.67%, and 39.40% ± 10.53% (P≤.01). Conclusions: This is a first demonstration of end-to-end planning and delivery implementation of a fully dynamic APBI TMAT. Compared with IMRT, TMAT resulted in marked reduction of the breast tissue volume irradiated at high doses.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Relations between environmental noise and electronic coupling for optimal exciton transfer in one- and two-dimensional homogeneous and inhomogeneous quantum systems

Recent studies have indicated that environmental noise may increase energy-transfer efficiency in quantum systems. For homogeneous networks of chromophores previous studies have primarily considered excitonic transport in one-dimensional (linear) networks. In our study, we expand previous research to a two-dimensional fully coupled topology of chromophore molecules. We demonstrate that not only does an optimal dephasing rate exist in both one- and two-dimensional networks but also that it increases in magnitude with increasing coupling strength between chromophores. Optimal transport occurs when the noise quenches the entanglement between local modes that prevent the exciton from moving efficiently to the target site. We find that these results are insensitive to minor site defects such as those found in realistic systems. We contrast these findings to systems with a high degree of inhomogeneity, in which the optimal dephasing rate is largely set by the system topology and does not vary significantly with respect to coupling strength. Our findings have potential applications to systems such as quantum dot arrays and carbon nanotube structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Decoding and optimized implementation of SECDED codes over GF(q)

A plurality of columns for a check matrix that implements a distance d linear error correcting code are populated by providing a set of vectors from which to populate the columns, and applying to the set of vectors a filter operation that reduces the set by eliminating therefrom all vectors that would, if used to populate the columns, prevent the check matrix from satisfying a column-wise linear independence requirement associated with check matrices of distance d linear codes. One of the vectors from the reduced set may then be selected to populate one of the columns. The filtering and selecting repeats iteratively until either all of the columns are populated or the number of currently unpopulated columns exceeds the number of vectors in the reduced set. Columns for the check matrix may be processed to reduce the amount of logic needed to implement the check matrix in circuit logic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Design, decoding and optimized implementation of SECDED codes over GF(q)

A plurality of columns for a check matrix that implements a distance d linear error correcting code are populated by providing a set of vectors from which to populate the columns, and applying to the set of vectors a filter operation that reduces the set by eliminating therefrom all vectors that would, if used to populate the columns, prevent the check matrix from satisfying a column-wise linear independence requirement associated with check matrices of distance d linear codes. One of the vectors from the reduced set may then be selected to populate one of the columns. The filtering and selecting repeats iteratively until either all of the columns are populated or the number of currently unpopulated columns exceeds the number of vectors in the reduced set. Columns for the check matrix may be processed to reduce the amount of logic needed to implement the check matrix in circuit logic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Decoding and optimized implementation of SECDED codes over GF(q)

A plurality of columns for a check matrix that implements a distance d linear error correcting code are populated by providing a set of vectors from which to populate the columns, and applying to the set of vectors a filter operation that reduces the set by eliminating therefrom all vectors that would, if used to populate the columns, prevent the check matrix from satisfying a column-wise linear independence requirement associated with check matrices of distance d linear codes. One of the vectors from the reduced set may then be selected to populate one of the columns. The filtering and selecting repeats iteratively until either all of the columns are populated or the number of currently unpopulated columns exceeds the number of vectors in the reduced set. Columns for the check matrix may be processed to reduce the amount of logic needed to implement the check matrix in circuit logic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Understanding the Impacts of Biomass Blending on the Uncertainty of Hydrolyzed Sugar Yield from a Stochastic Perspective

Feedstock price and availability are key challenges for biorefinery development and biomass blending has been suggested as a route to overcome these limitations. However, the impacts of feedstock blending on the variability of yields of hydrolyzed sugar remained unclear. This study quantifies the uncertainties in the sugar yields from the blends of corn stover, switchgrass, and grass clippings by considering both variable and model uncertainties. The results indicate that feedstock blending reduces the uncertainties in sugar yields and delivers feedstock of more uniform quality. A further investigation on xxx demonstrates that the reduction in overall uncertainty is achieved by reduced compositional uncertainties of blended feedstock. Based on these results, the optimization of blending ratios was performed with respect to three objectives: (1) probability of meeting the sugar yield target; (2) expected sugar yields; and (3) sugar yield per unit feedstock expense, while satisfying constraints of feedstock availability and cost. In conclusion, the optimization method developed in this study, is readily applied to other combinations of feedstocks, biofuels production processes, and constraints.

09 BIOMASS FUELS↗
DOE OSTI2023

First principles optimization of plutonium electrorefining

Herein this work presents a means of controlling plutonium electrorefining at a maximum rate regardless of equipment setup through the derivation of power supply current and potential governing equations for normal and off-normal operations. The governing equations are demonstrated by electrorefining surrogate materials. A simple linear current sweeping method was used to determine the maximum electrorefining current for the surrogate system. This method can be used to develop autonomous process optimization, real-time online processing monitoring, and real-time process endpoint detection. Ultimately, this research provides the foundation to optimize the liquid metal electrorefining rate to decrease the time needed to the physical limit for the process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Optimal Transfer Operators in Algebraic Two-Level Methods for Nonsymmetric and Indefinite Problems

Consider an algebraic two-level method applied to the 𝑛-dimensional linear system 𝐴⁢𝒙 = 𝒃 using fine-space preconditioner (i.e., “relaxation” or “smoother”) 𝑀, with 𝑀 ≈ 𝐴, restriction and interpolation 𝑅 and 𝑃, and algebraic coarse-space operator 𝐴 𝑐 : = 𝑅 ∗ ⁢𝐴⁢𝑃. Then, what are the best possible transfer operators 𝑅 and 𝑃 of a given dimension 𝑛 𝑐 < 𝑛? Brannick et al. [12] showed that when 𝐴 and 𝑀 are Hermitian positive definite (HPD), the optimal interpolation is such that its range contains the 𝑛 𝑐 smallest generalized eigenvectors of the matrix pencil (𝐴, 𝑀). Recently, in Ali et al. [5] we generalized this framework to the non-HPD setting, by considering both right (interpolation) and left (restriction) generalized eigenvectors of (𝐴, 𝑀) and defining corresponding nonsymmetric transfer operators {𝑅#, 𝑃#}. Tight convergence bounds for {𝑅#, 𝑃#} are derived in spectral radius, as well as a proof of pseudo-optimality. Note, {𝑅#, 𝑃#} are typically complex valued, which is not practical for real-valued problems. Here, in this work, we build on [5], first characterizing all inner products in which the coarse-space correction defined by {𝑅#, 𝑃#} is orthogonal. We then develop tight two-level convergence bounds in these norms, and prove that the underlying transfer operators {𝑅#, 𝑃#} are genuinely optimal. As a special case, our theory both recovers and extends the HPD results from [12]. Finally, we show how to construct optimal, real-valued transfer operators in the case of that 𝐴 and 𝑀 are real valued, but are not HPD. Numerical examples arising from a discretized advection-reaction equation, wave-equation, and Stokes equations are used to verify and illustrate the theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

TH-C-17A-07: Visualizing and Quantifying Radiation Therapy in Real-Time Using a Novel Beam Imaging Technique

Purpose: Radiation therapy uses invisible high energy X-rays to treat an invisible tumor. Proper positioning of the treatment beam relative to the patient's anatomy during dose delivery is critically important to the success of treatment. We develop and characterize a novel radiation therapy beam visualization technique for real-time monitoring of patient treatment. Methods: Custom made flexible scintillator sheets were fabricated from gadolinium oxysulfide (GOS) particles that had been doped with terbium within a silicone elastomer matrix. Sheets of several thicknesses ranging from 0.3 to 1mm were prepared and tested. Sheets were exposed to megavoltage X-ray and electron beams from a Varian linac and the resulting optical signal was collected by multiple CMOS cameras placed in the treatment room. Real-time images were collected for different beam energies and dose rates. Signal intensity and SNR were calculated by processing the acquired images. Results: All signals were detectable in the presence of full room lighting and at an integration time of 45ms. Average signal intensity and SNR increased with both sheet thickness and dose rate and decreased with beam energy and incident light. For a given sheet thickness and beam energy the correlation between dose rate and signal intensity was highly linear. Increased sheet thickness or dose rate results in a linear increase in the detected signal. All results are consistent with analytical approximations. Conclusion: The technique offers a means of accurately visualizing a radiation therapy beam shape and fluence in real time. The effects of salient parameters have been characterized and will enable further optimization of the technique as it is implemented into the clinical workflow. The project described was supported by the National Center for Advancing Translational Sciences of the National Institutes of Health through UL1 TR001085.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

The value of continuity: Refined isogeometric analysis and fast direct solvers

Here, we propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers. Starting from a highly continuous Isogeometric Analysis (IGA) discretization, we introduce C0-separators to reduce the interconnection between degrees of freedom in the mesh. By doing so, both the solution time and best approximation errors are simultaneously improved. We call the resulting method “refined Isogeometric Analysis (rIGA)”. To illustrate the impact of the continuity reduction, we analyze the number of Floating Point Operations (FLOPs), computational times, and memory required to solve the linear system obtained by discretizing the Laplace problem with structured meshes and uniform polynomial orders. Theoretical estimates demonstrate that an optimal continuity reduction may decrease the total computational time by a factor between p 2 and p 3 , with pp being the polynomial order of the discretization. Numerical results indicate that our proposed refined isogeometric analysis delivers a speed-up factor proportional to p 2 . In a 2D mesh with four million elements and p=5, the linear system resulting from rIGA is solved 22 times faster than the one from highly continuous IGA. In a 3D mesh with one million elements and p=3, the linear system is solved 15 times faster for the refined than the maximum continuity isogeometric analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Inverse compton light source: a compact design proposal

In the last decade, there has been an increasing demand for a compact Inverse Compton Light Source (ICLS) which is capable of producing high-quality X-rays by colliding an electron beam and a high-quality laser. It is only in recent years when both SRF and laser technology have advanced enough that compact sources can approach the quality found at large installations such as the Advanced Photon Source at Argonne National Laboratory. Previously, X-ray sources were either high flux and brilliance at a large facility or many orders of magnitude lesser when produced by a bremsstrahlung source. A recent compact source was constructed by Lyncean Technologies using a storage ring to produce the electron beam used to scatter the incident laser beam. By instead using a linear accelerator system for the electron beam, a significant increase in X-ray beam quality is possible, though even subsequent designs also featuring a storage ring offer improvement. Preceding the linear accelerator with an SRF reentrant gun allows for an extremely small transverse emittance, increasing the brilliance of the resulting X-ray source. In order to achieve sufficiently small emittances, optimization was done regarding both the geometry of the gun and the initial electron bunch distribution produced off the cathode. Using double-spoke SRF cavities to comprise the linear accelerator allows for an electron beam of reasonable size to be focused at the interaction point, while preserving the low emittance that was generated by the gun. An aggressive final focusing section following the electron beam's exit from the accelerator produces the small spot size at the interaction point which results in an X-ray beam of high flux and brilliance. Taking all of these advancements together, a world class compact X-ray source has been designed. It is anticipated that this source would far outperform the conventional bremsstrahlung and many other compact ICLSs, while coming closer to performing at the levels found at large facilities than ever before. The design process, including the development between subsequent iterations, is presented here in detail, with the simulation results for this groundbreaking X-ray source.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Thermal effects on geologic carbon storage

One of the most promising ways to significantly reduce greenhouse gases emissions, while carbon-free energy sources are developed, is Carbon Capture and Storage (CCS). Non-isothermal effects play a major role in all stages of CCS. In this paper, we review the literature on thermal effects related to CCS, which is receiving an increasing interest as a result of the awareness that the comprehension of non-isothermal processes is crucial for a successful deployment of CCS projects. We start by reviewing CO 2 transport, which connects the regions where CO 2 is captured with suitable geostorage sites. The optimal conditions for CO 2 transport, both onshore (through pipelines) and offshore (through pipelines or ships), are such that CO 2 stays in liquid state. To minimize costs, CO 2 should ideally be injected at the wellhead in similar pressure and temperature conditions as it is delivered by transport. To optimize the injection conditions, coupled wellbore and reservoir simulators that solve the strongly non-linear problem of CO 2 pressure, temperature and density within the wellbore and non-isothermal two-phase flow within the storage formation have been developed. CO 2 in its way down the injection well heats up due to compression and friction at a lower rate than the geothermal gradient, and thus, reaches the storage formation at a lower temperature than that of the rock. Inside the storage formation, CO 2 injection induces temperature changes due to the advection of the cool injected CO 2 , the Joule-Thomson cooling effect, endothermic water vaporization and exothermic CO 2 dissolution. These thermal effects lead to thermo-hydro-mechanical-chemical coupled processes with non-trivial interpretations. These coupled processes also play a relevant role in “Utilization” options that may provide an added value to the injected CO 2 , such as Enhanced Oil Recovery (EOR), Enhanced Coal Bed Methane (ECBM) and geothermal energy extraction combined with CO 2 storage. If the injected CO 2 leaks through faults, the caprock or wellbores, strong cooling will occur due to the expansion of CO 2 as pressure decreases with depth. Finally, we conclude by identifying research gaps and challenges of thermal effects related to CCS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A neural-network-enhanced parameter-varying framework for multi-objective model predictive control applied to buildings

Management of the electrical grid is becoming more complex due to the increased penetration of alternative energy generation technologies and a broadening diversity of electric loads. This complexity creates challenges in balancing demand and generation that can increase the potential for grid instabilities. One effective way to address this issue is to leverage previously unexploited demand flexibility through advanced control strategies. In this work, we propose an advanced control method, called adaptive neural parameter-varying model predictive control (ANPV-MPC), to control the temperature and energy consumption of a building via its Heating, Ventilation, and Air Conditioning system. ANPV-MPC combines key ideas in parameter-varying control, adaptive control, and online learning strategies to bridge the gap between computationally efficient linear model predictive control and more accurate nonlinear model predictive control. The novelty in ANPV-MPC is the use of a physics-inspired Bayesian neural network to estimate the coefficients of the parameter-varying linear control model. The Bayesian neural network additionally provides uncertainty estimates, triggering online training to capture evolving building system conditions. We show that ANPV-MPC can approximate the building system dynamics with a 28.39% higher accuracy than traditional linear model predictive control, resulting in 36.23% better control performance without increasing complexity of the optimal control problem. ANPV-MPC also adapts in real time to previously unseen conditions using online learning, further improving its performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗