[Nek5000-users] Rayleigh Benard with symmetry BCs

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Feb 18 11:53:35 CST 2013


My apologies upfront for any simple misconceptions as I'm new to nek.

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).

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:

     1   100 **ERROR**: Failed in HMHOLTZ: TEMP   1.4999E-02   3.7064E+01   2.4836E-03
          1   2.0000E-02  6.0370E-03 Heat done

     1   100 **ERROR**: Failed in HMHOLTZ: VELX   2.8743E-21   1.7086E-16   2.3416E-23
  0.000000000000000E+000  p22            1           1
 New CG1-tolerance (RINIT*epsm) =   1.010379228335448E-010
  2.341589846098494E-023


     1   100 **ERROR**: Failed in HMHOLTZ: VELY   1.0837E-01   1.0104E+03   1.0104E-10
 New CG2-tolerance (RINIT*10-5/10-10) =   2.795810823528297E-003
  3.472834773690102E-025



Do I need to adjust the tolerances because the mesh/BCs are different?


I don't think it matters, but I ran this on 16 procs (there are 16 per node on the machine I'm using).

                                               -Jared

CCS2/CNLS Postdoctoral Research Associate
Los Alamos National Laboratory
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130218/8550233e/attachment-0001.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: NEK.map
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130218/8550233e/attachment-0001.ksh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NEK.rea
Type: application/octet-stream
Size: 30722 bytes
Desc: NEK.rea
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130218/8550233e/attachment-0001.obj>


More information about the Nek5000-users mailing list