Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 595 records · Page 33

Low-Rate Cycling Regeneration Sustains NMC811–Graphite Battery Performance through Extreme Ultrafast Charging

Urban air mobility (UAM) systems, such as electric vertical takeoff and landing (eVTOL) aircraft, require batteries capable of extreme ultrafast charging to support high-frequency operations. Here, in this work, we evaluate the electrochemical stability of NMC811–graphite full cells using a complex mixed-modal protocol designed to simulate the aggressive turnaround demands of UAM missions. We demonstrate that incorporating a periodic low-rate (0.3C) “regeneration” cycle enables cells to maintain 96% capacity retention over 900 cycles, despite repeated exposure to ultrafast pulses. To validate the efficacy of this regeneration mechanism, we performed simplified control experiments comparing continuous 10C charging against a 10C–0.3C interleaved protocol. Advanced characterization reveals that while mechanical particle fracture is pervasive at these rates, the regeneration cycles decelerate the electrochemical degradation by resolving spatially heterogeneous Ni oxidation states and relaxing lithium-ion concentration gradients. These findings provide a vital mechanistic framework for battery operation, proving that restorative charging can effectively delay and mitigate failure in demanding UAM applications.

batteries↗

A scaling law of the neutral penetration length and Balmer- α wing shape in high-temperature plasmas

Hydrogen atoms penetrating deep inside high-temperature magnetically confined plasmas by repetitive charge-exchange collisions result in a particle source of the plasma, which affects the plasma transport significantly. In this paper, we present an approximate solution of the fluid equations for neutral transport and an analytical representation of the neutral penetration length, in a simplified plasma geometry. This analysis predicts a power-law decay in the Balmer-α line wings which reflects the velocity distribution of the neutral atoms, with the power-law exponent analytically represented as well. These scaling laws are compared with a simple Monte–Carlo simulation and spectroscopic observations of Large Helical Device plasmas. Since the Balmer-α line wings are experimentally accessible, our formulation opens the possibility to quickly estimate the neutral penetration length from spectroscopic observations.

neutral opacity↗

Novel Results Visualization for Dynamic PSA and New Modeling Features in EMRALD

The Event Modeling Risk Assessment Linked Diagram (EMRALD) tool, developed at the Idaho National Laboratory (INL), was designed to simplify the creation of dynamic models and support various research projects. One of the primary goals of EMRALD was to provide visual methods for modeling. EMRALD consists of two main components: a web-based user interface for model development and a solve engine for running model simulations. Over time, it has evolved to meet the diverse needs of its users. Initially, EMRALD's results were simple text outputs with final key state percentages and uncertainty bounds. However, because EMRALD utilizes a three-phase discrete event simulation and tracks the paths of each simulation run leading to a key state, there is significant potential to analyze large sets of path results data, including state paths, events, and timing. Visualizing this data meaningfully posed a challenge. To address this, a novel time-based Sankey diagram was developed. EMRALD exports results data in a format that can be opened in this Sankey viewer, allowing users to visualize paths, occurrences, events, and probability data for the entire simulation run in a single diagram. Moreover, when EMRALD was first created, there were limited tools capable of meeting its graphical requirements, many of which are no longer supported. In 2024, a new web-based interface was developed using modern graphing tools, enabling additional modeling features. This paper discusses the new dynamic PSA results visualization capability and the enhanced modeling tools available in EMRALD.

97 - MATHEMATICS AND COMPUTING↗

GeoThermalCloud for EGS – An Open-source, User-friendly, Scalable AI Workflow for Modeling Enhanced Geothermal Systems

Enhanced Geothermal Systems (EGS) offer a vast potential to expand the use of geothermal energy. Heat is extracted from this engineered system by injecting relatively cold water into subsurface fractures, which are in contact with hot dry rock, and brought back to surface through production wells. Creating EGS requires improving the natural permeability of hot crystalline rocks. In this short conference paper, we present a reproducible workflow for modeling EGS. Our workflow called the GeoThermalCloud (GTC) for EGS, leverages recent advances in machine learning, deep learning, and high-performance computing. This GTC framework is currently being made open-source, user-friendly, and reproducible through python scripts as well as Google Colab/Jupyter Notebooks. This GTC for EGS modeling scripts are made available at https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS and will constantly be updated to cater for geothermal community. Current GTC framework provides scripts to train deep learning (DL) models for techno-economics and data worth analysis. The Geothermal Design Tool (https://github.com/GeoDesignTool/GeoDT.git), a fast and simplified multi-physics solver, is used to develop a database for training DL models. This short paper provides details on the scripts to curate, process, and train DL models. The scripts can easily be modified to train on databases generated by other popular open-source simulators such as PFLOTRAN, STOMP, TOUGH, and GEOSX or commercial software such as ResFrac and COMSOL.

15 GEOTHERMAL ENERGY↗

Insert Modeling in UNF ST&DARDS

The Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS) is a software tool that integrates a used nuclear fuel (UNF) or spent nuclear fuel (SNF) relational database and key analysis capabilities to simplify and automate numerous UNF management and fuel cycle–related activities. UNF-ST&DARDS is being developed for the US Department of Energy’s Office of Nuclear Energy Spent Fuel and Waste Disposition program. UNF-ST&DARDS provides an integrated framework that uses advanced modeling and simulation to predict the behavior of SNF over the timescales associated with permanent disposal in a geologic repository. After leaving the spent fuel pool, SNF is transferred to dry storage in a dual-purpose canister (DPC). DPCs are considered “dual purpose” because they are designed for both storage and transportation, removing the need to transfer the fuel to a separate transportation cask. However, much research has been conducted investigating the feasibility of directly disposing of DPCs in geologic repositories. Direct disposal of DPCs could reduce worker exposure during repackaging, reducing the amount of low-level waste from the discarded DPCs and potentially saving billions of dollars. Therefore, direct disposal of as-loaded DPCs is desirable if it can be done safely.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Model for Pair Production Limit Cycles in Pulsar Magnetospheres

Abstract It was recently proposed that the electric field oscillation as a result of self-consistent e ± pair production may be the source of coherent radio emission from pulsars. Direct particle-in-cell simulations of this process have shown that the screening of the parallel electric field by this pair cascade manifests as a limit cycle, as the parallel electric field is recurrently induced when pairs produced in the cascade escape from the gap region. In this work, we develop a simplified time-dependent kinetic model of e ± pair cascades in pulsar magnetospheres that can reproduce the limit-cycle behavior of pair production and electric field screening. This model includes the effects of a magnetospheric current, the escape of e ± , as well as the dynamic dependence of pair production rate on the plasma density and energy. Using this simple theoretical model, we show that the power spectrum of electric field oscillations averaged over many limit cycles is compatible with the observed pulsar radio spectrum.

Astronomy & Astrophysics↗

MicroBooNE Public Data Sets: a Collaborative Tool for LArTPC Software Development

Among liquid argon time projection chamber (LArTPC) experiments MicroBooNE is the one that continually took physics data for the longest time (2015-2021), and represents the state of the art for reconstruction and analysis with this detector. Recently published analyses include oscillation physics results, searches for anomalies and other BSM signatures, and cross section measurements. LArTPC detectors are being used in current experiments such as ICARUS and SBND, and being planned for future experiments such as DUNE. MicroBooNE has recently released to the public two of its data sets, with the goal of enabling collaborative software developments with other LArTPC experiments and with AI or computing experts. These data sets simulate neutrino interactions on top of off-beam data, which include cosmic ray background and noise. The data sets are released in two formats: the native art/ROOT format used internally by the collaboration and familiar to other LArTPC experts, and the HDF5 format which contains reduced and simplified content and is suitable for usage by the broader community. This contribution presents the open data sets, discusses their motivation, the technical implementation, and the extensive documentation -- all inspired by FAIR principles. Finally, opportunities for collaborations are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Robust multi‐loop control of a floating wind turbine

Abstract A principal challenge facing the control of floating offshore wind turbines (FOWTs) is the problem of instability, or “negative damping,” when using blade pitch feedback to control generator speed. This closed‐loop instability can be attributed to non‐minimum phase zeros in the transfer function from blade pitch to generator speed. Standard approaches to improving stability and performance include robust tuning of control gains and introducing multiple feedback loops to respond to platform motion. Combining these approaches is nontrivial because multiple control loops complicate the impact of coupling in the system dynamics. The single‐loop approach to analyzing stability robustness neglects inter‐loop coupling, while a simplistic multi‐loop approach is highly sensitive to dimensional scaling and overestimates the robustness of the single‐loop controller. This work proposes a sensitivity representation that separates some of the natural FOWT dynamic coupling into a parallel feedback loop in the sensitivity function loop to address both of these concerns. The modified robustness measure is used with a simplified linear FOWT model to optimize scheduled multi‐loop control parameters in an automated tuning procedure. This controller is implemented for the 10‐MW Ultraflexible Smart FLoating Offshore Wind Turbine (USFLOWT) and compared against conventional single‐ and multi‐loop controllers tuned using frequency‐domain analysis and high‐fidelity OpenFAST simulations. The multi‐loop robust controller shows the highest overall performance in generator speed regulation and tower load reduction, though consideration of power quality, actuator usage, and other structural loading leads to additional trade‐offs.

Stockhouse, David↗

A Statistics-Based Threshold for the RMS-Energy Oscillation Detector

System operators and reliability coordinators currently rely on extensive baselining studies to set thresholds for their oscillation detectors. The resulting thresholds are based largely on engineering judgment and may vary significantly between organizations. In this paper, statistical distribution theory is used to derive a detection threshold for the widely deployed root-mean-square (RMS)-energy detector. This expression provides a theoretical basis for the detector's configuration, simplifies the process of selecting the threshold, and enables improved consistency among organizations that need to coordinate during system-wide events. Three methods for calculating the threshold using synchrophasor measurements are also proposed. These methods ensure that the threshold can be calculated reliably for various applications. Tests with simulated and field-measured data demonstrate that the statistics-based threshold provides consistent detection of grid disturbances while maintaining a low probability of false alarm.

Follum, James D.↗

SCALE depletion capabilities for molten salt reactors and other liquid-fueled systems

Nuclear reactor systems that use fuel dissolved in a liquid have the potential for enhanced safety characteristics, improved fuel-cycle outcomes, and more efficient isotope-production configurations. In these reactor systems, the fueled liquid may simultaneously undergo irradiation, physical and chemical removal processes, and fueling. The modeling and simulation of this transmutation and decay with material additions and removals is an ongoing research area. An accurate simulation tool is critical to the reactor and fuel-cycle design, reactor deployment, and source-term characterization for these advanced reactor systems. The work described herein involved implementing, testing, and applying the capability to perform reactor physics simulations within the Oak Ridge National Laboratory-developed SCALE suite for nuclear systems analyses and design, leveraging much of its pedigree in quality-assurance and reactor-analysis capabilities. The functionalities to simulate irradiation with material feeds and removals had been added in ORIGEN, and the TRITON reactor physics sequence was extended to calculate the total removed material and track external nonirradiated mixtures to estimate separate processing or waste streams. Results from these capabilities align with analytical expectations obtained from ORIGEN for simplified test cases and with expectations for a molten salt reactor application. This implementation, available with the SCALE 6.3 release, provides for a more efficient and accurate material accountability methodology, allowing for the characterization, design, and analysis of the complete isotopic material inventory of advanced liquid-fueled systems for a variety of applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Scattering in the Ising model with the quantum Lanczos algorithm *

Abstract Time evolution and scattering simulation in phenomenological models are of great interest for testing and validating the potential for near-term quantum computers to simulate quantum field theories. Here, we simulate one-particle propagation and two-particle scattering in the one-dimensional transverse Ising model for 3 and 4 spatial sites with periodic boundary conditions on a quantum computer. We use the quantum Lanczos algorithm to obtain all energy levels and corresponding eigenstates of the system. We simplify the quantum computation by taking advantage of the symmetries of the system. These results enable us to compute one- and two-particle transition amplitudes, particle numbers for spatial sites, and the transverse magnetization as functions of time. The quantum circuits were executed on various IBM Q superconducting hardware. The experimental results are in very good agreement with the values obtained using exact diagonalization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proxy Applications for Converged Workloads: DMC LDRD Initiative

Modern scientific applications are complicated and require coordination of several components. Proxy application driven software-hardware co-design plays a vital role in driving innovation among the developments of applications, software infrastructure and hardware architecture. Proxy applications are self-contained and simplified codes that are intended to model the performance-critical computations within applications. Applications executing on modern High Performance Computing (HPC) systems are susceptible to network congestion, insufficient memory bandwidth within and across compute nodes, and inadvertent loss of performance due to bugs and unoptimized programming models. Modern numerical simulations and machine learning models play a critical role in studying physical phenomenon under myriad uncertainties. Such applications often exhibit irregular computation and memory accesses at specific regions of the application code, which can contribute to various performance bottlenecks at scale. To mitigate such issues and prepare the next generation hardware for a variety of computation and data movement contingencies, a well-known practice is to consider "proxy" applications as representative motifs for various classes of scientific applications. While there is disagreement in the HPC community on the mechanisms of construction of the proxy applications, there is a strong consensus on their positive impact in co-design. Proxy Applications for Converged Workloads (PACER) is about facilitating software-hardware co-design through proxy applications with the goal of improving the performance of converged science workflows on heterogeneous systems.

97 MATHEMATICS AND COMPUTING↗

Cross sections for neutron-induced reactions from surrogate data: Reexamining the Weisskopf-Ewing approximation for ( n , n ' ) and ( n , 2 n ) reactions

Background: Modeling nuclear reaction networks for nuclear science applications and for simulations of astrophysical environments relies on cross section data for a vast number of reactions, many of which have never been measured. Cross sections for neutron-induced reactions on unstable nuclei are particularly scarce, since they are the most difficult to measure. Consequently, we must rely on theoretical predictions or indirect measurements to obtain the requisite reaction data. For compound nuclear reactions, the surrogate reaction method can be used to determine many cross sections of interest. Purpose: Earlier work has demonstrated that cross sections for neutron-induced fission and radiative neutron capture can be determined from a combination of surrogate reaction data and theory. For the fission case, it was shown that the Weisskopf-Ewing approximation, which significantly simplifies the implementation of the surrogate method, can be employed. Capture cross sections cannot be obtained, and require a detailed description of the surrogate reaction process. Here, we examine the validity of the Weisskopf-Ewing approximation for determining unknown (n, n') and (n, 2n) cross sections from surrogate data. Methods: Using statistical reaction calculations with realistic parametrizations, we investigate first whether the assumptions underlying the Weisskopf-Ewing approximation are valid for (n, n') and (n, 2n) reactions on representative target nuclei. We then produce simulated surrogate reaction data and assess the impact of applying the Weisskopf-Ewing approximation when extracting (n, n') and (n, 2n) cross sections in situations where the approximation is not strictly justified. Results: We find that peak cross sections can be estimated using the Weisskopf-Ewing approximation, but the shape of the (n, n') and (n, 2n) cross sections, especially for low neutron energies, cannot be reliably determined without accounting for the angular-momentum differences between the neutron-induced and surrogate reaction. Conclusions: To obtain reliable (n, n') and (n, 2n) cross sections from surrogate reaction data, a detailed description of the surrogate reaction mechanisms is required. To do so for the compound-nucleus energies and decay channels relevant to these reactions, it becomes necessary to extend current modeling capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Framework for idealized climate simulations with spatiotemporal stochastic clouds and planetary-scale circulations

In climate predictions, clouds are the leading source of uncertainty. This is partly because, to simulate the fluid dynamics of climate over the entire globe, a large grid spacing must be used, so clouds are a subgrid-scale parametrization rather than a resolved feature. Here, a framework is investigated with finer grid spacing of O(1) or O(10)km so that some clouds are not subgrid-scale; instead, clouds evolve on the numerical grid. This cloud evolution is achieved using stochastic modeling. Hence, the framework is idealized in the sense that the full fluid dynamics of cloud circulations is still not resolved, and simplified vertical structures are used. Nevertheless, the fluid dynamics model includes evolving clouds that interactively adjust in size, shape, lifetime, and regional coverage. In addition, different cloud types are included with different roles in the climate system, including deep convective clouds and also boundary-layer clouds such as shallow cumulus and stratocumulus clouds. Other basic aspects of the idealized climate system are planetary-scale circulations (e.g., Walker circulation) and radiation. With these ingredients (evolving clouds, planetary-scale circulations, and radiation), the framework has the potential for idealized investigations of climate change with interactive cloud–radiative feedback of individual clouds. Here, the formulation of the model equations is presented, and numerical simulations are shown to illustrate the model dynamics and climate change.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗

A Robust Neural Network for Extracting Dynamics from Electrostatic Force Microscopy Data

Advances in scanning probe microscopy (SPM) methods such as time-resolved electrostatic force microscopy (trEFM) now permit the mapping of fast local dynamic processes with high resolution in both space and time, but such methods can be time-consuming to analyze and calibrate. Here, we design and train a regression neural network (NN) that accelerates and simplifies the extraction of local dynamics from SPM data directly in a cantilever-independent manner, allowing the network to process data taken with different cantilevers. We validate the NN’s ability to recover local dynamics with a fidelity equal to or surpassing conventional, more time-consuming, calibrations using both simulated and real microscopy data. We apply this method to extract accurate photoinduced carrier dynamics on n = 1 butylammonium lead iodide, a halide perovskite semiconductor film that is of interest for applications in both solar photovoltaics and quantum light sources. Lastly, we use SHapley Additive exPlanations to evaluate the robustness of the trained model, confirm its cantilever-independence, and explore which parts of the trEFM signal are important to the network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗