本文来源于我本科毕设,在利用COMSOL进行仿真模拟吸放氢的过程中,也遇到了很多困难,同时也了解到有很多同行,正在开展这方面的研究,所以把利用COMSOL进行仿真模拟储氢合金吸氢过程的建模方法整理分享了出来,希望有助于后续学者的储氢罐结构设计和换热设计。
一、模型与假设
上图给出了本研究中仿真所采用的储氢罐结构示意图。储氢罐为圆柱形,底部直径为d,罐体高度为L。储氢罐中填充着金属氢化物(MH),其孔隙率为ε。氢气从顶部恒压流入储氢罐。储氢罐底部和外表面均与恒定温度为Tf 的换热流体接触并通过换热流体进行对流换热,有效换热系数为h。
在本文的模拟过程中,采用以下假设:
(1) 氢气为理想气体;
(2) 粉末床合金颗粒均为球体且分布均匀;
(3) 罐体内部局部热平衡,并忽略辐射散热;
(4) 金属氢化物的渗透率在吸氢过程中保持不变。
二、建模过程
1.开始界面
选择模型向导,二维。
添加物理场接口:多孔介质传热、达西定律、域常微分和微分代数方程模块。
注意修改域常微分方程模块的单位。
选择瞬态。
2.参数与变量
在参数界面导入LANI5-data文件,获取基本物性参数。
选择几何1,添加矩形。
输入宽度d,高度H,位置基准x为d/2,点击构建选定对象。
形成联合体,点击构建选定对象。
右键定义,插入变量
导入变量文件-LANI5-variable。
3.多孔介质传热
选择多孔介质传热,将参考温度改为参数表中的T_ref。
多孔介质传热-多孔介质-流体:修改导热系数、密度、恒压热容为用户定义,输入对应参数的名称。
多孔介质传热-多孔介质-多孔基体:修改孔隙率、干本体导热系数、干本体密度和干本体恒压热容为用户定义,输入对应参数的名称。
多孔介质传热-初始值:温度修改为起始温度T_ref
右键多孔介质传热,添加热源和热通量。
多孔介质传热-热源:域选择-选择右侧图像1,修改下方热源项Q0,设定为S。
多孔介质传热-热通量:修改热通量类型为对流热通量,传热系数设置为300,外部温度设置为T_ref1,选择下方边界和左右边界。
4.达西定律
达西定律:因变量改为P。
达西定律-多孔介质-流体:修改密度和动力黏度为用户定义,输入对应参数的名称。
达西定律-多孔介质-多孔基体:修改孔隙率为用户定义,输入对应参数的名称,设置渗透率为1e-8。
达西定律-初始值:修改压力为P_ref。
右键达西定律,添加质量源和压力。
达西定律-质量源:域选择,选择右侧整体,下方质量源设置为-m。
达西定律-压力1:选择上边界,压力设置为P_ref。
5.域常微分和微分代数方程
域常微分和微分代数方程-分布式常微分方程:修改源项为m,阻尼或质量系数为1-e_p,质量系数为0。
域常微分和微分代数方程-初始值:
修改rhos为8410,修改初始时间导数为:Ca*exp(-Ea/(R_const*T))*log(P_ref/200000[Pa])*(8521-8400)[kg/m^3]。
6.探针、网格与计算
右键变量,添加探针-域探针
表达式,设置为反应分数f。
网格-全部构建
研究-瞬态:修改时步为0,5,3000,计算。
储氢合金吸氢仿真结果如上图。
三、不同合金与变工况
改变合金种类,则可在参数与变量中,修改合金的热物性参数,平衡压方程与域常微分的初始值等参数。