Search NASA⌕ Search

SEARCH · Search NASA

Results for “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,747 records · Page 97

Accelerating matrix-centric graph processing on GPUs through bit-level optimizations

Even though it is well known that binary values are common in graph applications (e.g., adjacency matrix), how to leverage the phenomenon for efficiency has not yet been adequately explored. This paper presents a systematic study on how to unlock the potential of the bit-level optimizations of graph computations that involve binary values. It proposes a two-level representation named Bit-Block Compressed Sparse Row (B2SR) and presents a series of optimizations to the graph operations on B2SR by the intrinsics of modern GPUs. It additionally introduces Deep Reinforcement Learning (DRL) as an efficient way to best configure the bit-level optimizations on the fly. Additionally, the DQN-based adaptive tile size selector with dedicated model training can reach 68% prediction accuracy. Evaluations on NVIDIA Pascal and Volta GPUs show that the optimizations bring up to 40× and 6555× for essential GraphBLAS kernels SpMV and SpGEMM, respectively, making GraphBLAS-based BFS accelerate up to 433×, SSSP, PR, and CC up to 35×, and TC up to 52×.

79 ASTRONOMY AND ASTROPHYSICS↗

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗

Li-ion battery design through microstructural optimization using generative AI

Lithium-ion batteries are used across various applications, necessitating tailored cell designs to enhance performance. Optimizing electrode manufacturing parameters is a key route to achieving this, as these parameters directly influence the microstructure and performance of the cells. However, linking process parameters to performance is complex, and experimental or modeling campaigns are often slow and expensive. This study introduces a fast computational optimization framework for electrode manufacturing parameters. A generative model, trained on a small dataset of microstructural images associated with different manufacturing parameters, efficiently generates representative microstructures for new parameters. This model is integrated into a Bayesian optimization loop that includes microstructure generation, characterization, and simulation, aiming to find optimal manufacturing parameters for a particular application. Significant improvement in the energy density of a 4680 cell is achieved through bespoke cell design, highlighting the importance of cell-scale normalization. The framework’s modularity allows its application to various advanced materials manufacturing scenarios.

batteries↗

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Capacity assessment and co-optimization of CO 2 storage and enhanced oil recovery in residual oil zones

Residual oil zones (ROZs) are increasingly being commercially exploited using CO 2 -enhanced oil recovery (CO 2 -EOR) method. In this study, CO 2 storage potential, long-term CO 2 fate and oil recovery potential in ROZs are characterized based on a reservoir model for Goldsmith-Landreth San Andres Unit in the Permian Basin. The effects of CO 2 injection rates, well patterns (five-spot and line-drive), well spacings, injection modes (continuous CO 2 injection and water-alternatinggas injection) on the CO 2 retention in the reservoir and the oil production are investigated. After the preliminary assessment of CO 2 storage and EOR potentials in ROZs, we next develop a novel approach based on a newly developed optimization algorithm-Stochastic Simplex Approximate Gradient (StoSAG) and predictive empirical models constructed using machine learning technique to co-optimize CO 2 storage and oil recovery in ROZs. The performance of co-optimization of CO 2 storage and oil recovery is compared with the performance of optimization of only CO 2 storage.

58 GEOSCIENCES↗

Economic dispatch of offshore renewable energy resources for islanded communities with optimal storage sizing

Coastal or isolated microgrids depend on diesel generators and could benefit from renewable energy resources, especially offshore wind and wave energy. Integrating these resources into microgrids is complicated by their high intermittency, which requires optimal economic dispatch to effectively evaluate. This study considers three coastal or islanded sites, and uses mid-fidelity models of wind and wave energy technologies, and local demand data to solve the optimal economic dispatch problem. An optimal storage sizing method is developed that finds the smallest capacity of energy storage required to meet the microgrid load during each season. The storage capacity decreases by a factor of two at most when adding wave energy converters to a system. Adding wave energy converters to a farm decreases cost by about 30%. Furthermore, the required storage size varies by two to three times from summer to winter. Compared with the state-of-the-art approaches that often overlook realistic offshore renewable energy technology in microgrid economic dispatch and optimal storage sizing, the proposed solution introduced in this study allows for better site selection, microgrid design, converter selection, and storage sizing considerations for isolated microgrids.

16 TIDAL AND WAVE POWER↗

Joint Optimization of Multimodal Transit Frequency and Shared Autonomous Vehicle Fleet Size with Hybrid Metaheuristic and Nonlinear Programming

Shared autonomous vehicles (SAVs) bring competition to traditional transit services but redesigning multimodal transit network can utilize SAVs as feeders to enhance service efficiency and coverage. This paper presents an optimization framework for the joint multimodal transit frequency and SAV fleet size problem, a variant of the transit network frequency setting problem. The objective is to maximize total transit ridership (including SAV-fed trips and subtracting boarding rejections) across multiple time periods under budget constraints, considering endogenous mode choice (transit, point-to-point SAVs, driving) and route selection, while allowing for strategic route removal by setting frequencies to zero. Due to the problem’s non-linear, non-convex nature and the computational challenges of large-scale networks, we develop a hybrid solution approach that combines a metaheuristic approach (particle swarm optimization) with nonlinear programming for local solution refinement. To ensure computational tractability, the framework integrates analytical approximation models for SAV waiting times based on fleet utilization, multimodal network assignment for route choice, and multinomial logit mode choice behavior, bypassing the need for computationally intensive simulations within the main optimization loop. Applied to the Chicago metropolitan area’s multimodal network, our method illustrates a 33.3% increase in transit ridership through optimized transit route frequencies and SAV integration, particularly enhancing off-peak service accessibility and strategically reallocating resources.

Ng, Max↗

Optimization of the IPP-bypass mevalonate pathway and fed-batch fermentation for the production of isoprenol in Escherichia coli

Isoprenol (3-methyl-3-buten-1-ol) is a drop-in biofuel and a precursor for commodity chemicals. Biological production of isoprenol via the mevalonate pathway has been developed and optimized extensively in Escherichia coli, but high ATP requirements and isopentenyl diphosphate (IPP) toxicity have made it difficult to achieve high titer, yield, and large-scale production. To overcome these limitations, an IPP-bypass pathway was previously developed using the promiscuous activity of diphosphomevalonate decarboxylase, and enabled the production of isoprenol at a comparable yield and titer to the original pathway. In this study, we optimized this pathway, substantially improving isoprenol production. A titer of 3.7 g/L (0.14 g isoprenol per g glucose) was achieved in batch conditions using minimal medium by pathway optimization, and a further optimization of the fed-batch fermentation process enabled an isoprenol titer of 10.8 g/L (yield of 0.105 g/g and maximum productivity of 0.157 g L -1 h -1 ), which is the highest reported titer for this compound. The substantial increase in isoprenol titer via the IPP-bypass pathway in this study will facilitate progress toward commercialization.

59 BASIC BIOLOGICAL SCIENCES↗

Direct construction of optimized stellarator shapes. Part 2. Numerical quasisymmetric solutions

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding-centre particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field strength $B$ . Here, instead, we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in the distance from the magnetic axis, without using optimization. The approach involves solution of a one-dimensional nonlinear ordinary differential equation, originally derived by Garren & Boozer ( Phys. Fluids B, vol. 3, 1991, p. 2805). We illustrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cost Optimization of Osmotically Assisted Reverse Osmosis

We develop a nonlinear optimization model to identify minimum cost designs for osmotically assisted reverse osmosis (OARO), a multistaged membrane-based process for desalinating high-salinity brines. The optimization model enables comprehensive evaluation of a complex process configuration and operational decision space that includes nonlinear process performance and implicit relationships among membrane stages, saline sweep cycles, and makeup, purge, and recycle streams. Here, the objective function minimizes cost, rather than energy or capital expenditures, to accurately account for the trade-offs in capital and operational expenses inherent in multistaged membrane processes. Generally, we find that cost-optimal OARO processes minimize the number of stages, eliminate the use of saline makeup streams, purge from the first sweep cycle, and successively decrease stage membrane area and sweep flow rates. The optimal OARO configuration for treating feed salinities of 50–125 g/L total dissolved solids with water recoveries between 30–70% results in costs less than or equal to $6 per m3 of product water. Sensitivity analysis suggests that future research to minimize OARO costs should focus on minimizing the membrane structural parameter while maximizing the membrane burst pressure and reducing the membrane unit cost.

42 ENGINEERING↗

A Julia Framework for Graph-Structured Nonlinear Optimization

Graph theory provides a convenient framework for modeling and solving structured optimization problems. Under this framework, the modeler can arrange/assemble the components of an optimization model (variables, constraints, objective functions, and data) within nodes and edges of a graph, and this representation can be used to visualize, manipulate, and solve the problem. In this work, we present a Julia framework for modeling and solving graph-structured nonlinear optimization problems. Our framework integrates the modeling package Plasmo.jl (which facilitates the construction and manipulation of graph models) and the nonlinear optimization solver MadNLP.jl (which provides capabilities for exploiting graph structures to accelerate solution). We illustrate with a simple example how model construction and manipulation can be performed in an intuitive manner using Plasmo.jl and how the model structure can be exploited by MadNLP.jl. We also demonstrate the scalability of the framework by targeting a large-scale, stochastic gas network problem that contains over 1.7 million variables.

Cole, David↗

Tracking Excited States in Wave Function Optimization Using Density Matrices and Variational Principles

We introduce an approach for finding individual excited states' energy stationary points in complete active space self-consistent field theory that is compatible with standard optimization methods and highly effective at overcoming difficulties due to root flipping and near-degeneracies. Inspired by both the maximum overlap method and recent progress in excited-state variational principles, our approach combines these ideas in order to track individual excited states throughout the orbital optimization process. In a series of tests involving root flipping, near-degeneracies, charge transfers, and double excitations, we show that this approach is more effective for state-specific optimization than either the naive selection of roots on the basis of energy ordering or a more direct generalization of the maximum overlap method. We offer evidence that this state-specific approach improves the performance of complete active space perturbation theory for vertical excitation energies. Furthermore, we find that the state-specific optimization can help avoid state-averaging-induced discontinuities on potential energy surfaces. With a simple implementation, a low cost, and compatibility with large active space methods, the approach is designed to be useful in a wide range of excited-state investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Engineering Antifreeze Proteins to Optimally Resist Engulfment by Ice

Antifreeze proteins (AFPs) facilitate the survival of organisms in cold climates by inhibiting the growth and/or recrystallization of ice. To function, AFPs must first bind to ice crystals; bound AFPs must then resist engulfment by using their nonbinding side (NBS) to pin the ice–water interface. Here, we seek to understand how the molecular characteristics of an NBS, such as its ice-phobicity or shape, influence its ability to resist engulfment. By characterizing the free energy barriers that impede the engulfment of model AFPs, we find that the critical supercooling ΔT*, above which an AFP is engulfed, is dictated by an optimal pinning site on the NBS. We further find that the optimal pinning site is determined by an interplay between the contact line perimeter P and a pinning efficiency η, with ΔT* ∝ ηP at the optimal pinning site. For a hemispherical AFP, which displays progressively inward tapering, we find that η increases during engulfment, whereas P decreases; conversely, an NBS with outward tapering can achieve high P, but it suffers from low η. Because the product of η and P determines ΔT*, the inverse correlation between them limits ΔT. To circumvent such limiting behavior, we propose an NBS shape with an outward bulge; by initially tapering outward, a bulged NBS permits higher P, and by subsequently tapering inward, it promotes high η as well. Importantly, we find that ΔT* is enhanced by more than a factor of 2 with an outward bulge of only 1 nm. We also find that the more ice-phobic an NBS is, the more efficiently it pins the ice–water interface, resulting in a higher ΔT. Furthermore, our findings shed light on how the NBS molecular characteristics influence ΔT* and suggest strategies for engineering the NBS to optimally resist engulfment by ice.

Antifreeze↗

Evolution of a Novel Ribbon Phase in Optimally Doped Bi 2 Sr 2 CaCu 2 O 8+δ at High Pressure and Its Implication to High-T C Superconductivity

One challenge in studying high-temperature superconductivity (HTSC) stems from a lack of direct experimental evidence linking lattice inhomogeneity and superconductivity. Here in this paper, we apply synchrotron hard X-ray nanoimaging and small-angle scattering to reveal a novel micron-scaled ribbon phase in optimally doped Bi 2 Sr 2 CaCu 2 O 8+δ (Bi-2212, with δ = 0.1). The morphology of the ribbon-like phase evolves simultaneously with the dome-shaped T C behavior under pressure. X-ray absorption studies show that the increasing of T C is associated with oxygen-hole redistribution in the CuO 2 plan, while T C starts to decrease with pressure when oxygen holes become immobile. Additional X-ray irradiation experiments reveal that nanoscaled short-range ordering of oxygen vacancies could further lower TC, which indicates that the optimal T C is affected not only by an optimal morphology of the ribbon phase, but also an optimal distribution of oxygen vacancies. Our studies thereby provide for the first time compelling experimental evidence correlating the T C with micron to nanoscale inhomogeneity.

36 MATERIALS SCIENCE↗

Investigating and Optimizing the Lysate-Based Expression of Nonribosomal Peptide Synthetases Using a Reporter System

Lysate-based cell-free expression (CFE) systems are accessible platforms for expressing proteins that are difficult to synthesize in vivo, such as nonribosomal peptide synthetases (NRPSs). NRPSs are large (>100 kDa), modular enzyme complexes that synthesize bioactive peptide natural products. This synthetic process is analogous to transcription/translation (TX/TL) in lysates, resulting in potential resource competition between NRPS expression and NRPS activity in cell-free environments. Moreover, CFE conditions depend on the size and structure of the protein. Here, a reporter system for rapidly investigating and optimizing reaction environments for NRPS CFE is described. This strategy is demonstrated in E. coli lysate reactions using blue pigment synthetase A (BpsA), a model NRPS, carrying a C-terminal tetracysteine (TC) tag which forms a fluorescent complex with the biarsenical dye, FlAsH. A colorimetric assay was adapted for lysate reactions to detect the blue pigment product, indigoidine, of cell-free expressed BpsA-TC, confirming that the tagged enzyme is catalytically active. An optimized protocol for end point TC/FlAsH complex measurements in reactions enables quick comparisons of full-length BpsA-TC expressed under different reaction conditions, defining unique requirements for NRPS expression that are related to the protein’s catalytic activity and size. Importantly, these protein-dependent CFE conditions enable higher indigoidine titer and improve the expression of other monomodular NRPSs. Notably, these conditions differ from those used for the expression of superfolder GFP (sfGFP), a common reporter for optimizing lysate-based CFE systems, indicating the necessity for tailored reporters to optimize expression for specific enzyme classes. In conclusion, the reporter system is anticipated to advance lysate-based CFE systems for complex enzyme synthesis, enabling natural product discovery.

59 BASIC BIOLOGICAL SCIENCES↗

Optimizing Carbon Cycle Parameters Drastically Improves Terrestrial Biosphere Model Underestimates of Dryland Mean Net CO 2 Flux and its Inter-Annual Variability

Drylands occupy ~40% of the land surface and are thought to dominate global carbon (C) cycle inter-annual variability (IAV). Therefore, it is imperative that global terrestrial biosphere models (TBMs), which form the land component of IPCC earth system models, are able to accurately simulate dryland vegetation and biogeochemical processes. However, compared to more mesic ecosystems, TBMs have not been widely tested or optimized using in situ dryland CO 2 fluxes. Here, for this work, we address this gap using a Bayesian data assimilation system and 89 site-years of daily net ecosystem exchange (NEE) data from 12 southwest US Ameriflux sites to optimize the C cycle parameters of the ORCHIDEE TBM. The sites span high elevation forest ecosystems, which are a mean sink of C, and low elevation shrub and grass ecosystems that are either a mean C sink or “pivot” between an annual C sink and source. We find that using the default (prior) model parameters drastically underestimates both the mean annual NEE at the forested mean C sink sites and the NEE IAV across all sites. Our analysis demonstrated that optimizing phenology parameters are particularly useful in improving the model's ability to capture both the magnitude and sign of the NEE IAV. At the forest sites, optimizing C allocation, respiration, and biomass and soil C turnover parameters reduces the underestimate in simulated mean annual NEE. Our study demonstrates that all TBMs need to be calibrated for dryland ecosystems before they are used to determine dryland contributions to global C cycle variability and long-term carbon-climate feedbacks.

59 BASIC BIOLOGICAL SCIENCES↗