第三类边界条件 UDF 望指点

han_wenjing

各位大侠，我现在在做流固耦合传热计算，几何体是一个竖直圆柱，竖直方向是Z方向，圆柱里面流流体，圆柱边是固体壁面，现在想在固体壁面上加第三类边界条件，即q=h(T1-Tw). 我想选择convection边界条件，其中h已知，就是常数，Tw是固体温度，而T1是流体温度，T1则是随着Z轴（即高度方向）变化的，他们的关系是T1=563.51626+0.0192*z+0.69236*z*z-0.05431*z*z*z+0.00119*z*z*z*z，我自己用UDF尝试的编写了一段程序，想把它嵌进去，因为对UDF不是特别懂，只是仿照帮助文件写了一段，大侠们能不能帮我看看到底是哪里错了吗，计算的结果总是不对，真心希望高手指点！

#include "udf.h"
DEFINE_PROFILE(temperature, thread, position)
{        real x[ND_ND];
        real z;
        face_t f;
        begin_f_loop(f,thread)
        {        F_CENTROID(x,f,thread);
                z=x[2];
                F_PROFILE(f,thread,position)=563.51626+0.0192*z+0.69236*z*z-0.05431*z*z*z+0.00119*z*z*z*z;
        }
        end_f_loop(f,loop)
}