Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase Field”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theory and phase-field simulations of electrical control of spin cycloids in a multiferroic

Cycloidal spin orders are common in multiferroics. One of the prototypical examples is BiFeO 3 (BFO) which shows a large polarization and a cycloidal antiferromagnetic (AFM) order at room temperature. Here, we employ Landau theory and phase-field simulations to analyze the coupled switching dynamics of polarization and cycloidal AFM orders in BFO. We identify 14 types of transitional spin structures between two cycloids and nine electric-field-induced spin switching paths. We demonstrate the electric-field-induced rotation of wave vectors of the cycloidal spins and discover two types of cycloidal spin switching dynamics: fast local spin flips and slow rotation of wave vectors. Also, we construct roadmaps to achieve the switching between any two spin cycloids through multistep applications of electric fields. Furthermore, this paper provides a theoretical framework for the phenomenological description of spin cycloids and a fundamental understanding of the switching mechanisms to achieve electrical control of magnetic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2025

Impact of position and density of nanoscale voids on fracture initiation in iron from phase field fracture simulation

Understanding the impact of these bubbles on crack propagation, like that of helium bubble-induced cracking in irradiated materials is incredibly complex. A useful first study towards understanding bubble effects on fracture is to examine how voids impact fracture first. In this work, we used phase-field fracture simulations to examine the influence of voids and their distribution on Mode I fracture in Fe. Assuming brittle fracture, two simulation configurations were considered: (1) nanoscale systems with one or two voids, and (2) nanoscale systems with an experimentally relevant distribution of voids, with up to 20% void area. Results from simulations with one and two voids showed that voids within 10 nm of a crack tip reduce the stress required for crack growth, with the magnitude of reduction depending on void-to-crack orientation. Comparisons with linear elastic fracture mechanics and evaluation of one versus two void systems revealed deviations from linear superposition, implying complex interactions between void and crack tip stress fields. In multi-void simulations, as void sizes increase, the nearest void to the crack tip exerts a greater influence on fracture stress than the overall porosity. Furthermore this study provides valuable insights into the relationship between void size and concentration, and the stress necessary for crack growth, marking a step forward towards understanding He bubble-induced fracture in ferrous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of a grain growth model for U 3 Si 2 using experimental data, phase field simulation and molecular dynamics

The purpose of this work is to develop a model for normal grain growth in U 3 Si 2 . The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. Here, the mobility was determined by comparing phase field grain growth simulations to the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boundary conditions manipulation of polar vortex domains in BiFeO 3 membranes via phase-field simulations

Polar vortex domains have recently become an emergent research field due to the abundant physical phenomena and potential applications in high-density memories. Here, we explore the mechanisms of creating polar vortex domains in the BiFeO3 (BFO) membranes subjected to different boundary conditions using phase-field simulations. A major difference is that the vortex in membrane can be stabilized even under short-circuit electrical boundary conditions compared to vortex in other systems, such as thin film or superlattice. We found that (a) the formation of polar vortex domains at the membrane interior under bending is mainly driven by the reduction of elastic energy under short-circuit boundary condition, and the vortex chirality (namely, clockwise and counterclockwise) could be identified by n-shape and u-shape bending; (b) in the unbent open-circuit BFO membrane case, exotic trapezoid-shaped vortex nanodomains form at the terminations of 109 degrees domain walls (DWs) and partially charged 71 degrees DWs, which is driven by the local depolarization field and the interplay among electrostatic, elastic, and gradient and Landau energies. We also examine Kittel's law by establishing the dependence of vortex periods on the membrane thickness. Furthermore, these results give further understanding of the effect of boundary conditions on the formation of polar vortex domains, guiding experimental designs of vortex-based high-density memories.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Domain patterns and super-elasticity of freestanding BiFeO 3 membranes via phase-field simulations

Super-elasticity of functional ferroelectric oxides offers promises for integrating ferroelectric films into flexible electronics. However, super-elastic deformation is a complex phenomenon related to possibly multiple concurrent mechanisms. Fundamentally understanding how multiple mechanisms contribute to the super-elasticity of ferroelectric oxides is crucial to realizing their potential flexible electronic applications. In this study, we employ phase-field simulations to model the dynamics of ferroelectric domain patterns of freestanding BiFeO 3 membranes to understand the origin of their super-elasticity under substantial bending deformation (5% strain). It is demonstrated that both a reversible Rhombohedral-Tetragonal (R-T) phase transition and a nearly reversible domain evolution of BiFeO 3 membranes contribute to accommodating the large deformation and thus their super-elasticity. The dynamics of domain evolution also reveal the formation of an exotic ferroelectric vortex and polarization rotation before the phase transition. We constructed a diagram of phases and domain patterns as a function of the membrane thickness and bending angle, which allows one to readily predict the emergence of T phase and ferroelectric vortex in bent BFO membranes. These results not only provide fundamental understanding of mesoscale super-elastic mechanisms but also reveal exotic domain states of ferroelectric membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of four strain energy decomposition methods for phase field fracture models using quasi-static and dynamic benchmark cases

Abstract Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.

Zhang, Shuaifang↗
DOE OSTI2024

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2022

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Phase-field simulation of domain size effect on dielectric and piezoelectric responses in K 0.5 Na 0.5 NbO 3 epitaxial thin films with superdomain structures

Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field simulations to inform nuclear fuel performance modeling

Software tools to simulate nuclear fuel performance at the engineering scale, such as Idaho National Laboratory (INL)’s BISON code, are increasingly relied upon in regulatory and economic decision-making. However, accurate results from these tools depends on the availability of materials parameters that are used as input. In recent years, atomistic and mesoscale simulation methods have emerged as a cost-effective, expedient means to obtain such input parameters. Phase-field simulations using INL’s Marmot application have been used to obtain microstructure-level parameters and to improve material models for fuel performance modeling using BISON. In this talk, recent examples of this process are given, including applications in UO2, U3Si2, and UZr fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Phase field modeling of hydrogen release in nitrile-butadiene rubber composites after high-pressure hydrogen exposure

Experimental data showed that a simple diffusion model doesn’t suffice to explain the hydrogen (H 2 ) release kinetics in nitrile butadiene rubber (NBR) composites (one with no filler or plasticizer, one with plasticizer only, one with fillers and plasticizer, and one with fillers only) after exposed to hydrogen gas at a pressure of 27.6 MPa. In this work, we developed a phase field model that considers the effect of H 2 gas bubble evolution on H 2 release kinetics. With the model, we simulated the effect of thermodynamic and kinetics properties on the nucleation, growth, and shrinkage of gas bubbles and H 2 release kinetics during and after decompression. Here, the results demonstrated that 1) the model built upon certain thermodynamic and kinetics properties of a given material can well describe its H 2 release kinetics measured in experiments, and 2) the model can predict the gas bubble evolution, which is associated with the material property degradation and failures, with accurate thermodynamic and kinetics properties or H 2 release data from experiments.

08 HYDROGEN↗