<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Hi Ruth,<br>
    <br>
    I'm able to get a volume mesh for some objects, but I'm finding that
    for others I'm getting:<br>
    <br>
    Error   : Self intersecting surface mesh, computing intersections
    (this could take a while)<br>
    Error   : Surface mesh is wrong, cannot do the 3D mesh<br>
    <br>
    and<br>
    <br>
    Error:  Invalid PLC.<br>
      Facet (41, 172, 188) (3749) is not a valid polygon.<br>
    <br>
    An example is attached - command line is<br>
    <br>
    gmsh try3.geo -3 -format msh -saveall -smooth 3 -o try3-out.gmsh<br>
    <br>
    If I change the geometry somewhat it works fine, which makes me
    think there is some problem with the edges of the surfaces.  In
    general I can't force the triangulation on the geometry surfaces to
    match on curved edges where they meet, and I can see there are small
    gaps between the triangulations (small compared to the mesh I want).<br>
    <br>
    Is this causing the problem?  Would it help if I specified what
    edges are adjacent to what surfaces, or provided some kind of
    tolerance? <br>
    <br>
    Thanks,<br>
    <br>
        Wayne<br>
    <br>
    On 3/25/12 12:05 AM, Ruth V. Sabariego wrote:
    <blockquote
      cite="mid:73AC802C-48CC-4F74-BD24-25D1A5CFE610@ulg.ac.be"
      type="cite">Hi Wayne,
      <div>For that you have to reclasify the elements of your geometry
        and create the volume with the new surfaces.</div>
      <div><br>
      </div>
      <div>Try this:</div>
      <div>**********************</div>
      <div>
        <div>Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2)
          split only with metis</div>
        <div><br>
        </div>
        <div>Merge "try1-geom.gmsh";</div>
        <div>RefineMesh;</div>
        <div>CreateTopology;</div>
        <div><br>
        </div>
        <div>ll[] = Line "*";</div>
        <div>For j In {0 : #ll[]-1}</div>
        <div>  Compound Line(newl) = ll[j];</div>
        <div>EndFor</div>
        <div><br>
        </div>
        <div>ss[] = Surface "*";</div>
        <div>For j In {0 : #ss[]-1}</div>
        <div>  ssn[] += news ;</div>
        <div>  Compound Surface(news) = ss[j];</div>
        <div>EndFor</div>
        <div><br>
        </div>
        <div><br>
        </div>
        <div>Surface Loop(10) = { ssn[] };</div>
        <div>Volume(11) = {10};</div>
        <div><br>
        </div>
        <div>Physical Volume("Cylinder.2") = {11};</div>
      </div>
      <div>**********************</div>
      <div><br>
      </div>
      <div><br>
      </div>
      <div>Regards,</div>
      <div>Ruth</div>
      <div><br>
        <div apple-content-edited="true">
          <div style="word-wrap: break-word; -webkit-nbsp-mode: space;
            -webkit-line-break: after-white-space; "><span
              class="Apple-style-span" style="border-collapse: separate;
              color: rgb(0, 0, 0); font-family: Helvetica; font-style:
              normal; font-variant: normal; font-weight: normal;
              letter-spacing: normal; line-height: normal; orphans: 2;
              text-align: -webkit-auto; text-indent: 0px;
              text-transform: none; white-space: normal; widows: 2;
              word-spacing: 0px; -webkit-border-horizontal-spacing: 0px;
              -webkit-border-vertical-spacing: 0px;
              -webkit-text-decorations-in-effect: none;
              -webkit-text-size-adjust: auto; -webkit-text-stroke-width:
              0px; font-size: medium; ">
              <div style="word-wrap: break-word; -webkit-nbsp-mode:
                space; -webkit-line-break: after-white-space; ">
                <div>--</div>
                <div>Dr. Ir. Ruth V. Sabariego</div>
                <div>University of Liege, Electrical Engineering &
                  Computer Science, </div>
                <div>Applied & Computational Electromagnetics (ACE),</div>
                <div>phone: +32-4-3663737 - fax: +32-4-3662910 - <a
                    moz-do-not-send="true"
                    href="http://ace.montefiore.ulg.ac.be/">http://ace.montefiore.ulg.ac.be/</a></div>
                <br class="Apple-interchange-newline">
                <br>
              </div>
            </span></div>
        </div>
        <br>
        <div>
          <div>On 24 Mar 2012, at 23:41, Wayne Christopher wrote:</div>
          <br class="Apple-interchange-newline">
          <blockquote type="cite">
            <meta content="text/html; charset=ISO-8859-1"
              http-equiv="Content-Type">
            <div bgcolor="#FFFFFF" text="#000000"> Hi Ruth,<br>
              <br>
              Thanks - that fixes my volume problem.<br>
              <br>
              My other problem remains, though.  I want to define
              geometry using a msh file, and then create a mesh with new
              surface elements, not the ones I am using for geometry
              surfaces.    Is there a way to do this?<br>
              <br>
                  Wayne<br>
              <br>
              On 3/23/12 2:09 PM, Ruth V. Sabariego wrote:
              <blockquote
                cite="mid:D56BD075-2972-4FE0-818A-0AD5C97536B1@ulg.ac.be"
                type="cite">Hi Wayne, 
                <div><br>
                </div>
                <div>You need to define the geometrical Volume:</div>
                <div>Volume(11) = {10};</div>
                <div>the Physical Volume is used just for tagging the
                  elements and their further use in a e.g. FE solver.</div>
                <div><br>
                </div>
                <div>Regards,</div>
                <div>Ruth</div>
                <div><br>
                  <div apple-content-edited="true"> <span
                      class="Apple-style-span" style="border-collapse:
                      separate; font-family: Helvetica; font-style:
                      normal; font-variant: normal; font-weight: normal;
                      letter-spacing: normal; line-height: normal;
                      orphans: 2; text-align: -webkit-auto; text-indent:
                      0px; text-transform: none; white-space: normal;
                      widows: 2; word-spacing: 0px;
                      -webkit-border-horizontal-spacing: 0px;
                      -webkit-border-vertical-spacing: 0px;
                      -webkit-text-decorations-in-effect: none;
                      -webkit-text-size-adjust: auto;
                      -webkit-text-stroke-width: 0px; font-size: medium;
                      "><span class="Apple-style-span"
                        style="border-collapse: separate; font-family:
                        Helvetica; font-style: normal; font-variant:
                        normal; font-weight: normal; letter-spacing:
                        normal; line-height: normal; orphans: 2;
                        text-align: -webkit-auto; text-indent: 0px;
                        text-transform: none; white-space: normal;
                        widows: 2; word-spacing: 0px;
                        -webkit-border-horizontal-spacing: 0px;
                        -webkit-border-vertical-spacing: 0px;
                        -webkit-text-decorations-in-effect: none;
                        -webkit-text-size-adjust: auto;
                        -webkit-text-stroke-width: 0px; font-size:
                        medium; ">
                        <div style="word-wrap: break-word;
                          -webkit-nbsp-mode: space; -webkit-line-break:
                          after-white-space; "><span
                            class="Apple-style-span"
                            style="border-collapse: separate;
                            font-family: Helvetica; font-style: normal;
                            font-variant: normal; font-weight: normal;
                            letter-spacing: normal; line-height: normal;
                            orphans: 2; text-align: -webkit-auto;
                            text-indent: 0px; text-transform: none;
                            white-space: normal; widows: 2;
                            word-spacing: 0px;
                            -webkit-border-horizontal-spacing: 0px;
                            -webkit-border-vertical-spacing: 0px;
                            -webkit-text-decorations-in-effect: none;
                            -webkit-text-size-adjust: auto;
                            -webkit-text-stroke-width: 0px; font-size:
                            medium; ">
                            <div style="word-wrap: break-word;
                              -webkit-nbsp-mode: space;
                              -webkit-line-break: after-white-space; ">
                              <div>--</div>
                              <div>Dr. Ir. Ruth V. Sabariego</div>
                              <div>University of Liege, Electrical
                                Engineering & Computer Science, </div>
                              <div>Applied & Computational
                                Electromagnetics (ACE),</div>
                              <div>phone: +32-4-3663737 - fax:
                                +32-4-3662910 - <a
                                  moz-do-not-send="true"
                                  href="http://ace.montefiore.ulg.ac.be/">http://ace.montefiore.ulg.ac.be/</a></div>
                              <br class="Apple-interchange-newline">
                              <br>
                            </div>
                          </span></div>
                      </span></span> </div>
                  <br>
                  <div>
                    <div>On 23 Mar 2012, at 18:41, Wayne Christopher
                      wrote:</div>
                    <br class="Apple-interchange-newline">
                    <blockquote type="cite">
                      <div>I'm trying to mesh a volume defined by
                        surfaces that are defined by tessellations.  I
                        don't want the mesh to be based on the surface
                        triangles - those are purely geometry and in
                        most cases finer than the volume mesh that I
                        want.<br>
                        <br>
                        When I run gmsh like this:<br>
                        <br>
                        gmsh try1.geo -3 -format msh -saveall -smooth 3
                        -o try1-out.gmsh<br>
                        <br>
                        (files attached) I do not get any tetras in the
                        output, and it appears that the output mesh is a
                        refinement of the input mesh.<br>
                        <br>
                        What am I doing wrong?  Thanks,<br>
                        <br>
                           Wayne<br>
                        <br>
                        <span><try1.geo></span><span><try1-geom.gmsh></span><span><try1-out.gmsh></span>_______________________________________________<br>
                        gmsh mailing list<br>
                        <a moz-do-not-send="true"
                          href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>
                        <a moz-do-not-send="true"
                          class="moz-txt-link-freetext"
                          href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
                      </div>
                    </blockquote>
                  </div>
                  <br>
                </div>
              </blockquote>
            </div>
            _______________________________________________<br>
            gmsh mailing list<br>
            <a moz-do-not-send="true" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>
            <a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
          </blockquote>
        </div>
        <br>
      </div>
    </blockquote>
  </body>
</html>