Hi Paul,<div><br></div><div>I tried calling fix_geom in usrdat2 and it did not fix the problem. </div><div><br></div><div>Jérémie<br><br>On Sunday, October 26, 2014,  <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div lang="EN-CA" link="#0563C1" vlink="#954F72">
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<div><br>
</div>
Hi Jeremie,
<div><br>
</div>
<div>Try calling "fix_geom"</div>
<div><br>
</div>
<div>in usrdat2 after you modify the sides... this might provide the requisite mesh clean up.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<hr>
<div style="direction:ltr"><font face="Tahoma" color="#000000"><b>From:</b> <a href="javascript:_e(%7B%7D,'cvml','nek5000-users-bounces@lists.mcs.anl.gov');" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="javascript:_e(%7B%7D,'cvml','nek5000-users-bounces@lists.mcs.anl.gov');" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of <a href="javascript:_e(%7B%7D,'cvml','nek5000-users@lists.mcs.anl.gov');" target="_blank">nek5000-users@lists.mcs.anl.gov</a> [<a href="javascript:_e(%7B%7D,'cvml','nek5000-users@lists.mcs.anl.gov');" target="_blank">nek5000-users@lists.mcs.anl.gov</a>]<br>
<b>Sent:</b> Sunday, October 26, 2014 4:10 PM<br>
<b>To:</b> <a href="javascript:_e(%7B%7D,'cvml','nek5000-users@lists.mcs.anl.gov');" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
<b>Subject:</b> [Nek5000-users] Curved GLL Points on Airfoil<br>
</font><br>
</div>
<div></div>
<div>
<div>
<p class="MsoNormal">Hello Neks,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am trying to simulate the flow around an interpolated airfoil with curved boundaries. To do so I am generating the collocation points separately and writing them into an input file, which is read in the usrdat2 subroutine before geometry
 re-initialization. The only curved elements in the mesh are located at the airfoil boundary. I am not getting any errors during geometry re-initialization. However, the calculation starts off with a large courant number (around 7), which ends up diverging
 after 5 times steps. Along the element edges near the airfoil, there seems to be discontinuities in the flow which cause the solution to diverge. However, when inputting non-curvilinear collocation points, written in the exact same order to the input file,
 the simulation runs flawlessly. </p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Are there any switches that must be activated in the rea file for curved elements?</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">If not, are there any general guidelines to using curved elements in this manner?</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Thanks,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Jeremie</p>
</div>
</div>
</div>
</div>
</div>
</div>

</blockquote></div>