Engineering Papers⌕ Search

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.

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Xyce™ Parallel Electronic Simulator Users’ Guide, Version 7.10

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

MCOR User Guide

User guide for the MCOR software package which is currently publicly hosted on Github (https://github.com/pnnl/MCOR). The Microgrid Component Optimization for Resilience (MCOR) tool simulates the operation of a renewable energy, battery, and back-up generator microgrid under a large range of outage conditions to understand how a potential system can meet the resilience goals of a particular site. It is an open-source, command line, Python-based tool that produces an output Excel spreadsheet as well as several types of plots to enable a user to compare different microgrid system sizes and costs. It is intended for high-level system planning and opportunity identification, and not for detailed electric system modeling and design. The tool includes a range of input parameters that can be adjusted or tuned to provide a more custom analysis as needed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present the development of a machine learning (ML) based regulation system for third-order resonant beam extraction in the Mu2e experiment at Fermilab. Classical and ML-based controllers have been optimized using semi-analytic simulations and evaluated in terms of regulation performance and training efficiency. We compare several controller architectures and discuss the integration of neural control into an adaptive framework. We also present progress on surrogate models that predict the controller response given a spill intensity and controller action history. To enable real-time deployment, we report progress on implementing low-latency, edge-based inference suitable for hardware-constrained environments. Our results demonstrate the feasibility and advantages of ML-based control in managing complex, time-varying physical systems, with broader implications for accelerator operations and other domains requiring fast, adaptive regulation.

Berlioz, Jose Rene [Fermilab]↗
DOE OSTI2020

Crop Modeling Application to Improve Irrigation Efficiency in Year-Round Vegetable Production in the Texas Winter Garden Region

Given a rising demand for quality assurance, rather than solely yield, supplemental irrigation plays an important role to ensure the viability and profitability of vegetable crops from unpredictable changes in weather. However, under drought conditions, agricultural irrigation is often given low priority for water allocation. This reduced water availability for agriculture calls for techniques with greater irrigation efficiency, that do not compromise crop quality and yield, and that provide economic benefit for producers. This study developed vegetable growing models for eight different vegetable crops (bush bean, green bean, cabbage, peppermint, spearmint, yellow straight neck squash, zucchini, and bell pepper) based on data from several years of field research. The ALMANAC model accurately simulated yields and water use efficiency (WUE) of all eight vegetables. The developed vegetable models were used to evaluate the effects of various irrigation regimes on vegetable growth and production in several locations in the Winter Garden Region of Texas, under variable weather conditions. Based on our simulation results from 960 scenarios, optimal irrigation amounts that produce high yield as well as reasonable economic profit to producers were determined for each vegetable crop. Overall, yields for all vegetables increased as irrigation amounts increased. However, irrigation amounts did not have a sustainable impact on vegetable yield at high irrigation treatments, and the WUEs of most vegetables were not significantly different among various irrigation regimes. When vegetable yields were compared with water cost, the rate decreased as irrigation amounts increased. Thus, producers will not receive economic benefits when vegetable irrigation water demand is too high.

Kim, Sumin↗
DOE OSTI2023

Ducted Fuel Injection Provides Consistently Lower Soot Emissions in Sweep to Full-Load Conditions

Earlier studies have proven how ducted fuel injection (DFI) substantially reduces soot for low- and mid-load conditions in heavy-duty engines, without significant adverse effects on other emissions. Nevertheless, no comprehensive DFI study exists showing soot reductions at high- and full-load conditions. This study investigated DFI in a single-cylinder, 1.7-L, optical engine from low- to full-load conditions with a low-net-carbon fuel consisting of 80% renewable diesel and 20% biodiesel. Over the tested load range, DFI reduced engine-out soot by 38.1–63.1% compared to conventional diesel combustion (CDC). This soot reduction occurred without significant detrimental effects on other emission types. Thus, DFI reduced the severity of the soot–NO x tradeoff at all tested conditions. While DFI delivered considerable soot reductions in thepresent study, previous DFI studies at low- and mid-load conditions delivered larger soot reductions (>90%) compared to CDC operation at the same conditions. Therefore, the DFI configuration used here has been deemed nonoptimal (in terms of parameters such as the injector-spray and piston geometries), and several improvements are recommended for future studies with high-load DFI. Further, these improvements include employing better spray-duct alignment, a deeper piston bowl with a smaller injector umbrella angle, and a fuel injector that opens and closes faster. The study also suggests future research to make DFI ready for commercialization, such as metal-engine tests to ensure desirable DFI performance over an engine’s complete speed/load map. Overall, this study supports the continued development and commercialization of DFI to meet upcoming emissions regulations for heavy-duty vehicles. Specifically, multicylinder engine experiments and CFD simulations should be utilized to optimize the performance and clarify the full potential of DFI.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

A Spectral Approach for the Design of Experiments: Design, Analysis and Algorithms

This paper proposes a new approach to construct high quality space-filling sample designs. First, we propose a novel technique to quantify the space-filling property and optimally trade-off uniformity and randomness in sample designs in arbitrary dimensions. Second, we connect the proposed metric (defined in the spatial domain) to the quality metric of the design performance (defined in the spectral domain). This connection serves as an analytic framework for evaluating the qualitative properties of space-filling designs in general. Using the theoretical insights provided by this spatial-spectral analysis, we derive the notion of optimal space-filling designs, which we refer to as space-filling spectral designs. Third, we propose an efficient estimator to evaluate the space-filling properties of sample designs in arbitrary dimensions and use it to develop an optimization framework for generating high quality space-filling designs. Finally, we carry out a detailed performance comparison on two different applications in varying dimensions: a) image reconstruction and b) surrogate modeling for several benchmark optimization functions and a physics simulation code for inertial confinement fusion (ICF). Our results clearly evidence the superiority of the proposed space-filling designs over existing approaches, particularly in high dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantum computing hardware for HEP algorithms and sensing

Quantum information science harnesses the principles of quantum mechanics to realize computational algorithms with complexities vastly intractable by current computer platforms. Typical applications range from quantum chemistry to optimization problems and also include simulations for high energy physics. The recent maturing of quantum hardware has triggered preliminary explorations by several institutions (including Fermilab) of quantum hardware capable of demonstrating quantum advantage in multiple domains, from quantum computing to communications, to sensing. The Superconducting Quantum Materials and Systems (SQMS) Center, led by Fermilab, is dedicated to providing breakthroughs in quantum computing and sensing, mediating quantum engineering and HEP based material science. The main goal of the Center is to deploy quantum systems with superior performance tailored to the algorithms used in high energy physics. In this Snowmass paper, we discuss the two most promising superconducting quantum architectures for HEP algorithms, i.e. three-level systems (qutrits) supported by transmon devices coupled to planar devices and multi-level systems (qudits with arbitrary N energy levels) supported by superconducting 3D cavities. For each architecture, we demonstrate exemplary HEP algorithms and identify the current challenges, ongoing work and future opportunities. Furthermore, we discuss the prospects and complexities of interconnecting the different architectures and individual computational nodes. Finally, we review several different strategies of error protection and correction and discuss their potential to improve the performance of the two architectures. This whitepaper seeks to reach out to the HEP community and drive progress in both HEP research and QIS hardware.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗
DOE OSTI2020

Methodology and Application of Physical Security Effectiveness Based on Dynamic Force-on-Force Modeling

This report describes the research and development being performed at INL towards a dynamic modeling and simulation framework to enable physical security optimization at commercial nuclear power plants. The framework is based on the dynamic modeling tool EMRALD and is demonstrated for applications that can result in physical security optimization. Two main applications are presented: 1. Integrating FLEX portable equipment performance with FOF models of a plant’s physical security posture, and 2. Location optimization of bullet resistant enclosure. The generic framework for modeling FLEX portable equipment is described in detail, followed by a case study modeling an adversarial attack aimed at causing a radiological release by sabotaging the plant’s power supply and its ultimate heat sink capabilities at a hypothetical PWR. Two distinct FLEX deployment strategies, series and parallel, are modeled with distinct timelines. The results of the adversarial attack modeled in a commercial FOF tool, AVERT, are integrated with the FLEX deployment model in EMRALD. Monte Carlo simulation is used to model the distribution of the timeline in FLEX deployment strategies. Thermal-hydraulic analysis of FLEX performance is performed in RELAP5 and integrated with the EMRALD simulations to provide more realistic timelines in the models. The results demonstrate that, even in the extreme case of a successful adversarial attack, deployment of FLEX equipment can result in a significantly high likelihood of preventing radiological release. The modeling and simulation framework of integrating FLEX equipment with FOF models enables the NPPs to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security. The objective of location optimization of BRE is to determine the best location in the plant for a new BRE being planned by the plant to enhance their physical security effectiveness. The plant physical security FOF model is integrated with EMRALD that performs Monte Carlo simulation to run different attack scenarios and a discrete set of potential BRE locations. Sensitivity analysis is used to determine the most effective location for the BRE. The optimization approach can be extended to wide applications such as location optimization of remotely operated weapons and other strategic fixed assets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗
DOE OSTI2019

Calibration of Car-Following Models Using a Dual Genetic Algorithm with Central Composite Design

We report the calibration of microscopic traffic simulation models has received much attention in the simulation field. Although no standard has been established for it, a genetic algorithm (GA) has been widely employed in recent literature because of its high efficiency to find solutions in such optimization problems. However, the performance still falls short in simulation analyses to support fast decision making. This paper proposes a new calibration procedure using a dual GA and central composite design (CCD) in order to improve the efficiency. The calibration exercise goes through three major sequential steps: (1) experimental design using CCD for a quadratic response surface model (RSM) estimation, (2) 1st GA procedure using the RSM with CCD to find a near-optimal initial population for a next step, and (3) 2nd GA procedure to find a final solution. The proposed method was applied in calibrating the Gipps car-following model with respect to maximizing the likelihood of a spacing distribution between a lead and following vehicle. In order to evaluate the performance of the proposed method, a conventional calibration approach using a single GA was compared under both simulated and real vehicle trajectory data. It was found that the proposed approach enhances the optimization speed by starting to search from an initial population that is closer to the optimum than that of the other approach. This result implies the proposed approach has benefits for a large-scale traffic network simulation analysis. Finally, this method can be extended to other optimization tasks using GA in transportation studies.

42 ENGINEERING↗
DOE OSTI2022

Machine-learning-assisted high-temperature reservoir thermal energy storage optimization

High-temperature reservoir thermal energy storage (HT-RTES) has the potential to become an indispensable component in achieving the goal of the net-zero carbon economy, given its capability to balance the intermittent nature of renewable energy generation. In this study, a machine-learning-assisted computational framework is presented to co-optimize the performance metrics of HT-RTES by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Further, neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Non-trivial symmetries in quantum landscapes and their resilience to quantum noise

Very little is known about the cost landscape for parametrized Quantum Circuits (PQCs). Nevertheless, PQCs are employed in Quantum Neural Networks and Variational Quantum Algorithms, which may allow for near-term quantum advantage. Such applications require good optimizers to train PQCs. Recent works have focused on quantum-aware optimizers specifically tailored for PQCs. However, ignorance of the cost landscape could hinder progress towards such optimizers. In this work, we analytically prove two results for PQCs: (1) We find an exponentially large symmetry in PQCs, yielding an exponentially large degeneracy of the minima in the cost landscape. Alternatively, this can be cast as an exponential reduction in the volume of relevant hyperparameter space. (2) We study the resilience of the symmetries under noise, and show that while it is conserved under unital noise, non-unital channels can break these symmetries and lift the degeneracy of minima, leading to multiple new local minima. Based on these results, we introduce an optimization method called Symmetry-based Minima Hopping (SYMH), which exploits the underlying symmetries in PQCs. Our numerical simulations show that SYMH improves the overall optimizer performance in the presence of non-unital noise at a level comparable to current hardware. Overall, this work derives large-scale circuit symmetries from local gate transformations, and uses them to construct a noise-aware optimization method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Production Through Simulation: Using Simulation Technologies to Create and Evaluate Nuclear Fuel Fabrication Facility Designs

This project explores how physics-based and discrete-event simulation technologies can be jointly employed to model a suite of nuclear fuel fabrication facility designs with variable radiological environments, optimize their design to lower operational costs subject to a set of design constraints, and evaluate their performance. Initially, three variants of U-20Pu-10Zr metallic nuclear fuel in three specific geometric configurations are modeled and studied across more than 50 physics simulations in nuclear physics software packages SCALE and Monte Carlo N-Particle Transport (MCNP6.2). From these tests, values for effective dose rate and criticality are recorded for each of the nine alloy-geometry configurations. These values are then incorporated into discrete-event models of fuel fabrication in the simulation program ExtendSim Pro 10 as process attributes. This allows one to track the dose to facility personnel as fuel is fabricated, dually enabling us to evaluate the safety of a given facility design in terms of absorbed dose per worker per year and to craft operational guidelines that adhere to federal and local safety regulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Controlling Ligand Excimer Formation with Dipole Changes in Emissive Rare-Earth/Phosphonic Acid Complexes

The interactions between substituted arylvinyl phosphonic acid (AVPA) ligands within a Eu-AVPA complex are shown to influence the outcomes of excited state evolution after photoexcitation. Compared with unfunctionalized AVPAs, pairs of ligands functionalized with CF 3 in the para position preassociate in the ground state of complexes with Eu 3+ according to calculated geometry optimizations. The CF 3 -substituted AVPA complexes show evidence of red-shifted optical absorption and undergo more efficient excimer formation, as revealed by transient absorption spectroscopy. We rationalize this behavior through simulations of excited-state geometry optimizations that reveal evolution toward interligand phenyl-phenyl planarity for specific excited states. Emission from complexed Eu 3+ after energy transfer from the ligand is found to be weaker with CF 3 substitution, which we hypothesize is due to intracomplex, interligand aggregates with excimer-promoting geometries. These observations point to the need to consider ground-state geometries as well as dynamic excited-state processes to understand the flow of energy in rare earth coordination complexes.

09 BIOMASS FUELS↗
DOE OSTI2018

Initial Studies into Longitudinal Ionization Cooling for the Muon g-2 Experiment

Fermilab's Muon g-2 experiment aims to measure the anomalous magnetic moment of the muon to an unprecedented precision of 140 ppb. It relies on large numbers of muons surviving many turns in the storage ring without colliding with the sides, at least long enough for the muons to decay. Longitudinal ionization cooling is introduced with respect to Fermilab's Muon g-2 experiment in an attempt to increase storage and through this the statistics and quality of results. The ionization cooling is introduced to the beam through a material wedge, an initial simulation study is made into the positioning, material, and geometrical parameters of this wedge using G4Beamline. Results suggest a significant increase of 20 - 30% in the number of stored muons when the optimal wedge is included in the simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Extraction of Rare Earth Elements from Chloride Media with Tetrabutyl Diglycolamide in 1-Octanol Modified Carbon Dioxide

Rare earth elements (REEs) are critical to our modern world. Recycling REEs from used products could help with potential supply issues. Extracting REEs from chloride media with tetrabutyl diglycolamide (TBDGA) in carbon dioxide could help recycle REEs with less waste than traditional solvents. Carbon dioxide as a solvent is inexpensive, inert, and reusable. Conditions for extraction of Eu from aqueous chloride media were optimized by varying moles percent of 1-octanol modifier, temperature, pressure, Eu concentration, TBDGA concentration, Cl– concentration, and HCl concentration. These optimized conditions were tested on a Y, Ce, Eu, Tb simulant material, REEs containing NdFeB magnets, and lighting phosphor material. The optimized conditions were found to be 23 °C, 24.1 MPa, 0.5 mol% 1-octanol, with an excess of TBDGA. At these conditions 95 ± 2% Eu was extracted from 8 M (mol/m3) HCl. Extraction from the mixed REE simulate material resulted in separation of Y, Eu, and Tb from the Ce which remained in the aqueous solution. The extraction on NdFeB magnet dissolved into 8 M HCl resulted in extraction of Pr, Nd, Dy, and Fe >97%. This ends in a separation from B, Al, and Ni. Extraction from a trichromatic lighting phosphor leachate resulted in extraction of Y and Eu >93% and no extraction of Ba, Mg, and Al.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗