找回密码
 注册
查看: 1813|回复: 2

求高手帮忙 已知方程求指点粒子轨迹udf

[复制链接]
发表于 2011-5-27 20:32:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
粒子被放在流体域中预先定好的位置,他们的位置通过积分下面公式:
dx(t)/dt=v(x,y,z)
采用四阶Runge--Kutta方法积分得

K1l=vl(t0,x0,y0,z0)
K2l=vl(t0+h/2,x1,y1,z1)
K3l=vl(t0+h/2,x2,y2,z2)
K4l=vl(t+h,x3,y3,z3)
(l=1,2,3)
x(i+1)=x(i)+h*(k[1][1]+2*k[2][1]+2*k[3][1]+k[4][1])/6
y(i+1)=y(i)+h*(k[1][2]+2*k[2][2]+2*k[3][2]+k[4][2])/6
z(i+1)=z(i)+h*(k[1][3]+2*k[2][3]+2*k[3][3]+k[4][3])/6
其中
x[j+1]=x[j]+k[j][1]*h/2
y[j+1]=y[j]+k[j][2]*h/2
z[j+1]=z[j]+k[j][3]*h/2
(j=1,2,3,4)
求指点该用哪些函数,编程思路改怎么样呢,小妹不胜感激啊
 楼主| 发表于 2011-5-29 21:25:18 | 显示全部楼层

回复 1# ltzzq 的帖子

咋就没人指点一下啊
发表于 2017-12-10 22:31:58 | 显示全部楼层
ltzzq 发表于 2011-5-29 21:25
咋就没人指点一下啊

楼主,我碰到你一样的问题,这个有办法解决吗
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表