Adaptive finite element methods are a powerful tool to obtain numerical simulation results in a reasonable time. Due to complex chemical and mechanical couplings in lithium-ion batteries, numerical simulations are very helpful to investigate promising new battery active materials such as amorphous silicon featuring a higher energy density than graphite. Based on a thermodynamically consistent continuum model with large deformation and chemo-mechanically coupled approach, we compare three different spatial adaptive refinement strategies: Kelly-, gradient recovery- and residual based error estimation. For the residual based case, the strong formulation of the residual is explicitly derived. With amorphous silicon as example material, we investigate two 3D representative host particle geometries, reduced with symmetry assumptions to a 1D unit interval and a 2D elliptical domain. Our numerical studies show that the Kelly estimator overestimates the error, whereas the gradient recovery estimator leads to lower refinement levels and a good capture of the change of the lithium flux. The residual based error estimator reveals a strong dependency on the cell error part which can be improved by a more suitable choice of constants to be more efficient. In a 2D domain, the concentration has a larger influence on the mesh distribution than the Cauchy stress.
翻译:自适应有限元方法是一种在合理时间内获得数值模拟结果的强大工具。由于锂离子电池中复杂的化学与力学耦合作用,数值模拟对于研究有前景的新型电池活性材料(如具有比石墨更高能量密度的非晶硅)非常有帮助。基于一个包含大变形和化学-力学耦合方法的热力学一致连续介质模型,我们比较了三种不同的空间自适应细化策略:Kelly估计法、梯度恢复估计法和基于残差的误差估计法。对于基于残差的情况,我们明确推导了残差的强形式。以非晶硅作为示例材料,我们研究了两种具有代表性的3D宿主颗粒几何形状,通过对称性假设将其简化为1D单位区间和2D椭圆域。我们的数值研究表明,Kelly估计法高估了误差,而梯度恢复估计法导致较低的细化水平,并能较好地捕捉锂通量的变化。基于残差的误差估计器显示出对单元误差部分的强烈依赖性,这可以通过更合适的常数选择来提高效率。在2D域中,浓度对网格分布的影响大于柯西应力。