[Nek5000-users] userchk in turbChannel

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

Hi David,

The routine is doing just some sanity checks on your geometry so you shouldn't be concerned. However if you want to import a lot more from that example file it could be that certain routines/operations are hard wired to a geometry generated by genbox (which is the case here). So keep wary and if something doesn't work you need to adapt that to your case.

Oana
________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [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: Friday, September 18, 2015 1:49 PM
To: Nek Mailing List
Subject: Re: [Nek5000-users] userchk in turbChannel


Oana,
Okay, I will keep trying to get around it. Can you tell me what exactly this function is doing? Is it checking for mesh or input errors?
David

On Sep 16, 2015 3:53 PM, <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
Hi David,

I think that piece of code applies to a perfectly rectangular geometry where both nelx, nely are clearly defined as provided in e.g. genbox. For pipe it is indeed a tricky if not impossible matter.

However the check performed there seems to me not too relevant for your problem. You could in fact just compute
nelx*nely*nelz>lelx*lely*lelz

Oana
________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov> [nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov>] on behalf of nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> [nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>]
Sent: Wednesday, September 16, 2015 2:14 PM
To: Nek Mailing List
Subject: [Nek5000-users] userchk in turbChannel

Hi all,

I'm trying to use a round pipe mesh with this LES formulation but this
section in userchk is tripping me up:

      nelx = 8       ! Number of elements in x,y, and z directions.
      nely = 8       ! NOTE, this may vary from one mesh to the next.
      nelz = 8       !

      ntot = nx1*ny1*nz1*nelv

      ! do some checks
      if (istep.eq.0) then
         if(mod(nely,2).ne.0) then
           if(nid.eq.0) write(6,*) 'ABORT: nely has to be even!'
           call exitt
         endif
         if(nelx.gt.lelx .or. nely.gt.lely .or. nelz.gt.lelz) then
           if(nid.eq.0) write(6,*) 'ABORT: nel_xyz > lel_xyz!'
           call exitt
         endif
      endif

How should I determine nelx and nely for the circular area at the end
of the pipe?  I have 56 elements in this area (x-y plane) extruded to
10 layers with n2to3 (560 total elements).

Thanks,

David
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto: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<mailto: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/20150918/2d0abbf5/attachment.html>