SEARCH · Search NASA

Results for “Cluster Dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Irradiation enhanced diffusion and diffusional creep in U₃Si₂

U₃Si₂ is a candidate as an advanced fuel for light water reactors (LWRs). Compared to UO₂, there are significant data gaps for the thermophysical and thermomechanical properties of U₃Si₂. Atomic scale simulations can help fill in those gaps, while also providing a more mechanistic understanding of the material properties. Creep is an important property for fuel performance as it influences the pellet-cladding gap, which in turn affects the temperature of the fuel (a key parameter in many other aspects of fuel behavior). In this work, we employed a DFT-informed cluster dynamics framework to predict point defect concentrations in U₃Si₂. The concentrations, diffusivities, and volumes of these defects were used to develop a diffusional creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusion processes. It was found that irradiation enhanced Nabarro-Herring creep is dominant at low temperatures, while Coble creep dominates for high temperatures. The model compares well against available experimental data and has been implemented in BISON. A demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U₃Si₂, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Implementation and testing of physics-based pulverization model in BISON

This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Development of Mechanistic Fission Gas Release and Swelling Models for UN Fuels in BISON

This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Non-Fermi Liquid Behavior in the Three-Dimensional Hubbard Model

We present a numerical study on the non-Fermi liquid behavior of a three-dimensional strongly correlated system. The Hubbard model in a simple cubic lattice is simulated by the dynamical cluster approximation; in particular, the quasi-particle weight is calculated at finite dopings for a range of temperatures. By fitting the quasi-particle weight to the marginal Fermi liquid form at finite doping near the putative quantum critical point, we find evidence of a separatrix between Fermi liquid and non-Fermi liquid regions. Our results suggest that a marginal Fermi liquid and possibly a quantum critical point exist in the non-symmetry broken solution of the three-dimensional interacting electron systems. We also calculate the spectral function, close to the half-filling, and we obtain evidence of pseudogap.

Kellar, Samuel↗

DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

NASA NTRS1977

On the fragmentation of cosmic gas clouds. III - The initial stellar mass function

Radiation from the first opaque protostellar fragments to form in a nonhomologously collapsing molecular cloud can provide a significant heat source in regions where the matter remains relatively diffuse over the initial free-fall time scale. This energy input provides a negative feedback that gradually inhibits fragmentation as larger protostellar fragments form that are also more effective radiators. A highly simplified model is described that yields an expression for the initial mass function in stellar clusters. Dynamical dissipation, driven by the more massive protostars, tends to decouple the gas and grain temperatures and steepen the mass function. Fragmentation of less massive stars appears to proceed more efficiently when the heavy-element abundance is decreased, indicating the possibility of an inverse correlation between gradients in Z and M/L. Other possible observational tests of fragmentation are also suggested.

Silk, J.↗

NASA NTRS1993

A study of the x ray environment of radio galaxies

We are currently working on a program to use extensive x-ray and radio databases to investigate the relationship between extended radio emission and environment in clusters of galaxies. The radio galaxy morphology is determined using VLA imaging and the x-ray properties are determined from Einstein IPC images. This study is motivated by the hypothesis that the key to understanding radio galaxies lies in the local environment. To test this hypothesis we have studied the detailed relationship between galaxy radio emission and the x-ray morphology of their parent clusters. In this pilot study we have used 35 radio sources found in 27 clusters. We have determined the position angle of the x-ray and radio emission, and x-ray and radio luminosities. The x-ray position was taken to be the position of peak flux of the subclump containing the radio galaxy. The radio position was taken to be the position of the galaxy. We do not find a correlation between the x-ray and radio source position angle. This remains true when the sample is divided into subsamples according to radio morphology (wide angle tail, twin jet, narrow angle tail galaxies). We find a weak correlation between the radio source luminosity and the x-ray luminosity. We have computed the distance from the radio galaxy position to the center of the x-ray clump. We find a mean distance from the x-ray clump center of 0.16 Mpc for the radio galaxies in this sample. The mean distance to the nearest clump of x-ray emission is typically half the distance to the optical cluster center. We thus find strong evidence that radio galaxies are located very close to clumps of x-ray emission. These subclumps are not always affiliated with the central cluster x-ray emission. This supports our hypothesis that x-ray emission may provide a key to understanding radio galaxy morphology. We find evidence that radio galaxies occur in clusters that contain prominent substructures. Radio galaxies may thus provide an added diagnostic of the cluster dynamical state.

Rhee, George F.↗

NASA NTRS2017

Preliminary Analysis of Two Years of the Massive Collision Monitoring Activity

It is hypothesized that the interactions between many of the most massive derelicts in low Earth orbit are more frequent than modeled by the traditional combination of kinetic theory of gases and Poisson probability distribution function. This is suggested by the fact that there are clusters of derelicts where members' inclinations are nearly identical and their apogees/perigees overlap significantly resulting in periodic synchronization of the objects' orbits. In order to address this proposition, an experiment was designed and conducted over the last two years. Results from this monitoring and characterization experiment are presented with implications for proposed debris remediation strategies. Four separate clusters of massive derelicts were examined that are centered around 775km, 850km, 975km, and 1500km, respectively. In aggregate, the constituents of these clusters contain around 500 objects and about 800,000kg of mass; this equates to a third of all derelict mass in LEO. Preliminary analysis indicates that encounter rates over this time period for these objects are greater than is estimated by traditional techniques. Hypothesized dependencies between latitude of encounter, relative velocity, frequency of encounters, inclination, and differential semi-major axis were established and verified. This experiment also identified specific repeatable cluster dynamics that may reduce the cost/risk and enhance the effectiveness of debris remediation activities and also enable new operational debris remediation options.

McKnight, Darren↗

DOE OSTI2021

A Case Study of LLVM-Based Analysis for Optimizing SIMD Code Generation

This paper presents a methodology for using LLVM-based tools to tune the DCA++ (dynamical cluster approximation) application that targets the new ARM A64FX processor. The goal is to describe the changes required for the new architecture and generate efficient single instruction/multiple data (SIMD) instructions that target the new Scalable Vector Extension instruction set. During manual tuning, the authors used the LLVM tools to improve code parallelization by using OpenMP SIMD, refactored the code and applied transformation that enabled SIMD optimizations, and ensured that the correct libraries were used to achieve optimal performance. By applying these code changes, code speed was increased by 1.98× and 78 GFlops were achieved on the A64FX processor. The authors aim to automatize parts of the efforts in the OpenMP Advisor tool, which is built on top of existing and newly introduced LLVM tooling.

Huber, Joseph↗

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↗

DOE OSTI2024

Multi-scale modeling of fuel fragmentation and microstructural evolution

An overview of NEAMS program activities to model fuel fragmentation and microstructual evolution to a group of industry attendees at a meeting hosted by EPRI. Topics include atomistic methods to determine grain boundary strength, phase-field fracture modeling to determine critical pressure at grain boundaries, and coupled phase-field-cluster dynamics to model microstructural evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗

DOE OSTI2026

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

DOE OSTI2021

Examination of Radiation Belt Dynamics During Substorm Clusters: Activity Drivers and Dependencies of Trapped Flux Enhancements

Here, dynamical variations of radiation belt trapped electron fluxes are examined to better understand the variability of enhancements linked to substorm clusters. Analysis is undertaken using the Substorm Onsets and Phases from Indices of the Electrojet substorm cluster algorithm for event detection. Observations from low earth orbit are complemented by additional measurements from medium earth orbit to allow a major expansion in the energy range considered, from medium energy energetic electrons up to ultra-relativistic electrons. The number of substorms identified inside a cluster does not depend strongly on solar wind drivers or geomagnetic indices either before, during, or after the cluster start time. Clusters of substorms linked to moderate (100 nT < AE ≤ 300 nT) or strong AE (AE ≥ 300 nT) disturbances are associated with radiation belt flux enhancements, including up to ultra-relativistic energies by the strongest substorms (as measured by strong southward Bz and high AE). These clusters reliably occur during times of high speed solar winds streams with associated increased magnetospheric convection. However, substorm clusters associated with quiet AE disturbances (AE ≤ 100 nT) lead to no significant chorus whistler mode intensity enhancements, or increases in energetic, relativistic, or ultra-relativistic electron flux in the outer radiation belts. In these cases the solar wind speed is low, and the geomagnetic Kp index indicates a lack of magnetospheric convection. Our study clearly indicates that clusters of substorms occurring outside of high speed wind streams are not by themselves sufficient to drive acceleration, which may be due to the lack of pre-cluster convection.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

Einstein imaging observations of clusters of galaxies

X-ray observations suggest a broad division between clusters having a centrally-located, X-ray dominant galaxy characterized by a centrally-peaked X-ray emission, and clusters with a broader X-ray profile. Within these categories, such evolutionary indicators as X-ray temperature, vorticity dispersion, spiral fraction and central galaxy density may indicate the relative stages of evolution. X-ray imaging observations are presented for (1) the Virgo cluster, (2) other dynamically unevolved clusters, (3) double clusters which may be in an intermediate evolutionary stage, and (4) dynamically evolved clusters. The data are analyzed under the assumption of an isothermal X-ray gas which is gravitationally bound by galaxies and clusters.

Forman, W.↗

NASA NTRS1988

The kinematics and dynamics of the rich cluster of galaxies Abell 539

A sample of 289 redshifts within 10 deg of the center of the cluster A539 have been obtained in order to investigate the cluster kinematics and dynamics. Within 1 Mpc of the center, the physical parameters of A539 are found to be typical of those of rich clusters. It is shown that early-type galaxies are more concentrated toward the cluster center and that the velocity distributions of early-type and late-type galaxies differ marginally.

Ostriker, Eve C.↗

NASA NTRS1984

Constraints on the anisotropy of the velocity dispersion of the Coma cluster

The dynamics of the Coma cluster are probed through the study of spiral galaxies that have been stripped of gas and ellipticals that have been tidally stripped. Since such stripping is most effective in regions where intergalaxy or cluster gas density is highest, the affected galaxies probably have orbits that take them through the center of the cluster. The Bothun et al. (1982, 1983) data analyzed indicate that the spirals in the cluster center are deficient in H I relative to other spiral members. Observations by Strom and Strom (1978) may indicate the tidal stripping of elliptical galaxies near the cluster center. For both data sets, the most likely models are everywhere isotropic.

Pryor, C.↗