| 查看: 2094 | 回复: 4 | ||
潇潇雨歇99新虫 (小有名气)
|
[求助]
请问怎么计算Ni在Ni-Al合金中的扩散系数
|
|
我想计算Ni在Ni-Al合金中的扩散系数，还有Al在Ni-Al合金中的扩散系数. 可是得到的msd 斜率全是一样的。 units metal boundary p p p atom_style atomic read_data ni1al9.data group Ni type 1 group Al type 2 pair_style eam/alloy pair_coeff * * NiAl.eam.alloy Ni Al timestep 0.001 thermo 50 thermo_modify lost warn neighbor 0.6 bin neigh_modify every 10 delay 0 check yes variable N equal step variable pote equal pe variable Etotal equal etotal variable T equal temp variable Press equal press variable V equal vol variable Pxx equal pxx variable Pyy equal pyy variable Pzz equal pzz variable Pxy equal pxy variable Pxz equal pxz variable Pyz equal pyz compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 4.0 velocity all create 2500 825577 dist gaussian fix extra all print 50 "${N} ${T} ${pote} ${Etotal} ${Press} ${V}" file 0414ni1al9.dat screen no dump 1 all cfg 20000 a.*.cfg id type xs ys zs c_3 c_4 c_5 dump 2 all atom 20000 dump.2500 fix 1 all npt temp 2500 2500 0.2 iso 0 0 0.2 drag 1 run 100000 unfix 1 write_restart restart.$N undump 1 undump 2 label loop variable x index 2500 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 500 400 variable y index 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 500 400 300 # Cooling steps: fix 1 all npt temp $x $y 0.1 iso 0 0 0.2 drag 1 run 100000 unfix 1 write_restart restart.$y.$N # Equilibrium steps: fix 2 all npt temp $y $y 0.2 iso 0 0 0.2 drag 1 run 100000 unfix 2 write_restart restart.$y.$N # Collection steps: dump 1 all cfg 5000 a$y.*.cfg id type xs ys zs c_3 c_4 c_5 dump 2 all custom 5000 dump.$y.lammpstrj id type x y z xs ys zs compute mymsd all msd fix fix_msd all ave/time 5 10 50 c_mymsd[1] c_mymsd[4] file 0414ni1al9msd.$y.dat compute msd1 Ni msd fix fix_msd1 Ni ave/time 5 10 50 c_msd1[1] c_msd1[4] file 0414nimsd.$y.dat compute msd2 Al msd fix fix_msd2 Al ave/time 5 10 50 c_msd2[1] c_msd2[4] file 0414almsd.$y.dat fix 3 all npt temp $y $y 0.2 iso 0 0 0.2 drag 2 run 100000 unfix 3 unfix fix_msd unfix fix_msd1 unfix fix_msd2 uncompute msd1 uncompute msd2 uncompute mymsd undump 1 undump 2 write_restart restart.$N.$y next x next y jump ni1al9.txt loop 以上是我算MSD的in文件  我听师兄说用group分开1和2就可以 ，但这样算出的结果， msd1 msd2 根合金总的mymsd 三个值一模一样。 所以求问到底怎么才能只算一种原子的MSD @oxox6085@smutao |
回复此楼» 猜你喜欢
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
潇潇雨歇99
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 446.6
- 散金: 130
- 帖子: 127
- 在线: 28.4小时
- 虫号: 3708895
- 注册: 2015-03-04
- 性别: MM
- 专业: 特种医学(航空、航天、航海
2楼2017-04-18 09:39:52
mengfan
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7931.9
- 红花: 3
- 帖子: 639
- 在线: 92.8小时
- 虫号: 390430
- 注册: 2007-06-02
- 性别: GG
- 专业: 金属结构材料
3楼2017-04-18 10:27:06
杨and王
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1710.1
- 帖子: 36
- 在线: 28.5小时
- 虫号: 4659848
- 注册: 2016-05-05
- 性别: GG
- 专业: 金属材料的微观结构
4楼2017-05-26 09:04:52
杨and王
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1710.1
- 帖子: 36
- 在线: 28.5小时
- 虫号: 4659848
- 注册: 2016-05-05
- 性别: GG
- 专业: 金属材料的微观结构
5楼2017-06-30 20:02:27
20
