| 查看: 1915 | 回复: 8 | ||
[求助]
万事开头难！请问一下lLAMMPS可以模拟石墨在铁水中的溶解吗？ 已有3人参与
|
| 需要掌握哪方面的知识？求大神指点方向！ |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
完全是火: 金币+30, ★★★很有帮助, 指明了方向！ 2016-01-07 20:03:10
完全是火: 金币+10, ★★★★★最佳答案, 后面的回复非常详细！ 2016-01-08 21:21:50
感谢参与，应助指数 +1
完全是火: 金币+30, ★★★很有帮助, 指明了方向！ 2016-01-07 20:03:10
完全是火: 金币+10, ★★★★★最佳答案, 后面的回复非常详细！ 2016-01-08 21:21:50
|
能不能模拟取决于你能否找到以下势函数： 1，描述高温液态铁中铁原子之间的相互作用的势函数（EAM可以么？MEAM可以么？需要调研才行） 2，描述高温液态铁水中铁原子和碳原子的相互作用的势函数（Fe已经是液态了，能否用LJ？目测不行） 所以，你应该做的第一步就是先调研相关文献，看目前有没有现成的势函数。 如果有，开始学习LAMMPS 如果没有，你有两个选择： 1，放弃这个课题 2，先学习第一性原理，然后用ab initio的方法自己造一个Fe和C的势函数出来，不过很有难度，慎重选择这条路！ |
2楼2016-01-07 17:33:09
|
感谢您的回复！我原先看LAMPPS有很多计算Cu、Al熔化的例子，觉得Fe应该也有。经您分析，发现这些势函数并不通用，需要找到合适的势函数。搜索了一阵，发现还是找不到。想再请教您几个问题： 1. 是不是分别需要找到熔体Fe 原子间的势函数，C 原子间的势函数，Fe和C原子间的势函数，才有可能用LAMMPS计算下去？ 2. 用ab initio构造势函数的话需要哪些基础？ -------------------------------------------------------------- 课题没得换了，不过有看到文献中有人用蒙特卡罗法做出来，不知道您对蒙特卡罗法有涉猎吗？看蒙特卡罗法好像有很多软件，不知道哪个适合我？望您抽空回复！谢谢！ |
3楼2016-01-07 20:02:46
4楼2016-01-08 14:14:38
5楼2016-01-08 15:29:04
wolfshow
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 126.4
- 散金: 200
- 帖子: 28
- 在线: 17.9小时
- 虫号: 575281
- 注册: 2008-06-17
- 性别: GG
- 专业: 固体力学
6楼2016-01-08 15:31:47
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
|
1. 是不是分别需要找到熔体Fe 原子间的势函数，C 原子间的势函数，Fe和C原子间的势函数，才有可能用LAMMPS计算下去？ 没错，需要找到以上三个势函数才行！这是先决条件。 首先是怎么找： C 原子间的势函数：上次看到你这个帖子我想到的是AIREBO势函数，但现在怀疑高温下这个势函数还适不适用？ Fe和C原子间的势函数：Fe3C也许可用，你去找到那篇文献，仔细看他们是在什么条件下适用的。看他们在创造这个势函数的时候测试了哪些参数：是偏重于力学参数还是热学参数。 熔体Fe 原子间的势函数：如果Fe3C可用，那么应该也可以用Fe3C势函数分别描述Fe-Fe和C-C之间的作用。 其次是找到之后，如果这些势函数在LAMMPS有现成的格式可用（例如EAM，MEAM等形式），那么只需要设置参数即可。如果LAMMPS中没有现成的格式，那么需要手动修改LAMMPS源程序并编译链接——工作量将会很大。。。 2. 用ab initio构造势函数的话需要哪些基础？ 你需要学习量子化学为基础，然后可以VASP实际操作。 如果你周围没有人用过VASP，那么你将陷入无数的细节不可自拔。 不知道你LAMMPS现在熟练程度怎么样，但我只能跟你讲VASP比LAMMPS复杂度和容易出错程度高一两个数量级。。。 3，我不会蒙特卡洛，所以无法给你这方面建议，sorry... |
7楼2016-01-08 19:18:21
8楼2016-01-08 21:20:14
9楼2016-01-08 21:26:27
100