[Nek5000-users] Generating a 2D circular mesh

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

Hi David,

It's been awhile since I tested this option in genbox.

I would suggest the alternative route of simply
wrapping the box geometry in usrdat2, as is done
for the 2D Taylor Couette example in

    nek5_svn/examples/Taylor

Please let me know if that approach works for you.

Paul



On Sat, 10 Sep 2011, nek5000-users at lists.mcs.anl.gov wrote:

> Hi Neks,
>
> I'm having a little trouble with my first foray into circular geometries. I
> tried to construct a 2D annular mesh using genbox with the following .box
> file:
>
> -----------------------------------------------------------------------
> an.rea
> 2                     spatial dimension
> 2                     number of fields
> Y

vi
> -4  -12          nelr,nel_theta
> 0 0              x0 y0
> ccbb
> 1 2 1
> 0 1 1
> SYM,SYM,         V bc's  ! NB:  3 characters each !
> f  ,f  ,         T bc's  !      You must have 2 spaces!!
> -----------------------------------------------------------------------
>
> genbox seems to work properly, with the exception of the fact that nothing
> after the boundary data gets copied over from the old .rea file.  Nonethess,
> I copy that chunk over manually and run nek5000, which exits with the
> following error:
>
>  WARNINGa: Detected non-right-handed element.
>  Number       1  C1-4:  0.6667E+00  0.8333E+00 -0.2083E+00 -0.4167E-01
>
>
> Any guidance would be much appreciated.
>
> Best,
>
> David
>