想在超音速流中计算自己定义的广义力，如何获得表面压强

ranjhzhl

最近做一个数值试验，想计算跨超音速流情况下的广义气动力，问题就出在物面压强的调用上，在fleunt6.3....26的版本中，针对压力远场边界条件是可以用F_P(f,t)来获得face center的压强的，但是在物面上却不能用该函数来获取物面上face center的压强，从而导致广义气动力不能计算。

解算器必须选基于密度的解算器。

有人建议用紧邻单元c0的压强值来代表，可是我以现翼型完全对称的机翼在0度攻角下，计算出来的广义力不为零。理论上说，对称翼型0度攻角M=0.9时升力系数是零，上下两表面压强完全相同，方向相反，积分后，升力系数所以是零，广义力也应该是零，可我现在算出来的结果不为零，原因未知。

当前所用的广义力计算公式：压强*面积矢量*振型(即面中心坐标变形前后得到的矢量差)，然后整个物面积分，主要代码有：F_AREA(A,f,t);  c_p(c0,t)*(A[0]*x[0]+A[1]*x[1]+A[2]*x[3]),请高手指条出路，如何解决当前的问题。谢谢

[ 本帖最后由 ranjhzhl 于 2009-6-20 13:30 编辑 ]

gbd20022002

F_P(f,t)可以获得壁面压强