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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A guide to the design of the virtual element methods for second- and fourth-order partial differential equations

Here we discuss the design and implementation details of two conforming virtual element methods for the numerical approximation of two partial differential equations that emerge in phase-field modeling of fracture propagation in elastic material. The two partial differential equations are: (i) a linear hyperbolic equation describing the momentum balance and (ii) a fourth-order elliptic equation modeling the damage of the material. Inspired by, we develop a new conforming VEM for the discretization of the two equations, which is implementation-friendly, i.e., different terms can be implemented by exploiting a single projection operator. We use C 0 and C 1 virtual elements for the second-and fourth-order partial differential equation, respectively. For both equations, we review the formulation of the virtual element approximation and discuss the details pertaining the implementation.

42 ENGINEERING↗
DOE OSTI2020

Phase field simulation of anode microstructure evolution of solid oxide fuel cell through Ni(OH) 2 diffusion

Microstructure evolution in a solid oxide fuel cell anode is strongly affected by operating condition. The detailed mesoscale evolution mechanism under operating condition is still under debate. Here we develop a phase-field model to simulate microstructure evolution through formation and diffusion of gaseous Ni(OH)2. We studied the coarsening kinetics and redistribution of Ni under different steam distributions and compare it to that under pure hydrogen condition. Overall, the results suggest that although the presence of a steady gradient of steam leads to redistribution of Ni, Ni(OH) 2 formation and diffusion do not significantly change the Ni coarsening rate under strictly humid conditions, contrary to commonly reported hypotheses. It is concluded that competing mechanisms other than Ni(OH) 2 diffusion must be responsible for the experimentally observed Ni redistribution and enhanced coarsening under humid conditions.

25 ENERGY STORAGE↗
DOE OSTI2017

Integrated Predictive Tools for Customizing Microstructure and Material Properties of Additively Manufactured Aerospace Components

Additive Manufacturing (AM) refers to a process by which digital three-dimensional (3-D) design data is converted to build up a component by depositing material layer-by-layer. United Technologies Corporation (UTC) is currently involved in fabrication and certification of several AM aerospace structural components made from aerospace materials. This is accomplished by using optimized process parameters determined through numerous design-of-experiments (DOE)-based studies. Certification of these components is broadly recognized as a significant challenge, with long lead times, very expensive new product development cycles and very high energy consumption. Because of these challenges, United Technologies Research Center (UTRC), together with UTC business units have been developing and validating an advanced physics-based process model. The specific goal is to develop a physics-based framework of an AM process and reliably predict fatigue properties of built-up structures as based on detailed solidification microstructures. Microstructures are predicted using process control parameters including energy source power, scan velocity, deposition pattern, and powder properties. The multi-scale multi-physics model requires solution and coupling of governing physics that will allow prediction of the thermal field and enable solution at the microstructural scale. The state-of-the-art approach to solve these problems requires a huge computational framework and this kind of resource is only available within academia and national laboratories. The project utilized the parallel phase-fields codes at Oak Ridge National Laboratory (ORNL) and Lawrence Livermore National Laboratory (LLNL), along with the high-performance computing (HPC) capabilities existing at the two labs to demonstrate the simulation of multiple dendrite growth in threedimensions (3-D). The LLNL code AMPE was used to implement the UTRC phase field model that was previously developed for a model binary alloy, and the simulation results were compared against the UTRC simulation results, followed by extension of the UTRC model to simulate multiple dendrite growth in 3-D. The ORNL MEUMAPPS code was used to simulate dendritic growth in a model ternary alloy with the same equilibrium solidification range as the Ni-base alloy 718 using realistic model parameters, including thermodynamic integration with a Calphad based model for the ternary alloy. Implementation of the UTRC model in AMPE met with several numerical and parametric issues that were resolved and good comparison between the simulation results obtained by the two codes was demonstrated for two dimensional (2-D) dendrites. 3-D dendrite growth was then demonstrated with the AMPE code using nondimensional parameters obtained in 2-D simulations. Multiple dendrite growth in 2-D and 3-D were demonstrated using ORNL’s MEUMAPPS code using simple thermal boundary conditions. MEUMAPPS was then modified to incorporate the complex, time-dependent thermal boundary conditions obtained by UTRC’s thermal modeling of single track AM experiments to drive the phase field simulations. The results were in good agreement with UTRC’s experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A combined kinetic Monte Carlo and phase field approach to model thermally activated dislocation motion

In this work, a kinetic Monte Carlo (KMC) scheme is integrated into a phase field dislocation dynamics (PFDD) model to account for thermal activation. While the implementation is general and could be applied to many different thermally activated processes, we apply the model to simulate a dislocation transmitting through an interface. Two cases are analyzed: a dislocation transmitting through a strong interface, and a dislocation transmitting through a more realistic grain boundary in tungsten. In conclusion, results and trends for both cases are discussed in detail, and ultimately show that the combined KMC-PFDD algorithm is effective in capturing thermally activated dislocation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Kinetics of Domain Switching by Mechanical and Electrical Stimulation in Relaxor-Based Ferroelectrics

Ferroelectric materials have been extensively explored for applications in high-density nonvolatile memory devices because of their ferroelectric-ferroelastic domain-switching behavior under electric loading or mechanical stress. However, the existence of ferroelectric and ferroelastic backswitching would cause significant data loss, which affects the reliability of data storage. In this work, we apply in situ transmission electron microscopy and phase-field modeling to explore the unique ferroelastic domain-switching kinetics and the origin of this in relaxor-based Pb(Mg 1/3 Nb 2/3 )O 3 -33%PbTiO 3 single-crystal pillars under electrical and mechanical stimulations. Results showed that the electric-mechanical hysteresis loop shifted for relaxor-based single-crystal pillars because of the low energy levels of domains in the material and the constraint on the pillars, resulting in various mechanically reversible and irreversible domain-switching states. The phenomenon can potentially be used for advanced bit writing and reading in nonvolatile memories, which effectively overcomes the backswitching problem and broadens the types of ferroelectric materials for nonvolatile memory applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Thermodynamics of Light-Induced Nanoscale Polar Structures in Ferroelectric Superlattices

Here we study the thermodynamics of nanoscale polar structures in PbTiO 3 /SrTiO 3 ferroelectric superlattices induced by above-bandgap optical excitation using a phase-field model explicitly considering both structural and electronic processes. We demonstrate that the light-excited carriers provide the charge compensation of polarization bound charges and the lattice thermal energy, both of which are key to the thermodynamic stabilization of a previously observed supercrystal, a three dimensionally periodic nanostructure, within a window of substrate strains, while different mechanical and electrical boundary conditions can stabilize a number of other nanoscale polar structures by balancing the competing short-range exchange interactions responsible for the domain wall energy and longrange electrostatic and elastic interactions. The insights into the light-induced formation and richness of nanoscale structures from this work offer theoretical guidance for exploring and manipulating the thermodynamic stability of nanoscale polar structures employing a combination of thermal, mechanical, and electrical stimuli as well as light.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Predicting the morphologies of γ' precipitates in cobalt-based superalloys

Cobalt-based alloys with γ/γ' microstructures have the potential to become the next generation of superalloys, but alloy compositions and processing steps must be optimized to improve coarsening, creep, and rafting behavior. While these behaviors are different than in nickel-based superalloys, alloy development can be accelerated by understanding the thermodynamic factors influencing microstructure evolution. In this work, we develop a phase field model informed by first-principles density functional theory and experimental data to predict the equilibrium shapes of Co-Al-W γ' precipitates. Three-dimensional simulations of single and multiple precipitates are performed to understand the effect of elastic and interfacial energy on coarsened and rafted microstructures; the elastic energy is dependent on the elastic stiffnesses, misfit strain, precipitate size, applied stress, and precipitate spatial distribution. We observe characteristic microstructures dependent on the type of applied stress that have the same γ' morphology and orientation seen in experiments, indicating that the elastic stresses arising from coherent γ/γ' interfaces are important for morphological evolution during creep. Here, the results also indicate that the narrow γ channels between γ' precipitates are energetically favored, and provide an explanation for the experimentally observed directional coarsening that occurs without any applied stress.

36 MATERIALS SCIENCE↗