Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster Dynamics 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Comparison of dynamical approximation schemes for non-linear gravitational clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.

Melott, Adrian L.↗
DOE OSTI2023

Multiscale modeling for high burnup structure formation and associated pulverization

This report summarizes the lower length scale modeling work performed in the fiscal year 2023 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstrutural evolution and associated pulverization criteria in the high burnup regions. This year, the resolution mechanisms within the cluster dynamics code has been updated and it influence on the bubble growth in HBS regions at the mesoscale is studied. To improve the accuracy of phase-field models of fission gas bubble growth, a new Helmholtz free energy for high-density Xe gas is derived from a virial equation of state. Two previously used strategies for representing net vacancy production in phase-field models of fission gas bubble growth are compared with each other and with analytical models of bubble growth. The vacancy source only model is found to be more convenient to parameterize realistically compared with the vacancy source+sink model. Furthermore, the phase-field-fracture simulation have been performed using MD- informed failure stress values and realistic HBS structure obtained from the phase-field simulations. We also present the uncertainty bands on prediction of the critical stress to account for the effect of the lower length scale variabilities on the failure criteria at the mesoscale. It is observed that pulverization may occur in partially restructured regions with restructuring fraction as low as 17%. In light of this, an update to the BISON’s pulverization criteria is recommended.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

MEDPRAT Treatment Clusters: Improving Representation of Mission Medical Risk

INTRODUCTION The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) implements a computational model that aims to quantify spaceflight medical risk by utilizing probabilistic techniques to simulate critical event incidence and outcomes over thousands of simulated mission trials. The goal of MEDPRAT is to characterize mission medical risk and provide insight into medical resource utilization. In order to analyze the medical resource space, treatment must be mapped from each simulated condition, and resources consumed as a result of this treatment must be tracked throughout the course of the mission. A new MEDPRAT feature, ‘treatment clusters’, provide a more sophisticated method of defining the structure and interaction between resources, more closely mimicking the way treatment is carried out clinically. METHODS Treatment clusters expand on the two existing treatment groupings (combination and alternate) adding a new grouping: bundled treatment. Treatment clusters may be combined to any depth, giving users the ability to specify complex treatment trees whose behavior is governed by several user-specified parameters. This approach emphasizes reusability, as treatment clusters, once defined, can be used to create more complex treatment trees or applied to many conditions. By configuring parameters for contribution, efficacy, necessity, primacy, and equivalence, resource relationships and dependencies can be more accurately represented, thereby allowing users to build capabilities with desired treatment properties, for example an intravenous capability for conditions such as anaphylaxis, acute radiation syndrome, etc. MEDPRAT v1.0 remains backward compatible with existing treatment structures, giving users the ability to define new treatment clusters as evidence becomes available, without having to recode their existing treatment databases. In addition to facilitating the representation of more complex treatment options, by pairing treatment clusters with the internal optimization routine, the MEDPRAT set selector, medical resources can be identified as organized in bundles, where appropriate, so that optimized resource sets include groups of highly-dependent resources only when all resources of the group are together. For example, it would be wasteful to include ultrasound gel but not an ultrasound machine, since the gel provides no benefit as a treatment without the ultrasound machine. With treatment clusters, the user may require that both resources are available to provide any benefit as treatment, so that if one resource is optimized out of the set, the other resource will be optimized out as well. RESULTS AND CONCLUSIONS We will report on MEDPRAT treatment clusters used in a bundling study under the IMPACT project of the ExMC element. We will discuss an example of a complex treatment tree. Through the implementation of this feature MEDPRAT enables treatment to be defined and applied in a way that is more representative of the real world, providing more accurate insight into mission medical risk and the medical resource space.

Lawrence Leinweber↗
DOE OSTI2025

Structure of human MUTYH and functional profiling of cancer-associated variants reveal an allosteric network between its [4Fe-4S] cluster cofactor and active site required for DNA repair

Abstract MUTYH is a clinically important DNA glycosylase that thwarts mutations by initiating base-excision repair at 8-oxoguanine (OG):A lesions. The roles for its [4Fe-4S] cofactor in DNA repair remain enigmatic. Functional profiling of cancer-associated variants near the [4Fe-4S] cofactor reveals that most variations abrogate both retention of the cofactor and enzyme activity. Surprisingly, R241Q and N238S retained the metal cluster and bound substrate DNA tightly, but were completely inactive. We determine the crystal structure of human MUTYH bound to a transition state mimic and this shows that Arg241 and Asn238 build an H-bond network connecting the [4Fe-4S] cluster to the catalytic Asp236 that mediates base excision. The structure of the bacterial MutY variant R149Q, along with molecular dynamics simulations of the human enzyme, support a model in which the cofactor functions to position and activate the catalytic Asp. These results suggest that allosteric cross-talk between the DNA binding [4Fe-4S] cofactor and the base excision site of MUTYH regulate its DNA repair function.

Science & Technology - Other Topics↗
DOE OSTI2024

Molecular insights into the in situ early-stage assembly of metal–organic frameworks on cellulose nanofibrils

The integration of metal–organic frameworks (MOFs) into renewable matrices, including those derived from cellulose, to create cellulose/MOF hybrids has attracted significant interest due to the synergistic combination of cellulosic biopolymer properties and MOFs' multifunctional features. The interfacial interactions between cellulose and MOF not only ensure stable bonding, but ultimately also determine the bulk physical properties of the macroscopic composites. However, the mechanistic understanding of the in situ assembly between the two materials remains unclear. In this study, we employ a combination study of synthesis, experimental characterization, and molecular dynamics simulations to explore the early-stage assembly of a cellulose/MOF hybrid. Our results revealed that the growth of MOF clusters on the 2,2,6,6-tetramethyl-1-piperidinyloxy oxidized cellulose nanofibrils (TOCNF) follows an inhomogeneous sequential transformation pathway. The carboxylates of TOCNF form coordination-like bonds with the metal ions, while the hydroxyl groups of TOCNF form hydrogen bonds with MOF ligands. These interactions provide the initial nucleation sites that mediate the growth of MOF clusters and also guide the assembly of larger MOF clusters onto the TOCNF substrate. In conclusion, the fundamental insights into the in situ assembly of MOF nanoparticles on cellulosic substrates are essential for the rational design of high-performance materials with tailored morphology and optimized properties.

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

Searching for the Most Harmful Field Errors in the HSR IR Superconducting Magnets

In this project, we improve beam stability for the Electron-Ion Collider. Magnetic field errors can reduce beam stability, making it essential to identify the field errors that have the greatest impact on accelerator performance. However, this is particularly challenging because beam stability depends on the complex interactions of many magnetic field errors, resulting in a high-dimensional and nonlinear optimization problem. We determine which field errors are the most influential for the large physical aperture superconducting magnet B2PF, a critical magnet in the Interaction Region (IR) in the Hadron Storage Ring (HSR). We complete and analyze nearly 30,000 simulations on the Brookhaven National Laboratory Linux Cluster by varying 18 nonlinear magnetic field errors. We evaluate beam stability using the dynamic aperture and the tune diffusion. We identify the field errors that most strongly influence beam stability and establish quantitative field error tolerances that improve accelerator performance.

43 PARTICLE ACCELERATORS↗
NASA NTRS1991

Gravitational clustering in the expanding universe - Controlled high-resolution studies in two dimensions

Results are presented from a series of gravitational clustering simulations in two dimensions. These simulations are a significant departure from previous work, since in two dimensions one can have large dynamic range in both length scale and mass using present computer technology. Controlled experiments were conducted by varying the slope of power-law initial density fluctuation spectra and varying cutoffs at large k, while holding constant the phases of individual Fourier components and the scale of nonlinearity. Filaments are found in many different simulations, even with pure power-law initial conditions. By direct comparison, filaments, called 'second-generation pancakes' are shown to arise as a consequence of mild nonlinearity on scales much larger than the correlation length and are not relics of an initial lattice or due to sparse sampling of the Fourier components. Bumps of low amplitude in the two-point correlation are found to be generic but usually only statistical fluctuations. Power spectra are much easier to relate to initial conditions, and seem to follow a simple triangular shape (on log-log plot) in the nonlinear regime. The rms density fluctuation with Gaussian smoothing is the most stable indicator of nonlinearity.

Beacom, John Francis↗
NASA NTRS2001

A Hardware-in-the-Loop Testbed for Spacecraft Formation Flying Applications

The Formation Flying Test Bed (FFTB) at NASA Goddard Space Flight Center (GSFC) is being developed as a modular, hybrid dynamic simulation facility employed for end-to-end guidance, navigation, and control (GN&C) analysis and design for formation flying clusters and constellations of satellites. The FFTB will support critical hardware and software technology development to enable current and future missions for NASA, other government agencies, and external customers for a wide range of missions, particularly those involving distributed spacecraft operations. The initial capabilities of the FFTB are based upon an integration of high fidelity hardware and software simulation, emulation, and test platforms developed at GSFC in recent years; including a high-fidelity GPS simulator which has been a fundamental component of the Guidance, Navigation, and Control Center's GPS Test Facility. The FFTB will be continuously evolving over the next several years from a too[ with initial capabilities in GPS navigation hardware/software- in-the- loop analysis and closed loop GPS-based orbit control algorithm assessment to one with cross-link communications and relative navigation analysis and simulation capability. Eventually the FFT13 will provide full capability to support all aspects of multi-sensor, absolute and relative position determination and control, in all (attitude and orbit) degrees of freedom, as well as information management for satellite clusters and constellations. In this paper we focus on the architecture for the FFT13 as a general GN&C analysis environment for the spacecraft formation flying community inside and outside of NASA GSFC and we briefly reference some current and future activities which will drive the requirements and development.

Leitner, Jesse↗
DOE OSTI2023

Insights into radiation resistance of titanium alloys from displacement cascade simulations

Radiation damage in beam window materials limits the use of high power proton beams in high energy physics research. The alloy Ti-6Al-4V is presently used as a beam window material but a prospective alternative, Ti-15V-3Cr-3Sn-3Al has been proposed. Since both these alloys contain dual phases at room temperature, we compare the radiation damage in the α and ß-phases of these two materials via primary knock-on atom (PKA) cascade simulations in the 10-40keV energy range. At PKA energies 30 and 40keV, the number of Frenkel pairs in the ballistic stage is higher in the ß-phase of Ti-6Al-4V than that in the ß-phase of Ti-15V-3Cr-3Sn-3Al almost by a factor of 2. The α-phase, of both these alloys, by far outperforms the ß-phases of the two alloys, both in terms of damage during the ballistic stage and in terms of the surviving defects. The average displacement threshold energy (E d ) in the α-phase of both alloys was found to be 66eV while that in the ß-phase of Ti-15-3 was 55 eV and in the ß-phase of Ti-6-4 was 46. Further, while the number of surviving defects is almost equal in both alloys, the vacancy and interstitial clustering mechanisms differ notably, which can impact the degree of radiation hardening and loss of ductility. Our simulations show larger vacancy and interstitial clusters form in Ti-6Al-4V as compared to that in Ti-15-3-3-3 alloy. These results indicate that Ti-15-3-3-3 alloys may be a promising candidate for next generation beam window material with a higher radiation tolerance than the existing Ti-6-4 alloy.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Unsteady evaporation and combustion of a drop cluster inside a vortex

A model has been developed which describes the evaporation, ignition and combustion of a drop cluster embedded in a large vortex. The purpose of this model is to simulate the behavior of drops in large coherent vortices produced in the shear layer of a jet. The model treats the dynamic interactions between the drops and the vortex, and also takes into account the drop proximity to calculate the heat and mass transfers between drops and ambient gas. The gas phase outside the cluster is treated as an unsteady, reacting phase, whereas quasi-steadiness is assumed between the drops and surrounding gas inside the cluster. The results show a very complex interaction between the dynamics of the drop-loaded vortex, the flame, and the evaporation process.

Fichot, F.↗
DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The cosmology dependence of galaxy clustering and lensing from a hybrid N -body–perturbation theory model

ABSTRACT We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from N-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of N-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these non-linear spectra from redshifts z = 0 to z = 2. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h \,{\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift z ≃ 0.4 to scales of $k_{\rm max} \approx 0.6\, h\,\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy–galaxy lensing analyses with data drawn from an independent N-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming data sets from galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Using Atom Dynamics to Map the Defect Structure Around an Impurity in Nano-Hematite

The geochemical cycling of actinides such as U is strongly controlled by Fe (oxyhydr)oxides. Indeed, the strong affinity of these minerals for U has factored into the design of geologic repositories for the long-term storage of nuclear waste. Many decades of work have focused on detailing the local coordination environment U associated with hematite (a-Fe2O3), a common and stable Fe(III) oxide, in order to better predict the fate and transport of U in the sub-surface. Here, for the first time, molecular dynamic simulations and atomically resolved electron microscopy of uranium atom mobility under the beam were used to map the topology of defect clusters surrounding structurally incorporated U in hematite. The ability to observe vacancies by means of metal atom hopping can be used to directly probe defects in other materials, with promising applications in geochemistry and materials science.

Ilton, Eugene S.↗
DOE OSTI2021

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Dynamic Simulation and Stability Analysis

Dynamic Interaction Simulation of Controls and Structure (DISCOS) program developed for dynamic simulation and stability analysis of passive and actively controlled spacecraft. In use of DISCOS, physical system undergoing analysis generally described as cluster of contiguous flexible structures (bodies) that comprise mechanical system, such as spacecraft.

Frisch, H. P.↗