[求助]一个非常简单的问题，就是弄不明白

zhoubuaa

我算一个三维的椭圆旋成体的低速无粘绕流问题，理论上壁面上前后两个端点的压力系数应该都为１，因为无粘性损失．（壁面附近的网格较密），但是每次算出来的结果后面的端点压力系数都远小于１，请高手指点一下．多谢．

周华

要看你用的什么格式，控制方程是什么。即使是欧拉方程，也会存在数值粘性，在结果中体现出“粘性”影响。不过，如果前后差的很多，可能就是不是这么简单了。

zhoubuaa

多谢！
我用的是一阶迎风格式．因为是无粘，肯定是欧拉方程了．前后两点的绝对压力相差５０ｐａ左右．另外您说的数值粘性指的是计算误差吗？
不知道周华是不是您的真名字，两个字中间加上一个字，就是我的名字．

周华

是真名，我也是北航毕业的。数值粘性的来源很多，有些格式里面为了格式稳定要增加一些人工粘性项，另外误差有时也可以体现为粘性。虽然用的是无粘的欧拉方程，但是如果壁面边界条件用无滑移条件的话，软件会在壁面与外流之间做一个差值，仔细观察可以看到一个“边界层”。由于这样的粘性，流动就会有机械能损失，在前后驻点就会形成一个压差。

周华

北航4系、5系有很多做CFD的老师，如果离的不远，可以找他们当面请教一下。

小凡

如果我在北京就好了,我这里做CFD的人太少了

zhoubuaa

多谢多谢．

沈沉

最近也发现这个问题，算的达朗贝尔疑题，算不对，也怀疑是人工粘性的问题。
不同格式数值粘性不同，我感觉Quick、二阶、三阶MUSCL粘性明显比较小。

[ 本帖最后由 沈沉 于 2011-10-9 11:26 编辑 ]

onebyte

Quick、二阶等格式本身就是为了减少数值粘性做了改进