基于龙格库塔四阶积分的流线可视化方法.docVIP

基于龙格库塔四阶积分的流线可视化方法 　　摘要：龙格库塔积分的基本思想是通过多次点定位和速度插值的方式抑制误差，提高精度，是工程中常用到的单步但有不失精度的算法，其中，四阶龙格库塔性能最好文章基于CUDA平台对四阶龙格库塔积分做了实现，借助GPU强大的并行计算能力，提高积分计算的速度，使得使用四阶龙格库塔积分计算生成的流线可视化达到实时性的要求。 中国论文网 /8/view-7218478.htm 　　关键词：龙格库塔积分；点定位；速度插值 　　1概述 　　基于几何形状的三维流场可视化技术相比于基于纹理合成的可视化技术是一种较为传统的流场可视化方式。采用几何形状来进行三维流场的可视化，首先要确定使用何种几何形，主要有字形的如线段和箭头，线形如流线、迹线，面形如流面、流管、流带等。选定几何形状后，接下来需要选取合适的积分器来完成流场可视化中积分工作。 　　积分器中最具代表性的两个就是欧拉积分器和龙格库塔积分器。作为数值积分呈现三维流场可视化技术中的两大经典方法，欧拉积分和龙格库塔积分有着各自的特点和优势。欧拉积分的优势在公式简单，易推导，从流线的积分生成公式中可以直接得来，在对步长进行选取后，可以快速地从初始点积分得到下一个点的位置信息，计算量不大，在积分前对由模拟器得到的原始数据进行网格化的划分后，得到规整的流场数据，然后用这些数据作为输入来进行数值积分操作，虽然在此同时会丢失掉流体的质量信息，不再满足流体质量守恒定律，使得最终结果和真实值有所偏差，但误差仍在可接受范围内，不会出现现实不可能现象的发生。相比于欧拉积分，龙格库塔积分显得更为复杂一些，但带来的是误差的减小和精度的提升，但同时计算能力也要求得到提升以满足龙格库塔积分所要求的大量计算任务。 　　2数据预处理 　　要进行数据预处理首先是获得流场原始数据，由于原始流场数据复杂多变，数据量及其庞大，所以先要对得到的流场原始数据进行预处理以得到较为清晰简洁的输入数据。根据需要进行预处理的实际原始数据类型采用不同的预处理应对方案，主要应用的技术包括：挑选原始数据中真正能够在程序的计算过程中需要使用的，减少荣冗余信息量，提取关键信息，可以通过对原始数据进行过滤，流场特征的检测提取等方法来进行；针对分布稀疏的流场数据，由于可能会导致最后绘制的流线都一样，体现不出流场的特征来，需要对这种影响视觉效果的数据进行插值处理，常用的预处理技术还有很多如噪声的消除，计算速度矢量的法矢量，还有变换参数域等。在进行完上述原始数据的过滤精炼后，还要进行的一项就是数据组织。 　　本文采用体元（即体素）划分的方式。体元划分的目的是让流场中各个场值点的分布相互影响，以达到促进流场分布均匀化和各点影响的局部化，以此来调整流场拓扑结构，避免绘制结果中出现狭长体元。体元划分其实就是在选定的区域内对体元顶点进行线性插值，使区域内点对流场的影响局限在相邻的体元间。体元划分采用的体元主要有四面体和六面体，由于精度要求的原因，本文采用六面体来进行实际的体元划分。基于六面体体元的数据组织方法是目前流场可视化中应用得很多的数据组织方法之一，通过将矢量场数据依次分割为相同大小的六面体体元，广泛应用在规则数据网格和虽不规则但已进行过结构化的网格上。根据有限元理论可以完成一个六面体体元对应于一个雅可比矩阵的求解。 　　3龙格库塔四阶积分 　　龙格库塔积分家族中最常用的就是四阶龙格库塔积分法，简称为“RK4”。四阶龙格库塔积分利用给定的矢量场初值信息，省去求解微分方程的复杂过程，通过积分求得矢量场的特征。相比欧拉积分法，减小了误差，提高了求解的精度。其主要的公式内容如下所示： 　　其中xn是当前积分位置，xn+1是积分所得的下一点位置，dt是选定的步长，v（xn）表示矢量场中位置xn的速度矢量大小，0（dt2）表示计算值与真实值之间的误差大小。图1展示了当给定步长dt为1.0时，四阶龙格库塔积分法的计算过程，与欧拉积分的图例相似，其中绿线表示真实值（groundtruth），蓝线是由龙格库塔积分计算得到，红线为计算值与真实值之间的误差，虚线和红线黑箭头的线表示对每次求得的部分结果分截所得，如图1所示。 　　4可视化映射及结果的显示与绘制 　　在此基础上，最后需要将经过预处理和体元划分后得到的规则的流场数据通过数值积分的方式生成具有形状、亮度、颜色等属性的几何元素（流线，流管，流面），如图2所示。通过利用目前已经很成熟的计算机图形学的方法和理论，使用多种绘制算法实现显示结果的多样化，包括几何变换，颜色材料属性，光照设置，纹理贴图，遮挡，隐藏面消除，深度纹理和阴影，反走样等绘制技术将可视化映射计算得到的几何数据和属性解释成图像处理数据输出到要最终的显示设备上，完成流场的可视化工作。