|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
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)
求指点该用哪些函数,编程思路改怎么样呢,小妹不胜感激啊 |
|