Search NASASearch

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 127 records · Page 7

A brief review on strain engineering of ferroelectric K x Na 1- x NbO 3 epitaxial thin films: Insights from phase-field simulations

Strains play a pivotal role in determining the phase equilibrium, domain configuration, and functional properties of the low-dimensional ferroelectrics. There is growing interest in the strain engineering of ferroelectric K x Na 1- x NbO 3 (KNN) epitaxial thin films, which exhibit excellent physical properties and promise as eco-friendly alternatives to lead-based ferroelectrics for microdevice applications. Further, advances have been made in understanding the phase equilibria and transitions, domains and domain walls, and their relations to the physical properties of KNN epitaxial thin films using a combination of experiments and theoretical modeling, particularly phase-field simulations. Here, we review recent progress in these aspects and showcase the phase-field method for establishing strain phase diagrams, elucidating the domain and domain wall structures at equilibrium, and predicting the structure–property relationships in ferroelectric KNN thin films. We also discuss challenges and opportunities to further advance our understanding of KNN thin films and potentially unlock new functionalities by leveraging phase-field simulations.

36 MATERIALS SCIENCE

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

Close Proximity to a Quantum Phase Transition in TmZn2⁢GaO5

TmZn2⁢GaO5 is a newly synthesized triangular lattice magnet that exhibits a unique quantum phase characterized by strong Ising anisotropy, a pseudodoublet crystal electric field ground state, and a low-energy gapped excitation at the 𝐾 point. Unlike its well-known counterparts, TmMgGa⁢O4 and YbMgGa⁢O4, this material crystallizes in a distinct hexagonal structure, leading to a platform for investigating frustrated magnetism. Magnetic susceptibility, heat capacity, and inelastic neutron scattering measurements confirm the absence of long-range magnetic order down to 50 mK, placing TmZn2⁢GaO5 in a distinct region of the transverse-field Ising model phase diagram. Theoretical calculations based on spin-wave theory and mean-field modeling reproduce key experimental observations, reinforcing the material’s placement in a quantum disordered or multipolar state near a quantum critical point that separates it from the low-energy phase of TmMgGa⁢O4. These results highlight its potential for exploring quantum disordered states, anisotropic excitations, and exotic quantum phases in frustrated spin systems.

Ennis, Matthew [Department of Physics, Duke Univer

Analytical prediction of the interior noise for cylindrical models of aircraft fuselages for prescribed exterior noise fields. Phase 1: Development and validation of preliminary analytical models

The basic theoretical work required to understand sound transmission into an enclosed space (that is, one closed by the transmitting structure) is developed for random pressure fields and for harmonic (tonal) excitation. The analysis is used to predict the noise reducton of unpressurized unstiffened cylinder, and also the interior response of the cylinder given a tonal (plane wave) excitation. Predictions and measurements are compared and the transmission is analyzed. In addition, results for tonal (harmonic) mechanical excitation are considered.

Pope, L. D.

Resolution requirements for high-speed droplet impact along a rigid wall

The mechanical and thermal loads that are produced on high-speed projectiles from a droplet are key in the projectile development.To accurately predict the pressure produced along the wall surface during impact, it is important to correctly model the air/water interface of the droplet. Novel modeling approaches such as the Phase-Field enable control of the (numerical) interface thickness. However, the interplay between the grid resolution (including potential adaptive mesh refinement) and the interface thickness is unclear. In this work, the resolution requirements for numerical simulations of a two-dimensional cylindrical water droplet impacting a rigid wall at Mach four are computationally investigated. Numerical simulations are performed using a second-order accurate scheme with adaptive mesh refinement and a Phase-Field interface model. The choice of numerical interface thickness and resolution parameters affects the maximum pressure and location of the maximum pressure along the wall. The overall maximum pressure does not converge even at high resolutions. The number of levels and thickness (number of cells across the interface) play a role in obtaining a well behaved solution.

Burrell, E.

Assessment of effective elastic constants of U-10Mo fuel: A multiscale modeling and homogenization study

The significant microstructural changes that U-Mo fuel undergoes during operation degrades its mechanical properties and structural integrity. Microstructural evolution entails the formation, evolution, and redistribution of porosity in conjunction with grain refinement. In the present paper, we employ numerical approaches to assess the impact of the various microstructural features—grains, nanoscale intragranular fission gas bubbles, and mesoscale intergranular voids—on the degradation of elastic constants. Phase-field microstructure models are combined with the asymptotic expansion homogenization technique in order to derive the effective elastic constants as a function of porosity and fission density. Here the results are verified and compared against theoretical bounds. Using this approach, elastic degradation in operating nuclear fuels can be quantified when the distributions of microstructural features are known.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Part-scale microstructure prediction for laser powder bed fusion Ti-6Al-4V using a hybrid mechanistic and machine learning model

Laser powder bed fusion (LPBF) Ti-6Al-4V is widely studied for use in structural applications in aerospace and medical industries, but mechanical anisotropy and microstructural inhomogeneity prohibits its wider adoption. Although successful microstructure prediction models have been developed, a remaining challenge is their limited integration across length/time scales and validation by experimental studies. Here, this work proposes a physics-augmented machine learning surrogate model to unite predictions of LPBF temperature, β phase morphology and texture, and α/α’ formation into a single framework that is calibrated and validated with experiments. First, a phase field (PF) model of the martensitic β→α’ transformation is developed and calibrated using data from in-situ synchrotron cyclic heating/cooling studies quantifying the variation of α phase fraction with time. In parallel, an established finite difference-Monte Carlo (FDMC) model predicts the part-scale temperature profile and β grain formation during solidification. A dataset is developed using LPBF cyclic temperature descriptors from the FDMC model as inputs and corresponding α/α’ phase fraction and width from the PF model as outputs. Five machine learning (ML) regression models are tested and optimized, having mean absolute error in testing ≤ 4 %, and the k-nearest neighbors (KNN) model is selected as the best performing. The KNN model is called at the nodal level during post-processing of the FDMC model to replace and downscale the response of the PF model. The combined agility and accuracy of the hybrid FDMC-ML model enables part-scale microstructure predictions that can be further used for property predictions to accelerate AM process optimization.

36 MATERIALS SCIENCE

Summary of nozzle-exhaust plume flowfield analyses performed in support of space shuttle applications

Results are summarized of various analyses of the exhaust plumes of the space shuttle. Specific topics discussed included: the development of the two-phase plume flow field model, including finite rate chemistry and free molecular effects; analyses of exhaust plume simulation studies; verification of the analytical two-phase plume flow field model; and complete documentation of the two-phase plume computer code.

Smith, S. D.

Improving Additive Manufactured Component Performance through Multi-Scale Microstructure Simulation and Process Optimization

The purpose of this project was to utilize computational tools to understand the relationships between processing, microstructure, and properties for additively manufactured (AM) aluminum alloys for automotive applications, and to provide an engineering solution for helping to optimize process conditions. The project leverages ORNL developments in computational modeling, including AM process modeling, phase-field based microstructure evolution predictions, and data analytics techniques for mapping process conditions to material outcomes. The project utilized an Al-Cu-Mn-Zr alloy as a model material for studying formation of defects and microstructural features in response to variations in process conditions. Based on both pre-existing experimental data and simulation results, statistical process maps were constructed to identify regions of process space with minimal defect formation and advantageous microstructures and properties. The software tools used for this purpose were successful disseminated to GM, who were able to successful compile the relevant HPC codes within their own computing ecosystem and perform initial calculations to reproduce ORNL results.

36 MATERIALS SCIENCE

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS

3D interface size effects on slip transfer in Ti/Nb nanolaminates

Two-phase nanolaminates are well-renowned for achieving extraordinarily high strengths but at the sacrifice of reduced toughness and strain to failure. Recently ”thick” interfaces, or so called 3D interfaces, in Cu/Nb nanolaminates were experimentally shown to improve both of these mechanical properties. Here, in this work, we study the effect of 3D interfaces in the hexagonal close packed (HCP)/body centered cubic (BCC) Ti/Nb nanolaminate system. Nanoindentation hardness testing suggests increased strength with the introduction of a 3D Ti–Nb interface and a positive size effect with increases in 3D interface thickness from 5 nm to 20 nm. To understand this effect from a single dislocation perspective, we present a phase-field dislocation dynamics (PFDD) model for multi-phase HCP/BCC systems. We employ the model to simulate stress-driven transfer of single dislocations across 3D Ti/Nb interfaces of various thicknesses. Our results show that the critical stress for slip transfer increases with the thickness of the interface. This positive size effect is stronger for transfer from basal or prismatic dislocations in the Ti layer to 110$\langle$111$\rangle$ dislocations in the Nb layer than the reverse. For this Ti/Nb system, a critical thickness of 2 nm is identified at which the asymmetry in slip transfer is minimized. This work showcases 3D interfaces as a beneficial microstructure modification to strengthen as well as reduce anisotropy in nanocrystalline materials containing HCP phases.

Dislocations

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE

Mesoscale Modeling of Dislocation Cell Structure Evolution and Radiation-Induced Segregation in Additively Manufactured Austenitic Stainless Steel

Structural alloys under irradiation develop radiation-induced segregation (RIS) at point defect sinks, which lead to undesired changes in the alloy's properties. Additively manufactured austenitic stainless steels are expected to yield a distinct response to irradiation damage owing to their unique as-printed dislocation cell structure. In this talk, we present the development of a mesoscale model for dislocation cell structure evolution coupled with RIS using the MOOSE framework. The evolution of network dislocations within spatially distinct cells and cell walls are modeled using climb-mediated edge annihilation and generation processes. Starting with initial microstructures of as-printed non-equilibrium segregations at dislocation cell walls, our 2D and 3D simulations explore the effects of irradiation temperature, dose rate, dislocation sink bias, and specimen thickness on microstructure evolution. The results will be compared against experimental characterizations of in-situ and ex-situ ion-irradiated samples, and the implications of evolving cell structure on irradiation damage response will be discussed.

additive manufacturing

In situ irradiation of spent nuclear fuels

To improve the economics of commercial nuclear reactors, nuclear vendors and utilities are seeking approval for increased burnup usage of the existing nuclear fleet. This is especially critical for meeting the clean energy initiative by the United States Government, Department of Energy (DOE). However, one of the key challenges the nuclear industry faces in this regard is that the regions exposed to high burnup and low temperatures exhibit a fine-grained microstructure with large bubbles known as high-burnup structure (HBS) [1]. The formation of HBS has been correlated to the diminished performance of the reactor, as well as fuel fragmentation and pulverization during transient and accidental conditions [2]. Therefore, it is paramount to understand the mechanisms for HBS formation along with its impact on the properties and performance of nuclear fuels. While existing programs, such as Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Fuel Campaign (AFC) are focusing on evaluating the performance impact of HBS, the physical mechanisms contributing to HBS formation are still not fully understood. In addition, having predictive capabilities and sound understanding of the microstructural evolution of nuclear fuel is essential for accelerated development, qualification, and deployment of new nuclear materials and novel reactor designs for advanced nuclear reactors. There is a lack of consensus among researchers regarding the mechanisms leading to such restructuring observed in HBS. Grain subdivision due to polygonization versus recrystallization, continuous versus discrete recrystallization occurring in tandem or conjunction, etc., have been proposed and debated. In general, it is hypothesized that defect accumulation and dislocation interaction within the grains cause the realignment of dislocations into grain boundaries, leading to the new subgrain formation, which over time transforms into new grains. However, due to the lack of transient data, the importance of fission rate, irradiation, thermal, and stress history of the fuel on the restructuring could not be assessed. In situ microstructural evolution under various irradiation conditions is desired to bridge this gap. Alternatively, phase-field-based models have been developed to capture HBS formation via discrete recrystallization utilizing the classical nucleation approach [3–5]. However, in these models, the grain nucleation criteria are often defined based on empirical relations for burnup and fission gas density leading to dislocation density change. A mechanistic approach to capture the dislocation interaction with the microstructural features leading to grain subdivision is lacking.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Thermodynamic origin of nonvolatility in resistive memory

Electronic switches based on the migration of high-density point defects, or memristors, are poised to revolutionize post-digital electronics. Despite significant research, key mechanisms for filament formation and oxygen transport remain unresolved, hindering our ability to predict and design device properties. For example, experiments have achieved 10 orders of magnitude longer retention times than predicted by current models. Here, using electrical measurements, scanning probe microscopy, and first-principles calculations on tantalum oxide memristors, we reveal that the formation and stability of conductive filaments crucially depend on the thermodynamic stability of the amorphous oxygen-rich and oxygen-poor compounds, which undergo composition phase separation. Including the previously neglected effects of this amorphous phase separation reconciles unexplained discrepancies in retention and enables predictive design of key performance indicators such as retention stability. Furthermore, this result emphasizes non-ideal thermodynamic interactions as key design criteria in post-digital devices with defect densities substantially exceeding those of today’s covalent semiconductors.

36 MATERIALS SCIENCE

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE