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 181 records · Page 10

Optimal Power Flow Derived Sparse Linear Solver Benchmarks

Due to the changing nature of the power grid, it is increasingly important to be able to solve a high-fidelity optimal power-flow models on large power networks. This high-fidelity problem, called AC Optimal Power Flow (ACOPF), is a nonlinear, nonconvex optimization problem. One of the few reliable ways of solving such a problem is interior point methods. These methods result in sparse linear systems where the coefficient matrix is symmetric, indefinite and nearly always ill-conditioned. As such, they are particularly challenging for sparse linear solvers and represent a considerable computational bottleneck in solving the ACOPF problem. In this paper, we introduce a repository of linear systems captured from ACOPF problems when solved by the open-source optimizer IPOPT. These matrices are meant to be used as a test suite for sparse linear solver development.

97 MATHEMATICS AND COMPUTING↗

Singular linear-quadratic control problem for systems with linear delay

A singular linear-quadratic optimization problem on the trajectories of non-autonomous linear differential equations with linear delay is considered. The peculiarity of this problem is the fact that this problem has no solution in the class of integrable controls. To ensure the existence of solutions is required to expand the class of controls including controls with impulse components. Dynamical systems with linear delay are used to describe the motion of pantograph from the current collector with electric traction, biology, etc. It should be noted that for practical problems fact singularity criterion of quality is quite commonly occurring, and therefore the study of these problems is surely important. For the problem under discussion optimal programming control contained impulse components at the initial and final moments of time is constructed under certain assumptions on the functional and the right side of the control system.

97 MATHEMATICS AND COMPUTING↗

Linear Energy Transfer-Guided Optimization in Intensity Modulated Proton Therapy: Feasibility Study and Clinical Potential

Purpose: To investigate the feasibility and potential clinical benefit of linear energy transfer (LET) guided plan optimization in intensity modulated proton therapy (IMPT). Methods and Materials: A multicriteria optimization (MCO) module was used to generate a series of Pareto-optimal IMPT base plans (BPs), corresponding to defined objectives, for 5 patients with head-and-neck cancer and 2 with pancreatic cancer. A Monte Carlo platform was used to calculate dose and LET distributions for each BP. A custom-designed MCO navigation module allowed the user to interpolate between BPs to produce deliverable Pareto-optimal solutions. Differences among the BPs were evaluated for each patient, based on dose–volume and LET–volume histograms and 3-dimensional distributions. An LET-based relative biological effectiveness (RBE) model was used to evaluate the potential clinical benefit when navigating the space of Pareto-optimal BPs. Results: The mean LET values for the target varied up to 30% among the BPs for the head-and-neck patients and up to 14% for the pancreatic cancer patients. Variations were more prominent in organs at risk (OARs), where mean LET values differed by a factor of up to 2 among the BPs for the same patient. An inverse relation between dose and LET distributions for the OARs was typically observed. Accounting for LET-dependent variable RBE values, a potential improvement on RBE-weighted dose of up to 40%, averaged over several structures under study, was noticed during MCO navigation. Conclusions: We present a novel strategy for optimizing proton therapy to maximize dose-averaged LET in tumor targets while simultaneously minimizing dose-averaged LET in normal tissue structures. MCO BPs show substantial LET variations, leading to potentially significant differences in RBE-weighted doses. Pareto-surface navigation, using both dose and LET distributions for guidance, provides the means for evaluating a large variety of deliverable plans and aids in identifying the clinically optimal solution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

New adaptive method to optimize the secondary reflector of linear Fresnel collectors

Performance of linear Fresnel collectors may largely depend on the secondary-reflector profile design when small-aperture absorbers are used. Optimization of the secondary-reflector profile is an extremely challenging task because there is no established theory to ensure superior performance of derived profiles. In this work, an innovative optimization method is proposed to optimize the secondary-reflector profile of a generic linear Fresnel configuration. The method correctly and accurately captures impacts of both geometric and optical aspects of a linear Fresnel collector to secondary-reflector design. The proposed method is an adaptive approach that does not assume a secondary shape of any particular form, but rather, starts at a single edge point and adaptively constructs the next surface point to maximize the reflected power to be reflected to absorber(s). As a test case, the proposed optimization method is applied to an industrial linear Fresnel configuration, and the results show that the derived optimal secondary reflector is able to redirect more than 90% of the power to the absorber in a wide range of incidence angles. Here, the proposed method can be naturally extended to other types of solar collectors as well, and it will be a valuable tool for solar-collector designs with a secondary reflector.

14 SOLAR ENERGY↗

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Wave Tank Testing Report for Controls Validation of a Heaving Point Absorber

The core objectives of this project is to improve the power capture of three different wave energy conversion (WEC) devices by more than 50% using an advanced control system and validate the attained improvements using wave tank and full scale testing. In parallel, we will bring along the development of a wave prediction system that is required to enable effective control and test it at full scale. The purposes of this report are to: 1. Plan and document the 1/25th scale device testing at the wave-tank facility; 2. Document the test article, setup and methodology, sensor and instrumentation, mooring, electronics, wiring, and data flow and quality assurance; 3. Communicate the testing results between the associated members; 4. Facilitate reviews that will help to ensure all aspects (risk, safety, testing procedures, etc.); 5. Provide a systematic guide to setting up, executing and decommissioning the experiment.

16 TIDAL AND WAVE POWER↗

Linear-in-temperature resistivity for optimally superconducting (Nd,Sr)NiO 2

The occurrence of superconductivity in proximity to various strongly correlated phases of matter has drawn extensive focus on their normal state properties, to develop an understanding of the state from which superconductivity emerges. The recent finding of superconductivity in layered nickelates raises similar interests. However, transport measurements of doped infinite-layer nickelate thin films have been hampered by materials limitations of these metastable compounds: in particular, a high density of extended defects. Here, by moving to a substrate (LaAlO 3 ) 0.3 (Sr 2 TaAlO 6 ) 0.7 that better stabilizes the growth and reduction conditions, we can synthesize the doping series of Nd 1– x Sr x NiO 2 essentially free from extended defects. In their absence, the normal state resistivity shows a low-temperature upturn in the underdoped regime, linear behaviour near optimal doping and quadratic temperature dependence for overdoping. This is phenomenologically similar to the copper oxides despite key distinctions—namely, the absence of an insulating parent compound, multiband electronic structure and a Mott–Hubbard orbital alignment rather than the charge-transfer insulator of the copper oxides. In this study, we further observe an enhancement of superconductivity, both in terms of transition temperature and range of doping. These results indicate a convergence in the electronic properties of both superconducting families as the scale of disorder in the nickelates is reduced.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Systems Analysis Approach to Polyethylene Terephthalate and Olefin Plastics Supply Chains in the Circular Economy: A Review of Data Sets and Models

The environmental and economic impacts of implementing a circular economy in plastic waste supply chains are not well understood. The proposed systems analysis framework assesses environmental, social, and economic impacts of plastic waste supply chains in a circular economy. The first objective of this article is to identify datasets, models and knowledge gaps associated with waste plastic supply chain processes, mainly in the U.S. Our literature review indicated that the best datasets exist for virgin plastic resin production, mechanical recycling, landfilling, and incineration, with materials recovery facility being intermediate, and with chemical recycling the lowest. The second objective of this article is to develop an illustrative application of the framework by conducting a preliminary systems analysis of PET bottles with closed-loop recycling. Here, the preliminary systems analysis of PET bottles utilized a linear programming optimization method. Our optimization model indicated that both chemical and mechanical recycling processes are needed to achieve a true circular economy of PET bottles with the least greenhouse gas emissions, specifically reductions of 24% when compared with the linear economy. Good quality and standardized life cycle assessment and techno-economic analysis studies are needed to better understand the environmental, economic, and social impacts of advanced sorting and chemical recycling technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cherenkov imaging method for rapid optimization of clinical treatment geometry in total skin electron beam therapy

Purpose: A method was developed utilizing Cherenkov imaging for rapid and thorough determination of the two gantry angles that produce the most uniform treatment plane during dual-field total skin electron beam therapy (TSET). Methods: Cherenkov imaging was implemented to gather 2D measurements of relative surface dose from 6 MeV electron beams on a white polyethylene sheet. An intensified charge-coupled device camera time-gated to the Linac was used for Cherenkov emission imaging at sixty-two different gantry angles (1° increments, from 239.5° to 300.5°). Following a modified Stanford TSET technique, which uses two fields per patient position for full body coverage, composite images were created as the sum of two beam images on the sheet; each angle pair was evaluated for minimum variation across the patient region of interest. Cherenkov versus dose correlation was verified with ionization chamber measurements. The process was repeated at source to surface distance (SSD) = 441, 370.5, and 300 cm to determine optimal angle spread for varying room geometries. In addition, three patients receiving TSET using a modified Stanford six-dual field technique with 6 MeV electron beams at SSD = 441 cm were imaged during treatment. Results: As in previous studies, Cherenkov intensity was shown to directly correlate with dose for homogenous flat phantoms (R{sup 2} = 0.93), making Cherenkov imaging an appropriate candidate to assess and optimize TSET setup geometry. This method provided dense 2D images allowing 1891 possible treatment geometries to be comprehensively analyzed from one data set of 62 single images. Gantry angles historically used for TSET at their institution were 255.5° and 284.5° at SSD = 441 cm; however, the angles optimized for maximum homogeneity were found to be 252.5° and 287.5° (+6° increase in angle spread). Ionization chamber measurements confirmed improvement in dose homogeneity across the treatment field from a range of 24.4% at the initial angles, to only 9.8% with the angles optimized. A linear relationship between angle spread and SSD was observed, ranging from 35° at 441 cm, to 39° at 300 cm, with no significant variation in percent-depth dose at midline (R{sup 2} = 0.998). For patient studies, factors influencing in vivo correlation between Cherenkov intensity and measured surface dose are still being investigated. Conclusions: Cherenkov intensity correlates to relative dose measured at depth of maximum dose in a uniform, flat phantom. Imaging of phantoms can thus be used to analyze and optimize TSET treatment geometry more extensively and rapidly than thermoluminescent dosimeters or ionization chambers. This work suggests that there could be an expanded role for Cherenkov imaging as a tool to efficiently improve treatment protocols and as a potential verification tool for routine monitoring of unique patient treatments.

61 RADIATION PROTECTION AND DOSIMETRY↗

Multiphase Optimal and Non-Singular Power Flow by Successive Linear Approximations

We propose an optimal power flow method for multiphase power systems. The method works for both radial and meshed networks, and is able to accommodate different types of load/source connections such as wye, delta, and combination thereof. Moreover, the method generates a sequence of nonsingular points (i.e., associated with non-singular power-flow Jacobian) in the electrical state space and a corresponding sequence of power injections, which can be used in real-time control settings. Under certain conditions, the sequences attain the local minimum. The method has been tested using several typical networks, including the IEEE 37-bus and 123-bus test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Large-scale quasi-Newton trust-region methods with low-dimensional linear equality constraints

We propose two limited-memory BFGS (L-BFGS) trust-region methods for large-scale optimization with linear equality constraints. Here, the methods are intended for problems where the number of equality constraints is small. By exploiting the structure of the quasi-Newton compact representation, both proposed methods solve the trust-region subproblems nearly exactly, even for large problems. We derive theoretical global convergence results of the proposed algorithms, and compare their numerical effectiveness and performance on a variety of large-scale problems.

97 MATHEMATICS AND COMPUTING↗

On the Efficacy and High-Performance Implementation of Quaternion Matrix Multiplication

Quaternion symmetry is ubiquitous in the physical sciences. As such, much work has been afforded over the years to the development of efficient schemes to exploit this symmetry using real and complex linear algebra. Recent years have also seen many advances in the formal theoretical development of explicitly quaternion linear algebra with promising applications in image processing and machine learning. Despite these advances, there do not currently exist optimized software implementations of quaternion linear algebra. The leverage of optimized linear algebra software is crucial in the achievement of high levels of performance on modern computing architectures, and thus provides a central tool in the development of high-performance scientific software. In this work, a case will be made for the efficacy of high-performance quaternion linear algebra software for appropriate problems. In this pursuit, an optimized software implementation of quaternion matrix multiplication will be presented and will be shown to outperform a vendor tuned implementation for the analogous complex matrix operation. The results of this work pave the path for further development of high-performance quaternion linear algebra software which will improve the performance of the next generation of applicable scientific applications.

97 MATHEMATICS AND COMPUTING↗

Algebraic multigrid preconditioning of the Hessian in optimization constrained by a partial differential equation

Summary We construct an algebraic multigrid (AMG) based preconditioner for the reduced Hessian of a linear‐quadratic optimization problem constrained by an elliptic partial differential equation. While the preconditioner generalizes a geometric multigrid preconditioner introduced in earlier works, its construction relies entirely on a standard AMG infrastructure built for solving the forward elliptic equation, thus allowing for it to be implemented using a variety of AMG methods and standard packages. Our analysis establishes a clear connection between the quality of the preconditioner and the AMG method used. The proposed strategy has a broad and robust applicability to problems with unstructured grids, complex geometry, and varying coefficients. The method is implemented using the Hypre package and several numerical examples are presented.

Barker, Andrew T.↗

Cascade energy optimization for waste heat recovery in distributed energy systems

The efficiency of distributed energy systems can be significantly increased through waste heat recovery from industry or power generation. The technologies used for this process are typically dependent on the quality and temperature grades of waste heat. To maximize the efficiency of cascade heat utilization, it is important to optimize the choice of waste heat recovery technologies and their operation. Here, a detailed mixed integer linear programming optimization model is proposed for waste heat recovery in a district-scale microgrid. The model can distinguish waste heat quality for planning and operation optimization of distributed energy systems. Heat utilization technologies are formulated in this developed model and categorized in different temperature grades. The developed model is validated using four typical cases under different settings of system operation and business models. It is found that the optimization model, by distinguishing waste heat temperature, can increase energy cost savings by around 5%, compared to models that do not consider waste heat temperature grades. Additionally, the results indicate that the developed model can provide more realistic configuration and technologies dispatch.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Depot Charging Schedule Optimization for Medium- and Heavy-Duty Battery-Electric Trucks

Charge management, which lowers charging costs for fleets and prevents straining the electrical grid, is critical to the successful deployment of medium- and heavy-duty battery-electric trucks (MHD BETs). This study introduces an energy demand and cost management framework that optimizes depot charging for MHD BETs by combining an energy consumption machine learning model and a linear program optimization model. The framework considers key factors impacting real-world MHD BET operations, including vehicle and charger configurations, duty cycles, use cases, geographic and climate conditions, operation schedules, and utilities’ time-of-use (TOU) rates and demand charges. The framework was applied to a hypothetical fleet of 100 MHD BETs in California under three different utilities for 365 days, with results compared to unmanaged charging. The optimized charging solution avoided more than 90% of on-peak charging, reduced fleet charging peak load by 64–75%, and lowered fleet energy variable costs by 54–64%. This study concluded that the proposed charge management framework significantly reduces energy costs and peak loads for MHD BET fleets while making recommendations for fleet electrification infrastructure planning and the design of utility TOU rates and demand charges.

Song, Shuhan↗

Continuous-Time Mean-Variance Portfolio Selection with Random Horizon

This paper examines the continuous-time mean-variance optimal portfolio selection problem with random market parameters and random time horizon. Treating this problem as a linearly constrained stochastic linear-quadratic optimal control problem, I explicitly derive the efficient portfolios and efficient frontier in closed forms based on the solutions of two backward stochastic differential equations. Some related issues such as a minimum variance portfolio and a mutual fund theorem are also addressed. All the results are markedly different from those in the problem with deterministic exit time. A key part of my analysis involves proving the global solvability of a stochastic Riccati equation, which is interesting in its own right.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct microstability optimization of stellarator devices

Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Scientific/Technical Report

In this work, we have built upon our results from previous DOE funding (DEFG 0204ER25655), where we developed new and more efficient methods for solving certain optimization problems with many inequality constraints. This past work resulted in efficient algorithms (and analysis of their convergence) for linear programming, convex quadratic programming, and the training of support vector machines. The algorithms are based on using constraint reduction in interior point methods: at each iteration we consider only a smaller subset of the inequality constraints, focusing on the constraints that are close enough to be relevant. Surprisingly, we have been able to show theoretically that such algorithms are globally convergent and to demonstrate experimentally that they are much more efficient than standard interior point methods.

97 MATHEMATICS AND COMPUTING↗