Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Phase-field-guided design of record-high piezoelectricity and discovery of simultaneous high light transparency and high piezoelectricity in relaxor ferroelectrics

Here, the phase-field method has been extensively applied to predicting the domain structures and their responses to external fields and understanding experimentally observed domain states under different electromechanical conditions in ferroelectric heterostructures. This article highlights the successful examples of phase-field applications in guiding the design of relaxor ferroelectric ceramics and crystals with record-high piezoelectricity and the discovery of simultaneous high light transparency and piezoelectricity of relaxor ferroelectric crystals for a wide range of biomedical, robotics, and optoelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing diffraction patterns of microstructures from phase-field simulations

The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field modeling of aging-induced microstructure evolution in pentaerythritol tetranitrate thin films and ramifications for shock initiation

Aging of energetic materials may change performance and affect their safety and reliability, but the relationship between microstructure changes induced by aging and consequent performance changes has not been fully established. This work presents results of phase-field method simulations used to model microstructure evolution of vapor-deposited pentaerythritol tetranitrate (PETN) thin films. Simulated aging is shown to induce grain coarsening and substantial changes of the configuration of porosity in the film: Specifically, we show that porosity tends to concentrate in large pores to a greater degree in aged films, a state that is arrived at by closure or consolidation of small pores. To evaluate the performance of the as-deposited and aged films, we perform two-dimensional hydrocode flyer-film impact simulations that incorporate the phase-field output microstructures directly, permitting us to connect features therein to changes in reactivity, a key metric of energy output for shock initiation. The results demonstrate that declining sensitivity obtained for the simulated aged films can be correlated with the loss of fine-structured pores relatively early in the aging process, while long-term microstructure evolution that gradually alters the shape of large, branching pores is less impactful. Finally, we discuss commonalities and discrepancies between our simulation results and high-throughput initiation experiments on shock initiation of aged PETN thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerate microstructure evolution simulation using graph neural networks with adaptive spatiotemporal resolution

Abstract Surrogate models driven by sizeable datasets and scientific machine-learning methods have emerged as an attractive microstructure simulation tool with the potential to deliver predictive microstructure evolution dynamics with huge savings in computational costs. Taking 2D and 3D grain growth simulations as an example, we present a completely overhauled computational framework based on graph neural networks with not only excellent agreement to both the ground truth phase-field methods and theoretical predictions, but enhanced accuracy and efficiency compared to previous works based on convolutional neural networks. These improvements can be attributed to the graph representation, both improved predictive power and a more flexible data structure amenable to adaptive mesh refinement. As the simulated microstructures coarsen, our method can adaptively adopt remeshed grids and larger timesteps to achieve further speedup. The data-to-model pipeline with training procedures together with the source codes are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Marmot V2

MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.

Aagesen, LarryK.↗
DOE OSTI2023

Phase Field Sintering Simulations of Tagged UO 2

Isotopic taggants are being considered as an additive for UO 2 manufacturing to improve nuclear security. Before taggants can be adopted, it must be determined whether taggants have any effect on the fuel microstructure or performance. As part of this effort, a grand potential formulation of the phase field method was used to model untagged, Cr-tagged, and Fe-tagged UO 2 . The simulations were compared to determine whether the taggants had any effect on fuel density or grain size. The results suggest that the taggants will not affect the density or grain size, but these findings may not be reliable because of high uncertainty. Causes of the uncertainty were discussed, and future work is proposed to reduce the model uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Multiscale, mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, the effective diffusion coefficient of the fission product Ag through the grain boundary (GB) network is calculated using a combination of atomistic and phase-field methods. Atomistic calculations of Ag diffusivity in SiC bulk and GBs are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Since GBs serve as pathways for Ag diffusion, Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Understanding the effect of irradiation and silicon carbide microstructure on silver diffusion in TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, a combination of atomistic and phase-field methods focuses on quantifying the effect of irradiation and SiC microstructure on Ag diffusion. Atomistic calculations provide the Ag diffusivity in SiC bulk and GBs, also providing an estimate of the effect of irradiation defects on Ag diffusivity in bulk SiC. These lower length scale efforts are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 and AGR-2 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Accelerating phase-field-based microstructure evolution predictions via surrogate models trained by machine learning methods

Abstract The phase-field method is a powerful and versatile computational approach for modeling the evolution of microstructures and associated properties for a wide variety of physical, chemical, and biological systems. However, existing high-fidelity phase-field models are inherently computationally expensive, requiring high-performance computing resources and sophisticated numerical integration schemes to achieve a useful degree of accuracy. In this paper, we present a computationally inexpensive, accurate, data-driven surrogate model that directly learns the microstructural evolution of targeted systems by combining phase-field and history-dependent machine-learning techniques. We integrate a statistically representative, low-dimensional description of the microstructure, obtained directly from phase-field simulations, with either a time-series multivariate adaptive regression splines autoregressive algorithm or a long short-term memory neural network. The neural-network-trained surrogate model shows the best performance and accurately predicts the nonlinear microstructure evolution of a two-phase mixture during spinodal decomposition in seconds, without the need for “on-the-fly” solutions of the phase-field equations of motion. We also show that the predictions from our machine-learned surrogate model can be fed directly as an input into a classical high-fidelity phase-field model in order to accelerate the high-fidelity phase-field simulations by leaping in time. Such machine-learned phase-field framework opens a promising path forward to use accelerated phase-field simulations for discovering, understanding, and predicting processing–microstructure–performance relationships.

36 MATERIALS SCIENCE↗
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 OSTI2024

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗