<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Dear George,<br>
    <br>
    You can also use the "Fast Curving" algorithm in the "High order
    tools" after converting your 3D mesh to high order. It is still
    work-in-progress and not guaranteed to work, specially for order
    > 2, but for the case you sent, it corrects the invalid 2nd-order
    3D mesh in less than 2 seconds.<br>
    <br>
    Best regards,<br>
    <br>
    Thomas<br>
    <br>
    <br>
    <div class="moz-cite-prefix">On 12/13/2013 02:14 PM, Jean-François
      Remacle wrote:<br>
    </div>
    <blockquote
      cite="mid:7042DA22-C077-4217-9039-2B36E76FA496@uclouvain.be"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
      <br>
      <div>
        <div>Le 11 déc. 2013 à 11:53, George Ntemos <<a
            moz-do-not-send="true"
            href="mailto:george.ntemos11@imperial.ac.uk">george.ntemos11@imperial.ac.uk</a>>
          a écrit :</div>
        <br class="Apple-interchange-newline">
        <blockquote type="cite">
          <meta http-equiv="content-type" content="text/html;
            charset=ISO-8859-1">
          <div bgcolor="#FFFFFF" text="#000000">
            <meta http-equiv="CONTENT-TYPE" content="text/html;
              charset=ISO-8859-1">
            <p style="margin-bottom: 0cm">Dear Gmsh Development Team,<br>
              <br>
              <br>
              My name is George Ntemos and I am a PhD student in the
              Aeronautics Department of Imperial College London. Our
              team is developing and applying CFD software, using gmsh
              for high-order meshing purposes.</p>
            <p style="margin-bottom: 0cm">I am currently studying a flow
              case around a rod-aerofoil set-up and I've been
              experiencing some problems with high-order meshing. More
              specifically, the meshing algorithm fails to mesh in high
              order around the (spline-defined) aerofoil profile. Please
              find attached a lightweight .geo file of my configuration
              just for demonstration purposes. If one tries meshing this
              file in 3d and for order > 1, gmsh will warn that, to
              my understanding, the algorithm fails to properly converge
              at certain points, visibly resulting in quite severe
              discrepancies for a few elements (curvilinear on). It also
              warns that it cannot orient the normal for the two
              extruded surfaces of the rod (cylinder) boundary layer. It
              is worth noting here that I am experiencing none of the
              above problems when I am meshing in 2d. Since my 3d mesh
              is just a simple extrusion, I would assume that extending
              to 3d is just a trivial case of translating every
              operation along the 3<sup>rd</sup> dimension. I am of
              course not particularly experienced in using gm</p>
          </div>
        </blockquote>
        <div><br>
        </div>
        <div>The 2D mesh is already wrong as is. If you use the high
          order tools with the current parameters, this</div>
        <div>is what it gives : </div>
        <div><br>
        </div>
        <div>
          <div>Info    : --- Optimization pass with initial jac. range
            (-0.528961, 1.76584), jacBar = -0.581857</div>
          <div>Info    : Reached jacobian (0.148494 1.77654)
            requirements, setting null gradient</div>
          <div>Info    : Optimization finalized after 11 iterations (41
            function evaluations),</div>
          <div>Info    : because gradient norm is no more than EpsG</div>
          <div>Info    : Optimization done Range (0.148494,1.77654)</div>
          <div>Info    : Optimization succeeded</div>
          <div>Info    : Done optimizing high order mesh (278.816 s)</div>
          <div><br>
          </div>
        </div>
        <div>It is quite long indeed (work clearly in progress for
          boundary layers) but it worked !</div>
        <div><br>
        </div>
        <div>JF</div>
        <br>
        <blockquote type="cite">
          <div bgcolor="#FFFFFF" text="#000000">
            <p style="margin-bottom: 0cm">sh however and I may very well
              be missing something here.</p>
            <p style="margin-bottom: 0cm">Any help regarding these
              issues would be greatly appreciated.</p>
            <p style="margin-bottom: 0cm">Thank you for your time.</p>
            <p style="margin-bottom: 0cm"><br>
            </p>
            <p style="margin-bottom: 0cm">Best regards,</p>
            <p style="margin-bottom: 0cm">George</p>
            <title></title>
            <meta name="GENERATOR" content="LibreOffice 4.0.4.2 (Linux)">
            <style type="text/css">
        <!--
                @page { margin: 2cm }
                P { margin-bottom: 0.21cm }
        -->
        </style> </div>
          <span><rod_aerofoil_3d_light.geo></span>_______________________________________________<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 apple-content-edited="true">
        <span class="Apple-style-span" style="border-collapse: separate;
          border-spacing: 0px; ">
          <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-size:
              12px; font-style: normal; font-variant: normal;
              font-weight: normal; letter-spacing: normal; line-height:
              normal; orphans: 2; 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; ">
              <div style="word-wrap: break-word; -webkit-nbsp-mode:
                space; -webkit-line-break: after-white-space; ">
                <div>
                  <div><span class="Apple-style-span" style="font-size:
                      medium; ">
                      <div>------------------------------------------------------------------</div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Prof. Jean-Francois Remacle</span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Universite catholique de Louvain
                            (UCL)</span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Ecole Polytechnique de Louvain (EPL)
                            - Louvain School of Engineering</span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Institute of Mechanics, Materials
                            and Civil Engineering (iMMC)</span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Center for Systems Engineering and
                            Applied Mechanics (CESAME)</span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; ">Tel : +32-10-472352 -- Mobile :
                            +32-473-909930 </span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; "><br>
                          </span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; "> </span></font></div>
                      <div><font class="Apple-style-span" size="3"><span
                            class="Apple-style-span" style="font-size:
                            12px; "><br>
                          </span></font></div>
                      <div><br>
                      </div>
                    </span></div>
                </div>
              </div>
            </span></div>
          <br class="Apple-interchange-newline">
        </span><br class="Apple-interchange-newline">
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
    </blockquote>
    <br>
    <pre class="moz-signature" cols="72">-- 
Thomas Toulorge
Université catholique de Louvain (UCL)
Institute of Mechanics, Materials and Civil Engineering (iMMC)
Applied Mechanics and Mathematics Division (MEMA)
Avenue George Lemaître 4-6 bte L4.05.02
1348 Louvain-la-Neuve, BELGIUM
Tel.: +32 10 47 23 54
</pre>
  </body>
</html>