Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Analysis of Vector Particle-In-Cell (VPIC) memory usage optimizations on cutting-edge computer architectures

Vector Particle-In-Cell (VPIC) is one of the fastest plasma simulation codes in the world, with particle numbers ranging from one trillion on the first petascale system, Roadrunner, to ten trillion particles on the more recent Blue Waters supercomputer. As supercomputers continue to grow rapidly in size, so too does the gap between computing capability and memory capability. Current memory systems limit VPIC simulations greatly as the maximum number of particles that can be simulated directly depends on the available memory. In this study, we present a suite of VPIC memory optimizations (i.e., particle weight, half-precision, and fixed-point optimizations) that enable a significant increase in the number of particles in VPIC simulations. Here, we assess the optimizations’ impact on memory and runtime performance for a suite of cutting-edge computer architectures such has the NVIDIA V100 GPU, the IBM Power9, and the Fujitsu A64FX architectures. Our optimizations enable a 31.25% reduction in memory usage and up to 40% increase in the number of particles. This paper extends our work on developing particle storage format optimizations Tan et al.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Optimization of metal-supported solid oxide fuel cells with a focus on mass transport

Performance of symmetric-architecture metal-supported solid oxide fuel cells was improved significantly by optimizing the catalyst infiltration process and metal support structure. Optimization of component structure and processing parameters was performed during tape-casting and fabrication of button cells. Mass transport of oxygen in the metal support was identified as a major limitation. To overcome this limitation, pore former loading and thickness of the metal support (130–250 μm) were optimized. The catalyst infiltration process was also improved by studying the impact of firing temperature (400 °C–900 °C) and infiltration cycle numbers (1–15). The maximum power density of the optimized cell was 0.9 W cm –2 at 700°C using hydrogen as a fuel, a three-fold increase over the baseline cell performance. The degradation rate of optimized cells at 550 °C, 600 °C, and 700 °C was 2%, 4.5%, and 5.5% per 100 h, respectively. Here, the phenomena of mass transport, catalyst coarsening, and chromium poisoning on the catalyst were analyzed by electrochemical impedance spectroscopy and scanning electron microscopy.

13 HYDRO ENERGY↗
DOE OSTI2025

Property optimized energy absorber for automotive bumpers utilizing multi-material and structural design strategies

This study proposes a novel design for automotive bumper using optimized lattice structures and multi-materials to balance low-speed collision and high-speed pedestrian impact performance. Different blends of 20 % carbon fiber-reinforced acrylonitrile butadiene styrene with thermoplastic polyurethane were used to tailor material properties. The energy absorber features lattice structures with customized mechanical responses, created by varying the incline angle θ from 0 to 180°. We conducted 576 finite element simulations on a half-scale model to optimize energy absorption and stiffness, leading to 66 optimized designs that met both low-speed and high-speed impact criteria. Two sub-scale optimized energy absorbers with different peak forces—both meeting low-speed impact requirements—were 3D printed and validated through drop-weight testing. The one with lower peak stress demonstrated a more compliant response, exhibiting approximately 90 % lower initial peak force and an increase in energy absorption of around 33 % (from 24 J to 32 J). Compared to the baseline triangular lattice, the optimized absorber increased energy absorption by 68 % from (19 J to 32 J) and reduced peak stress by 70 %. It also showed near-complete recovery with minimal fractures, making it suitable for repeated use. This design improves safety while offering a lightweight, durable, and cost-effective bumper system.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Optimal beam margins in linac-based VMAT stereotactic ablative body radiotherapy: a Pareto front analysis for liver metastases

We explored the Pareto fronts mathematical strategy to determine the optimal block margin and prescription isodose for stereotactic body radiotherapy (SBRT) treatments of liver metastases using the volumetric-modulated arc therapy (VMAT) technique. Three targets (planning target volumes [PTVs] = 20, 55, and 101 cc) were selected. A single fraction dose of 26 Gy was prescribed (prescription dose [PD]). VMAT plans were generated for 3 different beam energies. Pareto fronts based on (1) different multileaf collimator (MLC) block margin around PTV and (2) different prescription isodose lines (IDL) were produced. For each block margin, the greatest IDL fulfilling the criteria (95% of PTV reached 100%) was considered as providing the optimal clinical plan for PTV coverage. Liver D{sub mean}, V7Gy, and V12Gy were used against the PTV coverage to generate the fronts. Gradient indexes (GI and mGI), homogeneity index (HI), and healthy liver irradiation in terms of D{sub mean}, V7Gy, and V12Gy were calculated to compare different plans. In addition, each target was also optimized with a full-inverse planning engine to obtain a direct comparison with anatomy-based treatment planning system (TPS) results. About 900 plans were calculated to generate the fronts. GI and mGI show a U-shaped behavior as a function of beam margin with minimal values obtained with a +1 mm MLC margin. For these plans, the IDL ranges from 74% to 86%. GI and mGI show also a V-shaped behavior with respect to HI index, with minimum values at 1 mm for all metrics, independent of tumor dimensions and beam energy. Full-inversed optimized plans reported worse results with respect to Pareto plans. In conclusion, Pareto fronts provide a rigorous strategy to choose clinical optimal plans in SBRT treatments. We show that a 1-mm MLC block margin provides the best results with regard to healthy liver tissue irradiation and steepness of dose fallout.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Control co-design under uncertainty for offshore wind farms: Optimizing grid integration, energy storage, and market participation

Offshore wind farms (OWFs) are set to significantly contribute to global decarbonization efforts. Developers often use a sequential approach to optimize design variables and market participation for grid-integrated offshore wind farms. However, this method can lead to sub-optimal system performance, and uncertainties associated with renewable resources are often overlooked in decision-making. Here, this paper proposes a control co-design approach, optimizing design and control decisions for integrating OWFs into the power grid while considering energy market and primary frequency market participation. Additionally, we introduce optimal sizing solutions for energy storage systems deployed onshore to enhance revenue for OWF developers over time. This framework addresses uncertainties related to wind resources and energy prices. We analyze five U.S. west-coast offshore wind farm locations and potential interconnection points, as identified by the Bureau of Ocean Energy Management (BOEM). Results show that optimized control co-design solutions can increase market revenue by 3.2% and provide flexibility in managing wind resource uncertainties.

Control Co-design↗
DOE OSTI2019

An adjoint method for neoclassical stellarator optimization

Stellarators remain an enticing route to steady-state fusion power. However, to achieve the required confinement, the magnetic geometry must be highly optimized. This optimization requires navigating high-dimensional spaces, often necessitating the use of gradient-based methods. The gradient of the neoclassical fluxes is expensive to compute with classical methods, requiring $O(N)$ flux computations, where $N$ is the number of parameters. To reduce the cost of the gradient computation, we present an adjoint method for computing the derivatives of moments of the neoclassical distribution function for stellarator optimization. The linear adjoint method allows derivatives of quantities which depend on solutions of a linear system, such as moments of the distribution function, to be computed with respect to many parameters from the solution of only two linear systems. This reduces the cost of computing the gradient to the point that the finite-collisionality neoclassical fluxes can be used within an optimization loop. With the neoclassical adjoint method, we compute solutions of the drift kinetic equation and an adjoint drift kinetic equation to obtain derivatives of neoclassical quantities with respect to geometric parameters. When the number of parameters in the derivative is large ( $O(10^{2})$ ), this adjoint method provides up to a factor of 200 reduction in cost. We demonstrate adjoint-based optimization of the field strength to obtain minimal bootstrap current on a surface. With adjoint-based derivatives, we also compute the local sensitivity to magnetic perturbations on a flux surface and identify regions where tight tolerances on error fields are required for control of the bootstrap current or radial transport. Moreover, the solve for the ambipolar electric field is accelerated using a Newton method with derivatives obtained from the adjoint method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Physics↗
DOE OSTI2024

Stellarator optimization with constraints

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a quasi-axisymmetric stellarator for favourable physics properties while preventing strong shaping of the plasma boundary, which can be difficult to create with external current sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multiple aspects maintenance ontology-based intelligent maintenance optimization framework for safety-critical systems

Abstract Maintenance optimization is a process for improving the efficiency of maintenance strategies and activities, considering various aspects of the target system and components, such as the probabilities of system failures and the cost of repair and replacement of a failed component. The improvement of maintenance optimization algorithms generally requires information from various data sources. For example, it may require the system risk information derived from risk analysis tools or the residual lifetime of a component from fault prognosis tools. The requirements of data acquisition (DAQ) and aggregation pose new challenges for maintenance management systems (MMSs) that implement and use these maintenance optimization algorithms. This paper proposes a multiple aspects maintenance ontology-based framework to facilitate DAQ from MMSs, online monitoring systems, fault detection and discrimination tools, risk assessment tools, decision-making tools, and component identification tools, and accelerate the implementation and verification of contemporary maintenance optimization models and algorithms. The proposed framework consists of a multi-aspect maintenance ontology with critical information for maintenance optimization and application interfaces for collecting information from various data sources, such as fault prognosis tools, online monitoring tools, risk assessment tools, and decision-making algorithms. In addition, this paper proposes a heuristic method for integrating concepts and properties from other existing ontologies into the proposed framework when the existing ontology is not fully compatible with the ontology under construction. Finally, the paper verifies the proposed ontology framework using a feedwater system designed for nuclear power plants with valves and filters as the components under maintenance.

Diao, Xiaoxu (ORCID:0000000346726352)↗
DOE OSTI2018

An Ensemble Optimization Framework for Coupled Design of Hydropower Contracts and Real-Time Reservoir Operating Rules

Often, revenues from hydropower generation depend on the operator's ability to provide firm power in the presence of uncertain inflows. The primary options available for optimizing revenue are negotiation of a firm power contract before operations begin and adjustment of the reservoir release during operations. Contract and release strategy optimization are closely coupled and most appropriately analyzed with stochastic real-time control methods. Here we use an ensemble-based approach to stochastic optimization that provides a convenient way to construct nonparametric revenue probability distributions to explore the implications of uncertainty. The firm power contract is a simplified bilateral fixed price agreement that partially insulates operator and buyer from price fluctuations. The release control laws and firm energy target are jointly optimized to maximize the operator's expected revenue. Revenue probability distributions and related spill performance statistics indicate that predictive operating strategies such as stochastic dynamic programming and model predictive control can give significantly better performance than standard deterministic operating rules. The performance obtained from batch optimization with perfect inflow information establishes a convenient upper bound on potential revenue and provides a baseline for assessing the significance of differences between real-time operating strategies. Sensitivity analysis indicates that the benefits of predictive operational strategies are greatest for reservoirs with medium nondimensional residence times and less important for reservoirs with large residence times. Overall, probabilistic analysis of the coupled hydropower contract-operations problem provides a realistic way to assess revenue and risk for reservoirs that must provide firm power when inflows are uncertain.

13 HYDRO ENERGY↗
DOE OSTI2022

Evaluation and Optimization of Snow Albedo Scheme in Noah-MP Land Surface Model Using In Situ Spectral Observations in the Colorado Rockies

The Biosphere-Atmosphere Transfer Scheme (BATS) ground snow albedo algorithm is commonly used in land-surface models (LSM), weather forecasting and research applications. This study addresses key uncertainties in BATS simulated ground snow albedo within the Noah-MP LSM framework through evaluation and optimization of the Noah-MP BATS ground snow albedo formulation using 2-band (visible and near-infrared (NIR)) in situ albedo observations at Rocky Mountain field stations. The Noah-MP BATS ground snow albedo scheme is extremely sensitive to its input parameters. Namely, an ensemble generated by varying BATS input parameters within potentially plausible ranges provides an average daily range (maximum ensemble member minus minimum ensemble member) of ground snow albedo exceeding 0.45 in visible and NIR bands. Parameter optimization improves agreement between simulated and in situ observed ground snow albedo in visible, NIR and broadband spectrums. Importantly, optimized parameters result in reduced biases relative to observed fresh-snow albedo and better agreement with observed albedo decay. Our analysis across different sites supports that the optimized BATS ground snow albedo parameters are appropriate to transfer in space and time, at least within the region studied (the central-southern Rocky Mountains). The primary error source remaining after parameter optimization is that observed fresh-snow albedo is highly variable, particularly in the NIR spectrum, whereas BATS fresh-snow albedo is constant, an issue which requires further investigation. This study shows significant correlations between observed fresh-snow albedo and surface meteorological conditions (e.g., downward shortwave radiation and temperature) which can support future model development that attempts to include a time-varying formulation for fresh-snow albedo.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Constrained quantum optimization for extractive summarization on a trapped-ion quantum computer

Abstract Realizing the potential of near-term quantum computers to solve industry-relevant constrained-optimization problems is a promising path to quantum advantage. In this work, we consider the extractive summarization constrained-optimization problem and demonstrate the largest-to-date execution of a quantum optimization algorithm that natively preserves constraints on quantum hardware. We report results with the Quantum Alternating Operator Ansatz algorithm with a Hamming-weight-preserving XY mixer (XY-QAOA) on trapped-ion quantum computer. We successfully execute XY-QAOA circuits that restrict the quantum evolution to the in-constraint subspace, using up to 20 qubits and a two-qubit gate depth of up to 159. We demonstrate the necessity of directly encoding the constraints into the quantum circuit by showing the trade-off between the in-constraint probability and the quality of the solution that is implicit if unconstrained quantum optimization methods are used. We show that this trade-off makes choosing good parameters difficult in general. We compare XY-QAOA to the Layer Variational Quantum Eigensolver algorithm, which has a highly expressive constant-depth circuit, and the Quantum Approximate Optimization Algorithm. We discuss the respective trade-offs of the algorithms and implications for their execution on near-term quantum hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Computational design of thermoelectric alloys through optimization of transport and dopability

Alloying is a common technique to optimize the functional properties of materials for thermoelectrics, photovoltaics, energy storage etc. Designing thermoelectric (TE) alloys is especially challenging because it is a multi-property optimization problem, where the properties that contribute to high TE performance are interdependent. In this work, we develop a computational framework that combines first-principles calculations with alloy and point defect modeling to identify alloy compositions that optimize the electronic, thermal, and defect properties. We apply this framework to design n-type Ba 2(1–x) Sr 2x CdP 2 Zintl thermoelectric alloys. Our predictions of the crystallographic properties such as lattice parameters and site disorder are validated with experiments. To optimize the conduction band electronic structure, we perform band unfolding to sketch the effective band structures of alloys and find a range of compositions that facilitate band convergence and minimize alloy scattering of electrons. Here, we assess the n-type dopability of the alloys by extending the standard approach for computing point defect energetics in ordered structures. Through the application of this framework, we identify an optimal alloy composition range with the desired electronic and thermal transport properties, and n-type dopability. Such a computational framework can also be used to design alloys for other functional applications beyond TE.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Stochastic Optimization with Risk Aversion for Virtual Power Plant Operations: A Rolling Horizon Control

While the concept of aggregating and controlling renewable distributed energy resources (DERs) to provide grid services is not new, increasing policy support of DER market participation has driven research and development in algorithms to pool DERs for economically viable market participation. Sandia National Laboratories recently undertook a three-year research program to create the components of a real-world virtual power plant (VPP) that can simultaneously participate in multiple markets. Our research extends current state-of-the-art rolling horizon control through the application of stochastic programming with risk aversion at various time resolutions. Our rolling horizon control consists of (1) day-ahead optimization to produce an hourly aggregate schedule for the VPP operator and (2) sub-hourly optimization for real-time dispatch of each VPP subresource. Both optimization routines leverage a two-stage stochastic program (SP) with risk aversion, and integrate the most up-to-date forecasts to generate probabilistic scenarios in real operating time. Our results demonstrate the benefits to the VPP operator of constructing a stochastic solution regardless of the weather. In more extreme weather, applying risk optimization strategies can dramatically increase the financial viability of the VPP. As a result, the methodologies presented here can be further tailored for optimal control of any VPP asset fleet and its operational requirements.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Accelerating the convergence of auxiliary-field quantum Monte Carlo in solids with optimized Gaussian basis sets

We investigate the use of optimized correlation-consistent Gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second-order Møller–Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely, calculations in the Kohn–Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn–Sham basis and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si, and MgO and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for the systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.

Morales, Miguel A. (ORCID:0000000263893067)↗
DOE OSTI2024

A multifidelity Bayesian optimization method for inertial confinement fusion design

Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and effectively explore the numerous parameters required to describe a typical experiment. Therefore, traditionally, ICF design has relied on low-fidelity modeling to initially identify potentially interesting design regions, which are then subsequently explored via selected high-fidelity modeling. In this paper, we demonstrate that this two-step approach can be insufficient: even for simple design problems, a two-step optimization strategy can lead high-fidelity searching toward incorrect regions and consequently waste computational resources on parameter regimes far away from the true optimal solution. We reveal that a primary cause of this behavior in ICF design problems is the presence of low-fidelity optima in different regions of the parameter space far away from high-fidelity optima. To address this issue, we propose an iterative multifidelity Bayesian optimization method based on Gaussian Process Regression that leverages both low- and high-fidelity models simultaneously. We demonstrate, using both two- and eight-dimensional ICF test problems, that our algorithm can effectively utilize both low-fidelity and high-fidelity models to refine the designs. This approach proves to be more efficient than relying solely on high-fidelity modeling for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗