原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
中 国 科 学 (A 辑)
第 30 卷 第 2 期 SCIENCE IN CHINA ( Series A) 2000 年 2 月
可压混合层流动转捩到湍流的
直接数值模拟
① ①②
傅德薰 马延文
( ①中国科学院力学研究所 , LNM 开放实验室 ; ②LHD 开放实验室 北京 100080)
张林波
( 中国科学院计算数学与科学工程计算研究所,LSEC 开放实验室 ,北京 100080)
摘要 采用五阶精度迎风紧致格式和六阶对称紧致格式 ,结合三阶 RK 方法求解
可压缩三维NavierStokes 方程 ,直接数值模拟了时间发展平面混合层流. 给出了对
流 Mach 数 Mc = 0. 8 ,Reynolds 数 Re = 200 时流动从初始扰动为一对相等且相反的斜
波开始 ,经过失稳发展形成 ∧涡、马蹄涡和蘑菇云等拟序结构 ,然后破碎为小尺度、再
破碎为更小尺度的涡结构 ,形成以小尺度运动为主导的流动 ,最终达到湍流.
关键词 可压平面混合层 转捩机制 湍流 紧致格式
混合层流动是研究湍流机理的典型流动 ,也是解决很多实际问题的基础. 如可压缩平面
混合层流的研究是提高超声速冲压发动机燃烧效率的基础. 对不可压平面混合层流已有很多
学者通过实验、分析和直接数值模拟的方法进行了研究. 而对于可压缩平面混合层流的研究
是近年才开始的. 多数文章仅限于数值模拟涡形成的初期阶段. 很多问题尚不清楚.
[1 ]
早些时候 Sandham 等人 采用直接数值模拟的方法 ,研究了时间发展可压平面混合层流.
给出了对流 Mach 数 Mc = 0. 8 时 ,流动从失稳发展而形成 ∧涡的过程 ,并认为 ∧涡之后将产生
涡对并. 以后,Luo 等人[2 ] 的直接数值模拟结果指出 ,可压平面混合流的转捩机制与不可压混
合流不同 ,它是从初始不稳定直接转捩到湍流的. 没有二次失稳 ,没有涡对并产生. 但他们认
为在 Mc = 0. 8 时 ,转捩过程中的拟序结构是对称结构. 有关可压缩混合层流动的转捩机理、
可压缩效应及混合层增长率随对流 Mach 数增长而下降的原因等问题仍不清楚 ,有待进一步研
究.
本文采用文献[3 ] 中提出的高精度差分方法 ,求解三维可压缩 NavierStokes 方程 ,直接数
(
值模拟了时间发展的平面混合层流动. 对于对流 Mach 数 采用文献[ 1 ] 中的对流Mach 数定
( ) ( ) )
义 , 即 Mc = u1 - u2 / c1 + c2 , u1 , u2 和 c1 , c2 为对应上、下流动的来流速度和声速 Mc =
0. 8 , Reynolds 数 Re = 200 的流动 , 给出了流动从初始失稳、转捩到湍流的整个过程. 指出
收稿 收修改稿
( )
国家自然科学基金 批准号 : 及科学技术部“九五”攀登 A 预选项目“流体与空气动力学关键基础问题研究”
资助项目
162
文档评论（0）