Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Field Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improvement of the $\mathrm{BISON U_3Si_2}$ modeling capabilities based on multiscale developments to modeling fission gas behavior

Uranium silicide (U 3 Si 2 ) is a concept explored as a potential alternative to UO 2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U 3 Si 2 has been previously used in research reactors in the form of dispersion fuel, but operated at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U 3 Si 2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U 3 Si 2 in the BISON fuel performance code. These lower length scale informed models are then utilized in the validation of BISON U 3 Si 2 modeling capabilities to simulate the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the validation process to identify where further experiments and lower length scale modeling would be beneficial. Here, the multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Understanding the effect of crystal anisotropy on grain growth, texturing and transport via the orthorhombic ?-U system [Slides]

The alpha phase of uranium exhibits the orthorhombic crystal structure, resulting in significant physical property anisotropy. These complex properties make the material challenging to work with in engineering settings but provide a rich arena to investigate the fundamental, multi-scale effects of anisotropy on physical behaviors such as grain growth and microstructure evolution under irradiation. We apply molecular dynamics and phase field modeling to study mass transport and grain growth in alpha-uranium. We characterize the mobilities of crystalline defects and the interfacial energy, and We find that the grain boundary energy is highly variable depending on the interface structure and that the relative rates of defect transport change with temperature. We also find that anisotropic thermal expansion has a significant impact on grain growth kinetics and texture development. Our results are used to explain experimental observations and to provide a basis for further hypotheses into the complex physical behavior of alpha-uranium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Comparison between phase field simulations and experimental data from intragranular bubble growth in UO{sub 2}

In this work, we used the phase field method to simulate the post-irradiation annealing of UO{sub 2} described in the experimental work by Kashibe et al., 1993 [1]. The simulations were carried out in 2D and 3D using the MARMOT FEM-based phase-field modeling framework. The 2-D results compared fairly well with the experiments, in spite of the assumptions made in the model. The 3-D results compare even more favorably to experiments, indicating that diffusion in all three directions must be considered to accurate represent the bubble growth. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

The effects of stress, temperature and facet structure on growth of { 10 1 ¯ 2 } twins in Mg: a molecular dynamics and phase field study

The present study focuses on the roles of structure, stress and temperature on the mobility of facets bounding twin domains and on the overall impact on the growth kinetics of { 10 1 ¯ 2 } twins in Mg. We first use molecular dynamics simulations (MD) to study the relationship between interface structures, stress, temperature and mobility. As facets can exhibit distinct atomic scale structures, we also quantify the effect of stress relaxation mediated by misfit dislocations on the mobility of the basal-prismatic interface. The information gathered for ten different facets is used to calibrate an anisotropic phase field model predicting the effects of temperature and stress on the kinetics and shape of a growing { 10 1 ¯ 2 } 3D twin, by comparison with MD data. The comparative analysis suggests that, as compared to simulations of individual facets, the effective mobility of facets can be significantly different during growth of a 3D twin. This can be attributed to the complex three-dimensional internal stress state of the twin and to coupling in facet motion triggered by facet junctions. As such, by directly fitting both the twin growth kinetics and morphology as predicted by MD for different stresses and temperatures, we can identify a limited set of rate-limiting facets that control twin growth and morphology, and measure their effective mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Uncovering the re-distribution mechanism of Ni in a de-alloyed Ni-Cr alloy in molten fluorinated salts

A mechanism of Ni redeposition during dealloying corrosion of Ni-Cr is investigated. A model Ni20Cr (wt%) metal alloy was exposed to molten LiF-NaF-KF eutectic (FLiNaK) at 600 °C and at an applied potential of +2.1 VK+/K, above the critical potential for onset of dealloying. Upon extended exposure times (up to 12 h), prominent salt-filled corrosion channels emerge along grain boundaries. A unique grain boundary corrosion mechanism, with respect to the exposed faces of the grains, a central focus of this investigation, is intrinsically connected to the formation of high purity Ni-rich de-alloyed regions within the salt-filled channel. We implement microscale techniques such as energy dispersive spectroscopy (EDS) and electron backscatter diffraction (EBSD) to uncover morphological and compositional variations in relevance to the formation of bicontinuous porosity from corrosion dealloying. To rationalize our findings, a phase-field model is developed and discovers a mechanism in which dissolved Ni from one grain can be redeposited on an adjacent grain at a given difference in interfacial energies. The interaction of chemical and structural factors at the grain boundaries plays a central role in elucidating the dynamics of this phenomenon and its implications towards corrosion.

Mills, Sean H↗
DOE OSTI2018

Effect of Coherent Lattice Mismatch on the Morphology and Kinetics of Ordered Precipitates

The morphological evolution and precipitation kinetics of ordered precipitates under the influence of lattice mismatch were studied by using the phase field model. The variation of the average aspect ratio of the precipitates and their distributions was investigated quantitatively, and the relationship of particles shape and particles fraction with the variation of the aspect ratio was determined. As the lattice mismatch increases, the coarsening rate constants and the particles number density decrease, the particles size distribution (PSD) becomes narrower, and the peak value of the PSD becomes lower than the value obtained by the Lifshitz–Slyozov–Wagner theory, while it is close to the value obtained by the Brailsford–Wynblatt theory. The change rate of the particles number density shows four-stage characteristics and decreases as the lattice mismatch and time increase, indicating that the precipitation and coarsening are retarded for a large lattice mismatch and that a large lattice mismatch results in fewer particles being dispersed in the matrix phase. The morphology and the particles number density show that a suitable lattice mismatch can enhance the mechanical properties of the alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Numerical study of the effect of the shape of the phase diagram on the eutectic freezing temperature

To evaluate the reliability of metal-carbon eutectic systems as fixed points for the next generation of high-temperature standards the effect of thermodynamic properties related to the shape of eutectic phase diagram on the freezing temperature is investigated within the context of the numerical multi-phase-field model. The partition coefficient and liquidus slopes of the two solids involved in the eutectic reaction are varied deliberately and independently. The difference between the eutectic temperature and the freezing temperature is determined in dependence of the solid/liquid (s/l) interface shape and concentration. Where appropriate reference is made to the Jackson-Hunt analytical theory. It is shown that there are mainly two typical conditions to decrease the undercooling: 1) a small liquidus slope and 2) the associated difference between the eutectic composition and the liquid composition during solidification.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Intermittent dynamics in externally driven ferroelastics and strain glasses

The interplay of elastic anisotropy and disorder dictates many of the properties of ferroic materials, specifically martensites. Here we use a phase-field model for ferroelastic athermal materials to study their response to an increasing external stress that couples to the strain order parameter. We show that these systems evolve through avalanches and study the avalanche-size distribution for ferroelastic systems (large anisotropy and/or small disorder) and for the strain glass (small anisotropy and/or large disorder) using various statistical analysis techniques, including the maximum likelihood method. The model predicts that in the former case the distribution is subcritical or power law (in agreement with experimental observations), whereas in the latter case it becomes supercritical. Finally, our results are consistent with experiments on martensitic materials, and we predict specific avalanche behavior that can be tested and used as an alternative means to characterize strain glasses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

In-situ Monitoring of Dynamic Phenomena during Solidification and Phase Transformation Processing

The purpose of this project is to: (1) Directly observe phase transformations and microstructure evolution using proton (and synchrotron x-ray) radiography and tomography; (2) Constrain phase-field models for microstructure evolution; (3) Experimentally control microstructure evolution during processing to enable co-design; and (4) Advance toward the MaRIE vision. Understand microstructure evolution and chemical segregation during solidification {yields} solid-state transformations in Pu-Ga.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Orientation Dependence of Columnar Dendritic Growth with Sidebranching Behaviors in Directional Solidification: Insights from Phase-Field Simulations

In this study, a thin-interface phase-field model was employed to study the orientation dependence of the columnar dendritic growth with sidebranching behaviors in directional solidification. It was found that the dimensionless tip undercooling increases with the increase of misorientation angle for three pulling velocities. The primary spacing is found to be a function of misorientation angle, and the dimensionless primary spacing with respect to the misorientation angle follows the orientation correction given by Gandin and Rappaz (Acta. Metall. 42:2233–2246, 1994). For the analysis of the dendritic tip, the two-dimensional (2-D) form of the nonaxisymmetric needle crystal was used to determine the radius of the tilted columnar dendrite. Based on the definitions of open side and constrained side of the dendrite, the analysis of the width active sidebranches and the dendritic area in 2-D with respect to the distance from the dendritic tip was carried out to investigate the asymmetrical dendrite envelop and sidebranching behaviors on the two sides in directional solidification. The obtained prefactor and exponent with respect to misorientation angle are discussed, showing that the sidebranching behaviors of a tilted columnar dendritic array obey a similar power-law relationship with that of a free dendritic growth.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Nanoscale origins of nonlinear behavior in ferroic thin-films

The nonlinear response of a ferroic to an applied stimulus (e.g. electric field, mechanical stress) is a fundamental characteristic that underpins a number of technologically significant applications[1-3]. It is also the driving feature in numerous physical phenomena, such as interfacial motion[4,5], spin glasses[6], relaxors[7] and phase transitions[8]. In particular, nonlinearity associated with minor hysteresis loops is an extremely useful avenue to explore energy dissipation and losses in such systems. This knowledge is necessary for the design of future materials with enhanced low-field properties. Quantitatively, the macroscopic nonlinear response of ferroic systems at low to mid-range amplitudes of driving fields is given by the phenomenological Rayleigh law[9], first conceived in 1887 for magnetic materials. Yet, the applicability of the Rayleigh law at small length scales has not been extensively studied. Here, we show using a combination of scanning probe techniques and phase field modeling, that nanoscale response appears to follow a non-Rayleigh regime. However, through statistical analysis, we find that a distribution in the individual responses can lead to directly to Rayleigh-like behavior of the strain on a macroscale. The studies shed light on the nanoscale origins of nonlinear behavior in disordered ferroics.

ferroelectric↗
DOE OSTI2023

Nondestructively Visualizing and Understanding the Mechano-Electro-chemical Origins of “Soft Short” and “Creeping” in All-Solid-State Batteries

All-solid-state Li-metal batteries (ASLMBs) represent a significant breakthrough in the quest to overcome limitations associated with traditional Li-ion batteries, particularly in energy density and safety aspects. However, widespread implementation is stymied due to a lack of profound understanding of the complex mechano-electro-chemical behavior of Li metal in the ASLMBs. Herein, operando neutron imaging and X-ray computed tomography (XCT) are leveraged to nondestructively visualize Li behaviors within ASLMBs. This approach offers real-time observations of Li evolutions, both pre- and post- occurrence of a “soft short”. The coordination of 2D neutron radiography and 3D neutron tomography enables charting of the terrain of Li metal deformation operando. Concurrently, XCT offers a 3D insight into the internal structure of the battery following a “soft short”. Despite the manifestation of a “soft short”, the persistence of Faradaic processes is observed. To study the elusive “soft short”, phase field modeling is coupled with electrochemistry and solid mechanics theory. The research unravels how external pressure curbs dendrite growth, potentially leading to dendrite fractures and thus uncovering the origins of both “soft” and “hard” shorts in ASLMBs. Furthermore, by harnessing finite element modeling, it dive deeper into the mechanical deformation and the fluidity of Li metal.

25 ENERGY STORAGE↗
DOE OSTI2020

Unexpected Giant Microwave Conductivity in a Nominally Silent BiFeO 3 Domain Wall

Nanoelectronic devices based on ferroelectric domain walls (DWs), such as memories, transistors, and rectifiers, have been shown in recent years. Practical high-speed electronics, on the other hand, usually demand operation frequencies in the gigahertz (GHz) regime, where the effect of dipolar oscillation is important. Herein, an unexpected giant GHz conductivity on the order of 10 3 S m -1 is observed in certain BiFeO 3 DWs, which is about 100 000 times greater than the carrier-induced direct current (dc) conductivity of the same walls. Surprisingly, the nominal configuration of the DWs precludes the alternating current (ac) conduction under an excitation electric field perpendicular to the surface. Theoretical analysis reflects that the inclined DWs are stressed asymmetrically near the film surface, whereas the vertical walls in a control sample are not. The resultant imbalanced polarization profile can then couple to the out-of-plane microwave fields and induce power dissipation, which is confirmed by the phase-field modeling. Since the contributions from mobile-carrier conduction and bound-charge oscillation to the ac conductivity are equivalent in a microwave circuit, the research on local structural dynamics may open a new avenue to implement DW nano-devices for radio-frequency applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Chemomechanical origin of morphological disparity in lithium metal electrodeposition

The morphology of electrochemically deposited lithium (Li) critically governs the cycling stability and safety of Li metal batteries, yet the underlying controlling factors remain poorly understood. Even within the same coin cell, Li deposits can exhibit strikingly different morphologies, for example, sparse whiskers coexisting with particle-like deposits, indicting strong local variations in growth conditions. By combining cryogenic transmission electron microscopy with phase-field modeling, we here identify the root cause of this disparity. We reveal that the morphological divergence originates from the variations in structure and chemical composition of the initial SEI layer formed on the copper (Cu) substrate. These variations arise from substrate defects that induce localized electrochemical fields, driving selective dissociation of either solvent or salt, producing SEI layers with distinct chemomechanical properties. Solvent-derived organic SEI favors whisker growth, whereas salt-derived SEI promotes particle formation. This study establishes a direct mechanistic link between SEI chemomechanical properties and Li morphology, providing design principles for tailoring SEI layers to control Li deposition

energy storage↗
DOE OSTI2026

Decoding THz‐Driven Dynamic Fingerprints of Ferroelectric Nanotwin Networks

Ultrafast polarization dynamics in ferroelectrics are of considerable interest for high-speed tunable dielectrics and electro-optics. Extended domain wall networks formed in ferroelectric twin nanodomains can support collective dynamics in the terahertz regime but require techniques that track polarization and strain evolution driven by ultrafast stimulus. Here, we use multi-modal probing of THz-pulse-driven excitations in PbTiO 3 /SrTiO 3 superlattices by combining X-ray free electron laser measurements that directly tracks lattice changes, with optical second harmonic generation that tracks the electronic potential coupled with the lattice potential. Dynamical phase-field modeling enables fingerprinting of these collective modes as superpositions of domain “breathing” through wall oscillations and polarization “rotations” with still walls. Ultrafast domain wall motion at 0.1–0.5 THz is observed at practical fields of 100 kV/cm with wall velocities of >4000 m/s, approaching typical speed of sound in PbTiO 3 . A unique “charging” mode is discovered that can electrically charge and discharge domain walls on ∼4 ps time scale thus dynamically tuning wall conductivity. Integrated experimental and theoretical fingerprinting of the dynamical landscape presented here enables ultrafast control of ferroics for high-speed microelectronics and optical applications.

THz dynamics↗