Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theory and phase-field simulations of electrical control of spin cycloids in a multiferroic

Cycloidal spin orders are common in multiferroics. One of the prototypical examples is BiFeO 3 (BFO) which shows a large polarization and a cycloidal antiferromagnetic (AFM) order at room temperature. Here, we employ Landau theory and phase-field simulations to analyze the coupled switching dynamics of polarization and cycloidal AFM orders in BFO. We identify 14 types of transitional spin structures between two cycloids and nine electric-field-induced spin switching paths. We demonstrate the electric-field-induced rotation of wave vectors of the cycloidal spins and discover two types of cycloidal spin switching dynamics: fast local spin flips and slow rotation of wave vectors. Also, we construct roadmaps to achieve the switching between any two spin cycloids through multistep applications of electric fields. Furthermore, this paper provides a theoretical framework for the phenomenological description of spin cycloids and a fundamental understanding of the switching mechanisms to achieve electrical control of magnetic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2022

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase-field simulations of intergranular fission gas bubble behavior in U 3 Si 2 nuclear fuel

U 3 Si 2 is a potential accident-tolerant fuel that shows promise due to its high thermal conductivity and higher uranium density relative to UO 2 . However, its swelling and fission gas release behavior in light water reactor (LWR) conditions is relatively unknown. To provide mechanistic insight and determine parameters for engineering-scale fuel performance modeling of pellet-form U 3 Si 2 , phase-field simulations of the growth, interconnection, and venting of intergranular fission gas bubbles were performed. The fractional coverage of the grain boundary and the fraction of bubble area that is vented were calculated as a function of time. From the simulation data, the fractional grain boundary coverage at saturation, an important parameter needed in engineering-scale modeling of swelling and fission gas release, was determined. Multiple simulations were run to determine the uncertainty in the calculated value. Finally, the effect of model assumptions and input parameters that are not well known was evaluated. Simulation results are compared to related theoretical and computational work. Based on the simulation results, a value of 0.60 for the fractional grain boundary coverage at saturation is recommended for U 3 Si 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Optimizing Piezoelectric Nanocomposites by High‐Throughput Phase‐Field Simulation and Machine Learning

Abstract Piezoelectric nanocomposites with oxide fillers in a polymer matrix combine the merit of high piezoelectric response of the oxides and flexibility as well as biocompatibility of the polymers. Understanding the role of the choice of materials and the filler‐matrix architecture is critical to achieving desired functionality of a composite towards applications in flexible electronics and energy harvest devices. Herein, a high‐throughput phase‐field simulation is conducted to systematically reveal the influence of morphology and spatial orientation of an oxide filler on the piezoelectric, mechanical, and dielectric properties of the piezoelectric nanocomposites. It is discovered that with a constant filler volume fraction, a composite composed of vertical pillars exhibits superior piezoelectric response and electromechanical coupling coefficient as compared to the other geometric configurations. An analytical regression is established from a linear regression‐based machine learning model, which can be employed to predict the performance of nanocomposites filled with oxides with a given set of piezoelectric coefficient, dielectric permittivity, and stiffness. This work not only sheds light on the fundamental mechanism of piezoelectric nanocomposites, but also offers a promising material design strategy for developing high‐performance polymer/inorganic oxide composite‐based wearable electronics.

42 ENGINEERING↗
DOE OSTI2021

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dendrite-resolved, full-melt-pool phase-field simulations to reveal non-steady-state effects and to test an approximate model

In this report we study the epitaxial, columnar growth of (multiply oriented) dendrites/cells for a spot melt in a polycrystalline Al–Cu substrate using two-dimensional, phase-field, direct numerical simulations (DNS) at the full-melt-pool scale. Our main objective is to compare the expensive DNS model to a much cheaper but approximate “line” model in which a single-crystal phase-field simulation is confined to a narrow rectangular geometry. To perform this comparison, we develop algorithms that automatically extract quantities of interest (QoIs) from both DNS and line models. These QoIs allow us to quantitatively assess the assumptions in the line model and help us analyze its discrepancy with the DNS model. We consider four sets of heat source parameters, mimicking welding and additive manufacturing conditions, that create a combination shallow and deep melt pools. Our largest DNS simulation used 16K x 14K grid points in space. Our main findings can be summarized as follows. Under AM conditions, the QoIs of line models are in excellent agreement with the full DNS results for both shallow and deep melt pools. Under welding conditions, the primary spacing of the DNS model is smaller than the prediction of line model. We identify a geometric crowding effect that accounts for the discrepancies between the DNS and line models. We propose two potential mechanisms that determine the response of the microstructure to geometric crowding.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field simulations of fission gas bubbles in high burnup UO2 to inform engineering-scale fuel performance modeling

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). The cause of pulverization has been hypothesized to be the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted. Simulation results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of a grain growth model for U 3 Si 2 using experimental data, phase field simulation and molecular dynamics

The purpose of this work is to develop a model for normal grain growth in U 3 Si 2 . The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. Here, the mobility was determined by comparing phase field grain growth simulations to the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase field simulation of anode microstructure evolution of solid oxide fuel cell through Ni(OH) 2 diffusion

Microstructure evolution in a solid oxide fuel cell anode is strongly affected by operating condition. The detailed mesoscale evolution mechanism under operating condition is still under debate. Here we develop a phase-field model to simulate microstructure evolution through formation and diffusion of gaseous Ni(OH)2. We studied the coarsening kinetics and redistribution of Ni under different steam distributions and compare it to that under pure hydrogen condition. Overall, the results suggest that although the presence of a steady gradient of steam leads to redistribution of Ni, Ni(OH) 2 formation and diffusion do not significantly change the Ni coarsening rate under strictly humid conditions, contrary to commonly reported hypotheses. It is concluded that competing mechanisms other than Ni(OH) 2 diffusion must be responsible for the experimentally observed Ni redistribution and enhanced coarsening under humid conditions.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of HPC based phase field simulations tool for modification of alloy morphology to enhance material properties during additive manufacturing (AM) process

Currently, additive manufacturing (AM) technology is extensively used to manufacture aerospace components because of the unique advantages of the technique in producing custom designs, complex geometries in a near-net fashion that significantly reduces the manufacturing cost, and improved energy efficiency. However, the thermodynamic properties of existing aerospace alloys and the thermal conditions that prevail in most of the powder-bed based additive manufacturing processes lead to undesirable columnar / columnar-dendritic microstructures that promote solidification cracking during the process, and the development of anisotropy in the mechanical properties of the as-built components. Time consuming, energy-intensive post AM treatments are required to recover the mechanical properties of the wrought components for which the existing alloy compositions were designed for. However, novel alloy compositions, that exploit the AM thermal conditions to produce fine, equiaxed grain structures in the as-built condition, with mechanical properties equal to or exceeding those of the wrought work-horse alloy, Ti-6Al-4V has recently been demonstrated in Ti-Cu alloys. Significant energy savings up to 66% can be obtained for a typical aerospace component through digitally designing the AM process, process parameters, and alloy composition. The current project focused on developing a fundamental understanding of the evolution of the columnar to equiaxed transition (CET) occurring during solidification under AM thermal conditions in Ti-Cu and Ti-Cu-X alloys using high-fidelity phase field (PF) simulations using the leadership class computing facilities that exist in the national laboratories. The simulations were able to capture the effect of alloy and process conditions on CET, and clearly showed the beneficial effect of a ternary solute addition to Ti-Cu binary alloys on CET. The simulations indicated that bulk nucleation in the liquid ahead of the solidifying epitaxial front was promoted by a large, transient, thermal undercooling promoted by the mismatch between the growth rate of the dendrite tips in the epitaxial front and the imposed solidification rate. AM experiments performed at Raytheon Technologies Research Center using the same thermal conditions used in the PF simulations indicated the formation of equiaxed grains, although the grain size was bigger than the ones reported in the literature. On the other hand, the equiaxed grain size predicted by PF simulations were smaller than the grain size of Ti-Cu alloys reported in the literature. The industry will follow up the existing work in a future program that would involve further optimization of the process for the Ti-Cu-X alloys to demonstrate the application of the process and the alloy to a specific aerospace component.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The effects of dose, dose rate, and irradiation type and their equivalence on radiation-induced segregation in binary alloy systems via phase-field simulations

Radiation-induced segregation is a phenomenon commonly observed in many alloys which consists of the redistribution of elements (solute or interstitial impurities) under irradiation. The onset and development of radiation-induced segregation can only occur when a sufficient flux of defects is sustained and defect sinks are present. Irradiation dose, dose rate, and particle types all affect defect flux. In this work, we employ a phase-field model to examine the effects of dose, dose rate, and type of incident particles on radiation-induced segregation behavior in a model binary alloy. The phase-field model takes into account the formation and evolution of point defects as well as defect clusters, the diffusion and clustering of alloy species, the presence of additional extrinsic defect sinks in the form of dislocations, and two different methods of radiation-damage insertion, which are intended to simulate either light-ion/electron irradiation via Frenkel pairs or heavy-ion irradiation in the form of cascades. Our results show a dose-rate and particle-type dependence on the amount of solute segregation. We show that the material systems exposed to higher dose rates are less subjected to solute segregation at equivalent doses. We also show that such dose-rate-dependence behavior is due to a delay of the incubation dose at which radiation-induced segregation effectively starts. Particle type and the presence of dislocations can accentuate this behavior. Our model predictions correlate with many experimental observations made over the years on radiation-induced segregation providing credence to the simulation results. The methodology presented in this study allows for a first-order prediction of the dose rate at which proxy irradiation experiments could be performed to approximate radiation-induced segregation behaviors seen in targeted irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗