I've simplified this to a case the works a little better (no crash) but still not there.

I've started with the gmsh tutorial file 1: t1.geo and modified it a bit

------------------------------------------
lc = 0.009;

Point(1) = {0, 0, 0, lc};
Point(2) = {.1, 0,  0, lc} ;
Point(3) = {.1, .3, 0, lc} ;
Point(4) = {0,  .3, 0, lc} ;

Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;

Line Loop(5) = {4,1,-2,3} ;

Plane Surface(6) = {5} ;

Physical Point(1) = {1,2} ;
Physical Point(2) = {3,4} ;

Physical Surface(7) = {6};
-------------------------------

I'd like to set the voltage along line 1 to 0V, and along line 3 to 150V say. I tried 'Physical Line', but that added '2_2' elements to the vtk file that sfepy didn't like. So.. I tried 'Physical Point' and that seemed to allow my code to run, but its still not setting up the region correctly. I used:

regions = {
    'Omega' : ('all', {}),
    'Omega_Source' : ('nodes of group 1', {}),
    'Omega_Target' : ('nodes of group 2', {}),
}


which ran "OK", but didn't set the voltages with:

ebcs = {
    'u1' : ('Omega_Source', {'u.0' : 0.0}),
    'u2' : ('Omega_Target', {'u.0' : 150.0}),
}

I'll keep looking. Thanks for any insight.

-steve