[Howard818]

Hi Friends,

I want to use volume fraction gradient in define_source. but I get an error "Error: received a fatal signal (Segmentation fault)".

Could anyone help me resolve this problem?

Thank you.




////////////////////////////////////////////////////////////////
#include "udf.h"


DEFINE_EXECUTE_AT_END(calculate_value)

{ Thread *t;
Thread **pt;
cell_t c;
Domain *domain = Get_Domain(1);
real n1, n2, n3, n4;

mp_thread_loop_c (t,domain,pt)

{

begin_c_loop(c, t)

{

n1 = C_VOF_G (c, pt[0])[0];

n2 = C_VOF_G (c, pt[0])[1];

n3 = C_VOF_G (c, pt[1])[0];

n4 = C_VOF_G (c, pt[1])[1];

C_UDMI(c, t, 0) = n1;

C_UDMI(c, t, 1) = n2;

C_UDMI(c, t, 2) = n3;

C_UDMI(c, t, 3) = n4;

}end_c_loop(c, t)

} }


DEFINE_SOURCE(source_wx, cell, thread_l, dS, eqn)
{

real U_Dx_L, U_Dx_G;
real aa, bb, source;
Thread *mix_thread, *thread_g;

mix_thread=THREAD_SUPER_THREAD(thread_l);
thread_g=THREAD_SUB_THREAD(mix_thread, 0);

aa=C_UDMI(cell, thread_l, 1);

bb=C_UDMI(cell, thread_l, 0);


U_Dx_L=C_U(cell,thread_l)*aa/C_VOF(cell, thread_l);

U_Dx_G=C_U(cell,thread_g)*bb/C_VOF(cell,thread_g);

source=U_Dx_L+U_Dx_G;

dS[eqn]=0.0;
return source;
}