天天看点

lammps不规则界面建模技巧

大家好,我是小马老师。

本文介绍使用Atomsk软件建立不规则界面的方法。

在一些lammps模拟中会涉及到不规则界面对模拟结果的影响,如凹凸界面对摩擦力的影响、凹凸界面对纳米流体的影响等等。

Atomsk可根据设定条件重新调整原子的位置,从而生成不同的结构。

下面以Al板为例介绍不规则界面的建模代码。

(1)生成2层原子的Al板

atomsk --create fcc 4.046 Al -duplicate 20 20 1 Al_supercell.cfg      

结果如图所示:

lammps不规则界面建模技巧

(2)新建disp.txt文件,并输入以下内容:

# Atom displacements
displacement function
uz = 10*cos(4*pi*x/Hx)      

第一句是注释,运行时不起作用,仅说明文件功能。

第二句说明对原子施加displacement(移位)操作,原子位移受函数function控制

第三句设定原子z坐标调整为“10cos(4pi*x/Hx)”

(3)生成不规则界面

atomsk Al_supercell.cfg -properties disp.txt final.lmp      

运行结果如下:

lammps不规则界面建模技巧

可以看出,原子z坐标与x坐标有关,截面为一条余弦线,公式中的10和4分别控制幅度和周期。

disp.txt更改为:

# Atom displacements
displacement function
uz = 2*cos(4*pi*x/Hx) + 2*cos(4*pi*y/Hy)      

运行结果如下图所示:

lammps不规则界面建模技巧

更换不同的公式,可得到不同的界面结构,大家有兴趣的话可以尝试一下。

公众号案例代码以及学习交流,请加QQ群:754749935

扫描关注微信公众号:lammps加油站。