[Nek5000-users] RE : problem with closed internal boundaries

Tue May 11 18:20:06 CDT 2010

Hi Michael,

I just ran your case and all seems OK.  I built a new *.map but it
turned out to be the same as yours.  Everything is at:
http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/temp/a

Cheers,
Frank

On Wed, 2010-05-12 at 00:49 +0200, nek5000-users at lists.mcs.anl.gov
wrote:
> Hi Stefan, 
> 
> as compilers, I use mpif77 and mpicc for the makenek, and gfortran and mpicc for the maketools (assuming that prenek and/or genmap compilation could have an influence here).
> 
> Frank, indeed my whole geometry, as interpreted by nek5000, is messed up. That's the point but I have no clue where this problem is coming from.
> 
> Cheers,
> 
> 
> Michael 
> ________________________________________
> De : nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] de la part de nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
> Date d'envoi : mardi 11 mai 2010 21:10
> À : nek5000-users at lists.mcs.anl.gov
> Objet : Re: [Nek5000-users] RE :  problem with closed internal boundaries
> 
> Hi Michael,
> 
> It seems that perhaps your mesh and/or grid mapping is messed up.  The
> numbers concerning your mesh in your "out" file are huge:
> 
>  verify mesh topology
>   9.90000000000000105E+021 -9.90000000000000105E+021  Xrange
>   9.90000000000000105E+021 -9.90000000000000105E+021  Yrange
>   9.90000000000000105E+021 -9.90000000000000105E+021  Zrange
>  done :: verify mesh topology
> 
> 
> >From one of my simulations:
> 
> 
>  verify mesh topology
>   -1.51515151515152        1.51515151515152       Xrange
>   -1.51515151515152        1.51515151515152       Yrange
>  -0.500000000000000       0.500000000000000       Zrange
>  done :: verify mesh topology
> 
>  
> Cheers,
> Frank
> 
> 
> On Tue, 2010-05-11 at 20:57 +0200, nek5000-users at lists.mcs.anl.gov
> wrote:
> > Hello Frank,
> >
> > thanks for the quick answer. The Nek team should be relieved that the community is growing and that some users start to be able to help the newbies.
> > I tried what you recommended and now that you've said it, it looks obvious that all the z parameters should be 1 for a 2D case, but this isn't changing a thing. I still get those NaNs.
> >
> > Do you have any other idea?
> >
> > Best,
> >
> >
> > Michael
> > ________________________________________
> > De : nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] de la part de nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
> > Date d'envoi : mardi 11 mai 2010 20:31
> > À : nek5000-users at lists.mcs.anl.gov
> > Objet : Re: [Nek5000-users] problem with closed internal boundaries
> >
> > Hello Michael,
> >
> > I am by no means knowledgeable concerning NEK, but will put my two-cents
> > in by saying that I believe that all dimensions in the z-direction must
> > be 1.  However, you have "parameter (lz2=lz1-2)" and "parameter
> > (lz3=lz2)".
> >
> > Cheers,
> > Frank
> >
> >
> > On Tue, 2010-05-11 at 18:20 +0200, nek5000-users at lists.mcs.anl.gov
> > wrote:
> > > Hello Neks,
> > >
> > > I'm trying to model a flow past an obstacle in a channel. As a first step I would like to do a 2D simulation.
> > > I'm encountering some problems during the geometry data generation in the nek5000 execution : it looks like the program can only generate NaN values although genmap seems to work well on the .rea file.
> > > I can also reopen the rea file with prenek without any problems.
> > >
> > > Do you have any idea what is going wrong here?
> > > Is it problematic to have a closed internal boundary condition?
> > > Thanks in advance for any help.
> > >
> > > Best,
> > >
> > >
> > > Michael B
> > > _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> > --
> > Frank Herbert Muldoon, Ph.D. Mechanical Engineering
> > Technische Universität Wien (Technical University of Vienna)
> > Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
> > Mechanics and Heat Transfer)
> > Resselgasse 3
> > 1040 Wien
> > Tel: +4315880132232
> > Fax: +4315880132299
> > Cell:+436765203470
> > fmuldoo (skype)
> > http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html
> >
> > _______________________________________________
> > Nek5000-users mailing list
> > Nek5000-users at lists.mcs.anl.gov
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> > _______________________________________________
> > Nek5000-users mailing list
> > Nek5000-users at lists.mcs.anl.gov
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> --
> Frank Herbert Muldoon, Ph.D. Mechanical Engineering
> Technische Universität Wien (Technical University of Vienna)
> Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
> Mechanics and Heat Transfer)
> Resselgasse 3
> 1040 Wien
> Tel: +4315880132232
> Fax: +4315880132299
> Cell:+436765203470
> fmuldoo (skype)
> http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-- 
Frank Herbert Muldoon, Ph.D. Mechanical Engineering
Technische Universität Wien (Technical University of Vienna)
Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
Mechanics and Heat Transfer)
Resselgasse 3
1040 Wien
Tel: +4315880132232
Fax: +4315880132299 
Cell:+436765203470
fmuldoo (skype)
http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html