间断Galerkin方法计算叶型转捩流动

宋寅; 奉凡; 顾春伟

间断Galerkin方法计算叶型转捩流动

清华大学热能工程系热科学与动力工程教育部重点实验室, 北京 100084

基金项目: 国家自然科学重点基金(51136003); 高等学校博士学科点专项科研基金(20090002110052); 北京市自然科学基金(3112015)

计量
- 文章访问数: 1484
- HTML浏览量: 232
- PDF量: 908
- 被引次数: 0
出版历程
- 收稿日期: 2012-06-04
- 刊出日期: 2013-05-28

Transitional flow on turbomachine blade using discontinuous Galerkin method

Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Thermal Engineering, Tsinghua University, Beijing 100084, China

摘要

摘要: 使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.
- 叶型 /
- 转捩 /
- 间断Galerkin方法 /
- 大涡模拟 /
- 转捩模型
Abstract: Large eddy simulation (LES) and a correlation-based transition model were applied with discontinuous Galerkin method (DGM) to investigate transitional flow on turbomachine blade.The methodology was used to simulate the flowfield within T106 turbine cascade and Zierke-Deutsch compressor cascade.Comparison was made between the results of transition model and LES for T106 cascade.The pressure distribution and the locations of separation bubbles calculated by the transition model and LES both agreed well with experimental results. The axial locations of separation bubbles predicted by LES and the transition model are 0.145~0.165m and 0.150~0.168m, respectively.The complex flow details,including the vortex structures in the separation bubble and the dynamic process of laminar separation,vortex shedding and wake diffusion,were simulated by LES,which can help to understand the transition mechanism.As for the transition model,although less flow details were captured,the calculated pressure distribution and the location of separation bubble agreed well with experimental data and LES results.The advantage of the transition model lies in the low computational costs,and it is more suitable for engineering applications.The accuracy of the transition model was further validated by the calculation of Zierke-Deutsch compressor cascade.
- blade /
- transition /
- discontinuous Galerkin method /
- large eddy simulation /
- transition model

HTML

参考文献(24)

[1]	Tain L,Cumpsty N A.Compressor blade leading edges in subsonic compressible flow[J].Journal of Mechanical Engineering Science,2000,214(1):221-242.
[2]	Luo H,Baum J D,Lohner R.A discontinuous Galerkin method based on a Taylor basis for the compressible flows on arbitrary grids[J].Journal of Computational Physics,2008,227(20):8875-8893.
[3]	Feistauer M,Kucera V.On a robust discontinuous Galerkin technique for the solution of compressible flow[J].Journal of Computational Physics,2007,224(1):208-221.
[4]	Reed W H,Hill T R.Triangular mesh methods for the neutron transport equation[R].Los Alamos:Los Alamos Scientific Laboratory,LA-UR-73-479,1973.
[5]	Bassi F,Rebay S.A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations[J].Journal of Computational Physics,1997,131(2):267-279.
[6]	Bassi F,Crivellini A,Rebay S,et al.Discontinuous Galerkin solution of the Reynolds-averaged Navier-Stokes and k-ω turbulence model equations[J].Computers & Fluids,2005,34(4/5):507-540.
[7]	Collis S.The DG/VMS method for unified turbulence simulation[R].AIAA 2002-3124,2002.
[8]	Ramakrishnan S,Collis S.Multiscale modeling for turbulence simulation in complex geometries[R].AIAA 2004-241,2004.
[9]	Zhang X,StanescuD,Naughton J W.Development of a spectral element DNS/LES method for turbulent flow simulations[R].San Diego:ASME,FEDSM 2007-37443,2007.
[10]	Uranga A,Persson P O,Drela 1 M,et al.Implicit large eddy simulation of transition to turbulence at low Reynolds number using a discontinuous Galerkin method[J].International Journal for Numerical Methods in Engineering,2011,87(1/2/3/4/5):232-261.
[11]	苏铭德,张景明.Runge-Kutta间断Galerkin法在求解Navier-Stokes方程中的应用[J].计算物理,1999,16(2):167-176. SU Mingde,ZHANG Jingming.Application of RKDG in solution of N-S equation[J].Chinese Journal of Computation Physics,1999,16(2):167-176.(in Chinese)
[12]	卢义,袁新.隐式间断有限元方法的参数化边界修正[J].工程热物理学报,2009,30(7):1116-1118. LU Yi,YUAN Xin.Parametric boundary modification for implicit discontinuous Galerkin method[J].Journal of Engineering Thermophysics,2009,30(7):1116-1118.(in Chinese)
[13]	郝海兵,杨永,李喜乐.p型多重网格间断Galekin有限元方法研究[J].空气动力学学报,2010,12(6):715-719. HAO Haibing,YANG Yong,LI Xile.The research of p-multigrid solution for high-order discontinuous galerkin finite element method[J].Acta Aerodynamica Sinica,2010,12(6):715-719.(in Chinese)
[14]	贺立新,张来平,张涵信.任意单元间断 Galerkin 有限元计算方法研究[J].空气动力学学报,2007,25(2):157-162. HE Lixin,ZHANG Laiping,ZHANG Hanxin.Discontinuous Galerkin finite element method on 3D arbitrary elements[J].Acta Aerodynamica Sinica,2007,25(2):157-162.(in Chinese)
[15]	于剑,阎超.Navier-Stokes方程间断Galerkin有限元方法研究[J].力学学报,2010,42(5):962-970. YU Jian,YAN Chao.Study on discontinuous Galerkin method for Navier-Stokes equations[J].Chinese Journal of Theoretical and Applied Mechanics,2010,42(5):962-970.(in Chinese)
[16]	Langtry R,Menter F.Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes[J].AIAA Journal,2009,47(12):2894-2906.
[17]	Menter F.Two-equation eddy-viscosity turbulence models for engineering applications[J].AIAA Journal,1994,32(8):1598-1605.
[18]	Boris J P,Grinstein F F,Oran E S,et al.New insights into large eddy simulation[J].Fluid Dynamics Research,1992,10(4/5/6):199-228.
[19]	LI Xuesong,GU Chunwei.An all-speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour[J].Journal of Computational Physics,2008,227(10):5144-5159.
[20]	李雪松.基于可压缩流方法的大涡模拟及其工程化应用[D].北京:中国科学院,2006.LI Xuesong.Large eddy simulation based on the computational method of compressible flows and its application to engineering[D].Beijing:Chinese Academy of Sciences,2006.
[21]	Engber M,Fottner L.The effect of incoming wakes on boundary layer transition of a highly loaded turbine cascade[C]//The Proceeding of AGARD REP 85th Symposium on Loss Mechanisms and Unsteady Flows in Turbomachines.Derby,UK:AGARD,1996:21.
[22]	Raverdy B,Mary I,Sagaut P,et al.High-resolution large-eddy simulation of flow around low-pressure turbine blade[J].AIAA Journal,2003,41(3):390-397.
[23]	Zierke W,Deutsch S.The measurement of boundary layers on a compressor blade in cascade: Volume 1 experimental technique,analysis and results[R].NASA/CR-185118,1989.
[24]	Langtry R,Menter F,Likki S R,et al.A correlation-based transition model using local variables:Part Ⅱ test cases and industrial applications[R].ASME Paper GT2004-53454,2004.