Use the Internal Electrode Contact Resistance node as a subnode to the Electronic Conducting Phase to specify a Contact resistance at an interior boundary of the domain selection of the Electronic Conducting Phase. The contact resistance is modeled by adding slit boundary conditions. The slit boundary defines separate degrees of freedom for the electric potential at the up and down sides of the slit boundary, and then adds current density (Neumann) conditions at both the up and down side.