用mfix模拟流化床时边界压力需要注意的问题

没想到今天模拟一个冷态流化床都出现这么多问题。需要通入三种气体组成的混合物，这时入口边界的压力BC_P_g不能为零，否则会报错，但是，需要注意的是，收敛效果对这个压力边界非常敏感，我随意给了个30，结果发散，发散的很厉害，速度超限。然后我看了看官方给的一个流化床燃烧的例子，里面设的是3.04d7（30个大气压），设成这个值就没问题了。

所以，结论是：单种气体算流化床时，边界压力可以给0，但是通多种气体就必须严格考量这一值，可以参考官方案例，比如3.04d7。具体原因可能是用到了理想状态方程，PV=nRT，压力太小的时候，算出的解会异常。只是猜测，还得进一步了解里面的计算原理才行。