Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

doi:10.1007/s10483-017-2194-9

Applied Mathematics and Mechanics (English Edition) ›› 2017, Vol. 38 ›› Issue (5): 707-722.doi: https://doi.org/10.1007/s10483-017-2194-9

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Chenggong LI¹, J. P. Y. MAA²

1. National Engineering Laboratory for Methanol to Olefins, Dalian National Laboratory for Clean Energy, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning Province, China;
2. College of William and Mary, Virginia Institute of Marine Science, Virgina 23062, U.S.A

收稿日期:2016-10-10 修回日期:2016-10-10 出版日期:2017-05-01 发布日期:2017-05-01
通讯作者: Chenggong LI E-mail:chenggongli@dicp.ac.cn

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Chenggong LI¹, J. P. Y. MAA²

1. National Engineering Laboratory for Methanol to Olefins, Dalian National Laboratory for Clean Energy, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning Province, China;
2. College of William and Mary, Virginia Institute of Marine Science, Virgina 23062, U.S.A

Received:2016-10-10 Revised:2016-10-10 Online:2017-05-01 Published:2017-05-01
Contact: Chenggong LI E-mail:chenggongli@dicp.ac.cn

摘要/Abstract

摘要：

Large eddy simulation (LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components (D2Q9) lattice Boltzmann equation (LBE) with multi-relaxation-time (MRT) to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1×10⁷. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit (GPU) with a computing unified device architecture (CUDA) is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.

关键词: neutral equation, distribution of zeroes, estimate, lattice Boltzmann equation (LBE), large eddy simulation (LES), multi-relaxation-time (MRT), two-dimensional nine velocity components (D2Q9), graphic processing unit (GPU), computing unified device architecture (CUDA), Smagorinsky model

Abstract:

Key words: neutral equation, distribution of zeroes, estimate, lattice Boltzmann equation (LBE), large eddy simulation (LES), multi-relaxation-time (MRT), two-dimensional nine velocity components (D2Q9), graphic processing unit (GPU), computing unified device architecture (CUDA), Smagorinsky model

中图分类号:

O357

Chenggong LI, J. P. Y. MAA. Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit[J]. Applied Mathematics and Mechanics (English Edition), 2017, 38(5): 707-722.

参考文献

[1] McNamara, G. R. and Zanetti, G. Use of the Boltzmann equation to simulate lattice-gas automata. Physical Review Letters, 61, 2332-2335(1988)
[2] Chen, H. D., Chen, S. Y., and Matthaeus, W. H. Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method. Physical Review A, 45, 5339-5342(1992)
[3] Qian, Y. H., d'Humières, D., and Lallemand, P. Lattice BGK models for Navier-Stokes equation. Europhysics Letters, 17, 479-484(1992)
[4] Shan, X. W. and Chen, H. D. Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47, 1815-1820(1993)
[5] Dou, Z., Zhou, Z. F., Huang, Y., and Wu, W. Numerical study of non-aqueous phase liquid transport in a single filled fracture by lattice Boltzmann method. Journal of Hydrodynamics, 24, 130-137(2012)
[6] Guo, Z. L. and Zhao, T. S. Lattice Boltzmann model for incompressible flows though porous media. Physical Review E, 66, 036304(2002)
[7] Hou, S. L., Sterling, J., Chen, S. Y., and Doolen, G. D. A lattice Boltzmann subgrid model for high Reynolds number flows. Fields Institute Communications, 6, 151-166(1996)
[8] Yu, H. D., Girimaji, S. S., and Luo, L. S. DNS and LES of decaying isotropic turbulence with and without frame rotation using lattice Boltzmann method. Journal of Computational Physics, 209, 599-616(2005)
[9] Fernandino, M., Beronov, K., and Ytrehus, T. Large eddy simulation of turbulent open duct flow using a lattice Boltzmann approach. Mathematics and Computers in Simulation, 79, 1520-1526(2009)
[10] Krafczyk, M., Tölke, J., and Luo, L. S. Large-eddy simulations with a multiple-relaxation-time LBE model. International Journal of Modern Physics B, 17, 33-39(2003)
[11] Yu, H. D., Luo, L. S., and Girimaji, S. S. LES of turbulent square jet flow using an MRT lattice Boltzmann model. Computers and Fluids, 35, 957-965(2006)
[12] Premnath, K. N. and Pattison, M. J. Generalized lattice Boltzmann equation with forcing term for computation of wall-bounded turbulent flows. Physical Review E, 79, 026703(2009)
[13] Sukop, M. C. and Thorne, D. T. Lattice Boltzmann Modeling:an Introduction for Geoscientists and Engineers, Springer, Berlin (2006)
[14] d'Humières, D. Generalized lattice Boltzmann equations. Rarefied Gas Dynamics:Theory and Simulations (eds. Shizgal, B. D. and Weaver, D. P.), Progress in Astronautics and Aeronautics, American Instiute of Aeronautics and Astronauics, Reston (1992)
[15] Lallemand, P. and Luo, L. S. Theory of the lattice Boltzmann method:dispersion, dissipation, isotropy, Galilean invariance, and stability. Physical Review E, 61, 6546-6562(2000)
[16] d'Humières, D., Ginzburg, I., Krafczyk, M., Lallemand, P., and Luo, L. S. Multiple-relaxation-time lattice Boltzmann models in three dimensions. Philosophical Transactions of the Royal Society A, 360, 437-451(2002)
[17] Chai, Z. H., Shi, B. C., and Zheng, L. Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method. Chinese Physics, 15, 1855- 1863(2006)
[18] Hou, S. L., Chen, S. Y., Doolen, G. D., and Cogley, A. C. Simulation of cavity flow by the lattice Boltzmann method. Journal of Computational Physics, 118, 329-347(1995)
[19] Li, C. G., Maa, J. P. Y., and Kang, H. G. Solving generalized lattice Boltzmann model for 3-D cavity flows using CUDA-GPU. Science China Physics, Mechanics and Astronomy, 55, 1894-1904(2012)
[20] Li, W., Wei, X. M., and Kaufman, A. Implementing lattice Boltzmann computation on graphics hardware. Visual Computer, 19, 444-456(2003)
[21] NVIDIA Corporation. NVIDIA CUDA 4.1 Programming Guide (2011)
[22] Kirk, D. B. and Hwu, W. W. Programming Massively Parallel Processors:a Hands-on Approach, Morgan Kaufmann, Burlington (2010)
[23] Tölke, J. Implementation of lattice Boltzmann kernel using the compute unified device architecture developed by NVIDIA. Computing and Visualization in Science, 13, 29-39(2009)
[24] Tölke, J. and Krafczyk, M. TeraFLOP computing on a desktop PC with GPUs for 3D CFD. International Journal of Computational Fluid Dynamics, 22, 443-456(2008)
[25] Kuznik, F., Obrecht, C., Rusaouen, G., and Roux, J. LBM based flow simulation using GPU computing processor. Computers and Mathematics with Applications, 59, 2380-2392(2010)
[26] Pope, S. Turbulent Flows, Cambridge University Press, New York (2000)
[27] Smagorinsky, J. General circulation experiments with the primitive equations:I, the basic equations. Monthly Weather Review, 91, 99-164(1963)
[28] Ghia, U., Ghia, K. N., and Shin, C. T. High-Re solutions for incompressible flow using the NavierStokes equations and a multigrid method. Journal of Computational Physics, 48, 387-411(1982)
[29] Erturk, E., Corke, T. C., and Gökçöl, C. Numerical solutions of 2-D steady incompressible driven cavity flow at high Reynolds numbers. International Journal for Numerical Methods in Fluids, 48, 747-774(2005)
[30] Du, R., Shi, B. C., and Chen, X. W. Multi-relaxation-time lattice Boltzmann model for incompressible flow. Physical Letters A, 359, 564-572(2006)
[31] Erturk, E. Discussions on driven cavity flow. International Journal for Numerical Methods in Fluids, 60, 275-294(2009)

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	Y. ZHANG, M. VANIERSCHOT. Proper orthogonal decomposition analysis of coherent motions in a turbulent annular jet[J]. Applied Mathematics and Mechanics (English Edition), 2021, 42(9): 1297-1310.
[2]	Ying LIU, Yanping CHEN, Yunqing HUANG, Qingfeng LI. Analysis of a two-grid method for semiconductor device problem[J]. Applied Mathematics and Mechanics (English Edition), 2021, 42(1): 143-158.
[3]	Zhou YU, Hongda ZHANG, Taohong YE, Minming ZHU. Large eddy simulation of turbulent premixed piloted flame using artificial thickened flame model coupled with tabulated chemistry[J]. Applied Mathematics and Mechanics (English Edition), 2018, 39(9): 1277-1294.
[4]	Chao XU, Dongyang SHI, Xin LIAO. A new streamline diffusion finite element method for the generalized Oseen problem[J]. Applied Mathematics and Mechanics (English Edition), 2018, 39(2): 291-304.
[5]	Zhou YU, Hongda ZHANG, Taohong YE, Minming ZHU. Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method[J]. Applied Mathematics and Mechanics (English Edition), 2018, 39(12): 1719-1736.
[6]	Zhendong LUO, Fei TENG. Reduced-order proper orthogonal decomposition extrapolating finite volume element format for two-dimensional hyperbolic equations[J]. Applied Mathematics and Mechanics (English Edition), 2017, 38(2): 289-310.
[7]	Xiaodi WU, Fu CHEN, Huaping LIU. Combined immersed boundary method and multiple-relaxation-time lattice Boltzmann flux solver for numerical simulations of incompressible flows[J]. Applied Mathematics and Mechanics (English Edition), 2017, 38(12): 1679-1696.
[8]	Chao XU, Dongyang SHI, Xin LIAO. Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations[J]. Applied Mathematics and Mechanics (English Edition), 2016, 37(8): 1095-1112.
[9]	Jincun LIU, Hong LI, Yang LIU, Zhichao FANG. Reduced-order finite element method based on POD for fractional Tricomi-type equation[J]. Applied Mathematics and Mechanics (English Edition), 2016, 37(5): 647-658.
[10]	Jikun ZHAO, Shipeng MAO, Weiying ZHENG. Anisotropic adaptive finite element method for magnetohydrodynamic flow at high Hartmann numbers[J]. Applied Mathematics and Mechanics (English Edition), 2016, 37(11): 1479-1500.
[11]	Mengda LIN, Guixiang CUI, Zhaoshun ZHANG. Large eddy simulation of aircraft wake vortex with self-adaptive grid method[J]. Applied Mathematics and Mechanics (English Edition), 2016, 37(10): 1289-1304.
[12]	Hongbo GUAN, Dongyang SHI. P₁-nonconforming triangular finite element method for elliptic and parabolic interface problems[J]. Applied Mathematics and Mechanics (English Edition), 2015, 36(9): 1197-1212.
[13]	Yirang YUAN;Aijie CHENG;Danping YANG;Changfeng LI;Yunxin LIU. Theory and application of numerical simulation method of capillary force enhanced oil production[J]. Applied Mathematics and Mechanics (English Edition), 2015, 36(3): 379-400.
[14]	李想张卫国李正明. Shape analysis and damped oscillatory solutions for a class of nonlinear wave equation with quintic term[J]. Applied Mathematics and Mechanics (English Edition), 2014, 35(1): 117-132.
[15]	方志朝李宏. Numerical solutions to regularized long wave equation based on mixed covolume method[J]. Applied Mathematics and Mechanics (English Edition), 2013, 34(7): 907-920.