| 查看: 2091 | 回复: 18 | ||||
| 本帖产生 1 个 模拟EPI ,点击这里进行查看 | ||||
xiaowu787木虫 (正式写手)
|
[交流]
【讨论】高温高压的体系需要注意的环节已有4人参与
|
|||
|
希望LAMMPS模拟的模型中心产生一个高温高压(有时间寿命)的局部区,然后观测温度和压力的向周围扩展的情况。可是温度高压强高时会出现lost atoms,现在 第一个疑问是:周期边界条件下的image atom怎么不起作用了--周期性的边界条件不是可以保证粒子数不变? 第二: 局部升温的选择: 1. fix all nve+fix 局部 temp/rescale,可是压强无法控制 2. fix all nve+fix 局部 nvt or npt? 符合实际情况吗? 3. fix all npt+fix 局部 nvt or npt? 符合实际情况吗? 4. fix all nvt+fix 局部 nvt or npt? 符合实际情况吗? 那种选择更好一些? |
» 收录本帖的淘帖专辑推荐
动力学或热力学计算 |
» 猜你喜欢
三无产品还有机会吗
已经有6人回复
垃圾破二本职称评审标准
已经有7人回复
投稿返修后收到这样的回复,还有希望吗
已经有7人回复
博士申请都是内定的吗?
已经有14人回复
谈谈两天一夜的“延安行”
已经有13人回复
氨基封端PDMS和HDI反应快速固化
已经有11人回复
之前让一硕士生水了7个发明专利,现在这7个获批发明专利的维护费可从哪儿支出哈?
已经有11人回复
论文投稿求助
已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
高温高压反应釜防火措施
已经有5人回复
高温高压反应釜防火措施
已经有10人回复
求《安全血液和血液制品》原本或相关题库
已经有4人回复
【求助】一个高温高压催化试验台的若干问题
已经有10人回复
【转帖】大力推进煤矿瓦斯抽采利用
已经有7人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-09-29 18:02:49
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
3楼2010-09-29 20:37:00
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4):谢谢 2010-09-29 22:31:27
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4):谢谢 2010-09-29 22:31:27
|
记得先前有位版友也提到过 高温的时候 lost atoms。 我不是用lammps的,所以具体不清楚什么原因。但是据经验分析,lost atoms,是不是因为粒子速度过大,在步长之内的位移过大,以至于越出了边界? 温度的扩散容易模拟和控制,但是压强的扩散就麻烦一点了。一般来说,压强可以通过外加应力,或者改变晶格参数来实现。如果不是恒压,应该选取NPT的系综,而且体积大小、形状都允许发生变化。 建议,先做温度扩散,再做压力扩散,然后再综合考虑。 |
4楼2010-09-29 21:08:59
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
5楼2010-09-29 21:15:44
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
6楼2010-09-29 21:33:56
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
7楼2010-09-29 21:35:40
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
8楼2010-09-29 22:28:35
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
9楼2010-09-29 23:07:27
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
★ ★ ★ ★
zh1987hs(金币+4):谢谢,鼓励交流 2010-10-01 15:44:50
zh1987hs(金币+4):谢谢,鼓励交流 2010-10-01 15:44:50
|
时间步长设置 0.01fs 300K初始弛豫,时间步长改为 1fs 后开始升温2000K(整体NVE+部分temp/rescale(速度标度))也可以进行。但是有个问题--由于整体NVE,导致在部分进行速度标度是,温度向周围传递的很快(形成连续的温度梯度),希望的温度梯度最好不是很连续--中间快速升温时,两端温度上升的越慢越好,当随着时间的推移,温度真正传递到这里时才出现真正的升温过程。 [ Last edited by xiaowu787 on 2010-9-30 at 09:54 ] |
10楼2010-09-30 09:51:16













回复此楼
