multiphase velocities——VOF

David

hello everybody.

I'm using VOF to modeling my problems,program is list below,the question is: C_U(c,pri_ph),C_U(c,sec_ph)and C_U(c,mix_ph) returns same values, the relative velocity of the primary phase and secondary phase is zero.why?

DEFINE_SOURCE(Flg_x,cell,mix_ph,dS,eqn) { Thread **ph; Thread *pri_ph,*sec_ph; real NV_VEC(ul),NV_VEC(ug),NV_VEC(sub);

ph=THREAD_SUB_THREADS(mix_ph); pri_ph=ph[0]; sec_ph=ph[1]; if(C_VOF(cell,sec_ph)>0 && C_VOF(cell,sec_ph)<1) { ...... NV_D(ul, =, C_U(cell,sec_ph),C_V(cell,sec_ph),0); NV_D(ug, =, C_U(cell,pri_ph),C_V(cell,pri_ph),0); Message("velocity is : %g\n,%g\n,%g\n",ul[0],ug[0]); ...... } }

thanks for your help! David