| 查看: 1675 | 回复: 16 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
fannieyang金虫 (小有名气)
|
[交流]
分子动学模拟中的一个问题
|
||
|
做了一个孤立系统分子动力学模拟测试(即没有热浴、没有耗散力,只有带电粒子之间的库伦作用),发现系统的动能、势能随时间都在变大,直到模拟崩掉。感到有这么几个困惑: (1)程序应该没有问题,因为加上热浴后的计算结果就是正确的。 (2)假定程序没有问题,就从物理基本原理方面理解:系统是孤立的,那么根据熵增原理,系统状态应该趋向无序演化,直到平衡态。因此我能否这样理解:也许我的程序在崩掉的时候还没有达到平衡态。。。 可是系统的势能、动能随时间都增大显然不符合能量守恒,这怎么解释呢。。。 (3) 无奈,俺只能怀疑求解运动方程中数值误差造成的能量不守恒,于是减小了时间步长测试,果然可以计算的时间延长了,但是动能、势能依然较为缓慢的增加。是不是一定要使用精度更高的数值解法,才能满足能量守恒呢? 各位大侠怎么看呢?求指教。 |
» 猜你喜欢
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
现代”学阀”该如何界定
已经有14人回复
求合成方法
已经有8人回复
师弟论文见刊大半年才想起来申请专利,还能抢救一下吗?
已经有4人回复
有机合成以后会不会被AI改变?做科研的虫友怎么看
已经有8人回复
我的奶奶
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
fannieyang
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1429.9
- 散金: 72
- 红花: 1
- 帖子: 182
- 在线: 76.8小时
- 虫号: 445802
- 注册: 2007-10-28
- 专业: 基础物理学
12楼2015-01-27 08:20:16
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学

7楼2015-01-25 16:41:36
10楼2015-01-26 14:56:30
ktwyw
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 11 (小学生)
- 金币: 2729.4
- 散金: 500
- 红花: 42
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化

11楼2015-01-26 18:58:43










回复此楼