NACA0012失速点附近RANSLES数值模拟.pdfVIP

第十六届全国计算流体力学会议 NACA0012失速点附近RANS/LES 数值模拟 杨伟，邓小兵，赵晓慧 （中国空气动力研究与发展中心空气动力学国家重点实验室，四川绵阳 621000) 摘要：众所周知，许多翼型在失速点附近都表现出低频振荡现象，其Strouhal数约为0.02。 为了深入对这一现象的动力学机理和特征进行研究，对NACA0012翼型失速点附近绕流流场 6 进行了RANS/LES混合模拟，马赫数M∞= 0.2 ，雷诺数Re =6 ×10 ，攻角范围13°～20°。把 数值结果和试验验进行了比较。结果表明，低频振荡现象存在于失速点附近一定攻角范围 内，升力系数表现出较大的非定常低频振荡，流动在大范围分离流和附着流之间反复交替。 进一步考察了流场结构，研究分析了分离泡的形成、破裂的过程和及流动振荡机理。 关键词：临界点，低频振荡， RANS/LES混合算法，数值模拟 一 引言 某些情形下，飞行器的气动力会在某些临界点出现突变。采用传统的定常 （或时间平均） 的实验或计算流体力学方法来研究这些问题的时候，通常在临界点得到的气动力曲线是突变 的。然而，在气动力突变的临界点附近，事实上存在着两个稳定解，因而在非定常的实验或 数值模拟中，常常发现流场结构在两个稳定状态之间反复振荡，从而气动力在临界点附近表 现出较强的非定常特征。 Jones[1]很早就指出翼型的升阻力在最大升力附近会出现低频率的剧烈振荡。Zaman等 4 5 人[2]首先在研究低雷诺数条件下（Re在4×10 到4×10 之间）LRN(1)-1007翼型上流动对声 学激励的响应时发现了流动低频振荡现象，Strouhal数约为0.02。随后，Zaman等人[3]对这 种流动低频振荡进行了进一步的实验和计算研究，确认了这种现象的原因在于流体动力学， 而非设备等其他因素。他们指出流动在失速和非失速两种状态下呈拟周期反复切换，升力随 之低频振荡，振幅达50%。他们还发现NACA0012也具有类似现象。 其他研究人员也实验中报道了这种低频振荡现象。Bragg和Khodadoust[4]发现雷诺数为 1.5×106 时，NACA0012翼型上流动有低频振荡现象。随后，Bragg等人[5]又在更大的风洞 6 里对LRN(1)-1007翼型进行了实验，将翼型上流动低频振荡现象严拓到雷诺数1.4×10 。 Reda[6]在雷诺数106条件下对SAND0018/50翼型进行了液晶可视化实验，发现翼型上流动也 1 基金项目：国家自然科学基金资助项目（批准号 ￥邓小兵： dengxb7504@163.com 第十六届全国计算流体力学会议 有类似低频振荡发生。Bragg等人[7,8]进一步用表面流场可视化实验揭示，流动振荡过程中， 翼型前沿附近存在较大层流分离泡。另外，随着攻角、雷诺数的增加，Strouhal会不断变大。 Broeren和Bragg[9，10]用实验研究了翼型失速类型和非定常低频振荡之间的关系，发现薄 翼失速比尾沿失速和前沿失速的流动振荡更加剧烈。他们通过LDV实验手段展示了前沿分离 泡和后沿分离如何发展并最终导致翼型失速的基本过程。Tanaka[11]用PIV技术对NACA0012 失速点附近流动低频振荡进行了实验，发现这种不一定发生在某一个特定的攻角，它实际上 失速点附近一定攻角范围内的非定常流动现象。Rinoie和Takemura[12]也对NACA0012做了类 似的的实验，试图通过深入考察流动的分离-再附过程，来揭示流动低频振荡的机理。他们 指出，翼型上流动低频振荡问题，对预测飞行器失速有至关重要的意义。 在数值模拟研究方面，Mukai等人[13]首次用大涡模拟手段研究了失速点附近翼型上流 动低频振荡问题。他们用标准Smag