Engineering Papers⌕ Search

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.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mesoscale modeling of restructuring in high burnup UO 2 fuel

Here, this work aims to simulate the restructuring behavior observed in different regions of high burnup fuel, providing a first-of-its-kind restructuring model for the dark zone and rim region of high-burnup UO 2 fuel. We employed a grand-potential-based phase-field model to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of different initial conditions and different modeling parameters were systematically studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to begin around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Restructuring was demonstrated to be influenced by a combination of initial dislocation densities, burnup rate, subgrain formation rate, and temperature. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. Furthermore, bubble size and distribution do not significantly affect the rate of restructuring. The predicted microstructures are consistent with experimental observations of the restructured regions. Finally, we present a correlation demonstrating the evolution of the restructuring volume fraction as a function of local burnup.

UO2↗
DOE OSTI2019

A novel liquid-mediated nucleation mechanism for explosive crystallization in amorphous germanium

We report a novel mechanism for explosive crystallization in amorphous germanium (a-Ge), which operates through liquid-mediated nucleation occurring under extreme thermal gradient conditions. The crystallization kinetics of sputter-deposited films with thicknesses ranging from 30 to 150 nm were characterized using in situ movie-mode dynamic transmission electron microscopy (MM-DTEM). After localized heating from a short laser pulse, explosive liquid phase nucleation (LPN) was observed to occur during the early stage (<2 μs) of crystallization in the thicker (>50 nm) films deposited on silicon nitride substrates. The crystallization front propagated at ~12–15 m/s and produced nanocrystalline microstructure with ~50 nm grains. A mechanism involving the existence of a relatively thick (>100 nm) transient liquid layer and a high nucleation rate is proposed to explain the behavior. In conclusion, the key thermodynamic and kinetic features as well as the feasibility of the mechanism are further explored by employing parametric and systematic phase-field modeling and simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructure-Based Degradation Modeling in Solid Oxide Cells

The presentation is intended for a graduate class in Georgia Southern University taught by Dr. Hayri Sezer on electrochemistry, SOFC and Li ion battery. It summaries NETL's recent works on microstructure-based modeling of degradation in solid oxide cells. It covers the degradation mechanism of Ni coarsening, Ni migration and the build-up of oxygen partial pressure inside the electrolyte. The phase-field models developed for Ni coarsening and migration are reviewed and the results are compared to experiments in literature to examine different diffusion mechanisms (Ni self-diffusion, Ni(OH)2 diffusion in pores and Ni(OH) surface diffusion) and driving forces (Ni(OH)x concentration and Ni-YSZ contact angle change). It also demonstrates the developed model for oxygen partial pressure and shows that the predicted oxygen partial pressure can explain the Ni oxidation in the hydrogen electrode and crack growth on the electrolyte-oxygen electrode interface under certain operating conditions.

nickel coarsening↗
DOE OSTI2012

Effects of internal stresses and intermediate phases on the coarsening of coherent precipitates: A phase-field study

Phase stability, topology and size evolution of precipitates are important factors in determining the mechanical properties of crystalline materials. In this article, the Cahn-Hilliard type of phase-field model was coupled to elasticity equations within a mixed-order Galerkin finite element framework to study the coarsening morphology of coherent precipitates. The effects of capillarity, particle size and fraction, compositional strain, and inhomogeneous elasticity on the kinetics and kinematics of coherent precipitates in a binary dual phase crystal admitting a third intermediate stable/meta-stable phase were investigated. The results demonstrated the ability of the model to simulate coarsening under the concomitant action of Ostwald ripening and mismatch elastic strain mechanisms. Using a phenomenological coarsening power law, coarsening rates were determined to depend on precipitate size and volume fraction, compositional strain, and strain mismatch between precipitates and the matrix. Results also showed that the necking incubation time between two neighboring precipitates depends inversely on the precipitate's initial sizes; however, under fixed volume fraction of precipitates, any increase in the initial sizes of the precipitates mitigates the coarsening. Meanwhile, the compositional strain and the growth of the intermediate stable/meta-stable phase leads to substantial enhancements of precipitate coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A 3D phase field dislocation dynamics model for body-centered cubic crystals

In this work, we present a 3D Phase Field Dislocation Dynamics (PFDD) model for body-centered cubic (BCC) metals. The model formulation is extended to account for the dependence of the Peierls barrier on the line-character of the dislocation. Simulations of the expansion of a dislocation loop belonging to {110} < 111 > the slip system are presented with direct comparison to Molecular Statics (MS) simulations. The extended PFDD model is able to capture the salient features of dislocation loop expansion predicted by MS simulations. In conclusion, the model is also applied to simulate the motion of a straight screw dislocation through kink-pair motion.

36 MATERIALS SCIENCE↗
DOE OSTI2017

FY17 end of year MARMOT materials model development update

This document summarizes the lower length scale progress made in FY17 in developing materials models for fuel performance modeling of oxide fuels and their implementation in the MARMOT code. The progress covers the development of a phase field recrystallization model for high-burnup structure formation in oxide fuels, a phase field model for studying stage-3 fission gas release in UO 2 , a phase-field model for irradiation-enhanced densification in UO 2 , a quantitative model for hydride morphology in Zirconium based cladding, and the calculation of hydrogen diffusivity under stress in various Zirconium based alloys as cladding candidates. Further development in the MARMOT code has been carried out as well to improve the multi-phase, multi-component phase field model, the representation of interfacial energy in the Kim-Kim-Suzuki model, and the description of interfacial stress. For most of the above, separated reports have been generated with detailed description. Here a brief summary is provided, with references pointing to the detailed reports to guide interested readers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Integrated Simulation Framework for Additively Manufactured Ti-6Al-4V: Melt Pool Dynamics, Microstructure, Solid-State Phase Transformation, and Microelastic Response

In this work, to accelerate the establishment of fundamental understanding of the additive manufacturing (AM) process and its influence on microstructural evolution and related properties, we develop a multiphysics and multiscale modeling framework that integrates: (1) a high-fidelity powder-scale three-dimensional simulation of transient heat transfer and melt flow dynamics, (2) cellular automaton simulation of solidification grain structure and texture, (3) phase-field modeling of precipitation and dissolution of second-phase precipitate during repeated thermal cycles, and (4) microstructure-based micro- and mesoscopic elastic response calculation. Using Ti-6Al-4V as a model system, we demonstrate the application of the integrated framework to simulate complex microstructure evolution during a single-track laser powder bed fusion process and the associated mechanical response. Our modeling framework successfully captures the solidification β grain structure as a function of laser power and scanning speed, α precipitation upon subsequent cooling with different rates, and elastic response of the resulting (α + β) two-phase microstructure. The key features of solidification and second-phase precipitate microstructures, and their dependence on processing parameters, agree well with existing experimental observations. Finally, the established modeling framework is generally applicable to other metallic materials fabricated by AM.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Colossal Strain Tuning of Ferroelectric Transitions in KNbO 3 Thin Films

Strong coupling between polarization ( P ) and strain (ɛ) in ferroelectric complex oxides offers unique opportunities to dramatically tune their properties. Here colossal strain tuning of ferroelectricity in epitaxial KNbO 3 thin films grown by sub‐oxide molecular beam epitaxy is demonstrated. While bulk KNbO 3 exhibits three ferroelectric transitions and a Curie temperature ( T c ) of ≈676 K, phase‐field modeling predicts that a biaxial strain of as little as −0.6% pushes its T c > 975 K, its decomposition temperature in air, and for −1.4% strain, to T c > 1325 K, its melting point. Furthermore, a strain of −1.5% can stabilize a single phase throughout the entire temperature range of its stability. A combination of temperature‐dependent second harmonic generation measurements, synchrotron‐based X‐ray reciprocal space mapping, ferroelectric measurements, and transmission electron microscopy reveal a single tetragonal phase from 10 K to 975 K, an enhancement of ≈46% in the tetragonal phase remanent polarization ( P r ), and a ≈200% enhancement in its optical second harmonic generation coefficients over bulk values. These properties in a lead‐free system, but with properties comparable or superior to lead‐based systems, make it an attractive candidate for applications ranging from high‐temperature ferroelectric memory to cryogenic temperature quantum computing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A phase field study of the thermal migration of gas bubbles in UO 2 nuclear fuel under temperature gradient

Phase field models are developed to study the gas bubble migration in uranium dioxide nuclear fuel in which a large temperature gradient exists during the operation. In this work, thermal diffusion mechanism for nanosized gas bubbles and vapor transport process for micron-sized gas bubbles are considered, respectively. In both cases, gas bubbles migrate to the high-temperature area. Due to the velocity difference between leading and trailing edges of the gas bubbles, nanosized gas bubbles are elongated along the temperature gradient direction when thermal diffusion is dominated. Micron-sized gas bubbles are either compressed along temperature gradient direction to form lenticular shape bubbles or elongated along temperature gradient direction, depending on the location of the gas bubbles within the fuel pellet. Initial gas bubble radius has no significant effect on the gas bubble migration velocity for both thermal diffusion and vapor transport mechanisms. We notice that the shape change of the gas bubble due to vapor transport mechanism has no significant effect on the migration velocity. Furthermore, the center cavity formation is also captured by our model which is due to the migration and accumulation of lenticular gas bubbles at the center of the fuel pellet. The modeling results compare well with experimental observations and theoretical analysis in the literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2011

Overpotential-Dependent Phase Transformation Pathways

An objective in battery development for higher storage energy density is the design of compounds that can accommodate maximum changes in ion concentration over useful electrochemical windows. Not surprisingly, many storage compounds undergo phase transitions in situ, including production of metastable phases. Unique to this environment is the frequent application of electrical over- and underpotentials, which are the electrical analogs to undercooling and superheating. Surprisingly, overpotential effects on phase stability and transformation mechanisms have not been studied in detail. Here we use synchrotron X-ray diffraction performed in situ during potentiostatic and galvanostatic cycling, combined with phase-field modeling, to reveal a remarkable dependence of phase transition pathway on overpotential in the model olivine Li{sub 1-x}FePO{sub 4}. For a sample of particle size {approx}113 nm, at both low (e.g., <20 mV) and high (>75 mV) overpotentials a crystal-to-crystal olivine transformation dominates, whereas at intermediate overpotentials a crystalline-to-amorphous phase transition is preferred. As particle sizes decrease to the nanoscale, amorphization is further emphasized. Implications for battery use and design are considered.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Compositionally-Driven Formation Mechanism of Hierarchical Morphologies in Co-Deposited Immiscible Alloy Thin Films

Co-deposited, immiscible alloy systems form hierarchical microstructures under specific deposition conditions that accentuate the difference in constituent element mobility. The mechanism leading to the formation of these unique hierarchical morphologies during the deposition process is difficult to identify, since the characterization of these microstructures is typically carried out post-deposition. We employ phase-field modeling to study the evolution of microstructures during deposition combined with microscopy characterization of experimentally deposited thin films to reveal the origin of the formation mechanism of hierarchical morphologies in co-deposited, immiscible alloy thin films. Our results trace this back to the significant influence of a local compositional driving force that occurs near the surface of the growing thin film. We show that local variations in the concentration of the vapor phase near the surface, resulting in nuclei (i.e., a cluster of atoms) on the film’s surface with an inhomogeneous composition, can trigger the simultaneous evolution of multiple concentration modulations across multiple length scales, leading to hierarchical morphologies. We show that locally, the concentration must be above a certain threshold value in order to generate distinct hierarchical morphologies in a single domain.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Synergistic Effect of Graphene Oxide for Impeding the Dendritic Plating of Li

Dendritic growth of lithium (Li) has severely impeded the practical application of Li-metal batteries. Herein, a 3D conformal graphene oxide nanosheet (GOn) coating, confined into the woven structure of a glass fiber separator, is reported, which permits facile transport of Li-ions thought its structure, meanwhile regulating the Li deposition. Electrochemical measurements illustrate a remarkably enhanced cycle life and stability of the Li-metal anode, which is explained by various microscopy and modeling results. Utilizing scanning electron microscopy, focused ion beam, and optical imaging, the formation of an uniform Li film on the electrode surface in the case of GO-modified samples is revealed. Ab initio molecular dynamics (AIMD) simulations suggest that Li-ions initially get adsorbed to the lithiophilic GOn and then diffuse through defect sites. This delayed Li transfer eliminates the “tip effect” leading to a more homogeneous Li nucleation. Meanwhile, C-C bonds rupture observed in the GO during AIMD simulations creates more pathways for faster Li-ions transport. In addition, phase-field modeling demonstrates that mechanically rigid GOn coating with proper defect size (smaller than 25 nm) can physically block the anisotropic growth of Li. This new understanding is a significant step toward the employment of 2D materials for regulating the Li deposition

25 ENERGY STORAGE↗
DOE OSTI2015

Miscibility Gap Closure, Interface Morphology, and Phase Microstructure of 3D Li x FePO 4 Nanoparticles from Surface Wetting and Coherency Strain

We study the mesoscopic effects which suppress phase-segregation in Li x FePO 4 nanoparticles using a multiphysics phase-field model implement on a high performance cluster. We simulate 3D spherical particles of radii from 3nm to 40nm and examine the equilibrium microstructure and voltage profiles as a they depend on size and overall lithiation. The model includes anisotropic, concentration-dependent elastic moduli, misfit strain, and facet dependent surface wetting within a Cahn-Hilliard formulation. Here, we find that the miscibility gap vanishes for particles of radius ~ 5 nm, and the solubility limits change with overall particle lithiation. The corresponding voltage plateau, indicative of phase-segregation, changes in extent and also vanishes. Surface wetting is found to have a strong effect on stabilizing a variety of microstructures, exaggerating the shifting of solubility limits, and shortening the voltage plateau.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Phase-Field Study on Internal to External Oxidation Transition in High-Temperature Structural Alloys

Structural alloys applied at high temperatures rely on an external dense layer of oxide scale for protection. As some reactive alloy components are selectively oxidized internally, understanding how these dispersed metal oxide particles reach the surface to form a protective scale—the so-called internal to external oxidation transition—is crucial for designing these alloys. While the literature is replete with experimental studies on oxidation of alloys, there is a lack of computational studies in this realm due to the complex nature of coupled reaction and diffusion processes in multicomponent multi-phase alloy systems. In this work, we apply a recently developed phase-field model to simulate the oxidation processes under different compositions and nucleation scenarios to gain insights into how a continuous oxide scale can be established. The results show that while alloy composition is critical for internal to external oxidation transition, the oxide nuclei size, shape and distribution also have significant impact on the transition kinetics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

On the unsteady Darcy–Forchheimer–Brinkman equation in local and nonlocal tumor growth models

A mathematical analysis of local and nonlocal phase-field models of tumor growth is presented that includes time-dependent Darcy–Forchheimer–Brinkman models of convective velocity fields and models of long-range cell interactions. A complete existence analysis is provided. In addition, a parameter-sensitivity analysis is described that quantifies the sensitivity of key quantities of interest to changes in parameter values. Two sensitivity analyses are examined; one employing statistical variances of model outputs and another employing the notion of active subspaces based on existing observational data. Remarkably, the two approaches yield very similar conclusions on sensitivity for certain quantities of interest. The work concludes with the presentation of numerical approximations of solutions of the governing equations and results of numerical experiments on tumor growth produced using finite element discretizations of the full tumor model for representative cases.

Mathematics↗
DOE OSTI2013

SISGR -- Domain Microstructures and Mechanisms for Large, Reversible and Anhysteretic Strain Behaviors in Phase Transforming Ferroelectric Materials

This four-year project (including one-year no-cost extension) aimed to advance fundamental understanding of field-induced strain behaviors of phase transforming ferroelectrics. We performed meso-scale phase field modeling and computer simulation to study domain evolutions, mechanisms and engineering techniques, and developed computational techniques for nanodomain diffraction analysis; to further support above originally planned tasks, we also carried out preliminary first-principles density functional theory calculations of point defects and domain walls to complement meso-scale computations as well as performed in-situ high-energy synchrotron X-ray single crystal diffraction experiments to guide theoretical development (both without extra cost to the project thanks to XSEDE supercomputers and DOE user facility Advanced Photon Source).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗