Search NASA⌕ Search

SEARCH · Search NASA

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 811 records · Page 45

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗

Tools for analysis of optimization models

The software is code for analyzing, debugging, and solving optimization models. The code implements several algorithms and provides convenient APIs to apply these algorithms to optimization models. The algorithms are primarily based on graph theory. They compute well-known partitions of graphs, and use these partitions to provide the user information about their optimization model, including diagnosing certain types of structural modeling errors. Some algorithms may be used as a subroutine to solve an optimization problem, and may call an optimization solver to facilitate this. Well-known linear algebra routines are called to provide further diagnostics. The software contains no data and no models other than toy models used for testing.

Parker, Robert↗

Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. Key to the efficiency of the method is a technique for fast computation of the product of the blob-blob mobility matrix and a vector. For unbounded suspensions, we rely on existing analytical expressions for the Rotne-Prager-Yamakawa tensor combined with a fast multipole method (FMM) to obtain linear scaling in the number of particles. For suspensions sedimented against a single no-slip boundary, we use a direct summation on a graphical processing unit (GPU), which gives quadratic asymptotic scaling with the number of particles. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently developed rigid-body immersed boundary method by B. Kallemov, A. P. S. Bhalla, B. E. Griffith, and A. Donev (Commun. Appl. Math. Comput. Sci. 11 (2016), no. 1, 79-141) to suspensions of freely moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid-body equations converges in a bounded number of iterations regardless of the system size. In our approach, each iteration only requires a few cycles of a geometric multigrid solver for the Poisson equation, and an application of the block-diagonal preconditioner, leading to linear scaling with the number of particles. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Embedded Sensors and Controls to Improve Component Performance and Reliability - System Dynamics Modeling and Control System Design

This report documents the current status of the modeling, control design, and embedded control research for the magnetic bearing canned rotor pump being used as a demonstration platform for deeply integrating instrumentation and controls (I&C) into nuclear power plant components. This pump is a highly inter-connected thermo/electro/mechanical system that requires an active control system to operate. Magnetic bearings are inherently unstable system and without active, moment by moment control, the rotor would contact fixed surfaces in the pump causing physical damage. This report details the modeling of the pump rotordynamics, fluid forces, electromagnetic properties of the protective cans, active magnetic bearings, power electronics, and interactions between different dynamical models. The system stability of the unforced and controlled rotor are investigated analytically. Additionally, controllers are designed using proportional derivative (PD) control, proportional integral derivative (PID) control, voltage control, and linear quadratic regulator (LQR) control. Finally, a design optimization problem that joins the electrical, mechanical, magnetic, and control system design into one problem to balance the opposing needs of various design criteria using the embedded system approach is presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nalu's Linear System Assembly using Tpetra

The Nalu Exascale Wind application assembles linear systems using data structures provided by the Tpetra package in Trilinos. This note describes the initialization and assembly process. The purpose of this note is to help Nalu developers and maintainers to understand the code surrounding linear system assembly, in order to facilitate debugging, optimizations, and maintenance.

17 WIND ENERGY↗

A Tutorial for Using an Open-Source Solver for the Regional Energy Deployment System (ReEDS) Model

ReEDS is a publicly available model developed at NREL that can be used to analyze the potential evolution of the U.S. electric power system into the future. ReEDS is formulated as a linear program, written in GAMS, and solved using a linear programming solver (e.g., CPLEX). This tutorial provides context and understanding for the implications of using an open-source solver for ReEDS.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

CRYSTAL CHEMISTRY OF THREE-COMPONENT WHITE DWARFS AND NEUTRON STAR CRUSTS: PHASE STABILITY, PHASE STRATIFICATION, AND PHYSICAL PROPERTIES

A systematic search for multicomponent crystal structures is carried out for five different ternary systems of nuclei in a polarizable background of electrons, representative of accreted neutron star crusts and some white dwarfs. Candidate structures are “bred” by a genetic algorithm and optimized at constant pressure under the assumption of linear response (Thomas–Fermi) charge screening. Subsequent phase equilibria calculations reveal eight distinct crystal structures in the T = 0 bulk phase diagrams, five of which are complicated multinary structures not previously predicted in the context of compact object astrophysics. Frequent instances of geometrically similar but compositionally distinct phases give insight into structural preferences of systems with pairwise Yukawa interactions, including and extending to the regime of low-density colloidal suspensions made in a laboratory. As an application of these main results, we self-consistently couple the phase stability problem to the equations for a self-gravitating, hydrostatically stable white dwarf, with fixed overall composition. To our knowledge, this is the first attempt to incorporate complex multinary phases into the equilibrium phase-layering diagram and mass–radius-composition dependence, both of which are reported for He–C–O and C–O–Ne white dwarfs. Finite thickness interfacial phases (“interphases”) show up at the boundaries between single-component body-centered cubic (bcc) crystalline regions, some of which have lower lattice symmetry than cubic. A second application—quasi-static settling of heavy nuclei in white dwarfs—builds on our equilibrium phase-layering method. Tests of this nonequilibrium method reveal extra phases that play the role of transient host phases for the settling species.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating Nationwide Supply Chain for Circularity of PET and Olefin Plastics

PET (Polyethylene Terephthalate, #1) and olefin plastics including HDPE (High Density Polyethylene, #2), LDPE/LLDPE (Low Density/Linear Low-Density Polyethylene, #4) and PP (Polypropylene, #5) together comprise nearly 80% of the U.S. plastic market. Due to their high market share and extensive application in packaging these polymers have better potential for circularity than other polymer types. To understand the potential of future scenarios with higher recycling rates, supply chain scenarios for PET and olefin plastic packaging need to be analyzed with increased availability and collection of plastics. Currently there exists a knowledge gap to understand the circular supply chain on a national level. Performing a nationwide analysis introduces certain challenges such as variability in the plastic mix of the recycling stream in different municipalities, lack of collection of some types of plastics, regional differences in cost and estimating a national supply curve based on state level access rate and participation rate. Currently, a limited number of Plastic Reclaimers recycle the nation’s collected plastic involving transportation over long distances surpassing state boundaries. Modeling a nationwide scenario instead of regional/state based scenario will facilitate transportation beyond state boundaries for the development of a circular economy. A Mixed-integer Linear Programming model was developed, to identify the optimal location and capacity of the Material Recovery Facilities (MRFs) nationwide subject to maximizing the profit margin of the industrial entities within the model. We considered three different scenarios including the base case scenario (S1) collecting 2.27 million metric tons/year as well as two additional scenarios with available plastic supply of 1.8 times (S2) and 2.2 times (S3) of the base case scenario. The model identified 166, 275 and 319 counties as potential MRF locations for scenarios S1, S2 and S3 respectively. Compared to the existing number of counties having MRFs in the US, the model results indicated a reduction of 38% of MRFs for S1. For S2 and S3, the number of counties with MRFs increased 3% and 20% respectively compared to counties currently with MRFs. The average profit remains constant between $180-$181/ton regardless of the increased plastic collection.

54 - ENVIRONMENTAL SCIENCES/GLOBAL CLIMATE CHANGE ↗

Online Primal-Dual Methods With Measurement Feedback for Time-Varying Convex Optimization

This paper addresses the design and analysis of feedback-based online algorithms to control systems or networked systems based on performance objectives and engineering constraints that may evolve over time. The emerging timevarying convex optimization formalism is leveraged to model optimal operational trajectories of the systems, as well as explicit local and network-level operational constraints. Departing from existing batch and feed-forward optimization approaches, the design of the algorithms capitalizes on an online implementation of primal-dual projected-gradient methods; the gradient steps are, however, suitably modified to accommodate feedback from the system in the form of measurements - hence, the term 'online optimization with feedback.' By virtue of this approach, the resultant algorithms can cope with model mismatches in the algebraic representation of the system states and outputs, they avoid pervasive measurements of exogenous inputs, and they naturally lend themselves to a distributed implementation. Under suitable assumptions, analytical convergence claims are established in terms of dynamic regret. Furthermore, when the synthesis of the feedback-based online algorithms is based on a regularized Lagrangian function, Q-linear convergence to solutions of the time-varying optimization problem is shown.

97 MATHEMATICS AND COMPUTING↗

A Linear Programming Approach to Backtracking for Single-Axis Trackers on Rolling Terrain

In this article, we present a computationally efficient method for determining optimal backtracking rotations for single-axis solar trackers on nonuniform terrain. The method allows for ganged tracking, mechanical rotation constraints, uneven row spacing, and arbitrary maximum allowable shaded fractions (to enable “fractional backtracking”). As with previous 2-D approaches, the method is suitable for terrain that varies in the transverse direction with respect to the rotation axis of the trackers. The novelty of the method lies in formulating the problem of shade avoidance as a linear problem, which is achieved by using the row interception width as the optimization variable instead of rotation angles. Formulating backtracking as a linear problem enables the use of extremely efficient linear programming algorithms, making the method highly scalable, requiring less than 1 min to compute optimal rotation schedules for hundreds of trackers. It also produces more effective backtracking rotations, reducing the frequency of shading by 4× and improving system energy output by 1%–2%.

Optimization↗

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

SU-E-T-197: Helical Cranial-Spinal Treatments with a Linear Accelerator

Purpose: Craniospinal irradiation (CSI) of systemic disease requires a high level of beam intensity modulation to reduce dose to bone marrow and other critical structures. Current helical delivery machines can take 30 minutes or more of beam-on time to complete these treatments. This pilot study aims to test the feasibility of performing helical treatments with a conventional linear accelerator using longitudinal couch travel during multiple gantry revolutions. Methods: The VMAT optimization package of the Eclipse 10.0 treatment planning system was used to optimize pseudo-helical CSI plans of 5 clinical patient scans. Each gantry revolution was divided into three 120° arcs with each isocenter shifted longitudinally. Treatments requiring more than the maximum 10 arcs used multiple plans with each plan after the first being optimized including the dose of the others (Figure 1). The beam pitch was varied between 0.2 and 0.9 (couch speed 5- 20cm/revolution and field width of 22cm) and dose-volume histograms of critical organs were compared to tomotherapy plans. Results: Viable pseudo-helical plans were achieved using Eclipse. Decreasing the pitch from 0.9 to 0.2 lowered the maximum lens dose by 40%, the mean bone marrow dose by 2.1% and the maximum esophagus dose by 17.5%. (Figure 2). Linac-based helical plans showed dose results comparable to tomotherapy delivery for both target coverage and critical organ sparing, with the D50 of bone marrow and esophagus respectively 12% and 31% lower in the helical linear accelerator plan (Figure 3). Total mean beam-on time for the linear accelerator plan was 8.3 minutes, 54% faster than the tomotherapy average for the same plans. Conclusions: This pilot study has demonstrated the feasibility of planning pseudo-helical treatments for CSI targets using a conventional linac and dynamic couch movement, and supports the ongoing development of true helical optimization and delivery.

60 APPLIED LIFE SCIENCES↗

Maximizing dynamic range and performance of anatase TiO 2 ECRAM through structure and programming

Here, in this study, we investigate the structure-dependent modulation characteristics of all-solid-state three-terminal electrochemical random-access memory (ECRAM) based on an anatase Li x TiO 2 channel. By directly comparing “asymmetric” and “symmetric” ECRAM device architectures, we reveal significant insight into the impact of a non-zero gate-drain open-circuit voltage and its influence on voltage vs. current-controlled gating. We also explore the impact of potentiation/depression write parameters on the symmetry, linearity, and dynamic range of the device response. Together, initial results from optimizing structure and programming approaches yielded unprecedented G max /G min ratios of >1,000 for ECRAM and hundreds of tunable memory states with excellent linearity and symmetry. Simulations based on these ECRAM devices further illustrate the promise of this analog memory technology, achieving near 2% classification error in the MNIST digit recognition benchmark for a range of training parameters compared to a theoretical best of 1.66% and outperforming other device models extracted from the literature.

AIHWKit↗

Multi-output multilevel best linear unbiased estimators via semidefinite programming

Multifidelity forward uncertainty quantification (UQ) problems often involve multiple quantities of interest and heterogeneous models (e.g., different grids, equations, dimensions, physics, surrogate and reduced-order models). While computational efficiency is key in this context, multi-output strategies in multilevel/multifidelity methods are either sub-optimal or non-existent. In this paper we extend multilevel best linear unbiased estimators (MLBLUE) to multi-output forward UQ problems and we present new semidefinite programming formulations for their optimal setup. Not only do these formulations yield the optimal number of samples required, but also the optimal selection of low-fidelity models to use. While existing MLBLUE approaches are single-output only and require a non-trivial nonlinear optimization procedure, the new multi-output formulations can be solved reliably and efficiently. Here, we demonstrate the efficacy of the new methods and formulations in practical UQ problems with model heterogeneity.

97 MATHEMATICS AND COMPUTING↗

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants: Model-Predictive Control Library Development (Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.

01 COAL, LIGNITE, AND PEAT↗

An optimal guarding scheme for thermal conductivity measurement using a guarded cut-bar technique, part 1 experimental study

In the guarded cut-bar technique, a guard surrounding the measured sample and reference (meter) bars is temperature controlled to carefully regulate heat losses from the sample and reference bars. Guarding is typically carried out by matching the temperature profiles between the guard and the test stack of sample and meter bars. Problems arise in matching the profiles, especially when the thermal conductivitiesof the meter bars and of the sample differ, as is usually the case. In a previous numerical study, the applied guarding condition (guard temperature profile) was found to be an important factor in measurement accuracy. Different from the linear-matched or isothermal schemes recommended in literature, the optimal guarding condition is dependent on the system geometry and thermal conductivity ratio of sample to meter bar. To validate the numerical results, an experimental study was performed to investigate the resulting error under different guarding conditions using stainless steel 304 as both the sample and meter bars. The optimal guarding condition was further verified on a certified reference material, pyroceram 9606, and 99.95% pure iron whose thermal conductivities are much smaller and much larger, respectively, than that of the stainless steel meter bars. Additionally, measurements are performed using three different inert gases to show the effect of the insulation effective thermal conductivity on measurement error, revealing low conductivity, argon gas, gives the lowest error sensitivity when deviating from the optimal condition. The result of this study provides a general guideline for the specific measurement method and for methods requiring optimal guarding or insulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hydrodynamic Design Optimization and Wave Tank Testing of a Self-Reacting Two-Body Wave Energy Converter

As worldwide energy consumption continues to increase, so does the demand for renewable energy sources. The total available wave energy resource for the United States alone is 2,640 TWh/yr; nearly two thirds of the 4,000 TWh of electricity used in the United States each year. It is estimated that nearly half of that available energy is recoverable through wave energy conversion techniques. In this thesis, a two-body ‘point absorber’ type wave energy converter with a mechanical power-takeoff is investigated. The two-body wave energy converter extracts energy through the relative motion of a floating buoy and a neutrally buoyant submerged body. Using a linear frequency-domain model, analytical solutions of the optimal power and the corresponding power-takeoff components are derived for the two-body wave energy converter. Using these solutions, a case study is conducted to investigate the influence of the submerged body size on the absorbed power of the device in regular and irregular waves. Here it is found that an optimal mass ratio between the submerged body and floating buoy exists where the device will achieve resonance. Furthermore, a case study to investigate the influence of the submerged body shape on the absorbed power is conducted using a time-domain numerical model. Here it is found that the submerged body should be designed to reduce the effects of drag, but to maintain relatively large hydrodynamic added mass and excitation force. To validate the analytical and numerical models, a 1/30th scale model of a two-body wave energy converter is tested in a wave tank. The results of the wave tank tests show that the two-body wave energy converter can absorb nearly twice the energy of a single-body ‘point absorber’ type wave energy converter.

16 TIDAL AND WAVE POWER↗

Development and clinical evaluation of an ionization chamber array with 3.5 mm pixel pitch for quality assurance in advanced radiotherapy techniques

Purpose: To characterize a new air vented ionization chamber technology, suitable to build detector arrays with small pixel pitch and independence of sensitivity on dose per pulse. Methods: The prototype under test is a linear array of air vented ionization chambers, consisting of 80 pixels with 3.5 mm pixel pitch distance and a sensitive volume of about 4 mm{sup 3}. The detector has been characterized with {sup 60}Co radiation and MV x rays from different linear accelerators (with flattened and unflattened beam qualities). Sensitivity dependence on dose per pulse has been evaluated under MV x rays by changing both the source to detector distance and the beam quality. Bias voltage has been varied in order to evaluate the charge collection efficiency in the most critical conditions. Relative dose profiles have been measured for both flattened and unflattened distributions with different field sizes. The reference detectors were a commercial array of ionization chambers and an amorphous silicon flat panel in direct conversion configuration. Profiles of dose distribution have been measured also with intensity modulated radiation therapy (IMRT), stereotactic radiosurgery (SRS), and volumetric modulated arc therapy (VMAT) patient plans. Comparison has been done with a commercial diode array and with Gafchromic EBT3 films. Results: Repeatability and stability under continuous gamma irradiation are within 0.3%, in spite of low active volume and sensitivity (∼200 pC/Gy). Deviation from linearity is in the range [0.3%, −0.9%] for a dose of at least 20 cGy, while a worsening of linearity is observed below 10 cGy. Charge collection efficiency with 2.67 mGy/pulse is higher than 99%, leading to a ±0.9% sensitivity change in the range 0.09–2.67 mGy/pulse (covering all flattened and unflattened beam qualities). Tissue to phantom ratios show an agreement within 0.6% with the reference detector up to 34 cm depth. For field sizes in the range 2 × 2 to 15 × 15 cm{sup 2}, the output factors are in agreement with a thimble chamber within 2%, while with 25 × 25 cm{sup 2} field size, an underestimation of 4.0% was found. Agreement of field and penumbra width measurements with the flat panel is of the order of 1 mm down to 1 × 1 cm{sup 2} field size. Flatness and symmetry values measured with the 1D array and the reference detectors are comparable, and differences are always smaller than 1%. Angular dependence of the detector, when compared to measurements taken with a cylindrical chamber in the same phantom, is as large as 16%. This includes inhomogeneity and asymmetry of the design, which during plan verification are accounted for by the treatment planning system (TPS). The detector is capable to reproduce the dose distributions of IMRT and VMAT plans with a maximum deviation from TPS of 3.0% in the target region. In the case of VMAT and SRS plans, an average (maximum) deviation of the order of 1% (4%) from films has been measured. Conclusions: The investigated technology appears to be useful both for Linac QA and patient plan verification, especially in treatments with steep dose gradients and nonuniform dose rates such as VMAT and SRS. Major limitations of the present prototype are the linearity at low dose, which can be solved by optimizing the readout electronics, and the underestimation of output factors with large field sizes. The latter problem is presently not completely understood and will require further investigations.

60 APPLIED LIFE SCIENCES↗