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