multiphase velocities——VOF
 

hello everybody.

I have a question when using VOF,and program is listed below,the C_U(c,pri_ph), C_U(c,sec_ph) and C_U(c,mix_ph) returns same values, so the relative velocity between primary phase and secondary phase is zero.Why?

DEFINE_SOURCE(Flg_x,c,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(c,sec_ph)>0 && C_VOF(c,sec_ph)<1) { ...... NV_D(ul, =, C_U(c,sec_ph),C_V(c,sec_ph),0); NV_D(ug, =, C_U(c,pri_ph),C_V(c,pri_ph),0); NV_VV(sub, =, ul, -,ug); Message("velocity is : %g\n,%g\n",ul[0],ug[0]); ...... } } thanks for your help! David.

Re: multiphase velocities——VOF
 

By design VOF uses a common velocity field.

Re: multiphase velocities——VOF
 

Thank you, Joe. Then how can i get the velocity of primary and secondary phase respectively?