[Nek5000-users] userchk in turbChannel

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

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 [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, 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
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users