Hi all,<br><br>I am trying to produce a tetrahedral mesh of half of a  thick-walled prolate spheroid. To do this, I've first generated the geometry by adapting an example I found in the gmsh wiki. My script is included below. Then I run the 3D meshing algorithm based on MeshAdapt for 2D and Delaunay for 3D.  <br>

The problem I encounter is that a few nodes that are used to generate the triangular faces are then not used to generate any of the tetrahedral elements. Is there a way to avoid this? It's fundamental for the later use of the mesh within another tool that all the nodes defining triangular faces are part of one of the tetrahedra of the 3D mesh.<br>

<br>Thanks<br><br>Valentina <br><br><br>Thick-Walled Prolate Spheroid Script<br><br>lc = 0.09; //this value changes the resolution of the mesh<br>lz_o = 1; // outer long axis<br>lx_o = 0.5; // outer radius<br>wt=0.15;<br>

lx_i=lx_o-wt; //inner radius<br>lz_i=lz_o-wt; //inner long axis<br><br>Point(1) = {0,0,0,lc};<br>Point(2) = {lx_o,0,0,lc};<br>Point(3) = {0,0,-lz_o,lc};<br>Point(4) = {lx_i,0,0,lc};<br>Point(5) = {0,0,-lz_i,lc};<br>Ellipse(1) = {2,1,3,3};<br>

Ellipse(2) = {4,1,5,5};<br><br>//OUTER ELLIPSOID<br><br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{1};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{3};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{6};<br>

}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{9};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{12};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{15};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{18};<br>

}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{21};<br>}<br><br>//INNER ELLIPSOID<br><br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{2};<br>}<br><br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{27};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>

  Line{30};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{33};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{36};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{39};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>

  Line{42};<br>}<br>Extrude {{0,0,1}, {0,0,0}, Pi/4} {<br>  Line{45};<br>}<br><br>//Upper lines connecting inner and outer half spheroid and the related surfaces<br><br><br>Line(51) = {16, 9};<br>Line(52) = {15, 8};<br>Line(53) = {14, 7};<br>

Line(54) = {13, 6};<br>Line(55) = {4, 2};<br>Line(56) = {19, 12};<br>Line(57) = {18, 11};<br>Line(58) = {17, 10};<br>Line Loop(59) = {40, 58, -16, -51};<br>Plane Surface(60) = {59};<br>Line Loop(61) = {51, -13, -52, 37};<br>

Plane Surface(62) = {61};<br>Line Loop(63) = {52, -10, -53, 34};<br>Plane Surface(64) = {63};<br>Line Loop(65) = {53, -7, -54, 31};<br>Plane Surface(66) = {65};<br>Line Loop(67) = {54, -4, -55, 28};<br>Plane Surface(68) = {67};<br>

Line Loop(69) = {55, -25, -56, 49};<br>Plane Surface(70) = {69};<br>Line Loop(71) = {56, -22, -57, 46};<br>Plane Surface(72) = {71};<br>Line Loop(73) = {57, -19, -58, 43};<br>Plane Surface(74) = {73};<br>Plane Surface(75) = {59};<br>

Surface Loop(76) = {14, 17, 20, 23, 26, 5, 8, 11, 64, 62, 38, 41, 44, 47, 50, 29, 32, 35, 66, 68, 70, 72, 74, 60};<br>Volume(77) = {76};<br><br>