Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field models”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Phase-field modeling of void anisotropic growth behavior in irradiated zirconium

A three-dimensional (3D) phase field model was developed to study the effects of surface energy and diffusivity anisotropy on void growth behavior in irradiated Zr. The gamma surface energy function, which is used in the phase field model, was developed with the surface energy anisotropy calculated from the molecular dynamics (MD) simulations. It is assumed that vacancies have much larger mobility in c-axis than a- and b- axes while interstitials have much larger mobility in basal plane then that in c-axis. With the model, the equilibrium void morphology and the effect of defect concentrations and defect mobility anisotropy on void growth behavior were simulated. The simulations demonstrated that 1) The developed phase-field model can correctly reproduce the faceted void morphology predicted by the Wullf construction. 2) With isotropic diffusivity the void prefers to grow on the basal plane. 3) When the vacancy has large mobility along c-axis and interstitial has a large mobility on the basal plane of hexagonal closed packed (hcp) Zr alloys a platelet void grows in c-direction and shrinks on the basal plane, which is in agreement with the experimental observation of void growth behavior in irradiated Zr.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2020

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models

In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the summation constraint, the conservation constraint, and the consistency of reduction. The first application of the volume distribution algorithm is to determine the Lagrange multipliers in multiphase Phase-Field models that enforce the mass conservation, and a multiphase conservative Allen-Cahn model that satisfies the consistency of reduction is developed. A corresponding consistent and conservative numerical scheme is developed for the model. The multiphase conservative Allen-Cahn model has a better ability than the multiphase Cahn-Hilliard model to preserve under-resolved structures. The second application is to develop a numerical procedure, called the boundedness mapping, to map the order parameters, obtained numerically from a multiphase model, into their physical interval, and at the same time to preserve the physical properties of the order parameters. Along with the consistent and conservative schemes for the multiphase Phase-Field models, the numerical solutions of the order parameters are reduction consistent, conservative, and bounded, which are theoretically analyzed and numerically validated. Then, the multiphase Phase-Field models are coupled with the momentum equation by satisfying the consistency of mass conservation and the consistency of mass and momentum transport, thanks to the consistent formulation. Finally, it is demonstrated that the proposed model and scheme converge to the sharp-interface solution and are capable of capturing the complicated multiphase dynamics even when there is a large density and/or viscosity ratio.

42 ENGINEERING↗
DOE OSTI2022

A phase-field model for void and gas bubble superlattice formation in irradiated solids

A phase-field model to simulate the formation of both void and gas bubble superlattices is derived from a grand potential functional, assuming 1D diffusion of self-interstitial atoms. The model is capable of accounting for superlattice formation by either a nucleation and growth or spinodal decomposition mechanism; in this work, we focus on the nucleation and growth mechanism, using a discrete nucleation approach in the phase-field model. In simulations of void formation, short aligned rows of voids were initially formed, followed by growth in the size of the aligned rows, and finally leading to superlattice formation, consistent with experimental observations. Void superlattice spacing as a function of model parameters was quantified. Increased nucleation rates led to smaller superlattice spacing, while increased vacancy diffusivity led to larger superlattice spacing. In simulations of gas bubble superlattice formation, superlattice spacing decreased with increasing gas atom flux, consistent with experimental observations. Here, simulations with increasing ratio of gas atom flux to vacancy-interstitial pair production led to increased superlattice spacing, which is inconsistent with experimental observations at high gas atom flux rates; this discrepancy is attributed to the fact that the model does not account for bound gas atom-vacancy pairs that may dominate at high gas flux rates.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Phase-field model of pitting corrosion kinetics in metallic materials

This study presents a thermodynamically consistent phase field model for the quantitative prediction of the pitting corrosion kinetics in metallic materials. An order parameter is introduced to represent the metal’s physical state at each point in the entire domain. The free energy of the metal-electrolyte system is defined in terms of its metal ion concentration and the order parameter. The mass transport in the electrolyte and the electrochemical reactions in the electrolyte/metal interface are taken into consideration in the model to simulate metal corrosion in a corrosive environment. The governing equations for mass transport and the order parameter are solved in such a manner that the system’s free energy is reduced as a result of diffusion and migration processes, and the distribution of electrostatic potential is governed by Poisson’s equation. A calibration study is performed to couple the kinetic interface parameter with the corrosion current density, which results in a direct relationship between overpotential and the kinetic interface parameter. A comparative study is performed to validate the phase field model against the experimental results. Several case studies are also examined to understand the corrosion behavior of closely located pits, stressed material, ceramic particles-reinforced steel, and crystallographic orientation dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2021

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A phase-field model for hydride formation in polycrystalline metals: Application to δ-hydride in zirconium alloys

We report a phase-field model for simulating metal hydride formation involving large volume expansion in single- and polycrystals. As an example, we consider δ-hydride formation in α-zirconium (Zr), which involves both displacive crystallographic structural change and hydrogen diffusion process. Thermodynamic Gibbs energy functions are extracted from the available thermodynamic database based on the sublattice model for the interstitial solid solutions. Solute-grain boundary interactions and inhomogeneous elasticity of polycrystals are taken into consideration within the context of diffuse-interface description. The stress-free transformation strains of multiple variants for hcp -Zr (α) to fcc -hydride (δ) transformation are derived based on the well-established orientation relationship between the α and δ phases as well as the corresponding temperature-dependent lattice parameters. In particular, to account for the large volume expansion, we introduced the mixed interfacial coherency concept between those phases—basal planes are coherent and prismatic planes are semi-coherent in computing the strain energy contribution to the thermodynamics. We analyzed the morphological characteristics of hydrides involving multiple structural variants and their interactions with grain boundaries. Moreover, our simulation study allows for the exploration of the possible hydride re-orientation mechanisms when precipitating under applied tensile load, taking into account the variation in the interfacial coherency between hydrides and matrix, their elastic interactions with the applied stress, as well as their morphology-dependent interactions with grain boundaries. Here, the phase-field model presented here is generally applicable to hydride formation in any binary metal–hydrogen systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effects of introducing elasticity using different interpolation schemes to the grand potential phase field model

Introducing elastic energy in the phase field method has been shown to influence interfacial energy, depending on the elastic interpolation scheme. This study investigates the impact of the elastic energy when using a grand potential-based phase field method, comparing the result of Khachaturyan’s strain interpolation scheme (KHS) and Voight-Taylor’s elastic energy interpolation scheme (VTS). The KHS model leads to a decrease in the interfacial energy, while the VTS model leads to an increase. The change in interfacial energy is greater with the VTS model than the KHS model, which suggests that the KHS model is more appropriate to limit the artificial impact of the elastic energy on the interfacial energy. When the contribution at the interface is not negligible, it is shown that both the microstructure evolution kinetics and the equilibrium microstructure can be influenced by the choice of the elastic scheme being used. In addition, this paper shows that the grand potential model might not be appropriate when the system requires the introduction of a composition-dependent term in the elastic energy contribution. This limitation is due to the need for an explicit and invertible relation between the total potential and the composition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2014

Phase field modelling of stressed grain growth: Analytical study and the effect of microstructural length scale

We establish the correlation between the diffuse interface and sharp interface descriptions for stressed grain boundary migration by presenting analytical solutions for stressed migration of a circular grain boundary in a bicrystalline phase field domain. The validity and accuracy of the phase field model is investigated by comparing the phase field simulation results against analytical solutions. The phase field model can reproduce precise boundary kinetics and stress evolution provided that a thermodynamically consistent theory and proper expressions for model parameters in terms of physical material properties are employed. Quantitative phase field simulations are then employed to investigate the effect of microstructural length scale on microstructure and texture evolution by stressed grain growth in an elastically deformed polycrystalline aggregate. The simulation results reveal a transitional behaviour from normal to abnormal grain growth by increasing the microstructural length scale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Visualizing Uranium Crystallization from Melt: Experiment-Informed Phase Field Modeling and Machine Learning

The focus of this project was to observe and simulate the solidification of uranium metal at the crystallographic level from its molten state. Melting experiments were conducted at two different scales to observe microstructural evolution using either a laboratory-scale induction furnace (hundreds of grams of metal) or a microscope heating stage (hundreds of milligrams of metal), respectively. Experimental parameters and characterization data were then used to inform a phase field model of gamma-U crystal growth as dendrites with or without secondary phase impurities in the form of uranium carbide particles. Finally, training datasets were generated by the phase field model as inputs to a neural network, developed with the aim of providing a faster, cheaper surrogate model for microstructural simulations within a given parameter space. Progress is reported herein for each of these task areas. Ultimately, 1) an optical microscope heating stage capability has been stood-up for uranium metal solidification studies, 2) a phase field model was advanced to simulate multiple uranium grains growing in the presence of carbide impurity particles and 3) a neural network was constructed and optimized to predict the microstructure features of individually growing uranium crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Phase field modeling of recrystallization phenomena

In this work, we studied the recrystallization phenomena in copper using phase field modeling technique. EBSD information is successfully translated on to MOOSE and quantitative phase field simulations for Copper are performed on MOOSE. DREAM 3D software is utilized to transform EBSD data into a format which is then used to conduct phase field simulations, with initial grain structure obtained from EBSD data. Several simulations are performed and a comparison is made against various cases: isotropic grain growth without any deformed grains, deformed grains with constant dislocation density, and deformed grains with orientation dependent dislocation density.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Multicomponent phase-field model for extremely large partition coefficients

We develop a multicomponent phase-field model specially formulated to robustly simulate concentration variations from molar to atomic magnitudes across an interlace, i.e., partition coefficients in excess of 10±23 such as may be the case with species which are predominant in one phase and insoluble in the other. Substitutional interdiffusion on a normal lattice and concurrent interstitial diffusion are included. The composition in the interlace follows the approach of Kim. Kim, and Suzuki [Phys. Rev. E 60, 7186 (1999)] and is compared to that of Wheeler, Boettinger, and McFadden [Phys. Rev. A 45, 7424 (1992)] in the context of large partitioning. The model successfully reproduces analytical solutions for binary diffusion couples and solute trapping for the demonstrated cases of extremely large partitioning.

36 MATERIALS SCIENCE↗