openfoam SA模型计算空翼出现Cl和Cd算不准的问题?



  • 用SA湍流模型计算的空翼,算是很常规的问题了,但是结果总是不对,请大家指教一下。

    NACA0015,攻角5度,速度17m/s,不可压缩流动,simpleFoam求解器。

    远场边界条件

    网格改了很多次,y+有低于1的,30以上的

    但是最后我计算的Cl=0.57(实验值0.5438),Cd=0.4(实验值0.0141)

    目前我认为的是调试湍流边界条件,nut和nuTilda,但是收敛很差,如果nut或者nuTilda给的值太低就会发散,震荡等等。

    我调试的依据是

    3_1482504095163_WechatIMG14.jpeg 2_1482504095163_WechatIMG13.jpeg 1_1482504095163_WechatIMG12.jpeg 0_1482504095162_WechatIMG11.jpeg

    如果取比较低的湍流度,无限远 nut约等于nuTilda

    我用nut=0.01~0.2的值都是不行的,要么震荡要么发散



    1. 对流项div(rhoPhi, U)div(phi, nuTilda)使用高阶格式;
    2. 网格,做一下网格尺寸无关性实验。


  • 你好,我换了很多套网格,结果都出奇的一致,跟实验对不上,至于scheme也换成了其他格式试过了,也不行



  • @yhdthu 您好,请问这个问题您解决了么,还有nut和nuTilda是怎么计算的呢



  • @yhdthu

    你无量纲化对不?可以找别的软件算算试试,会不会是SA算这种就不容易对。

    速度没有意义,得看Reynolds数。

    Cl和Cd有个参考长度,你取对了没?

    simpleFoam传说中因为Rhie-Chow插值问题导致收敛结果不是和松弛因子有关系么?你把relaxation关了试试?


  • 版主

    @yhdthuopenfoam SA模型计算空翼出现Cl和Cd算不准的问题? 中说:

    y+有低于1的,30以上的

    这样不行,机翼表面 y+至少要小于5。
    不知道你网格质量如何,震荡多半是网格质量不好引起的。还有,nut 的边界条件你用的什么?



  • 版主,请问您y+小于1的例子算好了么,我现在要计算一个算例 Re=6000000, 很多套网格,y+分布0.03到3,但是现在完全算不出来。用同样的边界条件和y+试了圆柱绕流,没有问题。但是机翼总是跑不了,不知道为什么。我是用point wise画的网格,在point wise里检测网格正交没问题,但是在openFoam里面检测总有50几个面非正交。不知道为什么? @yhdthu @xpqiu @lv1995


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待