[Nek5000-users] Error in generalev with MHD

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Jul 1 07:51:02 CDT 2014


Hi Jan,

On line 5912, velocity BC is written with capital letter V instead of 
small letter v (as in 5928).

If you write velocity with small letter, it'll read the BC in the .usr 
file (userbc subroutine), otherwise you have to specify the values of 
velocity in the .rea file.

In addition, in the elements (E) BC you have to specify

E     n1    n2   n3   n4   0  0  0
n1: Number of element
n2: side of element
n3: element connected
n4: side of the element connected.

However, you've set n3 and n4 to 0.

I hope it helps you.

Regards,
SL



El 01-07-2014 13:56, nek5000-users at lists.mcs.anl.gov escribió:
> Hello Neks,
> I' trying to solve the MHD equations and encounter the same error 
> reported here:
> http://lists.mcs.anl.gov/pipermail/nek5000-users/2011-June/001365.html.
> For testing, I have set up a small grid (384 elements) with B=0 at the
> boundaries and B=0 as initial conditions, and the error appears before
> the first timestep. As mentioned in the documentation, I have used the
> same BC-types for the B Field as for the velocity, so in the .rea file
> it reads:
> 
> 5909   ***** PASSIVE SCALAR 1  BOUNDARY CONDITIONS *****
> 5910  E    1  1   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5911  E    1  2   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5912  V    1  3   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5913  E    1  4   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5914  E    1  5   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5915  E    1  6   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5916  E    2  1   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5917  E    2  2   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5918  E    2  3   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5919  E    2  4   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5920  E    2  5   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5921  E    2  6   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5922  E    3  1   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5923  E    3  2   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5924  E    3  3   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5925  E    3  4   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5926  E    3  5   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5927  E    3  6   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5928  v    4  1   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5929  E    4  2   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> 5930  E    4  3   0.00000       0.00000       0.00000 0.00000       
> 0.00000
> etc.
> 
> The thread above mentioned that it was due to incorrect BCs, but
> unfortunately no more details were given.
> 
> Cheers,
> Jan
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users




More information about the Nek5000-users mailing list