24小时热门版块排行榜    

查看: 1628  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +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
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见