lammps教程：2D正弦粗糙界面建模建模教程

大家好，我是小马老师。
在前面的推文中，介绍了一维正弦粗糙界面的建模方法。

本文给出二维正弦粗糙界面建模代码。
二维曲面和一维曲面建模原理一样，只要能够构造出二维曲面的方程即可。
一维曲面沿y方向没有变化，而二维曲面y方向同样按照正弦曲线变化，因此，需要在方程中加入y的变化。
二维曲面方程形式为：

z=Asin(ωx)*sin(ωy)+k

对应的代码为：

variable        v equal "v_zz < (15 * sin(v_xx * 2.0*PI/50)*sin(v_yy * 2.0*PI/50) +30 ) "
create_atoms    1 box var v set x xx set y yy set z zz

最终效果为：

公众号：lammps加油站

微信：sunnyfirst888
QQ群：227321877