<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Hello,</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span>Yes, i will try, thank you very much for you help.</span></div><div><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Best regards,</div><div style="color: rgb(0, 0,
0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Dorin</div> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> Christophe Geuzaine <cgeuzaine@ulg.ac.be><br> <b><span style="font-weight: bold;">To:</span></b> Dorin Niculescu <niculescu_dorin85@yahoo.com> <br><b><span style="font-weight: bold;">Cc:</span></b> F. A. Portela <felipe.alves.portela@gmail.com>; "gmsh@geuz.org" <gmsh@geuz.org> <br> <b><span style="font-weight: bold;">Sent:</span></b>
Friday, 5 April 2013, 9:57<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [Gmsh] 3D Periodic Boundary Conditions<br> </font> </div> <br>
<br>On 22 Mar 2013, at 11:33, Dorin Niculescu <<a ymailto="mailto:niculescu_dorin85@yahoo.com" href="mailto:niculescu_dorin85@yahoo.com">niculescu_dorin85@yahoo.com</a>> wrote:<br><br>> Hello Filipe,<br>> <br>> Thank you for your quick answer. Yes, i've tried to use Periodic Surface too but the result is even worst, the mesh seems to be the same on all the faces but when you are testing the coordinates you can see that there is a difference (not really big, but it can introduce errors in the solution). To find the corresponding mesh points i had to code a subroutine for testing numbers which allows me to set a tolerance. When i use Extrude and Transfinite Surface i find perfectly corresponding mesh just on the first direction in which i extrude (i've tried to extrude in the other two directions hoping that it will work, but with no result). Is this a bug (i'm using Gmsh v 2.5.1) or is like this that it is supposed to work?<br>>
<br><br>The periodic surface code has been rewritten in version 2.7: can you give it a try?<br><br><br>> Thank you very much for your time!<br>> <br>> Dorin<br>> <br>> From: F. A. Portela <<a ymailto="mailto:felipe.alves.portela@gmail.com" href="mailto:felipe.alves.portela@gmail.com">felipe.alves.portela@gmail.com</a>><br>> To: Dorin Niculescu <<a ymailto="mailto:niculescu_dorin85@yahoo.com" href="mailto:niculescu_dorin85@yahoo.com">niculescu_dorin85@yahoo.com</a>> <br>> Cc: "<a ymailto="mailto:gmsh@geuz.org" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>" <<a ymailto="mailto:gmsh@geuz.org" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>> <br>> Sent: Thursday, 21 March 2013, 23:34<br>> Subject: Re: [Gmsh] 3D Periodic Boundary Conditions<br>> <br>> Hi Dorin,<br>> <br>> You can add this lines at the end of your .geo to force the opposing faces to have the same mesh:<br>> <br>> Periodic
Surface 27 {4,13,11,22} = 19 {-2,14,-9,18};<br>> Periodic Surface 6 {2,3,4,1} = 28 {9,10,11,8};<br>> Periodic Surface 15 {14,8,13,1} = 23 {18,-10,22,-3};<br>> <br>> However it seems to me that by making the faces transfinite you should not require this...<br>> <br>> Cheers,<br>> Felipe<br>> <br>> <br>> --<br>> Felipe Alves Portela<br>> MSc student in Aerospace Eng. at TU Delft<br>> http://www.linkedin.com/in/felipealvesportela <br>> <br>> <br>> On 21 March 2013 11:28, Dorin Niculescu <<a ymailto="mailto:niculescu_dorin85@yahoo.com" href="mailto:niculescu_dorin85@yahoo.com">niculescu_dorin85@yahoo.com</a>> wrote:<br>> Hello all,<br>> <br>> I am trying to implement periodic boundary conditions for my 3D CFD code and i have troubles generating matching points on the opposite faces of a simple cube. I've used Extrude and Transfinite Surface as you can see below but the points seems
to match perfectly only for the zy plan. Is there a way to make them match on the other two plans too?<br>> <br>> cl = 20.0;<br>> Point(1) = {-50, 50, 50, cl};<br>> Point(2) = {-50, -50, 50, cl};<br>> Point(3) = {-50, -50, -50, cl};<br>> Point(4) = {-50, 50, -50, cl};<br>> Line(1) = {1, 2};<br>> Line(2) = {2, 3};<br>> Line(3) = {3, 4};<br>> Line(4) = {4, 1};<br>> Line Loop(5) = {1, 2, 3, 4};<br>> Plane Surface(6) = {5};<br>> Extrude {100, 0, 0} {<br>> Surface{6};<br>> }<br>> Extrude {0, 0, -100} {<br>> Surface{15};<br>> }<br>> Extrude {0, 100, 0} {<br>> Surface{19};<br>> }<br>> Transfinite Surface {6} = {1, 2, 3, 4};<br>> Transfinite Surface {28} = {5, 6, 10, 14};<br>> Transfinite Surface {15} = {1, 2, 6, 5};<br>> Transfinite Surface {23} = {4, 3, 10, 14};<br>> Transfinite Surface {19} = {2, 6, 10, 3};<br>> Transfinite Surface {27} = {1, 5, 14,
4};<br>> <br>> <br>> Thank you verry much for your time.<br>> <br>> Dorin Niculescu<br>> <br>> _______________________________________________<br>> gmsh mailing list<br>> <a ymailto="mailto:gmsh@geuz.org" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>> http://www.geuz.org/mailman/listinfo/gmsh<br>> <br>> <br>> <br>> <br>> _______________________________________________<br>> gmsh mailing list<br>> <a ymailto="mailto:gmsh@geuz.org" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>> <a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a><br><br>-- <br>Prof. Christophe Geuzaine<br>University of Liege, Electrical Engineering and Computer Science <br>http://www.montefiore.ulg.ac.be/~geuzaine<br><br><br><br><br><br> </div> </div> </div></body></html>