4. 电堆热模型：从一维到二维，把热量算明白
各位工程师朋友，咱们今天聊聊电堆的热模型。说实话，我刚入行那会儿，觉得热管理不就是加个风扇、通个冷却液嘛，有什么难的？直到有一次，一个50kW的液流电池系统在满负荷测试时，电堆内部温度直接飙到了65℃，电压效率掉了将近8个百分点……嗯，从那以后，我再也不敢小看热模型了。
电堆的热行为，说白了就是“产热”和“散热”的博弈。产热算不准，散热就无从谈起。这一节，我带你手把手把电堆的一维和二维热模型建起来，把三个热源项——欧姆热、反应热、过电位热——一个一个算清楚。
核心观点：电堆热模型是热管理设计的“眼睛”。没有它，你就是在黑暗中摸索。
4.1 为什么需要热模型？
你想想看，液流电池电堆里，电解液流过电极表面，电化学反应在发生，电流在传导，离子在迁移——所有这些过程，都在产生热量。热量如果不及时带走，电堆温度就会上升。温度一高，问题就来了：
- 膜的选择性下降，交叉污染加剧
- 电极反应速率变化，副反应增多
- 电解液稳定性下降，可能析出沉淀
我在项目中遇到过最头疼的情况是：电堆进出口温差只有3℃，但内部局部热点温度比平均温度高了12℃。你说，如果只用平均温度来设计冷却系统，那个热点迟早会出问题。所以，热模型不是要不要建的问题，而是怎么建得准的问题。
4.2 一维热模型：快速估算的利器
一维模型，说白了就是沿着电解液流动方向，把电堆切成若干个小段，每一段假设温度均匀。这个模型简单、计算快，适合做初步设计和参数敏感性分析。
4.2.1 模型假设
- 沿流动方向（x方向）存在温度梯度
- 垂直于流动方向（y、z方向）温度均匀
- 电解液物性参数（比热容、密度）为常数
- 忽略辐射换热
4.2.2 控制方程
一维稳态能量平衡方程可以写成：
ρ·Cp·u·(dT/dx) = q_gen / H - (2·h / H)·(T - T_amb)其中：
- ρ — 电解液密度 (kg/m³)
- Cp — 电解液比热容 (J/(kg·K))
- u — 流速 (m/s)
- q_gen — 单位面积产热率 (W/m²)
- H — 流道高度 (m)
- h — 对流传热系数 (W/(m²·K))
- T_amb — 环境温度 (K)
我的经验：一维模型算出来的进出口温差，和实测值偏差通常在15%以内。如果你只需要做方案比选，一维模型完全够用。但如果你要定位热点，那就得上二维甚至三维了。
4.3 二维热模型：看见温度的“横截面”
二维模型比一维多了一个维度——通常是沿着流动方向（x）和垂直于膜的方向（y）。这样就能看到电堆内部的温度分布，包括电极、膜、双极板各自的温度。
4.3.1 模型假设
- 沿x和y方向存在温度梯度
- z方向（宽度方向）温度均匀
- 各层材料（电极、膜、双极板）的热物性已知
- 接触热阻可以等效为界面热阻
4.3.2 控制方程
二维稳态导热方程：
∂/∂x(k_x·∂T/∂x) + ∂/∂y(k_y·∂T/∂y) + q_gen = 0对于不同区域，k_x和k_y取值不同：
- 电极区域：k_x ≈ 0.5-1.5 W/(m·K)，k_y ≈ 0.3-0.8 W/(m·K)
- 膜区域：k ≈ 0.2-0.4 W/(m·K)
- 双极板区域：k ≈ 20-40 W/(m·K)（石墨板）
注意：电极的导热系数是各向异性的！我见过不少同行直接用各向同性假设，结果算出来的热点位置偏差了5mm以上。对于多孔电极，面内和厚度方向的导热系数可能差3-5倍。
4.4 热源项计算：三个“发热大户”
电堆里的热量，主要来自三个地方。我习惯把它们叫做“三兄弟”：欧姆热、反应热、过电位热。下面一个一个拆解。
4.4.1 欧姆热（Joule Heat）
欧姆热是电流通过电阻产生的热量。说白了，就是电子和离子在“挤来挤去”的过程中摩擦生热。
计算公式：
q_ohm = i² · R_eff其中：
- i — 电流密度 (A/m²)
- R_eff — 有效面电阻 (Ω·m²)
R_eff包括：
- 电解液欧姆电阻：R_elyte = δ / σ_elyte
- 膜电阻：R_mem = δ_mem / σ_mem
- 电极电阻：R_elec = δ_elec / σ_elec
- 接触电阻：R_contact（经验值，通常0.01-0.05 Ω·cm²）
避坑指南：我曾经在计算欧姆热时，忘了把接触电阻算进去。结果算出来的总产热量比实测低了8%。后来排查才发现，电堆压紧力不够，接触电阻偏大。从那以后，我每次都会问一句：“接触电阻测了吗？”
4.4.2 反应热（Reaction Heat）
反应热是电化学反应本身的热效应。对于全钒液流电池，正负极反应都是可逆的，反应热等于熵变乘以温度。
计算公式：
q_rxn = (T · ΔS) / (n · F) · i其中：
- ΔS — 反应熵变 (J/(mol·K))
- n — 电子转移数（全钒体系n=1）
- F — 法拉第常数 (96485 C/mol)
- i — 电流密度 (A/m²)
对于全钒液流电池：
- 正极：VO²⁺/VO₂⁺，ΔS ≈ -100 J/(mol·K)（放热）
- 负极：V²⁺/V³⁺，ΔS ≈ +80 J/(mol·K)（吸热）
有意思的是，负极反应是吸热的。所以整个电堆的净反应热，其实是正极放热和负极吸热的代数和。我刚开始做这个的时候，还纳闷为什么算出来的反应热这么小——原来正负极的热效应互相抵消了一部分。
4.4.3 过电位热（Overpotential Heat）
过电位热是“不完美”的代价。实际电位偏离平衡电位的那部分能量，最终都变成了热。
计算公式：
q_act = (η_act,pos + η_act,neg) · i其中η_act由Butler-Volmer方程给出：
i = i₀ · [exp(α·n·F·η_act/(R·T)) - exp(-(1-α)·n·F·η_act/(R·T))]在高电流密度下，可以简化为Tafel形式：
η_act = (R·T)/(α·n·F) · ln(i/i₀)关键数据：在额定工况下（电流密度80-120 mA/cm²），三个热源项的典型占比为：
| 热源项 | 占比范围 | 主要影响因素 |
|---|---|---|
| 欧姆热 | 45-55% | 电解液电导率、膜电阻 |
| 反应热 | 15-25% | 电流密度、温度 |
| 过电位热 | 25-35% | 电极活性、催化剂性能 |
4.5 热模型求解流程
下面这张图，是我做热模型分析时常用的流程。你可以把它当作一个检查清单。
4.6 一个简单的计算示例
假设我们有一个全钒液流电池单电池，参数如下：
| 参数 | 数值 | 单位 |
|---|---|---|
| 电流密度 | 100 | mA/cm² |
| 有效面电阻 | 0.8 | Ω·cm² |
| 温度 | 313 | K |
| 正极过电位 | 0.12 | V |
| 负极过电位 | 0.08 | V |
计算各热源项：
欧姆热：q_ohm = (0.1 A/cm²)² × 0.8 Ω·cm² = 0.008 W/cm² = 80 W/m²
反应热：q_rxn = (313 × 20) / (1 × 96485) × 0.1 = 0.0065 W/cm² = 65 W/m²
（注：正负极净熵变取20 J/(mol·K)）
过电位热：q_act = (0.12 + 0.08) × 0.1 = 0.02 W/cm² = 200 W/m²
总产热：q_total = 80 + 65 + 200 = 345 W/m²你看，在这个工况下，过电位热占了将近58%。这说明电极的催化活性还有提升空间。如果能把过电位降低20%，总产热量就能减少约12%。
我的习惯：每次做完热源项计算，我都会把三个热源的占比画成饼图。一眼就能看出哪个是“发热大户”，然后针对性地优化。比如欧姆热占比高，就考虑换高电导率的电解液或减薄膜；过电位热占比高，就琢磨怎么改进电极。
4.7 模型验证：别让计算骗了你
模型建好了，算完了，就完事了吗？远着呢。我见过太多人模型算得漂漂亮亮，一实测就露馅。验证这一步，绝对不能省。
验证方法：
- 进出口温差验证：在电堆进出口布置热电偶，对比实测温差和模型预测值
- 内部温度验证：在电堆内部插入薄片热电偶（注意密封！），测量几个关键位置的温度
- 热平衡验证：计算总产热量 = 冷却液带走热量 + 表面散热 + 温升蓄热，看是否守恒
血的教训：我曾经在一个项目中，模型算出来的电堆温度比实测低了8℃。排查了三天，最后发现是热电偶的安装位置有问题——它插在了冷却液入口附近，测的是局部低温，不是电堆平均温度。从那以后，我每次布置测温点都会画一张详细的点位图，标注清楚每个测点代表什么。
好了，电堆热模型的核心内容就这些。一维模型帮你快速估算，二维模型让你看清内部温度分布，三个热源项的计算是基础中的基础。把这些搞明白了，热管理设计就有了“眼睛”。
下一节，我们会把这些热模型和冷却系统设计结合起来，看看怎么用模型指导工程实践。
公众号:蓝海资料掘金营，微信deep3321