关键词:
数值模拟
机器学习
大涡模拟
高精度数值模拟
摘要:
翼型绕流是流体力学研究重点之一,流动相对简单但广泛应用于各种流体机械及飞行器设计中。且翼型绕流可以通过圆柱绕流的势流变换得到,因此本文以圆柱及NACA0018翼型作为研究对象,采用不同湍流模型计算了研究对象在定常与非定常状态下的流场。通过对比分析研究了数值模拟结果,选用二维模型下SST k-ω湍流模型对研究对象的计算结果作为机器学习样本,对机器学习方法的流场复现方法进行了研究,选用多层前馈神经网络模型实现机器学习方法,并根据复现出的流场计算翼型的气动特性,主要方法及结论如下:1、采用S-A、SST k-ω与LES湍流模型,对NACA0018翼型在不同攻角下的流场进行数值模拟,研究表明,S-A模型在0-4°攻角下计算误差可以小于5%,当攻角超过8°时,S-A、SST k-ω方法开始出现较大误差。使用非定常方法进行计算时,SST k-ω与LES湍流模型在攻角α<8°时均能准确预测升力系数,但当α>8°时由于发生了流动分离,SST k-ω湍流模型难以准确预测,另外LES湍流模型可以较为准确的计算出翼型的阻力系数。最后对比了圆柱绕流LES方法与实验的对比,通过对比实验流场与数值计算流场发现SST k-ω方法无法得到圆柱绕流背压面旋涡被主流拖拽破裂的过程。综合对比三种数值模拟在圆柱绕流与翼型绕流中的计算结果可以发现LES在宏观外特性与流场特性均具有较高精度。2、本文采用数值模拟方法获取数据集,数据集的好坏直接影响预测准确性。圆柱绕流因计算简单,且可以通过势流变换得到翼型绕流,因此本文首先采用数值模拟方法对雷诺数为Re=3900的圆柱绕流进行计算,对比圆柱绕流试验压力系数与数值模拟计算压力系数可以发现LES湍流模型与SST k-ω湍流模型计算结果相近,但计算时长相差近35倍。综合考虑计算时间和准确性,选用二维模型下SST k-ω湍流模型建立直径为0.05m的圆柱在Re=3500-5500下压力场速度场数值模拟结果作为学习模型的数据库,研究不同深度、训练步数、学习率和不同网络结构的预测准确性。结果表明,保证训练步数在500-750步之间,学习率为0.1保持不变时预测结果很差,在减小学习率到0.01时可以获得较好的结果,但学习率继续减小会导致神经网络难以收敛,需要增大训练步数方能收敛。另外神经网络在层数增多时会有较好的预测结果,但当前馈网络超过6层时,预测准确性并没有很大的提高,认为使用6层前馈网络即可达到圆柱绕流预测的目的。最后加入0.005的正则化与1%权重衰减可以有效防止过拟合,提高预测准确性。3、将圆柱绕流优秀的预测模型应用于翼型绕流预测。首先利用数值模拟方法获得NACA0018翼型在攻角α=2°、4°、6°、8°,雷诺数Re=0.1-1.6×10范围内共计140组计算结果作为机器学习的样本数据库,建立基于前馈神经网络的机器预测模型,以绕流雷诺数与攻角作为自变量输入神经网络,通过神经网络的计算预测翼型附近压力场速度场。结果表明,采用八层全连接神经网络即可获得准确的计算结果,并可以通过预测得到的表面压力系数反算升阻力系数,发现机器学习方法能获得很好的准确性。对比Xfoil、数值模拟与神经网络计算得到的表面压力系数,发现SST k-ω湍流模型未算出翼型分离点,导致神经网络无法计算出分离点压力变动,表明机器学习方法准确度取决于数据库样本。4、通过计算矢量场速度之间的相关系数可以发现X、Y方向的速度矢量大小相关系数为0.3,相关性较弱,因此采用相互独立的神经网络预测速度矢量场。相比压力场、速度场的预测,因缺少数据间的约束,需要将神经网络训练步长增加至600步以上。并考虑到流场中驻点及部分区域Y方向速度接近于0,将计算结果中绝对值小于10的数据归零处理,最终可以将矢量场预测的平均误差降低至7%以下。