vector 赋值_【Fuent高阶教程】利用UDF给边界赋值

VS版本：Visual Studio 2010

Fluent版本：Fluent14.5

本文属于Fluent高阶教程，主要是研究采用最强IDE工具Visual Studio调试UDF代码。

(1)边界赋值介绍

边界赋值是指将已知信息通过插值的方式传输到边界上，包括查询和插值两部分。查询是指在通过位置信息中寻找到未知位置的K个近邻点。插值是指将K个近邻点上的信息插值到未知位置上，获得未知位置点上的信息。在fluent计算中有着较高频率的应用，如将实验结果作为初场插值到对应CFD计算面上，又如将与第三方软件(如EDEM)耦合计算时，将第三方软件计算结果传递到CFD计算边界中来。

(2)近邻查询^【1】

这里主要讲k-d树。k-d树是一种对k维空间中的实例点进行存储以便对其进行快速检索的树形数据结构，k-d树是二叉树。下面举例说明：

给定一个二维空间的数据集: T ={(2,3), (5,4), (9,6), (4,7), (8,1), (7,2)}

1)构造平衡k-d树

·        根结点对应包含数据集T的矩形选择x⁽¹⁾轴，6个数据点的x⁽¹⁾坐标的中位数是7，以超平面x⁽¹⁾ = 7将空间分为左右两个子矩形(子结点)

·        左矩形以x⁽²⁾ = 4为中位数分为两个子矩形

·        右矩形以x⁽²⁾ = 6 分为两个子矩形

·        如此递归，直到两个子区域没有实例存在时停止

   构造的k-d树如下：

2)利用k-d树搜索最近邻

输入：已构造的k-d树；目标点x;

输出：x的最近邻

·        在k-d树中找出包含目标点x的叶结点：从根结点出发，递归的向下访问k-d树，若目标点x的当前维的坐标小于切分点的坐标，则移动到左子结点，否则移动到右子结点，直到子结点为叶结点为止。

·        以此叶结点为“当前最近点”

·        递归地向上回退，在每个结点进行以下操作：(a)如果该结点保存的实例点比当前最近点距离目标点更近，则以该实例点为“当前最近点”;

(b)当前最近点一定存在于某结点一个子结点对应的区域，检查该子结点的父结点的另一子结点对应区域是否有更近的点(即检查另一子结点对应的区域是否与以目标点为球心、以目标点与“当前最近点”间的距离为半径的球体相交)；如果相交，可能在另一个子结点对应的区域内存在距目标点更近的点，移动到另一个子结点，接着递归进行最近邻搜索；如果不相交，向上回退

·        当回退到根结点时，搜索结束，最后的“当前最近点”即为x的最近邻点

(3)插值

空间插值常用于将离散点的测量数据转换为连续的数据曲面，以便与其它空间现象的分布模式进行比较。这里主要将反距离加权插值，其定义如下：

(4)UDF程序实现

这里以某入口边界添加入口流场为例，入口截面的半径为0.1m的圆面，位于xy平面，圆心为(0,0,0)

1)K-d树查询代码

#include

#include

#include

#include

#include

using namespace std;

struct KdTree{

   vector root;

   KdTree* parent;

   KdTree* leftChild;

   KdTree* rightChild;

   //默认构造函数

   KdTree(){parent = leftChild = rightChild = NULL;}

   //判断kd树是否为空

   bool isEmpty()

   {

       return root.empty();

   }

   //判断kd树是否只是一个叶子结点

   bool isLeaf()

   {

       return (!root.empty()) &&

            rightChild == NULL &&leftChild == NULL;

   }

   //判断是否是树的根结点

   bool isRoot()

   {

       return (!isEmpty()) && parent == NULL;

   }

   //判断该子kd树的根结点是否是其父kd树的左结点

   bool isLeft()

   {

       return parent->leftChild->root == root;

   }

   //判断该子kd树的根结点是否是其父kd树的右结点

   bool isRight()

   {

       return parent->rightChild->root == root;

   }

};

//int data[6][3] ={ {2,3,1},{5,4,5},{9,6,4},{4,7,8},{8,1,3},{7,2,7}};

//double data[6][2] ={ {2,3},{5,4},{7.2,6},{4,7},{8,1},{7,2}};

template

vector >Transpose(vector > Matrix)

{

   unsigned row = Matrix.size();

   unsigned col = Matrix[0].size();

   vector > Trans(col,vector(row,0));

   for (unsigned i = 0; i < col; ++i)

   {

       for (unsigned j = 0; j < row; ++j)