Hello everyone,

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”

/*

LABELS

FRONT 1

BACK 2

LEFT 3

RIGHT 4

TOP 5

BOTTOM 6

INTERFACE 7

*/

mesh3 cellunit=readmesh3(“Unitcell.mesh”);

savemesh(cellunit,“total.mesh”);

plot(cellunit);

func perio=[[1,x,z],[3,x,z],[2,y,z],[4,y,z],[5,x,y],[6,x,y]];

real Ce0=1.73e4;

real k1=0.003;

real k2=0.002;

real n=0.2;

real D=1.41e-14;

real l=0.5;

real Free=1;

real Symmetry=2;

fespace Vh(cellunit,P1, periodic=perio);

fespace Ph(cellunit,P0);

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.

Best Regards

Jacob