那要用密度定义宏

up

明白你的意思了，仍然需要定义比热，那一回事，温度改成压力

我用loop循环获取单元内的压力和温度。
但是标准初始化会卡死
DEFINE_SPECIFIC_HEAT(UDF_cp, T, Tref, h, yi)
{
real Tz, P, cp;
Domain *domain = Get_Domain(1);
cell_t c;
Thread *cthread;
thread_loop_c(cthread, domain)
{
begin_c_loop(c, cthread)
{
Tz = C_T(c, cthread);
P = C_P(c, cthread);
cp = /* 定义比热容与压力温度的函数 */
*h = cp * (T - Tref);
}
end_c_loop(c, cthread)
}
return cp;
}

up

你这个程序明显不对，比热的udf很难写的，它是最特殊的一个函数了。。我会写！但是还是要看你的函数形式，有些还是没法写的