[Nek5000-users] Sym boundary condition in mshconvert file

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

Rif,

I would try in your .usr file:

      subroutine usrdat
      :
      :

      integer e,f

      nface = 2*ndim
      do e=1,nelt
      do f=1,nface
          if (cbc(e,f,1).eq.'M  ') cbc(e,f,1) = 'SYM'
      enddo
      enddo

(this is assuming that the "M" gets loaded into the first position of the 3-character
string cbc; you could add:

          if (cbc(e,f,1).eq.' M ') cbc(e,f,1)='SYM'
          if (cbc(e,f,1).eq.'  M') cbc(e,f,1)='SYM'

if you want to cover all bases.

One important point --- "SYM" requires that the face in question to be planar with a normal
pointing in either X Y or Z direction...   If it it pointing somewhere in between you would have
to run with the stress formulation.

Paul

________________________________
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: Thursday, December 31, 2015 1:27 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Sym boundary condition in mshconvert file

Hi Neks,
I have downloaded the program mshconvert file to convert my Fluent mesh to nek5000 format. It works well with normal boundary conditions like W,v and O where W is wall, v is inlet velocity and O is outlet. However, the I am unable to prescribe symmetry, SYM as the code only takes the last letter (for this case, M). Any of you facing this problem and is there a way around this?

Rif
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20151231/b0634e7d/attachment.html>