[Gmsh] 3D Periodic Boundary Conditions

[Dorin Niculescu]

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 zyplan. 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130321/d13e2384/attachment.html>