<div dir="ltr">Hi Jared,<br><div><div class="gmail_extra"><br></div><div class="gmail_extra">If you send me your .box file, I will check that just in case.<br><br></div><div class="gmail_extra">Also, I see you are using variable time steps, so try reducing the target Courant number to 1/2.<br>

</div><div class="gmail_extra"><br>Another thing to be aware of with variable time stepping: if the Rayleigh number is large and the simulation starts 
from a static state, the exponentially growing velocities might still cause 
the CFL condition to be violated. The solution to this is to decrease 
the maximum allowed time step, DT. This is unlikely to be your problem at 
the current Rayleigh number, however.<br><br></div><div class="gmail_extra">-David<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 18, 2013 at 11:53 AM,  <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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">




<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">My apologies upfront for any simple misconceptions as I'm new to nek.<br>
<br>
I have successfully run the included 2d RB example and thought things 
looked reasonable.  I then changed the mesh to an aspect ratio 2, 4x8 
element mesh and modified the boundary conditions to be 'SYM' sidewalls 
and 'SYM' (free-slip) for the top and bottom
 velocity (see attached rea and map file).<br>
<br>
Without changing anything else from the example (I did not modify SIZE) 
the velocities appear to blow up, and I get a series of errors like the 
following in the logfile:<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: TEMP   1.4999E-02   3.7064E+01   2.4836E-03<br>
          1   2.0000E-02  6.0370E-03 Heat done<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: VELX   2.8743E-21   1.7086E-16   2.3416E-23<br>
  0.000000000000000E+000  p22            1           1<br>
 New CG1-tolerance (RINIT*epsm) =   1.010379228335448E-010<br>
  2.341589846098494E-023<br>
<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: VELY   1.0837E-01   1.0104E+03   1.0104E-10<br>
 New CG2-tolerance (RINIT*10-5/10-10) =   2.795810823528297E-003<br>
  3.472834773690102E-025<br>
<br>
<br>
<br>
Do I need to adjust the tolerances because the mesh/BCs are different?<br>
<br>
<br>
I don't think it matters, but I ran this on 16 procs (there are 16 per node on the machine I'm using).<br>
<br>
                                               -Jared<br>
<div><br>
<div style="font-family:Tahoma;font-size:13px">CCS2/CNLS Postdoctoral Research Associate
<div>Los Alamos National Laboratory</div>
</div>
</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></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Feb 18, 2013 at 11:53 AM,  <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>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">My apologies upfront for any simple misconceptions as I'm new to nek.<br>
<br>
I have successfully run the included 2d RB example and thought things looked reasonable.  I then changed the mesh to an aspect ratio 2, 4x8 element mesh and modified the boundary conditions to be 'SYM' sidewalls and 'SYM' (free-slip) for the top and bottom
 velocity (see attached rea and map file).<br>
<br>
Without changing anything else from the example (I did not modify SIZE) the velocities appear to blow up, and I get a series of errors like the following in the logfile:<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: TEMP   1.4999E-02   3.7064E+01   2.4836E-03<br>
          1   2.0000E-02  6.0370E-03 Heat done<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: VELX   2.8743E-21   1.7086E-16   2.3416E-23<br>
  0.000000000000000E+000  p22            1           1<br>
 New CG1-tolerance (RINIT*epsm) =   1.010379228335448E-010<br>
  2.341589846098494E-023<br>
<br>
<br>
     1   100 **ERROR**: Failed in HMHOLTZ: VELY   1.0837E-01   1.0104E+03   1.0104E-10<br>
 New CG2-tolerance (RINIT*10-5/10-10) =   2.795810823528297E-003<br>
  3.472834773690102E-025<br>
<br>
<br>
<br>
Do I need to adjust the tolerances because the mesh/BCs are different?<br>
<br>
<br>
I don't think it matters, but I ran this on 16 procs (there are 16 per node on the machine I'm using).<br>
<br>
                                               -Jared<br>
<div><br>
<div style="font-family:Tahoma;font-size:13px">CCS2/CNLS Postdoctoral Research Associate
<div>Los Alamos National Laboratory</div>
</div>
</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></div>