[Ludaer]

Hello everybody!
I'm using the DEFINE_SR_RATE to specify a surface reaction rate in2d porous media.However,once I calculcate the case,the error appeared as shown below:
Error: received a fatal signal (Segmentation fault).

#include"udf.h"
#define H2O 0
#define H2 1
#define C_S 2
#define CO 3
real k1_rate(real temp)
{
return
2000*exp(-32963/temp);
}
real k2_rate(real temp)
{
return
110*exp(-8958/temp);
}
real k3_rate(real temp)
{
return
201*exp(-11524/temp);
}
real reaction_rate(cell_t c, Thread *cthread,real mw[],real yi[])
{
real concenH2O = C_R(c,cthread)*yi[H2O]/mw[H2O]*C_P(c,cthread);
real concenH2 = C_R(c,cthread)*yi[H2]/mw[H2]*C_P(c,cthread);
return
k1_rate(C_T(c,cthread))*concenH2O/(1+k2_rate(C_T(c,cthread))*pow(concenH2,0.75)+ k3_rate(C_T(c,cthread))*concenH2O)*150;
}
DEFINE_SR_RATE(myrate,f,fthread,r,mw,yi,rr)
{
*rr =
reaction_rate(F_C0(f,fthread),THREAD_T0(fthread),m w,yi);
}
I already know the macro F_C0(f,fthread),THREAD_T0(fthread) are wrong ,how to modify it,thank you all attention.

[twcp0104]

In this sort of situation what I do is fill the code with markers:

Message("1");
Message("2");

If you put these all over, between every line of code, you'll have an idea of what part is giving you problems... In general, these segmentation fault deals have something to do with using a UDMI that has not been initialized in Fluent, or a variable that is being used but doesnt have a value yet, or an array that is not long enough.... the list goes on.

[Ludaer]

Originally Posted by twcp0104

In this sort of situation what I do is fill the code with markers:

Message("1");
Message("2");

If you put these all over, between every line of code, you'll have an idea of what part is giving you problems... In general, these segmentation fault deals have something to do with using a UDMI that has not been initialized in Fluent, or a variable that is being used but doesnt have a value yet, or an array that is not long enough.... the list goes on.

thank you,i have done as you say, i already know the macro F_C0(f,fthread),THREAD_T0(fthread) have problems, so i doubt the DEFINE_SR_RATE can't apply to the porous media zone