伪势模型
-
基于MATLAB的LBM两相流Shan-Chen模型数值模拟
本项目旨在利用MATLAB开发一套高效的二维格子Boltzmann方法（LBM）数值模拟程序，专门用于解决不混溶两相流体的动力学问题。核心算法基于Shan-Chen伪势模型（Pseudo-potential Model），采用经典的D2Q9离散速度模型和单弛豫时间（SRT-BGK）碰撞算子。程序主要功能涵盖：1. 精确构建LBM演化方程，包括平衡态分布函数的计算、碰撞步与流式步的迭代求解；2. 实现多相流相互作用机制，通过计算流体粒子间的非局部相互作用力来模拟表面张力和相界面自动捕获，无需预设界面位置；3. 多样化边界条件处理，支持全周期边界、标准反弹格式（Bounce-back）及半程反弹格式，以模拟封闭容器或无限域流动；4. 复杂的流体现象模拟，能够重现自发相分离（Spinodal Decomposition）、液滴在剪切流中的变形与破裂、液滴聚并以及不同润湿性壁面上的接触角现象；5. 高效的数据后处理与可视化，利用MATLAB强大的绘图引擎实时渲染密度场云图、速度场矢量图，并支持导出流场演化的动态视频，自动统计模拟过程中的平均密度、最大速度与相分离度等关键物理参数。
