|
[Sponsors] |
No base point for face... produces a valid tet decomposition gmsh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 18, 2023, 15:05 |
No base point for face... produces a valid tet decomposition gmsh
|
#1 |
Member
Join Date: Aug 2012
Posts: 40
Rep Power: 13 |
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}; //+ |
|
|
|
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 |