UDF问题解析.doc

1. real xND_ND 相当于声明一个数组，对于二维等价于real x2，对于三维等价于real x3，ND_ND是宏，有fluent自动确定其取值，二维取2，三维取3。这句的用处一般是用来存放坐标或和坐标有关的矢量的，三维情况下，x0代表x轴方向，x1代表y轴方向，x2代表z轴方向。2.fabs原型：在TC中原型是extern float fabs(float x);，在VC6.0中原型是double fabs( double x );。 用法：#include 功能：求浮点数x的绝对值 说明：计算|x|, 当x不为负时返回x，否则返回-x 类似函数：abs3. C_U(c,t): C_U(cell,thread),笛卡尔（直角）坐标系x轴方向速度。 单元线索（cell thread）在其中分配了材料数据和源项的单元组Thread 数据类型是 FLUENT 中的数据结构。它充当了一个与它描述的单元或面的组合相关的数据容器。4. 请问C_R(c,t)和C_VOF(c,t)的返回值是什么啊？fluent 的udf manual 5.2.3存取材料性质的宏中提到C_R(c,t)参数类型为cell_t c,Thread *t,返回值，密度；在5.6多相宏中提到C_VOF(c,pt0) 参数类型 cell_t c,Thread *pt 返回值 主要相的体积分数，C_VOF(c,ptn) 参数类型 cell_t c,Thread *pt 返回值 第n个辅助相的体积分数。5.6.定义源项不能定义离散坐标辐射模型定义 连续性方程 动量方程能量方程k, 方程种类质量分数方程UDS宏：DEFINE_SOURCE ( name, c, t, dS, eqn)参数类型：cell_t c c is cell thread variableThread *t c is a pointer to threadreal dS 定义dS为实数int eqn 定义方程为整型函数返回值为实数。对于DEFINE_SOURCE宏，有五个参数：name, c, t, dS, and eqn.Name 是UDF的名字，会在FLUENT图形用户界面中变得可见， 且可被选择。C是由t指向线程的单元，在这个单元中源项得到应用。数组dS指定了源项的衍生物（微分），是由传输方程中的变量得到的。这些微分如果能增强求解器的稳定性则可以用来线性化源项6. 关键在于找到那些不能包括到控制方程的非稳态项、对流项与扩散项中所有项，写出源项的表达式。局部线形化的时候注意Sp要小于等于0.7. Fluent UDF定义源项吴子牛编著的计算流体力学基本原理中提到以下几种情况会有源项。1.在非惯性坐标系中，N-S方程右端需添加反映惯性力影响的源项；2.考虑重力影响时，需在动量方程中添加重力项；3.考虑两相流时，在N-S方程右端需条件反映各种流体相互干扰的源项（包括质量交换，动量交换，能量交换）；4.考虑两方程湍流模型时，湍流模型方程右端含有反映湍流生成与消灭的源项；5.考虑化学反应时，需在组分浓度方程和能量方程添加由化学反应引起的质量变化与能量变化；6.将高维问题化为低维问题时，也会出现源项。流体力学中有点源，点汇的概念，我想可能这个与上述提到的应该有些相同吧。 陶文铨老师的书中对源项的描述是这样的：数值传热学书中所指的源项是广义量，它代表那些不能包括到控制方程的非稳态项、对流项与扩散项中所有其它各项之和。书中提到将源项局部线形化，S=Sc+SpTp,这种方法比将源项处理为常数会更合理。7. 温度梯度：具有连续温度场的物体内，过任意一点P温度变化率最大的方向位于等温线的法线方向上，称过点P的最大温度变化率为温度梯度，用gradt表示。#include udf.h#define C2 100.0DEFINE_SOURCE(xmom_source,c,t,dS,eqn) real xND_ND; real con,source; C_CENTROID(x,c,t); con=C2*0.5*C_R(c,