| 查看: 1173 | 回复: 16 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
fannieyang金虫 (小有名气)
|
[交流]
分子动学模拟中的一个问题
|
||
|
做了一个孤立系统分子动力学模拟测试（即没有热浴、没有耗散力，只有带电粒子之间的库伦作用），发现系统的动能、势能随时间都在变大，直到模拟崩掉。感到有这么几个困惑： （1）程序应该没有问题，因为加上热浴后的计算结果就是正确的。 （2）假定程序没有问题，就从物理基本原理方面理解：系统是孤立的，那么根据熵增原理，系统状态应该趋向无序演化，直到平衡态。因此我能否这样理解：也许我的程序在崩掉的时候还没有达到平衡态。。。 可是系统的势能、动能随时间都增大显然不符合能量守恒，这怎么解释呢。。。 (3) 无奈，俺只能怀疑求解运动方程中数值误差造成的能量不守恒，于是减小了时间步长测试，果然可以计算的时间延长了，但是动能、势能依然较为缓慢的增加。是不是一定要使用精度更高的数值解法，才能满足能量守恒呢？ 各位大侠怎么看呢？求指教。 |
回复此楼» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
ktwyw
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 11 (小学生)
- 金币: 2726.8
- 散金: 500
- 红花: 41
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化
16楼2015-01-27 10:28:41
阿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 (小学生)
- 金币: 2726.8
- 散金: 500
- 红花: 41
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化
11楼2015-01-26 18:58:43