[求助]如何根据速度场求解流线

LEVEN

这里的速度场是用三个三位数组Bx By Bz 表示的。其中Bx[x,y,z]表示在(x,y,z)点的速度的x方向的分量，期望得到三个一维数组lx，ly，lz。其中lx表示流线第i个点的x轴坐标。当然，流线的起点是知道的，即lx[0],ly[0],lz[0]已知。
流线问题理论上可以通过求解 dx / Bx = dy / By = dz / Bz （即 dr X B(x,y,z) = 0 )来得到流线方程。但是由于Bx，By，Bz只是用数组表示，而非解析的函数。所以，没法用现成的数值计算软件。还往那位高手能指点一下，只需说一下算法思路即可。
先谢谢了…

mycar07

可以查看一下VTK可视化软件包中的源码,vtkStreamer.h,vtkStreamer.cxx,vtkStreamLine.h,vtkStreamLine.cxx的源码,使用四阶龙格-库塔积分方法,由起始点出发,以一个小步长积分获得下一个空间坐标点,由该点周围的已知点线性插值,活得新空间点的速度,循环操作,直到达到退出条件~祝好运~

leven

谢谢啦 我先看看。。。

mysyy

可以参照MATLAB中有关插值的算法