lammps后处理技巧：速度云图的绘制

大家好，我是小马老师。
本文介绍lammps模拟后处理中速度云图的绘制方法。
在lammps中，云图的绘制应用较多，如应力、密度、温度、速度等等，绘制方法基本一致。

ovito可以根据单原子的参数绘制云图，但是总体效果不是很理想，因为原子之间的量相对离散，绘制出来之后得到的图像过度不平滑。
参考有限元的思想，把整个体系划分为若干个小的单元，一个单元可包含多个原子，计算单元的量（如温度、速度），以此进行绘图，得到的图像就会想的比较平滑。

compute chunk/atom命令可把体系划分为若干个单元（chunk），配合fix ave/chunk命令对每个chunk内的参数进行平均即可实现上述功能。
下面以流体模拟为例，这个案例在之前的推文中已经介绍过。

下面添加速度计算代码：
compute 1 all chunk/atom bin/2d x lower 2 y lower 2 units box
fix 01 all ave/chunk 1 200 200 1 vx file velocity.txt
模拟完成后，使用origin绘制等高线即可得速度云图：

本文介绍的方法，同样可以应用到密度分布、温度分布、应力分布等云图的绘制。