x

## 留言板

 引用本文:
 Citation:

## Three-dimensional 12-velocity multiple-relaxation-time lattice Boltzmann model of incompressible flows

Hu Jia-Yi, Zhang Wen-Huan, Chai Zhen-Hua, Shi Bao-Chang, Wang Yi-Hang
PDF
HTML
• #### 摘要

为提高多松弛(MRT)格子Boltzmann模型的计算效率, 运用反演法提出了一个求解三维不可压缩流的12速MRT格子Boltzmann模型(iD3Q12 MRT模型). 这个模型比通常使用的D3Q13 MRT模型具有更高的计算效率. 在数值模拟部分我们把iD3Q12 MRT模型与可压缩性较小的一个13速多松弛模型(He-Luo D3Q13 MRT模型)在精确性和稳定性方面作比较. 通过模拟不同的流动, 包括压力驱动的稳态泊肃叶流、周期变化的压力驱动的非稳态脉动流、顶盖驱动的方腔流, 可以发现iD3Q12 MRT模型模拟以上三种流动时得到的数值解与解析解或与已有的结果符合很好, 这说明我们提出的iD3Q12 MRT模型是准确的. 在模拟稳态的泊肃叶流时, 两个模型计算的速度场的全局相对误差完全相同, 且两个模型都具有二阶的空间精度. 在模拟非稳态脉动流时, 大多情况下是12速模型的计算误差更小, 但在脉动流的最大压降增大时, iD3Q12 MRT模型先发散, 这说明He-Luo D3Q13 MRT模型具有更好的稳定性. 在模拟不同雷诺数下的顶盖驱动的方腔流时, He-Luo D3Q13 MRT模型也比iD3Q12 MRT模型更稳定.

#### Abstract

In order to improve the computational efficiency of multiple-relaxation-time lattice Boltzmann model (MRT), a 12-velocity multiple-relaxation-time lattice Boltzmann model (iD3Q12 MRT model) for three-dimensional incompressible flows is proposed in this work by using an inversion method. This model has higher computational efficiency than the commonly used D3Q13 MRT model in principle. In numerical simulations, the accuracy and stability of iD3Q12 MRT model are validated by simulating different flows, including steady Poiseuille flow driven by pressure, unsteady pulsatile flow driven by periodic pressure and lid-driven cavity flow. We also compare the iD3Q12 MRT model with the 13-velocity multiple-relaxation-time lattice Boltzmann model(He-Luo D3Q13 MRT model).For the Poiseuille flow and pulsatile flow, the numerical solutions of the iD3Q12 MRT model agree well with the analytical solutions. In terms of accuracy, the iD3Q12 MRT model and He-Luo D3Q13 MRT model are used to simulate Poiseuille flow with different parameters. The global relative errors of the two models are identical. Similarly, we also simulate the pulsatile flow to calculate the global relative errors of flow fields at different times and different lattice spacing. It is found that the global relative errors of the iD3Q12 MRT model are smaller than those of the He-Luo D3Q13 MRT model, and both models have the second-order spatial accuracy. Furthermore, we also simulate the pulsatile flow by changing the lattice spacing or relaxation time when the maximal pressure drop of the channel is increased, and it is found that the global relative errors calculated by the iD3Q12 MRT model are smaller than those by the He-Luo D3Q13 MRT model in most cases, but the iD3Q12 MRT model diverges when the maximal pressure drop of the channel is large. This indicates that the iD3Q12 MRT model is more accurate than the He-Luo D3Q13 MRT model in simulating unsteady pulsatile flow, but less stable. For the lid-driven cavity flow, the results show that the numerical results of the iD3Q12 MRT model agree well with those given by Ku et al [Ku H C, Hirsh R S, Taylor T D 1987 J. Comput. Phys. 70 439]. In terms of stability, the iD3Q12 MRT model is quantitatively less stable than He-Luo D3Q13 MRT model.

#### 作者及机构信息

###### 通信作者: 张文欢, zhangwenhuan@nbu.edu.cn
• 基金项目: 浙江省自然科学基金(批准号: LQ16A020001)、浙江省教育厅科研基金(批准号: Y201533808)、宁波市自然科学基金(批准号: 2016A610075)和宁波大学王宽诚幸福基金资助的课题

#### Authors and contacts

###### Corresponding author: Zhang Wen-Huan, zhangwenhuan@nbu.edu.cn
• Funds: Project supported by the Natural Science Foundation of Zhejiang Province, China (Grant No. LQ16A020001), the Scientific Research Foundation of the Education Department of Zhejiang Province, China (Grant No. Y201533808), the Natural Science Foundation of Ningbo, China (Grant No. 2016A610075), and the K.C. Wong Magna Fund in Ningbo University, China.

#### 施引文献

• 图 1  三维泊肃叶流示意图

Fig. 1.  The schematic of three-dimensional Poiseuille flow.

图 2  泊肃叶流数值解与解析解的对比　(a) 泊肃叶流在$x=1$截面处z取不同的值时水平速度$u_{x}$y变化的函数图像; (b) 在截面$z=0$y取不同的值时压力px变化的函数图像; 直线: 解析解; 符号: 数值解; 松弛因子$\lambda_{\nu}=1.3$

Fig. 2.  Comparison between numerical and analytical solutions of Poiseuille flow: (a) The variation of $u_{x}$ with y for different locations of z at section $x=1$ for Poiseuille flow; (b) the variation of pressure with x for different locations of y at section $z=0$ for Poiseuille flow. Lines, analytical solutions; symbols, numerical results; the relaxation parameter ${\lambda}_{\nu}=1.3$.

图 3  不同的$\lambda_{\nu}$下, 模拟泊肃叶流得到的速度场的全局相对误差${\rm {GRE}}_u$随空间步长$\text{δ}{x}$的变化, 符号代表数值解, 连线表示拟合直线

Fig. 3.  The variation of ${\rm {GRE}}_u$ of velocity field with the lattice spacing $\text{δ}{x}$ at different $\lambda_{\nu}$ for Poiseuille flow. Symbols represent numerical solutions, lines represent fitting line.

图 4  $\eta=2.8285$时脉动流在$x=1$, $z=0$处水平速度uxy变化的函数. 直线: 解析解; 符号: 数值解

Fig. 4.  The variation of horizontal velocity ux with y for pulsatile flow at the location $x=1$, $z=0$, $\eta=2.8285$. Line, analytical solutions; symbols, numerical solutions.

图 5  同一周期四个不同时刻下变量${\rm {GRE}}_u$随空间步长$\text{δ}{x}$的变化

Fig. 5.  The variation of ${\rm {GRE}}_u$ with the lattice spacing at four different times in a period for pulsatile flow.

图 6  三维顶盖驱动的方腔流示意图

Fig. 6.  The schematic of three-dimensional lid-driven cavity flow

图 7  不同的雷诺数下模拟方腔流, 在截面$z=0.5$处竖直和水平中心线的速度分布　(a) Re = 100; (b) $Re=400$; (c) $Re=1000$

Fig. 7.  The velocity distribution in the vertical and horizontal center lines at section $z=0.5$ for cavity flows at different $Re$: (a) $Re=100$; (b) $Re=400$; (c) $Re=1000$.

•  [1] 马聪, 刘斌, 梁宏. 耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟. 物理学报, 2022, 71(4): 044701. doi: 10.7498/aps.71.20212061 [2] 陈效鹏, 冯君鹏, 胡海豹, 杜鹏, 王体康. 基于格子Boltzmann方法的二维气泡群熟化过程模拟. 物理学报, 2022, 71(11): 110504. doi: 10.7498/aps.70.20212183 [3] 陈效鹏, 冯君鹏, 胡海豹, 杜鹏, 王体康. 基于格子Boltzmann方法的二维汽泡群熟化过程模拟. 物理学报, 2022, (): . doi: 10.7498/aps.71.20212183 [4] 张恒, 任峰, 胡海豹. 基于格子Boltzmann方法的幂律流体二维顶盖驱动流转捩研究. 物理学报, 2021, 70(18): 184703. doi: 10.7498/aps.70.20210451 [5] 李洋, 苏婷, 梁宏, 徐江荣. 耦合界面力的两相流相场格子Boltzmann模型. 物理学报, 2018, 67(22): 224701. doi: 10.7498/aps.67.20181230 [6] 王佐, 张家忠, 王恒. 非正交多松弛系数轴对称热格子Boltzmann方法. 物理学报, 2017, 66(4): 044701. doi: 10.7498/aps.66.044701 [7] 任晟, 张家忠, 张亚苗, 卫丁. 零质量射流激励下诱发液体相变及其格子Boltzmann方法模拟. 物理学报, 2014, 63(2): 024702. doi: 10.7498/aps.63.024702 [8] 解文军, 滕鹏飞. 声悬浮过程的格子Boltzmann方法研究. 物理学报, 2014, 63(16): 164301. doi: 10.7498/aps.63.164301 [9] 陶实, 王亮, 郭照立. 微尺度振荡Couette流的格子Boltzmann模拟. 物理学报, 2014, 63(21): 214703. doi: 10.7498/aps.63.214703 [10] 史冬岩, 王志凯, 张阿漫. 任意复杂流-固边界的格子Boltzmann处理方法. 物理学报, 2014, 63(7): 074703. doi: 10.7498/aps.63.074703 [11] 刘邱祖, 寇子明, 韩振南, 高贵军. 基于格子Boltzmann方法的液滴沿固壁铺展动态过程模拟. 物理学报, 2013, 62(23): 234701. doi: 10.7498/aps.62.234701 [12] 郭亚丽, 徐鹤函, 沈胜强, 魏兰. 利用格子Boltzmann方法模拟矩形腔内纳米流体Raleigh-Benard对流 . 物理学报, 2013, 62(14): 144704. doi: 10.7498/aps.62.144704 [13] 曾建邦, 李隆键, 蒋方明. 气泡成核过程的格子Boltzmann方法模拟. 物理学报, 2013, 62(17): 176401. doi: 10.7498/aps.62.176401 [14] 苏进, 欧阳洁, 王晓东. 耦合不可压流场输运方程的格子Boltzmann方法研究. 物理学报, 2012, 61(10): 104702. doi: 10.7498/aps.61.104702 [15] 曾建邦, 李隆键, 廖全, 蒋方明. 池沸腾中气泡生长过程的格子Boltzmann方法模拟. 物理学报, 2011, 60(6): 066401. doi: 10.7498/aps.60.066401 [16] 曾建邦, 李隆键, 廖全, 陈清华, 崔文智, 潘良明. 格子Boltzmann方法在相变过程中的应用. 物理学报, 2010, 59(1): 178-185. doi: 10.7498/aps.59.178 [17] 张新明, 周超英, Islam Shams, 刘家琦. 用格子Boltzmann方法数值模拟三维空化现象. 物理学报, 2009, 58(12): 8406-8414. doi: 10.7498/aps.58.8406 [18] 卢玉华, 詹杰民. 三维方腔温盐双扩散的格子Boltzmann方法数值模拟. 物理学报, 2006, 55(9): 4774-4782. doi: 10.7498/aps.55.4774 [19] 李华兵, 黄乒花, 刘慕仁, 孔令江. 用格子Boltzmann方法模拟MKDV方程. 物理学报, 2001, 50(5): 837-840. doi: 10.7498/aps.50.837 [20] 吕晓阳, 李华兵. 用格子Boltzmann方法模拟高雷诺数下的热空腔黏性流. 物理学报, 2001, 50(3): 422-427. doi: 10.7498/aps.50.422
• 文章访问数:  4506
• PDF下载量:  47
• 被引次数: 0
##### 出版历程
• 收稿日期:  2019-06-26
• 修回日期:  2019-09-14
• 上网日期:  2019-11-27
• 刊出日期:  2019-12-05

## 三维不可压缩流的12速多松弛格子Boltzmann模型

• 1. 宁波大学数学与统计学院, 浙江　315211
• 2. 华中科技大学数学与统计学院, 湖北　430074
• ###### 通信作者: 张文欢, zhangwenhuan@nbu.edu.cn
基金项目: 浙江省自然科学基金(批准号: LQ16A020001)、浙江省教育厅科研基金(批准号: Y201533808)、宁波市自然科学基金(批准号: 2016A610075)和宁波大学王宽诚幸福基金资助的课题

/