应用于翼型绕流的线性%2f非线性湍流模式地研究.pdfVIP

第十二届全国计算流体力学会议论文 西安·2004年8月17日～20曰 应用于翼型绕流的线性／非线性 湍流模式的研究 褚洪杰一马晖扬 (中国科学院耐}究生院，北京， 100039) 摘要本文选取了四个线性湍流模式、叫个非线性}呙粘性湍流横式和～个鼬式代数庶_：|了模式对绕翼型的不可压缩分离 流动进彳亍r数值模拟。心计算鲁襻姓豹需要，其中部分模式在壁面酣近勰台了·方程模式。通过与实验结槊的比较， 对翼型在犬攻角情瓤下流动产生分离的气韵特性避静了评体。计算结粱亵蟒，菲线憾搂式能够较簿瓣反驶端流的吾蠢 辩性和曲；棼影响。 冀键词不可压缩分离流动；嚣篷绕流：端流；湍流模式 O善l言 ◆一 现阶段对作战飞机提出了越来越高的机动性及敏捷性的要求，这就使得飞机经常在火攻角和非 定常的状态下飞行。因此，翼型绕流的气动力计算及分祈就显褥尤为重要了。特潮是在大攻兔F出臻 背风面分离的典型流动情形，包含了诸如转捩、流线弯曲、逆压梯度、分离和迹流等具有大餐流动机 琏熬瑷象，其中，绕翼型低速演淀流场戆数馥摸拟，茏其是对失速玫角燃避漉场躲数建摸拟}研究分i ’ 离现象以及分离点附近流动性态，对于航空]二程具有煎要的意义。‘ 随蔫计算流体力学的发展，数俊计算在机翼空气动力学研究中起糟越来越重要的作矧。“对于湍’ 漉的数值≥}算来说，数值馁燕隧蓑冀法的不断发震褥到了稷犬稷菠静捧期，。=诗葵缨蘩‘嚣努耀翟禳丈群 度上依赖于所选择的湍流封闭的计算方法。由于直接数值模拟(DNs)和大涡模拟(LEs)的计算代价过高 }{1，就舅游来看。应琢较为广泛麴还是溃漉模式方法。该方法对雷诺平均N．S方程(萸ANsj中的环封． 闭量雷诺麻力进行合理的模化，从而住方程鲲封闭。湍流模式可以分为积犬袭j一桊是器用B鲫ss{oesq 涡粘性假设的框架，包括零方程模式，一方程模式和两方程檬式：另_誉是雷诺廨尢模式，氢埔霸诺，、 赢力赣运模式(二除矩模式)番{季℃数应力模式。在过去鹣卡年篓，国外学者蠢绕麓獒婺绕潺湍滚数慧 模拟做了大量的工作，如Queutey[2]}tlLesehziner[3]等，所采用的湍流模式为结构较为简龅触线性涡i 静，性模式。然而，对于绕翼型的湍流运动来说，当攻角增大到一定程度，翼型背风面上的池界屡曲于 受翻较强逆压梯发的影响迅速变荐。在靠近溉缘楚，滴流边界艨产生势簿，分离褥斑豸V群隽躐蛩，# 阻力增大，形成一个局部回流区。遮一复杂流动过程中，湍流的备向异性表现较豫，“问流隧的湍流结。 鞠强黧影镌流魂匏菇辫行为。线蛙潺旗牲摸式建立茬嚣毽窝蛙骰凌上，对鼍二失玫强龚翌绕波包含熬流 动机理，如各向异性、曲率敛应、压力梯度等。都不能很好的反映。菲线性涡牿性模式在模化过程中 包含了廊变张量和涡量张最的高阶项，因此可以较好的反映出各向异性、曲率及旋转等流动特征。二二 阶矩模式[4】直接觚K．s方程求得，先取平均德瓢雷诺方程，纛飙N—S穷程中藏去雷诺方禚褥§l器凌z 精方程，推得雷诺应力的微分方程，因此能较好地反映湍流的艇杂流动机理，但是要求解“_金庞火的 煽镦分方程组，计算费时缀大。为了缮到稳定购数擅职，二阶簸摸式对数值方法要求也I#锩苛妻瓠掰， 以在工程应用方面受到了限制。代数应力模式在二阶矩模式萋础上进行了简化，彬据雷诺唐为方程中! 对流项和扩散项的组合与湍动能方程相应组合的某种比例关系，代入压力一应变率关联项和湍流耗散 rt一‘，，o。’|，i， 心用于翼型绕流的线忡／非线性湍流模式的研究 项的模型得到。Pope在此基础上提出了显式求解雷诺府力的思想，建立丁显式代数应力模武，这种 模式和非线性涡粘性模式一样具有较高的工科应tL}I价值。本文研究的目的是评估几种较为先进的湍流 模式在翼型绕流这一外流问题中的应用情况，且这几种模式较雷诺应力模式更存易在工程问题中得到 实现。 本文其