[Nek5000-users] Error in generalev with MHD

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


Hey,
thanks for the suggestions.
I want to change some of the boundarie values to values =!0 for the 
actual simulations, that's why there is a mix of small and capital V's. 
However, userbc sets the bc to zero. I will try to use only one of the 
two, in case there is a problem with mixing v and V.
As far as I know, the specification for the Element-Element-Boundary (E) 
you mentioned is outdated and handled by the solver. I have run several 
hydrodynamic simulations with 0 for the options and never encountered 
any problems.

Am 01.07.14 14:51, schrieb nek5000-users at lists.mcs.anl.gov:
> 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
>
>
> _______________________________________________
> 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