基于高阶WENO格式的旋翼非定常涡流场数值模拟
收稿日期: 2015-08-18
修回日期: 2015-10-16
网络出版日期: 2015-11-30
基金资助
国家自然科学基金（11272150）；江苏高校优势学科建设工程基金
Numerical simulations for unsteady vortex flowfield of rotors based on high-order WENO scheme
Received date: 2015-08-18
Revised date: 2015-10-16
Online published: 2015-11-30
Supported by
National Natural Science Foundation of China (11272150);Priority Academic Program Development of Jiangsu Higher Education Institutions
为提高直升机旋翼黏性涡流场计算流体力学（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格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。
印智昭 , 招启军 , 王博 . 基于高阶WENO格式的旋翼非定常涡流场数值模拟[J]. 航空学报, 2016 , 37(8) : 2552 -2564 . DOI: 10.7527/S1000-6893.2015.0281
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.
[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.
/
| 〈 |
|
〉 |
CJA