基于小波奇异分析的流场计算方法及应用

王保国; 吴俊宏; 朱俊强

基于小波奇异分析的流场计算方法及应用

1.
北京理工大学宇航学院,北京 100081
2.
中国科学院工程热物理研究所,北京 100190

基金项目: 国家自然科学基金(50376004); 高等学校博士学科点专项基金(20030007028)

计量
- 文章访问数: 1636
- HTML浏览量: 171
- PDF量: 423
- 被引次数: 0
出版历程
- 收稿日期: 2010-06-04
- 修回日期: 2010-09-08
- 刊出日期: 2010-12-28

Method based on wavelet sigularity analysis for complicated flow and its application

1.
School of Aerospace Engineering, Beijing Institute of Technology,Beijing 100081,China
2.
Institute of Engineering Thermophysics, Chinese Academy of Sciences,Beijing 100190,China

摘要

摘要: 首次将三维小波以及二维小波奇异性分析的思想引进三维以及二维复杂流场的数值计算,发展了一种高效率、高精度、高分辨率的方法.该算法的核心是获取流场中物理量在不同空间点的Hlder指数α,而该指数α的获取又依赖于小波变换以及高维(即二维或三维)小波分析技术.在三维与二维欧氏空间中,为进行小波多分辨分析,需要在尺度空间与小波空间分别引进尺度基与小波基.对二维问题,尺度基与小波基的基底要由1个尺度函数与3个小波函数组成,而三维时要由1个尺度函数和7个小波函数组成.借助于小波奇异分析所找到的流场中奇异点区域与光滑区域,便可分别选用高分辨率、高精度的WENO(weighted essential non-oscillatory)格式与高精度中心差分格式进行流场的离散求解.一系列二维(即,①二维前台阶问题的Euler流;②二维双马赫反射的Euler流; ③著名的二维Riemann问题;④跨声速RAE2822翼型二维绕流; ⑤跨声速VKI-LS 59二维涡轮叶栅绕流)与三维(即,⑥跨声速轴流压气机转子NASA(National Aeronautics and Space Administration) Rotor 37三维黏性绕流;⑦跨声速风扇转子NASA Rotor 67三维黏性绕流)算例表明:该方法的计算效率比传统的WENO格式有较大的提高,大部分典型算例能够提高3~5倍,而且可以获取复杂流场中高分辨率的激波结构以及涡系分布,可以得到与有关实验数据较为吻合的数值结果.
- 三维小波 /
- 二维小波 /
- 奇异分析 /
- 多分辨分析 /
- 高精度WENO(weighted essential non-oscillatory)格式
Abstract: Three-dimensional and two-dimensional wavelet sigularity analysis was introduced in complicated flow and a high-order,high-resolution numerical method was developed.The core of this new method is the analysis of Hlder exponent α of complicated flow.The Hlder exponent α depends on wavelet transform and multi-dimensional wavelet analysis.In order to complete the multi-resolution analysis in two-dimensional and three-dimensional Euclidean space,scale basis and wavelet basis have been introduced in scale space and wavelet space.For two-dimensional and three-dimensional cases,scale basis and wavelet basis were constructed by one scale function and by three wavelet functions and seven wavelet functions respectively.High-resolution,high-order WENO (weighted essential non-oscillatory) scheme and high-order central difference scheme can be chosen adaptively according to the results of wavelet sigularity analysis.A series of examples in two-dimensional ( (1) forward facing step problem, (2) double Mach reflection of a strong shock problem, (3) two-dimensional Riemann problem, (4) two-dimensional transonic flow past RAE2822 airfoil problem, (5) two-dimensional transonic VKI-LS 59 turbine cascade flow problem) and three-dimensional ( (6) transonic axial compressor,NASA (National Aeronautics and Space Administration) Rotor 37 three-dimensional viscous flow and (7) NASA Rotor 67 three-dimensional viscous flow) were computed.The computation results have been compared with the experiment results in related references and showed that this new numerical method can not only improve the computation efficiency significantly(the efficiency can be increased by 3~5 times for most examples),but also capture the shock with high-resolution and vortices with an acceptable order.
- three-dimensional wavelet /
- two-dimensional wavelet /
- sigularity analysis /

HTML

参考文献(24)

[1]	Jaffard S,Meyer Y,Ryan R D.Wavelets:tools for science & technology[M].Philadelphia:SIAM,2001.
[2]	Mallat S G.A Wavelet tour of image processing[M].New York:Academic Press,1998.
[3]	Harten A,Engquist B,Osher S,et al.Uniformly high order accurate essentially non-oscillatory schemes,Ⅲ[J].Journal of Computational Physics,1987,71(2):231-303.
[4]	王保国,黄虹宾.叶轮机械跨声速及亚声速流场的计算方法[M].北京:国防工业出版社,2000.
[5]	Shu C W.Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws .NASA/CR-97-206253,ICASE Report No.97-65,1997.
[6]	王保国,刘淑艳,黄伟光.气体动力学[M].北京:国防科工委所属五所院校(北京理工大学、北京航空航天大学、西北工业大学、哈尔滨工程大学、哈尔滨工业大学)出版社,2005.
[7]	Farge M,Hunt J C R,Vassilicos J C.Wavelet,fractals,and Fourier transforms[M].Oxford:The Clarendon Press,1993.
[8]	Antoniadis A,Oppenheim G.Wavelets and statistics[M].Lecture Notes in Statistics 103.New York:Springer-Verlag,1995.
[9]	WANG Baoguo,WU Junhong.Construction of daubechies wavelet and its application in shock capturing .GyeongJu:9th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows (ISAIF) ,2009-3D-1,2009.
[10]	WU Junhong,WANG Baoguo.A new numerical method based on wavelet detection and its application .GyeongJu:9th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows (ISAIF),2009-2D-3,2009.
[11]	王保国,卞荫贵.转动坐标系中跨声速欧拉流的有限体积-TVD格式[J].空气动力学学报,1992,10(4):472-481. WANG Baoguo,BIAN Yingui.Finite volume-TVD scheme for 3-D Euler transonic flows computations in rotating curvilinear coordinates[J].Acta Aerodynamica Sinica,1992,10(4):472-481.(in Chinese)
[12]	WANG Baoguo,CHEN Naixing.An improved SIP scheme for numerical solutions of transonic stream-function equation[J].International Journal for Numerical Methods in Fluids,1990,10(5):591-602.
[13]	王保国,沈孟育.高速粘性内流的高分辨率高精度迎风型杂交格式[J].空气动力学学报,1995,13(4):365-373. WANG Baoguo,SHEN Mengyu.High resolution and high order accurate upwind TVD algorithm of high speed viscous internal flows[J].Acta Aerodynamica Sinica,1995,13(4):365-373.(in Chinese)
[14]	王保国,刘淑艳,张雅,等.双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用[J].航空动力学报,2005,20(4):534-539. WANG Baoguo,LIU Shuyan,ZHANG Ya,et al.Pseudo- time method of weighted ENO-strongly compact high resolution schemes and its application to turbomachinery unsteady flow[J].Journal of Aerospace Power,2005,20(4):534-539.(in Chinese)
[15]	王保国,李翔.多工况下高超声速飞行器再入时流场的计算[J].西安交通大学学报,2010,44(1):71-76. WANG Baoguo,LI Xiang.CFD prediction of hypersonic vehicle reentry flowfield of multiple cases[J].Journal of Xian Jiaotong University,2010,44(1):71-76.(in Chinese)
[16]	王保国,李翔,黄伟光.激波后高温高速流场中的传热特性研究[J].航空动力学报,2010,25(5):963-980. WANG Baoguo,LI Xiang,HUANG Weiguang.Analytical study of heat transfer of high-temperature and high-velocity flow field behind shock waves[J].Journal of Aerospace Power,2010,25(5):963-980.(in Chinese)
[17]	Woodward P R,Colella P.The numerical simulation of two-dimensional fluid flow with strong shocks[J].Journal of Computational Physics,1984,54(1):115-173.
[18]	Lax P D,Liu X D.Solution of two-dimensional Riemann problems of gas dynamics by positive schemes[J].SIAM J Scientific Computing,1998,19(2):319-340.
[19]	Tatum K E,Slater J H.The validation archive of the NPARC alliance .AIAA Paper 99-0747,1999.
[20]	Deconinck H,Hirsch C H.Finite element methods for transonic blade-to-blade calculation in turbomachines[J].ASME Journal of Engineering for Power,1981,103:665-677.
[21]	Reid L,Moore R D.Performance of a single-stage axial-flow transonic compressor with rotor and stator aspect ratios of 1.19 and 1.26,respectively,and with design pressure ratio of 1.82 .NASA/TP-1338,1978.
[22]	CHEN Naixing.Aerothermodynamics of turbomachinery:analysis and design[M].Singapore:John Wiley & Sons,2010.
[23]	陈懋章.粘性流体动力学基础[M].北京:高等教育出版社,2002:282-312.
[24]	Strazisar A J,Wood J R,Hathaway M D,et al.Laser anemometer measurements in a transonic axial-flow fan rotor .NASA/TP-2879,1989.