24小时热门版块排行榜    

查看: 1171  |  回复: 16

fannieyang

金虫 (小有名气)

[交流] 分子动学模拟中的一个问题

做了一个孤立系统分子动力学模拟测试(即没有热浴、没有耗散力,只有带电粒子之间的库伦作用),发现系统的动能、势能随时间都在变大,直到模拟崩掉。感到有这么几个困惑:

(1)程序应该没有问题,因为加上热浴后的计算结果就是正确的。
(2)假定程序没有问题,就从物理基本原理方面理解:系统是孤立的,那么根据熵增原理,系统状态应该趋向无序演化,直到平衡态。因此我能否这样理解:也许我的程序在崩掉的时候还没有达到平衡态。。。       可是系统的势能、动能随时间都增大显然不符合能量守恒,这怎么解释呢。。。
  (3)  无奈,俺只能怀疑求解运动方程中数值误差造成的能量不守恒,于是减小了时间步长测试,果然可以计算的时间延长了,但是动能、势能依然较为缓慢的增加。是不是一定要使用精度更高的数值解法,才能满足能量守恒呢?

各位大侠怎么看呢?求指教。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lx_PICO

铁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
崩掉是指什么情况?
10楼2015-01-26 14:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fannieyang(金币+1): 谢谢参与
祝福
2楼2015-01-25 16:16:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fannieyang(金币+1): 谢谢参与
3楼2015-01-25 16:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fannieyang(金币+1): 谢谢参与
4楼2015-01-25 16:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2015-01-25 16:36:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
6楼2015-01-25 16:37:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿Q~~

至尊木虫 (文坛精英)

路过的看了一下,帮顶~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
7楼2015-01-25 16:41:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
8楼2015-01-25 16:52:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
9楼2015-01-25 17:04:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fannieyang 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见