| 查看: 2762 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
lijun_255001金虫 (小有名气)
|
[求助]
在LAMMPS中怎么控制温度变化速率
|
|
| 在LAMMPS中怎么控制温度变化速率，我想按照指数下降的公式（公式已知）来控制温度，请问怎么实现，fix nvt之类的都是线性下降吧。有谁有这方面的经验。 |
回复此楼» 猜你喜欢
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 薄膜在lammps条件下的谐振
已经有5人回复
- 【分享】尝试lammps中, 分享中...
已经有15人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】lammps中怎样检验模拟系统在NVT系综下弛豫是否达到平衡？
已经有6人回复
- 【讨论】lammps中利用velocity命令无法控制住速度问题
已经有11人回复
- 【讨论】lammps中region和group的问题
已经有4人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【讨论】lammps进行弛豫平衡时调解不了温度
已经有3人回复
- 【讨论】lammps运行eam/fs出错
已经有38人回复
5楼2018-01-07 21:31:37
【答案】应助回帖
★ ★ ★
lijun_255001(金币+5): 公式为T=8000exp（-cdt），c可以修改下降速度。用循环语句fix nvt $ $这样您觉得可行吗 2011-11-19 09:42:45
ghcacj(金币+3): 谢谢 2011-11-19 11:31:13
lijun_255001(金币+5): 公式为T=8000exp（-cdt），c可以修改下降速度。用循环语句fix nvt $ $这样您觉得可行吗 2011-11-19 09:42:45
ghcacj(金币+3): 谢谢 2011-11-19 11:31:13
|
分段控制应该可以，你写一个脚本，先从T1以某个速度降到T2，这次得到的restart文件作为初始构型，通过脚本生成新的控制文件，新的控制文件使用脚本自动生成，只需要更改温度控制的那一行，用shell脚本的替换命令，替换掉旧输入文件中温度控制行，再重新执行计算。如此往复，就能任意控制温度。 我不知道lammps有没有其它控制方式，但是这种方式肯定可行。 |
2楼2011-11-18 17:24:24
xunbao15
木虫 (正式写手)
- 模拟EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2243.3
- 散金: 5
- 红花: 1
- 帖子: 457
- 在线: 57.3小时
- 虫号: 292130
- 注册: 2006-10-31
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2011-12-06 02:54:58
rush7
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 633.8
- 散金: 39
- 帖子: 286
- 在线: 75.9小时
- 虫号: 1556717
- 注册: 2011-12-30
- 专业: 金属材料的微观结构
4楼2012-04-25 18:49:53
