COMSOL金属氢化物-放氢过程

在此记录下放氢过程的软件设置思路

1、采用的是"达西定律"+"层流" 物理场，其中"层流"物理场选择了”弱可压缩流动“，这里主要是选择”可压缩流动“的话，算出来的瞬时流量值跟实测差距太大了。

2、设置"达西定律"的压力初始值为P0，"层流"的压力初始值也为P0，"层流"的出口设置压降函数，解决“找不到一致的初始值”问题，主要是一开始设置出口为1atm的话，压力场过渡太剧烈，不好找到初始解，要么就先算个稳态流场，这种也能求解。注：P0是根据实测PCT的插值曲线(函数名称为：p_iso)得到的，根据初始质量储氢量wtw0，和初始温度条件T0，然后根据函数p_iso(wtw0)*exp(-deltHd/R_const*(1/T0-1/303.15[K]))求解获得。

3、 瞬态求解器注意"初始步长"选择框勾上，设置初始步长为0.001s，最大步长可以设大一点，后续求解的收敛性仍然较好。若"初始步长"选择框不勾上，前面0~5s内的求解容易发散，然后求解失败。

4、最终CFD算出的瞬时流量和累计流量，跟实测的对比几乎吻合，拟合效果OK。