SEARCH · Search NASA

Results for “phase field simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Universal phase dynamics in VO 2 switches revealed by ultrafast operando diffraction

Understanding the pathways and time scales underlying electrically driven insulator-metal transitions is crucial for uncovering the fundamental limits of device operation. Using stroboscopic electron diffraction, we perform synchronized time-resolved measurements of atomic motions and electronic transport in operating vanadium dioxide (VO 2 ) switches. We discover an electrically triggered, isostructural state that forms transiently on microsecond time scales, which is shown by phase-field simulations to be stabilized by local heterogeneities and interfacial interactions between the equilibrium phases. This metastable phase is similar to that formed under photoexcitation within picoseconds, suggesting a universal transformation pathway. Our results establish electrical excitation as a route for uncovering nonequilibrium and metastable phases in correlated materials, opening avenues for engineering dynamical behavior in nanoelectronics.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Precipitate growth kinetics under inhomogeneous concentration fields using a phase-field model

We investigate precipitation dynamics in the presence of a local solute gradient using phase-field simulations. During the homogenization heat treatment of the solidified Inconel 718 alloy, high Nb concentration within the Laves phases or at the core of the secondary arms results in Nb diffusion into the γ matrix. The volume fraction and spatial distribution of precipitation during subsequent annealing can be controlled by tailoring the Nb concentration gradient in the matrix during homogenization. We use a surrogate Ni-Fe-Nb alloy for Inconel 718 to explore the growth dynamics of δ precipitates related to the local Nb concentration levels. The simulations indicate that in the presence of a Nb concentration gradient the growth rate of δ precipitates is higher than in a matrix of uniform average Nb concentration. The higher growth rate is a result of the higher local thermodynamic driving force at the interface between the solute-rich matrix and the δ interface. Here, we propose a phenomenological model to describe the diffusion-controlled growth kinetics of the δ phase under a solute concentration gradient.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coarsening dynamics of ternary polymer solutions with mobility and viscosity contrasts

Using phase-field simulations, we investigate the bulk coarsening dynamics of ternary polymer solutions undergoing a glass transition for two models of phase separation: diffusion only and with hydrodynamics. The glass transition is incorporated in both models by imposing mobility and viscosity contrasts between the polymer-rich and polymer-poor phases of the evolving microstructure. For microstructures composed of polymer-poor clusters in a polymer-rich matrix, the mobility and viscosity contrasts significantly hinder coarsening, effectively leading to structural arrest. For microstructures composed of polymer-rich clusters in a polymer-poor matrix, the mobility and viscosity contrasts do not impede domain growth; rather, they change the transient concentration of the polymer-rich phase, altering the shape of the discrete domains. This effect introduces several complexities to the coarsening process, including percolation inversion of the polymer-rich and polymer-poor phases—a phenomenon normally attributed to viscoelastic phase separation.

Chemistry↗

DOE OSTI2020

Lower length scale informed improvements to Bison U-Pu-Zr fuel swelling model

Due to renewed interest in metallic fuels from the U-(Pu)-Zr material system, significant improvements have recently been made to BISON's capability to simulate metallic fuel, particularly in the area of swelling. In this report, lower length scale simulations that have been conducted to inform the engineering-scale models in BISON during FY20 are described. For the high-temperature regions of the fuel that are dominated by the $\gamma$ phase of U-(Pu)-Zr, two new models for swelling have implemented in BISON. Phase-field simulations were conducted to calculate parameters for these and other BISON models that control when gaseous swelling ceases and fission gas release begins. For the low-temperature regions of the fuel dominated by the phase of uranium, microstructurally resolved simulations of pore growth and crystal deformation were performed to understand the effect of irradiation-generated crystal shape changes. Conformally meshed microstructures with pores were simulated with a crystal plasticity model and with a burnup eigenstrain model. Plastic deformation will not cause pore size to increase; burnup eigenstrain (irradiation-induced crystal shape changes) will increase the volume of existing porosity and cause it to become anisotropic in shape. Additional model development work for early-stage fission gas cluster/bubble formation in the $\gamma$ phase were performed, by coupling with defect clustering within cluster dynamics framework. Sensitivity analysis were conducted to identify the critical thermo-kinetic parameters needed for rigorous assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Hierarchical microstructures: a potential route to enhanced stability in structural materials for advanced nuclear reactors

The drive to increase efficiency in nuclear energy systems is leading to the need for materials that operate at higher temperatures and stress levels for extended periods, while maintaining stable microstructures to ensure their performance is not compromised. A novel route to producing materials that can perform well in such environments is the creation of hierarchical microstructures. A hierarchical microstructure is a microstructure in which features are present at multiple length scales simultaneously. In this work, a hierarchical microstructure is fabricated in a nickel-base superalloy, featuring nanometer-size gamma precipitates inside larger gamma-prime particles, which are in turn embedded in the gamma matrix phase. The hierarchical features of the microstructure lead to enhanced stability of the gamma-prime precipitates during annealing; the particle size does not follow the expected t^(1/3) growth law predicted by the classic LSW theory. Phase-field simulations are used to understand the unexpected stability of the gamma-prime precipitates.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Strain effects on stability of topological ferroelectric polar configurations in (PbTiO3) n /(SrTiO3) n superlattices

The (PbTiO3)n/(SrTiO3)n (PTO/STO) superlattice system has been shown to exhibit interesting topological phases (e.g., vortices and skyrmions) in addition to normal ferroelectric domain states. Existing studies are mostly focused on the dependence of topological polar distributions and properties of PTO/STO superlattice on its periodicity. Here, we study the strain effect on the topological phase transitions and ferroelectric domain structures employing phase-field simulations. We summarized in an isotropic strain (in-plane misfit strain along the x direction is equal to that along the y direction) periodicity phase diagram displaying the stability regions of different polar topological states, including normal ferroelectric twins, vortices, skyrmions, and mixtures of vortices and twins. We also analyzed the polarization configurations under anisotropic in-plane strains (in-plane misfit strain along the x direction is not equal to that along the y direction) and demonstrated that the strain anisotropy can be used to tune the directions of vortex arrays along either the [100] pc or [010] pc directions or labyrinth vortex arrays. This work offers guidance to manipulating polar structures in the PTO/STO superlattices via strain engineering.

Physics↗

DOE OSTI2024

Strain phase equilibria and phase‐field method of ferroelectric polydomain: A case study of monoclinic K x Na 1 − x NbO 3 thin films

Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Hierarchical microstructures: a potential route to enhanced stability in structural materials for advanced nuclear reactors

36 - MATERIALS SCIENCE↗

DOE OSTI2020

Induced ferroelectric phases in SrTiO 3 by a nanocomposite approach

Inducing new phases in thick films via vertical lattice strain is one of the critical advantages of vertically aligned nanocomposites (VANs). In SrTiO 3 (STO), the ground state is ferroelastic, and the ferroelectricity in STO is suppressed by the orthorhombic transition. In this study, we explore whether vertical lattice strain in three-dimensional VANs can be used to induce new ferroelectric phases in SrTiO 3 :MgO (STO:MgO) VAN thin films. The STO:MgO system incorporates ordered, vertically aligned MgO nanopillars into a STO film matrix. Strong lattice coupling between STO and MgO imposes a large lattice strain in the STO film. We have investigated ferroelectricity in the STO phase, existing up to room temperature, using piezoresponse force microscopy, phase field simulation and second harmonic generation. We also serendipitously discovered the formation of metastable TiO nanocores in MgO nanopillars embedded in the STO film matrix. Our results emphasize the design of new phases via vertical epitaxial strain in VAN thin films.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Austenite Formation and Manganese Partitioning during Double Soaking of an Ultralow Carbon Medium‐Manganese Steel

Double soaking (DS) is a thermal processing route intended to produce austenite–martensite microstructures in steels containing austenite‐stabilizing additions and consists of intercritical annealing (primary soaking), followed by heating and brief isothermal holding at an increased temperature (secondary soaking), and quenching. Herein, experimental dilatometry during DS of a medium‐manganese (Mn) steel with nominally 7 wt% Mn and an ultralow residual carbon concentration, in combination with phase‐field simulations of austenite formation during secondary soaking, is presented. The feasibility of maintaining heterogeneous Mn distributions during DS is demonstrated and insight is provided on the effects of the secondary soaking temperature and prior Mn distribution on the ferrite‐to‐austenite phase transformation during the secondary soaking portion of the DS treatment.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Asymmetry of the Ferroelectric Phase Transition in BaTiO 3

Symmetry changes during phase transformations fundamentally determine the behavior of thermodynamic systems, governing phenomena as diverse as water evaporation, fermion condensation, and epidemic spreading. Phase transitions are conventionally divided into two classes, and this classification is typically assumed to remain invariant upon reversing the transition. Here, an asymmetric phase transformation is uncovered in the ferroelectric–paraelectric transition of single-crystal BaTiO 3 , a model system for first-order transitions. Under slow temperature variation (≤0.1 °C min −1 ), thermodynamic, dielectric, and domain-structure measurements reveal that the ferroelectric-to-paraelectric transition exhibits latent heat, phase coexistence, and a discontinuous order parameter, while these signatures are absent upon cooling. Complementary phase-field simulations demonstrate similar behavior, attributing it to distinct elastic strain energy accumulation and release during heating and cooling. These findings reveal a first-order character upon heating but a second-order-like behavior upon cooling, challenging the conventional paradigm of symmetric phase-transition classification and suggesting new possibilities for ferroelectric-based energy storage.

Hershkovitz, Asaf [Technion-Israel Institute of Te↗

DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Spinodal electronic phase separation during insulator-metal transitions

Electronic phase transitions such as insulator-metal transitions are common in strongly correlated systems. Here, using a combination of thermodynamic linear-stability analysis and phase-field simulations and employing VO 2 as a prototypical example, we predict that an insulator-metal transition driven by photoexcitation may involve an intermediate, modulated charge density state with a temperature-dependent characteristic wavelength. It is shown that such an intermediate two-phase electronic state is formed through a spinodal mechanism and that its formation can be generic for insulator-metal transitions driven by fast stimuli. This transient electronic phase separation is expected to stimulate future experimental and computational efforts.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Phase field modeling of shearing processes of a dual-lobed γ"|γ'|γ" coprecipitate

The Ni-based superalloy IN718 is one of the most widely used commercial alloy in the aerospace industry since its development in 1960s. The excellent mechanical properties have been attributed in a great deal to strengthening by coherent γ' and γ" precipitates. The deformation mechanisms of these two phases have been well studied individually. Recent experimental characterization has shown coprecipitates of these two phases in a variety of morphologies and it was argued that these coprecipitates may lead to improved strengthening as compared to their monolithic counterparts. However, the deformation mechanisms of these coprecipitates are still not well understood. In this study, we performed microscopic phase field simulations, with generalized-stacking-fault (GSF) energy surfaces from ab initio calculations as inputs, to systematically study the shearing processes of a periodical array of dual-lobed coprecipitates as well as monolithic precipitates. We found that the coupling between the γ' and γ" phases in the coprecipitates forces dislocations to take high energy shearing pathways in both phases that would not occur if they were in monolithic forms. The coupling also creates stacking fault configurations in the coprecipitates that require high stress to form. Thus, the presence of coprecipitates in general should increase the resistance to dislocation shearing and lead to higher strength levels. Various fault configurations observed during the shearing process are documented as a reference for future comparison with experimental observations. The link between stacking fault shearing and microtwinning is also discussed. Here, the mechanisms analyzed in this study deepens our understanding of coprecipitation effects on alloy strength and may form a cornerstone for multi-precipitate strengthened alloy design strategies.

36 MATERIALS SCIENCE↗