[ImmersedBoundary] Problems when running potentialFoam in foam-extend 4.1 by using IBM

YatengQ · December 19, 2022, 04:28

Hi everyone,

I'm going to simulation a 3D helix rotating in the fluid domain by using IBM in foam-extend 4.1.

I used movingCylinderInChannelIcoFine case for reference, but in the Z direction the mesh cell is set to be 1 in blockMeshDict. When I set 100 cells in Z direction. The error appears in the log.potentialFoam as follows:

Create mesh for time = 0
153
154
155
156 --> FOAM FATAL ERROR:
157 Error in point ordering: mixed used and unused points at the end of point list.158 Last used point: 0 (-0.02 -0.02 0.005)
159 First unused point: 1 (-0.0198667 -0.02 0.005)
160 and point 1 (-0.0198667 -0.02 0.005) is used by a live face.
161 Face 0 4(1 302 60803 60502) with points 4((-0.0198667 -0.02 0.005) (-0.0198667
162 Done.
163
164 From function void polyMesh::initMesh()
165 in file meshes/polyMesh/polyMeshInitMesh.C at line 166.
166
167 FOAM aborting

I'm not sure if the IBM in foam-extend can simulate the 3D case. I am struggled with this error for several days.

hjasak · December 28, 2022, 22:34

Looks like a broken mesh; can you try to run something without IBM?

Hrv

YatengQ · February 11, 2023, 23:04

Hi hjasak,

Thanks for your reply. I run blockMesh, when I set the Z direction mesh to 1, the potentialFoam run correctly, when I set the Z direction mesh >1, the same error appears.

Thanks!

December 19, 2022, 04:28	Problems when running potentialFoam in foam-extend 4.1 by using IBM	#1
YatengQ New Member Yateng Qiao Join Date: Nov 2022 Posts: 13 Rep Power: 4	Hi everyone, I'm going to simulation a 3D helix rotating in the fluid domain by using IBM in foam-extend 4.1. I used movingCylinderInChannelIcoFine case for reference, but in the Z direction the mesh cell is set to be 1 in blockMeshDict. When I set 100 cells in Z direction. The error appears in the log.potentialFoam as follows: Create mesh for time = 0 153 154 155 156 --> FOAM FATAL ERROR: 157 Error in point ordering: mixed used and unused points at the end of point list.158 Last used point: 0 (-0.02 -0.02 0.005) 159 First unused point: 1 (-0.0198667 -0.02 0.005) 160 and point 1 (-0.0198667 -0.02 0.005) is used by a live face. 161 Face 0 4(1 302 60803 60502) with points 4((-0.0198667 -0.02 0.005) (-0.0198667 162 Done. 163 164 From function void polyMesh::initMesh() 165 in file meshes/polyMesh/polyMeshInitMesh.C at line 166. 166 167 FOAM aborting I'm not sure if the IBM in foam-extend can simulate the 3D case. I am struggled with this error for several days.

December 28, 2022, 22:34		#2
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,907 Rep Power: 33	Looks like a broken mesh; can you try to run something without IBM? Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Periodic B.C. in Foam extend IBM	Jiashun	OpenFOAM	0	March 12, 2022 07:58
[CGNS] CGNS converters available	mbeaudoin	OpenFOAM Meshing & Mesh Conversion	137	December 14, 2018 05:20
[swak4Foam] undefined symbols in .so libraries when compiling swak4foam 0.4.2 on foam extend 4.1	mAlletto	OpenFOAM Community Contributions	3	August 29, 2018 07:02
[blockMesh] BlockMesh FOAM warning	gaottino	OpenFOAM Meshing & Mesh Conversion	7	July 19, 2010 15:11
Statically Compiling OpenFOAM Issues	herzfeldd	OpenFOAM Installation	21	January 6, 2009 10:38

February 11, 2023, 23:04		#3
YatengQ New Member Yateng Qiao Join Date: Nov 2022 Posts: 13 Rep Power: 4	Hi hjasak, Thanks for your reply. I run blockMesh, when I set the Z direction mesh to 1, the potentialFoam run correctly, when I set the Z direction mesh >1, the same error appears. Thanks!