[lig]

Hi,All,

I am simulating particle-laden air flow through a porous zone. The source term in the transport equation need the velocity and concentration of the air when it flow into a cell. In another word, only two sides of a cell are needed in the calculation in a 2D case---- x, y velocity/concentration when flow direction is opposite to the normal direction of a face.

The current code I wrote as followed (part of the code) has 4 source output values (2D). Could you help me to correct it? Thank a lot.

c_face_loop(c,t,i)
{

ft=C_FACE_THREAD(c,t,i);
f=C_FACE(c,t,i);
c0=F_C0(f,ft);
tc0=THREAD_T0(ft);
c1=F_C1(f,ft);
tc1=THREAD_T1(ft);
scalar=C_UDSI(c0,tc0,7);
F_AREA(Af,f,ft);


NV_D(vel,=,C_U(c0,tc0),C_V(c0,tc0),C_W(c0,tc0));
flowrate=NV_DOT(vel,Af);

source = scalar*CEcell*flowrate/C_VOLUME(c,t);
dS[eqn] = CEcell*flowrate/C_VOLUME(c,t);
Message ("source_7=%g,%g\n",source);

}
return source;