<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div><style style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style></div><div><br></div><div>I generated again it didnt work.<br></div><div>I will explain what I did so far.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>I made a cube fluid mesh in ICEM.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Then I gave B.Cs to the cube as follows : outlet = Sideset no : 1, inlet = Sideset no : 5 , Wall_up = Sideset no:4 ,Wall_down = Sideset no:3 , P_Wall_Right =Side_set 6, P_wall_left = Sideset no: 2<br data-mce-bogus="1"></div><div>After that I exported the mesh from ICEM to exodus.exo format file.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Then using exo2nek tool from NEK5000, I converted the above file to exodus.re2 file. (Then I tried to convert into rea file but thrown some errors).<br data-mce-bogus="1"></div><div>Then I mapped the file into exodus.ma2 file.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Then from example Turbulent channel I took a .par file to give parameters.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Now I also given the B.Cs through the .usr file.<br data-mce-bogus="1"></div><div> as follows:<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>c-----------------------------------------------------------------------<br> subroutine usrdat<br> include 'SIZE'<br> include 'TOTAL'<br><br><br> do i=1,nelt<br> do j=1,2*ndim<br> if(bc(5,j,i,1).eq.1) then<br> cbc(j,i,1)='O '<br> elseif(bc(5,j,i,1).eq.2) then<br> cbc(j,i,1)='W '<br> elseif(bc(5,j,i,1).eq.3) then<br> cbc(j,i,1)='W '<br> elseif(bc(5,j,i,1).eq.4) then<br> cbc(j,i,1)='W '<br> elseif(bc(5,j,i,1).eq.5) then<br> cbc(j,i,1)='V '<br> elseif(bc(5,j,i,1).eq.6) then<br> cbc(j,i,1)='W '<br> else<br> cbc(j,i,1)='E '<br> endif<br> enddo<br> enddo<br><br><br> return<br> end<br>c-----------------------------------------------------------------------<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Then in my folder I have .ma2,.re2,.par,.usr files to run this case and when i run the case I get an error shown before.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Could you please tell me, if therer is anything wrong?<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Thank you<br data-mce-bogus="1"></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>De: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br><b>À: </b>"nek5000-users" <nek5000-users@lists.mcs.anl.gov><br><b>Envoyé: </b>Jeudi 31 Mai 2018 15:22:49<br><b>Objet: </b>Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file<br></div><div><br></div><div data-marker="__QUOTED_TEXT__">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Did you generate a new map file after changing the BCs?</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">This matters when switching periodic bcs on or off...</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Paul</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="divRplyFwdMsg" dir="ltr"><span data-mce-style="font-size: 11pt; color: #000000; font-family: Calibri,sans-serif;" style="font-size: 11pt; color: #000000; font-family: Calibri,sans-serif;" color="#000000" face="Calibri, sans-serif"><b>From:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Thursday, May 31, 2018 8:12:55 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file</span>
<div> </div>
</div>
<div class="" style="word-wrap:break-word; line-break:after-white-space">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.<div class=""><br class="">
</div>
<div class="">Marco<br class="">
<div><br class="">
<blockquote class="">
<div class="">On May 31, 2018, at 9:11 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div class="">
<div class="" style="font-family:arial,helvetica,sans-serif; font-size:12pt">
<div class=""><br class="">
</div>
<div class="">Hi Marco,<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I tried and it compiled . But before starting the time step it failed.
<br class="">
</div>
<div class="">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.<br class="">
</div>
<div class="">I get this error:<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Matrix: 0 aa 6 6 6<br class="">
0 aa 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
0 aa -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00<br class="">
0 aa 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00<br class="">
0 aa 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00<br class="">
0 aa 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity<br class="">
0 aa 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity<br class="">
<br class="">
Matrix: 0 bb 6 6 6<br class="">
0 bb 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
0 bb 0.00000E+00 1.66667E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
0 bb 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 4.16667E-02 0.00000E+00 0.00000E+00<br class="">
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 8.33333E-03 0.00000E+00<br class="">
0 bb 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
<br class="">
Matrix: 0 Aeig 6 6 6<br class="">
0 Aeig 8.33333E+01 -4.87842E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
0 Aeig -4.87842E+01 8.66667E+01 -4.87842E+01 7.11752E+00 -1.66667E+00 0.00000E+00<br class="">
0 Aeig 0.00000E+00 -4.87842E+01 8.33333E+01 -4.16667E+01 7.11752E+00 0.00000E+00<br class="">
0 Aeig 0.00000E+00 7.11752E+00 -4.16667E+01 8.33333E+01 -4.87842E+01 0.00000E+00<br class="">
0 Aeig 0.00000E+00 -1.66667E+00 7.11752E+00 -4.87842E+01 Infinity -Infinity<br class="">
0 Aeig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -Infinity Infinity<br class="">
<br class="">
Matrix: 0 Deig 1 6 6<br class="">
0 Deig 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00<br class="">
Error in generalev, info= 12 6 -6<br class="">
<br class="">
an error occured: dying ...<br class="">
<br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Could you please tell me why?<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Thank you<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Sijo GEORGE<br class="">
</div>
<hr id="x_zwchr" class="">
<div class=""><b class="">De: </b>"nek5000-users" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="" target="_blank">nek5000-users@lists.mcs.anl.gov</a>><br class="">
<b class="">À: </b>"nek5000-users" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="" target="_blank">nek5000-users@lists.mcs.anl.gov</a>><br class="">
<b class="">Envoyé: </b>Jeudi 31 Mai 2018 14:52:40<br class="">
<b class="">Objet: </b>Re: [Nek5000-users] How to give NEK5000 boundary conditions to exodus mesh file<br class="">
</div>
<br class="">
<div class="">I specify the boundary conditions in the usrdat subroutine using the following piece of code
<div class=""><br class="">
</div>
<div class="">
<div class=""> do i=1,nelt</div>
<div class=""> do j=1,2*ndim</div>
<div class=""> if(bc(5,j,i,1).eq.1) then</div>
<div class=""> cbc(j,i,1)='mv '</div>
<div class=""> elseif(bc(5,j,i,1).eq.2) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> elseif(bc(5,j,i,1).eq.3) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> elseif(bc(5,j,i,1).eq.4) then</div>
<div class=""> cbc(j,i,1)='W '</div>
<div class=""> else</div>
<div class=""> cbc(j,i,1)='E '</div>
<div class=""> endif</div>
<div class=""> enddo</div>
<div class=""> enddo</div>
</div>
<div class=""><br class="">
</div>
<div class="">Marco</div>
<div class=""><br class="">
<blockquote class="">
<div class="">On May 31, 2018, at 8:50 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div class="">
<div class="" style="font-family:arial,helvetica,sans-serif; font-size:12pt">
<div class="">Hi Neks,<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I did the conversion of a mesh from ICEM to exodus file then to re2 file. </div>
<div class=""><br class="">
</div>
<div class="">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.</div>
<div class=""><br class="">
</div>
<div class="">So could you please tell me how you put the BC in NEK5000 format to a mesh file with exodus format.<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Thank you<br class="">
</div>
</div>
</div>
_______________________________________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br class="">
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" class="" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br class="">
</div>
</blockquote>
</div>
<br class="">
<br class="">
_______________________________________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br class="">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">
</div>
</div>
</div>
_______________________________________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br class="">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
<br>_______________________________________________<br>Nek5000-users mailing list<br>Nek5000-users@lists.mcs.anl.gov<br>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br></div></div></body></html>