buoyantPimpleFoam计算问题

[Afr1yne]

请教下各位老师
我在计算二维山地风场 同时考虑温度 目前buoyantPimpleFoam计算一直发散

我把重力g设置成0的话 计算也会发散 但是速度场看起来会比较正常 是压力的问题吗？

@李东岳 手动艾特下大佬

[李东岳]

你这个尺度比较大
换成of8自带的bouyantBoussinesqPimpleFoam你试下，应该没问题

[Afr1yne]

@李东岳 李老师 我用bouyantBoussinesqPimpleFoam算了一下 看起来更合理 是不是就说明跟是否可压有关？

[李东岳]

这个求解器有bug 目前我的方法是把bouyantBoussinesqPimpleFoam移植到openfoam10了 没有时间研究他们的bug也..

[yjc95]

@Afr1yne 您好，请问您发散的结果用的是of8算的吗？问题解决了吗？

[Afr1yne]

@yjc95 我用的v2112，目前发现是跟密度有关，密度设为常数计算结果看起来比较合理，但是具体原因还不清楚

[yjc95]

@Afr1yne 谢谢

[runtowhere]

出口压力边界可以换成zeroGradient，或者fixedFluxPressure或者waveTransmissive 试试；
如果有ph_rgh的话，把top设成fixedValue；其他设成fixedFluxPressure。
如果top边界U是slip之类的可能更好收敛，PIOV会更复杂点

[Afr1yne]

@runtowhere 老师您好 我现在p_rgh边界设置的都是fixedFluxPressure这样可以吗 之前是出口fixedValue其他的全是fixedFluxPressure

[runtowhere]

@Afr1yne 都可以试试啊

[yjc95]

@Afr1yne 大佬，您是用of哪个版本下的bouyantBoussinesqPimpleFoam呀？我用of8提示没有这个求解器

[Afr1yne]

@yjc95 你看看buoyantFoam或者buoyantPimpleFoam 现在是把boussinesq作为状态方程放在physicalProperties里面了

[yjc95]

@Afr1yne 这俩我都试过了，of8下没有buoyantFoam，buoyantPimpleFoam在改变膨胀系数beta值的时候会出现非物理现象（https://www.cfd-china.com/topic/7611/buoyantpimplefoam相关问题/63
这个帖子的最后几个回复有图），东子老师说of8的buoyantPimpleFoam有bug，建议我试试buoyantBoussinesqPimpleFoam，我看您在of8下可以用buoyantBoussinesqPimpleFoam，但我试了之后发现用不了。

[Afr1yne]

@yjc95 我用的是v2112

[yjc95]

@Afr1yne 东子老师的虚拟机里没有v2112呀，是自己装的嘛

[Afr1yne]

@李东岳 东岳老师 你说buoyantPimpleFoam有bug 是新版本的有问题 还是所以版本都有问题？
我发现速度场垂直方向有问题 考虑是浮力项或者是p_rgh = p - rho*gh这两方面出现了错误
经过我测试
1 重力设为（0 0 0）
2 状态方程用rhoConst
3 尺度缩小
这三种情况计算看起来都比较合理

[李东岳]

是的，以上3个都没问题。

是buoyantPimpleFoam+bousinessq近似的问题。两者结合的问题。我还没详细debug咋回事。只是发现了这个问题。

[Afr1yne]

@李东岳 我用buoyantPimpleFoam + perfectGas 不好使
我想算大尺度同时考虑温度湿度，东岳老师还有求解器推荐吗？

[李东岳]

https://www.jianguoyun.com/p/DdwnJOEQwr7rChjP-98FIAA

这是我之前移到of9的，你看看能用这个就用这个

[李东岳]

另外你看看你得算例跟我这个是不是类似的 http://dyfluid.com/openfoam.html#id4

如果是的话把你的压力出口改成fixedMean看看