<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<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%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><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</font>
<div> </div>
</div>
<meta content="text/html; charset=utf-8">
<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 type="cite" class="">
<div class="">On May 31, 2018, at 9:11 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">
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="">nek5000-users@lists.mcs.anl.gov</a>><br class="">
<b class="">À: </b>"nek5000-users" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">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="">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="">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="">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>
</body>
</html>