comsol燃料电池模型 等温和不等温的均有 下面是验证图以及参考文献。【燃料电池建模手札】最近在COMSOL折腾质子交换膜燃料电池模型发现等温模型和不等温模型就像冰火两重天。先说个真实案例某次用等温模型预测电流密度分布美如画结果实测数据直接打脸——原来催化层局部过热导致性能骤降这才想起温度场才是隐藏BOSS。先看等温模型怎么玩适合摸鱼版建模def proton_transport(): sigma_ion 1.3*exp(0.013*(T-303)) #离子电导率温度修正项 i_ion -sigma_ion * grad(phi_ion) source_term vol_reaction_current * (1e4*active_surface_area) return div(i_ion) source_term这段的关键在volreactioncurrent这个反应电流源项实测发现用Butler-Volmer方程时别直接用教科书参数得根据实际MEA厚度做归一化处理。有个坑当电压低于0.4V时记得打开氢气渗透修正项否则极化曲线尾巴会翘起来。不等温模型才是真刀真枪的战场看这个能量方程配置% 伪代码表示多物理场耦合 energy_eq rho*Cp*dT/dt - div(k_eff*gradT) Q_rev Q_irrev; Q_rev T*dE_rev/dT * i_local; //可逆热项 Q_irrev (eta_act eta_conc)*i_local; //不可逆热重点在k_eff的有效导热系数实测GDL各向异性导热系数设置不当会让温度场出现诡异条纹。有个骚操作把冷却流道参数化扫描配合批处理计算自动寻找最优流速比手动调参效率高十倍。验证环节最刺激的是对比文献[1]的V-I曲线此处应有但没给的对比图。当电流密度突破1.5A/cm²时等温模型预测电压0.25V实际考虑温度骤升后直接掉到0.18V这和J.P.Koza的实测数据误差3%。网格划分有个魔鬼细节催化层至少分3层单元不然活化损失计算会出锯齿状异常。comsol燃料电池模型 等温和不等温的均有 下面是验证图以及参考文献。最后来个压轴技巧在求解器配置里把电荷守恒和动量方程分步求解收敛速度比全耦合快得像开了外挂。碰到发散别慌把电解质相变潜热作为阻尼项临时加入收敛后再逐步移除这招救过我好几次项目进度。参考文献[1] Wang Y, Chen KS. J.Power Sources 2013 (验证曲线原始数据来源)[2] COMSOL官方燃料电池案例库v5.6网格划分基准
COMSOL燃料电池模型：等温和不等温的验证图及参考文献
comsol燃料电池模型 等温和不等温的均有 下面是验证图以及参考文献。【燃料电池建模手札】最近在COMSOL折腾质子交换膜燃料电池模型发现等温模型和不等温模型就像冰火两重天。先说个真实案例某次用等温模型预测电流密度分布美如画结果实测数据直接打脸——原来催化层局部过热导致性能骤降这才想起温度场才是隐藏BOSS。先看等温模型怎么玩适合摸鱼版建模def proton_transport(): sigma_ion 1.3*exp(0.013*(T-303)) #离子电导率温度修正项 i_ion -sigma_ion * grad(phi_ion) source_term vol_reaction_current * (1e4*active_surface_area) return div(i_ion) source_term这段的关键在volreactioncurrent这个反应电流源项实测发现用Butler-Volmer方程时别直接用教科书参数得根据实际MEA厚度做归一化处理。有个坑当电压低于0.4V时记得打开氢气渗透修正项否则极化曲线尾巴会翘起来。不等温模型才是真刀真枪的战场看这个能量方程配置% 伪代码表示多物理场耦合 energy_eq rho*Cp*dT/dt - div(k_eff*gradT) Q_rev Q_irrev; Q_rev T*dE_rev/dT * i_local; //可逆热项 Q_irrev (eta_act eta_conc)*i_local; //不可逆热重点在k_eff的有效导热系数实测GDL各向异性导热系数设置不当会让温度场出现诡异条纹。有个骚操作把冷却流道参数化扫描配合批处理计算自动寻找最优流速比手动调参效率高十倍。验证环节最刺激的是对比文献[1]的V-I曲线此处应有但没给的对比图。当电流密度突破1.5A/cm²时等温模型预测电压0.25V实际考虑温度骤升后直接掉到0.18V这和J.P.Koza的实测数据误差3%。网格划分有个魔鬼细节催化层至少分3层单元不然活化损失计算会出锯齿状异常。comsol燃料电池模型 等温和不等温的均有 下面是验证图以及参考文献。最后来个压轴技巧在求解器配置里把电荷守恒和动量方程分步求解收敛速度比全耦合快得像开了外挂。碰到发散别慌把电解质相变潜热作为阻尼项临时加入收敛后再逐步移除这招救过我好几次项目进度。参考文献[1] Wang Y, Chen KS. J.Power Sources 2013 (验证曲线原始数据来源)[2] COMSOL官方燃料电池案例库v5.6网格划分基准
