考虑浮力的计算中边界条件设置

[王慧博]

我想做一个管道流动，管道中间放置一个颗粒，温度高于来流。在考虑浮力的求解器，比如buoyantPimpleFoam和buoyantBoussinesqPimpleFoam中，边界条件设置时遇到了问题。我参考OpenFOAM自带的算例，发现里面与浮力相关的算例，速度边界为(0,0,0)，压力边界为fixedFluxPressure。但是我在不考虑浮力的算例中设置，入口边界条件为：速度固定值，压力零梯度；出口边界条件为：速度零梯度，压力固定值。但是把这些边界条件代入到考虑浮力的条件中，会出现问题。
所以想求助一下老师们，在考虑浮力时，入口和出口的速度、压力应该采用什么样的边界条件比较合适呢？

[李东岳]

http://www.dyfluid.com/theory.pdf 看下144页

[王慧博]

@李东岳 感谢李老师，我看了您的笔记，上面写到了压力边界出口为prghPressure条件。那么速度出口条件是还是用零梯度就可以么

[李东岳]

@王慧博 是的

[王慧博]

@李东岳 东岳老师您好，我后面设置了一个算例，Re=100，边界条件如下所示： ，如果把重力加速度设置为0，计算正常，但是加入重力后，计算到后面就会发散，从云图上的感受是速度堵在出口没办法出去，下面是发散前我截出来的某个时刻的云图：

能否请教一下老师，帮我分析一下问题可能出在哪呢？