[tengteng]

Hi!

I'm trying to access my phase-1 in a multiphase simulation. As it's written now, I get segmentation error. And I have allocated extra memory, so my guess is that I'm doing something wrong when I try to get the phase thread. I would be very grateful if someone could tell me what I'm doing wrong. The code is the following:

#define fluid_zone_ID 2

DEFINE_SOURCE(Stirrer_Fx_source_to_mom_x_eqn, c_x, thread_x, dS, eqn)
{
real source_x = 0.0;
real current_time;
real f_0_magn,f_0_times_v,v_trav_wave;

int phase_index = 2; /* Phase id in fluent */
/* The row that I think is wrong!!!*/
Thread *melt_sub_thread = THREAD_SUB_THREAD(thread_x,phase_index);

current_time = CURRENT_TIME;


/* Force compensation in moving fluid is:
f_i = f_0_i*(1-(f_0*v)/(|f_0|*v_trav_wave)) =
f_0_i*(1-(sum(f_0_j*v_j))/(sqrt(sum(f_0_j^2))*(2*tau*f))) */

/* Magnitude of stirring force at zero velocity */
f_0_magn = sqrt(pow(C_UDMI(c_x,thread_x,FORCE_X_STORAGE),2.0) +pow(C_UDMI(c_x,thread_x,FORCE_Y_STORAGE),2.0)+pow (C_UDMI(c_x,thread_x,FORCE_Z_STORAGE),2.0));
/* Scalar product of stirring force and velocity, f_0*v */
f_0_times_v = C_UDMI(c_x,thread_x,FORCE_X_STORAGE)*C_U(c_x,melt_ sub_thread)+C_UDMI(c_x,thread_x,FORCE_Y_STORAGE)*C _V(c_x,melt_sub_thread)+C_UDMI(c_x,thread_x,FORCE_ Z_STORAGE)*C_W(c_x,melt_sub_thread);
/* Travelling wave speed */
v_trav_wave = 2.0*tau*freq;

/* Complete source term S [N/m^3] to momentum equation */
source_x = C_UDMI(c_x,thread_x,FORCE_X_STORAGE)*
(1.0-f_0_times_v/(f_0_magn*v_trav_wave));
/* Source term derivative wrt dependent variable, dS/du_i */
dS[eqn]= -pow(C_UDMI(c_x,thread_x,FORCE_X_STORAGE),2.0)/
(f_0_magn*v_trav_wave);


C_UDMI(c_x,thread_x,APPLIED_FORCE_X_STORAGE)=sourc e_x;
return source_x;
}

/Elin