udf初学求助

学习交流

问题为，热水与冷水混合，要求混合后出口温度不能高于600k，如果温度高于600k，热水入口速度为0，ID为4的面是出口。请高手看看这个对吗？谢谢！
include "udf.h"
#define JK_ID 4 /*面ID*/
float avg_temp=300.0;

DEFINE_ADJUST(JK_avg_temp,domain)
{
        float NV_VEC(A);
        float sum_T_A=0.0,sum_A=0.0;
        face_t f;
        Thread*f_thread;
        f_thread = Lookup_Thread(domain,JK_ID);
        begin_f_loop(f,f_thread)
        {
                F_AREA(A,f,f_thread);
                sum_A+=NV_MAG(A);               
                sum_T_A+=NV_MAG(A)*F_T(f,f_thread);
        }
        end_f_loop(f, f_thread)

        avg_temp=sum_T_A/sum_A;
}

DEFINE_PROFILE(VELOCITY_L,t,position) /*速度入口边界条件*/
{
        face_t f;
        begin_f_loop(f,t)
        {
                        if (avg_temp<=600.00)/*如果该面平均温度低于600k，则速度入口为10m/s*/
                        {
                                F_PROFILE(f,t,position) = 10;
                        }
                        else
                                F_PROFILE(f,t,position) = 0;


                 }

        end_f_loop(f,t)
}

学习交流

请大家帮助