双方程k-ω SST湍流模型的显式耦合求解及其在叶轮机械中的应用
收稿日期: 2013-03-12
修回日期: 2013-09-11
网络出版日期: 2013-09-15
基金资助
国家自然科学基金（51076131）
Explicit Coupled Solution of Two-equation k-ω SST Turbulence Model and Its Application in Turbomachinery Flow Simulation
Received date: 2013-03-12
Revised date: 2013-09-11
Online published: 2013-09-15
Supported by
National Natural Science Foundation of China (51076131)
双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法，即通过点隐的方式来处理湍流源项的刚性，并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合，从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性，进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例，确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟，验证了SST模型的精度；进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比，发现三者都能正确地捕捉出口参数分布，且SST与SA模型的模拟结果比较一致；对于该算例，SST模型在总温模拟上更具优势，而BL模型在总压分布上与试验值更加接近。
杨金广 , 吴虎 . 双方程k-ω SST湍流模型的显式耦合求解及其在叶轮机械中的应用[J]. 航空学报, 2014 , 35(1) : 116 -124 . DOI: 10.7527/S1000-6893.2013.0374
The two-equation k-ω shear stress transport (SST) turbulence model is commonly integrated in an implicit coupled way, or in an explicit loosely coupled/decoupled way. An explicit coupled implementation is proposed in the paper. In the present application, a point-implicit method is adopted to treat the stiffness of the turbulence source terms. Combined with hybrid Runge-Kutta time marching and popular accelerating techniques, such as local time stepping, implicit residual smoothing, etc., the turbulence equations can be solved simultaneously with the flow equations. In order to strengthen the robustness of the solver, the turbulent variables are limited. The proposed method is first validated against a DLR 2D cascade test case, which demonstrates the physical validity of the results, and determines the origin of the stiffness. Further, an NASA Rotor 67 test case is used to verify the accuracy of the SST model. The Baldwin-Lomax (BL) algebraic model and the Spalart-Allmaras (SA) one-equation model are also used for comparison. Final results indicate that the SA model and SST model achieve consistently better results, and the SST model has the best accuracy among the three models.
[1] Wilcox D C. Turbulence model for CFD[M]. California: DCW Industries Inc., 1993: 124-128.
[2] Menter F R. Improved two-equation k-omega turbulence models for aerodynamic flows, NASA TM-103975[R]. 1993.
[3] Menter F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 269-289.
[4] Menter F, Ferreira J C, Esch T, et al. The SST turbulence model with improved wall treatment for heat transfer predictions in gas turbines[C]//Proceedings of the International Gas Turbine Congress-IGTC 2003-TS-059, 2003.
[5] Menter F R, Kuntz M, Langtry R. Ten years of industrial experience with the SST turbulence model[J]. Turbulence, Heat and Mass Transfer, 2003, 4: 625-632.
[6] Bardina J E, Huang P G, Coakley T J. Turbulence modeling validation, AIAA-1997-2121[R]. 1997.
[7] Tang S M. Application of SST turbulence model in prediction of flows with strong adverse pressure gradient[J]. Ship & Ocean Engineering, 2008, 37(6): 43-47. (in Chinese) 汤士敏. SST湍流模型在强逆压流模拟中的应用[J]. 航海工程, 2008, 37(6): 43-47.
[8] Yin S, Jin D H, Gui X M, et al. Application and comparison of SST model in numerical simulation of the axial compressors[J]. Journal of Thermal Science, 2010, 19(4): 300-309.
[9] Zhao R Y, Yang H. Numerical simulation of the internal flow field of transonic fan and the study on turbulence model[J]. Compressor Blower & Fan Technology, 2010(6): 7-13.(in Chinese) 赵瑞勇, 杨慧. 跨音风扇内部流场数值模拟及湍流模型研究[J]. 风机技术, 2010(6): 7-13.
[10] Liu F, Zheng X. A staggered finite-volume scheme for solving cascade flow with a two-equation model of turbulence, AIAA-1993-1912[R]. 1993.
[11] Chima R V. A k-ω turbulence model for quasi-three-dimensional turbomachinery flows, NASA TM-107051[R]. 1995.
[12] Liu F, Zheng X. A strongly coupled time-marching method for solving the Navier-Stokes and k-ω turbulence model equations with multigrid[J]. Journal of Computational Physics, 1996, 128(2): 289-300.
[13] Rumsey C. The Menter shear stress transport turbulence model[EB/OL]. (2011-7-10)[2012-12-20]. http://turbmodels.larc.nasa.gov/sst.html.
[14] Yang J, Huang X, Wu H. Multirow inverse method based on adjoint optimization, ASME Paper, GT-2011-46130[R]. 2011.
[15] Jameson A, Schmidt W, Turkel E. Numerical solutions of the Euler equations by finite volume methods using Runge-Kutta time-stepping schemes, AIAA-1981-1259[R]. 1981.
[16] Blazek J. Computational fluid dynamics: principles and applications[M]. Oxford: Elsevier, 2005: 187-188.
[17] Mavriplis D J, Jameson A. Multigrid solution of the Navier-Stokes equations on triangular meshes[J]. AIAA Journal, 1990, 28(8): 1415-1425.
[18] Kunz R F, Lakshminarayana B. Stability of explicit Navier-Stokes procedures using k-ε and k-ε/algebraic Reynolds stress turbulence models[J]. Journal of Computational Physics, 1992, 103(1): 141-159.
/
| 〈 |
|
〉 |
CJA