[vidyadhar]

Hi,


I am trying to store Source term in UDM (user defined memory) in locations where the volume fraction of secondary phase lies between 0 and 1, i.e., at the interface. (volume fraction =1 corresponds to secondary phase and 0 corresponds to primary phase)


/*if (NULL != THREAD_STORAGE(cell_thread,SV_VOF_G))*/
if(NNULLP(THREAD_STORAGE(cell_thread,SV_VOF_G)))
{
mag = NV_MAG(C_VOF_G(cell,cell_thread));
printf("VOF_mag= %f\n",mag);
if ((NV_MAG(C_VOF_G(cell,cell_thread)))>0)
{
source=2.; C_UDMI(cell,cell_thread,0) = source;
printf("source= %g\n",source);
}
else
{
source=0;

printf("source= %g\n",source);
}
}


But, what I observed is that the commands inside the if-else conditioner is not getting executed AND user defined memory in the entire domain is assigned a value of zero.


If I remove the if(NNULLP(THREAD_STORAGE(cell_thread,SV_VOF_G))) statement it results in segmentation fault. From the forums, I partially understood that C_VOF_G data is freed up in memory and should not be allowed to do so by using text user-interface commnad: /solve/set/expert



I request any one to help me in this regard.


Thanks in advance!
vidyadhar