细长体绕流流态及其侧向力分布.pdfVIP

首届全国航空航天领域中的力学问题学术研讨会 细长体绕流流态及其侧向力分布 王刚 梁新刚 (清华大学航天航空系，北京100084) 摘 要：通过求解不可压层流NS方程模拟r细长体05到904攻角的绕流过程。绕体流态及其侧向力分布随 攻角变化经历6种连续的状态，即附着流态、对称涡流态、定常非对称二涡流态、定常非对称多重涡流态、过 渡流态和卡门涡街(或非定常尾迹)。绕体涡系的运动足轴向运动、横向运动和展向运动的综合过程；涡系的 展向运动与侧向力分布特性密切相关。完整的截面对称涡流态由主涡对、二次涡对和tertiary涡对组成。每个 主涡在脱体之前都诱导出包括二次涡和tertiary涡在内的次生涡结构，而且不同涡之间有涡量传递和相互合并 的现象。过渡流态是流动与侧向力特征既不同于定常非对称、也不同于完全非定常的一种流动状态，它存在于 一个较宽的攻角范围内。流动分离、非对称以及非定常都是首先起始于后体，然后随攻角增大沿轴向上游推进。 关键词：细长体；非对称涡；数值模拟：大攻角空气动力学 0引言 告(Q)+告(F一只)+啬(G—Gv) +喜(日一H，，)=0 细长体有攻角绕流是一个典型的三维流 ub、 7 动问题，由于它在航空、航天、航海等领域的 式中Q=(，1)(p，pu，gv，pw，e)。，F、G和H 重要应用，所以该问题也是倍受国内外学者关 是无粘矢通量项，Fv、G，和H。是粘性通量项； 注的问题之一nJ。在影响绕体流态及其气动力 J是Jacobiana(考，叩，亭)／a(工，Y，z)，‰v和w 的诸多因素中，攻角是～个重要的主控因素。 分别是』、Y和z方向的Cartesian速度分量，p 攻角不同，绕过细长体的流动状态就不同．因 为空气密度，P为压力。e为总能，且 而导致不同的气动力特性。由于对圆柱绕流问 题的认识臼趋成熟，早期人们就借助于二维圆 用三阶精度格式离散，求解方法采用有限体积 柱绕流对三维旋成体横截面内的流动特性进 法。更为详细的过程请参考文献【5】、【6】和【7】。 行了研究”】。对照圆柱绕流流态随Re的变化 物理模型为尖头拱形头部后接圆柱体，分形面 过程，绕细长体流态也依据攻角变化而划分为 为切形过渡，其中头部长细比为LN／D=3，后 附着流、对称涡、非对称涡、亚临界Re下的 体长细比LA／D=9．14。计算域网格总数为100 万以上。实验表明，在头部设置一定的微扰动 卡门涡街(或超临界Re时的随机尾迹)等4 种流态”】。这种简单的比拟与划分易于我们理 可以得到确定性的流场结构口】。数值模拟中首 解流动现象，但是却忽视了细长体绕流的三维 先在几何模型头部顶端方位角0=130。的位置 特性。本文作者等通过实验对细长体亚临界 设置了长A／／D=I．8x10-3、宽Ab／D=3xl旷、高 Re流动进行了探讨，结合测压与流动显示对 Ah／D=-3x10-4的微型边条作为小扰动。雷诺数 绕体流态进行了分区分析【4】。本文通过数值模 为ReD=V**D／v=I．Oxl05，此雷诺数下细长体绕 拟．进一步研究了细长体在各个攻角范围的涡 流位于亚临界分离状态“…。 流动问题。文中给出了不同攻角范围典型的涡 结构。 1数值方法 图1物理模型 本文采用了由Thomas、Krish和 2结果及分析 Anderson[5]发展的CFL3D方法，该方法后来 被Hsieh和K．C．Wang[6]进一步应用验证，他 通过对22个攻角状态的模拟表明，随着 们的模拟结果与其实验结果符合得很好。控制