[Nek5000-users] pipe prenek (nek5000-users at lists.mcs.anl.gov)

Wed Oct 18 08:57:48 CDT 2017

Hi ZJ,

I have written a python script to generate a mesh for pipe flow which 
should follow the nek format of the rea files. Feel free to use it if 
you want to.

You can use "git clone 
https://git.scc.kit.edu/nek5000_scripts/preprocessing" to clone my 
repository.
Then go into the folder "GenPipeMesh" and follow the instructions in the 
README.md.

As I have just made it public, you might run into trouble when cloning 
the repo. Moreover, consider the code as a work in progress. So there 
might be some mesh settings that do not work correctly.
If you have problems cloning the repo or generating the mesh, just let 
me know and I will see what I can do (steffen.straub at kit.edu).

HTH,

Steffen

On 10/18/2017 02:42 PM, nek5000-users-request at lists.mcs.anl.gov wrote:
> Message: 2
> Date: Wed, 18 Oct 2017 12:18:45 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
> 	<nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] pipe prenek
> Message-ID:
> 	<mailman.673.1508329131.19102.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="us-ascii"
>
>
> Dear ZJ,
>
>
> In order to parse the characters in the line, nek uses formatted reads for
>
> curve sides and for BC info.
>
>
> The format changes as the number of elements increases in order to accommodate the larger integers (i.e., the element number) in the same line,
>
> so you have to check either the reader source or prenek source or a valid
>
> .rea file to see what that formatting looks like.
>
>
> I've fixed your pipe mesh here.  Note that there is a 2D pipe mesh (but with a different base template) in NekExamples/pipe/w2dcyl020a.rea
>
>
> 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: Wednesday, October 18, 2017 4:48:38 AM
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] pipe prenek
>
> Hi Nek,
>
> Is there anyone knows how to calculate the curvature in .rea file. Recently, I want to generate a pipe mesh but it doesn't seem like a pipe in which I guess there is something wrong in curvature part in .rea file. I have attached the mesh part below.
>
>    **MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
>          12         2        12 NEL,NDIM,NELV
>        ELEMENT          1 [    1A]    GROUP     0
>     0.000000      0.350000       0.250000       0.000000
>     0.000000      0.000000       0.250000       0.350000
>        ELEMENT          2 [    1B]    GROUP     0
>     0.500000      0.353554       0.250000       0.350000
>     0.000000      0.353554       0.250000       0.000000
>        ELEMENT          3 [    1C]    GROUP     0
>     0.000000      0.000000       0.250000       0.353554
>     0.500000      0.350000       0.250000       0.353554
>        ELEMENT          4 [    1A]    GROUP     0
>     0.000000      0.000000       0.250000       0.350000
>     0.000000     -0.350000      -0.250000       0.000000
>        ELEMENT          5 [    1B]    GROUP     0
>     0.000000      0.353554       0.250000       0.000000
>    -0.500000     -0.353554      -0.250000      -0.350000
>        ELEMENT          6 [    1C]    GROUP     0
>     0.500000      0.350000       0.250000       0.353554
>     0.000000      0.000000      -0.250000      -0.353554
>        ELEMENT          7 [    1A]    GROUP     0
>     0.000000     -0.350000      -0.250000       0.000000
>     0.000000      0.000000      -0.250000      -0.350000
>        ELEMENT          8 [    1B]    GROUP     0
>    -0.500000     -0.353554      -0.250000      -0.350000
>     0.000000     -0.353554      -0.250000       0.000000
>        ELEMENT          9 [    1C]    GROUP     0
>     0.000000      0.000000      -0.250000      -0.353554
>    -0.500000     -0.350000      -0.250000      -0.353554
>        ELEMENT         10 [    1A]    GROUP     0
>     0.000000      0.000000      -0.250000      -0.350000
>     0.000000      0.350000       0.250000       0.000000
>        ELEMENT         11 [    1B]    GROUP     0
>     0.000000     -0.353554      -0.250000       0.000000
>     0.500000      0.353554       0.250000       0.350000
>        ELEMENT         12 [    1C]    GROUP     0
>    -0.500000     -0.350000      -0.250000      -0.353554
>     0.000000      0.000000       0.250000       0.353554
>    ***** CURVED SIDE DATA *****
>           8 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
>    1  2  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    4  3  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    1  5  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    4  6  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    1  8  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    4  9  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    1 11  0.50000   0.00000       0.00000       0.00000       0.00000     C
>    4 12  0.50000   0.00000       0.00000       0.00000       0.00000     C
>
>
> Kind regards,
>
> ZJ
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171018/fbc775d5/attachment-0001.html>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: b.rea.gz
> Type: application/x-gzip
> Size: 2667 bytes
> Desc: b.rea.gz
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171018/fbc775d5/attachment-0001.gz>
>
> ------------------------------