| 查看: 1622 | 回复: 16 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
fannieyang金虫 (小有名气)
|
[交流]
分子动学模拟中的一个问题
|
||
|
做了一个孤立系统分子动力学模拟测试（即没有热浴、没有耗散力，只有带电粒子之间的库伦作用），发现系统的动能、势能随时间都在变大，直到模拟崩掉。感到有这么几个困惑： （1）程序应该没有问题，因为加上热浴后的计算结果就是正确的。 （2）假定程序没有问题，就从物理基本原理方面理解：系统是孤立的，那么根据熵增原理，系统状态应该趋向无序演化，直到平衡态。因此我能否这样理解：也许我的程序在崩掉的时候还没有达到平衡态。。。 可是系统的势能、动能随时间都增大显然不符合能量守恒，这怎么解释呢。。。 (3) 无奈，俺只能怀疑求解运动方程中数值误差造成的能量不守恒，于是减小了时间步长测试，果然可以计算的时间延长了，但是动能、势能依然较为缓慢的增加。是不是一定要使用精度更高的数值解法，才能满足能量守恒呢？ 各位大侠怎么看呢？求指教。 |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
ktwyw
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 11 (小学生)
- 金币: 2729.4
- 散金: 500
- 红花: 42
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化
11楼2015-01-26 18:58:43
阿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
fannieyang
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1429.9
- 散金: 72
- 红花: 1
- 帖子: 182
- 在线: 76.8小时
- 虫号: 445802
- 注册: 2007-10-28
- 专业: 基础物理学
12楼2015-01-27 08:20:16
投票:
30