[Nek5000-users] Periodicity with curvature
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sat Jan 13 14:34:58 CST 2018
Hi Paul,
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.
What would be the best method do this, when the geometry isn't simple
enough to use auto-periodic?
Thank you,
Jon
On Sat, Jan 13, 2018 at 2:28 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
> Hi Jon,
>
>
> It's been awhile since I've run the cyclic BCs that you're trying to use.
>
> (Here: cyclic = periodic + turning through an angle.)
>
>
> I can't recall the full workflow at the moment (and am not in a position
>
> to check), but I think you'll have better luck if you build a straight
> mesh,
>
> run genmap, and then morph the geometry into the desired curved
>
> shape in usrdat2().
>
>
> The reason for doing so is that you'll then be guaranteed to have
>
> the correct mesh connectivity.
>
>
> hth,
>
>
> Paul
>
>
> ------------------------------
> *From:* Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf
> of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
> *Sent:* Saturday, January 13, 2018 1:00:48 PM
> *To:* nek5000-users at lists.mcs.anl.gov
> *Subject:* Re: [Nek5000-users] Periodicity with curvature
>
> My apologies, can you clarify what you mean?
> 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.
>
> Thank you,
> Jon
>
> On Sat, Jan 13, 2018 at 12:49 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
> can you try again with at least 3 elements in other direction.
>
> On 13 Jan 2018, at 19:42, "nek5000-users at lists.mcs.anl.gov" <
> nek5000-users at lists.mcs.anl.gov> wrote:
>
> Thanks for the quick response. I made some global refinements but it seems
> I'm still getting the same error with genmap.
>
> Attached is the new .rea file with 4 elements in the periodic direction.
>
>
> On Sat, Jan 13, 2018 at 12:32 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
> You need to have at least three elements in the periodic direction.
>
> On 13 Jan 2018, at 19:19, "nek5000-users at lists.mcs.anl.gov" <
> nek5000-users at lists.mcs.anl.gov> wrote:
>
> Dear Nek Users,
>
> 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?
>
>
>
> >> genmap
> reading a.rea
> Input mesh tolerance (default 0.2):
> NOTE: smaller is better, but generous is more forgiving for bad meshes.
> 0.01
> reading mesh data ...
> start locglob_lexico: 4 2 8
> 1.0000000000000000E-002
> locglob: 1 1 8
> locglob: 2 5 8
> locglob: 1 6 8
> locglob: 2 6 8
> done locglob_lexico: 6 6 8 2
> start periodic vtx: 2 6
> 1 4 4 Matrix: x0
> 1 x0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324
> 2 x0 -3.5355E-01 3.5355E-01 6.9423-310 6.9423-310
> 3 x0 3.5355E-01 -3.5355E-01 6.9529-310 6.9529-310
> 4 x0 6.9423-310 6.9529-310 2.3477-313 3.1632-317
> 2 4 4 Matrix: x1
> 1 x1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00
> 2 x1 -3.5355E-01 3.5355E-01 0.0000E+00 6.9423-310
> 3 x1 -3.5355E-01 3.5355E-01 6.9423-310 6.9529-310
> 4 x1 6.9423-310 6.9423-310 6.9423-310 6.9529-310
> 1 4 4 Matrix: z0
> 1 z0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324
> 2 z0 -2.1213E+00 -1.4142E+00 6.9423-310 6.9423-310
> 3 z0 2.1213E+00 1.4142E+00 6.9529-310 6.9529-310
> 4 z0 6.9423-310 6.9529-310 2.3477-313 3.1632-317
> 2 4 4 Matrix: z1
> 1 z1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00
> 2 z1 1.4142E+00 2.1213E+00 0.0000E+00 6.9423-310
> 3 z1 1.4142E+00 2.1213E+00 6.9423-310 6.9529-310
> 4 z1 6.9423-310 6.9423-310 6.9423-310 6.9529-310
> 1 3 0 1.00000000E-03 2.12132000E+00 abort: FACE MATCH FAIL
> 2 4 3 2.12132000E-03 1.00000031E+00 abort: FACE MATCH FAIL
>
> 0 quit
>
>
>
>
>
> The mesh within my .rea file is as follows:
> 20.0000 20.0000 -8.00000 -4.00000
> XFAC,YFAC,XZERO,YZERO
> **MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
> 2 2 2 NEL,NDIM,NELV
> ELEMENT 1 [ 1A] GROUP 0
> -2.121320 -1.414213 -0.8742276E-07 -0.1311341E-06
> 2.121320 1.414213 2.000000 3.000000
> ELEMENT 2 [ 1B] GROUP 0
> -0.1311341E-06 -0.8742276E-07 1.414213 2.121320
> 3.000000 2.000000 1.414213 2.121320
> ***** CURVED SIDE DATA *****
> 4 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
> 2 1 -2.00000 0.00000 0.00000 0.00000 0.00000
> C
> 4 1 3.00000 0.00000 0.00000 0.00000 0.00000
> C
> 2 2 -2.00000 0.00000 0.00000 0.00000 0.00000
> C
> 4 2 3.00000 0.00000 0.00000 0.00000 0.00000
> C
> ***** BOUNDARY CONDITIONS *****
> ***** FLUID BOUNDARY CONDITIONS *****
> P 1 1 2.00000 3.00000 0.00000 0.00000
> 0.00000
> W 1 2 0.00000 0.00000 0.00000 0.00000
> 0.00000
> E 1 3 2.00000 1.00000 0.00000 0.00000
> 0.00000
> W 1 4 0.00000 0.00000 0.00000 0.00000
> 0.00000
> E 2 1 1.00000 3.00000 0.00000 0.00000
> 0.00000
> W 2 2 0.00000 0.00000 0.00000 0.00000
> 0.00000
> P 2 3 1.00000 1.00000 0.00000 0.00000
> 0.00000
> W 2 4 0.00000 0.00000 0.00000 0.00000
> 0.00000
> ***** NO THERMAL BOUNDARY CONDITIONS *****
>
> <mime-attachment.png>
>
> _______________________________________________
>
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <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=>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> <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=>
>
>
> _______________________________________________
>
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <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=>
>
> <a.rea>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> <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=>
>
>
>
> _______________________________________________
> 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/20180113/b9ce137b/attachment-0001.html>
More information about the Nek5000-users
mailing list