I made a .geo file with gmsh 2.2.3 until it produced a mesh a approximately
representing what I had in mind. But when I used the same geo file with gmsh
2.2.4 I got a totally different mesh ( well the geometry remains the same but
the mesh density globally and locally different. See these gmsh terminal output
for a first idea.
Both versions were compiled from source with GCC.
Gmsh version 2.2.3
gmsh -rand 1.e-8 circle2.geo -2
Info : 'gmsh -rand 1.e-8 circle2.geo -2' started on Thu Sep 25 19:02:26 2008
Info : Reading 'circle2.geo'
Info : Read 'circle2.geo'
Info : Meshing 1D...
Info : Meshing curve 1 (Line)
Info : Meshing curve 2 (Line)
Info : Meshing curve 3 (Line)
Info : Meshing curve 4 (Line)
Info : Meshing curve 5 (Line)
Info : Meshing curve 6 (Line)
Info : Meshing curve 7 (Line)
Info : Meshing curve 8 (Line)
Info : Meshing curve 9 (Line)
Info : Meshing curve 10 (Line)
Info : Meshing curve 11 (Line)
Info : Meshing curve 12 (Line)
Info : Meshing curve 13 (Line)
Info : Meshing curve 14 (Line)
Info : Meshing curve 15 (Line)
Info : Meshing curve 16 (Line)
Info : Meshing curve 100 (Circle)
Info : Meshing curve 200 (Circle)
Info : Mesh 1D complete (0.504031 s)
Info : Mesh
Info : Meshing 2D...
Info : Meshing surface 106 (Plane, MeshAdapt+Delaunay)
Info : Mesh 2D complete (65.8241 s)
Info : Mesh
Info : 107605 vertices 215214 elements
Info : Writing 'circle2.msh'
Info : Wrote 'circle2.msh'
Gmsh version 2.2.4:
gmsh -rand 1e-8 circle2.geo -2
Info : Reading 'circle2.geo'
Info : Read 'circle2.geo'
Info : Meshing 1D...
Info : Meshing curve 1 (Line)
Info : Meshing curve 2 (Line)
Info : Meshing curve 3 (Line)
Info : Meshing curve 4 (Line)
Info : Meshing curve 5 (Line)
Info : Meshing curve 6 (Line)
Info : Meshing curve 7 (Line)
Info : Meshing curve 8 (Line)
Info : Meshing curve 9 (Line)
Info : Meshing curve 10 (Line)
Info : Meshing curve 11 (Line)
Info : Meshing curve 12 (Line)
Info : Meshing curve 13 (Line)
Info : Meshing curve 14 (Line)
Info : Meshing curve 15 (Line)
Info : Meshing curve 16 (Line)
Info : Meshing curve 100 (Circle)
Info : Meshing curve 200 (Circle)
Info : Mesh 1D complete (0.43 s)
Info : Mesh
Info : Meshing 2D...
Info : Meshing surface 106 (Plane, MeshAdapt+Delaunay)
Info : Mesh 2D complete (365.69 s)
Info : Mesh
Info : 750256 vertices 1500535 elements
Info : Writing 'circle2.msh'
Info : Wrote 'circle2.msh'
My questions are:
Can somebody tell me the reason for this?
Did the syntax change so that what I was using is not longer valid in 2.2.4?
Is there a Bug responsible for this. And if in which version?
Can somebody reproduce my findings ( I attached the .geo file )?
Or must be something wrong with my installation. ( although I could reproduce
this already on 2 machines).
kind regards
Moritz
//Abmessungen
//Radius des Behälters
r = 6E-2; // war mal 9 cm
//Länge eines Flügels
l = 2.32E-2;
//Dicke eines Flügels
d = 2.5E-4;
//Abstand der Flügel vom Ursprung
a = 0.94E-2;
//Elementauflösung der Punkte
/*
vHighPr = 0.00002;
highPr = 0.00004;
midPr = 0.00015;
lowPr = 0.0006;
vLowPr = 0.0006;
*/
vHighPr = 0.000018;
highPr = 0.00004;
midPr = 0.0002;
lowPr = 0.0006;
vLowPr = 0.001;
Point(0) = {0,0,0}; // Ursprung
//rechter Flügel
Point(1) = {a,-0.5*d,0,midPr};
Point(2) = {a+l,-0.5*d,0,midPr};
Point(3) = {a+l,0.5*d,0,midPr};
Point(4) = {a,0.5*d,0,midPr};
Line(1) = {1,2}; Physical Line (1) = {1};
Line(2) = {2,3}; Physical Line (2) = {2};
Line(3) = {3,4}; Physical Line (3) = {3};
Line(4) = {4,1}; Physical Line (4) = {4};
Line Loop(101) = {1,2,3,4};
//linker Flügel
Point(5) = {-a,0.5*d,0,midPr};
Point(6) = {-(a+l),0.5*d,0,midPr};
Point(7) = {-(a+l),-0.5*d,0,midPr};
Point(8) = {-a,-0.5*d,0,midPr};
Line(5) = {5,6};
Line(6) = {6,7};
Line(7) = {7,8};
Line(8) = {8,5};
Physical Line (5) = {5};
Physical Line (6) = {6};
Physical Line (7) = {7};
Physical Line (8) = {8};
Line Loop(102) = {5,6,7,8};
//oberer Flügel
Point(9) = {-0.5*d,a,0,midPr};
Point(10) = {-0.5*d,a+l,0,midPr};
Point(11) = {0.5*d,a+l,0,midPr};
Point(12) = {0.5*d,a,0,midPr};
Line(9) = {9,10};
Line(10) = {10,11};
Line(11) = {11,12};
Line(12) = {12,9};
Physical Line (9) = {9};
Physical Line (10) = {10};
Physical Line (11) = {11};
Physical Line (12) = {12};
Line Loop(103) = {9,10,11,12};
//unterer Flügel
Point(13) = {0.5*d,-a,0,midPr};
Point(14) = {0.5*d,-(a+l),0,midPr};
Point(15) = {-0.5*d,-(a+l),0,midPr};
Point(16) = {-0.5*d,-a,0,midPr};
Line(13) = {13,14};
Line(14) = {14,15};
Line(15) = {15,16};
Line(16) = {16,13};
Physical Line (13) = {13};
Physical Line (14) = {14};
Physical Line (15) = {15};
Physical Line (16) = {16};
Line Loop(104) = {13,14,15,16};
//AuÃenwand
Point(100) = {-r,0,0,vLowPr};
Point(200) = {r,0,0,vLowPr};
Circle(100) = {100,0,200};
Circle(200) = {200,0,100};
Physical Line (17) = {100,200};
Line Loop(105) = {100,200};
Plane Surface(106) = {101,102,103,104,105};
//lables for Elmer
Physical Surface ( 1 ) = { 106 };
//lokale Auflösung beeinflussen
//feine Auflösung für die Flügelecken
Field[1] = Attractor;
Field[1].EdgesList = {2,4,6,8,10,12,14,16};
Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = vHighPr;
Field[2].LcMax = vLowPr;
Field[2].DistMin = 0.025E-2;
Field[2].DistMax = 0.25E-2;
Field[5] = Attractor;
Field[5].EdgesList = {2,4,6,8,10,12,14,16};
Field[6] = Threshold;
Field[6].IField = 5;
Field[6].LcMin = highPr;
Field[6].LcMax = vLowPr;
Field[6].DistMin = 0.1E-2;
Field[6].DistMax = 1E-2;
// Auflösung der Flügel
Field[3] = Attractor;
Field[3].EdgesList = {1,3,5,7,9,11,13,15};
Field[4] = Threshold;
Field[4].IField = 3;
Field[4].LcMin = midPr;
Field[4].LcMax = vLowPr;
Field[4].DistMin = 0.5E-2;
Field[4].DistMax = 2E-2;
//Backround
Field[10] = Min;
Field[10].FieldsList = {2,4,6};
Background Field = 10;
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
