Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 721 records · Page 40

E-field controlled phase transformation in bismuth ferrite thin films, and effect of laser energy density

Heterostructures of BiFeO 3 /Pb(Mg 1/3 Nb 2/3 )O 3 -30%PbTiO 3 (BFO/PMN-30PT) were fabricated here by pulsed laser deposition. Electric fields (E) were applied on the PMN-30PT substrates along the in-plane and out-of-plane directions, and the induced strain of BFO was measured as a function of E by X-ray diffraction. Reciprocal space maps along the (110) zone axis evidenced that a monoclinic M A phase was induced in both PMN-PT and BFO, where the dominate domain structures could be switched between biaxial directions under application of in-plane E. Studies were done for films prepared under different laser energy densities (LED), in order to optimize the in-plane strain of BFO. A relatively low LED of 1 J/cm 2 produced a better BFO thin film with larger in-plane strain.

36 MATERIALS SCIENCE↗

Study of phase relationships in the Sr{sub 3}(PO{sub 4}){sub 2}–CePO{sub 4} system. Phase diagram and thermal characteristics of phases

A diagram representing phase relationships in the Sr{sub 3}(PO{sub 4}){sub 2}–CePO{sub 4} phosphate system has been developed on the basis of results obtained by thermal analysis (DTA/DSC/TGA) and X-ray diffraction (XRD) methods. One intermediate compound with the formula Sr{sub 3}Ce(PO{sub 4}){sub 3} occurs in the Sr{sub 3}(PO{sub 4}){sub 2}–CePO{sub 4} system at temperatures exceeding 1045 °C. The compound has a eulytite structure with the following structural parameters: a=b=c=10.1655(8) Å, α=β=γ=90.00°, V=1050.46(6) Å{sup 3}. It's melting point exceeds 1950 °C. A limited solid solution exists in the system, which possesses the structure of a low-temperature form of Sr{sub 3}(PO{sub 4}){sub 2}. At 1000 °C the maximal concentration of CePO{sub 4} in the solid solution is below 20 mol%. The solid solution phase field narrows with increased temperature. There is a eutectic point in the (Sr{sub 3}(PO{sub 4}){sub 2}+Sr{sub 3}Ce(PO{sub 4}){sub 3}) phase field at 1765 °C and 15 mol% of CePO{sub 4}. The melting temperature of Sr{sub 3}(PO{sub 4}){sub 2} is 1882±15 °C. - Graphical abstract: The phase diagram of Sr{sub 3}(PO{sub 4}){sub 2}–CePO{sub 4} system showing the stability ranges of limited solid solution and Sr{sub 3}Ce(PO{sub 4}){sub 3} phases. - Highlights: • Sr{sub 3}(PO{sub 4}){sub 2} melts at 1882 °C. • Phase diagram of Sr{sub 3}(PO{sub 4}){sub 2}–CePO{sub 4} system has been proposed. • Limited solid solution of CePO{sub 4} in Sr{sub 3}(PO{sub 4}){sub 2} forms in the system. • The Sr{sub 3}Ce(PO{sub 4}){sub 2} phosphate is stable at temperatures above 1045 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗

Numerical Simulation of Nanostructure Growth

Nanoscale structures, such as nanowires and carbon nanotubes (CNTs), are often grown in gaseous or plasma environments. Successful growth of these structures is defined by achieving a specified crystallinity or chirality, size or diameter, alignment, etc., which in turn depend on gas mixture ratios. pressure, flow rate, substrate temperature, and other operating conditions. To date, there has not been a rigorous growth model that addresses the specific concerns of crystalline nanowire growth, while demonstrating the correct trends of the processing conditions on growth rates. Most crystal growth models are based on the Burton, Cabrera, and Frank (BCF) method, where adatoms are incorporated into a growing crystal at surface steps or spirals. When the supersaturation of the vapor is high, islands nucleate to form steps, and these steps subsequently spread (grow). The overall bulk growth rate is determined by solving for the evolving motion of the steps. Our approach is to use a phase field model to simulate the growth of finite sized nanowire crystals, linking the free energy equation with the diffusion equation of the adatoms. The phase field method solves for an order parameter that defines the evolving steps in a concentration field. This eliminates the need for explicit front tracking/location, or complicated shadowing routines, both of which can be computationally expensive, particularly in higher dimensions. We will present results demonstrating the effect of process conditions, such as substrate temperature, vapor supersaturation, etc. on the evolving morphologies and overall growth rates of the nanostructures.

Hwang, Helen H.↗

Field-induced phase coexistence in an artificial spin ice

Artificial spin-ice systems are magnetic metamaterials consisting of nanomagnet arrays that can be designed to study exotic magnetic states not found in natural materials. Typically, these arrays are modelled as interacting binary macrospins that can only be in an up or down state and are described by the Ising model. These materials have demonstrated ordering transitions, but only via a spontaneous symmetry-breaking mechanism. In this study, we have designed and studied a quadrupole artificial spin-ice system consisting of interacting plaquettes of coupled single-domain nanomagnets that can be interpreted as a composite, ternary variable. After annealing this system in an external magnetic field, we observe both a ferroquadrupolar and an antiferroquadrupolar phase, with an apparent first-order phase boundary and a coexistence regime. The phase diagram of this material is reminiscent of a model used to describe phase coexistence in the superfluid transition of 4 He with 3 He impurities. Finally, these results illustrate how composite magnetic objects realize exotic statistical physics models beyond the Ising model.

36 MATERIALS SCIENCE↗

Three-dimensional in Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions in polycrystalline samples containing several grains with a small misorientation with respect to the temperature gradient. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions that shed new light on the role of subgrain boundaries in the spatiotemporal evolution of the primary cellular spacing and the selection of dendritic array structures.

Kaihua Ji↗

Three-Dimensional In Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions, which shed new light on the selection of dendritic array structures, and on roles of macroscopic curvatures and subgrain boundaries between grains with a small misorientation with respect to the temperature gradient in the spatiotemporal evolution of the primary cellular spacing.

Kaihua Ji↗

Sampling and Reconstruction of the Pupil and Electric Field for Phase Retrieval

This technology is based on sampling considerations for a band-limited function, which has application to optical estimation generally, and to phase retrieval specifically. The analysis begins with the observation that the Fourier transform of an optical aperture function (pupil) can be implemented with minimal aliasing for Q values down to Q = 1. The sampling ratio, Q, is defined as the ratio of the sampling frequency to the band-limited cut-off frequency. The analytical results are given using a 1-d aperture function, and with the electric field defined by the band-limited sinc(x) function. Perfect reconstruction of the Fourier transform (electric field) is derived using the Whittaker-Shannon sampling theorem for 1<Q<2. The Fourier transform is constructed by periodic extension, i.e., by spacing copies of the transform in a definite way, recognizing that no aliasing occurs for values of the sampling ratio such that 1<Q<2, which can be used to advantage in the application of phase retrieval estimation. A method was developed for propagating the electro magnetic field with no aliasing, which has been extended to 2-d optical apertures.

Dean, Bruce↗

Fission gas bubbles and recrystallization-induced degradation of the effective thermal conductivity in U-7Mo fuels

We have developed a mesoscale model to calculate the degradation of the effective thermal conductivity in irradiated U-Mo alloys caused by the fission-induced gas bubbles and recrystallization. The phase-field approach is employed to generate the grain microstructures of U-7Mo fuels with intra- and inter-granular gas bubbles. Based on the phase-field microstructures, the thermal conductivities of U-7Mo as a function of the fission density can be predicted by the developed mesoscale model. The predicted values of effective thermal conductivities are consistent with available experimental data. Results show that the effective thermal conductivity decreases rapidly with recrystallization compared to the one prior to recrystallization, which can be attributed to the sudden increase of grain boundary densities and corresponding intergranular gas bubbles at high fission densities. Smaller grain size fuel structure has a lower thermal conductivity at the same fission density due to the increased grain boundary density. The current study can provide a better understanding of the fission-induced degradation mechanism of the thermal conductivity in U-Mo fuels.

36 MATERIALS SCIENCE↗

Ductile Fracture - Gradient Regularization

Nuclear weapons alteration (ALT) and life extension programs (LEP) are of primary interest to the mission of Sandia National Laboratories. These programs continue to require experimental exploration and computational simulation of ductile failure scenarios to address qualification. Therefore, we invest in generating understanding about ductile failure as demonstrated though experimental procedures and computational simulation of engineering environments. In particular, we study an approach to ductile failure that incorporates the notion of phase-field fracture into our models of inelasticity appropriate for structural alloys. This report covers the formulations of the constitutive model and fracture models used within the phase-field approach and provides some numerical examples highlighting features and the state of the capability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Nonequilibrium Thermodynamics of Hydrate Growth on a Gas-Liquid Interface

We develop a continuum-scale phase-field model to study gas-liquid-hydrate systems far from thermodynamic equilibrium. We design a Gibbs free energy functional for methane-water mixtures that recovers the isobaric temperature-composition phase diagram under thermodynamic equilibrium conditions. The proposed free energy is incorporated into a phase-field model to study the dynamics of hydrate formation on a gas-liquid interface. We elucidate the role of initial aqueous concentration in determining the direction of hydrate growth at the interface, in agreement with experimental observations. Our model also reveals two stages of hydrate growth at an interface—controlled by a crossover in how methane is supplied from the gas and liquid phases—which could explain the persistence of gas conduits in hydratebearing sediments and other nonequilibrium phenomena commonly observed in natural methane hydrate systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows

A mass-conserving lattice Boltzmann method (LBM) for multiphase flows is presented in this paper. The proposed LBM improves a previous model (Lee and Liu, 2010 [21]) in terms of mass conservation, speed-up, and efficiency, and also extends its capabilities for implementation on non-uniform grids. The presented model consists of a phase-field lattice Boltzmann equation (LBE) for tracking the interface between different fluids and a pressure-evolution LBM for recovering the hydrodynamic properties. In addition to the mass conservation property and the simplicity of the algorithm, the advantages of the current phase-field LBE are that it is an order of magnitude faster than the previous interface tracking LBE proposed by Lee and Liu (2010) [21] and it requires less memory resources for data storage. Meanwhile, the pressure-evolution LBM is equipped with a multi-relaxation-time (MRT) collision operator to facilitate attainability of small relaxation rates thereby allowing simulation of multiphase flows at higher Reynolds numbers. Additionally, we reformulate the presented MRT-LBM on nonuniform grids within an adaptive mesh refinement (AMR) framework. Various benchmark studies such as a rising bubble and a falling drop under buoyancy, droplet splashing on a wet surface, and droplet coalescence onto a fluid interface are conducted to examine the accuracy and versatility of the proposed AMR-LBM. The proposed model is further validated by comparing the results with other LB models on uniform grids. A factor of about 20 in savings of computational resources is achieved by using the proposed AMR-LBM. As a more demanding application, the Kelvin–Helmholtz instability (KHI) of a shear-layer flow is investigated for both density-matched and density-stratified binary fluids. The KHI results of the density-matched fluids are shown to be in good agreement with the benchmark AMR results based on the sharp-interface approach. When a density contrast between the two fluids exists, a typical chaotic structure in the flow field is observed at a Reynolds number of 10000, which indicates that the proposed model is a promising tool for direct numerical simulation of two-phase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of MARMOT: A Mesoscale Fuel Performance Code

MARMOT is the mesoscale fuel performance code under development as part of the US DOE Nuclear Energy Advanced Modeling and Simulation Program. In this report, we provide a high level summary of MARMOT, its capabilities, and its current state of validation. The purpose of MARMOT is to predict the coevolution of microstructure and material properties of nuclear fuel and cladding. It accomplished this using the phase field method coupled to solid mechanics and heat conduction. MARMOT is based on the Multiphysics Object-Oriented Simulation Environment (MOOSE), and much of its basic capability in the areas of the phase field method, mechanics, and heat conduction come directly from MOOSE modules. However, additional capability specific to fuel and cladding is available in MARMOT. While some validation of MARMOT has been completed in the areas of fission gas behavior and grain growth, much more validation needs to be conducted. However, new mesoscale data needs to be obtained in order to complete this validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation studies of nucleation of ferroelectric polarization reversal.

Electric field-induced reversal of spontaneous polarization is the defining characteristic of a ferroelectric material, but the process(es) and mechanism(s) associated with the initial nucleation of reverse-polarity domains are poorly understood. This report describes studies carried out using phase field modeling of LiTaO 3 , a relatively simple prototype ferroelectric material, in order to explore the effects of either mechanical deformation or optically-induced free charges on nucleation and resulting domain configuration during field-induced polarization reversal. Conditions were selected to approximate as closely as feasible those of accompanying experimental work in order to provide not only support for the experimental work but also ensure that additional experimental validation of the simulations could be carried out in the future. Phase field simulations strongly support surface mechanical damage/deformation as effective for dramatically reducing the overall coercive field (Ec) via local field enhancements. Further, optically-nucleated polarization reversal appears to occur via stabilization of latent nuclei via the charge screening effects of free charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Grain Growth and Bubble Evolution in U-Mo Alloy by Multiscale Simulations

Increased grain size in U-Mo dispersion fuel is believed to affect the fuel swelling at high fission density. In this work, a multiscale simulation approach combining first-principles calculation and phase-field modeling is used to investigate the grain growth behavior in U-Mo alloys. The material properties of U-Mo alloys predicted by first-principles calculations are incorporated into the mesoscale phase-field models to study the effect of annealing temperature, annealing time and the initial grain structures of fuel particles on the grain growth. The grain growth rate is evaluated and compared with experiment. Meanwhile, the gas bubble evolution kinetics in irradiated U-Mo alloy fuels is investigated to understand its effect on fuel swelling. We systematically examine the effect of Xe, vacancy, and SIA concentration, fission defect generation, and elastic interaction on the growth kinetics of gas bubble. The bubble size distribution and swelling of U-Mo are simulated and compared to experimental measurements.

Mei, Zhi-Gang↗

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗

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↗