Olivier Allain wrote:
Dear sir,

 I am a new user of gmsh. This tool is very smart and easy to use. I
thank you for your effort. I have a question with a problem of
coherence between two surfaces. It appears in 3D using a extruded zone
 and an unstructured zone. If you look at the boundary between the
center and the external circle it appears the mesh is not conform. It
seems the mesh generate 2 volumes without taking into account there is
 a common face. It runs with Netgen but it crashes with Tetgen, and
Netgen seems to be not able to run for a big case. Is there an option
for Tetgen in order to solve this problem ?

Unfortunately, no: a 3D Delaunay algorithm cannot be guaranteed to
recover the surface triangulation (it sometimes needs to swap edges or
insert new points on the boundary).

We could try to make the meshes match at the interface a posteriori, but
that would destroy the structure of the structured mesh... It's the main
reason why we keep Netgen around.

(PS: In your particular example you could extrude the whole mesh: cf.
attached file.)



 Best Regards,
 O. Allain


------------------------------------------------------------------------

_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh


--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine



lc = 0.06;
lc2 = 3;
rayon = 0.6;
rayon2 = 12;
Point(1) = {rayon,0.0,0.0,lc};
Point(2) = {0.0,0.0,0.0,lc};
Point(3) = {-rayon,0.0,0.0,lc};
Point(4) = {-0.0,rayon,0.0,lc};
Point(5) = {-0.0,-rayon,0.0,lc};
Circle(1) = {5,2,1};
Circle(2) = {1,2,4};
Circle(3) = {4,2,3};
Circle(4) = {3,2,5};
Line Loop(5) = {2,3,4,1};
Plane Surface(6) = {5};

Point(11) = {rayon2,0.0,0.0,lc2};
Point(12) = {-rayon2,0.0,0.0,lc2};
Point(19) = {-0.0,rayon2,0.0,lc2};
Point(14) = {-0.0,-rayon2,0.0,lc2};
Circle(33) = {14,2,11};
Circle(34) = {11,2,19};
Circle(35) = {19,2,12};
Circle(36) = {12,2,14};
Line Loop(37) = {34,35,36,33};
Plane Surface(38) = {37,5};

Extrude {0,0,rayon} { 
  Surface{6, 38}; Layers {10};
}



_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh

Reply via email to