基于预处理HLLEW格式的全速域数值算法

计算物理 ›› 2016, Vol. 33 ›› Issue (3): 273-282.

基于预处理HLLEW格式的全速域数值算法

刘中玉, 张明锋, 郑冠男, 杨国伟

中国科学院力学研究所流固耦合系统力学重点实验室, 北京 100190

收稿日期:2015-03-25 修回日期:2015-05-28 出版日期:2016-05-25 发布日期:2016-05-25
通讯作者: 杨国伟,研究员,E-mail:gwyang@imech.ac.cn
作者简介:刘中玉(1988-),男,博士生,主要从事计算流体力学研究,E-mail:liuzhongyu@imech.ac.cn
基金资助:
国家高技术研究发展计划(2014AA110501)资助项目

Preconditioning HLLEW Scheme for Flows at All Mach Numbers

LIU Zhongyu, ZHANG Mingfeng, ZHENG Guannan, YANG Guowei

Key Laboratory for Mechanics in Fluid Solid Coupling Systems, Institute of Mechanics, CAS, Beijing 100190, China

Received:2015-03-25 Revised:2015-05-28 Online:2016-05-25 Published:2016-05-25

摘要/Abstract

摘要： 基于HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)格式引入预处理技术发展适合求解全速域流场的三维Navier-Stokes求解器.引入低速预处理技术,重新构造HLLEW格式的耗散项,给出预处理后的HLLEW格式,并根据预处理后的雅克比矩阵构造相应的隐式时间推进方程.利用预处理方法求解NACA 4412低速不可压流动与RAE 2822跨声速可压缩流动,并与实验结果及原有方法的计算结果对比.结果表明:预处理HLLEW格式不仅提高低速不可压缩流动的计算效率和精度,也保持了对可压缩流动的处理能力,是一种适用于全速域流场数值模拟的有效方法.

关键词: 预处理, 全速域流动, 不可压缩流动, HLLEW格式, 隐式时间推进

Abstract: Based on HLLEW (Harten-Lax-Van Leer-Einfeldt-Wada) scheme, low speed preconditioning technology is introduced to develop a three-dimensional Navier-Stokes solver for flows at all Mach numbers. Low speed preconditioning techniques is introduced to reconstruct dissipative term in HLLEW scheme and preconditioning HLLEW scheme is proposed. Implicit time-marching method is constructed based on preconditioning Jacobian Matrix. Results of NACA 4412 incompressible flow and RAE 2822 transonic flow with preconditioning HLLEW scheme are compared with results by original method and experimental data. It shows that preconditioning HLLEW method improves accuracy and convergence rate for low speed flow. It can be applied for flows at all Mach numbers.

Key words: preconditioning, flow at all mach number, incompressible flow, HLLEW scheme, implicit time-marching

中图分类号:

O241.3

刘中玉, 张明锋, 郑冠男, 杨国伟. 基于预处理HLLEW格式的全速域数值算法[J]. 计算物理, 2016, 33(3): 273-282.

LIU Zhongyu, ZHANG Mingfeng, ZHENG Guannan, YANG Guowei. Preconditioning HLLEW Scheme for Flows at All Mach Numbers[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(3): 273-282.

参考文献

[1] ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43(2):357-372.
[2] LIOU M S, STEFFEN C J. A new flux splitting scheme[J]. Journal of Computational Physics, 1993, 107(1):23-39.
[3] CHANG C H, LIOU M S. A robust and accurate approach to computing compressible multiphase flow:Stratified flow model and AUSM(+)-up scheme[J]. Journal of Computational Physics, 2007, 225(1):840-873.
[4] LIOU M S. A sequel to AUSM:AUSM(+)[J]. Journal of Computational Physics, 1996, 129(2):364-382.
[5] LIOU M S. A sequel to AUSM, Part Ⅱ:AUSM(+)-up for all speeds[J]. Journal of Computational Physics, 2006, 214(1):137-170.
[6] ZHA G, SHEN Y, WANG B. An improved low diffusion E-CUSP upwind scheme[J]. Computers & Fluids, 2011, 48(1):214-220.
[7] SHEN Y, ZHA G. Low diffusion E-CUSP scheme with implicit high order WENO scheme for preconditioned Navier-Stokes equations[J]. Computers & Fluids, 2012, 55:13-23.
[8] HARTEN A, LAX P D, LEER B V. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws[J]. SIAM Review, 1983, 25(1):35-61.
[9] OBAYASHI S, GURUSWAMY G P. Convergence acceleration of a Navier-Stokes solver for efficient static aeroelastic computations[J]. AIAA Journal, 1995, 33(6):1134-1141.
[10] EINFELDT B. On Godunov-type methods for gas dynamics[J]. SIAM Journal on Numerical Analysis, 1988, 25(2):294-318.
[11] EINFELDT B, MUNZ C, ROE P L, et al. On Godunov-type methods near low densities[J]. Journal of Computational Physics, 1991, 92(2):273-295.
[12] OBAYASHI S, WADA Y. Practical formulation of a positively conservative scheme[J]. AIAA Journal, 1994, 32(5):1093-1095.
[13] ZHENG G, YANG G, QIAN W. Flutter analyses of complete aircraft based on hybrid grids and parallel computing[J]. Science China Technological Sciences, 2013, 56(2):398-404.
[14] YANG G, OBAYASHI S, NAKAMINCHI J. Aileron flutter calculation for a supersonic fuselage-wing configuration[C]. 23rd International Congress of Aeronautical Science (ICAS), Toronto, 2002.
[15] FUKUSHIMA Y, MISAKA T, OBAYASHI S, et al. CFD-CAA coupled computation of fan noise propagation from engine nacelle based on Cartesian mesh method[C]. 9th AIAA/CEAS Aeroacoustics Conference, Berlin, 2013.
[16] GUILLARD H, VIOZAT C. On the behaviour of upwind schemes in the low Mach number limit[J]. Computers & Fluids, 1999, 28(1):63-86.
[17] CHORIN A J. A numerical method for solving incompressible viscous flow problems[J]. Journal of Computational Physics, 1967, 2(1):12-26.
[18] TURKEL E. Preconditioned methods for solving the incompressible and low speed compressible equations[J]. Journal of Computational Physics, 1987, 72(2):277-298.
[19] ERIKSSON L. A preconditioned Navier-Stokes solver for low Mach number flows[C]. 3rd ECCOMAS Computational Fluid Dynamics Conference, Paris, 1996.
[20] CHOI Y, MERKLE C L. The application of preconditioning in viscous flows[J]. Journal of Computational Physics, 1993, 105(2):207-223.
[21] TURKEL E. Preconditioned methods for solving the incompressible and low speed compressible equations[J]. Journal of Computational Physics, 1987, 72(2):277-298.
[22] WEISS J M, SMITH W A. Preconditioning applied to variable and constant density flows[J]. AIAA Journal, 1995, 33(11):2050-2057.
[23] COLES D, WADCOCK A J. Flying-hot-wire study of flow past an NACA 4412 airfoil at maximum lift[J]. AIAA Journal, 1979, 17(4):321-329.