[chunchen1987]

Dear all, I want to create a between-cell airflow rate matrix using UDF. For example, if there are 4 cells in the domain, I want to create a 4x4 matrix P, where P[i][j] is the airflow rate from cell i to cell j. I wrote a UDF as below, but it doesn't work. The idea is to use F_FLUX(f,t) for each face. Is this the right direction to do that? Or is there any stupid mistakes in the UDF?

I am a new learner of Fluent UDF. Any suggestions would be highly appreciated!!



#include"udf.h"

int P[4][4];

DEFINE_ON_DEMAND(calculate_P)
{
Domain *domain;
Thread *f_thread;
face_t f = -1;
Thread *t = NULL;
cell_t c0, c1 = -1;
Thread *t0, *t1 = NULL;
domain=Get_Domain(1);


thread_loop_f(f_thread, domain)/* loops over all face threads in a domain*/
{
c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
c1 = F_C1(f,t);
t1 = F_C1_THREAD(f,t);


if(F_FLUX(f,t)>0)
P[c0][c1]=F_FLUX(f,t);
else
P[c1][c0]=-F_FLUX(f,t);

}
}