On interior boundaries, the Surface Charge Accumulation node implements the boundary condition
where ρs is the solution of the following distributed ODE on the boundary:
where n · Ji is the normal component of the total ion current density on the wall and
n · Je is the normal component of the total electron current density on the wall.