24小时热门版块排行榜    

查看: 1625  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 4/200 2026-09-01 06:38 by o3EMxHUAQDtS
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 4/200 2026-09-01 06:07 by o3EMxHUAQDtS
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 06:06 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:45 by o3EMxHUAQDtS
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见