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 OSTI2022

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↗
DOE OSTI2022

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↗
NASA NTRS1998

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.↗
DOE OSTI2016

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↗
DOE OSTI2025

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↗
DOE OSTI2017

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↗
DOE OSTI2014

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↗
DOE OSTI2017

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
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↗