|
发表于 2011-1-16 16:23:31
|
显示全部楼层
#include"udf.h" ;;
DEFINE_ADJUST(my_adjust,d) ;;名字为my_adjust 的宏,用于计算域内耗散之和
;;没用过这个udf,好像里面积分的结果是耗散率与体积的乘积,没明白目的是什么
{
Thread*t; ;;声明一个指针,用于标记
/*Integrate dissipation.*/
real sum_diss=0.; ;;声明一个存储变量,此变量最后输出
cell_t c; ;;声明一个指针,指向cell
thread_loop_c(t,d) ;;开始遍历d上的所有cell
{
begin_c_loop(c,t) ;;
sum_diss+=C_D(c,t)*C_VOLUME(c,t); ;;相当于积分作用,C_D是获得c上的耗散率,C_VOLUME是c的体积
end_c_loop(c,t) ;;结束对域的遍历
}
printf("Volume integral of turbulent dissipation:%g\n",sum_diss);;;打印遍历结果
} |
|