三角翼旋涡破裂流动特性初步探讨.pdfVIP

三角翼旋涡破裂流动特性初步探讨 葛立新 庞铁波 (北京应用物理与计算数学研究所计算物理实验室) 捕要： 本文采用Beam-Warming的隐式近似因子分解法，通过求解三维、非定常、可压缩、层流 Ⅳ一S方程，对760后掠三角翼的旋涡破裂流场进行了数值模拟．考察了旋涡破裂点前后速度、压强 以及涡核轨迹等参变量的变化规律，对旋涡破裂机理作了探讨性的研究。 1 引言 旋涡破裂现象及其控制不仅是流体力学学科中有待于掌握的～个重要课题，同时在 航空航天及民用领域中也有重要的应用意义。在过去的四十年间，已开展了大量的研究 工作[1．3】。其中，对管流中的旋涡破裂现象的研究已取得较大进展。虽然三角翼绕流的 旋涡破裂研究也取得了不少成果[4．8】，但由于三角翼背风区旋涡运动高度复杂，人们对 三角翼上旋涡破裂的认识还不充分，还有待于进一步提高。 本论文的研究目的是探讨三角翼背风区旋涡破裂的机制。由于旋涡破裂发生后，流 动由层流转捩为湍流，采用层流假设对旋涡破裂流场进行数值模拟，不能完全反映实际 的流动现象。但由于计算条件的限制以及目前的湍流模型在反映实际湍流流动时存在不 确定性，而采用层流N-s方程模拟旋涡破裂流场，能够反映旋涡破裂后流场的基本结构， 避免了湍流模型所带来的不确定性，还可简化计算，对旋涡破裂的定性分析及机理性研 究不失为一种好的手段。文献[9]采用Euler方程、层流Ⅳ一S方程以及湍流Ⅳ一S方程， 对三角翼绕流流动进行了数值模拟，并将计算结果同实验测量值进行了对比。结果表明， 采用层流Ⅳ一s方程计算的结果与实验测量的结果最相近。因此，本问的计算模型采用了 层流假设。 计算方法及程序的验证见文献[1l】，这里不再赘述 3计算结果及分析 3．1速度轴向分量的变化特征 为考察旋涡破裂点前后的流动特性，本文计算了来流马赫数M。=0．2，雷诺数 Re=5×104，迎角a=36．5”的三角翼绕流流场。 图3给出了涡核处速度的轴向分量-y=“…肛，沿轴向的变化曲线。在破裂点的远 前方，订由零增加到最大值玎=2．4，说明此区间涡核内的流动呈现加速趋势；在破裂点 ’ 前的一小段距离内，订沿轴向急剧下降；至覃=x／c=0．5(c为弦长)时．玎已降为零， 在i=0．7处，订=-0．5，表明此破裂点下游区间涡核内的流动减速，并出现了回流； 之后“逐渐回升，在距离尾缘u．8倍根弦的尾流区处，酉回升到同来流相近的值。图4显 450 示了U’=“／“。在轴向几个站位点上沿z方向变化的曲线，对这些曲线的分析，可以得到 以下的结论：旋涡破裂前涡核内的u’呈射流形态(图4a)。随着流动向下游发展，u’逐 ● 站位点，’出现负值，如图4c中所示，此时出现了回流区。而且回流区的尺度沿轴向逐 渐扩大(图4d)。图5给出了速度的轴向分量”在旋涡破裂前、破裂点及破裂后几个轴 向站位上的等值线。可以看到：旋涡破裂前，涡心处的“值最大，同时，截面上u的分 布随其与涡心的距离的增加而递减，流动呈射流形态；“只有一个峰值，在涡心处，涡 核内的“呈对称分布(图5a)。旋涡破裂时，涡核内“的分布出现多个峰值：流动不再 对称(图5b)。旋涡破裂后，涡心上的u值变为最小，主涡内u的分布随其与涡心的距 离的增加而递增，流动呈尾流形态(图5c)。 3．2压强在破裂前后的变化 图6是压强在破裂前、破裂点及破裂后几个典型轴向站位上的等值线。在破裂点的 远前方，涡心的压强沿轴向由大变小，涡核内的流动受顺压控制(图6a及图6b)。接近破 裂点时，涡心的压强急剧增加，说明涡核内的流动受逆压控制，而且逆压梯度很大(图6c 及图6d)。该变化正是导致旋涡破裂的原因。图7给出了主涡涡核内的压力P沿涡轴的分 布，由该图也可以看出是逆压梯度B}起了旋涡的破裂。 3．3由涡核位置观察破裂形态