首页/文章/ 详情

如何准确地测量活性材料的固相扩散系数?

10月前浏览552
 

在电池仿真领域,P2D模型是我们模拟电池性能的强大工具。而在模型的一堆动力学参数中,固相扩散系数 (Ds) 无疑是核心之一。它直接决定了锂离子在活性材料颗粒内部的迁移速度,这几乎就定义了电池的倍率性能。试想一下,如果这个如此重要的参数,从源头输入时就是错的,那么你跑出来的仿真结果还能信吗 ?为了避免“垃圾进,垃圾出”的窘境,我们首先要搞清楚,这个系数该怎么测,以及我们测到的到底是什么?

测量固相扩散系数的“三板斧”

在实验室里,大家测量固相扩散系数常用的无非就是“三板斧”:GITT、PITT 和 EIS。

1. GITT (恒电流间歇滴定法)

像挤牙膏一样,先用一小段恒流充/放电,然后让电池“歇一会儿”弛豫,再重复这个过程 。通过分析弛豫过程中电压的变化,我们就可以计算出固相扩散系数。

2. PITT (恒电位间歇滴定法)

这个方法和GITT有些像,但它控制的是电压的阶跃,然后观察电流如何随时间衰减 。通过分析电流衰减的快慢,便可以反推出固相扩散系数(Ds)。

3. EIS (电化学阻抗谱)

通过给电池施加一个微小的正弦波扰动,来观察体系的响应。通过分析其在低频区的韦伯阻抗(Warburg Impedance)的斜率,也同样可以得到固相扩散系数。

表观固相扩散系数

从理论上看,上述这些方法都设计精巧,公式完备,似乎可以准确地得到我们想要的扩散系数。但问题真的这么简单吗?问题的关键在于,我们测量的对象是“极片”,而非纯粹的活性材料本身。一个真实的极片(图1),是由活性材料、导电剂和粘结剂混合而成的复合体,其内部的孔隙中还填充着电解液。

图1. 正极极片典型微观结构示意图

因此,当你对整个极片进行电化学测试时,你所记录到的信号,必然是一个包含了所有组分和过程影响的“综合信号”。这意味着,你得到的固相扩散系数,并不是材料“本征”的扩散系数,而是一个“表观”的固相扩散系数。 这个“表观值”包含了除固相扩散本身以外的多种过程的影响,例如 :

  • 电荷转移动力学
  • SEI/CEI 界面膜的阻抗
  • 液相离子传输
  • 固相电子传输等

因此,极片中的粘结剂、导电剂,以及极片的孔隙率、厚度等物理参数,都不可避免地会影响最终的测量结果。

“表观”与“本征”的差距

你可能会问,这些额外因素造成的影响到底有多大?我们来看一个最经典的例子——磷酸铁锂(LFP)。LFP是一种典型的橄榄石结构材料,其锂离子被认为是准一维扩散(主要沿b轴方向)。在早期的研究中,许多团队采用上述的常规方法,在极片上测得的LFP扩散系数普遍在1013至1016 cm2/s 这个数量级。这个值非常低,以至于当时学界普遍认为,LFP倍率性能的瓶颈就是缓慢的固相扩散。 然而,后来有研究者通过更精密的实验方法,例如使用微电极技术(图2)直接测量单颗粒的扩散行为,从而排除了导电剂、粘结剂、粒径分布等一系列的干扰因素 。结果你们猜怎么着?他们测出的LFP本征固相扩散系数,高达 108至109 cm2/s。

图2. 单颗粒微电极装置示意图

看到了吗?从极片测出的表观值到单颗粒测出的本征值,差距足足有5个数量级!甚至有学者发现,LFP单颗粒能在900C超高倍率下放电依然能保持其全容量的50%(图3)。这些结果是颠覆性的,它们告诉我们,对于性能优化的LFP材料而言,真正限制其倍率性能的,可能并非固相扩散本身,而是电荷转移或其它界面过程。如果你拿着从极片上测出的“表观”数据去指导材料改性,或是直接带入仿真模型,那么努力的方向可能从一开始就跑偏了。

图3. LFP单颗粒不同倍率下的放电曲线(750 nA相当于900 C)

结尾

读到这里,相信大家已经明白,我们日常从极片中测量到的固相扩散系数,本质上是一个受多种因素影响的“综合动力学指标”。但这并不意味着它没有价值,关键在于我们如何看待和使用它。正如文中所说,想要获得更接近材料本征的扩散系数,可以采用单颗粒测试等更为前沿的方法。 那么问题来了,想和各位同行与读者们深入探讨一下:

  • 除了进行更复杂的单颗粒测试,在常规的极片测试中,我们还能通过哪些实验设计或数据分析的技巧,来尽可能地排除干扰,让结果更具参考价值?
  • 对于构建靠谱的P2D模型,是必须追求“本征扩散系数”么?

欢迎在评论区留下您的真知灼见,分享您的经验与思考。如果大家对单颗粒测试技术感兴趣的人数较多,我们后续会考虑单独出一期内容,专门深入聊聊这个话题。

参考文献:

1. Journal of Power Sources 217 (2012) 444e448

2. J. Electrochem. 2015, 21(6): 566-571


来源:锂电芯动
化学电子理论材料控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-10-19
最近编辑:10月前
锂电芯动
博士 中科院博士,电芯仿真高级工程师
获赞 9粉丝 21文章 71课程 0
点赞
收藏
作者推荐

Matlab与Comsol联合仿真(三):多参数自动拟合

我相信很多仿真工程师都有过“建模一小时,调参一整天”的经历 。当模型建立后,参数拟合,尤其是涉及三个以上的多参数拟合,是一项非常不容易且耗时的工作。如果能借助软件自动实现多参数拟合,那将大大节省仿真工程师的宝贵时间,让他们能把时间放在更有意义的事情上面。COMSOL优化模块的局限性虽然 COMSOL 软件自带了优化模块,但个人感觉其参数拟合能力较弱 。这主要是因为其内置的优化算法多为常规的局部搜索方法(比如SNOPT、BOBYQA和Levenberg–Marquardt等),缺乏如遗传算法、粒子群算法、贝叶斯算法等更灵活的全局优化算法 。因此,当需要拟合的参数较多时,COMSOL 自带的优化算法很容易陷入局部最优状态 。我的实际感受常常是,即使花费了很长时间进行优化,参数的改动也不明显。因此,有必要将 COMSOL 与 MATLAB 联合使用,借助 MATLAB 强大的全局优化功能来弥补 COMSOL 在这方面的不足。MATLAB与COMSOL联合调参优势1. 调参效率显著提升传统的复杂参数调优过程需要工程师手动修改模型参数、运行仿真、导出数据、分析结果,然后重复这一循环,过程繁琐且效率低下。而利用MATLAB强大的脚本编程能力,用户可以轻松编写循环或批处理脚本,系统性地、自动化地改变COMSOL模型中的任意参数(如几何尺寸、材料属性、边界条件等),让计算机在夜间或周末自动执行,从而最大化利用计算资源,缩短研发周期。2. 可以融合高级优化算法工程师可以调用MATLAB的Optimization Toolbox™和Global Optimization Toolbox™,利用其中丰富的优化器,如遗传算法(ga)、粒子群算法(particleswarm)、模式搜索(patternsearch)以及各种梯度下降算法(fmincon)来解决COMSOL中的优化问题。这使得寻找非线性、非凸或多峰值复杂问题(在工程领域很常见)的全局最优解成为可能。3. 强大的数据后处理与可视化能力通过LiveLink™接口,仿真结果可以直接无缝地导入MATLAB工作空间。用户可以利用MATLAB强大的矩阵运算能力、信号处理工具箱、统计分析工具箱等对数据进行深入分析,例如执行傅里叶变换、数据拟合、相关性分析等。联合调参案例演示本次案例采用的是一个一维的 P2D(Pseudo-2-Dimension)模型,任务是拟合一条实验测的放电曲线,需要调整七个动力学参数,最终实现仿真的放电曲线与实验测的重合得很好。联合调参采用的是贝叶斯优化算法。选择该算法的原因是它具有极高寻优效率,能够在较少的仿真次数下找到最优的解,非常适合在计算资源受限或仿真周期较短的情况下使用。案例的关键结果展示如下:1. 寻优前的模拟与实验曲线 2. 待寻优的7个动力学参数 3. 最优参数及寻优时间 4. 最优参数的验证效果 补充说明本文展示的多参数自动拟合方法不仅仅局限于锂电池模型,只要在这个案例的代码基础上稍作修改,就可以适用于任何其他模型。为了方便大家根据自己的需求修改我的代码,我录制了一个视频详细地讲解了每一行代码的功能及注意事项。来源:锂电芯动

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈