基于超声速边界层转捩的WENO格式对比分析

刘朝阳; 王振国; 孙明波; 汪洪波

doi:10.13224/j.cnki.jasp.2015.09.018

基于超声速边界层转捩的WENO格式对比分析

doi: 10.13224/j.cnki.jasp.2015.09.018

国防科学技术大学航天科学与工程学院, 长沙 410073

基金项目:

霍英东教育基金高等院校青年教师基金(131055)

详细信息

作者简介:
刘朝阳(1990-),男,吉林公主岭人,硕士生,主要从事高超声速推进技术的研究.

中图分类号: V435.12
计量
- 文章访问数: 1236
- HTML浏览量: 69
- PDF量: 586
- 被引次数: 0
出版历程
- 收稿日期: 2014-03-17
- 刊出日期: 2015-09-28

Comparison and analysis of WENO schemes base on supersonic boundary layer transition

College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China

摘要

摘要: 基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.
- 边界层转捩 /
- 大涡模拟 /
- WENO格式 /
- 谐波扰动 /
- 湍流
Abstract: Large eddy simulations of the transition process in a supersonic flat-plate boundary layer at free-stream Mach number of 4.5 and Reynolds number of 10000 were carried out using high precision weighted essentially non-oscillation (WENO) schemes. The inviscid flux was discretized by fifth-order, seventh-order, ninth-order WENO schemes separately and the viscous flux was discretized by fourth-order center difference scheme. The temporal discretization was made by means of third-order Runge-Kutta with total variation diminishing(TVD). By imposing a pair of the most unstable oblique first-mode disturbances on the inflow boundary, the transition progress of the flow from laminar, transition process to full turbulence was compared using three WENO schemes. Results indicate that the fifth-order WENO scheme is more dissipative than the seventh-order and ninth-order schemes, so it's difficult to capture eddy and pulse in the flow field. The dissipation of ninth-order WENO scheme is small, enabling to capture smaller eddy and high frequency. To study a problem with high frequency fluctuation or capture fine eddy, it's advisable to use the scheme higher than seventh-order.
- boundary layer transition /
- large eddy simulation /
- weighted essentially non-oscillation scheme /
- oblique disturbance /
- turbulence

HTML

参考文献(18)

[1]	Gerolymos G A,Sénéchal D,Vallet I.Very-high-order WENO schemes[J].Journal of Computational Physics,2009,228(23):8481-8524.
[2]	Yee H C,Sjogreen B,Hadjadj A.Comparative study of three high order schemes for LES of temporally evolving mixing layers[R].Nashville:50th AIAA Aerospace Sciences Meeting,2012.
[3]	ZHANG Shuhai,SHU Chiwang.A new smoothness indicator for the WENO schemes and its effect on the convergence to steady state solutions[J].Journal of Scientific Computing,2007,31(1/2):273-305.
[4]	QIU Jingmei,SHU Chiwang,LIU Jiren,et al.A WENO algorithm for the growth of ionized regions at the reionization epoch[J].New Astronomy,2008,13(1):1-11.
[5]	XIONG Tao,SHU Chiwang,ZHANG Mengping.WENO scheme with subcell resolution for computing nonconservative Euler equations with applications to one-dimensional compressible two-medium flows[J].Journal of Scientific Computing,2012,53(1):222-247.
[6]	Toh H T.Large eddy simulation of supersonic twin-jet impingement using a fifth-order WENO scheme[D].Blacksburg,Virginia:Faculty of Virginia Polytechnic Institute and State University,2003.
[7]	Choi J Y,Yang V.Numerical study on the stability of the supersonic turbulent lifted flame using high resolution schemes[R].San Jose,CA:49th AIAA/ASME/SAE/ASEE Joint Propulsion Conference,2013.
[8]	孙明波.超声速来流稳焰凹腔的流动及火焰稳定机制研究[D].长沙:国防科学技术大学,2008. SUN Mingbo.Flow andflameholding mechanisms of cavity flameholders in supersonic flows[D].Changsha:National University of Defense Technology,2008.(in Chinese).
[9]	董平,黄洪雁,冯国泰.带有压力梯度的平板边界层转捩数值模拟[J].航空动力学报,2007,22(11):1909-1904. DONG Ping,HUANG Hongyan,FENG Guotai.Numerical simulation of boundary layer transition under pressure gradient[J].Journal of Aerospace Power,2007,22(11):1909-1914.(in Chinese)
[10]	潘宏禄,马汉东,王强.超声速边界层转捩拟序结构大涡模拟[J].宇航学报,2006,27(3):498-502. PAN Honglu,MA Hanwen,WANG Qiang.Large eddy simulation of transition coherent strucures in a supersonic boundary layer[J].Journal of Astronautics,2006,27(3):498-502.(in Chinese)
[11]	WANG Suozhu,GAO Zhenxun,LEE Chunhian,et al.Large eddy simulation of supersonic boundary layer transition over a flat-plate based on the spatial mode[J].Advances in Mechanical Engineering,2014,6:1-10.
[12]	黄章峰.超声速平板边界层从层流到湍流的转捩机理及湍流特性[D].天津:天津大学,2006. HUANG Zhangfeng.The mechanism of laminar-turbulent transition and the characteristics of turbulence in a supersonic boundary layer on a flat plate[D].Tianjin:Tianjin University,2006.(in Chinese).
[13]	SHEN Jun,SHU Chiwang,ZHANG Mengping A high order WENO scheme for a hierarchical size-structured population model[J].Journal of Scientific Computing,2007,33(3):279-291.
[14]	Dinshaw,Balsara S,SHU Chiwang.Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy[J].Journal of Computational Physics,2000,160(2):405-452.
[15]	Barth J,Deconinck H.High-order methods for computational physics[M].Berlin:Lecture Notes in Computational Science and Engineering,1999.
[16]	马汉东,潘宏禄,王强.超声速平板边界层斜波失稳转披过程研究[J].力学学报,2007,39(2):153-157. MA Handong,PAN Honglu,WANG Qiang.The study of flow transition process induced by oblique wave instability in a supersonic flat-plate boundary layer[J].Chinese Journal of Thoeretical and Applied Mechanics,2007,39(2):153-157.(in Chinese).
[17]	Terzi D A,Sandberg R D,Fasel H F.Identification of large coherent structures in supersonic axisymmetric wakes[J].Computers and Fluids,2009,38:1638-1650.
[18]	曹伟,黄章峰,周恒.超音速平板边界层转捩中层流突变为湍流的机理研究[J].应用数学和力学,2006,27(4):379-386. CAO Wei,HUANG Zhangfeng,ZHOU Heng.Study of themechanism of breakdown in laminar turbulent transition of a supersonic boundary layer on a flat plate[J].Applied Mathematics and Mechanics,2006,27(4):379-386.(in Chinese)