Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase Field”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantification and prediction of solidification textures under additive manufacturing conditions

Crystallographic textures are a major determinant of the macroscale anisotropic properties of polycrystalline metallic alloys produced in a wide range of additive manufacturing (AM) processes. Here, we introduce a statistical method that can accurately quantify the degree of orientational order of textures despite the large random fluctuations in the orientation of individual grains inherent in AM processes. The method, demonstrated for laser and resolidification of AlSi thin films, extends Z-scoring to a dynamical regime to assess the statistical significance of observed textures compared to randomly generated ones at different stages of solidification. We further show that, combined with phase-field modeling, this method can be used to infer fundamental anisotropic properties of the solid-liquid interface that are essential for texture prediction, and are compared here to the results of atomistic simulations. In addition, phase-field modeling reveals that, even at rapid AM solidification rates, the observed 〈110〉-dominated textures in the AlSi thin films are controlled predominantly by the anisotropy of the interface free-energy and sheds light on the physical mechanism of grain competition. These results significantly enhance both the existing tools for the quantification and prediction of AM crystallographic textures and our basic understanding of their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Observation of Unconventional Dynamics of Domain Walls in Uniaxial Ferroelectric Lead Germanate

Application of scanning probe microscopy techniques such as piezoresponse force microscopy (PFM) opens the possibility to re-visit the ferroelectrics previously studied by the macroscopic electrical testing methods and establish a link between their local nanoscale characteristics and integral response. It is reported that the nanoscale PFM studies and phase-field modeling of the static and dynamic behavior of the domain structure in the well-known ferroelectric material lead germanate Pb 5 Ge 3 O 11 . Here, several unusual phenomena are revealed: 1) domain formation during the paraelectric-to-ferroelectric phase transition, which exhibits an atypical cooling rate dependence; 2) unexpected electrically induced formation of the oblate domains due to the preferential domain walls motion in the directions perpendicular to the polar axis, contrary to the typical domain growth behavior observed so far; 3) absence of the bound charges at the 180 degrees head-to-head (H-H) and tail-totail (T-T) domain walls, which typically exhibit a significant charge density in other ferroelectrics due to the polarization discontinuity. This strikingly different behavior is rationalized by the phase field modeling of the dynamics of uncharged H-H and T-T domain walls. The results provide a new insight into the emergent physics of the ferroelectric domain boundaries, revealing unusual properties not exhibited by conventional Ising-type walls.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Nanoscale strain wave generation by a piezoelectric grating from polar vortices

Nanostructures formed by spontaneously broken symmetry have provided new ways to manipulate quantum states. Specifically, topological structures with periodic spatial ordering, such as polar vortices and skyrmions, can be ideal hosts for creating engineered responses in both spatial and frequency domains. So far, however, only a few examples of such hierarchical engineering have been reported in the literature. Here we demonstrate that the spatially modulated piezoelectric response of a polar vortex structure can create strain waves with a characteristic nanoscale wavefront. Using time-resolved pump–probe resonant X-ray scattering and diffraction measurements, coupled with dynamical phase-field simulations, we show that the piezoelectric modulation of the spontaneously formed polar vortex crystal functions as an acoustic diffraction grating. This system converts incoming laterally uniform strain waves into outgoing waves with a characteristic sub-terahertz frequency, driven by an intrinsic excitation of the polar vortex crystal. Moreover, our phase-field simulations suggest that the dynamic mechanical displacements exhibiting vortex textures are generated from both space- and time-varying piezoelectric responses. Our findings illustrate a new method for generating nanoscale strain waves with unique spatial textures by tuning the hierarchical order of polar topologies to engineer new collective modes, allowing for a wide range of control through the topological lattice.

ferroelectrics↗
DOE OSTI2023

Coded-mask-based X-ray phase-contrast and dark-field imaging

Phase contrast and dark-field X-ray imaging enable imaging of objects that absorb or reflect very little X-ray light. Disclosed is a method and systems for performing coded-mask-based multi-contrast imaging (CMMI). The method includes providing radiation to a coded mask that has a known phase and absorption profile according to a pre-determined pattern. The radiation is then impingent upon a sample, and the radiation is detected to perform phase-reconstruction and image processing. The method and associated systems allow for the use of maximum-likelihood and machine learning methods for reconstruction images of the sample from the detected radiation.

Shi, Xianbo↗
DOE OSTI2022

Learning two-phase microstructure evolution using neural operators and autoencoder architectures

Abstract Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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.↗
DOE OSTI2020

Controlled Nucleation and Stabilization of Ferroelectric Domain Wall Patterns in Epitaxial (110) Bismuth Ferrite Heterostructures

Ferroelectric domain walls, topological entities separating domains of uniform polarization, are promising candidates as active elements for nanoscale memories. In such applications, controlled nucleation and stabilization of domain walls are critical. In this work, using in situ transmission electron microscopy and phase-field simulations, a controlled nucleation of vertically oriented 109° domain walls in (110)-oriented BiFeO 3 (BFO) thin films is reported. In the switching experiment, reversed domains that are nucleated preferentially at the nanoscale edges of the “crest and sag” pattern-like electrode under external bias subsequently grow into a stable stripe configuration. In addition, when triangular pockets (with an in-plane polarization component) are present, these domain walls are pinned to form stable flux-closure domains. Phase field simulations show that i) field enhancement at the edges of the electrode causes site-specific domain nucleation, and ii) the local electrostatics at the domain walls drives the formation of flux closure domains, thus stabilizing the striped pattern, irrespective of the initial configuration. The results demonstrate how flux closure pinning can be exploited in conjunction with electrode patterning and substrate orientation to achieve a desired topological defect configuration. These insights constitute critical advancements in exploiting domain walls in next generation ferroelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fragmentation of faceted crystalline wires

The reduction of surface energy is a well-known driving force for the fragmentation of nanowires. Fragmentation can limit the service life of those small low-dimensional building blocks or can be controlled to induce beneficial shape changes. For isotropic surface energy, fragmentation is analogous to the classic Rayleigh-Plateau instability of liquid jets. However, commonly synthesized crystalline nanowires have strongly anisotropic surface energies and exhibit facets coinciding with cusps in the Wulff plot. Depending on growth orientation, different fragmentation behaviors have been seen in nanowires. Using phase-field simulations, we show that fragmentation of faceted nanowires with cubic crystal symmetry may occur by a finite-amplitude nonlinear instability, as opposed to a Rayleigh-Plateau-like linear instability, depending on nanowire growth orientation. Here, we carry out a weakly nonlinear analysis based on sharp interface theory to characterize the faceted nanowire shape corresponding to a surface-energy saddle point. The analysis predicts that the minimum amplitude of a periodic shape perturbation to trigger fragmentation increases with cusp strength but decreases inversely proportionally to the perturbation wavelength for long wavelengths, in good quantitative agreement with phase-field predictions. The results provide the theoretical foundation to predict nanowire stability as a function of length and surface energy in diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Néel-type antiferromagnetic order and magnetic field–temperature phase diagram in the spin-1/2 rare-earth honeycomb compound YbCl 3

Most of the searches for Kitaev materials deal with 4 d / 5 d magnets with spin-orbit-coupled J = 1 / 2 local moments such as iridates and α - RuCl 3 . Here we propose the monoclinic YbCl 3 with a Yb 3 + honeycomb lattice for the exploration of Kitaev physics. We perform thermodynamic, a c susceptibility, angle-dependent magnetic torque, and neutron diffraction measurements on YbCl 3 single crystal. We find that the Yb 3 + ion exhibits a Kramers doublet ground state that gives rise to an effective spin J eff = 1 / 2 local moment. Additionally, the compound exhibits short-range magnetic order below 1.20 K, followed by a long-range Néel-type antiferromagnetic order at 0.60 K, below which the ordered Yb 3 + spins lie in the a c plane with an angle of 16(11) ° away from the a axis. These orders can be suppressed by in-plane and out-of-plane magnetic fields at around 6 and 10 T, respectively. Moreover, the Néel temperature varies nonmonotonically under the out-of-plane magnetic fields, suggesting a reduced spin dimensionality. Finally, together with the strong in-plane magnetic anisotropy and the reduced order moment 0.8(1) μ B at 0.25 K, all indicate that YbCl 3 could be a two-dimensional spin system to proximate the Kitaev physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

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↗
DOE OSTI2023

Mechanistic FCCI model for BISON

In this report, multi-scale simulations to develop a mechanistic model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY23 are described. The model assumes that cladding wastage is dominated by the intermetallic phase (Fe,Cr)$_{17}$Nd$_2$, and calculates production and diffusion of Nd in the fuel, diffusion into the cladding, and formation and growth of the intermetallic wastage layer. Atomistic simulations are used to calculate diffusivity of Nd in the (Fe,Cr) system. Mesoscale simulations, informed by previous atomistic calculations, are used to determine the effective diffusivity of Nd through the fuel, accounting for empty and sodium-filled porosity within the fuel. The BISON model incorporates this data and calculates the wastage layer thickness as a function of time, assuming parabolic growth kinetics based on the best available data. The results of the new mechanistic BISON model are compared to existing assessment cases for the X447 experiment in EBR-II, and show good agreement with experimental data and the existing BISON empirical model. To improve predictions of kinetics for future work, a phase-field model that includes fuel, wastage, and cladding phases is developed. Preliminary testing indicates that the kinetics may differ significantly from the assumed parabolic growth law, suggesting direction for future improvement of the model that will allow it to be applied to other fuel designs with greater confidence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗