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]