换热系数UDF求助

Roaming

想给一段变截面的管子（厚壁+薄壁两段）外表面设置不同的随温度变化的换热系数，UDF应该怎么来写，只知道应该是用DEFINE_PROFILE宏，求助大佬给些指点

Roaming

我的换热系数的公式是h=3.15*10^6 W_f^0.616 [700+(T_s-700)/(exp((T_s-700)/10)+1)]^(-2.455) [1-1/(exp((T_s-250)/40)+1)]

#include "udf.h"

#define Wf 1.2  // 假设的质量流量，使用时替换为实际值

DEFINE_PROFILE(htc_profile, thread, position)
{
    face_t f;
    real Ts, h;

    begin_f_loop(f, thread)
    {
        Ts = F_T(f, thread);  // 获取壁面温度（单位：K）

        // 计算第二项：[700 + (T_s-700)/(exp((T_s-700)/10)+1)]
        real term2 = 700.0 + (Ts - 700.0)/(exp((Ts - 700.0)/10.0) + 1.0);
        
        // 计算第三项：[1 - 1/(exp((T_s-250)/40)+1)]
        real term3 = 1.0 - 1.0/(exp((Ts - 250.0)/40.0) + 1.0);

        // 完整公式计算
        h = 3.15e6 * pow(Wf, 0.616) * pow(term2, -2.455) * term3;

        F_PROFILE(f, thread, position) = h;
    }
    end_f_loop(f, thread)
}

还请大佬帮忙看看这个UDF有什么问题