August 18, 2023, 15:05
|
No base point for face... produces a valid tet decomposition gmsh
|
#1
|
Member
Join Date: Aug 2012
Posts: 41
Rep Power: 14
|
I've trying to run this case for an OpenFOAM application problem. I generated the mesh by using gmsh. However, I'm still getting this warning which seems to be a mesh problem and can not get any results. Not sure about it.
I would like to know if someone could give any clue on how I could solve it? Any help will be appreciated. Thanks.
--> FOAM Warning : From function Foam::triFace Foam::tetIndices::faceTriIs(const Foam:  olyMesh&) const in file /apps/easybuild/software/tinkercliffs-rome/OpenFOAM/9-foss-2021a/OpenFOAM-9/src/OpenFOAM/lnInclude/tetIndicesI.H at line 76 No base point for face 7468, 3(16 40 200), produces a valid tet decomposition.
Code:
//inputs
dx=0.0015/1000;
L1=0.04/1000;
L2=0.04/1000;
L3=dx/2;
L4=2e-5;
L5=5e-6;
L6=2.5e-5;
L7=2e-6;
L8=2e-6;
//Geometry
//inlet
Point(1) = {-L1/2,L2/2,L3};
Point(2) = {-L4, L5, L3};
Point(3) = {-L6, L5, L3};
Point(4) = {-L6, -L8, L3};
Point(5) = {-L4, -L8, L3};
Point(6) = {-L4, -L4, L3};
Point(7) = {-L4, -L4, -L3};
Point(8) = {-L4, -L7, -L3};
Point(9) = {-L6, -L7, -L3};
Point(10) = {-L6, L5, -L3};
Point(11) = {-L4, L5, -L3};
Point(12) = {-L4, L4, -L3};
Point(13) = {L1/2,-L2/2,L3};
Point(14) = {L1/2,-L2/2,-L3};
//outlet
Point(15) = {L1/2,L2/2,L3};
Point(16) = {L1/2,L2/2,-L3};
//y-dir
Line(1) = {1, 2, 3, 4, 5, 6};
Line(2) = {12, 11, 10, 9, 8, 7};
Line(8) = {13, 15};
Line(9) = {14, 16};
//z-dir
Line(3) = {1, 12};
Line(4) = {6, 7};
Line(7) = {13, 14};
Line(10) = {15, 16};
//x-dir
Line(5) = {6, 13};
Line(6) = {14, 7};
Line(11) = {15, 1};
Line(12) = {16, 12};
//X
Transfinite Curve {11, 5, 12, 6} = L1/dx Using Progression 1;
//Y
Transfinite Curve {1, 2, 8, 9} = L2/dx Using Progression 1;
//Z
Transfinite Curve {3, 4, 7, 10} = 2 Using Progression 1;
//inlet
Curve Loop(1) = {3, 2, -4, -1};
//+
Surface(1) = {1};
//outlet
Curve Loop(2) = {10, -9, -7, 8};
//+
Surface(2) = {2};
//empty y-dir //front back
Curve Loop(3) = {11, 1, 5, 8};
//+
Plane Surface(3) = {3};
//+
Curve Loop(4) = {12, 2, -6, 9};
//+
Plane Surface(4) = {4};
//empty z-dir //upperlower
Curve Loop(5) = {3, -12, -10, 11};
//+
Surface(5) = {5};
//+
Curve Loop(6) = {4, -6, -7, -5};
//+
Surface(6) = {6};
//Transfinite Surface {1} Alternated;
//+
//Transfinite Surface {2} Alternated;
//+
//Transfinite Surface {5} Alternated;
//+
//Transfinite Surface {6} Alternated;
//Recombine Surface {1, 2, 5, 6};
//Transfinite
//For j In {1:2}
//Transfinite Surface {j} = {};
//EndFor
//For j In {1:2}
//Recombine Surface {j};
//EndFor
//For j In {5:6}
//Transfinite Surface {j} = {};
//EndFor
//For j In {5:6}
//Recombine Surface {j};
//EndFor
//Coherence;
//Mesh.Smoothing = 100;
//volume
Surface Loop(1) = {3, 5, 1, 4, 6, 2};
//+
Volume(1) = {1};
TransfQuadTri{1};
Physical Surface("leftBoundary") = {1};
//+
Physical Surface("rightBoundary") = {2};
//+
Physical Surface("upperBoundary") = {5};
//+
Physical Surface("lowerBoundary") = {6};
//+
Physical Surface("emptyBoundary") = {3, 4, 5, 6};
//+
Physical Volume("internalField") = {1};
//+
|
|
|