27、流体流动模拟与多智能体通信研究

流体流动模拟与多智能体通信研究

流体在曲折通道中流动的格子气模型

在模拟流体在多孔材料孔隙水平的流动时，需要考虑孔隙壁的不同润湿性。
- 碰撞规则
- 无滑移边界条件 ：对于无滑移边界条件，边界层单元格的碰撞规则为 $\theta_{n - s} : (v, (i, j)) \to (v’, (i.j))$，其中 $v’ k = v {k + 3}$。此规则使壁面附近的流体速度为零。若要增强疏水性，可沿边界添加一些虚拟壁面单元格，这些单元格需替换孔隙单元格 $(i, j) \in \Omega_P$，且满足 $|T(i, j) \cap \Omega_P| \geq 4$，选择概率范围从 $p = 0$ 到 $p = 0.5$。
- 滑移碰撞规则 ：粒子进入边界层后，应尽量沿通道壁继续移动且改变方向最小。在规则直通道的多孔介质中，可通过镜面反射规则实现；在具有任意边界曲率的曲折多孔介质中，规则需根据边界单元格与其壁面相邻单元格的相对位置以及单元格邻域内的壁面单元格数量进行转换。若 $T(i, j) \cap \Omega_P = {\varphi_k(i, j), \varphi_{k + 1}(i, j)}$，则 $v’ k = v {k + 4}$，$v’ {k + 1} = v {k + 3}$，$k = 1, \cdots, 6$；否则过渡函数与无滑移情况相同。若要增强亲水性，可添加一些与边界相邻的虚拟壁面单元格来平滑边界，需满足 $|T(i, j) \cap \Omega_P| \leq 2$，选择概率