流体力学与飞行力学

新型单方程湍流模型构造及其应用

  • 白俊强 ,
  • 张扬 ,
  • 徐晶磊 ,
  • 华俊
展开
  • 1. 西北工业大学 航空学院, 陕西 西安 710072;
    2. 北京航空航天大学 能源与动力工程学院, 北京 100191;
    3. 中国航空研究院, 北京 100012
白俊强男,博士,教授,博士生导师。主要研究方向:设计空气动力学,非定常空气动力学,工程湍流模拟。Tel:029-88492694E-mail:junqiang@nwpu.edu.cn张扬男,博士研究生。主要研究方向:工程湍流模拟,流动控制。E-mail:iamvip2@163.com徐晶磊男,博士,讲师。主要研究方向:湍流稳定性,湍流模拟。Tel:010-82317403E-mail:xujl@buaa.edu.cn

收稿日期: 2013-08-29

  修回日期: 2013-12-23

  网络出版日期: 2014-01-08

基金资助

国家自然科学基金(11002014)

Construction and Its Application of a New One-equation Turbulence Model

  • BAI Junqiang ,
  • ZHANG Yang ,
  • XU Jinglei ,
  • HUA Jun
Expand
  • 1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China;
    2. School of Energy and Power Engineering, Beihang University, Beijing 100191, China;
    3. Chinese Aeronautics Establishment, Beijing 100012, China

Received date: 2013-08-29

  Revised date: 2013-12-23

  Online published: 2014-01-08

Supported by

National Natural Science Foundation of China

摘要

为了提高雷诺应力本构关系式对于非平衡湍流的预测精度并且兼顾求解效率,发展了一种基于湍动能k的单方程(KDO)湍流模型。其主要思路为:采用平板直接数值模拟(DNS)数据对原始Bradshaw假设进行重新标定,使得当地湍动能和雷诺主应力之比能够根据当地流动条件进行自适应调节;同时,对标准k-ε模型中的湍流耗散率输运方程采用代数形式进行模化,进而形成一种一方程湍流模型。算例结果表明:KDO湍流模型对于对数率能够准确反馈,而在带有激波或部件干扰等流动现象的RAE-2822、ONERA-M6和DLR-F6算例中,KDO湍流模型能够准确控制湍动能的增长和衰减,相比于Spalart-Allmaras和Menter k-ω剪切应力输运(SST)模型,KDO湍流模型的计算结果有了较为明显的改善。

本文引用格式

白俊强 , 张扬 , 徐晶磊 , 华俊 . 新型单方程湍流模型构造及其应用[J]. 航空学报, 2014 , 35(7) : 1804 -1814 . DOI: 10.7527/S1000-6893.2013.0502

Abstract

In order to enhance the predictive precision of the influence of Reynolds stress constitutive relation on non-equilibrium turbulent flow and preserve computational efficiency as well, a turbulent kinetic energy k based one-equation (KDO) turbulence model is proposed in this paper. The basic idea of the turbulence model is to recalibrate the original Bradshaw assumption results using the flat plate direct numerical simulation (DNS) data, which makes the local turbulence kinetic energy adaptively adjustable according to the local flow conditions. Meanwhile, the turbulence dissipation rate equation in the standard k-ε model is modeled by using the algebraic type and then a one-equation turbulence model has been completely constructed. The results of test cases showed that the KDO turbulence model could acceptably reflect the log-law. In the cases RAE-2822, ONERA-M6 and DLR-F6 with shock waves or component interferences, the KDO turbulence model is able to control the augmentations and decay of turbulence kinetic energy. Furthermore, comparing with the Spalart-Allmaras and Menter k-ω shear stress transport (SST) models, the KDO turbulence model evidently improved the calculation results.

参考文献

[1] Hoyas S, Jiménez J. Scaling of the velocity fluctuations in turbulent channels up to Retau=2003[J]. Physics of Fluids, 2006, 18(011702): 1-4.

[2] Spalart P R. Strategies for turbulence modelling and simulations[J]. International Journal of Heat and Fluid Flow, 2000, 21(3): 252-263.

[3] Prandtl L. Bericht uber untersuchungen zur ausgebildeten turbulenz[J]. Zeitschrift Für angewandte Mathematik und Mechanik, 1925, 5(1): 136-139.

[4] Spalart P, Allmaras S. A one-equation turbulence model for aerodynamic flows[J]. La Recherche Aerospatiale, 1994, 1(2): 5-21.

[5] Menter F R. Eddy viscosity transport equations and their relation to the k-ε model[J]. Journal of Fluids Engineering, 1997, 119: 876-884.

[6] Edwards J, Roy C, Blottner F, et al. Development of a one-equation transition/turbulence model[J]. AIAA Journal, 2001, 39(9): 1691-1698.

[7] Wilcox D. Reassessment of the scale-determining equation for advanced turbulence models[J]. AIAA Journal, 1988, 26(11): 1299-1310.

[8] Abid R. Evaluation of two-equation turbulence models for predicting transitional flows[J]. International Journal of Engineering Science, 1993, 31(6): 831-840.

[9] Rotta J. Statistische theorie nicht-homogener turbulenz I and II [J]. Journal of Physics, 1951, 129(6): 547-572.

[10] Launder B, Sharma B. Application of the energy-dissipation model of turbulence to the calculation of flow near spinning disc[J]. Letter Heat Mass Transfer, 1974(1): 131-138.

[11] Patel V C, Rodi W, Scheuerer G. Turbulence models for near-wall and low Reynolds number flows-A review[J]. AIAA Journal, 1985, 23(9): 1308-1319.

[12] Menter F. Improved two-equation k-ω turbulence models for aerodynamic flows, NASA-TM 103975. Washington D.C.: NASA, 1992.

[13] Bradshaw P, Ferriss D, Atwell N. Calculation of boundary-layer using the turbulent energy equation[J]. Journal of Fluid Mechanics, 1967, 28(3): 593-616.

[14] Wen X Q, Liu Y W, Fang L, et al. Improving the capability of k-ω SST turbulence model for predicting stall characteristics of airfoil[J]. Journal of Beijing University of Aeronautics and Astronautics, 2013, 31(8): 1127-1132. (in Chinese) 文晓庆, 柳阳威, 方乐, 等. 提高k-ω SST模型对翼型失速特性的模拟能力[J]. 北京航空航天大学学报, 2013, 31(8): 1127-1132.

[15] Liu J Y. An improved SST turbulence model for hypersonic flows[J]. Acta Aeronautica et Astronautica Snica, 2012, 33(12): 2193-2201. (in Chinese) 刘景源. SST湍流模型在高超声速绕流中的改进[J]. 航空学报, 2012, 33(12): 2193-2201.

[16] She Z S, Chen X, Wu Y, et al. New perspective in statistical modeling of wall-bounded turbulence[J]. Acta Mechanica Sinica, 2010, 26: 847-861.

[17] She Z S, Leveque E. Universal scaling laws in fully developed turbulence[J]. Physics Review Letter, 1997, 72(3): 336-339.

[18] Menter F, Kuntz M, Langtry R. Ten years of industrial experience with the SST turbulence model[J]. Turbulence Heat and Mass Transfer, 2003, 4(1): 625-632.

[19] Spalart P, Jou W, Strelets M, et al. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach//Liu C Q, Liu Z N. Advances in DNS/LES, 1st AFOSR International Conference on DNS/LES. Los Angles: Greyden Press, 1997: 137-147.

[20] Menter F. The scale-adaptive simulation method for unsteady turbulent flow predictions. Part 1: theory and model description[J]. Flow Turbulence Combust, 2010, 85(5): 113-138.

[21] Wilcox D. Turbulence modeling for CFD[M]. California: DCW Industries, Inc., 1994: 90.

[22] Xu J L, Yan C. A one-equation scale-adaptive simulation model[J]. Physics of Gases, 2010, 5(1): 79-82.(in Chinese) 徐晶磊, 阎超. 一个一方程Scale-Adaptive simulation模型的构造[J]. 气体物理, 2010, 5(1): 79-82.

[23] Zhang Y, Bai J Q, Hua J, et al. Improvement and assessment of the SST equation based on karman scale and filter method[J]. Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(3): 442-446. (in Chinese) 张扬, 白俊强, 华俊, 等.基于卡门尺度和滤波方法的SST方程改进和评估[J]. 力学学报, 2013, 45(3): 442-446.

[24] Xu J L, Zhang Y, Bai J Q, et al. A transition model based only on turbulence kinetic equation//Hu H Y. Chinese Congress of Theoretical and Applied Mechanics 2013. Xi'an: Chinese Journal of Theoretical and Applied Mechanics, 2013: 142. (in Chinese) 徐晶磊, 张扬, 白俊强, 等.一个方程驱动转捩的湍动能一方程模型//胡海岩. 中国力学大会文集2013. 西安: 中国力学学会, 2013: 142.

[25] White F. Viscous fluid flow[M]. New York: McGraw-Hill, Inc., 1991: 430.

[26] Spalding D B. A single formula for the "law of the wall" [J]. Journal of Applied Mechanics, 1961, 28(3): 455.

[27] Cook P, McDonald M, Firmin M. Aerofoil RAE-2822 pressure distributions, and boundary layer and wake measurements, AR-138, AGARD. : Experimental Data Base for Computer Program Assessment, 1979.

[28] Schmitt V, Charpin F. Pressure distributions on the ONERA-M6 wing at transonic mach numbers, AR-138, AGARD. : Experimental Data Base for Computer Program Assessment, 1979.

[29] Laflin R, Klausmeyer M, Zickuhr T, et al. Data summary of data from the second AIAA computational fluid dynamics drag prediction workshop[J]. Journal of Aircraft, 2005, 42(5): 1165-1178.

[30] Rumsey C L, Rivers S M, Morrison J H, et al. Study of CFD variation on transport configurations for the second drag prediction workshop[J]. Computers and Fluids, 2005, 34(7): 785-816.

[31] Wang Y T, Wang G X, Zhang Y L. Validation of TRIP 2.0: numerical simulation of DPW II complex configuration[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(1): 34-40. (in Chinese) 王运涛, 王光学, 张玉伦. TRIP 2.0软件的确认: DPW II复杂组合体的数值模拟[J]. 航空学报, 2008, 29(1): 34-40.

文章导航

/