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