[Gmsh] 3D Periodic Boundary Conditions

[F. A. Portela]

Hi Dorin,

You can add this lines at the end of your .geo to force the opposing faces
to have the same mesh:

  Periodic Surface 27 {4,13,11,22} = 19 {-2,14,-9,18};
  Periodic Surface 6 {2,3,4,1} = 28 {9,10,11,8};
  Periodic Surface 15 {14,8,13,1} = 23 {18,-10,22,-3};

However it seems to me that by making the faces transfinite you should not
require this...

Cheers,
Felipe


--
Felipe Alves Portela
MSc student in Aerospace Eng. at TU Delft
http://www.linkedin.com/in/felipealvesportela


On 21 March 2013 11:28, Dorin Niculescu <niculescu_dorin85 at yahoo.com> wrote:

> Hello all,
>
> 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?
>
> cl = 20.0;
> Point(1) = {-50, 50, 50, cl};
> Point(2) = {-50, -50, 50, cl};
> Point(3) = {-50, -50, -50, cl};
> Point(4) = {-50, 50, -50, cl};
> Line(1) = {1, 2};
> Line(2) = {2, 3};
> Line(3) = {3, 4};
> Line(4) = {4, 1};
> Line Loop(5) = {1, 2, 3, 4};
> Plane Surface(6) = {5};
> Extrude {100, 0, 0} {
>   Surface{6};
> }
> Extrude {0, 0, -100} {
>   Surface{15};
> }
> Extrude {0, 100, 0} {
>   Surface{19};
> }
> Transfinite Surface {6} = {1, 2, 3, 4};
> Transfinite Surface {28} = {5, 6, 10, 14};
> Transfinite Surface {15} = {1, 2, 6, 5};
> Transfinite Surface {23} = {4, 3, 10, 14};
> Transfinite Surface {19} = {2, 6, 10, 3};
> Transfinite Surface {27} = {1, 5, 14, 4};
>
>
> Thank you verry much for your time.
>
> Dorin Niculescu
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130321/3c7b5c5c/attachment.html>