高手看看这个传质的UDF有啥问题

jesseliuxl

做的液液传质，用VOF+species transport模拟液液不相容体系界面之间的传质，在两相作用里的mass transfer里加的UDF，迭代计算之后在display—counters—species里没有发现有组分传到另一相中，大家看看问题在哪？


#include "udf.h"
#include "sg_mphase.h"
#define kl 0.00002

DEFINE_MASS_TRANSFER

(oil_water_source,cell,thread,from_index,from_species_index,to_index,to_species_index)
{
  real m_l;
   
  Thread *oil = THREAD_SUB_THREAD(thread,from_index);
  Thread *water = THREAD_SUB_THREAD(thread,to_index);
  if (C_VOF(cell,oil)==1)
{
  m_l=0;

}
else if ((C_VOF(cell,oil)>0)&&(C_VOF(cell,oil)<1))

{
  m_l=kl*(C_YI(cell,oil,0)-C_YI(cell,water,0));

  }
return m_l;
}