基于格子Boltzmann方法的方通道湍流的大涡模拟

吴宏; 王蛟

基于格子Boltzmann方法的方通道湍流的大涡模拟

吴宏,
王蛟

北京航空航天大学能源与动力工程学院航空发动机气动热力国家级重点实验室,北京 100191

基金项目: 国家自然科学基金(50976008)

计量
- 文章访问数: 1994
- HTML浏览量: 242
- PDF量: 449
- 被引次数: 0
出版历程
- 收稿日期: 2011-03-04
- 修回日期: 2011-07-06
- 刊出日期: 2012-01-28

Large eddy simulation of square duct turbulence based on the lattice Boltzmann method

WU Hong,
WANG Jiao

National Key Laboratory of Science and Technology on Aero-Engine Aero-thermodynamics, School of Jet Propulsion, Beijing University of Aeronautics and Astronautics, Beijing 100191, China

摘要

摘要: 基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.
- 格子Boltzmann方法 /
- 大涡模拟 /
- 方通道湍流 /
- 多松弛时间模型 /
- 切应力改善Smagorinsky模型
Abstract: Large eddy simulation(LES),based on the lattice Boltzmann equation (LBE),was used to study the square duct turbulent flow at a shear Reynolds number equal to 300 based on the hydraulic diameter and the friction velocity.The flow field was simulated by the multiple relaxation time lattice Boltzmann equation.The non-closure term after filtering was modeled by the shear improved Smagorinsky sub-grid stress model.To combine the lattice Boltzmann and large eddy simulation,the subgrid scales turbulent viscosity was related to the subgrid-scale relaxation time.Then the effective relaxation times were computed locally.The turbulent statistics quantities,including the mean streamwise velocity,mean second velocity,the root-mean-square velocity and the instantaneous streamwise voriticity on different sections,were computed.The comparison result,with the DNS and experimental data,shows that the LES based on the LBE can be used to study the bounded turbulent flow.
- lattice Boltzmann method /
- large eddy simulation /
- square duct turbulence /
- multiple relaxation time /
- shear improved Smagorinsky model

HTML

参考文献(19)

[1]	Demuren A O,Rodi W.Calculation of turbulence-driven secondary motion in non-circular ducts[J].Journal of Fluid Mechanics,1984,140:189-222.
[2]	Gavrilakis S.Numerical simulation of low Reynolds number turbulent flow through a straight square duct[J].Journal of Fluid Mechanics,1992,244:101-129.
[3]	Huser A,Biringen S.Direct numerical simulation of turbulent flow in a square duct[J].Journal of Fluid Mechanics,1993,257:65-95.
[4]	Madabhushi R K, Vanka S P.Large eddy simulation of turbulence-driven secondary flow in a square duct[J].Physics of Fluids A:Fluid Dynamics,1991,3(11):2734-2745.
[5]	Kajishima T,Miyake Y. A discussion on eddy viscosity models on the basis of the large eddy simulation of turbulent flow in a square duct[J].Computers & Fluids,1992,21(2):151-161.
[6]	Su M D,Friedrich R.Investigation of fully developed turbulent flow in a straight duct with large eddy simulation[J].Journal of Fluids Engineering,1994,116(4):677-684.
[7]	Aidun C K,Clausen J R. Lattice-Boltzmann method for complex flows[J].Annual Review of Fluid Mechanics,2010,42:439-472.
[8]	Lammers P,Beronov K N,Volkert R,et al.Lattice BGK direct numerical simulation of fully developed turbulence in incompressible plane channel flow[J].Computers & Fluids,2006,35(10):1137-1153.
[9]	Teixeira C M. Incorporating turbulence models into the lattice-Boltzmann method[J].International journal of Modern Physics C,1998,9(8):1159-1175.
[10]	Premnath K N,Pattison M J,Banerjee S.Generalized lattice Boltzmann equation with forcing term for computation of wall-bounded turbulent flows[J].Physical Review E,2009,79(2):026703(19 pages).
[11]	Pattison M J,Premnath K N,Banerjee S.Computation of turbulent flow and secondary motions in a square duct using a forced generalized lattice Boltzmann equation[J].Physical Review E,2009,79(2):026704(13 pages).
[12]	Lévêque E,Toschi F,Shao L,et al.Shear-improved Smagorinsky model for large-eddy simulation of wall-bounded turbulent flows[J].Journal of Fluid Mechanics,2007,570:491-502.
[13]	Germano M,Piomelli U,Moin P,et al.A dynamic subgrid-scale eddy viscosity model[J].Physics of Fluids A:Fluid Dynamics,1991,3(7):1760-1765.
[14]	Lilly D K.A proposed modification of the Germano subgrid-scale closure method[J].Physics of Fluids A:Fluid Dynamics,1992,4(3):633-635.
[15]	HE Xiaoyi,LUO Lishi.Theory of the lattice Boltzmann method: from the Boltzmann equation to the lattice Boltzmann equation[J].Physical Review E,1997,56(6):6811-6817.
[16]	Lallemand P,Luo Lishi.Theory of the lattice Boltzmann method:dispersion,dissipation,isotropy,Galilean invariance,and stability[J].Physical Review E,2000,61(6):6546-6562.
[17]	Premnath K N,Abraham J. Three-dimensional multi-relaxation time (MRT) lattice-Boltzmann models for multiphase flow[J].Journal of Computational Physics,2007,224(2):539-559.
[18]	Guo Z,Zheng C,Shi B.Discrete lattice effects on the forcing term in the lattice Boltzmann method [J].Physical Review E,2002,65(4):046308(6 pages).
[19]	Lam K,Banerjee S.On the condition of streak formation in a bounded turbulent flow [J].Physics of Fluids A:Fluid Dynamics,1992,4(2):306-320.