[cfdmatteo]

hi there,

I am learning to use DEFINE_UDS_FLUX macro, and I follow the UDF Manual
example but somehow there's a problem :

Here's my code:

#include "udf.h"

DEFINE_UDS_FLUX(ws_flux,f,t,i)
{
cell_t cO, c1=-1;
Thread *t0, *t1 = NULL;

real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A,f,t);

if(BOUNDARY_FACE_THREAD_P(T))
{

real dens;
if(NNULP(THREAD_STORAGE(t,SV_DENSITY)))
dens=F_R(f,t);
else
dens=C_R(c0,t0);

NV_DS(psi_vec,=,F_U(f,t),F_V(f,t),F_W(f,t),*,dens) ;
flux=NV_DOT(psi_vec,A);
}
else
{
c1=F_C1(f,t);
t1=F_C_THREAD(f,t);

NV_DS(psi_vec,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),* ,C_R(c0,t0));
NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,C_R(c1,t1));
flux=NV_DOT(psi_vec,A)/2.0;
}
return flux;
}

And after I interpreting this macro, the command shows this error:

Error: D:\ryan\source_code\ws_flux.c: line 11: THREAD_STORE: undeclared variable

Can anybody help me? Any suggestion will be appreciated!

best regard, Ryan

[MMMSim]

Hi Ryan,

I think you wrote a typo in line 11:
t0 = F_C0_THEREAD(f,t);
try
t0 = F_C0_THREAD(f,t);
instead

good luck!

[cfdmatteo]

Originally Posted by MMMSim

Hi Ryan,

I think you wrote a typo in line 11:
t0 = F_C0_THEREAD(f,t);
try
t0 = F_C0_THREAD(f,t);
instead

good luck!

Dear Andreas K,
Thanks for your suggestion. : )
I'v fixed the code but still getting with the same problem.
It seems the DEFINE_UDS_FLUX example of the Manual is of defect.
Since I directly copy the code of example and somehow the same
error comes up.

[cfdmatteo]

Originally Posted by cfdmatteo

hi there,

I am learning to use DEFINE_UDS_FLUX macro, and I follow the UDF Manual
example but somehow there's a problem :

Here's my code:

#include "udf.h"

DEFINE_UDS_FLUX(ws_flux,f,t,i)
{
cell_t cO, c1=-1;
Thread *t0, *t1 = NULL;

real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A,f,t);

if(BOUNDARY_FACE_THREAD_P(T))
{

real dens;
if(NNULP(THREAD_STORAGE(t,SV_DENSITY)))
dens=F_R(f,t);
else
dens=C_R(c0,t0);

NV_DS(psi_vec,=,F_U(f,t),F_V(f,t),F_W(f,t),*,dens) ;
flux=NV_DOT(psi_vec,A);
}
else
{
c1=F_C1(f,t);
t1=F_C_THREAD(f,t);

NV_DS(psi_vec,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),* ,C_R(c0,t0));
NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,C_R(c1,t1));
flux=NV_DOT(psi_vec,A)/2.0;
}
return flux;
}

And after I interpreting this macro, the command shows this error:

Error: D:\ryan\source_code\ws_flux.c: line 11: THREAD_STORE: undeclared variable

Can anybody help me? Any suggestion will be appreciated!

best regard, Ryan

Well, the first problem seems to be solved since I subtitute
t0 = F_C0_THREAD(f,t) with t0 = THREAD_T0(t).

But here comes other problems :

for line 14 if(BOUNDARY_FACE_THREAD_P(T)) :
Error: D:\ryan\source_code\ws_flux.c: line 14: parse error.

and for line 20 dens = F_R(f,t); :
Error: D:\ryan\source_code\ws_flux.c: line 20: dens: undeclared variable

Any suggestion?
Ryan

[chenxizh]

Hi, you have serval mistakes in your udf, here the correct one is following, any way do you used it in multiphase?

#include "udf.h"

DEFINE_UDS_FLUX(ws_flux,f,t,i)
{
cell_t c0, c1=-1;
Thread *t0, *t1 = NULL;

real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A,f,t);

if(BOUNDARY_FACE_THREAD_P(t))
{

real dens;
if(NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
dens=F_R(f,t);
else
dens=C_R(c0,t0);

NV_DS(psi_vec,=,F_U(f,t),F_V(f,t),F_W(f,t),*,dens) ;
flux=NV_DOT(psi_vec,A);
}
else
{
c1=F_C1(f,t);
t1=F_C1_THREAD(f,t);

NV_DS(psi_vec,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),* ,C_R(c0,t0));
NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,C_R(c1,t1));
flux=NV_DOT(psi_vec,A)/2.0;
}
return flux;
}

Originally Posted by cfdmatteo

Well, the first problem seems to be solved since I subtitute
t0 = F_C0_THREAD(f,t) with t0 = THREAD_T0(t).

But here comes other problems :

for line 14 if(BOUNDARY_FACE_THREAD_P(T)) :
Error: D:\ryan\source_code\ws_flux.c: line 14: parse error.

and for line 20 dens = F_R(f,t); :
Error: D:\ryan\source_code\ws_flux.c: line 20: dens: undeclared variable

Any suggestion?
Ryan

[cfdmatteo]

Originally Posted by chenxizh

Hi, you have serval mistakes in your udf, here the correct one is following, any way do you used it in multiphase?

#include "udf.h"

DEFINE_UDS_FLUX(ws_flux,f,t,i)
{
cell_t c0, c1=-1;
Thread *t0, *t1 = NULL;

real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A,f,t);

if(BOUNDARY_FACE_THREAD_P(t))
{

real dens;
if(NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
dens=F_R(f,t);
else
dens=C_R(c0,t0);

NV_DS(psi_vec,=,F_U(f,t),F_V(f,t),F_W(f,t),*,dens) ;
flux=NV_DOT(psi_vec,A);
}
else
{
c1=F_C1(f,t);
t1=F_C1_THREAD(f,t);

NV_DS(psi_vec,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),* ,C_R(c0,t0));
NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,C_R(c1,t1));
flux=NV_DOT(psi_vec,A)/2.0;
}
return flux;
}

hi
Thanks a lot for your help.
I am trying to calculate a multiphase problem. And as following your fixed code, it's still get some problem. Anyway, some say that might be a "interpreted or compiled" problem. Did you read this code by "interpreted" or by "compiled" in Fluent?

I was interpreting it and always got errors as :

Error: D:\ryan\paper comfirm\cfd_help_ws_flux.c: line 11: THREAD_STORE:undeclared variable

[chenxizh]

I use complie, there is no problem in my code
notice,the scalar transport in multiphase contain the volume fraction

Originally Posted by cfdmatteo

hi
Thanks a lot for your help.
I am trying to calculate a multiphase problem. And as following your fixed code, it's still get some problem. Anyway, some say that might be a "interpreted or compiled" problem. Did you read this code by "interpreted" or by "compiled" in Fluent?

I was interpreting it and always got errors as :

Error: D:\ryan\paper comfirm\cfd_help_ws_flux.c: line 11: THREAD_STORE:undeclared variable

[cfdmatteo]

Originally Posted by chenxizh

I use complie, there is no problem in my code
notice,the scalar transport in multiphase contain the volume fraction

hi,
I appreciate for your reply, it's really helpful for me.

[chenxizh]

what is the problem you are studying? I aslo used uds in multiphase, and I find that the uds is unbound.

Originally Posted by cfdmatteo

hi,
I appreciate for your reply, it's really helpful for me.

[cfdmatteo]

Originally Posted by chenxizh

what is the problem you are studying? I aslo used uds in multiphase, and I find that the uds is unbound.

hi,
I study on fuel cell calculating, and I use the PEMFC module of the Fluent but somehow the model is of defect.
I think the multiphase model of this module is UFT model,but the flux term is neglected, so I try to fixed the problem by UDF.

[juzer_700]

Is there a predefined macro under which I can write a new model (combustion model). I dont plan to use any combustion model provided in FLUENT.

Please help. I want to link the advance combustion model to fluent.

Under which DEFINE macro should I write my code???

Any help will be appreciate.

[sunjian]

Originally Posted by chenxizh

Hi, you have serval mistakes in your udf, here the correct one is following, any way do you used it in multiphase?

#include "udf.h"

DEFINE_UDS_FLUX(ws_flux,f,t,i)
{
cell_t c0, c1=-1;
Thread *t0, *t1 = NULL;

real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A,f,t);

if(BOUNDARY_FACE_THREAD_P(t))
{

real dens;
if(NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
dens=F_R(f,t);
else
dens=C_R(c0,t0);

NV_DS(psi_vec,=,F_U(f,t),F_V(f,t),F_W(f,t),*,dens) ;
flux=NV_DOT(psi_vec,A);
}
else
{
c1=F_C1(f,t);
t1=F_C1_THREAD(f,t);

NV_DS(psi_vec,=,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),* ,C_R(c0,t0));
NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,C_R(c1,t1));
flux=NV_DOT(psi_vec,A)/2.0;
}
return flux;
}

Hi, chenxizh. Have you tried to calculate with this code? This code could be compiled without any error. But when iteration was run, the error access_violation occured. What do you think the problem is? I am looking forward to your reply. Thank you.

[xuansy]

Hello, sunjian. I have the same problem, can you tell me how to solve this problem. Thank you.

[xuansy]

Originally Posted by sunjian

Hi, chenxizh. Have you tried to calculate with this code? This code could be compiled without any error. But when iteration was run, the error access_violation occured. What do you think the problem is? I am looking forward to your reply. Thank you.

Hello, sunjian. I have the same problem, can you tell me how to solve this problem. Thank you.

[toolpost]

Hi xuan. I am in same situation, but the error message is different here. My case is serial run on a Linux system. It compiles fine, but crashes in first iteration itself.

  Node 0: Process 8824: Received signal SIGSEGV.
  ==============================================================================
 MPI Application rank 0 exited before MPI_Finalize() with status 2
  The fluent process could not be started.

Did you able to solve it?

[qin chunqiu]

Dear Shi-yao,
Have you solved this problem. In my case, the iteration is OK but results in a result contradicing to physics.

[pakk]

Can you elaborate on the results you see? In which way do they contradict physics?

This information would help on how to fix the UDF...