Grain boundary diffusion in polycrystalline materials is a physical phenomenon of great fundamental interest and practical significance. Although accelerated atomic transport along grain boundaries has been known for decades, atomic-level understanding of diffusion mechanisms remains poor. Previous atomistic simulations focused on low temperatures where the grain boundary structure is ordered or high temperatures where it is highly disordered. Here, we conduct molecular dynamics simulations of grain boundary diffusion at intermediate temperatures most relevant to applications. A surprising result of this work is the observation of intermittent GB diffusion behavior and its strong system-size dependence unseen in previous work. Both effects are found to originate from thermally activated point-defect avalanches. We identify the length and time scales of the avalanches and link their formation to dynamic heterogeneity in partially disordered systems. Our findings have implications for future computer modeling of grain boundary diffusion and mass transport in nano-scale materials.
Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.
more » « less- NSF-PAR ID:
- 10512757
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Modelling and Simulation in Materials Science and Engineering
- Volume:
- 32
- Issue:
- 5
- ISSN:
- 0965-0393
- Format(s):
- Medium: X Size: Article No. 055027
- Size(s):
- Article No. 055027
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Grain boundaries critically limit the electronic performance of oxide perovskites. These interfaces lower the carrier mobilities of polycrystalline materials by several orders of magnitude compared to single crystals. Despite extensive effort, improving the mobility of polycrystalline materials (to meet the performance of single crystals) is still a severe challenge. In this work, the grain boundary effect is eliminated in perovskite strontium titanate (STO) by incorporating graphene into the polycrystalline microstructure. An effective mass model provides strong evidence that polycrystalline graphene/strontium titanate (G/STO) nanocomposites approach single crystal‐like charge transport. This phenomenological model reduces the complexity of analyzing charge transport properties so that a quantitative comparison can be made between the nanocomposites and STO single crystals. In other related works, graphene composites also optimize the thermal transport properties of thermoelectric materials. Therefore, decorating grain boundaries with graphene appears to be a robust strategy to achieve “phonon glass–electron crystal” behavior in oxide perovskites.
-
more » « less
The complexity of shear-induced grain boundary dynamics has been historically difficult to view at the atomic scale. Meanwhile, two-dimensional (2D) colloidal crystals have gained prominence as model systems to easily explore grain boundary dynamics at single-particle resolution but have fallen short at exploring these dynamics under shear. Here, we demonstrate how an inherent interfacial shear in 2D colloidal crystals drives microstructural evolution. By assembling paramagnetic particles into polycrystalline sheets using a rotating magnetic field, we generate a particle circulation at the interface of particle-free voids. This circulation shears the crystalline bulk, operating as both a source and sink for grain boundaries. Furthermore, we show that the Read-Shockley theory for hard-condensed matter predicts the misorientation angle and energy of shear-induced low-angle grain boundaries based on their regular defect spacing. Model systems containing shear provide an ideal platform to elucidate shear-induced grain boundary dynamics for use in engineering improved/advanced materials.
-
Efficient heat dissipation in batteries is important for thermal management against thermal runaway and chemical instability at elevated temperatures. Nevertheless, thermal transport processes in battery materials have not been well understood especially considering their complicated microstructures. In this study, lattice thermal transport in lithium cobalt oxide (LiCoO 2 ), a popular cathode material for lithium ion batteries, is investigated via molecular dynamics-based approaches and thermal resistance models. A LiCoO 2 single-crystal is shown to have thermal conductivities in the order of 100 W m −1 K −1 with strong anisotropy, temperature dependence, and size effects. By comparison, polycrystalline LiCoO 2 is more isotropic with much lower thermal conductivities. This difference is caused by random grain orientations, the thermal resistance of grain boundaries, and size-dependent intra-grain thermal conductivities that are unique to polycrystals. The grain boundary thermal conductance is calculated to be in the range of 7.16–25.21 GW m −2 K −1 . The size effects of the intra-grain thermal conductivities are described by two empirical equations. Considering all of these effects, two thermal resistance models are developed to predict the thermal conductivity of polycrystalline LiCoO 2 . The two models predict a consistent thermal conductivity–grain size relationship that agrees well with molecular dynamics simulation results. The insights revealed by this study may facilitate future efforts on battery materials design for improved thermal management.more » « less
-
more » « less
Interactions between an evolving solid and inviscid flow can result in substantial computational complexity, particularly in circumstances involving varied boundary conditions between the solid and fluid phases. Examples of such interactions include melting, sublimation, and deflagration, all of which exhibit bidirectional coupling, mass/heat transfer, and topological change of the solid–fluid interface. The diffuse interface method is a powerful technique that has been used to describe a wide range of solid-phase interface-driven phenomena. The implicit treatment of the interface eliminates the need for cumbersome interface tracking, and advances in adaptive mesh refinement have provided a way to sufficiently resolve diffuse interfaces without excessive computational cost. However, the general scale-invariant coupling of these techniques to flow solvers has been relatively unexplored. In this work, a robust method is presented for treating diffuse solid–fluid interfaces with arbitrary boundary conditions. Source terms defined over the diffuse region mimic boundary conditions at the solid–fluid interface, and it is demonstrated that the diffuse length scale has no adverse effects. To show the efficacy of the method, a one-dimensional implementation is introduced and tested for three types of boundaries: mass flux through the boundary, a moving boundary, and passive interaction of the boundary with an incident acoustic wave. Two-dimensional results are presented as well these demonstrate expected behavior in all cases. Convergence analysis is also performed and compared against the sharp-interface solution, and linear convergence is observed. This method lays the groundwork for the extension to viscous flow and the solution of problems involving time-varying mass-flux boundaries.
