lammps模拟激光加热案例

lammps模拟激光加热案例

大家好，我是小马老师。本文介绍lammps模拟激光加热的案例。 lammps模拟激光的原理是把激光的能量输入到被照射区域，用到的命令为fix heat。 在热导系数的模拟中，fix heat也被用来持续输入热量，两者原理是一样的，区别是激光能量高，持续时间短，热导模拟中，输入的热量相对较小，持续时间长。 fix heat命令的语法格式为：

fix ID group-ID heat N eflux （region）

group-ID：被加热的原子组N：输入热量的间隔eflux：热流，指单位时间内输入的能量，metal单位下为ev/psregion：激光照射区域，这个关键词是可选的，不写region表示激光照射group组的全部原子，设置region区域后，激光仅加热region区域内的原子。

下面以激光加热Cu球为例，简单介绍fix heat命令的使用，在本例中激光仅加热上部的Cu球。

我平时举例子以Cu、Al、Fe较多，有些粉丝以为我只会做这几种材料的模拟。用这几种材料较多，主要原因是我比较懒，这几种材料势文件较多，不用单独再去-，随便从一个文件夹复制过来就可以用。 我写推文基本以介绍lammps模拟技巧或者lammps命令为主，像本例，重点是fix heat命令的用法，用哪种材料并不是重点，Cu可以用，换另外一种材料同样可以用。 in文件代码：

units metalatom_style atomicboundary p p ptimestep 0.001neighbor 2 binneigh_modify every 1 delay 0 check yes#读取模型read_data cu.data#设置激光加热区域region 2 sphere 36 18 30 10 units boxgroup sphere type 2#力场设置pair_style eam/fspair_coeff * * Cu1.eam.fs Cu Cuthermo 100#温度初始化velocity all create 300 989#激光加热dump 1 all custom 100 cu.xyz id type x y zfix 1 all nve fix 2 sphere heat 10 2.0 region 2run 10000

模拟结果：

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。