potentialFoam计算的驻点压力系数为啥不是1？

[lwjetmann]

算垂直平板绕流势流，入流速度为1、压力接近0（出口压力为0）。前驻点压力应该为1吧？咋是1.427？

手机拍的图，搞不会paraview显示点数据时咋截图 →_→

[lwjetmann]

添加一个自带potentialFoam的cylinder算例结果：
压力分布

速度分布

p+0.5*U^2在前后驻点及上顶点分别为0.39，0.37，1.28，不是常数

[李东岳]

你用-writep参数了么

[lwjetmann]

@李东岳 感谢老师回复，用了的。

runApplication $(getApplication) -withFunctionObjects -writePhi -writep
runApplication postProcess -func grad"(Phi)"

不writep压力都是0

[李东岳]

我这主要搞方程离散求解的，流体力学基础比较差，还不知道驻点压力是啥

驻点压力就是$p+0.5|\bfU|^2$么

[lwjetmann]

@李东岳 是我想得不对吗？这难道不就是伯努利方程 ，势流情况下总压守恒？入口速度U=1，压力接近于0（因为给的是出口压力）；那么在流场中速度为零的点，压力系数为1。
网上没搜到类似的帖子，顶多说翼型计算时驻点压力系数不为1，那个有粘性的影响好理解。potentialFoam就是算势流的吧，怎么想也不觉得会出现这个问题。难道potentialFoam算的压力不是真实压力？

[李东岳]

难道potentialFoam算的压力不是真实压力？

在老版本的p是速度势的概念。但是新版已经改过来了。

[lwjetmann]

@李东岳 我用的v2006，Phi是速度势，用-grad(Phi)计算的速度和输出的U有些差异，尤其是在远场边界。