收稿日期: 2017-04-25
网络出版日期: 2019-05-05
基金资助
国家重点基础研究发展计划(973计划资助项目2014CB744100)
Lattice Boltzmann method for simulation of flow with high Reynolds number
Received date: 2017-04-25
Online published: 2019-05-05
针对贴体网格下广义形式的补充插值 LBM方法(generalized interpolation lattice Boltzmann method, GILBM),发展了一种当地时间步法,在改善高雷诺数下翼型绕流模拟中计算稳定性的同时,有效地提高了计算效率. 另外, 通过引入非平衡态外推边界处理,合理解决了 LBGK 模型中处理复杂边界困难的问题. 对 $5.0\times10^4\leqslant Re\leqslant 5.0\times10^6$ 的 NACA0012 翼型绕流进行了数值模拟,并对算法的计算效率和稳定性进行分析, 证明了发展方法的优越性.
关键词： 格子 Boltzmann 方法; 贴体网格; 当地时间步; 高雷诺数; 非平衡态外推
陈彦晓, 李孝伟, 丁珏, 翁培奋 . 高雷诺数流动模拟的 LBM 方法[J]. 上海大学学报(自然科学版), 2019 , 25(2) : 256 -265 . DOI: 10.12066/j.issn.1007-2861.1933
A local time step technique for the generalized interpolation lattice Boltzmann method (GILBM) with the body-fitted grid is developed to improve stability and accelerate convergence of airfoil flow simulation at a high Reynolds number. The difficulty in dealing with complex boundaries in the LBGK model is overcome by introducing the non-equilibrium extrapolation boundary treatment. Flows of NACA0012 airfoil within the range are simulated. Advancement of the method is verified by comparing efficiency and stability of the calculation with traditional GILBM.
| [1] | 乔治$\cdot$埃姆$\cdot$卡尼亚达克斯, 埃里$\cdot$柏斯考克. 微流动------基础与模拟 [M]. 中国科学院过程工程研究所多相反应重点实验室多相复杂系统与多尺度方法课题组, 译. 北京: 化学工业出版社, 2006: 5-37. |
| [2] | Qian Y H ,D'Humieres D,Lallemand P.Lattice BGK models for Navier-Stokes equation[J]. Europhys Lett, 1992,17(6):479-484. |
| [3] | 何雅玲, 王勇, 李庆 . 格子 BOLTZMANN 方法的理论及应用 [M]. 北京: 科学出版社, 2009: 145-151. |
| [4] | Imamura T,Suzuki K,Nakamura T, et al. Flow simulation around an airfoil using lattice Boltzmann method on generalized coordinates [R/OL]. ( 2016- 08- 30)[2017-02-20]. . |
| [5] | Fillipova O, Succi S . Multiscale lattice Boltzmann schemes with turbulence modeling[J]. Journal of Computational Physics, 2001,170(2):812-829. |
| [6] | 李明秀, 陶文铨, 何雅玲 , 等. 格子-BOLTZMANN 方法非均分网格的实施[J]. 工程热物理学报, 2003,24(1):73-75. |
| [7] | Imamura T, Suzuki K, Nakamura T , et al. Flow simulation around an airfoil by lattice Boltzmann method on generalized coordinates[J]. AIAA Journal, 2005,43(9):1968-1973. |
| [8] | He X Y, Doolen G D . Lattice Boltzmann method on curvilinear coordinates system: flow around a circular cylinder[J]. Journal of Computational Physics, 1997,134(2):306-315. |
| [9] | Yu H D, Girimaji S S . Near-field turbulent simulations of rectangular jets using lattice Boltzmann method[J]. Physics of Fluids, 2005,17(12):125106. |
| [10] | Imamura T, Suzuki K, Nakamura T , et al. Acceleration of steady-state lattice Boltzmann simulations on non-uniform grid using local time step method[J]. Journal of Computational Physics, 2005,202(2):645-663. |
| [11] | 应纯同 . 气体输运理论及应用 [M]. 北京: 清华大学出版社, 1990: 45-51. |
| [12] | 郭照立, 郑楚光, 李青 , 等. 流体动力学的格子 Boltzmann 方法 [M]. 武汉: 湖北科学技术出版社, 2002: 54-60. |
| [13] | Im H S ,Zhan G C. Delayed detached eddy simulation of astall flow over NACA0012 airfoil using high order schemes[R/OL].( 2016 -05-31) [2017-02-12]. . 2011-1297. |
| [14] | Gregory N ,O'Reilly C L. Low-speed aerodynamics of MACA0012 aerofoil section, including the effects of upper-surface roughness simulation hoar frost[R].National Physical Laboratory, NPL Aero Report 1308, 1970. |
| [15] | Harris C D . Two-dimensional aerodynamic characteristics of the NACA0012 airfoil in the Langley 8-foot transonic pressure tunnel[M]. Hampton: NASA Technical Memorandum, 1981: 55-76. |
/
| 〈 |
|
〉 |
