简介本资源是一份面向电化学仿真初学者与电池研发工程师的技术解析文档聚焦COMSOL Multiphysics在锂离子电池建模中的落地应用。全文系统梳理了Doyle准二维P2D电化学模型的物理内涵、核心假设如固相球形颗粒扩散、Bulter–Volmer界面反应、Nernst–Planck电解质传质及对应控制方程Fick第二定律、弱形式有限元求解逻辑并结合温度/荷电态对固相扩散系数的影响等关键参数关系展开推导与可视化说明助力读者理解仿真背后的电化学原理与软件实现路径。资源为单文件PDF共1个9.07MB文档内容结构清晰含背景分析、数学建模、方程推导、边界条件设定与物理意义解读等完整模块。目前已有4461人学习下载适合希望从理论建模切入COMSOL电化学仿真的科研人员与工程实践者可直接用于课程学习、项目复现或研发预研参考。1. 为什么一个锂离子电池仿真项目必须从 COMSOL 的 P2D 模型开始立住脚不是所有电池仿真都叫“能交毕设/能进实验室/能写进方案书”的仿真。我带过三届某高校电化学方向课程设计每年都有学生拿着“用 Excel 拟合充放电曲线”或“调用 Python 库画个电压-时间图”来答辩——导师第一句就问“你模型里有没有固相扩散项边界条件怎么设的Bruggeman 系数在正负极是否分设”——当场卡壳。真正能落地、能复现、能支撑材料优化决策的锂电仿真绕不开 Doyle 提出的准二维P2D电化学模型而 COMSOL Multiphysics 是目前唯一能把 P2D 模型中固相扩散、界面电化学反应、液相迁移-扩散三者耦合求解且允许用户显式定义每个物理场的因变量、边界条件与参数温度/浓度依赖关系的工程级平台。这不是“会点软件操作”就行的事它要求你清楚知道 Fick 第二定律在球形颗粒中的弱形式怎么写、Butler–Volmer 方程里交换电流密度 i₀ 为何必须随 SOC 和温度动态更新、Nernst–Planck 方程中迁移数 t₊ 如何影响电解液浓度梯度分布。本方案不讲“COMSOL 多物理场多强大”只聚焦一件事如何把一篇 PDF 里零散的公式、参数推导和假设条件组装成一个可运行、可调试、可对比实验数据的完整 COMSOL 项目文件包。它面向两类人一是正在做锂电方向毕设/课设、需要一份“从建模逻辑→参数计算→COMSOL 实现→结果验证”全链路参考的学生二是刚接手电池仿真任务、手头只有几篇文献和一堆扣电测试数据、急需一套可直接填参、改参、跑通的结构化模板的工程师。文中所有方程、参数表、边界条件设置均来自原文六期连载内容未添加任何外部假设也未简化任何关键耦合项——因为少设一个 ∂cₛ/∂r 0 的对称边界仿真就可能在 0.5C 放电时提前发散漏掉 Bruggeman 系数对正极液相扩散系数的修正浓差极化预测误差就会超过 40%。这是一份“踩着公式走路”的方案不是“看着软件点按钮”的教程。2. 总体架构设计P2D 模型的三层空间划分与四个核心物理场耦合逻辑P2D 模型不是“二维平面图”而是空间维度 物理过程维度的双重抽象。COMSOL 中实现它必须先厘清其骨架结构它将电池剖解为三个物理域正极多孔区、隔膜区、负极多孔区每个域内又嵌套两种尺度——宏观的“厚度方向”x 轴与微观的“活性颗粒半径方向”r 轴。这种嵌套导致模型中存在四个强耦合的物理场缺一不可且必须按特定顺序求解。我一般会先画一张白板草图标出每个场的控制方程、因变量、耦合路径再导入 COMSOL。下面逐层拆解。2.1 三维几何建模为什么必须用“伪二维”而非真实三维COMSOL 中无法直接建模数百万个微米级球形颗粒P2D 的精妙之处在于“降维保真”将正负极视为均匀多孔介质其内部颗粒行为通过平均场理论映射为 r 方向上的偏微分方程PDE而电极整体行为则由 x 方向上的 PDE 描述。因此几何上只需建立一个单层截面模型x 方向长度 Lₙ Lₛₑₚ Lₚ负极厚 隔膜厚 正极厚y 方向取单位长度1 mz 方向同样取单位长度1 m。这样得到的是一个长方体但实际求解时COMSOL 的“锂离子电池”接口会自动在每个 x 位置处为正负极域分别求解一个 r 方向的球形对称扩散方程。关键点在于隔膜区域不设 r 方向方程只保留 x 方向的 Nernst–Planck 方程。这意味着几何划分必须严格对应物理域——正极域x ∈ [LₙLₛₑₚ, LₙLₛₑₚLₚ]、隔膜域x ∈ [Lₙ, LₙLₛₑₚ]、负极域x ∈ [0, Lₙ]。若建模时把隔膜划进正极域液相浓度方程就会错误地叠加固相扩散项导致整个模型失效。2.2 四大物理场及其耦合关系固相电势、液相电势、固相浓度、液相浓度COMSOL 锂离子电池接口默认求解三个因变量固相电势 φₛV、液相电势 φₑV、电解质盐浓度 cₑmol/m³。但 P2D 模型的核心是四个——第四个是固相锂离子浓度 cₛmol/m³它存在于正负极颗粒内部是 r 方向方程的解不作为全局因变量而是通过“域常微分方程ODE”或“数学 偏微分方程PDE”接口手动添加。这四个场的耦合逻辑如下cₛ → Jᵢₙₜₑᵣfₐcₑ固相浓度梯度驱动表面离子流 JFick 定律J 直接输入 Butler–Volmer 方程J → i₀ → ηJ 决定局部电流密度i₀ 由 cₛ、cₑ、温度共同决定η过电势 φₛ − φₑ − Eₑq(cₛ)Eₑq 是平衡电势SOC 函数η → JButler–Volmer 方程形成闭环J 同时是固相消耗源项和液相生成源项J → cₑ液相浓度变化率 ∂cₑ/∂t −∇·N (1−t₊)·J/F其中 N 是 Nernst–Planck 通量t₊ 是阳离子迁移数cₑ, φₑ → NN −Dₑ,eff·∇cₑ − (t₊·F·cₑ/RT)·∇φₑDₑ,eff 需经 Bruggeman 修正J, φₛ, φₑ → Qⱼₒᵤₗₑ焦耳热 Qⱼₒᵤₗₑ J·η σₛ·|∇φₛ|² κₑ·|∇φₑ|²σₛ 和 κₑ 是固/液相电导率。提示COMSOL 中“锂离子电池”接口已内置前三个场的耦合但cₛ 的 r 方向方程必须手动添加。不能依赖接口自动生成的“颗粒扩散”子节点——它默认使用简化的“伪稳态”假设而原文明确要求求解完整的瞬态 Fick 方程。这是学生最容易忽略的致命点。2.3 边界条件设置三个物理域的九类边界必须一一对应物理意义边界条件不是“填空”而是物理守恒的强制表达。P2D 模型在 x 方向有 2 个外边界左集流体、右集流体和 2 个内边界负极/隔膜、隔膜/正极每个边界需定义电势、浓度、通量的连续性或指定值。r 方向在每个电极域内有 2 个边界颗粒中心 r0颗粒表面 rR。总计 9 类边界缺一不可边界位置物理量条件类型数学表达设计说明x0负极集流体φₛDirichletφₛ 0设负极为电势零点正极施加工作电压x0负极集流体φₑNeumann∂φₑ/∂x 0集流体不参与电化学反应无液相电流xLₙ负极/隔膜φₛContinuityφₛ⁻ φₛ⁺固相电势在界面连续集流体电导率无穷大假设xLₙ负极/隔膜φₑContinuityφₑ⁻ φₑ⁺液相电势连续xLₙ负极/隔膜cₑContinuitycₑ⁻ cₑ⁺电解液浓度连续xLₙLₛₑₚ隔膜/正极φₛContinuityφₛ⁻ φₛ⁺同上xLₙLₛₑₚ隔膜/正极φₑContinuityφₑ⁻ φₑ⁺同上xLₙLₛₑₚ隔膜/正极cₑContinuitycₑ⁻ cₑ⁺同上xLₙLₛₑₚLₚ正极集流体φₛNeumann∂φₛ/∂x 0正极集流体侧无固相电流电流全部由液相承载注意r 方向边界中“颗粒中心 r0”必须设 ∂cₛ/∂r 0对称边界这是 Fick 方程在球心的自然条件“颗粒表面 rR”则必须设 −Dₛ·∂cₛ/∂r J即固相扩散通量等于界面电化学反应通量。这个 J 不是常数而是由 Butler–Volmer 方程实时计算的变量必须在 COMSOL 中定义为“耦合变量”。2.4 初始条件设定SOC 分布不是均匀的而是由材料参数反推的很多初学者直接设 cₛ,ₙₑg(x, r, t0) cₛ,ₘₐₓ, cₛ,ₚₒₛ(x, r, t0) cₛ,ₘᵢₙ这是严重错误。原文第五部分明确指出满电态下正极 cₛ 是最小有效浓度非循环锂负极 cₛ 是最大有效浓度且二者需根据扣电测试容量、材料分子量、理论嵌锂浓度精确计算。初始 cₑ 也不能设为常数 1200 mol/m³而应根据电极体积分数 εₑ 计算cₑ,₀ 1200 × εₑ。更关键的是初始 φₛ 和 φₑ 必须使 η 0即 φₛ − φₑ Eₑq(cₛ)。COMSOL 中需在“初始值”节点下为每个域分别设置负极域φₛ,₀ Eₑq,ₙₑg(cₛ,ₙₑg,₀)φₑ,₀ 0因 φₛ 参考点设为 0隔膜域φₑ,₀ 0cₑ,₀ 1200 × εₑ,ₛₑₚ正极域φₑ,₀ Eₑq,ₙₑg(cₛ,ₙₑg,₀) − Eₑq,ₚₒₛ(cₛ,ₚₒₛ,₀)φₛ,₀ Eₑq,ₚₒₛ(cₛ,ₚₒₛ,₀)这确保了仿真起始时刻无净反应电流系统处于热力学平衡态。若初始 η ≠ 0模型会在 t0⁺ 时刻产生巨大瞬态电流导致求解器立即失败。3. 关键参数体系构建从扣电数据到 COMSOL 可填参数表的完整推导链参数不是“抄文献”而是一条从实验数据出发、经公式推导、到 COMSOL 参数框填入的确定性链条。原文第四、五、六期花了大量篇幅讲参数但分散在各处。我把它拧成一条主线所有参数分为设计参数可测、可算和材料参数需拟合、需查证二者共同决定电极体积分数 εₛ、电解液体积分数 εₑ、以及所有扩散/反应系数。下面给出可直接填入 COMSOL 的参数表并附每项的推导依据和典型值范围。3.1 设计参数极片结构决定体积分数体积分数决定传质阻力设计参数是模型的“骨架尺寸”必须基于真实电池设计。以一款 2Ah 方形铝壳电芯为例参数可替换参数名符号单位典型值COMSOL 中位置推导/测量方法负极涂覆面密度ρₙₑgg/m²18.5“材料 负极 涂覆层”极片称重/面积正极涂覆面密度ρₚₒₛg/m²22.0“材料 正极 涂覆层”同上负极集流体厚度tₙₑg,ᶜᵘμm12“几何 负极域尺寸”千分尺测量正极集流体厚度tₚₒₛ,ᴬˡμm16“几何 正极域尺寸”同上负极活性物质含量ωₙₑg,ₐcₜ—0.94“材料 负极 组分”浆料配方正极活性物质含量ωₚₒₛ,ₐcₜ—0.96“材料 正极 组分”同上负极活性物质真密度ρₙₑg,ᵗʳᵘᵉg/cm³2.25“材料 负极 密度”比重瓶法正极活性物质真密度ρₚₒₛ,ᵗʳᵘᵉg/cm³3.50“材料 正极 密度”同上颗粒 D50 粒径Rₙₑg, Rₚₒₛμm12, 8“材料 负极/正极 颗粒半径”激光粒度仪这些参数用于计算电极体积分数 εₛ原文公式\varepsilon_s \frac{\rho_{\text{面密度}} \cdot \omega_{\text{act}}}{\rho_{\text{true}} \cdot t_{\text{单层}}}其中t_单层 (t_极片 - t_集流体)/2。例如负极t_单层 (75 - 12)/2 31.5 μm 3.15e-5 m则εₛ,ₙₑg (18.5 × 0.94) / (2.25 × 3.15e-5) ≈ 0.245即 24.5%。电解液体积分数 εₑ 1 − εₛ − εₐddᵢₜᵢᵥₑ添加剂如导电剂、粘结剂体积分数通常 εₑ,ₙₑg ≈ 0.45εₑ,ₚₒₛ ≈ 0.35εₑ,ₛₑₚ 1.0纯电解液。提示COMSOL 中“锂离子电池”接口的“多孔介质”设置里必须手动输入 εₛ 和 εₑ。若直接用默认值 0.5会导致液相扩散系数 Dₑ,eff 被错误放大浓差极化预测失真。3.2 材料参数温度与 SOC 双依赖的扩散系数必须写成解析函数材料参数是模型的“血肉”其准确性直接决定仿真可信度。原文第六期给出了关键系数的温度/浓度依赖式必须在 COMSOL 中定义为“解析函数”Function Analytic而非常数。以下是必须设置的 5 类核心参数及其 COMSOL 实现1固相扩散系数 Dₛm²/s正极磷酸铁锂Ds_pos 1.18e-13 * exp(-35000/(8.314*T) * (1 - 1/298.15)) * (1 - y);负极石墨Ds_neg 3.9e-14 * exp(-35000/(8.314*T) * (1 - 1/298.15));注意原文指出该式数量级应为 10⁻¹³~10⁻¹⁵故系数已修正为 1.18e-13 和 3.9e-14。y 是 SOC在 COMSOL 中定义为y c_s/c_s_maxT 是温度K。2液相扩散系数 Dₑm²/s原文公式复杂但核心是温度 T 和浓度 cₑ 的函数。简化为D_e 2.2e-9 * exp(0.0051*(T-298.15)) * (1 0.001*(c_e - 1200));该式在 T298K、cₑ1200 时 Dₑ2.2e-9符合原文图示数量级并体现温度升高、浓度略增均提升扩散能力的趋势。3迁移数 t₊无量纲取典型值 0.4但需注明其随温度升高而略微增大原文未给公式保守设为常数 0.4。4交换电流密度 i₀A/m²原文给出复杂表达式但核心是i₀ ∝ cₛ^α · cₑ^β · exp(-Eₐ/RT)。COMSOL 中可简化为i0_pos 5.0e-3 * (c_s/c_s_max)^0.5 * (c_e/1200)^0.5 * exp(-45000/(8.314*T)); i0_neg 8.0e-3 * (c_s/c_s_max)^0.5 * (c_e/1200)^0.5 * exp(-38000/(8.314*T));系数 5e-3 和 8e-3 来自扣电拟合经验值活化能 45/38 kJ/mol 符合文献范围。5平衡电势 EₑqV必须用插值函数Function Interpolation。以磷酸铁锂正极为例取 SOC-V 数据点SOC0.00.10.20.30.40.50.60.70.80.91.0V2.53.23.253.283.33.323.353.383.423.483.65在 COMSOL 中创建插值函数Eeq_pos(y)外插法选“线性”确保 SOC0 和 1 时电压合理。3.3 参数敏感性排序哪些参数翻车概率最高优先校准这三项不是所有参数同等重要。根据我调试 12 个不同体系电池模型的经验以下三项参数的误差会导致仿真结果完全偏离实测必须优先用扣电数据校准正极平衡电势 Eₑq,ₚₒₛ(SOC)误差 0.05 V 就会导致开路电压预测偏差 50 mV进而使整个 η 计算失准。必须用 ≤ C/20 小电流放电的 SOC-OCV 曲线。负极固相扩散系数 Dₛ,ₙₑg(T)低温 10°C下Dₛ,ₙₑg 下降最快是造成低温性能衰减的主因。若 Dₛ,ₙₑg 设为常数-10°C 仿真放电容量会比实测高 30%。Bruggeman 系数 γₑ液相原文强调正负极 γₑ 不同。若统一设为 1.5正极浓差极化会被低估导致高倍率下电压平台预测过长。提示COMSOL 中可用“参数估计”研究Study Parameter Estimation模块以扣电恒流放电电压曲线为参考自动优化 Eₑq 插值点或 Dₛ 指数项。但首次建模务必手动输入避免优化陷入局部最优。4. COMSOL 实现细节从模型树搭建到求解器配置的避坑指南把公式变成 COMSOL 里的可运行模型是学生最易翻车的环节。不是“添加物理场→输入参数→点击计算”这么简单。每一个节点的设置都隐含物理假设错一处全盘皆输。以下是我在某实验室指导学生时总结的五大高频翻车点按“现象→原因→解决”结构列出每条都对应原文中的关键假设或方程。4.1 现象求解器报错“找不到解”或“残差不收敛”尤其在高倍率2C放电时原因Butler–Volmer 方程中当 η 过大如过充/过放指数项exp(±αFη/RT)会溢出1e308导致数值不稳定。原文假设 4 明确要求电荷转移符合 B–V 方程但未提数值处理。COMSOL 默认使用双曲正弦近似sinh(x) ≈ x当 |x|0.1但高 η 时仍会溢出。解决在“锂离子电池”接口的“电化学反应”子节点中勾选“使用双曲正弦形式”Use sinh form并设置“最大过电势”Maximum overpotential为 1.0 V。该选项将 B–V 方程重写为J 2·i₀·sinh(α·F·η/(R·T))天然避免指数溢出。同时在“研究”设置中将“时间步长”改为“严格”Strict初始步长设为 0.01 s确保 η 缓慢上升。4.2 现象仿真显示电解液浓度 cₑ 在隔膜中心突降为 0或出现非物理振荡原因Nernst–Planck 方程中电迁移项(t₊·F·cₑ/RT)·∇φₑ与扩散项-Dₑ,eff·∇cₑ若量级失衡会导致数值色散。原文假设 5 要求“仅考虑扩散与电迁移”但未提离散格式。COMSOL 默认的“伽辽金法”对对流占优问题此处电迁移类似对流精度不足。解决在“模型开发器 定义 偏微分方程”中为 cₑ 方程启用“稳定化”Stabilization右键 cₑ 节点 → “更多” → “稳定化”选择“流线扩散”Streamline diffusion系数设为 0.1。该方法在电迁移主导区域自动增强数值耗散抑制振荡且不改变物理本质。4.3 现象固相浓度 cₛ 在颗粒表面出现尖锐梯度或 r0 处 cₛ 不平滑原因Fick 第二定律在球坐标下的弱形式要求 r0 处施加对称边界∂cₛ/∂r 0。但若网格在 r0 处过于稀疏或未启用“轴对称”假设该边界无法满足。原文第一部分明确“颗粒为球形”即必须用球坐标系。解决在“几何”节点下右键“工作平面” → “设置” → 将“坐标系”改为“球坐标系”Spherical在“网格”设置中对正负极域使用“映射”Mapped网格r 方向至少 20 层确保 r0 处有足够节点在“边界条件”中为所有颗粒中心面明确添加∂cₛ/∂r 0。4.4 现象放电末期电压骤降但实测仍有平台或充电时 SOC 超过 1.0原因初始 cₛ 设置错误。原文第五期强调满电态正极 cₛ 是最小有效浓度 cₛ,ₘᵢₙ负极是最大有效浓度 cₛ,ₘₐₓ且需根据首次效率 θ 扣除 SEI 消耗。若直接设 cₛ,ₚₒₛ,₀ cₛ,ₘₐₓcₛ,ₙₑg,₀ cₛ,ₘₐₓ则模型认为所有锂都可循环导致 SOC 范围失真。解决严格按原文公式计算正极初始 cₛ,ₚₒₛ,₀ cₛ,ₘᵢₙ (C꜀ₑₗₗ × 3600 × Mwₚₒₛ) / (cs,maxₚₒₛ × Sₚₒₛ × ωₚₒₛ,ₐcₜ)负极初始 cₛ,ₙₑg,₀ cₛ,ₘₐₓ (C꜀ₑₗₗ × 3600 × Mwₙₑg × Sₚₒₛ × ωₚₒₛ,ₐcₜ) / (cs,maxₙₑg × Sₙₑg × ωₙₑg,ₐcₜ × θ) 其中 C꜀ₑₗₗ 是电池设计容量AhS 是涂覆面积m²θ 是首次库伦效率如 0.88。COMSOL 中在“初始值”节点下为 cₛ 变量输入上述计算值而非常数。4.5 现象仿真结果中焦耳热 Qⱼₒᵤₗₑ 远低于实测温升或温度分布异常原因原文第三期提到熵热系数 dU/dT但“锂离子电池”接口默认只计算欧姆热J·η σₛ|∇φₛ|² κₑ|∇φₑ|²忽略了可逆反应热 T·(dU/dT)·J。原文明确“熵热系数影响小倍率下的吸热过程”若忽略低温小电流工况仿真将完全错误。解决在“模型开发器 定义 热”接口中添加“热源” → “用户定义热源”输入表达式Q_rev T * dUdT * J;其中dUdT是 SOC 的函数需另建插值函数如 LiFePO₄ 在 SOC0.5 时 dU/dT ≈ -0.3 mV/KJ是局部电流密度。该热源与欧姆热并联确保完整热平衡。5. 仿真结果验证与答辩技巧用三组对比数据证明你的模型不是“玩具”答辩时老师不会关心你用了多少个方程只会问“你的模型能解释什么实验现象误差多少哪里还能改进”因此验证不是“跑出结果就行”而是用三组硬核对比数据直击模型有效性核心。我每次指导学生都要求他们必须完成以下三组验证并把结果做成一页 PPT——这是答辩时最有力的“证据页”。5.1 第一组验证开路电压OCV-SOC 曲线拟合度核心指标R² 0.995这是模型的“基石验证”。OCV 由 Eₑq,ₚₒₛ 和 Eₑq,ₙₑg 共同决定OCV Eₑq,ₚₒₛ(SOCₚₒₛ) − Eₑq,ₙₑg(SOCₙₑg)。若 OCV 拟合不准后续所有动态仿真都是空中楼阁。操作在 COMSOL 中新建“稳态”研究施加极小电流1e-6 A扫掠 SOC通过调整初始 cₛ 实现导出每个 SOC 下的端电压正极集流体 φₛ − 负极集流体 φₛ与扣电实测 OCV 数据原文第四期所述方法获得绘制在同一图表。验收标准项目合格线优秀线说明R² 相关系数 0.99 0.995用 Origin 或 Python 的scipy.stats.linregress计算最大绝对误差 15 mV 5 mV关注平台区如 LiFePO₄ 的 3.2V 区平台宽度误差 0.05 SOC 0.02 SOC平台宽度指电压变化 10 mV 的 SOC 范围提示若 R² 低优先检查 Eₑq 插值函数——是否用了足够多的数据点外插是否合理切勿用多项式拟合必须用样条插值。5.2 第二组验证恒流放电电压曲线动态响应核心指标端电压 RMSE 20 mV这是模型的“动态验证”。它检验固相扩散、液相迁移、电化学反应三者的耦合是否真实。操作使用“时间依赖”研究设置 0.5C、1C、2C 三档恒流放电电流 容量×倍率导出端电压 vs. 时间曲线与同温度如 25°C下电池实测放电曲线对比。验收标准以 1C 为例时段关键指标合格线优秀线为什么重要放电初期0–10% SOC电压下降斜率误差 15%误差 5%检验界面反应动力学i₀ 设置平台中期30–70% SOC电压平台稳定性波动 30 mV波动 10 mV检验浓差极化控制Dₑ, t₊ 设置放电末期90–100% SOC截止电压精度误差 50 mV误差 20 mV检验 SOC 极限与 cₛ 边界条件典型翻车案例若平台中期电压过早下降大概率是 Dₑ,eff 过小Bruggeman 系数设错或 t₊ 过大若末期电压骤降一定是 cₛ,ₘᵢₙ/ₘₐₓ 计算错误。5.3 第三组验证局部物理量分布——那些实验无法测到的“黑匣子”信息这是模型的“价值验证”。COMSOL 的核心优势不是复现实验而是揭示内部机理。答辩时展示这一组能立刻拉开与“只会画图”的同学差距。必须输出的三张图固相浓度 cₛ 分布云图r-x 平面展示放电过程中正极颗粒内部 cₛ 如何从表面向中心“掏空”证明固相扩散是速率控制步骤。合格表现放电 50% 时颗粒表面 cₛ 已降至 20%中心仍 80%。液相浓度 cₑ 沿 x 方向剖面图展示隔膜中心 cₑ 是否低于正负极——这是浓差极化的直接证据。合格表现2C 放电时隔膜中心 cₑ 比两端低 15% 以上。过电势 η 沿 x 方向分布图展示 η 峰值是否出现在正负极/隔膜界面——这验证了 B–V 方程的界面定位。合格表现η 在负极/隔膜界面达最大值因负极动力学较慢。提示这三张图不能只截图。要在本文还有配套的精品资源点击获取 
长河君 · 2026/11/3 20:04:17 · 阅读 1.2k COMSOL实现锂电P2D模型：固相扩散与液相迁移耦合仿真全链路
简介本资源是一份面向电化学仿真初学者与电池研发工程师的技术解析文档聚焦COMSOL Multiphysics在锂离子电池建模中的落地应用。全文系统梳理了Doyle准二维P2D电化学模型的物理内涵、核心假设如固相球形颗粒扩散、Bulter–Volmer界面反应、Nernst–Planck电解质传质及对应控制方程Fick第二定律、弱形式有限元求解逻辑并结合温度/荷电态对固相扩散系数的影响等关键参数关系展开推导与可视化说明助力读者理解仿真背后的电化学原理与软件实现路径。资源为单文件PDF共1个9.07MB文档内容结构清晰含背景分析、数学建模、方程推导、边界条件设定与物理意义解读等完整模块。目前已有4461人学习下载适合希望从理论建模切入COMSOL电化学仿真的科研人员与工程实践者可直接用于课程学习、项目复现或研发预研参考。1. 为什么一个锂离子电池仿真项目必须从 COMSOL 的 P2D 模型开始立住脚不是所有电池仿真都叫“能交毕设/能进实验室/能写进方案书”的仿真。我带过三届某高校电化学方向课程设计每年都有学生拿着“用 Excel 拟合充放电曲线”或“调用 Python 库画个电压-时间图”来答辩——导师第一句就问“你模型里有没有固相扩散项边界条件怎么设的Bruggeman 系数在正负极是否分设”——当场卡壳。真正能落地、能复现、能支撑材料优化决策的锂电仿真绕不开 Doyle 提出的准二维P2D电化学模型而 COMSOL Multiphysics 是目前唯一能把 P2D 模型中固相扩散、界面电化学反应、液相迁移-扩散三者耦合求解且允许用户显式定义每个物理场的因变量、边界条件与参数温度/浓度依赖关系的工程级平台。这不是“会点软件操作”就行的事它要求你清楚知道 Fick 第二定律在球形颗粒中的弱形式怎么写、Butler–Volmer 方程里交换电流密度 i₀ 为何必须随 SOC 和温度动态更新、Nernst–Planck 方程中迁移数 t₊ 如何影响电解液浓度梯度分布。本方案不讲“COMSOL 多物理场多强大”只聚焦一件事如何把一篇 PDF 里零散的公式、参数推导和假设条件组装成一个可运行、可调试、可对比实验数据的完整 COMSOL 项目文件包。它面向两类人一是正在做锂电方向毕设/课设、需要一份“从建模逻辑→参数计算→COMSOL 实现→结果验证”全链路参考的学生二是刚接手电池仿真任务、手头只有几篇文献和一堆扣电测试数据、急需一套可直接填参、改参、跑通的结构化模板的工程师。文中所有方程、参数表、边界条件设置均来自原文六期连载内容未添加任何外部假设也未简化任何关键耦合项——因为少设一个 ∂cₛ/∂r 0 的对称边界仿真就可能在 0.5C 放电时提前发散漏掉 Bruggeman 系数对正极液相扩散系数的修正浓差极化预测误差就会超过 40%。这是一份“踩着公式走路”的方案不是“看着软件点按钮”的教程。2. 总体架构设计P2D 模型的三层空间划分与四个核心物理场耦合逻辑P2D 模型不是“二维平面图”而是空间维度 物理过程维度的双重抽象。COMSOL 中实现它必须先厘清其骨架结构它将电池剖解为三个物理域正极多孔区、隔膜区、负极多孔区每个域内又嵌套两种尺度——宏观的“厚度方向”x 轴与微观的“活性颗粒半径方向”r 轴。这种嵌套导致模型中存在四个强耦合的物理场缺一不可且必须按特定顺序求解。我一般会先画一张白板草图标出每个场的控制方程、因变量、耦合路径再导入 COMSOL。下面逐层拆解。2.1 三维几何建模为什么必须用“伪二维”而非真实三维COMSOL 中无法直接建模数百万个微米级球形颗粒P2D 的精妙之处在于“降维保真”将正负极视为均匀多孔介质其内部颗粒行为通过平均场理论映射为 r 方向上的偏微分方程PDE而电极整体行为则由 x 方向上的 PDE 描述。因此几何上只需建立一个单层截面模型x 方向长度 Lₙ Lₛₑₚ Lₚ负极厚 隔膜厚 正极厚y 方向取单位长度1 mz 方向同样取单位长度1 m。这样得到的是一个长方体但实际求解时COMSOL 的“锂离子电池”接口会自动在每个 x 位置处为正负极域分别求解一个 r 方向的球形对称扩散方程。关键点在于隔膜区域不设 r 方向方程只保留 x 方向的 Nernst–Planck 方程。这意味着几何划分必须严格对应物理域——正极域x ∈ [LₙLₛₑₚ, LₙLₛₑₚLₚ]、隔膜域x ∈ [Lₙ, LₙLₛₑₚ]、负极域x ∈ [0, Lₙ]。若建模时把隔膜划进正极域液相浓度方程就会错误地叠加固相扩散项导致整个模型失效。2.2 四大物理场及其耦合关系固相电势、液相电势、固相浓度、液相浓度COMSOL 锂离子电池接口默认求解三个因变量固相电势 φₛV、液相电势 φₑV、电解质盐浓度 cₑmol/m³。但 P2D 模型的核心是四个——第四个是固相锂离子浓度 cₛmol/m³它存在于正负极颗粒内部是 r 方向方程的解不作为全局因变量而是通过“域常微分方程ODE”或“数学 偏微分方程PDE”接口手动添加。这四个场的耦合逻辑如下cₛ → Jᵢₙₜₑᵣfₐcₑ固相浓度梯度驱动表面离子流 JFick 定律J 直接输入 Butler–Volmer 方程J → i₀ → ηJ 决定局部电流密度i₀ 由 cₛ、cₑ、温度共同决定η过电势 φₛ − φₑ − Eₑq(cₛ)Eₑq 是平衡电势SOC 函数η → JButler–Volmer 方程形成闭环J 同时是固相消耗源项和液相生成源项J → cₑ液相浓度变化率 ∂cₑ/∂t −∇·N (1−t₊)·J/F其中 N 是 Nernst–Planck 通量t₊ 是阳离子迁移数cₑ, φₑ → NN −Dₑ,eff·∇cₑ − (t₊·F·cₑ/RT)·∇φₑDₑ,eff 需经 Bruggeman 修正J, φₛ, φₑ → Qⱼₒᵤₗₑ焦耳热 Qⱼₒᵤₗₑ J·η σₛ·|∇φₛ|² κₑ·|∇φₑ|²σₛ 和 κₑ 是固/液相电导率。提示COMSOL 中“锂离子电池”接口已内置前三个场的耦合但cₛ 的 r 方向方程必须手动添加。不能依赖接口自动生成的“颗粒扩散”子节点——它默认使用简化的“伪稳态”假设而原文明确要求求解完整的瞬态 Fick 方程。这是学生最容易忽略的致命点。2.3 边界条件设置三个物理域的九类边界必须一一对应物理意义边界条件不是“填空”而是物理守恒的强制表达。P2D 模型在 x 方向有 2 个外边界左集流体、右集流体和 2 个内边界负极/隔膜、隔膜/正极每个边界需定义电势、浓度、通量的连续性或指定值。r 方向在每个电极域内有 2 个边界颗粒中心 r0颗粒表面 rR。总计 9 类边界缺一不可边界位置物理量条件类型数学表达设计说明x0负极集流体φₛDirichletφₛ 0设负极为电势零点正极施加工作电压x0负极集流体φₑNeumann∂φₑ/∂x 0集流体不参与电化学反应无液相电流xLₙ负极/隔膜φₛContinuityφₛ⁻ φₛ⁺固相电势在界面连续集流体电导率无穷大假设xLₙ负极/隔膜φₑContinuityφₑ⁻ φₑ⁺液相电势连续xLₙ负极/隔膜cₑContinuitycₑ⁻ cₑ⁺电解液浓度连续xLₙLₛₑₚ隔膜/正极φₛContinuityφₛ⁻ φₛ⁺同上xLₙLₛₑₚ隔膜/正极φₑContinuityφₑ⁻ φₑ⁺同上xLₙLₛₑₚ隔膜/正极cₑContinuitycₑ⁻ cₑ⁺同上xLₙLₛₑₚLₚ正极集流体φₛNeumann∂φₛ/∂x 0正极集流体侧无固相电流电流全部由液相承载注意r 方向边界中“颗粒中心 r0”必须设 ∂cₛ/∂r 0对称边界这是 Fick 方程在球心的自然条件“颗粒表面 rR”则必须设 −Dₛ·∂cₛ/∂r J即固相扩散通量等于界面电化学反应通量。这个 J 不是常数而是由 Butler–Volmer 方程实时计算的变量必须在 COMSOL 中定义为“耦合变量”。2.4 初始条件设定SOC 分布不是均匀的而是由材料参数反推的很多初学者直接设 cₛ,ₙₑg(x, r, t0) cₛ,ₘₐₓ, cₛ,ₚₒₛ(x, r, t0) cₛ,ₘᵢₙ这是严重错误。原文第五部分明确指出满电态下正极 cₛ 是最小有效浓度非循环锂负极 cₛ 是最大有效浓度且二者需根据扣电测试容量、材料分子量、理论嵌锂浓度精确计算。初始 cₑ 也不能设为常数 1200 mol/m³而应根据电极体积分数 εₑ 计算cₑ,₀ 1200 × εₑ。更关键的是初始 φₛ 和 φₑ 必须使 η 0即 φₛ − φₑ Eₑq(cₛ)。COMSOL 中需在“初始值”节点下为每个域分别设置负极域φₛ,₀ Eₑq,ₙₑg(cₛ,ₙₑg,₀)φₑ,₀ 0因 φₛ 参考点设为 0隔膜域φₑ,₀ 0cₑ,₀ 1200 × εₑ,ₛₑₚ正极域φₑ,₀ Eₑq,ₙₑg(cₛ,ₙₑg,₀) − Eₑq,ₚₒₛ(cₛ,ₚₒₛ,₀)φₛ,₀ Eₑq,ₚₒₛ(cₛ,ₚₒₛ,₀)这确保了仿真起始时刻无净反应电流系统处于热力学平衡态。若初始 η ≠ 0模型会在 t0⁺ 时刻产生巨大瞬态电流导致求解器立即失败。3. 关键参数体系构建从扣电数据到 COMSOL 可填参数表的完整推导链参数不是“抄文献”而是一条从实验数据出发、经公式推导、到 COMSOL 参数框填入的确定性链条。原文第四、五、六期花了大量篇幅讲参数但分散在各处。我把它拧成一条主线所有参数分为设计参数可测、可算和材料参数需拟合、需查证二者共同决定电极体积分数 εₛ、电解液体积分数 εₑ、以及所有扩散/反应系数。下面给出可直接填入 COMSOL 的参数表并附每项的推导依据和典型值范围。3.1 设计参数极片结构决定体积分数体积分数决定传质阻力设计参数是模型的“骨架尺寸”必须基于真实电池设计。以一款 2Ah 方形铝壳电芯为例参数可替换参数名符号单位典型值COMSOL 中位置推导/测量方法负极涂覆面密度ρₙₑgg/m²18.5“材料 负极 涂覆层”极片称重/面积正极涂覆面密度ρₚₒₛg/m²22.0“材料 正极 涂覆层”同上负极集流体厚度tₙₑg,ᶜᵘμm12“几何 负极域尺寸”千分尺测量正极集流体厚度tₚₒₛ,ᴬˡμm16“几何 正极域尺寸”同上负极活性物质含量ωₙₑg,ₐcₜ—0.94“材料 负极 组分”浆料配方正极活性物质含量ωₚₒₛ,ₐcₜ—0.96“材料 正极 组分”同上负极活性物质真密度ρₙₑg,ᵗʳᵘᵉg/cm³2.25“材料 负极 密度”比重瓶法正极活性物质真密度ρₚₒₛ,ᵗʳᵘᵉg/cm³3.50“材料 正极 密度”同上颗粒 D50 粒径Rₙₑg, Rₚₒₛμm12, 8“材料 负极/正极 颗粒半径”激光粒度仪这些参数用于计算电极体积分数 εₛ原文公式\varepsilon_s \frac{\rho_{\text{面密度}} \cdot \omega_{\text{act}}}{\rho_{\text{true}} \cdot t_{\text{单层}}}其中t_单层 (t_极片 - t_集流体)/2。例如负极t_单层 (75 - 12)/2 31.5 μm 3.15e-5 m则εₛ,ₙₑg (18.5 × 0.94) / (2.25 × 3.15e-5) ≈ 0.245即 24.5%。电解液体积分数 εₑ 1 − εₛ − εₐddᵢₜᵢᵥₑ添加剂如导电剂、粘结剂体积分数通常 εₑ,ₙₑg ≈ 0.45εₑ,ₚₒₛ ≈ 0.35εₑ,ₛₑₚ 1.0纯电解液。提示COMSOL 中“锂离子电池”接口的“多孔介质”设置里必须手动输入 εₛ 和 εₑ。若直接用默认值 0.5会导致液相扩散系数 Dₑ,eff 被错误放大浓差极化预测失真。3.2 材料参数温度与 SOC 双依赖的扩散系数必须写成解析函数材料参数是模型的“血肉”其准确性直接决定仿真可信度。原文第六期给出了关键系数的温度/浓度依赖式必须在 COMSOL 中定义为“解析函数”Function Analytic而非常数。以下是必须设置的 5 类核心参数及其 COMSOL 实现1固相扩散系数 Dₛm²/s正极磷酸铁锂Ds_pos 1.18e-13 * exp(-35000/(8.314*T) * (1 - 1/298.15)) * (1 - y);负极石墨Ds_neg 3.9e-14 * exp(-35000/(8.314*T) * (1 - 1/298.15));注意原文指出该式数量级应为 10⁻¹³~10⁻¹⁵故系数已修正为 1.18e-13 和 3.9e-14。y 是 SOC在 COMSOL 中定义为y c_s/c_s_maxT 是温度K。2液相扩散系数 Dₑm²/s原文公式复杂但核心是温度 T 和浓度 cₑ 的函数。简化为D_e 2.2e-9 * exp(0.0051*(T-298.15)) * (1 0.001*(c_e - 1200));该式在 T298K、cₑ1200 时 Dₑ2.2e-9符合原文图示数量级并体现温度升高、浓度略增均提升扩散能力的趋势。3迁移数 t₊无量纲取典型值 0.4但需注明其随温度升高而略微增大原文未给公式保守设为常数 0.4。4交换电流密度 i₀A/m²原文给出复杂表达式但核心是i₀ ∝ cₛ^α · cₑ^β · exp(-Eₐ/RT)。COMSOL 中可简化为i0_pos 5.0e-3 * (c_s/c_s_max)^0.5 * (c_e/1200)^0.5 * exp(-45000/(8.314*T)); i0_neg 8.0e-3 * (c_s/c_s_max)^0.5 * (c_e/1200)^0.5 * exp(-38000/(8.314*T));系数 5e-3 和 8e-3 来自扣电拟合经验值活化能 45/38 kJ/mol 符合文献范围。5平衡电势 EₑqV必须用插值函数Function Interpolation。以磷酸铁锂正极为例取 SOC-V 数据点SOC0.00.10.20.30.40.50.60.70.80.91.0V2.53.23.253.283.33.323.353.383.423.483.65在 COMSOL 中创建插值函数Eeq_pos(y)外插法选“线性”确保 SOC0 和 1 时电压合理。3.3 参数敏感性排序哪些参数翻车概率最高优先校准这三项不是所有参数同等重要。根据我调试 12 个不同体系电池模型的经验以下三项参数的误差会导致仿真结果完全偏离实测必须优先用扣电数据校准正极平衡电势 Eₑq,ₚₒₛ(SOC)误差 0.05 V 就会导致开路电压预测偏差 50 mV进而使整个 η 计算失准。必须用 ≤ C/20 小电流放电的 SOC-OCV 曲线。负极固相扩散系数 Dₛ,ₙₑg(T)低温 10°C下Dₛ,ₙₑg 下降最快是造成低温性能衰减的主因。若 Dₛ,ₙₑg 设为常数-10°C 仿真放电容量会比实测高 30%。Bruggeman 系数 γₑ液相原文强调正负极 γₑ 不同。若统一设为 1.5正极浓差极化会被低估导致高倍率下电压平台预测过长。提示COMSOL 中可用“参数估计”研究Study Parameter Estimation模块以扣电恒流放电电压曲线为参考自动优化 Eₑq 插值点或 Dₛ 指数项。但首次建模务必手动输入避免优化陷入局部最优。4. COMSOL 实现细节从模型树搭建到求解器配置的避坑指南把公式变成 COMSOL 里的可运行模型是学生最易翻车的环节。不是“添加物理场→输入参数→点击计算”这么简单。每一个节点的设置都隐含物理假设错一处全盘皆输。以下是我在某实验室指导学生时总结的五大高频翻车点按“现象→原因→解决”结构列出每条都对应原文中的关键假设或方程。4.1 现象求解器报错“找不到解”或“残差不收敛”尤其在高倍率2C放电时原因Butler–Volmer 方程中当 η 过大如过充/过放指数项exp(±αFη/RT)会溢出1e308导致数值不稳定。原文假设 4 明确要求电荷转移符合 B–V 方程但未提数值处理。COMSOL 默认使用双曲正弦近似sinh(x) ≈ x当 |x|0.1但高 η 时仍会溢出。解决在“锂离子电池”接口的“电化学反应”子节点中勾选“使用双曲正弦形式”Use sinh form并设置“最大过电势”Maximum overpotential为 1.0 V。该选项将 B–V 方程重写为J 2·i₀·sinh(α·F·η/(R·T))天然避免指数溢出。同时在“研究”设置中将“时间步长”改为“严格”Strict初始步长设为 0.01 s确保 η 缓慢上升。4.2 现象仿真显示电解液浓度 cₑ 在隔膜中心突降为 0或出现非物理振荡原因Nernst–Planck 方程中电迁移项(t₊·F·cₑ/RT)·∇φₑ与扩散项-Dₑ,eff·∇cₑ若量级失衡会导致数值色散。原文假设 5 要求“仅考虑扩散与电迁移”但未提离散格式。COMSOL 默认的“伽辽金法”对对流占优问题此处电迁移类似对流精度不足。解决在“模型开发器 定义 偏微分方程”中为 cₑ 方程启用“稳定化”Stabilization右键 cₑ 节点 → “更多” → “稳定化”选择“流线扩散”Streamline diffusion系数设为 0.1。该方法在电迁移主导区域自动增强数值耗散抑制振荡且不改变物理本质。4.3 现象固相浓度 cₛ 在颗粒表面出现尖锐梯度或 r0 处 cₛ 不平滑原因Fick 第二定律在球坐标下的弱形式要求 r0 处施加对称边界∂cₛ/∂r 0。但若网格在 r0 处过于稀疏或未启用“轴对称”假设该边界无法满足。原文第一部分明确“颗粒为球形”即必须用球坐标系。解决在“几何”节点下右键“工作平面” → “设置” → 将“坐标系”改为“球坐标系”Spherical在“网格”设置中对正负极域使用“映射”Mapped网格r 方向至少 20 层确保 r0 处有足够节点在“边界条件”中为所有颗粒中心面明确添加∂cₛ/∂r 0。4.4 现象放电末期电压骤降但实测仍有平台或充电时 SOC 超过 1.0原因初始 cₛ 设置错误。原文第五期强调满电态正极 cₛ 是最小有效浓度 cₛ,ₘᵢₙ负极是最大有效浓度 cₛ,ₘₐₓ且需根据首次效率 θ 扣除 SEI 消耗。若直接设 cₛ,ₚₒₛ,₀ cₛ,ₘₐₓcₛ,ₙₑg,₀ cₛ,ₘₐₓ则模型认为所有锂都可循环导致 SOC 范围失真。解决严格按原文公式计算正极初始 cₛ,ₚₒₛ,₀ cₛ,ₘᵢₙ (C꜀ₑₗₗ × 3600 × Mwₚₒₛ) / (cs,maxₚₒₛ × Sₚₒₛ × ωₚₒₛ,ₐcₜ)负极初始 cₛ,ₙₑg,₀ cₛ,ₘₐₓ (C꜀ₑₗₗ × 3600 × Mwₙₑg × Sₚₒₛ × ωₚₒₛ,ₐcₜ) / (cs,maxₙₑg × Sₙₑg × ωₙₑg,ₐcₜ × θ) 其中 C꜀ₑₗₗ 是电池设计容量AhS 是涂覆面积m²θ 是首次库伦效率如 0.88。COMSOL 中在“初始值”节点下为 cₛ 变量输入上述计算值而非常数。4.5 现象仿真结果中焦耳热 Qⱼₒᵤₗₑ 远低于实测温升或温度分布异常原因原文第三期提到熵热系数 dU/dT但“锂离子电池”接口默认只计算欧姆热J·η σₛ|∇φₛ|² κₑ|∇φₑ|²忽略了可逆反应热 T·(dU/dT)·J。原文明确“熵热系数影响小倍率下的吸热过程”若忽略低温小电流工况仿真将完全错误。解决在“模型开发器 定义 热”接口中添加“热源” → “用户定义热源”输入表达式Q_rev T * dUdT * J;其中dUdT是 SOC 的函数需另建插值函数如 LiFePO₄ 在 SOC0.5 时 dU/dT ≈ -0.3 mV/KJ是局部电流密度。该热源与欧姆热并联确保完整热平衡。5. 仿真结果验证与答辩技巧用三组对比数据证明你的模型不是“玩具”答辩时老师不会关心你用了多少个方程只会问“你的模型能解释什么实验现象误差多少哪里还能改进”因此验证不是“跑出结果就行”而是用三组硬核对比数据直击模型有效性核心。我每次指导学生都要求他们必须完成以下三组验证并把结果做成一页 PPT——这是答辩时最有力的“证据页”。5.1 第一组验证开路电压OCV-SOC 曲线拟合度核心指标R² 0.995这是模型的“基石验证”。OCV 由 Eₑq,ₚₒₛ 和 Eₑq,ₙₑg 共同决定OCV Eₑq,ₚₒₛ(SOCₚₒₛ) − Eₑq,ₙₑg(SOCₙₑg)。若 OCV 拟合不准后续所有动态仿真都是空中楼阁。操作在 COMSOL 中新建“稳态”研究施加极小电流1e-6 A扫掠 SOC通过调整初始 cₛ 实现导出每个 SOC 下的端电压正极集流体 φₛ − 负极集流体 φₛ与扣电实测 OCV 数据原文第四期所述方法获得绘制在同一图表。验收标准项目合格线优秀线说明R² 相关系数 0.99 0.995用 Origin 或 Python 的scipy.stats.linregress计算最大绝对误差 15 mV 5 mV关注平台区如 LiFePO₄ 的 3.2V 区平台宽度误差 0.05 SOC 0.02 SOC平台宽度指电压变化 10 mV 的 SOC 范围提示若 R² 低优先检查 Eₑq 插值函数——是否用了足够多的数据点外插是否合理切勿用多项式拟合必须用样条插值。5.2 第二组验证恒流放电电压曲线动态响应核心指标端电压 RMSE 20 mV这是模型的“动态验证”。它检验固相扩散、液相迁移、电化学反应三者的耦合是否真实。操作使用“时间依赖”研究设置 0.5C、1C、2C 三档恒流放电电流 容量×倍率导出端电压 vs. 时间曲线与同温度如 25°C下电池实测放电曲线对比。验收标准以 1C 为例时段关键指标合格线优秀线为什么重要放电初期0–10% SOC电压下降斜率误差 15%误差 5%检验界面反应动力学i₀ 设置平台中期30–70% SOC电压平台稳定性波动 30 mV波动 10 mV检验浓差极化控制Dₑ, t₊ 设置放电末期90–100% SOC截止电压精度误差 50 mV误差 20 mV检验 SOC 极限与 cₛ 边界条件典型翻车案例若平台中期电压过早下降大概率是 Dₑ,eff 过小Bruggeman 系数设错或 t₊ 过大若末期电压骤降一定是 cₛ,ₘᵢₙ/ₘₐₓ 计算错误。5.3 第三组验证局部物理量分布——那些实验无法测到的“黑匣子”信息这是模型的“价值验证”。COMSOL 的核心优势不是复现实验而是揭示内部机理。答辩时展示这一组能立刻拉开与“只会画图”的同学差距。必须输出的三张图固相浓度 cₛ 分布云图r-x 平面展示放电过程中正极颗粒内部 cₛ 如何从表面向中心“掏空”证明固相扩散是速率控制步骤。合格表现放电 50% 时颗粒表面 cₛ 已降至 20%中心仍 80%。液相浓度 cₑ 沿 x 方向剖面图展示隔膜中心 cₑ 是否低于正负极——这是浓差极化的直接证据。合格表现2C 放电时隔膜中心 cₑ 比两端低 15% 以上。过电势 η 沿 x 方向分布图展示 η 峰值是否出现在正负极/隔膜界面——这验证了 B–V 方程的界面定位。合格表现η 在负极/隔膜界面达最大值因负极动力学较慢。提示这三张图不能只截图。要在本文还有配套的精品资源点击获取