储氢合金/金属氢化物LaNi5吸氢过程仿真模拟（COMSOL）

        本文来源于我本科毕设，在利用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，计算。

储氢合金吸氢仿真结果如上图。

三、不同合金与变工况

        改变合金种类，则可在参数与变量中，修改合金的热物性参数，平衡压方程与域常微分的初始值等参数。