[求助]谁能帮我看看我的udf程序什么地方错了！！

vion

下面是我用的程序降温的程序：
DEFINE_PROFILE(unsteady_temperature,  thread,position)
{
    face_t f;
    begin_f_loop(f, thread)
{
   real t = RP_Get_Real("flow-time");
   real temperature = C_T(f, thread);
   if (t <= 409.)       （在时间小于409s时，温度从310.15K降到98.15K,监测显示这部分能很好的实现）
    F_PROFILE(f, thread, position) = (-0.5183)*t+310.15;
   else if (409.<t <=1500.) （时间在409s到1500s期间，温度保持在98.15k，这部分也能实现）
    F_PROFILE(f, thread, position) = 98.15;
   else if (1500.<t <=1562.)   （温度在1500s到1562s之间从98.15k迅速升到348.15k，这就不能实现了，温度从1500k一直都是在98.15k不变）
    F_PROFILE(f, thread, position) =4.0322*t-5950.15;
   else if (1562.<t <=3000.)  （温度在1562s至3000s温度保持在348.15k，也就是75度，自然也就不能实现，一直到3000s温度都是在98.15k左右）
    F_PROFILE(f, thread, position) = 348.15;
}
end_f_loop(f, thread)
}

梅花雪剑

看你是考虑边界还是流体区域。

vion

我是给一个壁面程序降温，应该是边界

maxpayne

C_T(f, thread)；
应该是F_T

vion

我以前试过把C改成F，和没有改是一样的。

vion

有没有人知道怎么回事啊，急啊！！！！！

vion

无言了，郁闷中。。。。。。

tiandw

(409.<t <=1500.)
不能这么表达吧，应该用拆成两个用逻辑与连接起来