| 查看: 1425 | 回复: 6 | |||
| 当前主题已经存档。 | |||
superdirac木虫 (正式写手)
|
[交流]
【求助】想问个velocity-corrected verlet 算法的问题
|
||
|
v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/12 + O(dt**4) 这个式子是怎么推出来的?? 谢谢! 假如将 v(t+dt/2),v(t-dt/2),a(t+dt),a(t-dt)作泰勒展开,让后合并同类项: 得出的是 v(t)={v(t+dt/2)+v(t-dt/2}/2 - dt*{a(t+dt)-a(t-dt)}/16 + O(dt**4) 两者一直对不上! 找个高手指点一下! 呵呵!谢谢! [ Last edited by lei0736 on 2009-11-25 at 12:07 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2009-03-18 20:29:23
superdirac
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4269.5
- 散金: 33
- 红花: 2
- 帖子: 536
- 在线: 287.5小时
- 虫号: 522969
- 注册: 2008-03-12
- 专业: 理论和计算化学
3楼2009-03-18 22:04:05
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2009-03-19 10:17:10
superdirac
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4269.5
- 散金: 33
- 红花: 2
- 帖子: 536
- 在线: 287.5小时
- 虫号: 522969
- 注册: 2008-03-12
- 专业: 理论和计算化学
★ ★
xuefei06(金币+2,VIP+0):3Q！ 3-19 14:59
xuefei06(金币+2,VIP+0):3Q！ 3-19 14:59
|
呵呵! 这个不是 beeman http://en.wikipedia.org/wiki/Beeman's_algorithm |
5楼2009-03-19 11:15:21
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★
xuefei06(金币+3,VIP+0):thanks! 3-19 15:01
xuefei06(金币+3,VIP+0):thanks! 3-19 15:01
|
陈舜麟《计算材料科学》（化工，2005),p.78, 公式（4-76)，和你这个式子一样。说你这个叫改进的Beeman算法。这书上说，这叫Beeman算法中考虑到时间反演对称性的关于速度的改进式。 我用的程序基本上都是以Verlet算法为基础的，包括蛙跳和速度verlet。Lammps就是用的速度Verlet，DL_POLY可以选用蛙跳或速度verlet。好像Moldy是用的Beeman算法。但你这个式子，我以前真没见过。但陈的书上说这个式子降噪效果好，是个“值得采用的”算式。可惜没有推导。 [ Last edited by 老虎大王 on 2009-3-19 at 12:07 ] |
6楼2009-03-19 12:02:17
superdirac
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4269.5
- 散金: 33
- 红花: 2
- 帖子: 536
- 在线: 287.5小时
- 虫号: 522969
- 注册: 2008-03-12
- 专业: 理论和计算化学
7楼2009-03-19 18:30:33
