Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field modeling”

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

Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains

Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards understanding particle rigid-body motion during solid-state sintering

A quantitative understanding of particle rigid body (RB) motion that inherently accompanies grain boundary (GB) diffusion is highly desirable to understand and control the dynamic interplay between coarsening and densification during solid state sintering. In this work, by computer simulation using a multi-phase-field approach, we analyze systematically the roles played by each of these processes at different stages of the shrinkage of the internal pore in a three-particle green body as a function of particle size as well as thermodynamic and kinetic factors of interfaces. We demonstrate that particle RB translation promotes both neck growth, and pore rounding and shrinkage. Moreover, the forces acting at GBs and pulling neighboring particles towards one another dynamically evolve as particles fuse. In contrast, particle RB rotation has no contribution to pore shrinkage. The translational force acting on an individual particle varies with not only its size, but also the number and sizes of its neighboring particles.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring the relationship between solvent-assisted ball milling, particle size, and sintering temperature in garnet-type solid electrolytes

Garnet-type solid electrolytes, such as Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 (LLZTO), are promising materials for solid-state batteries, but processing remains a challenge, in part due to the high sintering temperature required for densification. Additionally, this temperature can be lowered by decreasing the initial particle size via solvent-assisted ball milling, but the relationship between solvent choice, particle properties, sintering behavior, and ionic conductivity is not well understood. In this work, we systematically explore these parameters, showing that milling in commonly used protic solvents, such as alcohols, effectively decreases the particle size but results in lithium loss (through Li + /H + exchange) that leads to poor sintering. By contrast, milling in aprotic solvents with surfactant reduces the particle size to ~220 nm without lithium loss, enabling the fabrication of dense samples (5.1 g/cm 3 ) with good ionic conductivity (0.43 mS/cm at 25 °C) at a lower sintering temperature (1000 °C). We compare ionic conductivities and activation energies for samples prepared with different particle sizes and sintering temperatures and use multiphase-field simulations to identify the mass transport and microstructural mechanisms responsible for the observed sintering dependence on particle size. These results further clarify the relationship between processing parameters and performance and represent important progress toward overcoming fabrication challenges for these materials.

25 ENERGY STORAGE↗
DOE OSTI2019

Radiation induced nanovoid shrinkage in Cu at room temperature: An in situ study

Radiation induced void swelling is widely observed in a variety of metallic materials. Here, by using sputtering deposition technique, we have introduced faceted nanovoids into Cu films. In-situ Kr ion irradiation was subsequently performed at room temperature to investigate the evolution of nanovoids. Most nanovoids found to shrink gradually with increasing dose. Irradiation induced high-density vacancy clusters exist mostly in the formof stacking fault tetrahedrons. Phase fieldmodeling reveales that void shrinkage arises from biased absorption of interstitials. These findings provide insights to the physical mechanisms of radiation response of nanovoids in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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

Dynamic Phase Engineering of Bendable Transition Metal Dichalcogenide Monolayers

Current interest in two-dimensional (2D) materials is driven in part by the ability to dramatically alter their optoelectronic properties through strain and phase engineering. A combination of these approaches can be applied in quasi-2D transition metal dichalcogenide (TMD) monolayers to induce displacive structural transformations between semiconducting (H) and metallic/semimetallic (T') phases. We classify such transformations in Group VI TMDs, and formulate a multiscale, first-principles-informed modeling framework to describe evolution of microstructural domain morphologies in elastically bendable 2D monolayers. We demonstrate that morphology and mechanical response can be controlled via application of strain either uniformly or through local probes to generate functionally patterned conductive T' domains. Such systems form dynamically programmable electromechanical 2D materials, capable of rapid local switching between domains with qualitatively different transport properties. Finally, this enables dynamic “drawing” of localized conducting regions in an otherwise semiconducting TMD monolayer, opening several interesting device-relevant functionalities such as the ability to dynamically “rewire” a device in real time.

2D materials↗
DOE OSTI2024

Fracture Dynamics in Silicon Anode Solid-State Batteries

Solid-state batteries (SSBs) with silicon anodes could enable improved safety and energy density compared to lithium-ion batteries. However, degradation arising from the massive volumetric changes of silicon anodes during cycling is not well understood in solid-state systems. Here, we use operando X-ray computed microtomography to reveal micro- to macro-scale chemo-mechanical degradation processes of silicon anodes in SSBs. Mud-type channel cracks driven by biaxial tensile stress form across the electrode during delithiation. We also find detrimental cracks at the silicon/ solid electrolyte interface that form due to local reaction competition between neighboring domains of different sizes. Continuum phase-field damage modeling quantifies stress-driven channel cracking and shows that the lithiated silicon stress state is critical for determining the extent of interfacial fracture. This work reveals the mechanisms that govern SSBs compared to conventional lithium-ion batteries and provides guidelines for engineering chemomechanically resilient electrodes for high-energy batteries.

25 ENERGY STORAGE↗
DOE OSTI2017

Flexible Multiferroic Bulk Heterojunction with Giant Magnetoelectric Coupling via van der Waals Epitaxy

Magnetoelectric nanocomposites have been a topic of intense research due to their profound potential in the applications of electronic devices based on spintronic technology. Nevertheless, in spite of significant progress made in the growth of high-quality nanocomposite thin films, the substrate clamping effect still remains a major hurdle in realizing the ultimate magnetoelectric coupling. To overcome this obstacle, an alternative strategy of fabricating a self-assembled ferroelectric–ferrimagnetic bulk heterojunction on a flexible muscovite via van der Waals epitaxy is adopted. In this paper, we investigated the magnetoelectric coupling in a self-assembled BiFeO 3 (BFO)–CoFe 2 O 4 (CFO) bulk heterojunction epitaxially grown on a flexible muscovite substrate. The obtained heterojunction is composed of vertically aligned multiferroic BFO nanopillars embedded in a ferrimagnetic CFO matrix. Moreover, due to the weak interaction between the flexible substrate and bulk heterojunction, the interface is incoherent and, hence, the substrate clamping effect is greatly reduced. The phase-field simulation model also complements our results. The magnetic and electrical characterizations highlight the improvement in magnetoelectric coupling of the BFO–CFO bulk heterojunction. A magnetoelectric coupling coefficient of 74 mV/cm·Oe of this bulk heterojunction is larger than the magnetoelectric coefficient reported earlier on flexible substrates. Finally and therefore, this study delivers a viable route of fabricating a remarkable magnetoelectric heterojunction and yet flexible electronic devices that are robust against extreme conditions with optimized performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Grain boundary formation through particle detachment during coarsening of nanoporous metals

Significance Nanoporous metals are prototypical bicontinuous structures that have a large surface area and many applications. We demonstrate that an in-grain orientation spread develops during coarsening in nanoporous metals, leading to the formation of a nanocrystalline–nanoporous structure. The process of particles detaching and reattaching as the structure coarsens is identified as a mechanism of grain boundary formation. Coarsening the structure results in significant particle detachment with detachment rates that vary drastically with changes in the volume fraction and volume fraction inhomogeneities in the nanostructure. These results demonstrate that particle detachment must be understood to have full control of the nanocrystalline–nanoporous structure and improve material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Interaction Dynamics Between Ferroelectric and Antiferroelectric Domains in a Pb Zr O 3 -Based Ceramic

The antiferroelectric-ferroelectric phase transition in PbZrO 3 -based oxides is of both fundamental and practical importance. In ceramics in which such a transition readily occurs, the antiferroelectric and the ferroelectric phases often coexist in individual grains with a coherent interphase interface. In this work, the electric biasing in situ transmission electron microscopy technique is employed to directly observe a unique microstructural dynamic when ferroelectric and antiferroelectric domains are driven by a moderate electric field to interact. It is found that, under monotonic loading, the ferroelectric domain grows until it is blocked by the ferroelectric-antiferroelectric interface. At the same time, a kink is formed on the interface at the contact point. Here, the interaction of the growing domain with the interface is interpreted in terms of depolarization field-assisted phase transition, which is supported by our phase-field simulation. Upon further bipolar cycling, the ferroelectric domain becomes less mobile and no longer reaches the ferroelectric-antiferroelectric interface, indicative of electric fatigue of the ferroelectric phase.

36 MATERIALS SCIENCE↗
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 OSTI2013

Structural short-range forces between solid-melt interfaces

We predict the structural interaction of crystalline solid-melt interfaces using amplitude equations which are derived from classical density functional theory or phase field crystal modeling. The solid ordering decays exponentially on the scale of the interface thickness at solid-melt interfaces; the overlap of two such profiles leads to a short-range interaction, which is mainly carried by the longest-range density waves, which can facilitate grain boundary premelting. Here, we calculate the tail of these interactions, depending on the relative translation of the two crystals, fully analytically and predict the interaction potential, and compare it to numerical simulations. For grain boundaries the interaction is predicted to decay twice as fast as for two crystals without misorientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Elastically mediated interactions between grain boundaries and precipitates in two-phase coherent solids

In this work, we investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B 81, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. 114, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed nonlinear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Cell dynamics simulations of coupled charge and magnetic phase transformation in correlated oxides

In this work, we present a comprehensive numerical study on the kinetics of phase transition that is characterized by two nonconserved scalar order parameters coupled by a special linear-quadratic interaction. This particular Ginzburg-Landau theory has been proposed to describe the coupled charge and magnetic transition in nickelates and the collinear stripe phase in cuprates. The inhomogeneous state of such systems at low temperatures consists of magnetic domains separated by quasimetallic domain walls where the charge order is reduced. By performing large-scale cell dynamics simulations, we find a two-stage phase-ordering process in which a short period of independent evolution of the two order parameters is followed by a correlated coarsening process. The long-time growth and coarsening of magnetic domains is shown to follow the Allen-Cahn power law. We further show that the nucleation-and-growth dynamics during phase transformation to the ordered states is well described by the Kolmogorov-Johnson-Mehl-Avrami theory in two dimensions. On the other hand, the presence of quasimetallic magnetic domain walls in the ordered states gives rise to a very different kinetics for transformation to the high-temperature paramagnetic phase. In this scenario, the phase transformation is initiated by the decay of magnetic domain walls into two insulator-metal boundaries, which subsequently move away from each other. Implications of our findings to recent nano-imaging experiments on nickelates are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗