基于高阶WENO格式的旋翼非定常涡流场数值模拟

doi:10.7527/S1000-6893.2015.0281

气动/运动耦合数值模拟研究

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

基于高阶WENO格式的旋翼非定常涡流场数值模拟

印智昭, 招启军, 王博

南京航空航天大学直升机旋翼动力学国家级重点实验室, 南京 210016

收稿日期:2015-08-18 修回日期:2015-10-16 出版日期:2016-08-15 发布日期:2015-11-30
通讯作者: 招启军,Tel.:025-84893753。E-mail:zhaoqijun@nuaa.edu.cn E-mail:zhaoqijun@nuaa.edu.cn
作者简介:印智昭,男,硕士研究生。主要研究方向:直升机计算流体力学、高精度CFD数值格式。E-mail:yinzhizhao@nuaa.edu.cn;招启军,男,博士,教授,博士生导师。主要研究方向:直升机计算流体力学、直升机空气动力学及流动控制、气动噪声、总体设计。Tel:025-84893753。E-mail:zhaoqijun@nuaa.edu.cn;王博,男,博士,讲师。主要研究方向:直升机计算流体力学、直升机空气动力学、新概念旋翼飞行器气动设计等。Tel:025-84893753。E-mail:wangbo@nuaa.edu.cn
基金资助:
国家自然科学基金（11272150）；江苏高校优势学科建设工程基金

Numerical simulations for unsteady vortex flowfield of rotors based on high-order WENO scheme

YIN Zhizhao, ZHAO Qijun, WANG Bo

National Key Laboratory of Science and Technology on Rotorcraft Aeromechanics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Received:2015-08-18 Revised:2015-10-16 Online:2016-08-15 Published:2015-11-30
Supported by:
National Natural Science Foundation of China (11272150);Priority Academic Program Development of Jiangsu Higher Education Institutions

摘要/Abstract

摘要：

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散，发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统，以Navier-Stokes方程作为主控方程，湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度，采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值，选取适当的加权因子，通过多个重构模板的凸组合来构造五阶WENO格式，在交界面附近获得具有五阶精度的状态变量。为了提高计算效率，采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后，运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟，并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下，五阶WENO格式能够捕捉到更精确的涡流场特性，反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。

关键词: 旋翼, 涡流, 加权本质无振荡格式, Navier-Stokes方程, 桨尖涡, 运动嵌套网格

Abstract:

Aimed at improving the computational fluid dynamics (CFD) computation accuracy and reducing the numerical dissipation of the vortex flowfield of helicopter rotor, the highly-accurate CFD method based on fifth-order weighted essentially non-oscillatory (WENO) scheme has been established for numerical simulation on the unsteady flowfield of rotors. In this method, the moving-embedded grid technique has been employed in the numerical simulation, the Navier-Stokes equations have been taken as governing equations, and one equation Spalart-Allmaras turbulence model has been employed. In order to improve the computation accuracy of the rotor vortex formation and evolution, the discretization of convective fluxes and the time marching are conducted by Roe-WENO scheme and the highly efficient implicit LU-SGS (Lower-Upper Symmetric Gauss-Seidel) scheme respectively, and the parallel strategy is employed to accelerate the calculation as well. The computation of the left and right flow quantities is based on the fifth-order WENO scheme, which uses the appropriate weights to combine all the candidate stencils to obtain the fifth-order approximation. Finally aerodynamic characteristics and blade-tip vortex of the several rotors (C-T and UH-60A rotor in hover, C-T rotor in non-lifting forward flight and SA349/2 rotor in forward flight) have been simulated by the present method. Comparing numerical results with experimental data and calculated results using the JST scheme and Roe-MUSCL scheme, it is demonstrated that the present method with only 10%~20% increase of computational consumption has high accuracy and low numerical dissipation in simulating the aerodynamic characteristics and blade-tip vortex.

Key words: rotor, vortex flow, weighted essentially non-oscillatory scheme, Navier-Stokes equations, blade-tip vortex, moving-embedded grid

中图分类号:

印智昭, 招启军, 王博. 基于高阶WENO格式的旋翼非定常涡流场数值模拟[J]. 航空学报, 2016, 37(8): 2552-2564.

YIN Zhizhao, ZHAO Qijun, WANG Bo. Numerical simulations for unsteady vortex flowfield of rotors based on high-order WENO scheme[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2016, 37(8): 2552-2564.

参考文献

[1] LEISHMAN J G. Rotorcraft aeromechanics:Getting through the dip[J]. Journal of the American Helicopter Society, 2010, 55(1):11001-1-11001-24.
[2] SRINIVASAN G R, BAEDER J D, OBAYASHI S, et al. Flowfield of a lifting rotor in hover-A Navier-Stokes simulation[J]. AIAA Journal, 1992, 30(10):2371-2378.
[3] POMIN H, WAGNER S. Navier-Stokes analysis of helicopter rotor aerodynamics in hover and forward flight[J]. Journal of Aircraft, 2002, 39(5):813-821.
[4] 杨爱明, 乔志德. 用运动嵌套网格数值模拟旋翼前飞非定常流场[J]. 空气动力学学报, 2000, 18(4):427-433. YANG A M, QIAO Z D. Numerical simulation of unsteady flow around a helicopter rotor in forward flight using moving overset grids[J]. Acta Aerodynamica Sinca, 2000, 18(4):427-433(in Chinese).
[5] ZHAO Q J, XU G H, ZHAO J G. New hybrid method for predicting the flowfields of helicopter rotors[J]. Journal of Aircraft, 2006, 43(2):372-380.
[6] HARIHARAN N. High order accurate numerical convection of vortices across overset interfaces:AIAA-2005-1263[R]. Reston:AIAA, 2005.
[7] MAVRIPLIS D J, DARMOFAL D, KEYES D, et al. Petaflops opportunities for the NASA fundamental aeronautics program:AIAA-2007-4084[R]. Reston:AIAA, 2007.
[8] 招启军, 徐国华. 使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟[J]. 航空动力学报, 2005, 20(2):186-191. ZHAO Q J, XU G H. Numerical simulations for the flowfield of helicopter rotors in hovering based on high-order upwind flux-difference splitting scheme[J]. Journal of Aerospace Power, 2005, 20(2):186-191(in Chinese).
[9] 徐广, 招启军, 王博, 等. 先进直升机旋翼悬停状态气动性能计算[J]. 航空学报, 2010, 31(9):1723-1732. XU G, ZHAO Q J, WANG B, et al. Prediction on aerodynamic performance of advanced helicopter rotor in hover[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(9):1723-1732(in Chinese).
[10] YU J, YAN C, JIANG Z H. A high resolution low dissipation hybrid scheme for compressible flows[J]. Chinese Journal of Aeronautics, 2011, 24(4):417-424.
[11] 李新亮. 高超声速湍流直接数值模拟技术[J]. 航空学报, 2015, 36(1):147-158. LI X L. Direct numerical simulation techniques for hypersonic turbulent flows[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1):147-158(in Chinese).
[12] XU L, WENG P F. High order accurate and low dissipation method for unsteady compressible viscous flow computation on helicopter rotor in forward flight[J]. Journal of Computational Physics, 2014, 258(2):470-488.
[13] HILGENSTOCK A. A fast method for the elliptic generation of three-dimensional grids with full boundary control[C]//Numerical Grid Generation in Computational Fluid Mechanics' 88. Swansea, Wales:Pineridge Press, 1988:137-146.
[14] MEAKIN R L. A new method for establishing intergrid communication among system of overset grids:AIAA-1991-1586[R]. Reston:AIAA, 1991.
[15] ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43(2):357-372.
[16] HARTEN A, HYMAN J M. Self adjusting grid methods for one-dimensional hyperbolic conservation laws[J]. Journal of Computational Physics, 1983, 50(2):235-269.
[17] 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.
[18] LIU X D, OSHER S, CHAN T. Weighted essentially non-oscillatory schemes[J]. Journal of Computational Physics, 1994, 115(1):200-212.
[19] JIANG G S, SHU C W. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126(1):202-228.
[20] BLAZEK J. Computational fluid dynamics:principles and applications[M]. Netherlands:Elsevier, 2001:204-208.
[21] SPALART P R, ALLMARAS S R. A one-equation turbulence model for aerodynamic flows:AIAA-1992-0439[R]. Reston:AIAA, 1992.
[22] CHANDRA R. Parallel programming in OpenMP[M]. San Francisco:Morgan Kaufmann, 2001:41-88.
[23] CARADONNA F X, TUNG C. Experimental and analytical studies of a model helicopter rotor in hover:NASA-TM-8122[R]. Cailifornia:Army Research and Technology Labs, 1980.
[24] LORBER P F. Aerodynamic results of a pressure-instrumented model rotor test at the DNW[J]. Journal of the American Helicopter Society, 1991, 36(4):66-76.
[25] CARADONNA F X, LAUB G H, TUNG C. An experimental investigation of the parallel blade-vortex interaction:NASA-TM-86005[R]. California:Army Research and Technology Labs., 1984.
[26] HEFFERNAN R M, GAUBERT M. Structural and aerodynamic loads and performance measurements of an SA 349/2 helicopter with an advanced geometry rotor:NASA-TM-88370[R]. California:Army Research and Technology Labs., 1986.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

基于高阶WENO格式的旋翼非定常涡流场数值模拟

Numerical simulations for unsteady vortex flowfield of rotors based on high-order WENO scheme

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	马锦超, 陆洋, 王亮权, 宋奎辉. 基于高阶谐波控制的倾转旋翼近场气动噪声主动控制试验[J]. 航空学报, 2024, 45(9): 528602-528602.
[2]	刘琦, 史勇杰, 胡志远, 徐国华. 共轴刚性旋翼气动及噪声特性的参数影响分析[J]. 航空学报, 2024, 45(9): 528856-528856.
[3]	柳家齐, 陈荣钱, 楼锦华, 韩旭, 吴昊, 尤延铖. 基于深度学习的高速直升机旋翼翼型气动优化设计[J]. 航空学报, 2024, 45(9): 529828-529828.
[4]	王畅, 何龙, 徐栋霞, 唐敏, 马率, 吴希明. 共轴刚性旋翼桨毂流动控制减阻研究[J]. 航空学报, 2024, 45(9): 529084-529084.
[5]	刘昌昊, 曹义华, 梅晓萌, 汪茂胜, 张广林. 高速直升机运输效能评估[J]. 航空学报, 2024, 45(9): 530182-530182.
[6]	张卫国, 唐敏, 武杰, 彭先敏, 章贵川, 聂博文, 王亮权, 李超群. 倾转旋翼机风洞试验综述[J]. 航空学报, 2024, 45(9): 530114-530114.
[7]	王宗辉, 杨云军, 赵弘睿, 王雪晨. 多飞行状态倾转旋翼气动优化设计[J]. 航空学报, 2024, 45(9): 529024-529024.
[8]	曾怡兰, 韩东, 刘壮壮, 周鑫. 高速飞行复合式直升机旋翼受上洗时的驱转特性[J]. 航空学报, 2024, 45(9): 529061-529061.
[9]	韩少强, 宋文萍, 韩忠华, 许建华. 高速共轴刚性旋翼非定常流动高精度数值模拟[J]. 航空学报, 2024, 45(9): 529064-529064.
[10]	王梓旭, 李攀, 鲁可, 朱振华, 陈仁良. 共轴刚性旋翼高速直升机配平策略优化设计[J]. 航空学报, 2024, 45(9): 529069-529069.
[11]	崔壮壮, 原昕, 赵国庆, 井思梦, 招启军. 共轴刚性旋翼高速直升机前飞性能操纵策略影响[J]. 航空学报, 2024, 45(9): 529256-529256.
[12]	孙朋朋, 刘平安, 樊枫, 曾伟. 悬停状态共轴刚性旋翼机身气动干扰特性[J]. 航空学报, 2024, 45(9): 529284-529284.
[13]	杨华, 陈树生, 高正红, 姜权峰, 张伟. 基于贝叶斯框架的旋翼气动力数据融合[J]. 航空学报, 2024, 45(8): 128960-128960.
[14]	邓景辉. 电动垂直起降飞行器的技术现状与发展[J]. 航空学报, 2024, 45(5): 529937-529937.
[15]	李国强, 宋奎辉, 覃晨, 赵光银, 吴霖鑫, 杨永东. 基于后缘小翼的翼型动态失速主动控制试验[J]. 航空学报, 2024, 45(3): 128699-128699.