lammps教程:python ovito编程绘制速度分布 |
您所在的位置:网站首页 › 编程语言python教程 › lammps教程:python ovito编程绘制速度分布 |
大家好,我是小马老师。 本文介绍使用python和ovito绘制速度分布曲线的方法。 在lammps流体模拟中,都会用到速度分布曲线,以应流体的流动性能。 有两种方法绘制速度分布曲线,一是在lammps in文件中使用分块计算的原理直接输出速度分布数据,二是使用编程的方法对数据进行后处理得到速度分布曲线。 本文介绍第二种方法,使用python和ovito绘制速度分布曲线。 ovito处理速度分布和lammps原理基本一致,也是先进行分块切割,得到每一块的原子数据后,再进行数据求和平均。 使用的命令为:SpatialBinningModifier 这个模块在ovito软件中也有,不过是pro版本才可以用。 如果没有ovito pro版本,则可以使用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加油站 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |