simpleFoam求解二维翼型Cd不变


  • OpenFOAM讲师

    @ycx1042341217 可以看一下你的数据吗?
    Cd系数基本不变?你是说绝对值吗,看看相对值有多大变化。
    至于阻力算不准的问题,多数是由于粘性力没算准。原因比较多,(1)什么湍流模型?(2)第一层网格法向是否足够密(y+)?(3)机翼前后缘的网格是否足够密,是否正确刻画了机翼的弯曲?



  • @dzw05 老师,我采用的是NACA0018翼型,弦长1米,速度1.463m/s,雷诺数100000,采用的是SA模型。不同攻角下升力系数较为吻合。阻力系数基本维持在14.多。之前划分网格的时候是跟着网上的一个教程划分的,没有考虑到y+。刚才我重新划分了网格(根据我的数据,按照y+等于1计算,我的壁面第一层网格高度大约为0.08mm。),但是在OpenFOAM里面checkMesh后,提示我High aspect ratio cells found, Max aspect ratio: 10536.4, number of cells 1020。老师我是不是需要重新划分网格?那为了避免这种情况,网格的纵横比最好不要超过多少呢?老师您看看还需要我的哪些数据吗?谢谢老师!


  • OpenFOAM讲师

    @ycx1042341217 你这个aspect ratio太大了,如果法向网格0.08mm,那么流向网格有800多mm?OpenFOAM中建议的最大aspect ratio为1000,需要重新划分一下。



  • @dzw05 老师您好,我的第一层网格高度大约为0.18mm,昨天跟您说错了。刚才重新划分网格后,checkMesh显示正常,但是进行计算后升力系数和阻力系数都与实验数据有偏差。老师我把我的网格文件,controlDict文件和初始条件传上来,您能帮我看看问题出在哪里吗?谢谢老师!4_1523608055301_U.pdf 3_1523608055301_p.pdf 2_1523608055301_nuTilda.pdf 1_1523608055301_nut.pdf 0_1523608055300_controlDict.pdf [网格](链接: https://pan.baidu.com/s/129VWe7pswJfTMWmG4oHxrg 密码: x7g4)



  • @dzw05 老师,我这两天重新划分网格的时候发现之前计算过程中郎库数非常大,现在我根据郎库数为1重新设定了时间步长,但阻力系数依旧不正确。我在gambit网格质量检查的时候,远场质量为0,翼型附近质量为0.1-0.2。不知道这样的网格质量可以吗,或者还有什么其他的可能问题吗(比如说我的nut,nutlida值会不会设定有问题),因为我现在修改了很多次网格,但是阻力系数总是维持在10以上。谢谢老师!


  • 管理员

    0_1523783724695_捕获.PNG

    建议你在这个方向上增加网格分辨率。

    但阻力系数依旧不正确。

    流场正确,但是计算不正确的话,不是特别好调节。网格上好好研究研究,看看能不能提高结果。另一方面是看看你用的湍流模型理论上是否能预测你的流场。



  • @东岳 老师您说的分辨率指的是gambit中划分网格时这个方向的successive ratio吗?


  • OpenFOAM讲师

    @ycx1042341217 边界条件可能有问题。
    (1)对于速度,你现在进口和出口都是freestream,而在openfoam中freestream其实是继承自inletOutlet条件,对于进口采用freestreamvalue,对于出口采用零梯度条件.
    (2)对于压力,你采用的是freestreampressure,其继承自零梯度条件。
    所以你现在并没有显式给定边界压力,这样的边界组合是不对的。你可以尝试在出口处将压力条件设为fixedValue,速度设为zeroGradient。
    我感觉openfoam在外流计算,或者说在空气动力学计算方面,做的不是特别好,不知道@东岳 老师怎么认为。



  • @dzw05 老师,我修改了边界条件,对网格也进行了一些修改,但是重新计算后阻力系数还是基本没变。会不会是采用的SA模型本身对于这种情况模拟就不是很好?


  • OpenFOAM讲师

    @ycx1042341217 从理论上来说,SA算这种外流机翼问题应该非常好,因为Spalart本人就是波音的,SA模式针对飞机的计算做了非常多的优化。但是,理论上的事谁能说得准呢,OpenFOAM中的SA我还没有用过,不知道是哪个版本,你可以试试其他的湍流模型。



  • @dzw05 老师我排查出来是我的湍流模型文件里面的参数设置出了问题,所以导致阻力系数一直存在问题,现在已经解决了。十分感谢老师这几天的帮助!


  • OpenFOAM讲师

    @ycx1042341217 能方便说一下哪个参数设置不对吗?让大家都知道怎么避免:)



  • @dzw05 湍流模型中有一个运动黏度nu,正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值,由于之后的计算我都用的是这个错误的湍流模型文件,导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。


登录后回复
 

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