UDF problem
 

I'm a user of fuel cell addon-module, and I'm trying to get volume averaged relative humidity. and I got an error form my customized code.
My code was compiled well but when it stared to calculate, error came up.

I think the reason of the error was C_YI(c,t,i_h2o). When I deleted this term, there was no error. Can't it be used below begin_c_loop macro??

or Is dimension wrong? please help me!

below is my code

real Membrane_Conductivity(real lam, cell_t c, Thread *t)
{
Domain *d;
real volume;
real vol_tot = 0.0;
real P_sat;
real C_h2o = 0.0;
real TT;
real RH_avg = 0.0;
int i=0;
real MC_h2o;
real MW_h2o = mw[i_h2o];

real RH_av;
real T = C_T(c,t);
real Ru = 8.314;
real x = 6.0; /* doping level*/
real Ea = -619.6*x+21750;
real a = (168*x*x*x)-(6324*x*x)+(65750*x)+8460;
real b;
real sigma;

int zt = get_fc_zone_type(t->id);
/*d = Get_Domain(ROOT_DOMAIN_ID);*/
d = Get_Domain(1);

thread_loop_c(t,d)
{

begin_c_loop(c,t)
{

volume = C_VOLUME(c,t);
TT = C_T(c,t);
P_sat = Get_P_sat(TT);
MC_h2o = C_YI(c,t,i_h2o);
/*C_h2o = C_YI(c,t,i_h2o)*C_R(c,t)/mw[i_h2o];*/

C_h2o = MC_h2o*C_R(c,t)/MW_h2o;

if (zt == aca_zone || zt == cca_zone)
{
vol_tot += volume;
RH_avg += (R_Gas*TT/P_sat*C_h2o)*volume;
}

}
end_c_loop(c,t)

RH_avg /= vol_tot;
}

RH_av = RH_avg;

b = 1+(0.01704*T-4.767)*RH_av;
sigma = a*b/T*exp(-Ea/Ru/T);

return sigma;

Hi,

When you started calculation, what error was thrown exactly?

error message was like this.

unable to read the cmd header on the pmi context, Error = -1

999999 (..\src\mpsystem.c@1173): mpt_read: failed: errno = 10054

999999: mpt_read: error: read failed trying to read 4 bytes: No such file or directory
unable to read the cmd header on the pmi context, Error = -1

received suspend command for a pmi context that doesn't exist: unmatched id = 7
received suspend command for a pmi context that doesn't exist: unmatched id = 6
received suspend command for a pmi context that doesn't exist: unmatched id = 5
received suspend command for a pmi context that doesn't exist: unmatched id = 4

job aborted:
rank: node: exit code[: error message]
0: x70-jof7ksxthyb: -1073741819: process 0 exited without calling finalize
1: x70-jof7ksxthyb: -1073741819: process 1 exited without calling finalize
2: x70-jof7ksxthyb: -1073741819: process 2 exited without calling finalize
3: x70-jof7ksxthyb: -1073741819: process 3 exited without calling finalize
4: x70-jof7ksxthyb: -1073741819: process 4 exited without calling finalize
5: x70-jof7ksxthyb: -1073741819: process 5 exited without calling finalize
6: x70-jof7ksxthyb: -1073741819: process 6 exited without calling finalize
7: x70-jof7ksxthyb: -1073741819: process 7 exited without calling finalize
received suspend command for a pmi context that doesn't exist: unmatched id = 3
received suspend command for a pmi context that doesn't exist: unmatched id = 2
received suspend command for a pmi context that doesn't exist: unmatched id = 1
The Parallel FLUENT process could not be started.

I am not sure exactly about the cause. I doubt about partitioning.
You may use:
1.In FLUENT select the "Parallel" menu
2.Choose "Auto Partition"
3.Unselect "Case File"
4.Select "Metis" as the partition method
5.Choose OK

I've changed parallel menu as you said , but it still gives error message.

Even serial processing case throws error message.

When C_YI term was deleted, the case ran well....

Hi,

Not sure, but the error seems related to parallel computing only.
http://stackoverflow.com/questions/1...ontext-error-1

Also try initializing C_YI terms with 0.0