应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【讨论】高温高压的体系需要注意的环节
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
192/2返回列表上一页12
查看: 2041  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 模拟EPI ，点击这里进行查看

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+5, 模拟EPI+1):老虎答疑这么久了，奖励EPI是必须的。 2010-09-30 14:23:02
1。初步判断，fix all nve＋选定部分速度标定，是比较合适的。如果你要研究的是内部原子加热后传递到外部的情况，那么这个内部原子＋外部原子所组成的整体，它与外界应是绝热的，否则就相当于这个整体又处于外部热浴之下。

2. 关于你说的传热很快的问题，这是你的系统的自然特性。我觉得不需要刻意影响它，你模拟工作本身就是设好一个条件，让它自然演化，你如果人为干预了，这个过程就不是自然的过程。

3. 如果你一定要干预这个过程，希望控制外部原子的温升，你也许可以用fix all nvt＋局部温度标定的形式，这个就强行把外部原子温度控住了。然后在你需要整体升温的时候，再切回到 fix all nve。管说明书上说temp/rescale 一般情况下不能与 fix nvt一起用，但你这个就属于abnormal 的情况。呵呵。
赞一下(2人) 回复此楼
11楼2010-09-30 13:30:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
ghcacj(金币+3):谢谢 2010-09-30 14:23:09
另外你说的 整体的NVE(或NVT或NPT）＋局部NVT（或NPT)是不行的，因为那样的话内部原子被积分了两次。
赞一下(1人) 回复此楼
12楼2010-09-30 13:35:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
★ ★
zh1987hs(金币+2):鼓励参与讨论 2010-10-01 22:26:54
引用回帖:
Originally posted by 老虎大王 at 2010-09-30 13:30:50:
1。初步判断，fix all nve＋选定部分速度标定，是比较合适的。如果你要研究的是内部原子加热后传递到外部的情况，那么这个内部原子＋外部原子所组成的整体，它与外界应是绝热的，否则就相当于这个整体又处于外部热 ...

开始一直用fix all nve＋fix temp/rescale进行升温。随着温度的上升，压强总是从一个正压强逐渐下降并且出现负压强。对于负压强总是感觉与实际情况不符，材料内部升温熔化应该压强不断上升才合理（感觉应该这样），解释的时候难以自圆其说，所以就想解决负压强的问题。
赞一下(1人) 回复此楼
13楼2010-10-01 15:36:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluewhale

铁杆木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+2):谢谢 2010-10-01 22:27:06
看起来你的力场可能不太适合于做高压模拟，不知具体是什么材料或使用的势？
赞一下(2人) 回复此楼
14楼2010-10-01 17:27:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
★ ★
zh1987hs(金币+2):谢谢 2010-10-01 22:27:14
引用回帖:
Originally posted by bluewhale at 2010-10-01 17:27:09:
看起来你的力场可能不太适合于做高压模拟，不知具体是什么材料或使用的势？

用EAM势，应该可以，在文献中看很多人都用过的
赞一下(1人) 回复此楼
15楼2010-10-01 20:56:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-10-02 11:28:57
负压强不是问题。照道理说，体积不变，升温了，压强确实应该增大，但你研究的不是压力的变化吧，这个问题可以暂时忽略它不管。如果不大放心，你可以做一个整体的加热，看看压力变化情况。这个问题有可能是势函数造成的。
赞一下(2人) 回复此楼
16楼2010-10-02 11:17:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
★ ★
zh1987hs(金币+2):谢谢 2010-10-02 20:03:28
引用回帖:
Originally posted by 老虎大王 at 2010-10-02 11:17:42:
负压强不是问题。照道理说，体积不变，升温了，压强确实应该增大，但你研究的不是压力的变化吧，这个问题可以暂时忽略它不管。如果不大放心，你可以做一个整体的加热，看看压力变化情况。这个问题有可能是势函数造 ...

看来是需要把压力暂时忽略了，不知道是不是我的模型建立出现了比较多的界面，并且界面处是无规则的（可能出现了有点地方界面间隙较大）。在温度上升的时候，界面出现的可能是个崩塌过程，而不是出现较高的压强直接传递到另一相。其实如果出现压强一直随着上升的话，那么观察压强的变化对我来说也将是很有意义的，因为也希望知道随着温度的变化压强是如何变化的。谢谢了！
赞一下(1人) 回复此楼
17楼2010-10-02 19:57:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3):谢谢老虎大王 2010-10-03 12:03:35
这个问题认真起来，可以分析每个原子的应力情况。当然这也是有意义的，就是在内部熔化的时候，压力是如何传递的，这个可能要用不同的势函数来试试。这个工作我以前也想过要做，可是一直没有做。不过你要研究的目标是什么，你先把你的目标完成，这些属于旁枝的，可以以后再说，压力的问题比较复杂，可以先回避它，只要你得到的结构是合理的。
赞一下(2人) 回复此楼
18楼2010-10-03 11:32:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
引用回帖:
Originally posted by 老虎大王 at 2010-10-03 11:32:15:
这个问题认真起来，可以分析每个原子的应力情况。当然这也是有意义的，就是在内部熔化的时候，压力是如何传递的，这个可能要用不同的势函数来试试。这个工作我以前也想过要做，可是一直没有做。不过你要研究的目标 ...

谢谢您的详细解答。先做着看看吧。
赞一下 回复此楼
19楼2010-10-03 21:56:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭