Working towards a modular, fully-coupled phase field fracture model integrating elasticity, plasticity, and damage
Explore the source record for details and available documents.
Nataraj, Chiraag M [Sandia National Laboratories (↗
SEARCH · Engineering Papers
Results for “Phase Field Model”
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.
Working towards a modular, fully-coupled phase field fracture model integrating elasticity, plasticity, and damage
Explore the source record for details and available documents.
Nataraj, Chiraag M [Sandia National Laboratories (↗
Niobium hydride formation in superconducting qubit thin films
The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.
36 - MATERIALS SCIENCE↗
Holistic Microstructure Control Strategies in Photopolymerization‐Induced Phase Separation of Acrylate Systems
Open porous materials, known for their large surface area and interconnected structures, are essential in various applications, including batteries, ion exchange, catalysis, filtration, and electronic waste recycling. A critical aspect of the functionality of porous membranes is the precise control of pore size and morphology. Photopolymerization-induced phase separation (photo-PIPS) offers a convenient and versatile methods for creating porous structures. However, controlling the porous morphology remains challenging due to the complex interplay between thermodynamics, polymerization kinetics, and monomer structures, which makes it difficult to establish the relationship between processing conditions and resulting morphology in photo-PIPS. Herein, a physics-based phase-field model capable of generating and characterizing the microstructures of porous materials based on both average and localized features is developed. Using the phase-field simulations as test bed, the effects of polarity, light intensity, and curing temperature, as well as the previously unexplored roles of chain transfer agents and substrates, on the morphology of the resulting porous microstructure are investigated. Experiments are performed to verify the results predicted by the simulations. This work lays out a comprehensive guide for designing PIPS-derived porous microstructures and offers practical engineering strategies for tailoring microstructure-level topology and size of pores for application-specific needs.
36 MATERIALS SCIENCE↗
Mechanistic FCCI model for BISON
In this report, multi-scale simulations to develop a mechanistic model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY23 are described. The model assumes that cladding wastage is dominated by the intermetallic phase (Fe,Cr)$_{17}$Nd$_2$, and calculates production and diffusion of Nd in the fuel, diffusion into the cladding, and formation and growth of the intermetallic wastage layer. Atomistic simulations are used to calculate diffusivity of Nd in the (Fe,Cr) system. Mesoscale simulations, informed by previous atomistic calculations, are used to determine the effective diffusivity of Nd through the fuel, accounting for empty and sodium-filled porosity within the fuel. The BISON model incorporates this data and calculates the wastage layer thickness as a function of time, assuming parabolic growth kinetics based on the best available data. The results of the new mechanistic BISON model are compared to existing assessment cases for the X447 experiment in EBR-II, and show good agreement with experimental data and the existing BISON empirical model. To improve predictions of kinetics for future work, a phase-field model that includes fuel, wastage, and cladding phases is developed. Preliminary testing indicates that the kinetics may differ significantly from the assumed parabolic growth law, suggesting direction for future improvement of the model that will allow it to be applied to other fuel designs with greater confidence.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Nanomechanics of Ferroelectric Thin Films and Heterostructures
The focus of this chapter is to provide basic concepts of how external strains/stresses altering ferroelectric property of a material and how to evaluate quantitatively the effect of strains/stresses on phase stability, domain structure, and material ferroelectric properties using the phase-field method. The chapter starts from a brief introduction of ferroelectrics and the Landau-Devinshire description of ferroelectric transitions and ferroelectric phases in a homogeneous ferroelectric single crystal. Due to the fact that ferroelectric transitions involve crystal structure change and domain formation, strains and stresses can be produced inside of the material if a ferroelectric transition occurs and it is confined. These strains and stresses affect in turn the domain structure and material ferroelectric properties. Therefore, ferroelectrics and strains/stresses are coupled to each other. The ferroelectric-mechanical coupling can be used to engineer the material ferroelectric properties by designing the phase and structure. The followed section elucidates calculations of the strains/stresses and elastic energy in a thin film containing a single domain, twinned domains to complicated multidomains constrained by its underlying substrate. Furthermore, a phase field model for predicting ferroelectric stable phases and domain structure in a thin film is presented. Examples of using substrate constraint and temperature to obtain interested ferroelectric domain structures in BaTiO3 films are demonstrated b phase field simulations.
Ferroelectrics↗
Effect of Recrystallization on Gas Bubble Swelling in UMo Fuels
The development of low-enriched uranium (LEU) fuels such as UMo dispersion fuels and monolithic UMo fuels is part of a global effort on nuclear nonproliferation. Although UMo fuels have the advantage of high density, excellent irradiation performance and good thermal conductivity, the volumetric swelling in irradiated metallic fuels is known to be an important design parameter because it affects not only the thermal conductivity of the fuels but also the mechanical integrity of fuel structures. Experiments also show that recrystallization dramatically affects the gas bubble swelling kinetics. A series of materials processes such as casting, mechanical rolling and thermal treatment are used in fuel fabrication. Resulting microstructure, such as grain size, porosity, and second phase, strongly depends on the material processes and process parameters. Therefore, a fundamental understanding how fabricated microstructures and irradiation conditions affect the nucleation of gas bubbles, recrystallization and gas bubble swelling kinetics is essential for optimizing material processes to obtain desired microstructure and fuel performance. Great effort has been made in the development of modeling capability for predicting gas bubble swelling kinetics and fission gas release kinetics in nuclear fuels. However, these models are based on the assumptions that the inter- and intra-granular gas bubble sizes are independent of position and only depend on time or fission density. Therefore, they are applicable to a material with uniform microstructure such as equiaxed grains with homogeneous defect distributions, homogeneous thermodynamic and kinetic properties of defects, and uniform radiation conditions. Very recently, we developed a three dimension (3D) microstructure dependent gas bubble swelling model for investigating the effect of heterogeneous microstructures and inhomogeneous thermodynamic properties of defects on gas bubble swelling in polycrystalline UMo fuels. However, the model doesn't take into account the effect of recrystallization on swelling kinetics. In this work, the spatially dependent rate theory is used to describe the evolution of interstitials and interstitial loop density. Recrystallization is described by a phase-field model with a concept of 'phase transition' when the local interstitial loop density is larger than a critical value. The intra- and inter- granular gas bubble evolution is described by the spatially dependent Booth model. A model integrating all the microstructure evolution enables one to study the effect of heterogeneous microstructures, recrystallization, and the metallic fuel thermodynamic and kinetic properties on gas bubble swelling kinetics.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Atomistic modeling of LiF microstructure ionic conductivity and its influence on nucleation and plating
We report the formation and degradation of the solid electrolyte interphase (SEI) and its underlying transport properties play an essential role in the overall performance of lithium-ion batteries. This paper presents classical molecular dynamics studies on polycrystalline inorganic lithium fluoride (LiF) layers to model and predict the SEI transport properties. The ionic conductivity is obtained from the lithium-ion diffusivity in polycrystalline structures of LiF using the Nernst-Einstein relation. The predicted molecular dynamics data are used in a continuum scale phase-field model to evaluate the plating kinetics under fast charging conditions. The analysis emphasizes that the SEI ionic conductivity properties impact the plating dynamics, where SEI's low ion conductivity value is prone to large plating and subsequent capacity degradation. The combination of atomic and continuum scale studies shown herein lays a foundation to tune in SEI transport properties to decrease the amount of lithium plating and improve the performance of fast-charging batteries.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Effect of substrate orientation and anisotropic strength on corrosion pits
The complexities of pitting corrosion and its associated pit morphologies are critical for predicting and preventing pit-related failure in materials. With the advances in computational tools, modeling of pitting corrosion has lately gained momentum. We recently developed a phase field model for 3-D pit morphology, and in the current work, the model is utilized to gain insight into the role of substrate orientation and anisotropic strength on crystallographic pits. The influence of these parameters on pit morphology and corrosion metrics such as pit volume, current, pit area, current density, and pit depth are studied. The crystallographic pit morphology changes with substrate orientation but not significantly with anisotropic strength. Finally, the anisotropic strength strongly influences the corrosion metrics whereas most corrosion metrics remain unaffected by substrate orientation.
36 MATERIALS SCIENCE↗
Overview of the Tusas Code for Simulation of Dendritic Solidification
The aim of this project is to conduct a parametric investigation into the modeling of two dimensional dendrite solidification, using the phase field model. Specifically, we use the Tusas code, which is for coupled heat and phase-field simulation of dendritic solidification. Dendritic solidification, which may occur in the presence of an unstable solidification interface, results in treelike microstructures that often grow perpendicular to the rest of the growth front. The interface may become unstable if the enthalpy of the solid material is less than that of the liquid material, or if the solute is less soluble in solid than it is in liquid, potentially causing a partition [1]. A key motivation behind this research is that a broadened understanding of phase-field formulation and microstructural developments can be utilized for macroscopic simulations of phase change. This may be directly implemented as a part of the Telluride project at Los Alamos National Laboratory (LANL), through which a computational additive manufacturing simulation tool is being developed, ultimately to become part of the Advanced Simulation and Computing Program within the U.S. Department of Energy [2].
97 MATHEMATICS AND COMPUTING↗
Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction
Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.
BISON↗
Linking Initial Microstructure to ORR Related Property Degradation in SOFC Cathode: A Phase Field Simulation
Microstructure evolution driven by thermal coarsening is an important factor for the loss of oxygen reduction reaction rates in SOFC cathode. In this work, the effect of an initial microstructure on the microstructure evolution in SOFC cathode is investigated using a recently developed phase field model. Specifically, we tune the phase fraction, the average grain size, the standard deviation of the grain size and the grain shape in the initial microstructure, and explore their effect on the evolution of the grain size, the density of triple phase boundary, the specific surface area and the effective conductivity in LSM-YSZ cathodes. It is found that the degradation rate of TPB density and SSA of LSM is lower with less LSM phase fraction (with constant porosity assumed) and greater average grain size, while the degradation rate of effective conductivity can also be tuned by adjusting the standard deviation of grain size distribution and grain aspect ratio. The implication of this study on the designing of an optimal initial microstructure of SOFC cathodes is discussed.
25 ENERGY STORAGE↗
Grizzly Staus Report
This report summarizes work during FY 2014 to develop capabilities to predict embrittlement of reactor pressure vessel steel, and to assess the response of embrittled reactor pressure vessels to postulated accident conditions. This work has been conducted a three length scales. At the engineering scale, 3D fracture mechanics capabilities have been developed to calculate stress intensities and fracture toughnesses, to perform a deterministic assessment of whether a crack would propagate at the location of an existing flaw. This capability has been demonstrated on several types of flaws in a generic reactor pressure vessel model. Models have been developed at the scale of fracture specimens to develop a capability to determine how irradiation affects the fracture toughness of material. Verification work has been performed on a previously-developed model to determine the sensitivity of the model to specimen geometry and size effects. The effects of irradiation on the parameters of this model has been investigated. At lower length scales, work has continued in an ongoing to understand how irradiation and thermal aging affect the microstructure and mechanical properties of reactor pressure vessel steel. Previously-developed atomistic kinetic monte carlo models have been further developed and benchmarked against experimental data. Initial work has been performed to develop models of nucleation in a phase field model. Additional modeling work has also been performed to improve the fundamental understanding of the formation mechanisms and stability of matrix defects caused.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Phase-field based Multiscale Modeling of Heterogeneous Solid Electrolytes: Applications to Nanoporous Li 3 PS 4
Modeling the effective ion conductivities of heterogeneous solid electrolytes typically involves the use of a computer-generated microstructure consisting of randomly or uniformly oriented fillers in a matrix. But, the structural features of the filler/matrix interface, which critically determine the interface ion conductivity and the microstructure morphology, have not been considered during the microstructure generation. In using nanoporous β-Li 3 PS 4 electrolyte as an example, we develop a phase-field model that enables generating nanoporous microstructures of different porosities and connectivity patterns based on the depth and the energy of the surface (pore/electrolyte interface), both of which are predicted through density functional theory (DFT) calculations. Room-temperature effective ion conductivities of the generated microstructures are then calculated numerically, using DFT-estimated surface Li-ion conductivity (3.14×10 -3 S/cm) and experimentally measured bulk Li-ion conductivity (8.93×10 -7 S/cm) of β-Li 3 PS 4 as the inputs. We also use the generated microstructures to inform effective medium theories to rapidly predict the effective ion conductivity via analytical calculations. Furthemore, when porosity approaches the percolation threshold, both the numerical and analytical methods predict a significantly enhanced Li-ion conductivity (1.74×10 -4 S/cm) that is in good agreement with experimental data (1.64×10 -4 S/cm). The present phase-field based multiscale model is generally applicable to predict both the microstructure patterns and the effective properties of heterogeneous solid electrolytes.
36 MATERIALS SCIENCE↗
Benchmark problems for the Me soscale M ultiphysics Ph ase F i eld S imulator (MEMPHIS)
This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Me soscale M ultiphysics Ph ase F i eld S imulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported 'classical' phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Irradiation induced void spheroidization, shrinkage and migration in Cu at elevated temperatures: An in situ study
We show understanding the void evolution in irradiation environment is of great interest and significance, as irradiation-induced voids typically lead to pronounced volumetric swelling and degradation of mechanical properties. In situ studies on the irradiation response of nanovoids at elevated temperature remain limited. In this work, we performed systematic in situ 1 MeV Kr ++ irradiations on Cu with nanovoids in a transmission electron microscope up to 350 °C. The in situ studies revealed intriguing void spheroidization, shrinkage and migration. Furthermore, the morphology evolution and migration of nanovoids showed a strong dependence on irradiation temperature and initial void size. Post-irradiation analyses identified defect clusters in the form of stacking fault tetrahedrons, and the remaining large faceted nanovoids. The underlying mechanisms of irradiation-induced void spheroidization and shrinkage were discussed based on phase-field modeling.
36 MATERIALS SCIENCE↗