Seta Takeshi, Rojas Roberto, Hayashi Kosuke, Tomiyama Akio
Graduate School of Science and Engineering for Research, University of Toyama, Toyama 930-8555, Japan.
Graduate School of Engineering, Kobe University, Hyogo 657-8501, Japan.
Phys Rev E Stat Nonlin Soft Matter Phys. 2014 Feb;89(2):023307. doi: 10.1103/PhysRevE.89.023307. Epub 2014 Feb 24.
In the present paper, we verify the effectiveness of the two-relaxation-time (TRT) collision operator in reducing boundary slip computed by the immersed boundary-lattice Boltzmann method (IB-LBM). In the linear collision operator of the TRT, we decompose the distribution function into symmetric and antisymmetric components and define the relaxation parameters for each part. The Chapman-Enskog expansion indicates that one relaxation time for the symmetric component is related to the kinematic viscosity. Rigorous analysis of the symmetric shear flows reveals that the relaxation time for the antisymmetric part controls the velocity gradient, the boundary velocity, and the boundary slip velocity computed by the IB-LBM. Simulation of the symmetric shear flows, the symmetric Poiseuille flows, and the cylindrical Couette flows indicates that the profiles of the numerical velocity calculated by the TRT collision operator under the IB-LBM framework exactly agree with those of the multirelaxation time (MRT). The TRT is as effective in removing the boundary slip as the MRT. We demonstrate analytically and numerically that the error of the boundary velocity is caused by the smoothing technique using the δ function used in the interpolation method. In the simulation of the flow past a circular cylinder, the IB-LBM based on the implicit correction method with the TRT succeeds in preventing the flow penetration through the solid surface as well as unphysical velocity distortion. The drag coefficient, the wake length, and the separation points calculated by the present IB-LBM agree well with previous studies at Re = 10, 20, and 40.
在本文中,我们验证了双弛豫时间(TRT)碰撞算子在减少通过浸入边界格子玻尔兹曼方法(IB-LBM)计算的边界滑移方面的有效性。在TRT的线性碰撞算子中,我们将分布函数分解为对称和反对称分量,并为每个部分定义弛豫参数。查普曼-恩斯科格展开表明,对称分量的一个弛豫时间与运动粘度有关。对对称剪切流的严格分析表明,反对称部分的弛豫时间控制着由IB-LBM计算的速度梯度、边界速度和边界滑移速度。对对称剪切流、对称泊肃叶流和圆柱库埃特流的模拟表明,在IB-LBM框架下由TRT碰撞算子计算的数值速度剖面与多弛豫时间(MRT)的速度剖面完全一致。TRT在消除边界滑移方面与MRT同样有效。我们通过解析和数值方法证明,边界速度的误差是由插值方法中使用的δ函数平滑技术引起的。在模拟绕圆柱流动时,基于带有TRT的隐式校正方法的IB-LBM成功地防止了流体穿透固体表面以及非物理速度畸变。由当前的IB-LBM计算的阻力系数、尾流长度和分离点在雷诺数Re = 10、20和40时与先前的研究结果吻合良好。