[mohsensalehi1993]

Hi ٍeveryone,
i'm going to simulate grains of solid rocket engine with udf.
i use DEFINE_PROFILE and DEFINE_CG MOTION macros to link moving boundary to pressure....
my udf's code As follows:


#include "udf.h"
#include "sg.h"

DEFINE_PROFILE(mflux,t,i)
{
int m = 0;
face_t f;
double a,nn,y1,y2;
double PRESSURE,PRESSURE2,press,rdot;
double x[ND_ND];
a=0.001499;
nn=0.4255;
PRESSURE=0;

begin_f_loop(f,t)
{
PRESSURE += F_P(f,t);
m = m + 1;
}
end_f_loop(f,t)
printf("sum P on boundary: %d = %f pas\n", PRESSURE);

PRESSURE2 = PRESSURE/m;
printf("average P on boundary: %d = %f pas\n", PRESSURE2);

press = (PRESSURE2)*1e-6;
rdot = (a * pow(press,nn));


begin_f_loop(f,t)
{

F_CENTROID(x,f,t);
y2=x[1];

if (y2<=0.05)
F_PROFILE(f,t,i) = 1605.*rdot + 1e-6; /* 1605=density*/
else
F_PROFILE(f,t,i) = 1e-6;
}
end_f_loop(f,t)

printf("burning rate: %f\n", rdot);

}


DEFINE_GRID_MOTION(moving_inlet,d,dt,time,dtime)
{
Thread *tf = DT_THREAD(dt);
face_t f;
Node *v;
real NV_VEC(disp);
real y,x,dis,dis_y,dis_x;
double PRESSURE,PRESSURE2,press,rdot;
int n,ig;
int m = 0;
real a=0.001499;
real nn=0.4255;
real k[ND_ND];

PRESSURE=0;

begin_f_loop(f,tf)
{
PRESSURE += F_P(f,tf);
m = m + 1;
}
end_f_loop(f,tf)

PRESSURE2 = PRESSURE/m;

press = (PRESSURE2)*1e-6;
rdot = (a * pow(press,nn));

SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));

dis_x = 0.0;
dis_y = 0.0;

begin_f_loop(f,tf)
{
f_node_loop(f,tf,n)
{
v = F_NODE(f,tf,n);
if(NODE_POS_NEED_UPDATE(v))
{
NODE_POS_UPDATED(v);
x = NODE_X(v);
y = NODE_Y(v);

if(y>=0.05)
{
dis_x = 0.;
dis_y = 0.;
}
else
{
dis_x = 0.;
dis_y = rdot*dtime;
}

NV_D(disp, = ,dis_x,dis_y,0);
NV_V(NODE_COORD(v), +=, disp);
}
}
}
end_f_loop(f,tf)

printf("dis y: %f\n", dis_y);

}


After I do the basic settings; while i choose "reference value"and compute "mass flow inlet" The following error is seen in fluent console:

Node 0: Process 34456: Received signal SIGSEGV.

================================================== ============================

================================================== ============================

Node 1: Process 29852: Received signal SIGSEGV.

================================================== ============================
MPI Application rank 0 exited before MPI_Finalize() with status 2
The fl process could not be started.


what is the problem?
So could anyone can help to solve this?
Thank you.