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 667 records · Page 37

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Microstructural Pattern Formation during Far-from-Equilibrium Alloy Solidification

We introduce a new phase-field formulation of rapid alloy solidification that quantitatively in- corporates non-equilibrium effects at the solid-liquid interface over a very wide range of interface velocities. Simulations identify a new dynamical instability of dendrite tip growth driven by solute trapping at velocities approaching the absolute stability limit. They also reproduce the formation of the widely observed banded microstructures, revealing how this instability triggers transitions between dendritic and microsegregation-free solidification. Predicted band spacings agree quantita- tively with observations in rapidly solidified Al-Cu thin films.

36 MATERIALS SCIENCE↗

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

Design and Evaluation of Environmental Barrier Coatings for Protection of Ceramic Matrix Composites in Hydrogen-Based Turbines

This presentation at the 2024 FECM Spring Project Review Meeting summarizes ongoing work at NETL aimed at enabling the use of ceramic matrix composite (CMC) materials in next-generation high-efficiency power generation turbines fueled by hydrogen, natural gas, or other fuels. The work broadly involves three aspects (1) phase-field damage modeling of CMCs and the environmental barrier coatings (EBCs) used to protect them (2) computational design of new and improved EBCs and (3) performance testing of EBCs using a novel gas turbine combustion simulation rig developed at NETL. An overview on current status and next steps for each aspect of the project are presented.

ceramic matrix composite (CMC)↗