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
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: [email protected] [[email protected]] on behalf of [email protected] [[email protected]] Sent: Thursday, December 31, 2015 1:27 AM To: [email protected] 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
participants (1)
-
nek5000-users@lists.mcs.anl.gov