[Nek5000-users] Error in generalev with MHD
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Wed Jul 2 07:49:14 CDT 2014
Hello again,
I changed all BC's for the magnetic field to 'v' and used userbc to set
the value to zero, but the error popped up again. The output reads:
216 Matrix: 0 aa 6 6 6
217 0 aa 7.80579E+01 -6.39449E+01 1.81958E+01 -7.90484E+00
2.16525E+00 0.00000E+00
218 0 aa -6.39449E+01 8.67803E+01 -4.74885E+01 1.84117E+01
-1.07792E+01 5.66872E+00
219 0 aa 1.81958E+01 -4.74885E+01 5.96643E+01 -4.03684E+01
3.02051E+01 -2.06952E+01
220 0 aa -7.90484E+00 1.84117E+01 -4.03684E+01 6.64609E+01
-7.64969E+01 4.76373E+01
221 0 aa 2.16525E+00 -1.07792E+01 3.02051E+01 -7.64969E+01
1.90127E+02 -1.67410E+02
222 0 aa 0.00000E+00 5.66872E+00 -2.06952E+01 4.76373E+01
-1.67410E+02 Infinity
223
224 Matrix: 0 bb 6 6 6
225 0 bb 7.19606E-02 9.84283E-03 -2.07087E-03 -1.02002E-03
7.34397E-04 0.00000E+00
226 0 bb 9.84283E-03 4.07143E-02 -4.22288E-03 7.31377E-03
-3.21920E-03 0.00000E+00
227 0 bb -2.07087E-03 -4.22288E-03 1.05630E-01 -1.61560E-02
6.44826E-03 0.00000E+00
228 0 bb -1.02002E-03 7.31377E-03 -1.61560E-02 1.05960E-01
-6.21764E-03 0.00000E+00
229 0 bb 7.34397E-04 -3.21920E-03 6.44826E-03 -6.21764E-03
5.05023E-02 0.00000E+00
230 0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 NaN
231
232 Matrix: 0 Aeig 6 6 6
233 0 Aeig NaN NaN NaN NaN NaN
NaN
234 0 Aeig NaN NaN NaN NaN NaN
NaN
235 0 Aeig NaN NaN NaN NaN NaN
NaN
236 0 Aeig NaN NaN NaN NaN NaN
NaN
237 0 Aeig NaN NaN NaN NaN NaN
NaN
238 0 Aeig NaN NaN NaN NaN NaN
NaN
239
240 Matrix: 0 Deig 1 6 6
241 0 Deig NaN NaN NaN NaN NaN
NaN
242 Error in generalev, info= 5 6 1
243
244 call exitt: dying ...
So the problem is probably the NaN and Inf entries. As a next step, I
took the gpf example as a starting point, copy pasted my own grid and bc
descriptions, set the bc's to zero in the usr-file and removed the
Temperaturefield. This produced the same error. Can this be linked to
the absence of the temperature field? If I leave out the B-Field the
case is running just fine, altough the BC's are the same in the *.rea file.
Jan
More information about the Nek5000-users
mailing list