lammps教程：python ovito编程绘制速度分布

大家好，我是小马老师。

本文介绍使用python和ovito绘制速度分布曲线的方法。

在lammps流体模拟中，都会用到速度分布曲线，以应流体的流动性能。

有两种方法绘制速度分布曲线，一是在lammps in文件中使用分块计算的原理直接输出速度分布数据，二是使用编程的方法对数据进行后处理得到速度分布曲线。

本文介绍第二种方法，使用python和ovito绘制速度分布曲线。

ovito处理速度分布和lammps原理基本一致，也是先进行分块切割，得到每一块的原子数据后，再进行数据求和平均。

使用的命令为：SpatialBinningModifier

这个模块在ovito软件中也有，不过是pro版本才可以用。

如果没有ovito pro版本，则可以使用python调用ovito库函数的方法实现同样的功能。
这种方法完全免费，只需要安装python和ovito库即可。
安装方法：
python ovito模块安装方法
安装完成后，新建flow.py，复制粘贴下面的代码。

该案例模拟了原子沿x方向的流动过程，需要计算的是原子沿y方向的速度分布。

from ovito.io import import_file, export_file
from ovito.modifiers import SpatialBinningModifier
#读入dump文件
pipeline = import_file("dump.xyz")
#分块处理速度数据
pipeline.modifiers.append(SpatialBinningModifier(
    property = 'Velocity.X',
    direction = SpatialBinningModifier.Direction.Y, 
    bin_count = 15,
    reduction_operation = SpatialBinningModifier.Operation.Mean,
))
#输出
export_file(pipeline, 'velocity_profile.txt', 'txt/table', key='binning')

运行python：

python flow.py

运行完成后，会生成一个文本文件：velocity_profile.txt，该文件一共两列数据，第一列是Y方向的每个分块bin的坐标，第二列为原子X方向的速度。
绘制曲线如下图所示：

-------------------------------
公众号：lammps加油站