[orchiders]

Hi. everyone.

I’m studying boiling phenomenon and carrying out the vof simulation on bubble growth when boiling.

A boundary between liquid and vapor is not clear.
So I assumed that
0< volume fraction <0.5 : water
0.5< volume fraction <1 : vapor

DEFINE_INIT(my_init_function, domain)
{
Thread *t;
Thread **pt;
Thread **st;
cell_t c;
Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,P_PHASE);
Domain *sDomain = DOMAIN_SUB_DOMAIN(domain,S_PHASE);

real xc[ND_ND], vof1, vof2;

mp_thread_loop_c (t,domain,pt)
if (FLUID_THREAD_P(t))
{
Thread *tp = pt[P_PHASE];

begin_c_loop (c,t)
{
C_CENTROID(xc,c,t);

vof1=C_VOF(c,tp);

if ( 0.5 < vof1 )
C_VOF(c,tp) = 1;

}
end_c_loop (c,t)
}

mp_thread_loop_c (t,domain,st)
if (FLUID_THREAD_P(t))
{
Thread *sp = st[P_PHASE];

begin_c_loop (c,t)
{
C_CENTROID(xc,c,t);

vof2=C_VOF(c,sp);

if ( 0.5 < vof2 )
C_VOF(c,sp) = 0;

}
end_c_loop (c,t)
}

}



and I wrote UDF. but It didn’t operate.

Would you please do me a favor? Please give me reply anytime you feel free. Thank you!