Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator Simulation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation

Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumvents assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

THIN CURRENT SHEETS AND ASSOCIATED ELECTRON HEATING IN TURBULENT SPACE PLASMA

Intermittent structures, such as thin current sheets, are abundant in turbulent plasmas. Numerical simulations indicate that such current sheets are important sites of energy dissipation and particle heating occurring at kinetic scales. However, direct evidence of dissipation and associated heating within current sheets is scarce. Here, we show a new statistical study of local electron heating within proton-scale current sheets by using high-resolution spacecraft data. Current sheets are detected using the Partial Variance of Increments (PVI) method which identifies regions of strong intermittency. We find that strong electron heating occurs in high PVI (>3) current sheets while no significant heating occurs in low PVI cases (<3), indicating that the former are dominant for energy dissipation. Current sheets corresponding to very high PVI (>5) show the strongest heating and most of the time are consistent with ongoing magnetic reconnection. This suggests that reconnection is important for electron heating and dissipation at kinetic scales in turbulent plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TH-CD-BRA-09: Towards Absolute Dose Measurement in MRI-Linac and Gamma-Knife: Design and Construction of An MR-Compatible Water Calorimeter

Purpose: The purpose of this work is to design and optimize a portable Water Calorimeter (WC) for use in a commercial MRI-linac and Gamma-knife in addition to conventional radiotherapy linacs. Water calorimeters determine absorbed dose to water at a point by measuring radiation-induced temperature rise of the volume (the two are related by the medium specific heat capacity). In this formalism, one important correction factor is heat transfer correction k-ht. It compensates for heat gain/loss due to conductive and convective effects, and is numerically calculated as ratio of temperature rise in the absence of heat loss to that in the presence of heat loss. Operating at 4°C ensures convection is minimal. Methods: A commercial finite element software was used to evaluate several WC designs with different insulation materials and thicknesses; channels allowing coolant to travel around WC (to sustain WC at 4°C) were modeled, and worst-case scenario variation in the temperature of the coolant was simulated for optimization purposes (2.6 mK/s). Additionally, several calorimeter vessel design parameters (front/back glass thickness/separation, diameter) were also simulated and optimized. Optimization is based on minimizing long term calorimeter drift (24h) as well as variation and magnitude of k-ht. Results: The final selected WC design reached a modest drift of 11µK/s after 15h for the worst-case coolant temperature variation. This design consists of coolant channels being encompassed on both sides by cryogel insulation. For the MRI-linac beam, glass thickness plays the largest effect on k-ht with variation of upto 0.6% in the first run for thicknesses ranging between 0.5–1.7mm. Subsequent runs vary only within 0.1% with glass thickness. Other factors such as vessel radius and top/bottom glass separation have sub 0.1% effects on k-ht. Conclusion: An MR-safe 4°C stagnant WC appropriate for dosimetry in MRI-linac and Gamma-Knife was designed, optimized, and construction is nearly completed. NSERC Discovery Grant RGPIN-435608.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

PARTICLE ACCELERATION AND THE ORIGIN OF X-RAY FLARES IN GRMHD SIMULATIONS OF SGR A*

Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform general relativistic magnetohydrodynamic simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

G4CMP: Condensed matter physics simulation using the G$\scriptsize{\mathrm{EANT}}$4 toolkit

G4CMP simulates phonon and charge transport in cryogenic semiconductor crystals using the Geant4 toolkit. The transport code is capable of simulating the propagation of acoustic phonons as well as electron and hole charge carriers. Herein processes for anisotropic phonon propagation, oblique charge-carrier propagation, and phonon emission by accelerated charge carriers are included. The simulation reproduces theoretical predictions and experimental observations such as phonon caustics, heat-pulse propagation times, and mean charge-carrier drift velocities. In addition to presenting the physics and features supported by G4CMP, this report outlines example applications from the dark matter and quantum information science communities. These communities are applying G4CMP to model and design devices for which the energy transported by phonons and charge carriers is germane to the performance of superconducting instruments and circuits placed on silicon and germanium substrates. The G4CMP package is available to download from GitHub: github.com/kelseymh/G4CMP.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Optimizing the emission, propagation, and focusing of an intense electron beam

Intense electron beams can be used to study the dynamical response of materials under shocks in order to adjust the models developed for hydrodynamics simulations. We present in this paper a characterization of beams produced in a field emission diode coupled to the generator RKA at CEA/CESTA. Cherenkov emission, produced by the beam interacting in a fused silica disk, was observed by fast optical cameras to estimate beam homogeneity. GEANT4 simulations were performed to estimate the transfer function of the silica target and to optimize the anode foil. First, we chose the best cathode material available among the most common materials used in field emission systems. In addition, we found that by optimization of the anode thickness, we could improve the spatial homogeneity of the beam which is of prime importance for computing the interaction of the beam with materials. Next, we changed the beam fluence by increasing the beam current and by reducing the beam radius. Finally, we studied the propagation and focusing of the electron beam in low pressure gases and observed that we could use self-magnetic field focusing in order to increase beam fluence at the target location. The experimental results are in good agreement with PIC simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Accelerating the Design of Solar Thermal Fuel Materials through High Throughput Simulations

Solar thermal fuels (STF) store the energy of sunlight, which can then be released later in the form of heat, offering an emission-free and renewable solution for both solar energy conversion and storage. However, this approach is currently limited by the lack of low-cost materials with high energy density and high stability. In this Letter, we present an ab initio high-throughput computational approach to accelerate the design process and allow for searches over a broad class of materials. The high-throughput screening platform we have developed can run through large numbers of molecules composed of earth-abundant elements and identifies possible metastable structures of a given material. Corresponding isomerization enthalpies associated with the metastable structures are then computed. Using this high-throughput simulation approach, we have discovered molecular structures with high isomerization enthalpies that have the potential to be new candidates for high-energy density STF. We have also discovered physical principles to guide further STF materials design through structural analysis. More broadly, our results illustrate the potential of using high-throughput ab initio simulations to design materials that undergo targeted structural transitions.

Liu, Y↗
DOE OSTI2022

Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays

Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Thin‐Film Ferroelectrics

Abstract Over the last 30 years, the study of ferroelectric oxides has been revolutionized by the implementation of epitaxial‐thin‐film‐based studies, which have driven many advances in the understanding of ferroelectric physics and the realization of novel polar structures and functionalities. New questions have motivated the development of advanced synthesis, characterization, and simulations of epitaxial thin films and, in turn, have provided new insights and applications across the micro‐, meso‐, and macroscopic length scales. This review traces the evolution of ferroelectric thin‐film research through the early days developing understanding of the roles of size and strain on ferroelectrics to the present day, where such understanding is used to create complex hierarchical domain structures, novel polar topologies, and controlled chemical and defect profiles. The extension of epitaxial techniques, coupled with advances in high‐throughput simulations, now stands to accelerate the discovery and study of new ferroelectric materials. Coming hand‐in‐hand with these new materials is new understanding and control of ferroelectric functionalities. Today, researchers are actively working to apply these lessons in a number of applications, including novel memory and logic architectures, as well as a host of energy conversion devices.

Fernandez, Abel↗
DOE OSTI2022

Probing the Kinetic Origin of Varying Oxidative Stability of Ethyl- vs. Propyl-spaced Amines for Direct Air Capture

Amine-based adsorbents are promising for direct air capture of CO 2 , yet oxidative degradation remains a key unmitigated risk hindering wide-scale deployment. Borrowing wisdom from the basic auto-oxidation scheme, insights are gained into the underlying degradation mechanisms of polyamines by quantum chemical, advanced sampling simulations, adsorbent synthesis, and accelerated degradation experiments. The reaction kinetics of polyamines are contrasted with that of typical aliphatic polymers and they elucidate for the first time the critical role of aminoalkyl hydroperoxide decomposition in the oxidative degradation of amino-oligomers. The experimentally observed variation in oxidative stability of polyamines with different backbone structures is explained by the relationship between the local chemical structure and the free energy barrier of aminoalkyl hydroperoxide decomposition, suggesting that its energetics can be used as a descriptor to screen and design new polyamines with improved stability. In conclusion, the developed computational capability sheds light on radical-induced degradation chemistry of other organic functional materials.

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

Mechanism of Nonequilibrium Formation of Novel Interface Structures in Explosive Welding

Novel structures are commonly observed in explosive welding products. This paper presents the nonequilibrium mechanism of the welding interface. A theoretical model was established first to deduce the motion equation of the rapid accelerated flyer plate. Numerical simulations of the welding interface were carried out with the smooth particle hydrodynamic method. Both the theoretical and numerical results were compared with experimental phenomena and showed that under the very rapid application of an impact load usually associated with high-strength plastic deformation, a significant strain rate, and temperature rise, the material volumes of the welding interface are in the viscoelastic state and prone to rotation, detrusion, and dilatation. In addition, the extreme values from the nonequilibrium development of the interface structure are a consequence of the competition between the loading and unloading stress waves during the explosive welding process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

FlowPM: Distributed TensorFlow implementation of the FastPM cosmological N-body solver

Here, we present FlowPM, a Particle-Mesh (PM) cosmological N-body code implemented in Mesh-TensorFlow for GPU-accelerated, distributed, and differentiable simulations. We implement and validate the accuracy of a novel multi-grid scheme based on multiresolution pyramids to compute large-scale forces efficiently on distributed platforms. We explore the scaling of the simulation on large-scale supercomputers and compare it with corresponding Python based PM code, finding on an average 10x speed-up in terms of wallclock time. We also demonstrate how this novel tool can be used for efficiently solving large scale cosmological inference problems, in particular reconstruction of cosmological fields in a forward model Bayesian framework with hybrid PM and neural network forward model. We provide skeleton code for these examples and the entire code is publicly available at https://github.com/modichirag/flowpm[Formula presented].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine learning accelerated discrete element modeling of granular flows

Granular flows are widely encountered in many industrial processes and natural phenomena. Discrete Element Modeling (DEM) is a useful tool for understanding and troubleshooting devices, which handle granular materials. However, its applicability is significantly limited by the huge computational cost associated with detecting and computing collisions. In this research, the computation speed of DEM was accelerated by orders of magnitude using a convolutional neural network to replace the direct calculation of particle-particle and particle-boundary collisions. The MFiX software was used to generate the training and testing dataset. Additionally, a GPU accelerated TensorFlow model was used to train the neural network and test the results. The model fluctuations caused by different training steps were reduced with a multi-scale loss function. The accuracy was improved with more frames within one training step. The modeling of a rotating drum and a hopper demonstrated the accuracy and efficiency of this machine learning accelerated DEM in the simulation of granular flows.

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

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2018

Re-additization of commercial biodiesel blends during long-term storage

Commercial biodiesel blends were aged at 43 degrees C while monitoring stability. The oxidation stability -- or oxidation reserve expressed as Rancimat induction period (IP) -- gradually decreased from its initial value. At a predetermined IP threshold, an antioxidant was used to restore IP to the ASTM D7467 specification minimum of 6 h, referred to as re-additization. At lower IP values, the amount of antioxidant required increased significantly, and the effectiveness tended to be reduced. Once IP fell to essentially zero, the acid content increased to above the allowable limit and insoluble material was also detected. Storage life was increased relative to the as-received fuels as evidenced by longer time to produce acids. Experience in the field may vary based on storage conditions; however, these results indicate re-additization can significantly increase storage life of biodiesel blends when used with regular monitoring of IP and acid number. An assessment of the storage stability of the as-received fuels showed that the initial IP did not predict storage behavior, although fuels above the specification minimum remained stable for >12 weeks accelerated aging (1 year simulated).

25 ENERGY STORAGE↗
DOE OSTI2026

Interpolation of computed gamma-ray detector response functions

Gamma-ray spectra measured by traditional detectors contain features that result from a combination of the effects of detector materials/geometry, the incident gamma-ray energy, and the angle of entry. The features, such as the full-energy photopeak, Compton continuum, annihilation peak, and escape peaks, are governed by simple relationships depending on incident energy and have been known for a long time. Monte Carlo computer simulations of gamma rays interacting with a detector will show these features, and with a resolution function applied, the results should look similar to real measurements. The traditional approach to creating a detector response function requires many separate simulations of monoenergetic gamma rays striking the detector. This paper presents a new approach to developing computed detector response functions. The new approach involves a much smaller number of monoenergetic gamma-ray simulations and uses interpolation to quickly generate the responses of gamma rays that were not simulated. During the interpolation process, the underlying physics equations are used to accurately compute the response of a given energy gamma ray from the small set of simulations. Such work enables accelerated generation of synthetic radiation detector data.

Detector response↗