[eng.masoud]

Hi fellows,
I am simulating a 2phase problem (liquid and air) that liquid phase comes from inlet and depositing on wall substrate. I am using heat source to apply and make it molten.
I defined the code as below, to select only second phase cells and apply on them but after starting Ive got the segmentation error. Would you please guide me on this regard.

Thanks


#include "udf.h"
DEFINE_SOURCE(heatsource, c, thsec, dS, eqn)
{
real x[ND_ND], source, current_time;
real xc; /* xc - x coordinate */
real yc; /* yc - y coordinate */

Thread *thmix;
thsec = THREAD_SUB_THREAD(thmix,3);

C_CENTROID(x,c,thsec);
xc=x[0]; /* the x coordinate of the cell - 0th component of vector x */
yc=x[1]; /* the y coordinate of the cell - 1st component of vector x */

if (xc>=0.0045 & xc<=0.0077 & yc>=-0.001 & yc<=0.0029)

{
source=10e11;}

else {
source=0;}

dS[eqn] = 0.0; /* the source does not depend on temperature so dS = 0 */

return source;

}