Verschaeve Joris C G
Department of Energy and Process Engineering, Norwegian University of Science and Technology, N-7491 Trondheim, Norway.
Phys Rev E Stat Nonlin Soft Matter Phys. 2009 Sep;80(3 Pt 2):036703. doi: 10.1103/PhysRevE.80.036703. Epub 2009 Sep 11.
An analytical and numerical analysis of the no-slip boundary condition at walls at rest for the lattice Boltzmann Bhatnagar-Gross-Krook method is performed. The main result of this analysis is an alternative formulation for the no-slip boundary condition at walls at rest. Numerical experiments assess the accuracy and stability of this formulation for Poiseuille and Womersley flows, flow over a backward facing step, and unsteady flow around a square cylinder. This no-slip boundary condition is compared analytically and numerically to the boundary conditions of Inamuro [Phys. Fluids 7, 2928 (1995)] and Zou and He [Phys. Fluids 9, 1591 (1997)] and it is found that all three make use of the same mechanism for the off-diagonal element of the stress tensor. Mass conservation, however, is only assured by the present one. In addition, our analysis points out which mechanism lies behind the instabilities also observed by Lätt [Phys. Rev. E 77, 056703 (2008)] for this kind of boundary conditions. We present a way to remove these instabilities, allowing one to reach relaxation frequencies considerably closer to 2.
针对格子玻尔兹曼 Bhatnagar-Gross-Krook 方法,对静止壁面处的无滑移边界条件进行了分析和数值分析。该分析的主要结果是给出了静止壁面处无滑移边界条件的一种替代形式。数值实验评估了该形式在泊肃叶流和沃默斯利流、后向台阶绕流以及方柱体周围非定常流中的准确性和稳定性。将这种无滑移边界条件与稻室 [《物理流体》7, 2928 (1995)] 以及邹和何 [《物理流体》9, 1591 (1997)] 的边界条件进行了分析和数值比较,发现这三种边界条件在应力张量非对角元素的处理机制上是相同的。然而,只有本文提出的边界条件能确保质量守恒。此外,我们的分析指出了拉特 [《物理评论 E》77, 056703 (2008)] 针对此类边界条件所观察到的不稳定性背后的机制。我们提出了一种消除这些不稳定性的方法,使得能够达到更接近 2 的松弛频率。