24小时热门版块排行榜    

查看: 586  |  回复: 3

dgwfly

新虫 (正式写手)

[求助] 求教下MD中Beeman算法 ?? 已有1人参与

想用分子动力学中的Beeman算法求解粒子的运动方程。查了文献,算法的形式如下,需要用t和t-\Delta t时刻的加速度来预测t+\Delta t时刻的位置。想问下在写程序时,可以根据初始位置计算初始时刻的力,但我如何表示在 初始时刻 前一\Delta t时刻的加速度呢,难道是0吗??

求教下MD中Beeman算法 ??


发自小木虫Android客户端
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果我没有记错的话, 这种积分方法无法自己启动, 在实际应用时, 开始几步是使用其他高精度的积分方法的, 如RK6之类, 然后再换用这个方法向下积分
2楼2016-07-26 09:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dgwfly

新虫 (正式写手)

引用回帖:
2楼: Originally posted by jerkwin at 2016-07-26 09:50:30
如果我没有记错的话, 这种积分方法无法自己启动, 在实际应用时, 开始几步是使用其他高精度的积分方法的, 如RK6之类, 然后再换用这个方法向下积分

谢谢,好主意。这样子的话,写程序的时候,我可以用velocity verlet等算法计算出前两步的力,然后再用beeman的算法算。

发自小木虫Android客户端
3楼2016-07-26 23:48:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福
4楼2016-07-27 00:58:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dgwfly 的主题更新
信息提示
请填处理意见