CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Open Source Meshers: Gmsh, Netgen, CGNS, ...

Hex with transfinite Volume

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 3, 2010, 10:00
Default Hex with transfinite Volume
  #1
New Member
 
Join Date: Sep 2009
Posts: 15
Rep Power: 0
Mjoelnir is on a distinguished road
Hi everybody,

i'm trying to mesh a section of sphere.

Code:
 
// Convert to meters
toM=1;
// Characteristic length
lc = 1*toM;
//Height and radius
//Aspect ratio Teta=D/h
D=0.29;
R=D*0.5*toM;
h=R*toM;
r=R/1.5*toM;
w=R*0.5;
W=R*0.6;
 
Ru=Sqrt(R*R-W*W);
Point(1) = {0,0,0,lc};
Point(2) = {R,0,0,lc};
Point(3) = {0,R,0,lc};
Point(5) = {r,0,0,lc};
Point(4) = {0,r,0,lc};
Point(6) = {Ru,0,-W,lc};
Point(7) = {0,Ru,-W,lc};
Point(8) = {r,0,-w,lc};
Point(9) = {0,r,-w,lc};
 
Circle (1) = {2,1,3};
Line (2) = {3,4};
Line (3) = {4,5};
Line (4) = {5,2};
Circle (5) = {2,1,6};
Circle (6) = {3,1,7};
Line (7) = {4,9};
Line (8) = {5,8};
Circle (9) = {6,1,7};
Line (10) = {7,9};
Line (11) = {9,8};
Line (12) = {8,6};
 
nodes1=11;
nodes2=5;
Transfinite Line {1,3,5,6,7,8,9,10,11} =nodes1;
Transfinite Line {2,4,10,12} =nodes2;
Line Loop (1) ={1,2,3,4};
Plane Surface (1)={1};
Line Loop (2) ={5,9,-6,-1};
Ruled Surface (2)={2};
Line Loop (3) ={6,10,-7,-2};
Plane Surface (3)={3};
Line Loop (4) ={7,11,-8,-3};
Plane Surface (4)={4};
Line Loop (5) ={8,12,-5,-4};
Plane Surface (5)={5};
Line Loop (6) ={9,10,11,12};
Plane Surface (6)={6};
Transfinite Surface {1} = {2,3,4,5};
Transfinite Surface {2} = {2,3,6,7};
Transfinite Surface {3} = {3,7,9,4};
Transfinite Surface {4} = {4,5,8,9};
Transfinite Surface {5} = {2,5,6,8};
Transfinite Surface {6} = {6,7,8,9};
Recombine Surface {1,2,3,4,5,6};
 
Surface Loop (1) = {1,2,3,4,5,-6} ;
Volume (1) = {1};
 
 
Transfinite Volume {1} = {2,3,4,5,6,7,8,9};
 
 
//Physical Surface("Wall") = {2,3,4,5,37,54,59,76,81,102,103,120,125} ;
 
//Physical Volume ("internalField") = {1};

i get following error Message:

Error : Incompatible surface 3 in transfinite volume 1
Error : Cannot tetrahedralize volume with quadrangles on boundary

i've no idea what could cause this problems


Thanks you
Mjoelnir is offline   Reply With Quote

Old   February 4, 2010, 00:59
Default
  #2
Senior Member
 
Eric Nutsch
Join Date: Sep 2009
Location: Eugene, Oregon USA
Posts: 112
Rep Power: 6
ericnutsch is on a distinguished road
Send a message via Skype™ to ericnutsch
The error message sounds fairly intuitive; Squares and triangles don't mix on surfaces, but hard to say without seeing it.

If you have tets' you may want to move their definitions to end of the script and put the hexes' at the beginning(long shot, but its worth a try) I will need some print screens of your problem to be of any further help.

You can attach photos with the "mange attachments" button below.


Best of Luck!
ericnutsch is offline   Reply With Quote

Old   February 9, 2010, 11:54
Default Pictures
  #3
New Member
 
Join Date: Sep 2009
Posts: 15
Rep Power: 0
Mjoelnir is on a distinguished road
If i don't use the command Recombine Surface, i can build a 3d Mesh
(see Mesh2dtirangle; Mesh3dquadangle)

Without this command i can't create an 3d Mesh.

(see Mesh3dquadangle, Mesh3dtriangle)

thanks

Mjoelnir
Attached Images
File Type: jpg Mesh2dtriangle.JPG (80.7 KB, 47 views)
File Type: jpg Mesh3dtriangle.JPG (95.7 KB, 46 views)
File Type: jpg Mesh2dquadangle.JPG (70.5 KB, 45 views)
File Type: jpg Mesh3dquadangle.JPG (87.5 KB, 43 views)
Mjoelnir is offline   Reply With Quote

Old   February 17, 2010, 10:55
Default Problem solved
  #4
New Member
 
Join Date: Sep 2009
Posts: 15
Rep Power: 0
Mjoelnir is on a distinguished road
i figured out the Problem. in the transfinite Volume command the order of the points ist important.
Mjoelnir is offline   Reply With Quote

Old   April 1, 2010, 19:00
Default
  #5
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 6
robingilbert is on a distinguished road
hi i have a similar problem!!
what order of points are we supposed to use??
thank you
robingilbert is offline   Reply With Quote

Old   April 8, 2010, 08:29
Default
  #6
New Member
 
Join Date: Sep 2009
Posts: 15
Rep Power: 0
Mjoelnir is on a distinguished road
its like in Blockmesh

http://www.openfoam.com/docs/user/me...23-1270005.1.1

for hexandron
Mjoelnir is offline   Reply With Quote

Old   April 8, 2010, 13:52
Default
  #7
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 6
robingilbert is on a distinguished road
oh ya!!! thanks a lot..
thank you so much for the help!!
robingilbert is offline   Reply With Quote

Old   May 12, 2014, 18:00
Default Incompatible surface in transfinite volume
  #8
New Member
 
Join Date: Sep 2013
Posts: 5
Rep Power: 2
menonshyam is on a distinguished road
Hello, I am getting a very similar error and was wondering if you could spot what is wrong in the numbering of the volume! The error i get is "Incompatible surface 22 in transfinite volume 28". Here is my code:

lc = 0.005;
length = 0.01; //chamber radius/sqrt(2)

Point(1) = {0,0,0,lc};
Point(2) = {length,length,-length,lc};
Point(3) = {-length,length,-length,lc};
Point(4) = {-length,-length,-length,lc};
Point(5) = {length,-length,-length,lc};

Circle(6) = {3,1,2};
Circle(7) = {2,1,5};
Circle(8) = {5,1,4};
Circle(9) = {4,1,3};

Line(13) = {4, 1};
Line(14) = {5, 1};
Line(15) = {2, 1};
Line(16) = {3, 1};

Line Loop(17) = {13, -16, -9};
Ruled Surface(18) = {17};
Line Loop(19) = {16, -15, -6};
Ruled Surface(20) = {19};
Line Loop(21) = {15, -14, -7};
Ruled Surface(22) = {21};
Line Loop(23) = {14, -13, -8};
Ruled Surface(24) = {23};
Line Loop(25) = {8, 9, 6, 7};
Ruled Surface(26) = {25};
Surface Loop(27) = {18, 22, 26, 24, 20};
Volume(28) = {27};
Transfinite Line {13, 14, 15, 16} = 10 Using Progression 1;

Woutlet = 26;
Wfront = 18;
Wback = 22;
Wtop = 20;
Wbottom = 24;
WinternalVol = 28;

Physical Volume("internalVol") = {WinternalVol};
Physical Surface("outlet") = {Woutlet};
Physical Surface("front") = {Wfront};
Physical Surface("back") = {Wback};
Physical Surface("top") = {Wtop};
Physical Surface("bottom") = {Wbottom};

Transfinite Surface "*";
Recombine Surface "*";
Transfinite Volume "*";


Thanks!
menonshyam is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
CheckMeshbs errors ivanyao OpenFOAM Running, Solving & CFD 2 March 11, 2009 03:34
On the damBreak4phaseFine cases paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 21:14
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 11:55
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00


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