FLuent与OpenFOAM中关于压力限制的区别

李东岳

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

打倒小怪兽

@李东岳老师，我还有一个疑问，pMin的值可以设置为负值吗？会有文件对pMin的值必须为非负值的限制吗？
这样是否可以避免直接取消pMin带来的风险？

李东岳

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

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

打倒小怪兽

@李东岳好的，谢谢老师指点迷津，我再自己研究一下。
刚刚尝试了一下，相同的算例使用compressibleInterFoam求解器计算，将pMin值设置为负值后很快就会报错负温度。看来不能简单的直接设置为负

李东岳

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

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

李东岳

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

打倒小怪兽

@李东岳好的，李老师，我去看一下。

我最近看了一下foam-extend版本的compressibleInterFoam求解器，没有耦合温度方程。他直接在pEqn.H文件里面给出的压力、密度的闭合方程——ρ=psi*p，如果是弱可压缩性的计算，是不是影响不大？也可以用来计算弱可压两相流动？

李东岳

是的，强弱都可以，你可以试试，那个不会出现负温度。

李东岳

这个帖子的绝对压力，fluent那面确实出现了负值。

我对负的绝对压力，目前理论上理解不上去。不过也没时间详细看了。我先看看能不能把org那面的compressibleInterFoam的能量方程去掉

打倒小怪兽

@李东岳是的，李老师。
我在使用Fluent计算相同设置的算例时，层流模型、没开能量方程、液体属性设置为可压缩流体。也会出现绝对负压（背压101325Pa，Solution中的限制压力1Pa），流场中会出现几倍于大气压的负值。

如果针对这个问题，在OpenFOAM中去掉能量方程的计算应该会稳定很多。不过，我能力有限，对org版本thermo那块复杂的，我看起来还很吃力，李老师辛苦研究一下！

foam-extend的compressibleInterFoam本身就没耦合能量方程，但是他的液体状态方程也没有Tait方程的形式，这两天我把这个方程整进去了，计算tut算例没啥问题，但是算我那个算例就发散，可能稳定性还是有点问题，还需要再调试一下。

李东岳

我用rhoPimpleFoam、rhoSimpleFoam测试了一下，直接把EEqn.H注释掉就好了。你试一下。在新版OpenFOAM中，原本的fluid模块，替换成isothermalFluid，等同于上面的效果。不过我测试NACA算例，注释掉能量方程之后收敛性比较差。流场结果出来倒是没有问题。

你试一下直接把compressibleInterFoam的TEqn.H注释掉

打倒小怪兽

@李东岳老师我尝试了一下将温度注释掉应该就可以了，不过也如您所说，注释掉能量方程后收敛性不是很好，可能要提高很多精度来计算，我再调试一下。
感谢李老师热心指导

CFD中文网