[Harry321]

Dear colleagues,

I have a problem with looping inside c_face_loop(c, t, i). At the beginning I should mention that I am working with two phases and the source is attached to the second phase.

I wrote simple UDF to find out where is a problem:
DEFINE_SOURCE(source, c, t, dS, eqn)
{
int i;
real a = 0;
real b = 0, cc = 0, d = 0, ee = 0, ff = 0, gg = 0;
face_t f;
Thread *t112;
c_face_loop(c, t, i)
{

f = C_FACE(c, t, i);
t112 = C_FACE_THREAD(c, t, i);
if (THREAD_TYPE(t112) != THREAD_F_WALL);
{
Message("ThreadID t112: %d\n", THREAD_ID(t112));
Message("F_W: %f\n", F_W(f, t));
}
}

dS[eqn] = cc-d;
return a -b;
}

And I have checked several cases and non of them seems to work good:

1) first case: Message("F_W: %f\n", F_W(f, t));
I run this UDF and I got several messages on the screen which shows values of F_W and after the while after fluent showed ID of a thread I got Segmentation fault.

2) second case: Message("F_W: %f\n", F_W(f, t112));

but in that case I got Segmentation fault after showing first ID.

3) third case: I added some more code befor If statement:
t11 = THREAD_SUB_THREADS(t112);
t12 = t11[0];
t13 = t11[1];
t14 = t11[2];
now Message is:
Message("F_W: %f\n", F_W(f, t13));
the same result as in case 2)

In this all cases ID of a t112 is interior so I am not a boundary.

Does anybody know how to fix this UDF? or maybe it is limitation of this macro and I cannot access face values using this one.

[Harry321]

I added also:



DEFINE_SOURCE(source, c, t, dS, eqn)
{
int i;
real a = 0;
real b = 0, cc = 0, d = 0, ee = 0, ff = 0, gg = 0;
face_t f;
Thread *t112;
Domain *superdomain= Get_Domain(1);
Domain *subdomain = DOMAIN_SUB_DOMAIN(superdomain, 1);
Thread *tm = THREAD_SUPER_THREAD(t);
c_face_loop(c, tm, i)
{
f = C_FACE(c, tm, i);
t112 = C_FACE_THREAD(c, tm, i);
Thread *t13 = Lookup_Thread(subdomain, THREAD_ID(t112));
if (THREAD_TYPE(t112) != THREAD_F_WALL);
{
Message("ThreadID t112: %d\n", THREAD_ID(t112));
Message("F_W: %f\n", F_W(f, t));
}
}
dS[eqn] = cc-d;
return a -b;

}



I run out of ideas what can I do :/

[upeksa]

Is your case 2D or 3D?

[Harry321]

The case is 3 D

[pakk]

Try what this (added Message) does in your code:

...
t112 = C_FACE_THREAD(c, tm, i);
Message("Type:%d\n",THREAD_TYPE(t112));
Thread *t13 = Lookup_Thread(subdomain, THREAD_ID(t112));
if (THREAD_TYPE(t112) != THREAD_F_WALL);
...

Check the last type printed before you get an error.
You now exclude types of THREAD_F_WALL, but maybe you should have exclude more thread types.

[Harry321]

Hi,
thank both of you for help.

Message("ThreadID t112: %d\n", THREAD_ID(t112));

This line shows me the ID of a thread appropriate for a face. It shows ID of a interior side so I am not at the any kind of BC.
I tested sever cases (including alocation of a second phase velocity in face UDS) when I try to retrieve UDS values for in c_face_loop(c, tm, i) {} in DEFINE_SOURCE by using Message("F_U: %f\n", F_UDSI(f, t13, 5)); where
Domain *superdomain = Get_Domain(1);
Domain *subdomain = DOMAIN_SUB_DOMAIN(superdomain, 1);
f = C_FACE(c, tm, i);
t112 = C_FACE_THREAD(c, tm, i);
t13 = Lookup_Thread(subdomain, THREAD_ID(t112));

I don't have a problem with cell values for two phases but it seems that FLUENT has some issues with this specific loop and especially with face values for phases.

[Harry321]

The answer to this problem is very simple. Fluent doesn't store face velocity at the interior only at the boundaries.