没想到今天模拟一个冷态流化床都出现这么多问题。需要通入三种气体组成的混合物，这时入口边界的压力BC_P_g不能为零，否则会报错，但是，需要注意的是，收敛效果对这个压力边界非常敏感，我随意给了个30，结果发散，发散的很厉害，速度超限。然后我看了看官方给的一个流化床燃烧的例子，里面设的是3.04d7（30个大气压），设成这个值就没问题了。
所以，结论是：单种气体算流化床时，边界压力可以给0，但是通多种气体就必须严格考量这一值，可以参考官方案例，比如3.04d7。具体原因可能是用到了理想状态方程，PV=nRT，压力太小的时候，算出的解会异常。只是猜测，还得进一步了解里面的计算原理才行。
还有一点，迭代步长最好习惯性的给定最大值DT_MAX，以免计算过程步长越来越大。