[Nek5000-users] userchk in turbChannel

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

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> 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 [
> 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
> _______________________________________________
> Nek5000-users mailing list
> 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/9b4bbabc/attachment.html>