[Gmsh] Small characteristic length -> Bad mesh ?
Christophe Geuzaine
c.geuzaine at ulg.ac.be
Thu Jan 20 17:41:13 CET 2005
Taras Danilchenko wrote:
> 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 ?
Hi Taras - Could you send the parameter values (R_of_pa, H_of_pa, etc.)?
>
> 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[]};
> _______
>
>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
>
>
--
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org