24小时热门版块排行榜    

查看: 1262  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
请填处理意见