《RAE2822翼型流场的fluent计算.docx》由会员分享，可在线阅读，更多相关《RAE2822翼型流场的fluent计算.docx（8页珍藏版）》请在三一办公上搜索。
1、word格式-可编辑-感谢下载支持RAE2822翼型跨声速绕流的CFD计算摘 要：使用ANSYS对RAE2822翼型进行网格划分，之后导入fluent中进行 计算。通过对不同的模型和不同边界层网格的计算，采用控制变量的方法分组 比较分析，并将计算结果中的压力系数与试验数据以及组内数据进行对比分 析，以验证FLUENT计算结果的准确性。关键词：RAE2822，控制变量，ANSYS，FLUENT引言：本文研究了速度场来流条件为Ma=0.729, a =2.31的情况下各种状态下的计算结 果。计算状态分别为无粘流动（欧拉方程、无附面层网格）；至少3种不同湍流模型计 算粘性绕流（同一带附面层网格，y+
2、=30）。采用S-A湍流模型，建立4种不同y+的网 格计算（y+1, y+=10, y+30, y+=50）。对y+=50的网格，额外采用流场求解网格自适 应功能（基于压力梯度）进行计算。以及在例8 = 0.75, Re=5X106，迎角a=-1, -0.5, 0, 0.5, 1的条件下对DLR-F6翼身组合体的绕流进行了数值求解，并将所得结果与实验 结果进行对比分析，对FLUENT软件计算三位复杂外形绕流的准确性进行验证。通过对不同状态下的计算结果分析飞行器的气动特性，并且将数值计算的结果与 相应的试验数据进行比较，从而对数值计算结果进行验证。一、RAE2822翼型1.1 RAE2822翼型
3、二维模型图1 SAE2822机翼翼型1.2计算初始参数流体介质：理想空气来流条件马赫数Ma=0.729攻角 a =2.31来流的压力和温度：P=101325Pa, T=300K声速 v= (KRT) 0.5=347.19m/s来流速度=Ma*v=253.10 m/s二、RAE2822翼型网格的划分和参数的设置2.1 RAE2882翼型网格的划分本文由机翼的特征长度L=1m、马赫数为Ma=0.729和雷诺数Re =6.5X106根据计 算可得到不同附面层的第一层厚度表1所示：根据表1的数据利用ICEM软件划分不同第一层厚度的边界层网格和无边界层厚 度的网格，以y+=30为例，网格图形如图2所示：
4、y+0.80.0000033100.0000041300.0001259500.0002099图2 y七30的网格表1不同附面层网格第一层厚度h WHHHnMmninilllii I II I III mill2.2参数的设置边界的温度T=300K；对于攻角的处理是通过在FLUENT计算时设置的来流速度的方向来体现，如图3所示:Turbulence图3边界条件的设置三、计算结果和数据处理分析3.1无粘模型和湍流模型计算结果的对比分析这里选用的无粘模型和用y+=30的边界层网格和用Spalart-Allmaras湍流模型的计算结 果分别与试验数据进行比较。比较结果分别如图4所示：图4无粘模型和S
5、palart-Allmaras湍流模型的压力系数与试验数据的比较由图4可知，湍流模型计算所得的压力系数比无粘模型计算所得压力系数与试验数据符合 的更好。3.2三种不同湍流模型的对比分析本文选用 y+=30 的边界层网格，分别采用 Spalart-Allmaras、k-epsilon 和 Scale-Adaptive Simulation这三种湍流模型计算所得的压力系数与试验数据的进行比较。结果如图5所示：图5三种不同湍流模型计算所得压力系数曲线与实验数据曲线的比较根据图5曲线对比分析，三种不同湍流模型的压力系数的计算结果与试验数据的曲线比较吻合，并且 三种湍流模型计算所得压力系数的偏差也很小。
6、3.3四种不同边界层网格的计算结果的对比分析本节选用y+=0.8、y+=10、y+=30和y+=50四种不同的边界层网格，均采用Spalart-Allmaras 湍流模型，分析不同边界层网格计算所得压力系数与试验数据的差别。四种计算结果和试验数据如图6所示：根据图6所示y+=10时与试验的数据差别最大，y+=50与试验数据的吻合度也比较差； y+=0.8和y+=30的计算数据与试验数据的吻合度比较好。所以单独对y+=0.8和y+=30的计 算数据与试验数据进行比较如图7所示：根据图7所示，对y+=0.8和y+=30的计算数据与试验数据吻合的都比较好，分析认为 是因为，y+=0.8和y+=30正
7、好是的边界层与网格线相重合，所以更加准确的计算了边界的流体流动状态。图6四种不同边界层网格的Spalart-Allmaras湍流模型计算所得压力系数与试验数据Q5Ewpcrtra frFlty=1ly=soy=50ErperlimiEnt y=0.8 .4。图7 y+=0.8和y+=30的Spalart-Allmaras湍流模型计算所得压力系数与试验数据word格式-可编辑-感谢下载支持3.4流场求解网格自适应功能计算结果的对比分析本节对于y+=50的网格额外采用流场求解网格自适应功能（基于压力梯度）进行计算。 将计算结果与试验结果和不用网格自适应功能的结果进行对比分析。网格自适应部分设置如
8、图8所示，设置结果如图9所示:图8压力梯度流程自适应的设置Dump U!SaqI 16315 CB11S+ 岫 F3DBS+ 101105311 Mia&S阮18415 c?llsT 岫 Fbe、nodesDinin(i usage:cells, 3ftiflr1d57 Faces,wdesMmptM / 1R9啊 / 38&UF / 19(15岫 /grid size C original /Wils (1S?如IS /faces t:搪/Interriuiptiviq -.图9流场自适应网格的结果对于流场求解自适应选用和y+=50相同的参数设置，计算流场结果与试验数据以及不 用流场求解自适
9、应的计算结果对比分析如图10所示：根据图10所示结果，在采用流场自适应功能（基于压力梯度后）机翼下表面积（即图 形上半部分）与试验数据的吻合更好。图10流场求解自适应功能与非自适应计算结果及试验数据对比分析图四结论本文通过对不同的模型和不同边界层网格的计算，采用控制变量的方法分组比较得 出湍流模型更加符合机翼绕流；对于相同的边界层网格，不同的湍流计算模型的计算结 果是一样的；对于相同的湍流模型，不同的边界层网格计算结果差异比较大，但是都能 反映出大体趋势；流场求解自适应功能的调整使得计算结果更加贴合试验数据。整体来 看，CFD技术能够比较准确的反映出流场的分布，可以帮助研发人员整体的判断流场 分布，并节省了大量的资金和时间，降低研发成本。
