关于OpenFOAM可压缩流求解压力方程的问题
-
在使用可压缩流求解器时,如果把密度仅仅视为温度的函数,不调用OpenFOAMD的热物理模型。在求解完能量守恒方程后直接用温度更新密度,那么在求解压力方程时是否可以把fvm::ddt(psi,p)改为fvc::ddt(rho)呢?
-
理论上可以,这时候rho作为源项出现在方程右手侧。但并没有试过这样处理。
-
-
@昝磊 这样改从方程上来说没问题,但你要注意求解之后的relax,以及更新flux的方式。我近期也在改Openfoam的方程,发现这些小地方很容易导致发散。
-
@昝磊 其实你这个问题和openfoam中的incompressiblePerfectGas的状态方程类似,理论上不需要改方程,直接用现有的buoyant类的求解器就行
-
@dzw05 好的,谢谢你的建议!我再试一试。