[Nek5000-users] n2to3 to create a tube

[nek5000-users at lists.mcs.anl.gov]

Hi Neks,

I am trying to extrude with a rotation a 2d mesh made with gmsh using 
n2to3. If the 2d mesh is a square, n2to3d performs the circular swept. 
However, if the mesh is a bit more complex, let's say two rectangles, 
one bigger than other, the routine n2to3 also works, but genmap gives 
the problem showed below (the 2d mesh perform 2d numerical simulations 
properly, in addition if the extrution is planar instead of circular, 
both n2to3d and genmap works properly and also Nek5000)

Is there anything I am missing about n2to3 and the circular swept? Can I 
rotate this mesh using a different routine and imposing 'P ' boundary 
conditions in the mesh sides?

I have seen the possibility to add something in userchk to extrude the 
mesh, however, what about the periodic boundary conditions? How can I do 
it?

Thank you in advance.
SL

El 12-04-2016 01:38, nek5000-users at lists.mcs.anl.gov escribió:
> Hola SL,
> 
> I met this problem just some days ago, and what I did to solve it via
> making sure the rotating axis is Y-axis !!!
> so you must be aware that your 2D mesh has a Y reflection symmetry.
> then the next step should be fine.
> 
> good luck!
> 
> saludo.
> 
> 2016-04-11 15:57 GMT+02:00 <nek5000-users at lists.mcs.anl.gov>:
> 
>> Hi Neks,
>> 
>> I am trying to use n2to3 to create an axisymmetric mesh. I have
>> created a 2D mesh using gmsh and it works. Now, the idea is to
>> extrude the 2D mesh to create a tube (circular sweep). However,
>> n2to3 doesn't seem to work.
>> 
>> In the 2d file the geometry is rectangular with SYM boundary
>> conditions in the x axis. Is it possible to use n2to3 to create the
>> tube, or half of the tube? If it is possible what are the files that
>> I need? I suppose that with the gemoetry_2d.rea file should be
>> enough, shouldn't it?
>> 
>> I realized that:
>> 
>> If I set the number of planes to, let's say -10, and
>>  input z min:
>> 0
>>  input z max:
>> 90
>> 
>> I have problems with genmap
>> 
>> NOTE: smaller is better, but generous is more forgiving for bad
>> meshes.
>> 0.2
>>  reading .rea file data ...
>>  start locglob_lexico:           8       47680     
>> 381440  0.20000000000000001
>>  locglob:           1           1      381440
>>  locglob:           2         666      381440
>>  locglob:           3        2536      381440
>>  locglob:           1       51469      381440
>>  locglob:           2       51644      381440
>>  locglob:           3       51644      381440
>>  locglob:           1       51644      381440
>>  locglob:           2       51644      381440
>>  locglob:           3       51644      381440
>> 
>>      1     2     4  Matrix:  SELF!!
>>        1 SELF!!    11362    11361    11362    11361
>>        2 SELF!!     9722     9721     9729     9728
>>  cont: SELF!!           1   ??
>> 
>>  ABORT: SELF-CHK            1           5       
>>    1           0
>>  Try to tighten the mesh tolerance!
>> 
>>            0  quit
>> 
>> However, if I increase the number of planes, let's say -80, it
>> seems that genmap also works. But, when I run the numerical
>> simulations and I try to plot some snapshots, I cannot see anything.
>> It seems that there is a problem with the connectivity. What I am
>> missing?
>> 
>> Thanks in advance.
>> Cheers
>> SL
>> 
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> 
> --
> 
> Wei HE 何偉
> 
> School of Aeronautics
> Universidad Politécnica de Madridhttp://www.cfm.upm.es/ [2]
> 
> 
> Links:
> ------
> [1] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> [2] http://www.cfm.upm.es/
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users