[Nek5000-users] Periodic B.C and Genmap Issue
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu Jun 24 13:29:04 CDT 2010
There should normally be two spaces after each bc in genbox.
Regarding format in the .rea file, this is a problem with
mixed characters + numbers on the same input line, particularly
when the input character string supports blanks.
I typically resolve the issue by looking at the source.
In this case, connect2.f in the nek/ repo indicates:
IF (NELGT.LT.1000) THEN
READ(9,50,ERR=500,END=500)
$ CHTEMP,
$ CBC(ISIDE,IEL,IFIELD),ID1,ID2,
$ (BC(II,ISIDE,IEL,IFIELD),II=1,NBCREA)
50 FORMAT(A1,A3,2I3,5G14.6)
ELSEIF (NELGT.LT.100000) THEN
READ(9,51,ERR=500,END=500)
$ CHTEMP,
$ CBC(ISIDE,IEL,IFIELD),ID1,ID2,
$ (BC(II,ISIDE,IEL,IFIELD),II=1,NBCREA)
51 FORMAT(A1,A3,I5,I1,5G14.6)
ELSEIF (NELGT.LT.1000000) THEN
READ(9,52,ERR=500,END=500)
$ CHTEMP,
$ CBC(ISIDE,IEL,IFIELD),ID1,
$ (BC(II,ISIDE,IEL,IFIELD),II=1,NBCREA)
52 FORMAT(A1,A3,I6,5G14.6)
ELSE
...
Paul
On Thu, 24 Jun 2010, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Paul,
>
> I am having a similar problem.
> I am trying to generate a mesh for a 3D periodic channel flow with periodic
> boundary conditions in x and z directions and walls on the y direction.
> I use genbox and then genmap. But i got the same error of Shriram.
> Note that i am having problems when i am using double periodic BCs :
> P ,P ,W ,W ,P ,P
> but single BCs works fine :
> P ,P ,W ,W ,W ,W
> W ,W ,W ,W ,P ,P
> W ,W ,P ,P ,W ,W
>
> Do you have a work around or a suggestion?
>
> Thank you
> francesco
>
>
>
> My error is
>
>
> reading .rea file data ...
> reading .re2 file data ...
> start locglob_lexico: 8 512 4096
> 0.2000000000000000
> locglob: 1 1 4096
> locglob: 2 9 4096
> locglob: 3 81 4096
> locglob: 1 729 4096
> locglob: 2 729 4096
> locglob: 3 729 4096
> locglob: 1 729 4096
> locglob: 2 729 4096
> locglob: 3 729 4096
> done locglob_lexico: 729 729 4096 8
> start periodic vtx: 512 729
> 1 1 5 3 0.00000000E+00 1 shift
>
> abort: PERIODIC MISMATCH 2:
> 1 5 P ie
> 449 6 P je
> 449.0000000000000 0.0000000000000000E+000 bc
>
>
> My .box file is :
>
>
> turbChannel.rea
> 3 spatial dimension
> 1 number of fields
> #
> # comments: This is the box immediately behind the
> # refined cylinder in Ugos cyl+b.l. run.
> #
> #
> #========================================================
> #
> Box 1 Pipe
> -8 -8 -8 Nelx Nely Nelz
> 0.0 6.28 1.0 x0 x1 ratio
> -1.0 1.0 1.0 y0 y1 ratio
> 0.0 3.14 1.0 z0 z1 ratio
> P ,P ,W ,W ,P ,P BCs: (cbx0, cbx1, cby0, cby1, cbz0, cbz1)
>
>
>
>
>
>
>
> -----Original Message-----
> From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov
> Sent: Thu 6/24/2010 6:22 PM
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] Periodic B.C and Genmap Issue
>
> Hi Paul,
>
> This post is in reference to the previous post by Michael (dated Jun 22).
> When I try to get a map file for a rea file that has periodic vertices, I
> get the following error. The rea file was generated the same way as before
> (Gambit -> Nek through Matlab). I checked on Matlab by plotting the periodic
> vertices and confirmed for a few faces that they are indeed periodic. Also,
> as you had earlier suggested (with the new genmap), I tried changing the
> tolerances fro 0.09 to 0.2 and I get the same periodic mismatch error. Could
> you provide any leads on where things are going wrong ? Please find below
> the error message. Thanks.
>
> <error>
> Input (.rea) file name:
> domain
> Input mesh tolerance (default 0.2):
> NOTE: smaller is better, but generous is more forgiving for bad meshes.
> 0.09
> reading .rea file data ...
> start locglob_lexico: 8 3912 31296
> 8.99999999999999967E-002
> locglob: 1 1 31296
> locglob: 2 36 31296
> locglob: 3 1224 31296
> locglob: 1 4851 31296
> locglob: 2 4851 31296
> locglob: 3 4851 31296
> locglob: 1 4851 31296
> locglob: 2 4851 31296
> locglob: 3 4851 31296
> done locglob_lexico: 4851 4851 31296 10
> start periodic vtx: 3912 4851
>
> abort: PERIODIC MISMATCH 2:
> 1093 1 P ie
> 0 1 je
> 4.14470000000000029E-003 0.0000000000000000 bc
>
>
> 8 quit
>
> </error>
>
>
> Regards
> Shriram
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
More information about the Nek5000-users
mailing list