[Nek5000-users] Problem in using .re2 file

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu May 18 09:25:18 CDT 2017


Please attached the .rea file. This is the smaller version (generated
with the .re2 file). Please let me know if you want the full version as 
well.

Thanks for your time.

NN.


On Thursday 18 May 2017 04:19 PM, nek5000-users at lists.mcs.anl.gov wrote:
> Can you please post your .rea file
>   
> -----Original message-----
>> From:nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
>> Sent: Thursday 18th May 2017 16:13
>> To: nek5000-users at lists.mcs.anl.gov
>> Subject: Re: [Nek5000-users] Problem in using .re2 file
>>
>> I have prepared the mesh using GAMBIT. Then using a converter I generated
>> the .rea file in ASCII format. This .rea file works fine with NEK5000.
>> However,
>> when I convert to .re2 the problem shows up...
>>
>>
>> On Thursday 18 May 2017 04:06 PM, nek5000-users at lists.mcs.anl.gov wrote:
>>> How do you generate your mesh/rea?
>>>    
>>> -----Original message-----
>>>> From:nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
>>>> Sent: Thursday 18th May 2017 16:02
>>>> To: nek5000-users at lists.mcs.anl.gov
>>>> Subject: [Nek5000-users] Problem in using .re2 file
>>>>
>>>> Hi Neks,
>>>>    I have encountered a problem while using .re2 file. When I use
>>>>         it,
>>>>         I get the following warning in the output file and the code hangs
>>>>         up.
>>>>       WARNINGb: Detected non-right-handed element.
>>>>               Number   14621  V1-8:  -.1290E+21   .0000E+00   .6991E+01
>>>>             .0000E+00
>>>>               WARNINGb: Detected non-right-handed element.
>>>>               WARNINGb: Detected non-right-handed element.
>>>> If, however, I use the full ASCII  .rea file, the code runs fine.
>>>>         What
>>>>         could be the problem in using .re2 file? Kindly suggest a remedy.
>>>> I would like to mention that first I create the full ASCII .rea
>>>>         file and generate
>>>>         the .map file from it. Then, I use the executable "reatore2" to
>>>>         make the
>>>>         .re2 file and smaller ASCII .rea file. I do not make any change in
>>>>         the smaller
>>>>         ASCII .rea file.
>>>> Many thanks,
>>>> NN.
>>>> _______________________________________________
>>>> Nek5000-users mailing list
>>>> Nek5000-users at lists.mcs.anl.gov
>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

-------------- next part --------------
 ****** PARAMETERS *****
    2.610000     NEKTON VERSION
            3 DIMENSIONAL RUN
          103 PARAMETERS FOLLOW
   1.00000         p1  DENSITY
   -100.00          p2  VISCOS
   0.
   0.
   0.
   0.
   1.00000         p7  RHOCP
   1.00000         p8  CONDUCT
   0.
   0.              p10 FINTIME
   80000.0          p11 NSTEPS
   -0.0015            p12 DT
   0.              p13 IOCOMM
   0.              p14 IOTIME
   500.000         p15 IOSTEP
   0.              p16 PSSOLVER
   0.
   0.250000E-01    p18 GRID
   -1.00000        p19 INTYPE
   7.0000          p20 NORDER
   0.100000E-10    p21 DIVERGENCE
   0.100000E-11    p22 HELMHOLTZ
   0               p23 NPSCAL
   0.100000E-09    p24 TOLREL
   0.100000E-09    p25 TOLABS
   5.00000         p26 COURANT/NTAU
   2.00000         p27 TORDER
   0.              p28 TORDER: mesh velocity (0: p28=p27)
   0.              p29 magnetic visc if > 0, = -1/Rm if < 0
   0.              p30 > 0 ==> properties set in uservp()
   -1.              p31 NPERT: #perturbation modes
   0.              p32 #BCs in re2 file, if > 0
   0.
   0.
   0.
   0.
   0.
   0.
   0.
   0.
   0.              p41 1-->multiplicative SEMG
   0.              p42 0=gmres/1=pcg
   0.00000         p43 0=semg/1=schwarz
   0.              p44 0=E-based/1=A-based prec.
   0.              p45 Relaxation factor for DTFS
   0.              p46 reserved
   0.              p47 vnu: mesh matieral prop
   0.
   0.
   0.
   0.
   0.              p52 IOHIS
   0.
   0.              p54 1,2,3-->fixed flow rate dir=x,y,z
   0.              p55 vol.flow rate (p54>0) or Ubar (p54<0)
   0.
   0.
   0.
   0.              p59 !=0 --> full Jac. eval. for each el.
   0.              p60 !=0 --> init. velocity to small nonzero
   0.
   0.              p62 >0 --> force byte_swap for output
   0.              p63 =8 --> force 8-byte output
   0.              p64 =1 --> perturbation restart
   1.              p65 #iofiles (eg, 0 or 64); <0 --> sep. dirs
   4.000         p66 output : <0=ascii, else binary
   4.000         p67 restart: <0=ascii, else binary
   0.              p68 iastep: freq for avg_all
   0.
   0.
   0.
   0.
   0.
   0.              p74 verbose Helmholtz
   0.
   0.
   0.
   0.
   0.
   0.
   0.
   0.
   0.
   0.               p84 !=0 --> sets initial timestep if p12>0
   0.               p85 dt ratio if p84 !=0, for timesteps>0
   0.               p86 reserved
   0.
   0.
   0.
   0.
   0.
   0.
   20.0000          p93 Number of previous pressure solns saved
   5.00000          p94 start projecting velocity after p94 step
   3.00000          p95 start projecting pressure after p95 step
   0.00000          p96 u0 translational velocity (in .usr)
   0.00000          p97 v0 translational velocity (in .usr)
   0.
   4.00000          p99   dealiasing: <0--> off/3--> old/4--> new
   0.
   0.               p101   No. of additional filter modes
   1.00000          p102   Dump out divergence at each time step
   .000             p103   weight of stabilizing filter (.01)
      4  Lines of passive scalar data follows2 CONDUCT; 2RHOCP
   1.00000       1.00000       1.00000       1.00000       1.00000
   1.00000       1.00000       1.00000       1.00000
   1.00000       1.00000       1.00000       1.00000       1.00000
   1.00000       1.00000       1.00000       1.00000
           13  LOGICAL SWITCHES FOLLOW
  T     IFFLOW
  F     IFHEAT
  T     IFTRAN
  T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
  F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
  F     IFAXIS
  F     IFSTRS
  F     IFSPLIT
  F     IFMGRID
  F     IFMODEL
  F     IFKEPS
  F     IFMVBD
  F     IFCHAR
   8.00000       8.00000     -0.500000      -4.00000     XFAC,YFAC,XZERO,YZERO
 **MESH DATA** 6 lines are X,Y,Z;X,Y,Z. Columns corners 1-4;5-8
      -21336  3         21336     NELT,NDIM,NELV
            0 PRESOLVE/RESTART OPTIONS  *****
            7         INITIAL CONDITIONS *****
C Default
C Default
C Default
C Default
C Default
C Default
C Default
  ***** DRIVE FORCE DATA ***** BODY FORCE, FLOW, Q
            4                 Lines of Drive force data follow
C
C
C
C
  ***** Variable Property Data ***** Overrrides Parameter data.
            1 Lines follow.
            0 PACKETS OF DATA FOLLOW
  ***** HISTORY AND INTEGRAL DATA *****
            0   POINTS.  Hcode, I,J,H,IEL
  ***** OUTPUT FIELD SPECIFICATION *****
            6 SPECIFICATIONS FOLLOW
  T      COORDINATES
  T      VELOCITY
  T      PRESSURE
  F      TEMPERATURE
  F      TEMPERATURE GRADIENT
            0      PASSIVE SCALARS
  ***** OBJECT SPECIFICATION *****
       0 Surface Objects
       0 Volume  Objects
       0 Edge    Objects
       0 Point   Object


More information about the Nek5000-users mailing list