<div dir="ltr">Thanks a lot Marco!<br></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Nov 14, 2017 at 11:27 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
It works for 2D and 3D meshes. Sideset ids should be set in the package you use to create the mesh. Then in the subroutine usrdat (inside the *.usr file), you will have to loop over the faces of each element and set the boundary conditions:
<div><br>
</div>
<div>
<div> do i=1,nelt</div>
<div> do j=1,2*ndim</div>
<div> !write(6,*) bc(5,j,i,1)</div>
<div> if(bc(5,j,i,1).eq.1) then</div>
<div> cbc(j,i,1)='mv '</div>
<div> !write(6,*) 'mv bc'</div>
<div> elseif(bc(5,j,i,1).eq.2) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.3) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.4) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.5) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.6) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.7) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.8) then</div>
<div> cbc(j,i,1)='W '</div>
<div> elseif(bc(5,j,i,1).eq.9) then</div>
<div> cbc(j,i,1)='W '</div>
<div> else</div>
<div> cbc(j,i,1)='E '</div>
<div> endif</div>
<div> enddo</div>
<div> enddo</div>
</div>
<div><br>
<div>The installation of the exo2nek script is not that complicated if you follow the instructions in the README file.</div>
<div><br>
</div>
<div>Marco</div></div></div><div style="word-wrap:break-word"><div>
<div><br>
</div>
<div><br>
<div>
<blockquote type="cite">
<div>On Nov 14, 2017, at 10:15 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="m_-154981083632584073Apple-interchange-newline">
<div>
<div dir="ltr">
<div>Does exo2nek support 3D conversion? Also, someone should consider updating the README for the exo2nek tools. Its not quite clear for someone who is starting out with nek5000, what exactly sideset ids are and how fluidic bcs are actually assigned.
I will however give it a go. <br>
<br>
</div>
Thanks for your help everyone.<br>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr">On Tue, Nov 14, 2017 at 9:39 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">can you generate the same mesh in a exodus format? If you can, use the exo2nek script to convert the mesh in a re2 file.
<div><br>
</div>
<div>Marco</div>
</div>
<div style="word-wrap:break-word">
<div><br>
<div>
<blockquote type="cite">
<div>On Nov 14, 2017, at 9:36 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="m_-154981083632584073m_-415075026871470738Apple-interchange-newline">
<div>
<div dir="auto">
<div></div>
<div>Note, this converter is not part of Nek5000 and not support by us. Yes this converter does not output the correct BC. </div>
<div><br>
On 14 Nov 2017, at 08:13, "<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>So your mesh has more than 1000 elements. I think that the way the python script is written the boundary conditions in a rea file is not correct. I had this issue in the past.<br>
<div>
<blockquote type="cite">
<div>On Nov 14, 2017, at 9:02 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="m_-154981083632584073m_-415075026871470738Apple-interchange-newline">
<div>
<div dir="ltr">Its a 25X25 square grid in 2D. Im trying to extrude this with 10 elements in the z direction.<br>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr">On Tue, Nov 14, 2017 at 8:13 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">What is the number of elements in your mesh?</div>
<div style="word-wrap:break-word"><br>
<div>
<blockquote type="cite">
<div>On Nov 13, 2017, at 11:40 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="m_-154981083632584073m_-415075026871470738m_-4113089898437144111Apple-interchange-newline">
<div>
<p dir="ltr">Gmsh to .rea using CPraveen's python mesh converter.</p>
<br>
<div class="gmail_quote">
<div dir="ltr">On Mon, Nov 13, 2017, 10:52 PM <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Looks like the boundary conditions in your .rea file are somehow messed up. How do you generate you .rea?<br>
<br>
-----Original message-----<br>
> <a href="mailto:From%3Anek5000-users@lists.mcs.anl.gov" target="_blank">
From:nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>><br>
> Sent: Tuesday 14th November 2017 4:47<br>
> To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
> Subject: [Nek5000-users] n2to3 mesh size error<br>
><br>
> Hi all,<br>
><br>
> I keep running into an error whenever I want to use n2to3 to extrude my mesh from 2 dimensions to 3. My workflow is as follows:<br>
><br>
> gmsh .geo file -> gmsh .msh file -> .rea file -> n2to3 > .rea (3D) file.<br>
><br>
> Now if the original 2D .rea file obtained from gmsh contains less than 1000 elements, n2to3 fails with the error<br>
><br>
> At line 651 of file n2to3.f (unit = 10, file = msh_test.rea)<br>
> Fortran runtime error: Bad value during integer read<br>
><br>
> However, should this not be the case (mesh size > 1000) its all smooth sailing. The trouble is with a mesh of that size my processor and memory requirements quickly blow up.<br>
><br>
> Incidentally, genmap fails under almost exactly the same conditions.<br>
><br>
> Is there a lower limit to mesh size issue for nek5000 that Im missing here?<br>
><br>
> Sincerely,<br>
> --<br>
> Amitvikram Dutta<br>
> Graduate Research Assistant<br>
> Fluid Mechanics Research Lab<br>
> University of Waterloo<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> <a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
</div>
<div dir="ltr">-- <br>
</div>
<div class="m_-154981083632584073m_-415075026871470738m_-4113089898437144111gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p>
<p dir="ltr">University of Waterloo</p>
</div>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div>
</blockquote>
</div>
<br>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
</div>
<div dir="ltr">-- <br>
</div>
<div class="m_-154981083632584073m_-415075026871470738gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p>
<p dir="ltr">University of Waterloo</p>
</div>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div>
</blockquote>
</div>
<br>
<pre>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</pre>
</div>
</blockquote>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div>
</blockquote>
</div>
<br>
</div>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote>
</div>
<div dir="ltr">-- <br>
</div>
<div class="m_-154981083632584073gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p>
<p dir="ltr">University of Waterloo</p>
</div>
</div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div>
</blockquote>
</div>
<br>
</div>
</div></div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><div dir="ltr">-- <br></div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p>
<p dir="ltr">University of Waterloo</p>
</div></div>