| 查看: 2552 | 回复: 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? 符合实际情况吗？ 那种选择更好一些？ |
回复此楼» 收录本帖的淘帖专辑推荐
 动力学或热力学计算 |
» 猜你喜欢
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高温高压反应釜防火措施
已经有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 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 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 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 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 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 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 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.4小时
- 虫号: 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
100