CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Special Topics > Mesh Generation & Pre-Processing

No base point for face... produces a valid tet decomposition gmsh

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 18, 2023, 15:05
Default No base point for face... produces a valid tet decomposition gmsh
  #1
Member
 
Join Date: Aug 2012
Posts: 40
Rep Power: 13
cesarjets is on a distinguished road
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};
//+
cesarjets is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] "no base point for face xxx - produces a valid tet decomposition" eugene.katsevman OpenFOAM Meshing & Mesh Conversion 18 April 7, 2017 11:35
[snappyHexMesh] How to define to right point for locationInMesh Mirage12 OpenFOAM Meshing & Mesh Conversion 7 March 13, 2016 14:07
rhoSImplecFoam still having issues! mihaipruna OpenFOAM Running, Solving & CFD 10 August 12, 2012 10:30
interFoam & decomposition method: scotch MacGyver OpenFOAM Running, Solving & CFD 2 May 23, 2012 07:00
OF 1.6 | Ubuntu 9.10 (64bit) | GLIBCXX_3.4.11 not found piprus OpenFOAM Installation 22 February 25, 2010 13:43


All times are GMT -4. The time now is 16:12.