Hi Nek's,<br><br>This question actually concerns the matlab code to transfer mesh from gambit to Nek. So far, I have used it to mesh pipe flow with a given inlet and oulet and it worked perfectly fine. I would like now to use periodic conditions between the inlet and oulet. As a matter of fact, I have remeshed my pipe with gambit and created 2 set of boundary conditions:<br>
<ol><li>the pipe wall</li><li>a set called 'periodic' that contains all the faces from the inlet and outlet</li></ol><p>I then export the mesh under the .neu format. The thing is when I use the matlab code to produce my pipe.rea file, here is what I get for any faces subject to the periodic condition:</p>
<p>P 1 5 0 0.00000 0.00000 0.00000 0.00000</p><p>The 5th face of element 1 should be connected to another element, but nothing comes out the matlab code. Have I missed something while meshing or is there any specific thing to do with the matlab code to be able to handle such periodic boundary conditions?</p>
<p>Sincerely yours,</p><p>JC<br></p><div class="gmail_quote">On 24 June 2010 20:29, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
There should normally be two spaces after each bc in genbox.<br>
<br>
Regarding format in the .rea file, this is a problem with<br>
mixed characters + numbers on the same input line, particularly<br>
when the input character string supports blanks.<br>
<br>
I typically resolve the issue by looking at the source.<br>
In this case, connect2.f in the nek/ repo indicates:<br>
<br>
<br>
IF (NELGT.LT.1000) THEN<br>
READ(9,50,ERR=500,END=500)<br>
$ CHTEMP,<br>
$ CBC(ISIDE,IEL,IFIELD),ID1,ID2,<br>
$ (BC(II,ISIDE,IEL,IFIELD),II=1,<u></u>NBCREA)<br>
50 FORMAT(A1,A3,2I3,5G14.6)<br>
ELSEIF (NELGT.LT.100000) THEN<br>
READ(9,51,ERR=500,END=500)<br>
$ CHTEMP,<br>
$ CBC(ISIDE,IEL,IFIELD),ID1,ID2,<br>
$ (BC(II,ISIDE,IEL,IFIELD),II=1,<u></u>NBCREA)<br>
51 FORMAT(A1,A3,I5,I1,5G14.6)<br>
ELSEIF (NELGT.LT.1000000) THEN<br>
READ(9,52,ERR=500,END=500)<br>
$ CHTEMP,<br>
$ CBC(ISIDE,IEL,IFIELD),ID1,<br>
$ (BC(II,ISIDE,IEL,IFIELD),II=1,<u></u>NBCREA)<br>
52 FORMAT(A1,A3,I6,5G14.6)<br>
ELSE<br>
<br>
<br>
...<br>
<br>
<br>
Paul<div class="im HOEnZb"><br>
<br>
<br>
<br>
On Thu, 24 Jun 2010, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
</div><div class="HOEnZb"><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Paul,<br>
<br>
I am having a similar problem.<br>
I am trying to generate a mesh for a 3D periodic channel flow with periodic<br>
boundary conditions in x and z directions and walls on the y direction.<br>
I use genbox and then genmap. But i got the same error of Shriram.<br>
Note that i am having problems when i am using double periodic BCs :<br>
P ,P ,W ,W ,P ,P<br>
but single BCs works fine :<br>
P ,P ,W ,W ,W ,W<br>
W ,W ,W ,W ,P ,P<br>
W ,W ,P ,P ,W ,W<br>
<br>
Do you have a work around or a suggestion?<br>
<br>
Thank you<br>
francesco<br>
<br>
<br>
<br>
My error is<br>
<br>
<br>
reading .rea file data ...<br>
reading .re2 file data ...<br>
start locglob_lexico: 8 512 4096<br>
0.2000000000000000<br>
locglob: 1 1 4096<br>
locglob: 2 9 4096<br>
locglob: 3 81 4096<br>
locglob: 1 729 4096<br>
locglob: 2 729 4096<br>
locglob: 3 729 4096<br>
locglob: 1 729 4096<br>
locglob: 2 729 4096<br>
locglob: 3 729 4096<br>
done locglob_lexico: 729 729 4096 8<br>
start periodic vtx: 512 729<br>
1 1 5 3 0.00000000E+00 1 shift<br>
<br>
abort: PERIODIC MISMATCH 2:<br>
1 5 P ie<br>
449 6 P je<br>
449.0000000000000 0.0000000000000000E+000 bc<br>
<br>
<br>
My .box file is :<br>
<br>
<br>
turbChannel.rea<br>
3 spatial dimension<br>
1 number of fields<br>
#<br>
# comments: This is the box immediately behind the<br>
# refined cylinder in Ugo’s cyl+b.l. run.<br>
#<br>
#<br>
#=============================<u></u>===========================<br>
#<br>
Box 1 Pipe<br>
-8 -8 -8 Nelx Nely Nelz<br>
0.0 6.28 1.0 x0 x1 ratio<br>
-1.0 1.0 1.0 y0 y1 ratio<br>
0.0 3.14 1.0 z0 z1 ratio<br>
P ,P ,W ,W ,P ,P BC’s: (cbx0, cbx1, cby0, cby1, cbz0, cbz1)<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.<u></u>mcs.anl.gov</a> on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
Sent: Thu 6/24/2010 6:22 PM<br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
Subject: [Nek5000-users] Periodic B.C and Genmap Issue<br>
<br>
Hi Paul,<br>
<br>
This post is in reference to the previous post by Michael (dated Jun 22).<br>
When I try to get a map file for a rea file that has periodic vertices, I<br>
get the following error. The rea file was generated the same way as before<br>
(Gambit -> Nek through Matlab). I checked on Matlab by plotting the periodic<br>
vertices and confirmed for a few faces that they are indeed periodic. Also,<br>
as you had earlier suggested (with the new genmap), I tried changing the<br>
tolerances fro 0.09 to 0.2 and I get the same periodic mismatch error. Could<br>
you provide any leads on where things are going wrong ? Please find below<br>
the error message. Thanks.<br>
<br>
<error><br>
Input (.rea) file name:<br>
domain<br>
Input mesh tolerance (default 0.2):<br>
NOTE: smaller is better, but generous is more forgiving for bad meshes.<br>
0.09<br>
reading .rea file data ...<br>
start locglob_lexico: 8 3912 31296<br>
8.99999999999999967E-002<br>
locglob: 1 1 31296<br>
locglob: 2 36 31296<br>
locglob: 3 1224 31296<br>
locglob: 1 4851 31296<br>
locglob: 2 4851 31296<br>
locglob: 3 4851 31296<br>
locglob: 1 4851 31296<br>
locglob: 2 4851 31296<br>
locglob: 3 4851 31296<br>
done locglob_lexico: 4851 4851 31296 10<br>
start periodic vtx: 3912 4851<br>
<br>
abort: PERIODIC MISMATCH 2:<br>
1093 1 P ie<br>
0 1 je<br>
4.14470000000000029E-003 0.0000000000000000 bc<br>
<br>
<br>
8 quit<br>
<br>
</error><br>
<br>
<br>
Regards<br>
Shriram<br>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote>
</div></div><br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Jean-Christophe<br>