Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase field 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 127 records · Page 7

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗

Evaluation of acoustic testing techniques for spacecraft systems

External acoustic environments, structural responses, noise reductions, and the internal acoustic environments have been predicted for a typical shroud/spacecraft system during lift-off and various critical stages of flight. Spacecraft responses caused by energy transmission from the shroud via mechanical and acoustic paths have been compared and the importance of the mechanical path has been evaluated. Theoretical predictions have been compared extensively with available laboratory and in-flight measurements. Equivalent laboratory acoustic fields for simulation of shroud response during the various phases of flight have been derived and compared in detail. Techniques for varying the time-space correlations of laboratory acoustic fields have been examined, together with methods for varying the time and spatial distribution of acoustic amplitudes. Possible acoustic testing configurations for shroud/spacecraft systems have been suggested and trade-off considerations have been reviewed. The problem of simulating the acoustic environments versus simulating the structural responses has been considered and techniques for testing without the shroud installed have been discussed.

Cockburn, J. A.↗

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗

Models of Neptune's smooth recurrent radio emission

The quantitative response of the Planetary Radio Astronomy (PRA) instrument to a wave with polarization ellipse of arbitrary shape and orientation, arriving at the antennas from any direction, can be determined. This capability is used to model the time variation of intensity and circular polarization over a range of radio frequencies for proposed radio-source locations and emission characteristics at Neptune. At frequencies below 400 kHz the observed variation of intensity, polarization, and phase are closely simulated in an offset tilted dipole magnetic field by conjugate sources at midlatitude with filled emission cones. The phase of emission at higher frequencies is reproduced by sources at lower latitude. Modeled wide-cone emission does not reach the spacecraft at the observed phase nor have the polarization sense observed before closest approach. Source-surface maps of apparent polarization for the period before closest approach when instrumental response is especially sensitive to source location is presented. The method is capable of extension to more realistic models of the magnetic field.

Sawyer, Constance↗

Focused Ion Beam Tomography of Alloy 617 Corroded in Molten Chloride Salt

Materials qualification of reactor structural materials is a critical step in rapid implementation of advanced nuclear reactor technologies, particularly to assess the corrosion performance in these designs. Accelerated qualification of reactor structural materials requires incorporating powerful computational toolsets, such as phase field modelling in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to predict the evolution of structural materials due to corrosion. Accordingly, computational toolsets will require experimental data generated at appropriate length scales to validate accuracy. Focused ion beam (FIB) provides a high degree of control over manipulation of materials for analytical purposes, including capturing data on the evolution in the microstructure and elemental composition of materials at the mesoscale, an appropriate length scale for phase field modelling of intergranular diffusion phenomena using the MOOSE framework. For instance, the FEI Helios G4 UX dual beam plasma FIB microscope at the Irradiated Materials Characterization Laboratory (IMCL) is capable of backscatter diffraction (EBSD) and energy-dispersive x-ray spectroscopy (EDS) documenting the evolution in the microstructure and elemental composition, respectively. The Helios can perform EDS and EBSD three-dimensionally (3D) using tomography, which is then combined using different software packages to visualize 3D volumes correlating elemental composition to microstructural data. The purpose of this investigation was to develop a streamlined characterization and data processing workflow for 3D tomography studies on the FEI Helios G4 plasma FIB. The investigation is segmented into three parts: 1) Optimizing the data collection workflow, 2) identifying appropriate data processing and visualization software (i.e. DREAM.3D, MIPAR, and VGStudioMax), and 3) establishing an infrastructure for public release. The optimization of the data collection workflow is in collaboration with members of the U220 department to setup formal training on the tomography operation of the G4, through ThermoFisher Scientific, and exploring DREAM.3D, MIPAR, and VGStudioMax data processing/visualization software packages. VGStudioMax currently demonstrates the most promise for future use. Optimization of the data collection and processing workflow is still ongoing. A collaboration with INL High Performance Computing (HPC) established an open-source license for expediting the public release of FIB tomography datasets through HPC. FIB tomography data generated by the G4 will provide comprehensive data for validating 3D phase field mesoscale modelling tools within the MOOSE framework for accelerated qualification of reactor structural materials.

Copeland-Johnson, Trishelle↗

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. The primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-110 km) for seasonal dates and locations during Phase I aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase I aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars' middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. ne primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-110 km) for seasonal dates and locations during Phase 1 aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase 1 aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars' middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. The primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-1 10 km) for seasonal dates and locations during Phase 1 aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase I aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗

Simulations of Large-Scale Zero Boiloff, Densification, and Solidification of Hydrogen

An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.

Computational Fluid Dynamics↗

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗

Three dimensional modeling of cirrus during the 1991 FIRE IFO 2: Detailed process study

A three-dimensional model of cirrus cloud formation and evolution, including microphysical, dynamical, and radiative processes, was used to simulate cirrus observed in the FIRE Phase 2 Cirrus field program (13 Nov. - 7 Dec. 1991). Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing of solution drops, growth by water vapor deposition, evaporation, aggregation, and vertical transport. Visible and infrared radiative fluxes, and radiative heating rates are calculated using the two-stream algorithm described by Toon et al. Wind velocities, diffusion coefficients, and temperatures were taken from the MAPS analyses and the MM4 mesoscale model simulations. Within the model, moisture is transported and converted to liquid or vapor by the microphysical processes. The simulated cloud bulk and microphysical properties are shown in detail for the Nov. 26 and Dec. 5 case studies. Comparisons with lidar, radar, and in situ data are used to determine how well the simulations reproduced the observed cirrus. The roles played by various processes in the model are described in detail. The potential modes of nucleation are evaluated, and the importance of small-scale variations in temperature and humidity are discussed. The importance of competing ice crystal growth mechanisms (water vapor deposition and aggregation) are evaluated based on model simulations. Finally, the importance of ice crystal shape for crystal growth and vertical transport of ice are discussed.

Jensen, Eric J.↗

Prediction of thermal-stress and deformations due to phase change in solidifying objects via flux/stress based finite element representations

The paper presents numerical simulations for the prediction of thermal-stress and deformation fields resulting from phase change in solidifying bodies employing new finite element representations. The formulations herein demonstrated provide different perspectives and physical interpretation for the modeling/analysis of thermo-mechanical problems and possess several inherent advantages. In comparison to traditional approaches for solving similar problems, the paper employs new flux/stress based representations to enhance the overall effectiveness. Comparative numerical applications validate applicability of the formulations for predicting the temperature induced deformations and stresses resulting from effects due to phase change.

Tamma, K. K.↗

Micro- to macroscale perspectives on space plasmas

The Earth's magnetosphere is the most accessible of natural collisionless plasma environments; an astrophysical plasma 'laboratory'. Magnetospheric physics has been in an exploration phase since its origin 35 years ago but new coordinated, multipoint observations, theory, modeling, and simulations are moving this highly interdisciplinary field of plasma science into a new phase of synthesis and understanding. Plasma systems are ones in which binary collisions are relatively negligible and collective behavior beyond the microscale emerges. Most readily accessible natural plasma systems are collisional and nearest-neighbor classical interactions compete with longer-range plasma effects. Except for stars, most space plasmas are collisionless, however, and the effects of electrodynamic coupling dominate. Basic physical processes in such collisionless plasmas occur at micro-, meso-, and macroscales that are not merely reducible to each other in certain crucial ways as illustrated for the global coupling of the Earth's magnetosphere and for the nonlinear dynamics of charged particle motion in the magnetotail. Such global coupling and coherence makes the geospace environment, the domain of solar-terrestrial science, the most highly coupled of all physical geospheres.

Eastman, Timothy E.↗

Phase-field modeling for restructuring in the dark zone of high burnup UO 2

This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗