buoyantPimpleFoam计算问题

[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看看

[Afr1yne]

@李东岳 东岳老师我试了下fixedMean 值给的101325跟其他边界和内部场一样 还是不行

[Afr1yne]

@李东岳 东岳老师 我想构建个压力梯度初始场，目前用codestream还不太行 您能看出来这篇论文是通过什么实现的吗

[Afr1yne]

@李东岳 东岳老师 这个case说的修改速度压力耦合策略的参数 是哪一个参数？
最近发现了大尺度大气流动计算不收敛的原因可能是压力不平衡 我也尝试了不加温度剖面 计算确实会收敛

[李东岳]

压力速度耦合那个我当时咋处理的有点忘了 ，可能是松弛因子，要不就是换成SIMPLE或者SIMPLEC

不过fixedMean看起来更加符合物理，我在讲课的时候都让他们改fixedMean