精选优质文档-倾情为你奉上浸入边界法1 IB-LBM的基本数学构造Hofler和Schwarzer在2000年提出了一种在对NS方程求解的有限差分方程中增加一项体积力的IBM来模拟充满粒子流场的固体粒子的运动。其粘性不可压流的控制方程表示为: (1) (2)f:作用在流体上的外力上述两个公式在IB-LBM中可有格子Boltzmann方程所代替,即:(3)(4)(5)F表示作用在物体表面上的力,其物理过程包括:碰撞过程:(6)迁移过程:(7)由此可知,碰撞仅仅是局部的,但迁移和邻近点皆有关系。平衡态分布函数为:(8)其中为声速,式中表示各离散速度方向的加权系数,其值分别为: (9)该模型中,流体的松弛时间和运动粘性系数的关系可表示为:(10)1.1 IB-LBM速度修正法在LBM中,同时满足NS方程的宏观量是通过Chapman-Enskog展开得到的。在速度修正法里其表达式方程有如下表示:(11) (12)其中,为中间层速度,为修正速度。结合方程(11),(1