在此记录下放氢过程的软件设置思路
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。
标签:P0,COMSOL,求解,初始值,步长,放氢,氢化物,初始,层流 From: https://blog.csdn.net/qq_24800941/article/details/140912964