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

[nek5000-users at lists.mcs.anl.gov]

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" <nek5000-users at lists.mcs.anl.gov> 
À: "nek5000-users" <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 




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 





_______________________________________________ 
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/9f21871d/attachment.html>