4. 电化学-热耦合建模基础
温升仿真的第一步不是打开软件画网格，而是决定“热从哪里来、电化学行为用什么描述”。在储能电池热管理中，常用的耦合路线有三条：以Bernardi产热模型为核心的“热源解析式”路线、以简化电化学模型（SPM/SPMe）为核心的“机理降阶”路线，以及以等效电路模型（ECM）为核心的“外特性拟合”路线。三者都能给出产热功率，但对参数、工况和精度的要求差异很大，选错模型往往导致温升趋势对、峰值错，或者计算量远超必要。
一、Bernardi产热模型：把产热拆成不可逆热与可逆熵热
Bernardi模型是锂离子电池产热计算的经典框架，其体积产热率常写成：
q = I·(U_ocv − U) / V + I·T·(∂U_ocv/∂T) / V其中第一项是不可逆焦耳热与极化热，第二项是可逆熵变热。符号约定需先明确：本文取充电电流为正、放电电流为负，因此放电时第一项为正（产热），充电时若极化很小，第二项可能为负（吸热）。实际仿真中，Bernardi模型通常不单独求解电化学场，而是把 U_ocv、U、∂U_ocv/∂T 作为SOC和温度的函数输入，再与三维热传导方程耦合。
它的适用边界很清晰：适合单体或模组级别的温升趋势分析、风冷/液冷方案对比、以及需要快速迭代的工况筛选。参数需求包括：开路电压曲线 U_ocv(SOC,T)、熵热系数曲线 ∂U_ocv/∂T(SOC)、直流内阻或极化电压 U_ocv−U 随电流和SOC的变化、以及电池体积或产热面积。若只给一个常数内阻，模型在低SOC大倍率下会明显低估产热。
二、简化电化学模型：SPM与SPMe的降阶逻辑
简化电化学模型从P2D模型降阶而来。单粒子模型（SPM）把正负极各视为一个球形活性颗粒，忽略电解液浓度梯度，只保留固相扩散和Butler-Volmer反应动力学。SPMe在此基础上增加电解液浓度和电位修正，能更好描述高倍率下的浓差极化。其核心方程包括固相扩散方程、反应动力学方程和电荷守恒方程，输出的是固相表面浓度、过电位和反应速率，再代入产热表达式得到空间分布更细的热源。
SPM/SPMe的适用边界是：中低倍率、需要看电极内部锂浓度和电位分布、以及希望产热沿厚度方向非均匀的场合。参数需求比Bernardi多得多：正负极颗粒半径、固相扩散系数、最大固相浓度、初始化学计量比、交换电流密度或反应速率常数、电解液电导率与扩散系数（SPMe）、以及电极厚度和面积。这些参数通常来自半电池测试或文献拟合，若扩散系数偏差一个数量级，高倍率下的端电压和产热都会明显偏移。
三、等效电路模型：用RC网络复现外特性
等效电路模型用开路电压源串联欧姆内阻和若干RC并联支路来模拟电池端电压。一阶RC模型包含 R0、R1、C1，二阶RC再增加一组。产热计算通常写成：
q = I²·R0 + I·U1 + I·U2其中 U1、U2 是RC支路端电压，代表极化热。ECM的优点是参数少、辨识容易、实时性好，适合BMS在线估算和系统级热仿真。它的适用边界是：不需要电极内部信息、只关心端电压和总产热、且工况在参数辨识范围内。参数需求包括：U_ocv(SOC,T)、R0(SOC,T,I)、R1、C1 等，通常通过HPPC脉冲测试辨识。若把 R0 固定为常数，低温大倍率下的产热会被严重低估。
四、三条路线的选择与参数准备
实际项目中，模型选择取决于三个问题：是否需要电极内部信息、工况倍率范围、以及可用参数。若只做模组风冷方案对比，Bernardi模型配合SOC相关的内阻曲线通常足够；若研究析锂风险或电极厚度方向温差，SPM/SPMe更合适；若要在BMS中实时估算产热并驱动热管理策略，ECM是首选。
参数准备上，Bernardi模型至少需要一张 U_ocv(SOC) 曲线和一张 R0(SOC,T) 或极化电压表。ECM需要HPPC测试数据，在10%、20%……90% SOC点分别做脉冲，用电压响应拟合 R0、R1、C1。SPM/SPMe则需要半电池或三电极测试获取正负极平衡电位和扩散系数。若参数不足，宁可降阶使用Bernardi或ECM，也不要用未辨识的SPM参数强行计算。
q ≈ I²R0/V，若 R0=1.2 mΩ、I=50 A、V=0.5 L，则焦耳热约 6 W/L；再叠加熵热项，总产热约 7~8 W/L。该数值仅为示例，用于说明量级，不代表实测结果。五、耦合方式与边界条件
产热模型给出 q(x,y,z,t) 后，需与热模型耦合。常见做法是：电化学/ECM在时间步内计算产热功率，作为热传导方程的内热源；热模型更新温度后，再把温度反馈给产热模型中的 R0(T)、U_ocv(T) 或扩散系数，形成双向耦合。若温度对参数影响很小，也可单向耦合以节省计算。边界条件需明确：对流换热系数 h、环境温度 T_amb、以及模组间的接触热阻。这些参数与产热模型同样重要，忽略接触热阻会导致模组内单体温差被低估。
R0 设为常数后直接用于低温大倍率仿真。低温下电解液电导率下降，R0 可能增大2~3倍，产热被严重低估，温升曲线会明显偏乐观。建议至少建立 R0(SOC,T) 二维查表。六、参数需求速查
| 模型 | 核心参数 | 典型获取方式 | 适用倍率 | 输出信息 |
|---|---|---|---|---|
| Bernardi | U_ocv(SOC,T)、∂U_ocv/∂T、极化电压、V | 小倍率充放电+熵热测试 | ≤2C 趋势分析 | 总产热功率 |
| SPM/SPMe | R_s、D_s、c_max、x0/y0、i0、κ、D_e | 半电池/三电极/文献拟合 | ≤3C 机理分析 | 固相浓度、过电位、空间产热 |
| ECM | U_ocv(SOC,T)、R0(SOC,T,I)、R1、C1 | HPPC脉冲测试 | 全倍率外特性 | 端电压、总产热 |
总结来看，Bernardi模型胜在参数少、物理意义清楚，适合快速温升趋势；SPM/SPMe胜在机理完整，适合电极级分析但参数门槛高；ECM胜在实时性和工程可落地，适合BMS与系统级仿真。实际建模时，可先用Bernardi或ECM确定产热量级和温升范围，再在关键工况下用SPM/SPMe校核产热分布，形成“粗筛—精算”的组合流程。