Microstructure Modeling of Transport and Kinetics in Solid State Batteries with LLZO Solid Electrolyte Fabricated Using the Freeze-Tape-Casting Method
Authors: Jordan Skylar McConnell , Yi Eongyu , Shen Hao , Doeff Marca M. , Roberts Scott A. , Nelson George J. , Shah Krishna 
Published: August 2024

In traditional Li-ion batteries with liquid electrolyte, the transport of lithium ions in the electrolyte and lithium in the active material within a porous electrode are studied extensively. The continuum-scale modeling of these transport phenomena considers a homogenized simulation domain consisting of liquid electrolyte and solid active material with effective transport properties accounting for the heterogeneity [1]. There have been fewer studies that consider digitized versions of the experimentally imaged microstructure in continuum-scale modeling studies, particularly to study mechanics [2]. Despite of the computational challenges, there are inherent advantages in using fully resolved microstructure modeling to improve fundamental understanding of the relevant physical and chemical phenomena in porous electrodes, such as localized degradation and stress concentration. In the case of solid electrolytes, specifically scaffold-type LLZO solid electrolytes, microstructure modeling is of even more relevance, particularly to study anisotropic transport, localized degradation phenomena, stress and deformation leading to loss of contact, and chemo-mechanics. Thus far, microstructure modeling studies on experimentally obtained microstructure consisting of solid electrolyte and active material have been quite limited [3]. The studies that have explored this realm have employed voxel-based mesh to discretize the computation domain [3]. While these methods are sufficient proof of concept, the lack of precision in defining the interface separating active material and electrolyte when using voxel-based mesh provides a source of potentially significant numerical error. Without a highly resolved interface mesh, the total interfacial area increases artificially which leads to inaccuracies in calculations of local reaction rate and overpotential. A highly resolved, smooth interface is also necessary when simulating mechanical deformation and stress field in the microstructure. In the present investigation, we model mass transport and charge transport in the active material/electrolyte microstructure to simulate concentration and potential fields with varying mesh types in order to compare their effectiveness and characterize the error that can be attributed to the mesh type chosen, particularly in the case of voxel-based meshes. Additionally, we study the effect of microstructure domain length chosen in the two directions perpendicular to the shortest length (electrode thickness) as well as boundary condition formulation on the boundaries in these axes on the simulation results. The microstructures studied in this work are derived from the image stacks of the bi/trilayer LLZO solid electrolyte fabricated using freeze-tape-casting method [4]. Overall, the present work aims to develop a better understanding of the best practices in microstructure modeling, particularly when using microstructure image stacks obtained from fabricated samples.

References

[1] M. Doyle, T.F. Fuller, J. Newman, J. Electrochem. Soc., 140, 1526 (1993).
[2] H. Mendoza, S.A. Roberts, V.E. Brunini, A.M. Grillet, Electrochim. Acta, 190, 1 (2016).
[3] M. Alabdali, F.M. Zanotto, V. Viallet, V. Seznec, A.A. Franco, Curr Opin Electrochem, 36, 101127 (2022).
[4] Hao Shen, Eongyu Yi, Stephen Heywood, Dilworth Y. Parkinson, Guoying Chen, Nobumichi Tamura, Stephen Sofie, Kai Chen, and Marca M. Doeff, ACS Appl Mater & Interfaces, 12, 3494 (2020).