SEARCH · Engineering Papers

Results for “Simulation 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Bayesian optimization for inverse calibration of expensive computer models: A case study for Johnson-Cook model in machining

Inverse model calibration for identifying the constitutive model parameters can be computationally demanding for expensive-to-evaluate simulation models. Here, this paper presents a modified Bayesian optimization (BO) method, denoted as BO-bound, that incorporates theoretical bounds on the quantity of interest. A case study for the inverse calibration of the Johnson Cook (J-C) flow stress model parameters is presented using machining (cutting) force data. The results show fast calibration of the five J-C parameters within 25 simulations. In general, the BO-bound method is applicable for inverse calibration of any expensive simulation models as well as optimization problems with known bounds.

Bayesian optimization↗

DOE OSTI2019

Bi-fidelity approximation for uncertainty quantification and sensitivity analysis of irradiated particle-laden turbulence

Particle-laden turbulent flows subject to radiative heating are relevant in many applications, for example concentrated solar power receivers. Efficient and accurate simulations provide valuable insights and enable optimization of such systems. However, as there are many uncertainties inherent in such flows, uncertainty quantification is fundamental to improve the predictive capabilities of the numerical simulations. For large-scale, multi-physics problems exhibiting high-dimensional uncertainty, characterizing the stochastic solution presents a significant computational challenge as most strategies require a large number of high-fidelity solves. This requirement might result in an infeasible number of simulations when a typical converged high-fidelity simulation requires intensive computational resources. To reduce the cost of quantifying high-dimensional uncertainties, we investigate the application of a non-intrusive, bi-fidelity approximation to estimate statistics of quantities of interest associated with an irradiated particle-laden turbulent flow. This method exploits the low-rank structure of the solution to accelerate the stochastic sampling and approximation processes by means of cheaper-to-run, lower fidelity representations. Here, the application of this bi-fidelity approximation results in accurate estimates of the quantities of interest statistics, while requiring a small number of high-fidelity model evaluations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Optimising response to an introduction of African swine fever in wild pigs

African swine fever virus (ASFv) is a virulent pathogen that threatens domestic swine industries globally and persists in wild boar populations in some countries. Persistence in wild boar can challenge elimination and prevent disease-free status, making it necessary to address wild swine in proactive response plans. In the United States, invasive wild pigs are abundant and found across a wide range of ecological conditions that could drive different epidemiological dynamics among populations. Information on the size of the control areas required to rapidly eliminate the ASFv in wild pigs and how this area should change with management constraints and local ecology is needed to optimize response planning. We developed a spatially explicit disease transmission model contrasting wild pig movement and contact ecology in two ecosystems in Southeastern United States. We simulated ASFv spread and determined the optimal response area (reported as the radius of a circle) for eliminating ASFv rapidly over a range of detection times (when ASFv was detected relative to the true date of introduction), culling capacities (proportion of wild pigs in the culling zone removed weekly) and wild pig densities. Large radii for response areas (14 km) were needed under most conditions but could be shortened with early detection (≤ 8 weeks) and high culling capacities (≥ 15% weekly). Under most conditions, the ASFv was eliminated in less than 22 weeks using optimal control radii, although ecological conditions with high rates of wild pig movement required higher culling capacities (≥ 10% weekly) for elimination within 1 year. The results highlight the importance of adjusting response plans based on local ecology and show that wild pig movement is a better predictor of the optimal response area than the number of ASFv cases early in the outbreak trajectory. Furthermore, our framework provides a tool for determining optimal control plans in different areas, guiding expectations of response impacts, and planning resources needed for rapid elimination.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Multi-scale simulations of wind farm performance with complex terrain and weather events

Predicting the response of wind farms to local weather is necessary for optimal design and operation. In this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine wakes, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with the GAD at10 m resolution. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed turbine wake structure, velocity deficit, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

DOE OSTI2019

Linear Regression Based Multi-fidelity Surrogate for Disturbance Amplification in Multi-phase Explosion

When simulations are very expensive and many are required, as for optimization or uncertainty quantification, a way to reduce cost is using surrogates. With multiple simulations to predict the quantity of interest, some being very expensive and accurate (high-fidelity simulations) and others cheaper but less accurate (low-fidelity simulations), it may be worthwhile to use multifidelity surrogates (MFSs). Moreover, if we can afford just a few high-fidelity simulations or experiments, MFS becomes necessary. Co-Kriging, which is probably the most popular MFS, replaces both low-fidelity and high-fidelity simulations by a single MFS. A recently proposed linear regression–based MFS (LR-MFS) offers the option to correct the LF simulations instead of correcting the LF surrogate in the MFS. When the low-fidelity simulation is cheap enough for use in an application, such as optimization, this may be an attractive option. Here, we explore the performance of LR-MFS using exact and surrogate-replaced low-fidelity simulations. The concern studied is a cylindrical dispersal of 100-μ m-diameter solid particles after detonation and the quantity of interest is a measure of the amplification of the departure from axisymmetry. We find very substantial accuracy improvements for this problem using the LR-MFS with exact low-fidelity simulations. Inspired by these results, we also compare the performance of co-Kriging to the use of Kriging to correct exact low-fidelity simulations and find a similar accuracy improvement when simulations are directly used. For this problem, further improvements in accuracy are achievable by taking advantage of inherent parametric symmetries. These findings may alert users of MFSs to the possible advantages of using exact low-fidelity simulations when this is affordable.

42 ENGINEERING↗

DOE OSTI2021

Development of high B s Fe–Ni-based metal amorphous nanocomposite by optimization of glass-forming ability

Thermocalc simulations identified compositions with good glass-forming ability (GFA) by locating minima in liquidus temperature, and solidification range in (Fe 70 Ni 30 ) x (B-Si-Nb) 100-x alloys with x = 82% and x = 85%, increased compared to previously developed x = 80% alloys. 3 compositions in the x = 82% system and 1 in the x = 85% system were successfully cast as amorphous magnetic ribbon. Magnetic properties of as-cast alloys showed improvements in Curie temperature and saturation induction, especially for the 85% alloy. Crystallization behavior in the 85% alloy was determined to occur as a 2-step process. Annealing studies showed (Fe 70 Ni 30 ) 82 Nb 2 Si 0 B 16 and (Fe 70 Ni 30 ) 82 Nb 3 Si 0 B 15 alloys to exhibit low coercivity and high B s with annealing. TEM results show significant crystallization in as-cast ribbon. The (Fe 70 Ni 30 ) 85 Nb 0.5 Si 0 B 14.5 alloy as-cast has small, uniformly distributed grains, allowing a useful nanocrystalline alloy to be produced as cast. Lastly, the parameter ΔT xg , the temperature range in which the material can be thermomechanically processed, was in the range of 12–23 °C.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Implicit neural representations for experimental steering of advanced experiments

Scattering measurements using electrons, neutrons, or photons are essential for obtaining microscopic insights into materials. However, limited facility availability and high-dimensional scattering data necessitate more efficient experimental steering techniques. Here, we report a machine learning method that guides scattering data collection and facilitates real-time estimation of model parameters, given a reliable forward model to simulate experimental signals. We employ implicit neural representations as efficient surrogates that link model parameters with simulated spectroscopies. This enables a Bayesian optimal experimental design framework to estimate the probability distributions of parameters from high-dimensional scattering data. We demonstrate the proposed method using inelastic neutron scattering with simulated and real experimental data, highlighting the method’s ability to provide real-time parameter estimation with quantified uncertainties and to deliver informed experimental guidance that reduces experimental time while maximizing scientific output. This approach paves the way for accelerated discoveries in condensed matter through scattering measurements.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Ultrafast Dynamics of Photoexcited Charge Carriers in In{sub 0.53}Ga{sub 0.47}As/In{sub 0.52}Al{sub 0.48}As Superlattices under Femtosecond Laser Excitation

The results of experimental studies of the time dynamics of photoexcited charge carriers in In{sub 0.53}Ga{sub 0.47}As/In{sub 0.52}Al{sub 0.48}As superlattices grown by molecular-beam epitaxy on a GaAs substrate with a metamorphic buffer are reported. On the basis of the results of the numerical simulation of band diagrams, the optimal thickness of the In{sub 0.52}Al{sub 0.48}As barrier layer (4 nm) is chosen. At this thickness, the electron wave functions in In{sub 0.53}Ga{sub 0.47}As substantially overlap the In{sub 0.52}Al{sub 0.48}As barriers. This makes it possible to attain a short lifetime of photoexcited charge carriers (τ ~ 3.4 ps) at the wavelength λ = 800 nm and the pumping power 50 mW without doping of the In{sub 0.53}Ga{sub 0.47}As layer with beryllium. It is shown that an increase in the wavelength to λ = 930 nm (at the same pumping power) yields a decrease in the lifetime of photoexcited charge carriers to τ ~ 2 ps. This effect is attributed to an increase in the capture cross section of trapping states for electrons with lower energies and to a decrease in the occupancy of traps at lower excitation densities.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Active Learning for Metamaterial Optimization on HPC and QC Integrated Systems

Active learning algorithms, integrating machine learning, quantum computing and optics simulation in an iterative loop, offer a promising approach to optimizing metamaterials. However, these algorithms can face difficulties in optimizing highly complex structures due to computational limitations. High-performance computing (HPC) and quantum computing (QC) integrated systems can address these issues by enabling parallel computing. In this study, we develop an active learning algorithm working on HPC-QC integrated systems. We evaluate the performance of optimization processes within active learning (i.e., training a machine learning model, problem-solving with quantum computing, and evaluating optical properties through wave-optics simulation) for highly complex metamaterial cases. Our results showcase that utilizing multiple cores on the integrated system can significantly reduce computational time, thereby enhancing the efficiency of optimization processes. Therefore, we expect that leveraging HPC-QC integrated systems helps effectively tackle large-scale optimization challenges in general.

Kim, Seongmin↗

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 OSTI2025

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗

DOE OSTI2018

Cogging torque minimization in transverse flux machines

This study presents the design considerations in cogging torque minimization in two types of transverse flux machines. The machines have a double stator-single rotor configuration with flux concentrating ferrite magnets. One of the machines has pole windings across each leg of an E-Core stator. Another machine has quasi-U-shaped stator cores and a ring winding. The flux in the stator back iron is transverse in both machines. Different methods of cogging torque minimization are investigated. Key methods of cogging torque minimization are identified and used as design variables for optimization using a design of experiments (DOE) based on the Taguchi method. A multi-level DOE is proposed as an optimization method to reach an optimum solution with minimum simulations. The case study is analyzed in a two level DOE optimization. Finite element analysis is used to study the different effects. Two prototypes are fabricated for validating the FEA results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Improved optical standing-wave beam splitters for dilute Bose–Einstein condensates

Bose–Einstein condensate (BEC)-based atom interferometry exploits low temperatures and long coherence lengths to facilitate high-precision measurements. Progress in atom interferometry promises improvements in navigational devices like gyroscopes and accelerometers, as well as applications in fundamental physics such as accurate determination of physical constants. Previous work demonstrates that beam splitters and mirrors for coherent manipulation of dilute BEC momentum in atom interferometers can be implemented with sequences of non-resonant standing-wave light pulses. While previous work focuses on the optimization of the optical pulses’ amplitude and duration to produce high-order momentum states with high fidelity, we explore how varying the shape of the optical pulses affects optimal beam-splitter performance, as well as the effect of pulse shape on the sensitivity of optimized parameters in achieving high fidelity in high-momentum states. Here, in simulations of two-pulse beam splitters utilizing optimized square, triangle, and sinc-squared pulse shapes applied to dilute BECs, we, in some cases, reduce parameter sensitivity by an order of magnitude while maintaining fidelity.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Beyond Price Taker: Conceptual Design and Optimization of Integrated Energy Systems Using Machine Learning Market Surrogates

Future electricity generation systems must be optimized to provide flexibility that counteracts the variability of non-dispatchable renewable energy sources and ensures the reliability and safety of critical infrastructure, including the electric grid. The current state-of-the-art is to co-optimize the design and operation of integrated energy systems (IES) treating historical or predicted time-series electricity prices as fixed parameters. Recent literature has shown the limitations of this price taker assumption, which neglects how IES optimization decisions influence market outcomes. As such, this paper proposes a new optimization formulation that uses machine learning surrogate models, trained from a library of annual market operation simulations, to embed IES market interactions into the co-optimization problem directly. Using a thermal generator example built in the open-source IDAES computational environment, we show that the price taker approach routinely over-predicts annual revenues by 8% or more compared to a validation simulation, where the proposed approach has a typical relative error of 1% or less.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2019

Evaluation of genomic selection and marker-assisted selection in Miscanthus and energycane

Abstract Although energycane ( Saccharum spp. hybrids) is widely used as a source of lignocellulosic biomass for bioethanol, breeding this crop for disease resistance is challenging due to its narrow genetic base. Therefore, efforts are underway to introgress novel sources of genetic resistance from Miscanthus into energycane. Given that disease resistance in energycane could be either qualitative or quantitative in nature, careful examination of a wide variety of genomic-enabled breeding approaches will be crucial to the success of such an undertaking. Here we examined the efficiency of both genomic selection (GS) and marker-assisted selection (MAS) for traits simulated under different genetic architectures in F 1 and BC 1 populations of Miscanthus × Miscanthus and sugarcane × sugarcane crosses. We observed that the performance of MAS was comparable and sometimes superior to GS for traits simulated with four quantitative trait nucleotides (QTNs). In contrast, as the number of simulated QTN increased, all four GS models that were evaluated tended to outperform MAS, select more phenotypically optimal F 1 individuals, and accurately predict simulated trait values in subsequent BC 1 generations. We therefore conclude that GS is preferable to MAS for introgressing genetic sources of horizontal disease resistance from Miscanthus to energycane, while MAS remains a suitable option for introgressing vertical disease resistance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Concurrent multi-peak Bragg coherent x-ray diffraction imaging of 3D nanocrystal lattice displacement via global optimization

Abstract In this paper we demonstrated a method to reconstruct vector-valued lattice distortion fields within nanoscale crystals by optimization of a forward model of multi-reflection Bragg coherent diffraction imaging (MR-BCDI) data. The method flexibly accounts for geometric factors that arise when making BCDI measurements, is amenable to efficient inversion with modern optimization toolkits, and allows for globally constraining a single image reconstruction to multiple Bragg peak measurements. This is enabled by a forward model that emulates the multiple Bragg peaks of a MR-BCDI experiment from a single estimate of the 3D crystal sample. We present this forward model, we implement it within the stochastic gradient descent optimization framework, and we demonstrate it with simulated and experimental data of nanocrystals with inhomogeneous internal lattice displacement. We find that utilizing a global optimization approach to MR-BCDI affords a reliable path to convergence of data which is otherwise challenging to reconstruct.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Linear Approximation Line Pack Model for Integrated Electricity and Natural Gas Systems OPF

With the increasing installation of renewable energy electricity generation capacity, the penetration of natural gas fired generation is also substantailly increasing in power systems due to their fast response capability to mitigate system variability and uncertainty. The dependencies between the electricity system and the natural gas system are becoming more important to maintain the reliability of both systems. In integrated electricity and natural gas systems (IEGSs) operations, the nonlinearity of the natural gas system brings challenges to the system operators. In this paper, the IEGSs optimal power flow (OPF) model is proposed that includes the line pack of the gas system via linear approximation. The simulation results demonstrate that the proposed model can obtain optimal solutions efficiently that are comparable in accuracy as the nonlinear line pack model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗