Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field modeling”

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
NASA NTRSDate not supplied

Evolution of Dendritic Extended 3D Patterns during Directional Solidification: Microgravity Experiments in DECLIC-DSI onboard ISS and Phase-field Modeling

To characterize the dynamical formation of three-dimensional (3D) arrays of dendrites under diffusive growth conditions, in situ monitoring of a series of experiments on a transparent succinonitrile – 0.46 wt% camphor model alloy was carried out under low gravity in the DECLIC Directional Solidification Insert onboard the International Space Station. These experiments offer continuous interface observation and enable the construction of space-time evolution maps of dendrite location and primary spacing during directional solidification. We present the evolution of dendritic extended 3D patterns for a large range of parameters displaying different levels of side branching in both experiments and phase-field simulations. The combined and separate effects of macroscopic interface curvature and crystal mis-orientation are evidenced and detailed. In addition, the aligned dendritic array caused by the invasion of dendrites formed at the edge of the crucible border is investigated.

Kaihua Ji↗
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 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 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 OSTI2024

Comparison of excess free energy at an interface according to the applied interpolation scheme for elasticity: A phase-field method

Phase-field modeling is an effective simulation technique for modeling microstructure evolution of elastically anisotropic systems. To introduce the elastic energy contribution in a phase field model, an interpolation scheme is used to define the mechanical properties within the phases and across the continuous interface. Several existing interpolation schemes introduce a potential excess elastic energy at the interface, which undesirable effect on microstructure evolution needs to be evaluated. In this study, we focused on three interpolation schemes including Khachaturyan’ scheme (KHS), Voigt–Taylor’s scheme (VTS), and Steinbach–Apel’s scheme (SAS). Comparisons of these schemes’ performances were performed in three configuration types using the MOOSE (Multiphysics Object-Oriented Simulation Environment) framework: bi-crystal, isotropic particle-matrix and anisotropic particle-matrix. The contribution of excess elastic energy on the interface energy as a function of interface width and the computational time to steady-state were evaluated in these three configurations. SAS introduces the lowest excess elastic energy contribution and the VTS has the biggest contribution amongst the considered schemes. Moreover, when modeling precipitation in an anisotropic elastic material, the SAS approach seems to predict more physical convex shapes during growth, making it preferable to KHS and VTS. Finally, as currently implemented, SAS requires the largest computational time and KHS requires the smallest time to reach steady-state amongst the considered schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗