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 451 records · Page 25

Evaluation of Structural Lattices for a Davis Gun Earth Penetrator Impact Experiment

The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.

36 MATERIALS SCIENCE↗

Modeling Codes Used for LAMP

LAMP team is planning on using multiple computer codes as tools for modeling beam dynamics, electromagnetics, and plasma environment in ion sources. Below is a description of the computer tools we have been using or planning to use in the future. Most of the codes chosen for the LAMP project are capable to model most of the beam dynamics in LAMP. A judicious choice of the appropriate code to use will be based on speed and fidelity of the modeling with faster calculations produced by simplifying some features.

43 PARTICLE ACCELERATORS↗

Web-Based Tools for Data-Informed Remedy Optimization: Software Theory and User Guide

This report documents the development and application of two web-based decision-support tools for pump-and-treat (P&T) groundwater remediation systems: PTOLEMY (Pump-and-Treat Optimized Location Evaluation to Maximize Yields) and OPTIMA (Optimization for Pump-and-Treat Implementation, Management, & Assessment). These tools enhance remedy design and management by leveraging advanced computational methods – specifically deep learning and multi-objective optimization – within a user-friendly platform. By integrating data-driven models with established hydrogeological knowledge, PTOLEMY and OPTIMA enable more efficient evaluation of well placement and operational strategies, helping site managers balance multiple remediation objectives under complex conditions. Both tools are implemented as modules within the SOCRATES (Suite Of Comprehensive Rapid Analysis Tools for Environmental Sites) web platform, which provides data access, visualization, and analytics to support remedy optimization across sites in the U.S. Department of Energy Office of Environmental Management complex. PTOLEMY is a rapid screening module designed to identify promising locations for new extraction wells. It employs a multi-channel three-dimensional convolutional neural network (MC3D-CNN) trained on high-fidelity simulation data to predict the relative performance (in terms of contaminant mass recovery) of potential well sites. Through an interactive web interface, PTOLEMY visualizes the probability of high performance across a site, highlighting areas where an extraction well is likely to yield above-threshold contaminant removal over a multi-year period. PTOLEMY’s map-based displays and exportable results support transparent communication of screening analyses. By focusing attention on the most favorable candidate locations, the tool augments traditional engineering judgment and physics-based modeling, providing a data informed basis for subsequent detailed evaluations. OPTIMA is a multi objective optimization module designed to find wellfield layouts and operating schedules that meet various cleanup goals. It quickly evaluates thousands of candidate setups – combinations of well locations, timing, and rates – and returns a small set of best trade-off options for comparison. At its core, OPTIMA uses a U-Net-based surrogate model – a deep-learning emulator of a groundwater flow and transport simulator – to dramatically accelerate scenario evaluations. Coupling this fast surrogate with the NSGA-II (Non-dominated Sorting Genetic Algorithm II) evolutionary algorithm, OPTIMA explores a wide decision space of well locations and schedules to identify Pareto-optimal solutions that trade off key objectives (e.g., minimizing cleanup time, maximizing contaminant mass removal, and minimizing plume extent). The tool outputs a family of optimal configurations and visualizes their trade-offs (Pareto frontiers of cleanup metrics and maps of optimized well placements). Site managers can use these results to understand the range of viable strategies and to select candidate designs for more detailed verification. OPTIMA is currently under active development and not yet fully released; this guide provides early documentation to support planning and gather user feedback.

54 ENVIRONMENTAL SCIENCES↗

Addressing GPU memory limitations for Graph Neural Networks in High-Energy Physics applications

Introduction Reconstructing low-level particle tracks in neutrino physics can address some of the most fundamental questions about the universe. However, processing petabytes of raw data using deep learning techniques poses a challenging problem in the field of High Energy Physics (HEP). In the Exa.TrkX Project, an illustrative HEP application, preprocessed simulation data is fed into a state-of-art Graph Neural Network (GNN) model, accelerated by GPUs. However, limited GPU memory often leads to Out-of-Memory (OOM) exceptions during training, due to the large size of models and datasets. This problem is exacerbated when deploying models on High-Performance Computing (HPC) systems designed for large-scale applications. Methods We observe a high workload imbalance issue during GNN model training caused by the irregular sizes of input graph samples in HEP datasets, contributing to OOM exceptions. We aim to scale GNNs on HPC systems, by prioritizing workload balance in graph inputs while maintaining model accuracy. Our paper introduces diverse balancing strategies aimed at decreasing the maximum GPU memory footprint and avoiding the OOM exception, across various datasets. Results Our experiments showcase memory reduction of up to 32.14% compared to the baseline. We also demonstrate the proposed strategies can avoid OOM in application. Additionally, we create a distributed multi-GPU implementation using these samplers to demonstrate the scalability of these techniques on the HEP dataset. Discussion By assessing the performance of these strategies as data loading samplers across multiple datasets, we can gauge their effectiveness in both single-GPU and distributed environments. Our experiments, conducted on datasets of varying sizes and across multiple GPUs, broaden the applicability of our work to various GNN applications that handle input datasets with irregular graph sizes.

Lee, Claire Songhyun↗

Low-speed, high-lift aerodynamic characteristics of slender, hypersonic accelerator-type configurations

Two investigations were conducted in the Langley 14 by 22 Foot Subsonic Tunnel to determine the low-speed aerodynamic characteristics of a generic hypersonic accelerator-type configuration. The model was a delta wing configuration incorporating a conical forebody, a simulated wrap-around engine package, and a truncated conical aftbody. Six-component force and moment data were obtained over a range of attack from -4 to 30 degrees and for a sideslip range of + or - 20 degrees. In addition to tests of the basic configuration, component build-up tests were conducted; and the effects of power, forebody nose geometry, canard surfaces, fuselage strakes, and engines on the lower surface alone were also determined. Control power available from deflections of wing flaps and aftbody flaps was also investigated and found to be significantly increased during power-on conditions. Large yawing moments resulted from asymmetric flow fields exhibited by the forebody as revealed by both surface pressure data and flow visualization. Increasing nose bluntness reduced the yawing-moment asymmetry, and the addition of a canard eliminated the yawing-moment asymmetry.

Gatlin, Gregory M.↗

Ion and relativistic electron transport in solar flares

Recent models for the ion and relativistic electron transport in solar flare magnetic loops are reviewed, focusing on the ways in which the models handle loop structure, particle acceleration, and particle transport. Results are presented from Monte Carlo simulations of ion and relativistic ion transport. Consideration is given to results concerning the depth distributions of gamma-ray production, the attenuation of gamma-ray lines, the time dependences of the various emissions and the angular distribution of the bremsstrahlung.

Ramaty, R.↗

Numerical Simulations of the Slingatron

The slingatron mass accelerator is described for several track configurations (shapes), and numerical simulations of this accelerating mass traversing a given track configuration are presented. The sled is modeled as a point mass that interacts with the slingatron track using both a conventional and a new empirical velocity dependent friction law. The closed loop circular slingatron was found to produce high maximum sled velocities provided the gyration angular speed is always increasing. In contrast several spiral shaped slingatron tracks reveal that high maximum sled velocities are obtainable with the gyration speed held constant. In fact, a slingatron constructed out of semi-circles is shown capable of generating high velocity sleds in such a way that no initial sled injection is necessary. Choosing the proper initial gyration phase with an empirically determined friction model allows the mass sled to gain ever-increasing velocities when placed in a semi-circle slingatron. The sled bearing pressure and its total acceleration are examined and presented.

Cooper, Gene R.↗

Hubble Space Telescope Degradation Data Used for Ground-Based Durability Projection of Insulation on the International Space Station

Ground-based environmental durability tests have indicated that exposing materials in accelerated tests to known spacecraft mission degradation sources predicted by an environmental model does not simulate the extent of damage that occurs in the space environment. The reasons for this may include the complex nature of the space environment, which is not simulated completely in any ground-based facility, the extreme differences in exposure rates in space and in ground tests, and inaccuracies in environmental models. One approach to overcoming the difficulties in simulating the space environment using ground-based testing is to calibrate the facility using data from actual space-exposed materials to determine exposure levels required to replicate degraded properties observed in space. Research was conducted at the NASA Glenn Research Center to develop a ground-to-space correlation method that determines the durability of Teflon-based insulation for the International Space Station (ISS) by using data obtained in a ground facility and degraded Teflon thermal insulation retrieved from the Hubble Space Telescope.

deGroh, Kim K.↗

Power Electronics Materials and Bonded Interfaces - Reliability and Lifetime

The emergence of wide-bandgap (WBG) devices in power electronics packages has brought significant attention to the reliability and thermal performance of the bonded interfaces. These interfaces, such as die-attach and substrate-attach, need to withstand the severe stress induced on them by the high-temperature (> 200°C) operation of the WBG devices. Sintered silver and transient liquid phase alloys are key bonded interfaces that can enable the high-temperature operation of WBG devices, however it is essential to evaluate their reliability and study the possible failure mechanisms before they can be used in power electronics packages. Sintered silver is already being adopted in commercial packages as a die-attach, but large-area substrate-attach offers additional reliability challenges. In this project, we investigate the reliability of these interface materials under accelerated thermal cycling and correlate these results with finite element simulations to develop lifetime prediction models. The failure mechanisms of these bonded interfaces are also studied.

47 OTHER INSTRUMENTATION↗

Engagement: Hyperparameter Optimization of Generative Adversarial Network Models for High-Energy Physics Simulations

We present our SciDAC FASTMath-HEP partnership results for tuning generative adversarial models (GANs) for high energy physics applications. The GANs are used in hybrid simulations to accelerate otherwise time-consuming computations. We optimize for both, prediction accuracy and variability with the goal to find GAN architectures that are reliable and robust.

high energy physics↗

Sensitivity to Snapshot Frequency in the POD-Based Reduced-Order Modeling of Flow over a Gaussian Bump

The flow over a Gaussian bump at a Reynolds number of 10^6, based on the bump length and freestream conditions, and a Mach number of 0.2, was simulated using wall-modeled large eddy simulation (WMLES). The incoming turbulent flow exhibits relaminarization upstream of the bump apex, due to the favourable pressure gradient acceleration induced by the bump curvature, as well as separation, reattachment and a detached shear layer downstream of the apex. A Galerkin reduced-order model (ROM) was then constructed using a proper orthogonal decomposition (POD) modal basis. The complex, three-dimensional flow features are captured by spatial POD modes differently based on the frequency of the high-fidelity snapshots originally provided. POD modes constructed from differently sampled snapshot matrices were then truncated to retain 90% of the turbulence kinetic energy in the domain of interest. A parametric study is conducted to evaluate the effect of the frequency of snapshots used to construct the POD-ROM.

turbulence↗

Multi-GPU porting of a phase-change cascaded lattice Boltzmann method for three-dimensional pool boiling simulations

The Lattice Boltzmann method (LBM) has proven effective in simulating phase-change phenomena, such as melting, solidification, evaporation, and boiling. In this work, we develop a highly parallelized multi-GPU implementation of LBM for three-dimensional pool boiling simulations. The code is based on the OpenACC programming model, which enables the code to be deployed efficiently on multi-core CPUs, GPUs, and potentially other accelerators, without the need for architecture-specific rewrites. To support large-scale simulations, the domain is decomposed and distributed across multiple compute nodes using MPI. We demonstrate that the code exhibits excellent scaling properties, with ideal strong-scaling running with up to 256 GPUs on the MareNostrum5 cluster.

97 MATHEMATICS AND COMPUTING↗

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗

Structural dynamic model obtained from flight use with piloted simulation and handling qualities analysis

The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted simulation. and handling qualities analysis has been developed. This technique was demonstrated using SR-71 flight test data. For the SR-71 aircraft, the most significant structural response is the longitudinal first-bending mode. This mode was modeled as a second-order system, and the other higher order modes were modeled as a time delay. The distribution of the modal response at various fuselage locations was developed using a uniform beam solution, which can be calibrated using flight data. This approach was compared to the mode shape obtained from the ground vibration test, and the general form of the uniform beam solution was found to be a good representation of the mode shape in the areas of interest. To calibrate the solution, pitch-rate and normal-acceleration instrumentation is required for at least two locations. With the resulting structural model incorporated into the simulation, a good representation of the flight characteristics was provided for handling qualities analysis and piloted simulation.

Powers, Bruce G.↗

An accelerated approach to inline non-LTE modeling

A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗