<div dir="ltr"><div><div><div>Hi NEKs,<br><br></div>I just started working on NEK5000 and am having trouble with my first real simulation (tried both PN/PN and PN/PN-2). I first created a .msh file of simple cylindrical geometry, with velocity-inlet (v), Outlet (O) and surrounding walls (W), in ICEM and used "mshconvert" to get a .rea file. In NEK, the simulation is dying before they actually begin. The error message is as follows:<br><br>*************************************************************************<br>*************************************************************************<br>Matrix:   0 aa         4       4       4<br>  0 aa        Infinity  -4.08542E+02   7.41342E+01   0.00000E+00<br>  0 aa    -4.08542E+02   4.65569E+02  -1.31161E+02   0.00000E+00<br>.<br>.<br>Matrix:   0 Aeig       4       4       4<br>  0 Aeig           NaN           NaN           NaN           NaN<br>  0 Aeig           NaN           NaN           NaN           NaN<br>  0 Aeig           NaN           NaN           NaN           NaN<br>  0 Aeig           NaN           NaN           NaN           NaN<br>.<br>.<br>Matrix:   0 Deig       1       4       4<br>  0 Deig           NaN           NaN           NaN           NaN<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br> Error in generalev, info=           3           4           1<br>  <br>call exitt: dying ...<br>*************************************************************************<br>*************************************************************************<br><br></div><div>I saw an older post regarding a similar problem and that was due to improper boundary conditions. I suspect I may be doing something similar but couldn't find anything until now. The .msh file is fine as I was able to import to "fluent" and run the simulation. Please let me know if you can help in some way. I am posting input files below (please pardon the length of this message):<br><br></div><div>************** FROM .rea FILE*********************<br></div><div>           0  objects found<br> 103   Parameters from file:<br>   1      1.00000         p1  DENSITY<br>   2     -500.000         p2  VISCOS<br>   7      1.00000         p7  RHOCP<br>   8      1.00000         p8  CONDUCT<br>  11      10000           p11 NSTEPS<br>  12      0.10000E-05     p12 DT<br>  15      5               p15 IOSTEP<br>  18     -20.00000        p18 GRID<br>  19     -1.00000         p19 INTYPE<br>  20      10.0000         p20 NORDER<br>  26      1.00000         p26 COURANT/NTAU<br>  27      2.00000         p27 TORDER<br>  65      1.00000         p65 #iofiles (eg, 0 or 64); <0 --> sep. dirs<br>  66      4.00000         p66 output : <0=ascii, else binary<br>  67      4.00000         p67 restart: <0=ascii, else binary<br>  93      20.0000          p93 Number of previous pressure solns saved<br>  94      3.00000          p94 start projecting velocity after p94 step<br>  95      5.00000          p95 start projecting pressure after p95 step<br>  99      0.               p99   dealiasing: <0--> off/3--> old/4--> new<br> 102      1.00000          p102   Dump out divergence at each time step<br> 103     -1.00000          p103   weight of stabilizing filter (.01)<br><br>           13  LOGICAL SWITCHES FOLLOW<br>  T     IFFLOW<br>  T     IFHEAT<br>  T     IFTRAN<br>  T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)<br>  F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)<br>  F     IFAXIS<br>  F     IFSTRS<br>  F     IFSPLIT<br>  F     IFMGRID<br>  F     IFMODEL<br>  F     IFKEPS<br>  F     IFMVBD<br>  F     IFCHAR<br>.<br>.<br> **MESH DATA**<br>       22620       3       22620      NEL,NDIM,NELV<br>      ELEMENT          1 [    1 ]    GROUP     0<br>  0.00000000e+00  0.00000000e+00  0.00000000e+00  0.00000000e+00<br> -7.07106729e-01 -6.82106736e-01 -5.66056759e-01 -5.87769989e-01<br>  7.07106833e-01  6.82106832e-01  7.76064003e-01  8.08965813e-01<br>  1.02564103e+00  1.02564103e+00  1.02564103e+00  1.02564103e+00<br> -7.07106715e-01 -6.82106722e-01 -5.65372395e-01 -5.87023410e-01<br>  7.07106818e-01  6.82106817e-01  7.75121255e-01  8.07937362e-01<br>      ELEMENT          2 [    1 ]    GROUP     0<br>  1.02564103e+00  1.02564103e+00  1.02564103e+00  1.02564103e+00<br> -6.82106722e-01 -5.65372395e-01 -5.43721379e-01 -6.57106729e-01<br>  6.82106817e-01  7.75121255e-01  7.42305149e-01  6.57106816e-01<br>  0.00000000e+00  0.00000000e+00  0.00000000e+00  0.00000000e+00<br> -6.82106736e-01 -5.66056759e-01 -5.44343529e-01 -6.57106742e-01<br>  6.82106832e-01  7.76064003e-01  7.43162193e-01  6.57106831e-01<br>.<br>.<br>  ***** FLUID   BOUNDARY CONDITIONS ***** <br> E      1  1 1.4149000e+04 1.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br> E      1  2 2.0000000e+00 1.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br> E      1  3 1.3000000e+01 1.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br> W      1  4 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br> v      1  5 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br>.<br> O  22620  6 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00<br>**************************************************************<br><br></div><div>************** .usr FILE**************<br></div><div>c All my user functions are empty except userbc to deal with the "v" condition<br><br>      subroutine userbc (ix,iy,iz,iside,ieg)<br>      include 'SIZE'<br>      include 'TSTEP'<br>      include 'INPUT'<br>      include 'NEKUSE'<br><br>      ux=0.005<br>      uy=0.0<br>      uz=0.0<br><br>      return<br>      end<br>****************************************<br></div><div><br></div><div>***********************FIRST FEW LINES OF .map FILE***********<br>       22620       24040          14       16384      180960       24040           0<br>        3597          14        4574        4561        4560        4562        4575        5164        5162<br>        3596        4575        5162        4581        5160        4574        4560        4580        4579<br>        3601        4581        5160        4572        5158        4580        4579        4571        4566<br>        3603        4572        5158        4570        5182        4571        4566        4569        5181<br>        3607        4570        5182        4693        5178        4569        5181        4692        5177<br>        3681        4693        5178        4645        5176        4692        5177        4644        5175<br>        3683        4645        5176        4649        5174        4644        5175        4648        5173<br>        3691        4649        5174        5277        5276        4648        5173        5275        5274<br>        3107        5277        5276        4035        4020        5275        5274        4034        4019<br>        3106        4035        4020        4033        4026        4034        4019        4032        4024<br>        3119        4033        4026        4025        4277        4032        4024        4023        4276<br>        3119        4025        4277        4055        4273        4023        4276        4022        4272<br>        3599        5164        4561        5172        5171        5162        4560        5170        5169<br>******************************************************************************<br><br></div><div></div>Thanks,<br></div>Ajay <br></div>