<div dir="ltr">Stefan,<div><br></div><div>When I say not in-line, picture an "L" shaped geometry with a cyclic patch on the upper most (north) face and the right most (east) face. Although they point different directions, it should be possible to make them cyclic. I am certain the lengths of my patches are the same size/area.</div><div><br></div><div>I tried adding "T IFCYC" in the rea logical switch but it doesn't appear to be the problem, as I get the same error when I try to use genmap (see initial post). When I use periodic in a prenek channel mesh, I don't get this error. Only when I use periodic with patches that are not in-line.</div><div><br></div><div>Thanks,</div><div>Jon</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jan 13, 2018 at 3:49 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div></div><div>Yes that’s the way to handle cyclic BCs. In addition you have to set IFCYC=true.</div><div><br></div><div>Note the surface mesh on the two cyclic patches have to match. What do you mean by „not in-line“?</div><div><br></div><div>Cheers</div><div>Stefan</div><div><div class="h5"><div><br>On 13 Jan 2018, at 21:35, "<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a>> wrote:<br><br></div><blockquote type="cite"><div><div dir="ltr">Hi Paul,<div>In the future, I will need a more complex geometry where cyclic BCs are not in-line. This simple example I presented here is simply a proof of concept for another mesh.</div><div>What would be the best method do this, when the geometry isn't simple enough to use auto-periodic?</div><div><br></div><div>Thank you,</div><div>Jon</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jan 13, 2018 at 2:28 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_5556649569104790821m_-5858926404772448990divtagdefaultwrapper" 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">Hi Jon,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">It's been awhile since I've run the cyclic BCs that you're trying to use.</p>
<p style="margin-top:0;margin-bottom:0">(Here: cyclic = periodic + turning through an angle.)</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I can't recall the full workflow at the moment (and am not in a position</p>
<p style="margin-top:0;margin-bottom:0">to check), but I think you'll have better luck if you build a straight mesh,</p>
<p style="margin-top:0;margin-bottom:0">run genmap, and then morph the geometry into the desired curved</p>
<p style="margin-top:0;margin-bottom:0">shape in usrdat2().</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">The reason for doing so is that you'll then be guaranteed to have</p>
<p style="margin-top:0;margin-bottom:0">the correct mesh connectivity.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">hth,</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="m_5556649569104790821m_-5858926404772448990divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.m<wbr>cs.anl.gov</a>> on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>><br>
<b>Sent:</b> Saturday, January 13, 2018 1:00:48 PM<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<b>Subject:</b> Re: [Nek5000-users] Periodicity with curvature</font>
<div> </div>
</div><div><div class="m_5556649569104790821h5">
<div>
<div dir="ltr">My apologies, can you clarify what you mean?
<div>I made global refinements one more time, and re-assigned the periodic BCs within prex. For this 2D case (attached), I have linked the 8 element sides with the right 8 sides. All other sides are walls.</div>
<div><br>
</div>
<div>Thank you,<br>
Jon</div>
</div>
<div class="m_5556649569104790821m_-5858926404772448990x_gmail_extra"><br>
<div class="m_5556649569104790821m_-5858926404772448990x_gmail_quote">On Sat, Jan 13, 2018 at 12:49 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>></span> wrote:<br>
<blockquote class="m_5556649569104790821m_-5858926404772448990x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="auto">
<div></div>
<div>can you try again with at least 3 elements in other direction.</div>
<div>
<div class="m_5556649569104790821m_-5858926404772448990x_h5">
<div><br>
On 13 Jan 2018, at 19:42, "<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<div dir="ltr">Thanks for the quick response. I made some global refinements but it seems I'm still getting the same error with genmap.
<div><br>
</div>
<div>Attached is the new .rea file with 4 elements in the periodic direction.</div>
<div><br>
</div>
</div>
<div class="m_5556649569104790821m_-5858926404772448990x_gmail_extra"><br>
<div class="m_5556649569104790821m_-5858926404772448990x_gmail_quote">On Sat, Jan 13, 2018 at 12:32 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>></span> wrote:<br>
<blockquote class="m_5556649569104790821m_-5858926404772448990x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="auto">
<div></div>
<div>You need to have at least three elements in the periodic direction. </div>
<div>
<div class="m_5556649569104790821m_-5858926404772448990x_m_2972450520427343467h5">
<div><br>
On 13 Jan 2018, at 19:19, "<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>> wrote:<br>
<br>
</div>
</div>
</div>
<blockquote type="cite">
<div>
<div dir="ltr">
<div>
<div class="m_5556649569104790821m_-5858926404772448990x_m_2972450520427343467h5">Dear Nek Users,
<div><br>
</div>
<div>I am trying to get periodicity to work in the lastest version of Nek. I've made an extremely simple mesh within prenek (prex) in polar coordinates and cannot seem to run this case. When I run genmap, I get the error below. Did I do something wrong in setting
up this mesh?</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>>> genmap</div>
<div>
<div> reading a.rea</div>
<div>Input mesh tolerance (default 0.2):</div>
<div>NOTE: smaller is better, but generous is more forgiving for bad meshes.</div>
<div>0.01</div>
<div> reading mesh data ...</div>
<div> start locglob_lexico: 4 2 8 1.0000000000000000E-002</div>
<div> locglob: 1 1 8</div>
<div> locglob: 2 5 8</div>
<div> locglob: 1 6 8</div>
<div> locglob: 2 6 8</div>
<div> done locglob_lexico: 6 6 8 2</div>
<div> start periodic vtx: 2 6</div>
<div> 1 4 4 Matrix: x0</div>
<div> 1 x0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324</div>
<div> 2 x0 -3.5355E-01 3.5355E-01 6.9423-310 6.9423-310</div>
<div> 3 x0 3.5355E-01 -3.5355E-01 6.9529-310 6.9529-310</div>
<div> 4 x0 6.9423-310 6.9529-310 2.3477-313 3.1632-317</div>
<div> 2 4 4 Matrix: x1</div>
<div> 1 x1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00</div>
<div> 2 x1 -3.5355E-01 3.5355E-01 0.0000E+00 6.9423-310</div>
<div> 3 x1 -3.5355E-01 3.5355E-01 6.9423-310 6.9529-310</div>
<div> 4 x1 6.9423-310 6.9423-310 6.9423-310 6.9529-310</div>
<div> 1 4 4 Matrix: z0</div>
<div> 1 z0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324</div>
<div> 2 z0 -2.1213E+00 -1.4142E+00 6.9423-310 6.9423-310</div>
<div> 3 z0 2.1213E+00 1.4142E+00 6.9529-310 6.9529-310</div>
<div> 4 z0 6.9423-310 6.9529-310 2.3477-313 3.1632-317</div>
<div> 2 4 4 Matrix: z1</div>
<div> 1 z1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00</div>
<div> 2 z1 1.4142E+00 2.1213E+00 0.0000E+00 6.9423-310</div>
<div> 3 z1 1.4142E+00 2.1213E+00 6.9423-310 6.9529-310</div>
<div> 4 z1 6.9423-310 6.9423-310 6.9423-310 6.9529-310</div>
<div> 1 3 0 1.00000000E-03 2.12132000E+00 abort: FACE MATCH FAIL</div>
<div> 2 4 3 2.12132000E-03 1.00000031E+00 abort: FACE MATCH FAIL</div>
<div><br>
</div>
<div> 0 quit</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>The mesh within my .rea file is as follows:</div>
<div>
<div> 20.0000 20.0000 -8.00000 -4.00000 XFAC,YFAC,XZERO,YZERO</div>
<div> **MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.</div>
<div> 2 2 2 NEL,NDIM,NELV</div>
<div> ELEMENT 1 [ 1A] GROUP 0</div>
<div> -2.121320 -1.414213 -0.8742276E-07 -0.1311341E-06</div>
<div> 2.121320 1.414213 2.000000 3.000000</div>
<div> ELEMENT 2 [ 1B] GROUP 0</div>
<div> -0.1311341E-06 -0.8742276E-07 1.414213 2.121320</div>
<div> 3.000000 2.000000 1.414213 2.121320</div>
<div> ***** CURVED SIDE DATA *****</div>
<div> 4 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE</div>
<div> 2 1 -2.00000 0.00000 0.00000 0.00000 0.00000 C</div>
<div> 4 1 3.00000 0.00000 0.00000 0.00000 0.00000 C</div>
<div> 2 2 -2.00000 0.00000 0.00000 0.00000 0.00000 C</div>
<div> 4 2 3.00000 0.00000 0.00000 0.00000 0.00000 C</div>
<div> ***** BOUNDARY CONDITIONS *****</div>
<div> ***** FLUID BOUNDARY CONDITIONS *****</div>
<div> P 1 1 2.00000 3.00000 0.00000 0.00000 0.00000</div>
<div> W 1 2 0.00000 0.00000 0.00000 0.00000 0.00000</div>
<div> E 1 3 2.00000 1.00000 0.00000 0.00000 0.00000</div>
<div> W 1 4 0.00000 0.00000 0.00000 0.00000 0.00000</div>
<div> E 2 1 1.00000 3.00000 0.00000 0.00000 0.00000</div>
<div> W 2 2 0.00000 0.00000 0.00000 0.00000 0.00000</div>
<div> P 2 3 1.00000 1.00000 0.00000 0.00000 0.00000</div>
<div> W 2 4 0.00000 0.00000 0.00000 0.00000 0.00000</div>
<div> ***** NO THERMAL BOUNDARY CONDITIONS *****</div>
</div>
<div><br>
</div>
</div>
</div>
<div><mime-attachment.png><br>
</div>
</div>
<pre>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=ncKzoPY2dwLIbAc1FITAEJsCsTeH2zRgDAqkjXYZe5A&s=nLM4v-0kL5jW9dMNEgiVqdVGDFSs6TFHQxJK53ePw04&e=" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
</pre>
</div>
</blockquote>
</div>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=BAFsct9orEUodBvhKrXssNMVdu-OcqgveZCrA11i8M8&s=DBQt47t5SSmsVmTy3GCoo7VqaiTms-bCOKWC0nrs9xo&e=" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
</blockquote>
</div>
<br>
</div>
<pre>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=BAFsct9orEUodBvhKrXssNMVdu-OcqgveZCrA11i8M8&s=DBQt47t5SSmsVmTy3GCoo7VqaiTms-bCOKWC0nrs9xo&e=" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
</pre>
</div>
</blockquote>
</div>
</div>
<blockquote type="cite">
<div><a.rea></div>
</blockquote>
</div>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFAg&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=sAFtxxjU2oId9i7igsGJ8tVajh1Q5GHP-CQMMPawlcs&s=NLpzCQaEmDTGEe8B503wySfumBz2BxWtkf6eEvhyGAU&e=" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div></div></div>
<br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=Bb2NwYJSr6FKn4V5QWP-feJQ7TDXc4Hz3edNRRRkb_s&s=yjFBkT1sIfkT1-SUIWsA1IJZsfv58lXi1UmuKM9S4Rw&e=" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>
<pre>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.mcs.anl.gov_mailman_listinfo_nek5000-2Dusers&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=gP2CkfxBOr64FG89YpFokg&m=Bb2NwYJSr6FKn4V5QWP-feJQ7TDXc4Hz3edNRRRkb_s&s=yjFBkT1sIfkT1-SUIWsA1IJZsfv58lXi1UmuKM9S4Rw&e=" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
</pre></div></blockquote></div></div></div><br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>