Hi to all,

is it possible to identify (by a UDF) whether or not a cell is a boundary cell, and which one of its faces belongs to the boundary?

Thank you!

Ale

You can directly check if a face belongs to a boundary looping over all faces and using the

THREAD_TYPE (t)

where t is the faces thread.

There are many values which THREAD_TYPE (t) may assume, according to which part of your domain the face belongs to.

I report some values.


Value of THREAD_TYPE (t) Part of the domain

THREAD_F_INTERIOR Face is internal to the
domain

THREAD_F_VINLET Face belongs to a velocity
inlet

THREAD_F_POUTLET Face belongs to a pressure
outlet

THREAD_F_OUTFLOW Face belongs to an outflow




You have to use a multiple if like:


if (THREAD_TYPE(t) == THREAD_F_INTERIOR) /* Interior
{
...
}

else if(THREAD_TYPE(t) == THREAD_F_POUTLET)
{
}

else if(THREAD_TYPE(t) == THREAD_F_VINLET)
{
}



Hope this helps

Hi

ap

Thanks a lot!

I found also the functions:

BOUNDARY_FACE_THREAD_P(tf)

BOUNDARY_FACE_GEOMETRY(f,tf,A,ds,es,A_by_es,dr0)

that can be useful to identify boundary cells and to get some information about them.

Regards,

Ale