numeca的翼型绕流计算问题

qazxsw

[这个贴子最后由qazxsw在 2004/02/01 10:34am 第 1 次编辑]

我对NACA0012做了绕流的欧拉计算
马赫数0.7，攻角3度
结果的-Cp与文献差别很大
而且升力系数阻力系数也差的特别多
不知怎么回事？

numeca

1. 对于FINE GUI 中输出的CP值并非为你所期望得到的值。该值可以在CFVIEW中进行计算。你可说出在CFVIEW中的参考压力、速度以及密度的取法
2. 对于EULER计算，阻力系数如何得出？
3. 必须注意，对于EULER计算，则需要用EULER网格而不建议用粘性网格，否则数值误差很大。
如可能，请将相关文件传给我们，可以一起研究一下。
support@numeca-beijing.com

qazxsw

谢谢，详细问题见email

qazxsw

对NS解
网格的影响太大了
我不断加密，细化，
越细阻力越小，升力变化幅度小一点
但升阻比变化太大了从30多到70多

numeca

通过对该算例的分析，答复如下：
1．对于该计算，首先存在的问题是采用粘性网格来进行无粘计算，这势必造成很大的数值粘性，并造成一定误差；
2．对于该算例的EULER计算，存在的最主要问题是边界条件设定不对。在计算中，来流静压为65286Pa,来流温度为182.05K，对应的来流音速应当为270m/s。在计算中边界条件中所设定的来流合速度为240.4m/s，这意味着来流马赫数不再是0.7，而是0.88左右（通过CFView也可以看出），这是导致结果与文献匹配不上的根本原因。
3．解决方法：
  a.调整网格，对于EULER计算，采用无粘网格；
  b.调整边界条件。附件Euler_NACA0012_To_WXF.rar压缩文件（上载至ftp）中包括了新网格以及计算结果，其中包括来流马赫数0.7(攻角3度)以及来流马赫数0.8(攻角1.25度)的计算结果。附图1给出了来流马赫数0.7时的Cp值分布，供参考。计算时取来流条件Ps = 101325,Ts=293K，音速343m/s。计算结束后，在CFVIEW中取出远场平均静压、密度以及平均速度，作为参考值进行无量刚化。此时升力系数CL = 0.491, CD = 0.01413，与文献中基本相同（文献中为计算值而非试验值，另请注意在文献中给出的阻力系数值为1.31´10-3是错误的，应当为1.31×10-2）。在马赫数为0.8（攻角1.25度时），计算得到的升力系数为CL=0.282，CD = 0.01489，与文献中给出的值亦基本吻合。

numeca

网格对于阻力的影响是不用质疑的，因为无论是激波阻力还是摩擦阻力，都和网格有很大的关系。
另外，在FINE GUI中，通过设定IDCDP为1,0,0所输出的阻力并非翼型的实际阻力而紧紧是合力在X轴上的投影。阻力应当是合力在来流方向上的投影。此问题已经通过详细的图形解释email给您了。

qazxsw

谢谢numeca！
请再讲讲无粘网格应该具有什么特征