[taosifme014]

Hi Everyone!

I am trying to model a moving heat source on the cells having the free surface (0.05 < VOF < 0.95) using an UDF. I think mainly my problem lies on finding the co-ordinates of the cells having the free surface of VOF value between 0.05 to 0.95. My simulation domain is 720 um * 270 um *210 um. And it has two fluids: air and liquid metal. When I try to run a simulation it shows an error like this in the console. :

" 999999: mpt_accept : error : accept failed : No error "

And the UDF is like this:

#include "udf.h"

#include "math.h"



#define vof_cutoff 0.05

#define P 250

#define eta 0.3
#define S 30e-6

#define d 30e-6

#define pi 3.1416

#define v_laser 300e-3

#define h 20

#define Te 300

#define sig 5.67e-8

#define emi 0.95



DEFINE_SOURCE (heat_source,c,t,dS,eqn)

{

/*Thread *tp = THREAD_SUB_THREAD (t,1);*/

real x[ND_ND];

real X, Y, Z;

real ti = CURRENT_TIME;

real Xs, Zs, a, H, T, I,q;



if ((C_VOF(c,t) > vof_cutoff) && (C_VOF(c,t) <(1- vof_cutoff)) )

{

C_CENTROID(x,c,t);

X = x[0];

Y = x[1];

Z = x[2];

T = C_T(c,t);

I = 1/0.75*(-2.25*pow((Y/S),2.0) +1.5 *(Y/S)+0.75);

Xs = v_laser*ti +0;

Zs = 105e-6;

a = 2*P/(pi*pow(d,2));

H = a * exp(-2*(pow((X-Xs),2.0)+pow((Z-Zs),2)/pow(d,2.0)));

q = (eta*H*I/S - h*(T-Te) - sig*emi*(pow(T,4.0)-pow(Te,4.0)));



}

else

{

q =0;

}

return q;

}



I would highly appreciate if somebody can help me in finding out my problem and give me suggestion or possible solution.

Thanks,

Taosif