[e_cfd]

Dear all,

I have tried to define source terms on a wall (which is a BC with zero thickness and located between two fluids) using F_C0 and F_C1 but the results are strange. I mean I want to have a fixed value of 500 on the top and -500 on the other side of the wall, but what I get is a range of [0 500] on the top and [-500 0] on the downside at two adjacent mesh cells (I need to define source terms on the very first mesh cells at the wall)!
Any help and suggestions are appreciated (or any alternative ways to define source terms).

My UDF code:

#include "udf.h"
DEFINE_ADJUST(a_source,domain)
{
Domain *d;

face_t f;
cell_t c;
Thread *t;
cell_t c0;
Thread *t0;
cell_t c1;
Thread *t1;

real source=500;
real sink=-500;

d = Get_Domain(1); /* it should be defined to point to the mixture level */

t = Lookup_Thread(d,3); /* the wall ID is 3 */


begin_f_loop(f,t)
{

c0 = F_C0(f,t); /* wall-up side */
t0 = THREAD_T0(t);

c1 = F_C1(f,t); /* wall-down side :
t1 = THREAD_T1(t);

C_UDMI(c0,t0,0)= source;
C_UDMI(c1,t1,1)= sink;

}


end_f_loop(f,t)

}