Check if face is Boundary

HorstvanGrass · October 31, 2018, 07:32

Hi,
I'm looping over cells and than of each cell I'm looping over the faces. Each time I got a face I want to check if it's a boundary zone face. Is there a macro, which is resolving this?
d = Get_Domain(1); /* mixture domain if multiphase */
thread_loop_c(t,d)
{

begin_c_loop(c,t)

check_if_boundary(f,t);
c_face_loop(c, t, n)
{
f = C_FACE(c,t,n);

faces[n]=f;
}

obscureed · October 31, 2018, 17:02

Hi Yann,

You won't get very far with "f=C_FACE(c,t,n);" without "ft=C_FACE_THREAD(c,t,n);" as well. When you have the thread ft (which is of type Thread*), you can test that using BOUNDARY_FACE_THREAD_P(ft) -- check it out in the UDF manual.

You can get other information about the thread type -- see for example In UDF, how to judge the type of boundary? . Remember that each face of a cell could be in a different thread.

Good luck!
Ed

HorstvanGrass · November 2, 2018, 10:26

Thx, it worked!

October 31, 2018, 07:32	Check if face is Boundary	#1
HorstvanGrass Member Horst van Gras Join Date: Oct 2018 Posts: 45 Rep Power: 8	Hi, I'm looping over cells and than of each cell I'm looping over the faces. Each time I got a face I want to check if it's a boundary zone face. Is there a macro, which is resolving this? d = Get_Domain(1); /* mixture domain if multiphase */ thread_loop_c(t,d) { begin_c_loop(c,t) check_if_boundary(f,t); c_face_loop(c, t, n) { f = C_FACE(c,t,n); faces[n]=f; }

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Radiation in semi-transparent media with surface-to-surface model?	mpeppels	CFX	11	August 22, 2019 08:30
Question about heat transfer coefficient setting for CFX	Anna Tian	CFX	1	June 16, 2013 07:28
[blockMesh] non-orthogonal faces and incorrect orientation?	nennbs	OpenFOAM Meshing & Mesh Conversion	7	April 17, 2013 06:42
[Netgen] Import netgen mesh to OpenFOAM	hsieh	OpenFOAM Meshing & Mesh Conversion	32	September 13, 2011 06:50
Problems in compiling paraview in Suse 10.3 platform	chiven	OpenFOAM Installation	3	December 1, 2009 08:21

October 31, 2018, 17:02		#2
obscureed Senior Member Join Date: Sep 2017 Posts: 246 Rep Power: 12	Hi Yann, You won't get very far with "f=C_FACE(c,t,n);" without "ft=C_FACE_THREAD(c,t,n);" as well. When you have the thread ft (which is of type Thread*), you can test that using BOUNDARY_FACE_THREAD_P(ft) -- check it out in the UDF manual. You can get other information about the thread type -- see for example In UDF, how to judge the type of boundary? . Remember that each face of a cell could be in a different thread. Good luck! Ed

November 2, 2018, 10:26		#3
HorstvanGrass Member Horst van Gras Join Date: Oct 2018 Posts: 45 Rep Power: 8	Thx, it worked!