CFD中文网

李

我在无痛苦ns方程笔记里面更新了密度、温度的更新流程，应该有点关系

xpqiu

@ChangranLv
你这个曲线有一个明显的规律：网格越多，中心速度就越小。很大一个可能性是你的计算根本没有收敛。用simpleFoam 求解稳态问题，网格越多收敛会越慢，需要越多迭代步。

李

是。对于理想气体，负的温度、压力、密度都不合理。不过你那个绝热理想液体，在p为负的时候，密度确实为正。在结合那个review说的绝对负压的存在，应该需要详细研究下。

不过如果没看到那个review，过去10多年我一直以为绝对压力一直为正。需要好好看看这个。

李

在你提到那个绝对负压之前，我从来没见过设置为负值的压力。从状态方程来看，基于可压缩性的算法psiThermo，p是负值是完全不可能的。rhoThermo中p和rho的关系可能更复杂，比如绝热理想液体，可能一个负的压力可以预测一个正的速度。但是听起来基本挑战之前的CFD数值常识。

你要想尝试绝热理想液体能否允许负压的存在，你需要做个算例做个测试。比如rhoPimpleFoam算一个特别简单的算例。

李

那个文章里面说确实可以存在负压。但是不表明OpenFOAM允许负压的存在。你把pMin取消掉，存在很大的风险。

李

multicomponentFluid模块

李

如果能满足这个要求，可以发上来我看看

李

那个密度公式，不会出现负的密度。
不过p倒是允许出现一定的负值，看起来没做限制，但要保证密度不是负值。

你得问题是什么？OpenFOAM里面的pMin的作用么？

李

HZDR那个确实感觉嘎嘎一嘎子，我还没细看算法。fluent那面我更不了解了

李

拉格朗日颗粒，还是连续相颗粒？