基于格子气元胞自动机的二维流体动力学模拟系统

资 源 简 介

本项目旨在利用格子气元胞自动机理论在MATLAB环境下实现对二维流体运动的数值模拟。系统通过构建离散的格点空间，并在格点间定义粒子的运动矢量，模拟流体微观粒子的平移与碰撞过程。核心功能包括实现FHP六角网格模型或HPP方格模型，能够模拟流体在受到外力作用或经过障碍物时的流动状态。具体实现过程分为迁移和碰撞两个基本步骤：迁移步模拟粒子沿特定方向移动至相邻格点，碰撞步则根据质量守恒和动量守恒原则重新分配格点内的粒子分布函数。该模拟系统支持自定义复杂的边界条件，如周期性边界、反弹边界（模拟无滑移固体壁面）以及压

详 情 说 明

基于格子气元胞自动机（LGCA）的流体动力学模拟系统

项目介绍

本项目是一个利用格子气元胞自动机（LGCA）理论，在MATLAB环境下实现的二维流体运动数值模拟系统。系统采用经典的FHP-I（Frisch-Hasslacher-Pomeau）模型，通过在六角格点上模拟大量离散微观粒子的迁移与碰撞，从而展现宏观流体力学现象。该系统能够模拟流体在管道内绕过圆形障碍物的复杂流动过程，是研究统计力学与流体力学底层联系的理想教学与实验平台。

功能特性

1. FHP-I 六角格点模型实现 系统通过在二维空间构建六角网格支架，定义了六个方向的粒子运动矢量。相比于普通的方格模型（如HPP），本系统采用的六角模型能够满足流体模拟所需的旋转不变性，真实还原纳维-斯托克斯方程的特性。

2. 核心演化算法 系统严格执行迁移（Streaming）和碰撞（Collision）两个基础步骤。迁移步处理粒子沿既定矢量的位移；碰撞步则根据质量守恒和动量守恒原则，处理二体和三体碰撞，实现动量的重新分布。

3. 复杂边界条件模拟 支持多种边界设置：

• 反弹边界（Bounce-back）： 当粒子遭遇障碍物或上下壁面时，其运动方向将发生180度反转，有效模拟了真实流体中的无滑移固体壁面。
• 入口驱动： 系统在左侧入口处通过控制粒子注入概率，模拟压力驱动或特定强度的流体注入。

4. 宏观物理量提取与平滑 系统提供了从微观布尔状态向宏观物理量转换的功能，包括局部密度、速度矢量场以及系统平均动能。为了消除LGCA模型固有的统计噪声，内置了基于卷积核的空间平均平滑算法。

5. 实时可视化监控 程序运行期间可动态展示密度分布云图、速度场云图、流速矢量场以及系统平均动能随时间的变化曲线，方便观测如压力波传播、障碍物后方流场扰动等物理过程。

使用方法

1. 参数配置 用户可以在程序开头修改模拟参数：

• 修改网格尺寸以改变模拟区域的大小。
• 调整初始密度以观察不同雷诺数下的流体状态。
• 设置入口扰动强度或障碍物的位置与大小。

2. 启动仿真 在MATLAB命令行窗口运行主体程序。系统将自动初始化格点并弹出可视化窗口，展示流体演化的完整过程。

3. 数据获取 仿真结束后，系统会自动将稳态下的密度场和速度场数据保存为MATLAB格式的数据文件，以便后续进行更深入的定量分析。

系统要求

• MATLAB R2016b 或更高版本。
• 无需特殊的附加工具箱，核心计算基于矩阵运算实现。

实现逻辑与细节分析

1. 六角网格的矩形映射逻辑 由于计算机内存采用矩形数组，本系统通过奇偶行平移的技术解决了六角网格的坐标映射。在奇数行和偶数行，粒子向斜向移动时会根据行数的奇偶性选择不同的横向偏移量，从而在高效率的矩阵位移（circshift）基础上保证了六角几何结构的物理准确性。

2. 迁移步（Streaming）的实现 系统利用矩阵位移操作，将六个方向的粒子分布分布层同步沿各自的运动矢量方向移动一个单位格点。这种实现方式最大限度地利用了MATLAB的矢量化运算能力，避免了慢速的循环遍历。

3. 碰撞步（Collision）的布尔逻辑 碰撞逻辑严格遵循FHP-I准则：

• 二体碰撞： 当两个粒子对向相撞且其余四个方向为空时，粒子将以50%的概率偏转60度或-60度。
• 三体碰撞： 当三个互成120度的粒子相撞时，粒子方向将反转。
• 这些逻辑通过纯粹的位运算和掩码计算实现，确保了在模拟大量粒子时仍能保持极高的计算速度。

4. 障碍物碰撞处理 在迁移步结束后，系统通过逻辑索引找出所有进入障碍物区域的粒子。利用反弹规则，将这些粒子的方向分量互换（如方向1与方向4互换），使粒子回到流场，模拟了固体壁面对流体的阻滞作用。

5. 统计与去噪算法 LGCA模型在单格点上是离散的0或1，直接观察会有强烈的噪声。系统采用5x5的均值卷积核对宏观量进行空间滤波，从而获得平滑的流速和密度连续分布，并基于该数据计算系统总动能，用于判断流场是否达到稳定状态。