lammps案例：流体绕过障碍物模拟代码讲解

大家好，我是小马老师。
本文介绍lammps官方案例中的一个流体案例，在模型中设置了两个球形障碍物，模拟流体原子绕过障碍物的过程。

本文对案例代码进行注释，以供初学者参考。
本例lammps in文件代码主要包含以下几个部分。
（1）模型基本参数设置
案例采用二维模型，边界条件为psp边界，y方向为自由边界，在模拟过程中也可以看到y方向尺寸会有变化。
dimension 2
boundary p s p
atom_style atomic
neighbor 0.3 bin
neigh_modify delay 5
（2）建模部分
晶格类型为hex类型，因为是二维模型，box在z方向尺寸只有0.5个晶格长度。
box内包含3种类型的原子。
lattice hex 0.7
region box block 0 40 0 10 -0.25 0.25
create_box 3 box
create_atoms 1 box
（3）力场设置
力场采用LJ势，截断半径为1.12246。