<div dir="ltr"><div><div>I check the first time step, the pressure is more or less correct. From the second step, the singular point starts to show up.  <br><br>I changed the Cartesian mesh of the example back, I have the same problem. So the problem looks like to produced during the generation of the mesh around the cylinder. <br>
<br>I used the import.rea from the ext_cyl example, so that I won't introduced any other problem. <br></div>After several tests, I isolated the problem to two potential issues:<br>1, genbox to generate the mesh around cylinder, the script is below<br>
2, prex, after importing the geometry and boundaries conditions, there is only the boundary condition of the cylinder left, which I imposed as 'wall' by using the script below<br><br>I feel that it is very probably caused by prex, but I do not see where is the bug? I attached all scripts in the mail as well. <br>
<br></div>Thank you! <br><br>Tubois<br><div><div><br>pretex << EOF<br>new<br>   1 READ PREVIOUS PARAMETERS <br>box<br>   1 BUILD FROM FILE          <br>box<br>  10 IMPORT MESH              <br>toto<br>y Would you like to displace existing elements in box?<br>
   1 END    ELEMENTS          <br>   1 ACCEPT MATL,QVOL         <br>   1 SET BCs                  <br>  11 SET ENTIRE LEVEL         <br>   3 WALL            <br>n Is this a lattice for hex-close-packed spheres?:<br>1 END  LEVEL               <br>
1 ACCEPT B.C.'s            <br>1 EXIT                     <br>EOF<br>genmap << EOF<br>new<br>.1<br>EOF<br><br><br><br><br>import.rea<br>#========================================================<br>2             spatial dimension<br>
1             number of fields<br>#<br># comments<br>#<br>#<br>#========================================================<br>#<br>Y             cylinder<br>4 -40 1         nelr,nel_theta,nelz<br>.0 .0             x0,y0 - center of cylinder<br>
cccbb             descriptors: c-cyl, o-oct, b-box (1 character + space)<br>.5 .6 .7 .8 1.          r0 r1 ... r_nelr<br>0 1 1             theta0/2pi theta1/2pi ratio<br>E  ,E  ,E  ,E  ,             bc’s (3 characters + comma)<br>
<br><br><br><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 5, 2014 at 7:29 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi nek users, <br><br></div>I am encountering a strange 
issue. Build on the example of cylinder flow, I write a matlab script to
 generate the Cartesian mesh and a genbox script to generate the mesh 
around the cylinder. When the domain is small, the code runs. When I 
increases the domain size and cell numbers, the calculation blows up 
after several steps?  <br>
<br>See the picture attached here, there is a singular point on the 
boundary for the pressure. But the velociy is ok, since I imposed a 
uniform flow for inlet and top bottom. <br><br></div>It also sometimes shows another error, the problem size is too large, even I run them in 32 procs instead of 4 as usual. <br>
</div><div><br></div><div>Here attached the output of nek. <br></div><div><br></div><div>Does someone have an idea? <br><br></div>Best regards<div><div><img src="https://mail.google.com/mail/u/0/images/cleardot.gif"><br>

</div><div>Tubois<br></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Université Pierre et Marie Curie<br>CNRS-LIMSI<br>Campus d'Orsay<br>91403 ORSAY
</div>