Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Self-accelerated corrosion of nuclear waste forms at material interfaces

The US plan for high-level nuclear waste includes the immobilization of long-lived radionuclides in glass or ceramic waste forms in stainless-steel canisters for disposal in deep geological repositories. Here we report that, under simulated repository conditions, corrosion could be significantly accelerated at the interfaces of different barrier materials, which has not been considered in the current safety and performance assessment models. Severe localized corrosion was found at the interfaces between stainless steel and a model nuclear waste glass and between stainless steel and a ceramic waste form. The accelerated corrosion can be attributed to changes of solution chemistry and local acidity/alkalinity within a confined space, which significantly alter the corrosion of both the waste-form materials and the metallic canisters. The corrosion that is accelerated by the interface interaction between dissimilar materials could profoundly impact the service life of the nuclear waste packages, which, therefore, should be carefully considered when evaluating the performance of waste forms and their packages. Moreover, compatible barriers should be selected to further optimize the performance of the geological repository system.

Guo, Xiaolei↗
DOE OSTI2023

Accelerating Battery Simulations by Using High Performance Computing and Opportunities with Machine Learning

The US Department of Energy (DOE) estimates that the battery pack cost of $60 per kilowatt-hour while increasing the driving range to over 300 miles and vehicle charging under 15 min or less would enable mass penetration of electric vehicles in the USA by 2030. These projections are based on the currently available high-density cell chemistry combined with a system level design and optimized electrodes. Electrodes for current state-of-the-art lithium-ion cell technology are fabricated from electrode slurry comprising of active material, polymeric binder, and conductive diluent such as carbon black that are coated on metal current collectors such as copper and aluminum. Given these challenging requirements for development of electrical energy storage devices for future transportation needs, a predictive simulation capability which can accelerate design by considering performance and safety implications of different geometry, materials, and chemistry choices is required. In this chapter, we discuss our state-of-the-art three-dimensional modeling framework, providing examples on battery performance and safety simulations. We also present approaches to use machine learning in the context of large datasets that are becoming available.

Allu, Srikanth↗
NASA NTRS2018

In-Flight Verification and Validation of Colloid Microthruster Performance

Colloid Micronewton Thrusters (CMNTs) use an electrospray to provide precision spacecraft position and pointing control. They were demonstrated in space for the first time as part of NASA’s Space Technology 7 (ST7) payload hosted by the European Space Agency’s (ESA’s) LISA Pathfinder (LPF) technology demonstration mission in January, 2016. CMNTs were the actuator in the disturbance reduction system (DRS) that provided drag-free operation of the LPF spacecraft, which will be necessary for future gravity wave observatories such as ESA’s Laser Interferometer Space Antenna (LISA) mission, currently in Phase A and scheduled for launch in 2034. The CMNT technology met performance requirements operating at 5-30 µN of thrust with 0.1 µN resolution and ≤0.1 µN/ÖHz thrust noise to deliver the required nanometer-level precision spacecraft control measured by the gravitational reference sensor (GRS) in the ESA LISA Technology Package (LTP). The performance of seven of the eight CMNTs in flight was consistent with ground test results, and as a system, all eight thrusters met mission-level performance requirements. The colloid microthruster performance model of thrust and thrust noise as a function of operational parameters (i.e. beam current, voltage, temperature, etc.) was validated in flight over a wide range of conditions. A model and simulation of the thruster control algorithm was developed and validated with flight data to predict thrust noise. This capability is important for future missions because it relates directly to the acceleration noise on the test masses, which provide the gravity wave measurements. The CMNT thruster model data and validation with LISA Pathfinder/ST7-DRS flight experiments are presented in this paper.

"O'Donnell, James"↗
DOE OSTI2024

Experimental Safe Extremum Seeking for Accelerators

Here, we demonstrate the recent designs of safe extremum seeking (Safe ES) on the 1-km-long charged particle accelerator at the Los Alamos Neutron Science Center (LANSCE). Safe ES is a modification of extremum seeking (ES) which, in addition to minimizing an analytically unknown cost, also employs a safety filter based on an analytically unknown control barrier function (CBF) safety metric. Tuning is necessitated by accelerators being large complex systems, with many drifting parameters due to thermal effects and degradation. At the same time, safe operation (the maintenance of state constraints) is crucial, as damage brings astronomical costs, both financially and in operation downtime. Our measured (but analytically unknown) safety metric is the beam current. We perform multivariable Safe ES on three accelerator applications, in which we adapt 4, 6, and 3 magnet strength parameters, respectively. Two of the three applications are for validated simulation models of beamlines at LANSCE: the first for the proton radiography (pRad) beamline of 800-MeV protons for spot size tuning; the second on a high-performance code, HPSim, for tuning the low-energy beam transport (LEBT) region that contains a beam of 750-keV protons. The third is an experimental tuning of the steering magnets in the LEBT at LANSCE.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Parallel quantum computing simulations via quantum accelerator platform virtualization

Quantum circuit execution is a central task in quantum computation. Due to inherent quantum-mechanical constraints, quantum computing workflows often involve a considerable number of independent measurements over a large set of slightly different quantum circuits. Here we discuss a simple model for parallelizing such quantum circuit executions that is based on introducing a large array of virtual quantum processing units (mapped to HPC nodes in our case) as a parallel quantum computing platform. Implemented within the XACC framework, the model can readily take advantage of its backend-agnostic features, enabling parallel quantum computing/simulation over any target backend supported by XACC. We illustrate the performance of this approach by demonstrating strong scaling in two pertinent domain science problems, namely in computing the gradients for the multi-contracted variational quantum eigensolver and in data-driven quantum circuit learning, where we vary the number of qubits and the number of circuit layers. Here, the latter simulation leverages the cuQuantum library to run efficiently on GPU-accelerated HPC platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗
DOE OSTI2023

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗
NASA NTRS2008

State of the Art Assessment of Simulation in Advanced Materials Development

Advances in both the underlying theory and in the practical implementation of molecular modeling techniques have increased their value in the advanced materials development process. The objective is to accelerate the maturation of emerging materials by tightly integrating modeling with the other critical processes: synthesis, processing, and characterization. The aims of this report are to summarize the state of the art of existing modeling tools and to highlight a number of areas in which additional development is required. In an effort to maintain focus and limit length, this survey is restricted to classical simulation techniques including molecular dynamics and Monte Carlo simulations.

Wise, Kristopher E.↗
DOE OSTI2020

Understanding the eco-geomorphologic feedback of coastal marsh under sea level rise: vegetation dynamic representations, processes interaction, and parametric sensitivity

A growing number of coastal eco-geomorphologic modeling studies have been conducted to understand coastal marsh evolution under sea-level rise (SLR). Although these models quantify marsh topographic change as a function of sedimentation and erosion, their representations of vegetation dynamics that control organic sedimentation differ. How vegetation dynamic schemes contribute to simulation outcomes is not well quantified. Additionally, the sensitivity of modeling outcomes to parameter selection in the available formulations has not been rigorously tested to date, especially under the influence of an accelerating SLR. In this paper, we used a coastal eco-geomorphologic model with different vegetation dynamic schemes to investigate the eco-geomorphologic feedbacks of coastal marshes and parametric sensitivity under SLR scenarios. We found that marsh platform relief increased with sea level rise rate. The simulations with different vegetation schemes exhibited different spatial-temporal variations in elevation and biomass. The nonlinear Spartina scheme presented the most resilient prediction with generally the highest marsh accretion and vegetation biomass, and the least elevation relief under SLR. But the linear Spartina scheme predicts the lowest unvegetated-vegetated ratio. We also found that vegetation-related parameters and sediment diffusivity, which were not well measured or discussed in previous studies, were identified as some of the most critical parameters. Additionally, the model sensitivity to vegetation-related parameters increased with SLR rates. The identified most sensitive parameters may inform how to appropriately choose modeling representations of key processes and parameters for different coastal marsh landscapes under SLR, and demonstrate the importance of future field measurements of these key parameters.

58 GEOSCIENCES↗
NASA NTRS2022

Megafires in a Warming World: What Wildfire Risk Factors Led to California’s Largest Recorded Wildfire

Massive wildfires and extreme fire behavior are becoming more frequent across the westernUnited States, creating a need to better understand how megafire behavior will evolve in our warmingworld. Here, the fire spread model Prometheus is used to simulate the initial explosive growth ofthe 2020 August Complex, which occurred in northern California (CA) mixed conifer forests. Hightemperatures, low relative humidity, and daytime southerly winds were all highly correlated withextreme rates of modeled spread. Fine fuels reached very dry levels, which accelerated simulationgrowth and heightened fire heat release (HR). Model sensitivity tests indicate that fire growth andHR are most sensitive to aridity and fuel moisture content. Despite the impressive early observedgrowth of the fire, shifting the simulation ignition to a very dry September 2020 heatwave predicted a>50% increase in growth and HR, as well as increased nighttime fire activity. Detailed model analysesof how extreme fire behavior develops can help fire personnel prepare for problematic ignitions.

Kevin Varga↗
NASA NTRS1979

Experiments in sensing transient rotational acceleration cues on a flight simulator

Results are presented for two transient motion sensing experiments which were motivated by the identification of an anomalous roll cue (a 'jerk' attributed to an acceleration spike) in a prior investigation of realistic fighter motion simulation. The experimental results suggest the consideration of several issues for motion washout and challenge current sensory system modeling efforts. Although no sensory modeling effort is made it is argued that such models must incorporate the ability to handle transient inputs of short duration (some of which are less than the accepted latency times for sensing), and must represent separate channels for rotational acceleration and velocity sensing.

Parrish, R. V.↗
DOE OSTI2021

3D Theoretical and simulation tools for microbunched cooling

The concept of microbunched electron cooling (MBEC, introduced in Ref. [1]) is currently under investigation as a candidate for the strong hadron cooling (SHC) component of EIC. Most of the theoretical treatments of this scheme are based on a simplified model for the longitudinal space charge [2–4]. In this model, the particle beams are treated as collections of thin charged disks, with the interaction function for two such disks being calculated by averaging the forces between constituent point charges. Although the disk model is flexible and captures a large part of the underlying physics, it lacks the fidelity and robustness of an analysis that explicitly treats the beams as collections of point charges. In what follows, we explore such a model, both from a theory-based and a simulation perspective.

43 PARTICLE ACCELERATORS↗
NASA NTRS2017

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗
DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗
DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 Mark-I Station Blackout: Part II – Uncertainty Analysis

Single case comparisons between severe accident simulations can provide detailed insights into severe accident model behavior, however, they cannot offer insights into model uncertainty, sensitivity to uncertain parameters, or underlying model biases.Here in this analysis, the single case benchmark comparison of the MELCOR material interaction models for a station blackout (SBO) scenario of a boiling water reactor (BWR) using representative Fukushima Daiichi Unit 1 boundary conditions is expanded to include an uncertainty analysis. As part of this uncertainty analysis, 1200 simulations are performed for each material interaction model (2400 total), with random sampling of 14 uncertain MELCOR input parameters. Input parameters are selected for their impact on models representing core degradation processes. These include candling, fuel rod failure, debris quenching and dryout. The analysis performed here is not a traditional “best-estimate” uncertainty analysis that uses best-estimate parameters or identifies best-estimate figure of merit distributions. Instead, it is an exploratory uncertainty analysis that identifies and interrogates underlying model form biases of the two material interaction models (eutectics and interactive materials models). Uniform distributions are applied to all uncertain parameters to ensure coverage of the model parameter uncertainty space. Key findings from this study include underlying model form biases exhibited by material interaction models, and notable differences in accident progression outcomes between the material interaction models. This uncertainty study extends and confirms the conclusions from the first part of this study, which compared the impact of material interaction modeling on simulation of a short-term station blackout scenario with representative Fukushima Daiichi Unit I boundary conditions. In particular, this study confirms that the eutectics model generally exhibits accelerated degradation and failure of fuel components, the core plate, and the lower head. The eutectics model also has a tendency to exhibit a greater degree of core degradation, greater debris mass formation, and larger debris mass ejection. Finally, the eutectics model exhibits higher maximum temperatures for fuel, cladding, particulate debris, oxidic molten pool, and metallic molten pool components than the interactive materials model; interactive materials model simulations exhibit a soft “limitation” on maximum temperatures that is related to the temperature at which material relocation occurs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Adsorption Hysteresis Under Control: Tuning Host–Guest Interactions via a Genetic Algorithm

Mesoporous adsorbent materials offer a large volumetric capacity; however, cyclic adsorption/desorption processes in these systems often suffer from hysteresis and may require a significant pressure swing to access this capacity. To mitigate hysteresis, a proposed strategy is to include nucleation sites on the walls of the mesoporous material to facilitate droplet and bubble formation, lowering the free energy barriers to the respective phase transitions. It is unclear, however, what combination of adsorbate− adsorbent interactions and spatial patterning would be beneficial for a given application, considering that improvements to some sorption properties may come at the expense of other attributes. To understand these interconnected observables, we examine two model systems, planar-slit and cylindrical pores with tunable interaction sites, using GPU-accelerated transition matrix Monte Carlo simulations. The simulations provide a free energy map of the pressure−adsorption space in a matter of minutes, which we use to track adsorption isotherm characteristics as a function of adsorbent properties. We then leverage the rapid acquisition of simulation data to construct a genetic algorithm to iteratively modify interaction sites of the slit-pore wall to minimize the hysteresis of this system without sacrificing uptake. We find that the adsorption branch of the isotherm is easily modulated via the average host−guest interaction strength, but desorption is only adjustable if there is a suitable bubble nucleation site. Within the context of a slit-pore system, we identify relative interaction strengths and patch sizes required to gain control over both branches of the hysteresis loop.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗