关于OpenFOAM可压缩流求解压力方程的问题

[昝磊]

在使用可压缩流求解器时，如果把密度仅仅视为温度的函数，不调用OpenFOAMD的热物理模型。在求解完能量守恒方程后直接用温度更新密度，那么在求解压力方程时是否可以把fvm::ddt(psi,p)改为fvc::ddt(rho)呢？

[dzw05]

理论上可以，这时候rho作为源项出现在方程右手侧。但并没有试过这样处理。

[昝磊]

@dzw05 谢谢！ 请问是改成这样吗

我改成这样之后压力的迭代数上千，发散了。

请问该怎么修改呢？ 谢谢！

[dzw05]

@昝磊 这样改从方程上来说没问题，但你要注意求解之后的relax，以及更新flux的方式。我近期也在改Openfoam的方程，发现这些小地方很容易导致发散。

[dzw05]

@昝磊 其实你这个问题和openfoam中的incompressiblePerfectGas的状态方程类似，理论上不需要改方程，直接用现有的buoyant类的求解器就行

[昝磊]

@dzw05 好的，谢谢你的建议！我再试一试。