I’m trying to define periodic boundary conditions on a cubic domain decomposed in two regions. In the paralell faces(labels [1,2], [3,4], [5,6]) were imposed identical boudary conditions. In the interface(label 7) was assumed as free.
load “msh3” load “iovtk” load “gmsh” load “MUMPS_seq”
fespace Vh(cellunit,P1, periodic=perio);
However when I execute the problem, I get the following error:
Exec error : Periodic 3d: the both number of face is not the same.
Below you have attached the .edp and the .mesh file.
Test.edp (534 Bytes)
Periodic 3d: the both number of face is not the same
Any suggestion will be very appreciated.