[aminhgn]

Hello everyone,

I have 2 phase flow which are a solution and humid-air. I want to define the humidity of water vapor in air at the interface. I know that I should use Define_Adjust(see the following code) but How can I access to only water vapor of humid-air at the interface? The following code pass through the whole domain and cells (even solution side(first species)) but I want to only define water vapor at the interface?

DEFINE_ADJUST(adjust_gradient,domain)

{

int phase_domain_index;

real T1, W;

Thread *t;

cell_t c;

Domain *subdomain;

sub_domain_loop(subdomain, domain, phase_domain_index)

{

if (DOMAIN_ID(subdomain) == 2){

thread_loop_c (t,subdomain)

{

begin_c_loop(c,t)

{

if(C_VOF(c,t)!=0. && C_VOF(c,t)!=1.){

T1 = C_T(c,t);

W = (3.56087706576757E-13*pow(T1,6))+(-6.49782616920242000000E-10*pow(T1,5))+(4.95600960068920000000E-07*pow(T1,4))+(-2.02032299336406000000E-04*pow(T1,3))+(4.63901242081371000000E-02*pow(T1,2))+(-5.68550924813065000000E+00*T1)+2.90434594274724000 000E+02;

C_YI(c,t,0)= W;

C_UDMI(c,t,0) = W;}

}

end_c_loop(c,t)

}

}

}

}