CFD中文网

@Andrew_Li 在 LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？ 中说：

dynamic Smagorinsky-Lilly

哪个openfoam版本的？

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

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

第三方写的求解器，如果你们自己还debug不出来毛病，一般只能找作者本人去debug

https://cfd-china.com/topic/677/fluent-空化两相流数值模拟结果绝对压力为负值

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

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

厉害厉害，这个图非常好

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

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

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

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

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

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

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

multicomponentFluid模块

https://cfd-china.com/topic/8018

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

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

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

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

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

const Thermo& thermo是什么变量？

是个Thermo变量。Thermo是一个类型。

也有相关领域的液相绝对负压的文献：CAUPIN, F. & HERBERT, E. 2006 Cavitation in water: a review. C. R. Phys. 7 (9–10), 1000–1017.

我第一次看到负压的讨论，值得深入研究一下

openfoam的adiabaticPerfectFluid模型，是更新密度的，跟你写的不一样。你确认一下，openfoam中的adiabaticPerfectFluid模型在压力为负的时候，密度可以是正的么

你这太小了，连残差的收敛标准也不过1e-6左右，你试试加到1米左右的速度看看

你得速度为什么这么小，e-8量级？