buoyantfoam计算发散求助

[ShaneHEEE]

有米有试试p_rgh除了出口fixedvalue，其他都给zeroGradient

[Miraitowa]

@ShaneHEEE 试过的，一样也是发散的，就很怪

[Afr1yne]

@Miraitowa 我看你说的是能算几步才发散是吧，可以把图放上来看一下。我之前用buoyantFoam发散是因为尺度问题，你把整个计算域缩小一下试试能不能跑

[Miraitowa]

@Afr1yne 您好，这个是计算域的图：

大小是19.5m9.6m9m的，下面是计算两个时间步后的结果：

我仔细看过之后发现速度是从出口开始发散，明显大于理论值，这个是残差图：

[Afr1yne]

@Miraitowa 把出口也设成fixedFluxPressure; 然后设置个pref试下

[Miraitowa]

@Afr1yne 好的谢谢，我试一下

[Miraitowa]

@Afr1yne 您好，我试过了，还是发散的，但是变成入口处速度明显大了

[李东岳]

你用单相流求解器不可以么？为什么要用buoyantFoam

[Miraitowa]

@李东岳 因为想要考虑重力分层，并且将温度和流场双向耦合起来

[李东岳]

你试试老版本的bouyantBounessneqPimpleFoam，或者自己写一个Bounessneq近似的求解器，直接求解T方程的，不需要求解能量方程，并且是不可压，不存在负温度的问题

[Miraitowa]

@李东岳 好的谢谢老师，我试一下

[李东岳]

buoyantBoussinesqSimpleFoam.tar.xz

给你做了一个openfoam10下的稳态求解器。你试试