高精度数值方法与非线性扰动研究.PDF

第 卷 第 期 空 气 动 力 学 学 报 ， ! # $% ’! (% ’ # 年 月 !#! !$%’()!*+! ,+)+! ， !))* )+ ,-. ’ !))* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） )!/01 20!/ !))* )#1)!+/1)/ 高精度数值方法与非线性扰动研究 2 ! # ! 夏治强 ，袁湘江 ，沈 清 ，涂国华 （ 防化研究院，北京 ； 中国空气动力研究与发展中心，四川 绵阳 ； 2 ’ 2)!!)/ ! ’ *!2))) # ’ 北京空气动力学研究所，北京 2)))3 ） 摘 要：通过数值模拟扰动波在平板边界层中的传播，考察了格式对 41, 波传播过程计算的精确性。本文还研究了 波 波与涡 涡两种典型的非线性扰动模型，结果表明：在自由剪切层中，波 波模型能加快扰动的增长速度，并 5 5 5 激发出了 结构，但仍然经过了线性失稳阶段；在钝体头部的超声速边界层中，涡 涡模型通过涡 涡相互作用， ! 5 5 产生“发卡”结构，进而出现湍流斑。在此过程中，没有线性失稳的迹象。 关键词：紧致格式；“ ”转捩；非线性扰动 67.899 中图分类号： 文献标识码： $!22 ’ # : 处产生湍流斑。湍流斑的持续出现和扩散，最终形成 - 引 言 湍流边界层。“ ”则是在有限幅值的非 67.899 8;9??%; 随着航空和航天事业的高度发展，人们在新型飞 线性扰动作用下，使扰动的演化过程跳过自然转捩的 行器的设计上，已不满足于仅仅对已有飞行器的修修 前两个阶段，直接形成湍流斑，从而迅速进入湍流。 补补，而是希望通过对流动性态变化的物理机制和对 但“ ”型转捩研究中的一个首要问题是，什么形 67.899 流场中拓扑结构发生、发展的物理规律的深入研究提 式的非线性扰动才能激发出这种“逾越”型转捩？换 出新的设计概念。这就给流场计算的数值方法提出 句话说，应该建立什么样的非线性扰动的数学物理模 了更高的要求，不仅需要较高的精度，还需要较高的