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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multi-scale modeling of the electric field assisted sintering process

The electric field assisted sintering (EFAS) process involves tightly coupled physics that influence microstructural evolution in the particles being compacted. It is also an inherently multi-scale phenomenon, with the microstructure of the compact influencing the subsequent engineering-scale response of the sintering system. To improve understanding of how processing parameters influence microstructural evolution, we have developed a multi-scale modeling approach that couples a continuum-level model of the sintering system with a phase-field model for microstructural evolution of particles within the compact. The phase-field model couples the effect of chemical and electrical driving forces on microstructural evolution and includes the effect of charged defect segregation to surfaces and grain boundaries; this segregation leads to enhanced defect transport and heat generation at these interfaces in response to applied electric field. The effect of enhanced heat generation on particle neck growth and the influence of microstructural evolution on the engineering-scale model are demonstrated.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.↗
DOE OSTI2015

Release 1.0 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. With this report, we release a new version of MARMOT, summarize its implemented capabilities, and document the ongoing development of new capabilities planned for the next release. The purpose of MARMOT is to predict the coevolution of microstructure and material properties of nuclear fuel and cladding due to stress, temperature, and irradiation damage. It accomplishes 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.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Interpreting Electrochemical and Chemical Sodiation Mechanisms and Kinetics in Tin Antimony Battery Anodes Using in Situ Transmission Electron Microscopy and Computational Methods

Intermetallic compounds such as SnSb are promising anode materials for sodium ion batteries; however, their nanoscale sodiation mechanisms are not well understood. Here, we used a combination of in situ transmission electron microscopy (TEM), first-principles electronic structure calculations, computational thermodynamic modeling, and phase-field simulations to reveal the sodiation mechanisms and to quantify microstructural effects contributing to the underlying reaction kinetics in SnSb electrodes. During in situ sodiation experiments, the nanocrystalline SnSb thin films underwent a rapid amorphous phase transformation upon sodiation, as determined by in situ TEM and electron diffraction experiments. The Na + diffusion coefficients were measured with and without an external electrical bias, and the data showed that an applied potential increased Na + diffusion by an order of magnitude compared to solid-state diffusion. Furthermore, there was a distinct decrease in sodium diffusion upon the formation of the amorphous phase that resulted from a change in the local structure and grain boundaries. To further understand how the Na + transport mechanism correlated with the changes observed in the SnSb thin films, phase-field modeling was used, which considered sodium diffusion within the grain boundaries together with their evolution and stress–strain state. As a result, these findings enhance our understanding of sodiation mechanisms within intermetallic anode materials for sodium ion battery applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Chiral properties of the zero-field spiral state and field-induced magnetic phases of the itinerant kagome metal YMn 6 Sn 6

Applying a magnetic field in the hexagonal plane of YMn 6 Sn 6 leads to a complex magnetic phase diagram of commensurate and incommensurate phases, one of which coexists with the topological Hall effect (THE) generated by a unique fluctuation-driven mechanism. Using unpolarized neutron diffraction, we report on the solved magnetic structure for two previously identified, but unknown, commensurate phases. These include a low-temperature, high-field fanlike phase and a room-temperature, low-field canted antiferromagnetic phase. Here, an intermediate incommensurate phase between the fanlike and forced ferromagnetic phases is also identified as the last known phase of the in-plane field-temperature diagram. Additional characterization using synchrotron powder diffraction reveals extremely high-quality, single-phase crystals, which suggests that the presence of two incommensurate magnetic structures throughout much of the phase diagram is an intrinsic property of the system. Interestingly, polarized neutron diffraction shows that the centrosymmetric system hosts preferential chirality in the zero-field double-flat-spiral phase, which, along with the THE, is a topologically nontrivial characteristic.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Competing Interface and Bulk Effect–Driven Magnetoelectric Coupling in Vertically Aligned Nanocomposites

Abstract Room‐temperature magnetoelectric (ME) coupling is developed in artificial multilayers and nanocomposites composed of magnetostrictive and electrostrictive materials. While the coupling mechanisms and strengths in multilayers are widely studied, they are largely unexplored in vertically aligned nanocomposites (VANs), even though theory has predicted that VANs exhibit much larger ME coupling coefficients than multilayer structures. Here, strong transverse and longitudinal ME coupling in epitaxial BaTiO 3 :CoFe 2 O 4 VANs measured by both optical second harmonic generation and piezoresponse force microscopy under magnetic fields is reported. Phase field simulations have shown that the ME coupling strength strongly depends on the vertical interfacial area which is ultimately controlled by pillar size. The ME coupling in VANs is determined by the competition between the vertical interface coupling effect and the bulk volume conservation effect. The revealed mechanisms shed light on the physical insights of vertical interface coupling in VANs in general, which can be applied to a variety of nanocomposites with different functionalities beyond the studied ME coupling effect.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Pd₂Cd 11–δ (0.21≤δ≤0.51)–a partly disordered γ-brass type phase and Pd₀.₂₃₈Cd₀.₇₆₂ -a γ-brass related incommensurate phase in the palladium–cadmium system

The Cd rich part of the Cd–Pd phase diagram was reassessed by means of synthesis and single crystal and powder X-ray diffraction. The region contains two phases that have been reported to have substantial compositional widths, Cd₁₁Pd₂ and Cd₄Pd. The phase Cd₁₁Pd₂ that has previously been reported to be a disordered γ-brass crystallizing is space group P4-bar 3m is here shown to crystallize in I4-bar 3m and the mechanism for compositional variation is explained. The phase Pd₄Cd has previously been shown to constitute a phase field or a phase bundle of modulated structures and here we determine the structure of a compound Pd₀.₂₃₈Cd₀.₇₆₂ which crystallizes in the orthorhombic superspace group mm(α00)0s0 (F=[(½, ½, 0, 0); (½, 0, ½, 0 ); (0, ½, ½, 0 )] with the fundamental cell dimensions a=4.687(2) Å, b=10.000(1) Å, c=14.140(2) Å, q=0.6432(6)a⁎. - Graphical abstract: The crystal structures of the partly disordered γ-brass type Pd₂Cd 11–δ (0.21≤δ≤0.51) and γ-brass related Pd₀.₂₃₈Cd₀.₇₆₂ in the palladium–cadmium binary system. Highlights: • Partly disordered γ-brass type Pd₂Cd 11–δ (0.21≤δ≤0.51) and γ-brass related Pd₀.₂₃₈Cd₀.₇₆₂ have been synthesized. • The Pd₂Cd 11–δ structure have been described by cluster concept. • Incommensurately modulated Pd₀.₂₃₈Cd₀.₇₆₂ have been described by (3+1) D space description. • The structure of Pd₀.₂₃₈Cd₀.₇₆₂ has two different chains of atomic subunits, each with their own translational periodicities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Atomic electrostatic maps of 1D channels in 2D semiconductors using 4D scanning transmission electron microscopy

Defects in materials give rise to fluctuations in electrostatic fields that reflect the local charge density, but imaging this with single atom sensitivity is challenging. However, if possible, this provides information about the energetics of adatom binding, localized conduction channels, molecular functionality and their relationship to individual bonds. Here, ultrastable electron-optics are combined with a high-speed 2D electron detector to map electrostatic fields around individual atoms in 2D monolayers using 4D scanning transmission electron microscopy. Simultaneous imaging of the electric field, phase, annular dark field and the total charge in 2D MoS 2 and WS 2 is demonstrated for pristine areas and regions with 1D wires. The in-gap states in sulphur line vacancies cause 1D electron-rich channels that are mapped experimentally and confirmed using density functional theory calculations. We show how electrostatic fields are sensitive in defective areas to changes of atomic bonding and structural determination beyond conventional imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Overview of lower length scale model development for accident tolerant fuels regarding U 3 Si 2 fuel and FeCrAl cladding

U 3 Si 2 and FeCrAl have been proposed as fuel and cladding concepts, respectively, for accident tolerance fuels with higher tolerance to accident scenarios compared to UO 2 . However, a lot of key physics and material properties regarding their in-pile performance are yet to be explored. To accelerate the understanding and reduce the cost of experimental studies, multiscale modeling and simulation are used to develop physics-based materials models to assist engineering scale fuel performance modeling. In this report, the lower-length-scale efforts in method and material model development supported by the Accident Tolerance Fuel (ATF) high-impact-problem (HIP) under the NEAMS program are summarized. Significant progresses have been made regarding interatomic potential, phase field models for phase decomposition and gas bubble formation, and thermal conductivity for U 3 Si 2 fuel, and precipitation in FeCrAl cladding. The accomplishments are very useful by providing atomistic and mesoscale tools, improving the current understanding, and delivering engineering scale models for these two ATF concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1985

Defect-phase-dynamics approach to statistical domain-growth problem of clock models

The growth of statistical domains in quenched Ising-like p-state clock models with p = 3 or more is investigated theoretically, reformulating the analysis of Ohta et al. (1982) in terms of a phase variable and studying the dynamics of defects introduced into the phase field when the phase variable becomes multivalued. The resulting defect/phase domain-growth equation is applied to the interpretation of Monte Carlo simulations in two dimensions (Kaski and Gunton, 1983; Grest and Srolovitz, 1984), and problems encountered in the analysis of related Potts models are discussed. In the two-dimensional case, the problem is essentially that of a purely dissipative Coulomb gas, with a sq rt t growth law complicated by vertex-pinning effects at small t.

Kawasaki, K.↗
DOE OSTI2015

Multi-scale modeling of inter-granular fracture in UO2

A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Microstructure-Level Modeling of Stage 3 Fission Gas Release in UO 2 Fuel

A new phase-field model of fission gas bubble evolution was developed and applied to gain an improved understanding of the microstructure-level processes leading to fission gas release from nuclear fuel and to inform improvements to the fission gas release model used in BISON. The phase-field model simulates the evolution of multiple fuel grains and intergranular bubbles simultaneously, tracks the local concentration of vacancies and gas atoms, and includes the effect of hydrostatic gas pressure on the surrounding fuel matrix. This represents the most realistic microstructure-level model of fission gas bubbles developed to date. The model was used to simulate the growth of grain face and grain edge bubbles in a hexagonally periodic 3D grain structure, which contains multiple grain boundaries and triple junctions. The coverage of grain boundaries and triple junctions by fission gas bubbles was calculated as a function of time. Preliminary results suggest that grain edge bubbles do not release a significant amount of gas prior to gas release by grain face bubbles, which supports the assumptions of the fission gas release model currently used in BISON. However, the simulations should be improved by including a physically realistic number of grain edge bubbles in the initial conditions before definitive conclusions can be drawn.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Draft of M2 Report on Integration of the Hybrid Hydride Model into INL's MBM Framework for Review

This report documents the development, demonstration and validation of a mesoscale, microstructural evolution model for simulation of zirconium hydride {delta}-ZrH{sub 1.5} precipitation in the cladding of used nuclear fuels that may occur during long-term dry storage. While the Zr-based claddings are manufactured free of any hydrogen, they absorb hydrogen during service, in the reactor by a process commonly termed ‘hydrogen pick-up’. The precipitation and growth of zirconium hydrides during dry storage is one of the most likely fuel rod integrity failure mechanisms either by embrittlement or delayed hydride cracking of the cladding. While the phenomenon is well documented and identified as a potential key failure mechanism during long-term dry storage (NUREG/CR-7116), the ability to actually predict the formation of hydrides is poor. The model being documented in this work is a computational capability for the prediction of hydride formation in different claddings of used nuclear fuels. This work supports the Used Fuel Disposition Research and Development Campaign in assessing the structural engineering performance of the cladding during and after long-term dry storage. This document demonstrates a basic hydride precipitation model that is built on a recently developed hybrid Potts-phase field model that combines elements of Potts-Monte Carlo and the phase-field models. The model capabilities are demonstrated along with the incorporation of the starting microstructure, thermodynamics of the Zr-H system and the hydride formation mechanism.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
DOE OSTI2015

Intrinsic space charge layers and field enhancement in ferroelectric nanojunctions

Conducting characteristics of topological defects in ferroelectric materials, such as charged domain walls, engendered a broad interest on their scientific merit and the possibility of novel applications utilizing domain engineering. At the same time, the problem of electron transport in ferroelectrics still remains full of unanswered questions and becomes yet more relevant over the growing interest in ferroelectric semiconductors and new improper ferroelectric materials. We have employed self-consistent phase-field modeling to investigate the physical properties of a local metal-ferroelectric (Pb(Zr{sub 0.2}Ti{sub 0.8})O{sub 3}) junction in applied electric field. We revealed an up to 10-fold local enhancement of electric field realized by large polarization gradient and over-polarization effects due to inherent non-linear dielectric properties of Pb(Zr{sub 0.2}Ti{sub 0.8})O{sub 3}. The effect is independent of bias polarity and maintains its strength prior, during and after ferroelectric switching. The observed field enhancement can be considered on similar grounds as increased doping level, giving rise to reduced switching bias and threshold voltages for charge injection, electrochemical and photoelectrochemical processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Quantum oscillations in strong magnetic fields, berry phase, and superconductivity in three-dimensional topological Bi{sub 2–x}Cu{sub x}Se{sub 3} insulators

Two-dimensional (2D) Shubnikov–de Haas oscillations and 2D Hall oscillations are observed in 3D copper-doped Bi{sub 2}Se{sub 3} single crystals in magnetic fields up to 19.5 T at temperatures down to 0.3 K. Three samples with a high bulk carrier concentration (n ≈ 10{sup 19}–10{sup 20} cm{sup –3}) are studied. The rotation of the samples in a magnetic field shows that these oscillations are related to numerous parallel 2D conducting channels 1–5 nm thick. Their basic kinetic parameters are found. Quantized Hall resistance R{sub xy} is detected in 1-nm-thick 2D conducting channels at high fields. The distance Δ(1/R{sub xy}) between the steps in the field dependence of 1/R{sub xy} is found to be constant for different Landau levels, 1.3e{sup 2}/h per 1-nm-thick layer. The constructed fan diagrams of 2D Landau levels for various angles of sample inclination with respect to the magnetic field direction allowed us to conclude that the Berry phase in the 2D conducting channels is γ ≈ π and independent of the magnetic field direction. When studying the angular dependence of upper resistive critical magnetic field H{sub c2} in one of the superconducting samples, we showed that it can be considered as a bulk superconductor consisting of superconducting layers with an effective thickness of about 50 nm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗