Hi Alex,<br><br>Lowering your DT (parameter 12) slightly will keep a more stable courant number and NEK won't blow up.<br><br><br>-Katie<br><br><div class="gmail_quote">On Thu, Aug 4, 2011 at 10:18 AM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Alex,<br>
<br>
I see that p26 is set to 3 -- this is the target courant number -- which should be 2, no?  You might even make it smaller than that to start with.<br>
<br>
Also, for my nondimensionalization, I had<br>
<br>
  1.000000E+00 p054: direction of flow rate -- 1-X, 2-Y, 3-Z<br>
  12.56637E+00 p055: volumetric flow rate<br>
<br>
whereas you have<br>
<br>
   -1.00000     p054 fixed flow rate dir: |p54|=1,2,3=x,y,z<br>
   1.00000     p055 vol.flow rate (p54>0) or Ubar (p54<0)<br>
<br>
finally,  I have<br>
<br>
  4.000000E+00     p099 dealiasing<br>
  0.050            p103 filter weight of last mode<br>
<br>
whereas you have<br>
<br>
   3.00000     p099 dealiasing: <0--> off/3--> old/4--> new<br>
  0.100000E-01 p103 weight of stabilizing filter (.01)<br>
<br>
<br>
Let me know if making those changes helps.<br>
<br>
--Mike<br>
<div><div></div><div class="h5"><br>
<br>
<br>
<br>
On Aug 3, 2011, at 3:41 AM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><mailto:<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>>> wrote:<br>


<br>
Hi neks,<br>
<br>
I'm trying to simulate a turbulent channel flow like in the presentation "A comparison<br>
of Nek5000 and OpenFOAM for DNS of turbulent channel flow" of Michael A. Sprague from the last<br>
user meeting (<a href="http://www.mcs.anl.gov/%7Efischer/nek5000/sprague_nek5000_dec2010.pdf" target="_blank">http://www.mcs.anl.gov/~fischer/nek5000/sprague_nek5000_dec2010.pdf</a> "mid"-case).<br>
In the presentation it is mentioned that the CFL number is approx. 2 so I have choosen the BDF3/OIFS3<br>
scheme. However the solution blows up and Nek5000 stops the computation after 49 steps.<br>
<br>
Are there any parameters I have to set in the rea-file to get it work? TCHAR is true and p026 is<br>
set to 3. Dealiasing is turned on as well, but I'm not sure if the weight of the filter may be wrong.<br>
<br>
The files of the case are attached.<br>
<br>
Thanks!<br>
Alex<br>
</div></div><sprague3D.tar.gz><ATT00001..txt><br>
<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>
</blockquote></div><br>