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

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗
DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2020

Phase-field model of deformation twin-grain boundary interactions in hexagonal systems

Here, a phase-field model for describing dynamic interactions between deformation twins and grain boundaries in a hexagonal close packed (HCP) metal is established. It is applied to simulating the coupled evolution mechanisms of grains and twins in magnesium (Mg) by probing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries. We analyze the effect of the strain relaxation near a grain boundary, which is related to the geometric compatibility arising from the misorientation angle between adjoining grains. Phase-field simulations demonstrate that twin transmission across a grain boundary into a neighboring grain leads to grain boundary migration towards the neighboring grain with a reduced grain boundary width. The preferred nucleation site of new twin variant in the neighboring grain is related to the elastic interaction energy distribution. These predicted transmission behaviors agree well with existing experimental observations and molecular dynamics simulations. By analyzing set of systematic phase-field simulation results, we establish a rotation-angle-related twin variant selection rule for analyzing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries in Mg and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effective parameterization of phase-field models of fission gas bubble growth

Fission gas bubbles are one of the most important microstructural features of ceramic nuclear fuels. As gas bubbles grow and interconnect, they allow release of gases, with important consequences for fuel performance. Phase-field modeling has been increasingly used to simulate the evolution of fission gas bubble microstructural because of its capability to capture complex microstructural features. However, computational performance limitations have made it difficult to simulate all the defects present in fuels during operation. For this reason, phase-field models have often simulated only vacancies and used multiple approaches to include the effect of vacancy-interstitial recombination and sinks in a simplified way. Here, we compare some of the most prevalent approaches, including source-only and source/sink. The kinetics of bubble growth using these approaches are analyzed analytically, and simulations with these approaches are compared to a full vacancy-interstitial model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Dynamical phase-field model of coupled electronic and structural processes

Abstract Many functional and quantum materials derive their functionality from the responses of both their electronic and lattice subsystems to thermal, electric, and mechanical stimuli or light. Here we propose a dynamical phase-field model for predicting and modeling the dynamics of simultaneous electronic and structural processes and the accompanying mesoscale pattern evolution under static or ultrafast external stimuli. As an illustrative example of application, we study the transient dynamic response of ferroelectric domain walls excited by an ultrafast above-bandgap light pulse. We discover a two-stage relaxational electronic carrier evolution and a structural evolution containing multiple oscillational and relaxational components across picosecond to nanosecond timescales. The phase-field model offers a general theoretical framework which can be applied to a wide range of functional and quantum materials with interactive electronic and lattice orders and phase transitions to understand, predict, and manipulate their ultrafast dynamics and rich mesoscale evolution dynamics of domains, domain walls, and charges.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase field modeling of discontinuous dynamic recrystallization in hot deformation of magnesium alloys

We have developed a quantitative phase field model of discontinuous dynamic recrystallization (PF-DDRX) for investigating grain structure evolution and its impact on mechanical response in magnesium alloy during hot deformation. A set of isothermal compression tests were conducted by Gleeble thermo-mechanical simulator and used to determine PF-DDRX model parameters such as critical stresses of recrystallization, activation energies of dislocation, grain boundaries and recrystallization grain nucleation. In the simulations, the experimentally measured average grain size and distribution were used to generate initial microstructures. The results show that the typical ‘necklace’ microstructures observed in DDRX forms at the earlier stage of recrystallization. The inhomogeneity of microstructures first increases rapidly and then decreases gradually and reaches a constant with the increase of dynamic recrystallization volume fraction. Both predicted grain structure and mechanical response are in good agreement with experimental results. This demonstrates that the model parameters determined by experimental data are reliable for the developed PF-DDRX model to predict the grain structure evolution and flow behavior of AZ80 magnesium alloy under hot deformation. This work presents a path to combine experiments and modeling to give quantitative prediction and minimize experimental cost.

magnesium alloy, Phase field method, dynamic recry↗
DOE OSTI2025

Phase-Field Modeling of Damage Evolution in Ceramic Matrix Composite (CMC) and Environmental Barrier Coating (EBC)

Ceramic matrix composites (CMCs) protected by environmental barrier coatings (EBCs) present a promising materials solution for next generation gas turbines. Developments of more robust and efficient EBCs and mechanically tougher CMCs are thus of significant technological importance. Here we develop a phase-field modeling framework that incorporates the thermally grown oxide (TGO), recognized as a critical factor for degradation and failure of EBCs. We simulate crack growth in the TGO and the potential extension into the bond coat / CMC substrate. The model efficiently takes account of the large inelastic deformation induced by the severe volume expansion of TGO, thanks to our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm. A phase-field model is built for damage evolution in CMCs including crack growth and interfacial sliding. The effects of fiber layout and interfacial sliding on the macroscopic toughness of CMCs are revealed by large-scale simulations and compared to experiments.

advanced energy systems and materials↗
NASA NTRS1999

Thin Interface Asymptotics for an Energy/Entropy Approach to Phase-Field Models with Unequal Conductivities

Karma and Rapped recently developed a new sharp interface asymptotic analysis of the phase-field equations that is especially appropriate for modeling dendritic growth at low undercoolings. Their approach relieves a stringent restriction on the interface thickness that applies in the conventional asymptotic analysis, and has the added advantage that interfacial kinetic effects can also be eliminated. However, their analysis focussed on the case of equal thermal conductivities in the solid and liquid phases; when applied to a standard phase-field model with unequal conductivities, anomalous terms arise in the limiting forms of the boundary conditions for the interfacial temperature that are not present in conventional sharp-interface solidification models, as discussed further by Almgren. In this paper we apply their asymptotic methodology to a generalized phase-field model which is derived using a thermodynamically consistent approach that is based on independent entropy and internal energy gradient functionals that include double wells in both the entropy and internal energy densities. The additional degrees of freedom associated with the generalized phased-field equations can be chosen to eliminate the anomalous terms that arise for unequal conductivities.

McFadden, G. B.↗
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2025

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗
DOE OSTI2025

Understanding grain refinement and intergranular gas bubble evolution in U-10Mo fuel using phase-field modeling

Monolithic U-10Mo fuel undergoes significant microstructural changes in the form of grain refinement and gas bubble formation during burnup, which degrades its mechanical properties. In this talk, I present a phase-field model for microstructure evolution in U-10Mo developed using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Simulations demonstrate that grain refinement initiates at pre-existing grain boundaries (GBs) to eliminate the lattice distortion energy caused by the accumulation of self-interstitial loops. By employing an equation of state for xenon gas, we simulate the evolution of gas bubbles in the polycrystal microstructure. Large, interconnected bubbles are found to form along the triple junctions. The effects of defect production rate, diffusivities and GB mobility on the microstructure evolution are systematically studied. Homogenization is employed on the microstructures to obtain effective elastic constants and diffusivity as a function of fission density. The simulations provide critical insights on microstructure and property degradation in U-10Mo fuel.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Phase-field model of coupled insulator-metal transitions and oxygen vacancy redox reactions

A predominant ubiquitous feature of strongly correlated oxides is the possible presence of oxygen vacancies, which has recently been shown to have profound effects on their electronic phase transitions. Here, we formulate a comprehensive phase-field model of intercoupled insulator-metal transitions and oxygen vacancy redox reactions, taking into account the valence electron state of oxygen vacancies. We use the model to study the voltage self-oscillation phenomenon in a prototypical strongly correlated oxide, VO 2 , and discover the mutual activation of the insulator-metal transition and oxygen vacancy redox reactions leading to systematic enhancement of the oscillation frequency. Furthermore, the established methodology and the mutual activation mechanism are generally applicable to understanding any insulator-metal transition dynamics in oxygen-deficient correlated oxides and improving the performance of voltage self-oscillation-based artificial neurons.

36 MATERIALS SCIENCE↗