[Gmsh] Small characteristic length -> Bad mesh ?
Taras Danilchenko
taras at ton.net.ru
Thu Jan 20 09:29:18 CET 2005
Hello,
I'm using GMSH which I find a very
powerful mesh generator. Actually with 2D meshing , I have not special
problems, but when I tried to make mesh for 2d surface with BSpline
boundary, I noticed that sometimes it does not work correctly,
especially when I reduce characteristic length.
AxiSym_Cyl_Jet.geo file give a good mesh only with h1 >= 0.002
Is it a bag, or the small characteristic length always leads to bad mesh ?
Thanks in advance.
Taras Danilchenko
SRI SET
St-Petersburg
Russia
____________
/* Gmsh AxiSym_Cyl_Jet.geo file. Written by Danilchenko T.Y. taras at ton.net.ru*/
p1=newp; Point(p1) = { 0 ,0 ,0, h};
p2=newp; Point(p2) = { R_of_pa ,0 ,0, h};
p3=newp; Point(p3) = { R_of_pa ,H_of_pa ,0, h};
p4=newp; Point(p4) = { 0 , H_of_pa ,0, h};
p5=newp; Point(p5) = { R_c ,0 ,0, h1};
p6=newp; Point(p6) = { R_c ,h_c ,0,h1 };
p7=newp; Point(p7) = { R_I*Cos(Angle_I) ,h_c+R_c+R_I*Sin(Angle_I) ,0,h1};
p8=newp; Point(p8) = { 0 ,h_c+R_c+R_I ,0, h1};
p9=newp; Point(p9) = { 0 ,h_c ,0, h};
p10=newp; Point(p10) = { 0 ,h_c+R_c ,0, h1};
p11=newp; Point(p11) = { 0 ,h_c+R_c+R_I/3 ,0, h1};
p12=newp; Point(p12) = { 0 ,h_c+R_c+R_I/3+R_I/4 ,0, h1};
p13=newp; Point(p13) = { 0 ,h_c+R_c+R_I/3-R_I/4 ,0, h1};
p14=newp; Point(p14) = { R_I/16 ,h_c+R_c+R_I/3 ,0, h1};
p15=newp; Point(p15) = {1.08*R_I*Cos(Pi/4+Angle_I/2) ,h_c+R_c+1.08*R_I*Sin(Pi/4+Angle_I/2) ,0,h1};
Line(1) = {1,5};
Line(2) = {5,2};
Line(3) = {2,3};
Line(4) = {3,4};
Line(5) = {4,8};
Line(6) = {8,12};
Ellipse(61) = {12,11,12,14};
Ellipse(62) = {14,11,12,13};
Line(63) = {13,10};
Line(64) = {13,12};
Line(7) = {10,9};
Line(8) = {9,1};
Line(10) = {6,5};
Circle(11) = {6,9,10};
BSpline(19) = {8,8,8,15,7,7,10,10,10};
Line Loop(65) = {2,3,4,5,19,-11,10};
Plane Surface(9000010) = {65};
Line Loop(9000011) = {-19,6,61,62,63};
Plane Surface(9000012) = {9000011};
Line Loop(9000013) = {62,64,61};
Plane Surface(9000014) = {9000013};
Tank_Surf[] = {2,3,4,11,10};
Diel1 = 9000001;Physical Surface(Diel1)={9000010};
Diel2 = 9000002;Physical Surface(Diel2)={9000012};
Diel3 = 9000003;Physical Surface(Diel3)={9000014};
Ground = 9000009;Physical Line(Ground)={Tank_Surf[]};
_______