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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

Modeling solidification microstructures has become an area of intense study in recent years. The properties of large scale cast products, ranging from automobile engine blocks to aircraft components and other industrial applications, are strongly dependent on the physics that occur at the mesoscopic and microscopic length scales during solidification. The predominant morphology found in solidification microstructures is the dendrite, a tree-like pattern of solid around which solidification proceeds. The microscopic properties of cast products are determined by the length scales of these dendrites, and their associated segregation profiles. For this reason understanding the mechanisms for pattern selection in dendritic growth has attracted a great deal of interest from the experimental and theoretical communities. In particular, a great deal of research has been undertaken to understand such issues as dendrite morphology, shape and growth speed. Experiments on dendrite evolution in pure materials by Glicksman and coworkers on succinonitrile (SCN), and more recently pivalic acid (PVA), as well as other transparent analogs of metals, have provided tests of theories for dendritic growth, and have stimulated considerable theoretical progress. These experiments have clearly demonstrated that in certain parameter ranges the physics of the dendrite tip can be characterized by a steady value for the dendrite tip velocity, radius of curvature and shape. Away from the tip, the time-dependent dendrite exhibits a characteristic sidebranching as it propagates, which is not yet well understood. These experiments are performed by observing individual dendrites growing into an undercooled melt. The experiments are characterized by the dimensionless undercooling. Most experiments are performed at low undercooling.

Dantzig, Jonathan A.↗
DOE OSTI2021

Non-orthogonal computational grids for studying dislocation motion in phase field approaches

In this work, new non-orthogonal computational grids are implemented into a phase field model called Phase Field Dislocation Dynamics (PFDD). We demonstrate that the new non-orthogonal grid can accommodate multiple slip planes in either the face centered cubic (FCC) or body centered cubic (BCC) crystallographic systems. Here, we show that they avoid numerical errors induced when modeling glide on inclined slip planes in an orthogonal grid. The Gibbs effect that arises in the orthogonal or rotated orthogonal grids is substantially diminished when a non-orthogonal grid is employed. A few test cases demonstrate the effectiveness of using non-orthogonal grids in solving systems with multiple non-planar slip systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗
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 OSTI2022

Phase-field dislocation modeling of cross-slip

The phase-field dislocation dynamics (PFDD) model is extended to simulate cross-slip in body-centered cubic materials. The model makes use of a non-orthogonal, body-centered cubic numerical grid such that all calculable points lie on dislocation glide planes. We demonstrate that the advanced PFDD model predicts a low-energy, non-planar core for screw dislocations and a planar core for edge dislocations. These core differences result in screw-edge character dependence in the critical stress to initiate glide and dislocation drag coefficients. Finally, we show that this model can be used to study the propensity of screw dislocations to cross-slip around or transmit through a coherent, crystallographic obstacle, depending on obstacle strength.

36 MATERIALS SCIENCE↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗