一、初始化对UDF运行的影响

初始化只会初始化网格上的物理数据、在UDF中常用的实际时间flow-time（CURRENT_TIME）、迭代步数N_ITER、UDM中的数据（其实也就是网格物理数据）。

初始化之后，UDF程序中的静态变量不会再初始化（程序变量初始化和fluent的初始化不同，请注意区分），也即fluent初始化后再计算不会让程序再回到第一次执行的状态，UDF程序依稀记得初始化前静态变量在运算过程中得到取值。

不知道各位读者理没理解上面这句话，如果不理解请再继续阅读下面的。

我举个例子。假设我们需要在UDF中实现写入数据，具体来说需要每次迭代结束时判断flow-time是否经历了一次写入周期（比如1秒写入一次），如果经历了一个周期就写入一次，而且第一次一定要写入一次。

为了实现上述特征，我编写了如下UDF代码。

DEFINE_EXECUTE_AT_END(iter_ending)
    {
        real End_iter_time=0.0;
        real curing_ratio_ave=0.0;  /* Volume weighted average curing degree */
        real heat_ave=0.0;  /* Volume weighted average heating energy */


        #if !RP_HOST
        Thread *sheet_ct;
        int sheet_id = sheet_zone_id;
        cell_t c;

        real curing_ratio_weighted_sum=0.0;
        real heat_weighted_sum=0.0;
        real volume_sum=0.0;
        #endif


        #if !RP_NODE
        FILE *fp = NULL;
        char filename[]="out.txt";
        static int remain_num = -1;  /* Static variables are initialized only once */
        #endif



        #if !RP_HOST
        sheet_ct = Lookup_Thread(Get_Domain(1), sheet_id);

        begin_c_loop_int(c, sheet_ct)
        {
            curing_ratio_weighted_sum += C_VOLUME(c, sheet_ct) * C_UDMI(c, sheet_ct, 0);
            heat_weighted_sum += C_VOLUME(c, sheet_ct) * C_UDMI(c, sheet_ct, 3);
            volume_sum += C_VOLUME(c, sheet_ct);
        }
        end_c_loop_int(c, sheet_ct)

        curing_ratio_weighted_sum = PRF_GRSUM1(curing_ratio_weighted_sum);
        heat_weighted_sum = PRF_GRSUM1(heat_weighted_sum);
        volume_sum = PRF_GRSUM1(volume_sum);
        curing_ratio_ave = curing_ratio_weighted_sum / volume_sum;
        heat_ave = heat_weighted_sum / volume_sum;
        End_iter_time = CURRENT_TIME;

        #endif

        node_to_host_real_4(End_iter_time, bou_temp, curing_ratio_ave, heat_ave);

        #if !RP_NODE
        if(N_ITER<1)  /* Clear the file when open the file for the first time */
        {
            if ((fp = fopen(filename, "w"))==NULL)
                Message("\nWarning: Unable to open %s for writing\n", filename);
            else
                Message("\nWriting info to %s...\n", filename);
                fprintf(fp, "real-time \t bou_temp \t  curing_ratio_ave \t heat_ave\n");
        }
        else
        {
            if ((fp = fopen(filename, "a"))==NULL)
                Message("\nWarning: Unable to open %s for writing\n", filename);
            else
                Message("\nIt is available to Write data to %s...\n", filename);
        }
        #endif

        #if !RP_NODE
        Message("real-time:%lf \t bou_temp:%lf \t  curing_ratio_ave:%lf \t heat_ave:%lf\n", End_iter_time, bou_temp, curing_ratio_ave, heat_ave);
        if(remain_num != (int) (End_iter_time / writig_frequency))
        {
            remain_num = (int) (End_iter_time / writig_frequency);
            Message("\nWritting data to %s...\n", filename);
            fprintf(fp, "%lf \t %lf \t %lf \t %lf \n", End_iter_time, bou_temp, curing_ratio_ave, heat_ave);
        }
        fclose(fp); 
        Message("Iter_Done\n");
        #endif
    }