如何解释计算喷嘴缩口处的负压问题？

[mohui]

如图所示 用interfoam算的结果这里出现很大的负压。主要疑惑是如果基于伯努利方程，流通面积减小是会导致压力降低，但是也不至于是负压吧，难道这个是符合实际情况还是求解本身的问题？
请问缩口出现负压正常吗？

[赵一铭]

你是说压力小于0？ 不可压缩模拟经常压力小于0啊。不可压缩模拟里面压力是相对的。你可以设置参考压力为50 或者500看看，应该就不是小于0了。这只是个相对压力。

[mohui]

@赵一铭 参考压力难道不是0吗？另外缩口处负压很大，达到-e7次方。这难道也是正常的吗？

[赵一铭]

那可能有点不正常了，我觉得不正常，可能你设置错了。

参考压力CFD界发过一个文章，你可以找找看。是可以随便设置的

[mohui]

@赵一铭 ! ，这是我见过为数不多一篇论文的截图，就我目前我所知的软件算的结果都会是这样子的结果，但是大部分的人貌似都不会讲这个问题或者自己修改一下，毕竟云图是最不能反应真实情况是什么样子的。

[赵一铭]

那还是相信论文比较好，哈哈。不过你那个速度也很大，600m/s了，确实是不可压缩？

[mohui]

@赵一铭 对啊，入口压力很高，加上缩口，速度确实很高，然后基于伯努利方程确实会得到负压，但是感觉怪怪的。另外拿可压缩求解器比如compressibleinterfoam算的话，其实会有一个更诡异的现象，气体会填充这个负压区域强行使得负压区域大于0.

[bfamq]

以前也遇到这个问题，入口流率一大，缩口的负压就会很高。

[mohui]

@bfamq 这个问题觉得可以这样子分析，首先得的负压计算结果是没问题，但是不代表着真实情况是负压，只是这个求解器本身就是不可压缩求解器，基于不可压缩的伯努利方程，缩口区域会出现负压。如果真要符合真实的情况，应该从两方面入手，一是是首先得用可压缩求解器，而是需要添加空化模型来修正压力，我认为只有这样子才能算出合理的情况。
另外如果用of中的compressibleinterfoam计算，缩口区域会出现空气相，alpha2=1-alpha1，由于缩口有很强的负压趋势，所以导致alpha1小于1，从而导致出现alpha2。单纯的用compressibleinterfoam算的也不对，正如我前面说的需要添加空化模型。
以上是我对这个问题的理解。

[bfamq]

@mohui 感谢你的解释