[Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file

Thu May 31 08:17:06 CDT 2018

Thanks . I tried W on both sides but stil i get the same error. 

I also tried to convert re2 to rea but the converter did not work well. 

Could you please tell me what is this error? 

Thank you 
Sijo 

De: "nek5000-users" <nek5000-users at lists.mcs.anl.gov> 
À: "nek5000-users" <nek5000-users at lists.mcs.anl.gov> 
Envoyé: Jeudi 31 Mai 2018 15:12:55 
Objet: Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file 

I do not think this method supports periodic boundary conditions. You will have to replace the periodic boundary conditions with wall boundary conditions and see it is runs. 

Marco 

On May 31, 2018, at 9:11 AM, [ mailto:nek5000-users at lists.mcs.anl.gov | 
nek5000-users at lists.mcs.anl.gov ] wrote: 

Hi Marco, 

I tried and it compiled . But before starting the time step it failed. 
Its a very simple case, the water flowing throgh a cube. So I have inlet, outlet, walls on upside and bottom and periodic conditions on side walls. 
I get this error: 

Matrix: 0 aa 6 6 6 
0 aa 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0 aa -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00 
0 aa 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00 
0 aa 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00 
0 aa 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity 
0 aa 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity 

Matrix: 0 bb 6 6 6 
0 bb 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0 bb 0.00000E+00 1.66667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0 bb 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00 
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00 
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 8.33333E-03 0.00000E+00 
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 

Matrix: 0 Aeig 6 6 6 
0 Aeig 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0 Aeig -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00 
0 Aeig 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00 
0 Aeig 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00 
0 Aeig 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity 
0 Aeig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity 

Matrix: 0 Deig 1 6 6 
0 Deig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
Error in generalev, info= 12 6 -6 

an error occured: dying ... 

Could you please tell me why? 

Thank you 

Sijo GEORGE 

De: "nek5000-users" < [ mailto:nek5000-users at lists.mcs.anl.gov | nek5000-users at lists.mcs.anl.gov ] > 
À: "nek5000-users" < [ mailto:nek5000-users at lists.mcs.anl.gov | nek5000-users at lists.mcs.anl.gov ] > 
Envoyé: Jeudi 31 Mai 2018 14:52:40 
Objet: Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file 

I specify the boundary conditions in the usrdat subroutine using the following piece of code 

do i=1,nelt 
do j=1,2*ndim 
if(bc(5,j,i,1).eq.1) then 
cbc(j,i,1)='mv ' 
elseif(bc(5,j,i,1).eq.2) then 
cbc(j,i,1)='W ' 
elseif(bc(5,j,i,1).eq.3) then 
cbc(j,i,1)='W ' 
elseif(bc(5,j,i,1).eq.4) then 
cbc(j,i,1)='W ' 
else 
cbc(j,i,1)='E ' 
endif 
enddo 
enddo 

Marco 

BQ_BEGIN

On May 31, 2018, at 8:50 AM, [ mailto:nek5000-users at lists.mcs.anl.gov | 
nek5000-users at lists.mcs.anl.gov ] wrote: 

Hi Neks, 

I did the conversion of a mesh from ICEM to exodus file then to re2 file. 

But I am bit confused on the boundary conditions. Because from ICEM before I export to exodus I give the B.Cs as side set numbers and after the conversion to re2 I wont be able to put the boundary conditions in NEK format. 

So could you please tell me how you put the BC in NEK5000 format to a mesh file with exodus format. 

Thank you 
_______________________________________________ 
Nek5000-users mailing list 
[ mailto:Nek5000-users at lists.mcs.anl.gov | Nek5000-users at lists.mcs.anl.gov ] 
[ https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users | https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users ] 

_______________________________________________ 
Nek5000-users mailing list 
[ mailto:Nek5000-users at lists.mcs.anl.gov | Nek5000-users at lists.mcs.anl.gov ] 
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users 
_______________________________________________ 
Nek5000-users mailing list 
[ mailto:Nek5000-users at lists.mcs.anl.gov | Nek5000-users at lists.mcs.anl.gov ] 
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users 

BQ_END

_______________________________________________ 
Nek5000-users mailing list 
Nek5000-users at lists.mcs.anl.gov 
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180531/9bc65cc0/attachment.html>