simpleFoam求解二维翼型Cd不变

[dzw05]

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

[ycx1042341217]

@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)

[ycx1042341217]

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

[李东岳]

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

但阻力系数依旧不正确。

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

[ycx1042341217]

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

[dzw05]

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

[ycx1042341217]

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

[dzw05]

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

[ycx1042341217]

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

[dzw05]

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

[ycx1042341217]

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