buoyantfoam计算发散求助

[Miraitowa]

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

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

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

[Afr1yne]

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

[Miraitowa]

@Afr1yne 好的谢谢，我试一下

[Miraitowa]

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

[李东岳]

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

[Miraitowa]

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

[李东岳]

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

[Miraitowa]

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

[李东岳]

buoyantBoussinesqSimpleFoam.tar.xz

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

[Miraitowa]

@李东岳 老师您好，我用这个求解器的时候会提示

FOAM FATAL ERROR: 
    request for thermophysicalTransport thermophysicalTransport from objectRegistry region0 failed
    available objects of type thermophysicalTransport are
0()

但是已经反复确认过在constant文件夹里设置了thermophysicalProperties，请问是什么原因呢

[李东岳]

问题在alphat，你去研究一下我那个算例，然后再用于你的算例

[李东岳]

buoyantBoussinesqSimpleFoam10-12.tar.xz

这个压缩包里面有两个笨笨的求解器，分别适用于

• OpenFOAM-10
• OpenFOAM-12 （模块化）

算法主要基于老版本的buoyantBoussinesqSimpleFoam。这个求解器不需要求解能量方程，同时将动量方程处理为不可压缩。求解器稳定的多，不会出现负温度。但是在OpenFOAM-7之后被删掉了。

压缩包里面同时含有1个测试算例，分别适用于两个版本的求解器。对模块化求解器感兴趣的可以研究研究。我在植入过程中，模块化求解器的植入相对麻烦，主要是成员函数比较分散。OpenFOAM-10作为一个顶层的求解器在写的时候一气呵成。二者体验还是挺大的。

如果想在OpenFOAM-12里面植入求解器，可以对比下我写的求解器在10下以及12下的区别。