麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

基于数值积分的流线生成方法及实现

简介

之前介绍过点图标在三维流场可视化中的一些缺陷,其中之一就是无法揭示数据内在的连续性。所以这里的流线生成算法正好解决这一缺陷,可以提前连续的流线。

流线的构造主要有两种:一种是基于数值积分的流线生成算法,另一种是基于流函数的流线生成算法。基于数值计算的流线构造具有计算简单、符合流体计算中的离散思想。所以这里以之为例。

数值积分的原理

我们假设某一流场中质点的运动规律以速度矢量来描述时可以表示为以下形式:

\[v=v(r,t)\]

其中r为点P的位置向量,为了求出该点P的位置关系随时间t的函数关系,则可以通过求解以下微分方程来得到:

\[\frac{dr(t)}{dt}=v(r(t))\]

求解上述方程,可得:

\[r(t+\Delta t)=r(t)+\int v(r(t))dt\]

对于上述方程所述的曲线,我们不可能用一个显式函数来表示,为了得到该曲线,我们引入数值积分。采用数值积分的思想,我们可以离散成许多离散点来表示曲线。只要我们给定了曲线的起始点,采用数值积分的方法就可以得到\(\Delta t,2\Delta t,3\Delta t …t+\Delta t\)时刻P点的位置,将这些离散点采用直线连接,就可以得到近似的曲线。如图:

2013-05-26_010233

数值积分方式:

一阶Euler方法

\[r(t+\Delta t)=r(t)+v(r(t))\cdot\Delta t\]

这种积分方法,逼近结果比较粗糙,精度也较低。

二阶Runge-Kutta方法

\[r*r(t+\Delta t)=r(t)+v(r(t))\cdot\Delta t\]

\[r(t+\Delta t)=r(t)+1/2[v(r(t))+v(r*(t+\Delta t))]\cdot\Delta t\]

该方法是利用t时刻P点的矢量值和P点在t+\(\Delta \)t时刻的矢量值之和的平均值,作为t时刻的矢量大小,再计算下一点。

四阶Runge-Kutta方法
\[a=v(r(t))\cdot\Delta t\]

\[b=v(r(t)+a/2)\cdot\Delta t\]

\[c=v(r(t)+b/2)\cdot\Delta t\]

\[d=v(r(t)+c)\cdot\Delta t\]

\[r(t+\Delta t)=r(t)+(a+2b+2c+d)/6\]

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。
微博关注
手机阅读
赞赏支持
手机阅读
微信捐助麻辣GIS
微信打赏
支付宝捐助麻辣GIS
支付宝打赏

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)

如您有疑问,可在文末留言,或到麻辣GIS QQ群中提问。

QQ群(一):291616564 加入QQ群

QQ群(二):166408035 加入QQ群

QQ群(三):627853279 加入QQ群

QQ群(四):436386604 加入QQ群

QQ群(五):606176554 加入QQ群

QQ群(六):946178380 加入QQ群

QQ群(七):861271808 加入QQ群

新疆互助群:910717627 加入QQ群

微信公众号:malagis,扫描右边二维码直接关注。

发表评论