If we choose electric potential as variable, at the material interface the boundary conditions would be ,
v1 = v2
epsilon1 (dv1/dn) - epsilon2 (dv2/dn) = sigma
How will these boundary conditions change if the materials at the interface have different conductivity ?