Search NASA⌕ Search

SEARCH · Search NASA

Results for “(non)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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Robust optimization of flexible diafiltration systems for critical mineral separations

This paper provides major contributions in expanding the literature for membrane process design with critical mineral recovery applications and showcasing the importance of robust design techniques for reducing risks of underperformance in such systems. Here, a membrane process flowsheet featuring PrOMMiS membrane models for recovering lithium/cobalt from spent batteries is showcased, uncertainty in membrane sieving and localized fouling are considered, and robust designs are obtained using the PyROS toolset. This paper is intended for a general audience of researchers working in critical minerals, membranes, and optimization related areas.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimization with Neural Network Feasibility Surrogates: Formulations and Application to Security-Constrained Optimal Power Flow

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Joint inversion of electromagnetic measurements for the determination of water saturation profiles in concrete structures

Highlights: • DC-electrical and dielectric data are combined to estimate concrete saturation degree. • A joint inversion approach of the electromagnetic measurements is proposed. • The sensitivity of the measurements to the saturation model parameters is analysed. • The new joint approach was developed and is applied to synthetic and real data. • The benefits of the joint approach over the inversion of one data type are highlighted. Water saturation profiles in concrete are essential to assess its durability and can be determined using non-destructive techniques, especially the electric and the capacitive methods. In this paper, we propose a new inversion scheme where both resistivity and permittivity measurements are inverted jointly to retrieve the saturation profile. The finite element method is used to model the measurements in 3D, the concrete having a saturation profile with depth, represented by a continuous model taking the form of a Weibull curve with four parameters. A non-linear least-squares optimization based on the Levenberg-Marquardt scheme is developed for the inversion of measurements. Results show that information gathered from both measurements enriches the reconstructed profile, leading to a more reliable estimation of saturation profiles. We believe that the joint inversion method herein developed could lead to the study of more complex phenomena, such as the coupled water-chloride ingress.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An analytical approach to evaluating magnetic-field closure and topological changes in FRC devices

Here, we describe mathematical methods based on optimizing a modified non-linear flux function (MFF) to evaluate whether odd-parity perturbations affect the local closure of magnetic field lines in field-reversed configurations. Using the MFF methodology, quantitative formulas are derived that provide the shift of the field minimum and the threshold for field-line opening, a discontinuous change in field topology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MULTI-LEADER: MULTI-source LEarning-Accelerated Design of high-Efficiency multi-stage compRessor (Final Technical Report)

The objective of MULTI-LEADER is to cut design costs by 80% while generating more energy-efficient designs of multi-stage compressors by developing and implementing novel machine learning (ML) techniques, which enable faster and fewer design iterations, improved solver performance, and concurrent multi-disciplinary design. Current industrial practices for the design of multi-stage compressors involve simulation-based design optimization with successive levels of model fidelity, iteratively evaluated between distinct disciplines, one stage at a time to tackle the high dimensional design variations. This project addresses these key design challenges: (1) concurrent optimization of multiple stages under many non-linear constraints; (2) multitude of evaluation of high-fidelity and expensive solvers and their gradients during optimization convergence in high-dimensional design; (3) multi-disciplinary design to maximize aerodynamic performance while guaranteeing structural integrity and additive manufacturability; (4) utilization of multiple fidelity of solvers with disparate parameterization and modeling assumptions. MULTI-LEADER achieved more than 5x speed up in detailed design of more energy-efficient compressors via these machine learning (ML) innovations: (i) rapid design surrogates by multi-source learning from diverse fidelities across multiple disciplines, (ii) physics-constrained data-augmented modeling for improved empiricism, (iii) generative manifold embedding for high dimensional concurrent design without gradient information; (iv) budget-constrained fidelity-adaptive sampling towards fewer design iterations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Proliferated Resilient Economical Half-Meter Aperture Space Telescopes (PREEMPT)

The PREEMPT LDRD was motivated by a major national need for better space-based imaging systems that are both high performing and affordable for the US Government. Current optical payloads for intelligence, surveillance and reconnaissance (ISR) and space domain awareness (SDA) missions can cost hundreds of millions of dollars and take years to develop, which makes it difficult to build the large constellations required for persistent, world-wide coverage. To address this, the team aimed to advance a different kind of large-aperture (>25 cm) telescope, called a monolithic telescope, in which key optical surfaces are built into a single piece of fused silica. This design greatly reduces payload and spacecraft complexity, the need for precision focus actuators, improves mechanical and thermal robustness, and lowers cost when compared with traditional Cassegrain telescopes that rely on many precisely aligned components. The project focused on five primary thrusts. The first thrust was to advance the concept of a V10, 25 cm, monolithic telescope forward from optical design to flight-ready stage. This was accomplished in partnership with Optimax, who delivered the first test unit in the early stages of the LDRD. The team developed several technologies necessary for this optic to be integrated into a flight demonstration. These include carbon fiber housings, highly detailed structural and thermal models and stress-reducing elastic averaging Hirth groove designs. These technologies resulted in a successful maturation of the optic, which is now slated to fly in late 2026/early 2027 for a demonstration mission. The second thrust was to advance the manufacturability of these optics. In collaboration with NIF’s optical manufacturing shop, we reduced polishing time from 480 hours to 65 hours through the implementation of optimized processes and new tools. The NIF team utilized a conceptual V8 (18 cm) optic to demonstrate this optimization, though it can be applied to the rest of the monolithic optic portfolio. Third, the team developed the first conceptual V20 (50 cm) payload, which is slated to be the next generation of LLNL optical payload systems. A set of structural, dynamic and thermal simulations were performed to identify potential challenges in the future development of this payload. Early-stage simulations suggest the payload is feasible, though thermal management will be the key focus area to maintain optimal performance. Fourth, a non-linear model of Viton was developed, to further enhance the reliability of our structural and dynamic models for future payloads. Viton acts as the primary interface material between the optic and its housing. Lastly, the team focused on successfully displaying the feasibility of using additively manufactured metal composites for optical space payloads. The team successfully demonstrated layer by layer deposition of Al-SiC composites, which have highly tunable structural and coefficient of thermal expansion (CTE) properties. These are crucial for optical payloads because CTE mismatch is one of the causes for degraded optical performance for telescopes in orbit. Overall, the work showed that monolithic telescopes could become a practical, lower-cost path to high-resolution space imaging for both national security and scientific missions.

42 ENGINEERING↗
DOE OSTI2023

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↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Peak Power Reduction for HVAC Operations in Multi-unit Commercial Buildings

The load profiles of most commercial consumers are characterized by brief periods of very high power consumption followed by intervals of relatively lower demand. In order to flatten commercial load profiles, several power utilities in addition to billing energy consumption, levy a demand charge (DC) on the monthly peak demand. In this work, we consider the problem of joint optimization of energy costs (EC) and DC incurred by a multi-unit building which follows a demand response (DR) program. Despite the non-linear structure of the problem, we show how the optimal solutions can be obtained efficiently using linear programming. We evaluate the performance of the proposed power control scheme for various climate zones in the US. We show that depending on the ambient conditions and the prescribed tariff structure, our strategy can result in savings of up to nearly 19% compared to the baseline.

Raza Naqvi, Syed A.↗
DOE OSTI2022

Peak Power Minimization for Commercial Thermostatically Controlled Loads in Multi-Unit Grid-Interactive Efficient Buildings

The load profiles of most commercial and industrial consumers are characterized by brief periods of very high power consumption followed by intervals of lower demand. To encourage such consumers to flatten their load profiles, power utilities in and around the world often levy a monthly demand charge (DC) on the peak demand measured over brief intervals. In this work, we consider the joint optimization of energy costs (EC) and the instantaneous peak power of a multi-unit building which uses a hydronic heating, ventilation and cooling (HVAC) system and responds to a demand response (DR) program. Despite the non-linear structure of the problem, we show how optimal solutions can be obtained efficiently using linear programming. Next, we study the power demand patterns resulting from our proposed strategy for thermostatically controlled loads (TCLs), and evaluate the strategy’s performance for various climate zones in the US, under both typical and atypical weather conditions. Finally, the results show that depending on the ambient conditions and the tariff structure, our strategy can result in utility bill savings of up to nearly 19% compared to the baseline. The results also indicate that our power control strategy can significantly reduce the instantaneous peak power consumption in commercial TCLs.

HVAC↗
DOE OSTI2021

Biochemical inhibition of acid phosphatase activity in two mountain spruce forest soils

The product inhibition of the 4-Methylumbelliferyl phosphate (MUB-P) decay rate, a measure of potential acid phosphatase activity, has not been considered in most of the published kinetic studies. The aim of this study was to determine the type and strength of the product inhibition in order to better define reaction conditions at which the acid phosphatase activity assay produces unbiased results. The MUB-P decay rate was measured in the forest floor and topsoil organic horizons of two spruce forest catchments at different initial MUB-P and P-PO4 3- concentrations. Furthermore, the type and strength of the inhibition was analyzed by non-linear regression. We found that MUB-P decay was competitively inhibited by the P-PO4 3- but also by dissolved organic P. By using estimated kinetic parameters, we calculated an underestimation of potential acid phosphatase activity by up to ~20% at concentrations of added MUB-P as high as 300 µmol g-1 dry weight and an incubation time one hour due to background inorganic and organic P concentrations in range of units of µmol per gram of dry soil. We discuss approaches that can be used to minimize the underestimation and show that the previously defined recommendations might not be always sufficient to avoid the bias. Therefore, we recommend to analyze progress curves at a wide range of initial MUB-P concentrations using non-linear methods when the enzyme assay is optimized for a given soil. When changes in inorganic and/or organic P concentrations are expected, we further recommend measuring background inorganic and organic P concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

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]↗
DOE OSTI2024

Third-integer Resonant Extraction Regulation System for Mu2e

A third-integer resonant slow extraction system is being developed for Fermilab's Delivery Ring to deliver protons to the upcoming Mu2e experiment. The timescale of the extraction (or spill) duration is 43 milliseconds, which is extremely short and unprecedented. Additionally, the experiment's strict and challenging requirements on the quality of the spill at this time scale has led to the development of a new Spill Regulation System (SRS) design. The SRS primarily consists of three components - slow regulation, fast regulation, and harmonic content suppressor. Contributions to the first two components of the SRS, i.e., Slow Regulation and Fast Regulation subsystems, will be presented in which new adaptive learning algorithm schemes for the slow regulation of the spill -- validated using particle tracking simulations -- shall be described. In addition to these novel methods for the enhancement of the spill regulation system, results of employing Machine Learning in enhancing the performance of the resonant extraction are also presented. At the forefront of applying ML techniques to solve non-linear accelerator control problems, this work includes optimizing the PID gains as well as the replacement of the traditional PID controller using Recurrent Neural Networks and Gated Recurrent Unit (GRU) ML models to achieve efficiencies greater than a PID controller. Cutting-edge on-going Reinforcement Learning efforts, including an actor-critic family of learning algorithms, to regulate the spill rate will be reviewed, as well as present analytical calculations pertaining the transit time of particles in a third-integer resonant extraction. Detailed numerical investigations and validations of such calculations, the model of which could be exported and reliably used in future analytical modeling of any resonant extraction, are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Online accelerator optimization with a machine learning-based stochastic algorithm

Abstract Online optimization is critical for realizing the design performance of accelerators. Highly efficient stochastic optimization algorithms are needed for many online accelerator optimization problems in order to find the global optimum in the non-linear, coupled parameter space. In this study, we propose to use the multi-generation Gaussian process optimizer for online accelerator optimization and demonstrate that the algorithm is significantly more efficient than other stochastic algorithms that are commonly used in the accelerator community.

Zhang, Zhe (ORCID:0000000281430381)↗
DOE OSTI2024

Uncertainty Quantification via Stable Distribution Propagation

We propose a new approach for propagating stable probability distributions through neural networks. Our method is based on local linearization, which we show to be an optimal approximation in terms of total variation distance for the ReLU non-linearity. This allows propagating Gaussian and Cauchy input uncertainties through neural networks to quantify their output uncertainties. To demonstrate the utility of propagating distributions, we apply the proposed method to predicting calibrated confidence intervals and selective prediction on out-of-distribution data. The results demonstrate a broad applicability of propagating distributions and show the advantages of our method over other approaches such as moment matching.

Artificial Intelligence (cs.AI)↗