| 查看: 1378 | 回复: 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 ] |
回复此楼» 猜你喜欢
2026年国自然面上资助率
已经有20人回复
-
一个有机合成实验室都需要哪些设备？
已经有5人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复
老虎大王
木虫 (著名写手)
- 模拟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