求助：fluent 子程序编写-热源和温度耦合问题！
一个传热问题，两种材料组成的一个物体，其中一种材料是一个体热源，体热源随；另一种材料的某个边界的平均温度变化，函数是source=1000000*(1-0.15*(1-0.0045*(avg_temp-298.0)))
热源受边界平均温度影响，同时热源也影响边界平均温度，是一个耦合的关系。按照一些例子，编了以下程序，但是热源和温度没有耦合在一起，不知哪里出了问题。
#include "udf.h"
real NV_VEC(A);
real avg_temp;
real sum_T_A=0.0;
real sum_A=0.0;
Thread *thread_out;
face_t f;
Domain *domain;
DEFINE_ADJUST(adjust,d)
{
domain=Get_Domain(2);
thread_out=Lookup_Thread(domain,6);
begin_f_loop(f,thread_out)
{
F_AREA(A,f,thread_out);
sum_A+=NV_MAG(A);
sum_T_A+=NV_MAG(A)*F_T(f, thread_out);
}
end_f_loop(f,thread_out)
avg_temp=sum_T_A/sum_A;
}
DEFINE_SOURCE(cell_source,c,t,ds,eqn)
{
real source;
begin_f_loop(f,t)
{
source=1000000*(1-0.15*(1-0.0045*(avg_temp-298.0)));
ds=0.0;
return source;}
end_f_loop(f,t)
}
页:
[1]