24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 2041  |  回复: 18
本帖产生 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 的情况。呵呵。
11楼2010-09-30 13:30:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★
ghcacj(金币+3):谢谢 2010-09-30 14:23:09
另外你说的 整体的NVE(或NVT或NPT)+局部NVT(或NPT)是不行的,因为那样的话内部原子被积分了两次。
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进行升温。随着温度的上升,压强总是从一个正压强逐渐下降并且出现负压强。对于负压强总是感觉与实际情况不符,材料内部升温熔化应该压强不断上升才合理(感觉应该这样),解释的时候难以自圆其说,所以就想解决负压强的问题。
13楼2010-10-01 15:36:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluewhale

铁杆木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+2):谢谢 2010-10-01 22:27:06
看起来你的力场可能不太适合于做高压模拟,不知具体是什么材料或使用的势?
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势,应该可以,在文献中看很多人都用过的
15楼2010-10-01 20:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-10-02 11:28:57
负压强不是问题。照道理说,体积不变,升温了,压强确实应该增大,但你研究的不是压力的变化吧,这个问题可以暂时忽略它不管。如果不大放心,你可以做一个整体的加热,看看压力变化情况。这个问题有可能是势函数造成的。
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:
负压强不是问题。照道理说,体积不变,升温了,压强确实应该增大,但你研究的不是压力的变化吧,这个问题可以暂时忽略它不管。如果不大放心,你可以做一个整体的加热,看看压力变化情况。这个问题有可能是势函数造 ...

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

老虎大王

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3):谢谢老虎大王 2010-10-03 12:03:35
这个问题认真起来,可以分析每个原子的应力情况。当然这也是有意义的,就是在内部熔化的时候,压力是如何传递的,这个可能要用不同的势函数来试试。这个工作我以前也想过要做,可是一直没有做。不过你要研究的目标是什么,你先把你的目标完成,这些属于旁枝的,可以以后再说,压力的问题比较复杂,可以先回避它,只要你得到的结构是合理的。
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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见