| 查看: 1264 | 回復(fù): 16 | |||
fannieyang金蟲 (小有名氣)
|
[交流]
分子動學(xué)模擬中的一個問題
|
||
|
做了一個孤立系統(tǒng)分子動力學(xué)模擬測試(即沒有熱浴、沒有耗散力,只有帶電粒子之間的庫倫作用),發(fā)現(xiàn)系統(tǒng)的動能、勢能隨時間都在變大,直到模擬崩掉。感到有這么幾個困惑: (1)程序應(yīng)該沒有問題,因為加上熱浴后的計算結(jié)果就是正確的。 (2)假定程序沒有問題,就從物理基本原理方面理解:系統(tǒng)是孤立的,那么根據(jù)熵增原理,系統(tǒng)狀態(tài)應(yīng)該趨向無序演化,直到平衡態(tài)。因此我能否這樣理解:也許我的程序在崩掉的時候還沒有達(dá)到平衡態(tài)。。。 可是系統(tǒng)的勢能、動能隨時間都增大顯然不符合能量守恒,這怎么解釋呢。。。 (3) 無奈,俺只能懷疑求解運動方程中數(shù)值誤差造成的能量不守恒,于是減小了時間步長測試,果然可以計算的時間延長了,但是動能、勢能依然較為緩慢的增加。是不是一定要使用精度更高的數(shù)值解法,才能滿足能量守恒呢? 各位大俠怎么看呢?求指教。 |
| 祝福 |
|
|
|
|
至尊木蟲 (文壇精英)

|
|
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6 | 墨墨漠 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|