多孔介质计算（非独立多孔介质）

多孔介质计算（非独立多孔介质）source = -con*fabs(C_U(c t))*C_U(c t); // C_U 是获取 x方向速度con = C2*0.5*C_R(c t)*x[1]; // C_R 是获取网格密度#include "udf.h" // 头文件#define C2 100.0 // 定义恒定常数DEFINE_SOURCE(xmom_source c t dS eqn) {real x[ND_ND]; // 声明一个实数型的数组，对于二维等价于real x[2]，对于三维等价于real x[3]，ND_ND是宏，有fluent自动确定其取值，二维取2，三维取3。这句的用处一般是用来存放坐标或和坐标有关的矢量的，三维情况下，x[0]代表x轴方向，x[1]代表y轴方向，x[2]代表z轴方向。real con source; // 声明两个实数类型的变量C_CENTRO

利用 DEFINE_SOURCE(name c t dS eqn) 宏指定源项

以下名为 xmom_source 的 UDF用于在 ANSYS fluent 中添加源项，源代码可以被解释或编译，该函数生成一个随y位置变化的x方向的动量源项。

图 1 计算模型

下面给是 UDF 源代码 ，用于在空间坐标系中指定 x 方向动量源项

#include "udf.h" // 头文件
#define C2 100.0 // 定义恒定常数
DEFINE_SOURCE(xmom_source c t dS eqn)
{
real x[ND_ND]; // 声明一个实数型的数组，对于二维等价于real x[2]，对于三维等价于real x[3]，ND_ND是宏，有fluent自动确定其取值，二维取2，三维取3。这句的用处一般是用来存放坐标或和坐标有关的矢量的，三维情况下，x[0]代表x轴方向，x[1]代表y轴方向，x[2]代表z轴方向。

real con source; // 声明两个实数类型的变量

C_CENTROID(x c t); // 获取网格中心坐标，动量源项表达式中需要用到 y 坐标值

con = C2*0.5*C_R(c t)*x[1]; // C_R 是获取网格密度

source = -con*fabs(C_U(c t))*C_U(c t); // C_U 是获取 x方向速度

dS[eqn] = -2.*con*fabs(C_U(c t)); // 是对源项求导，也就是对源项线性化

return source;
}