A numerical method for renal models that represent tubules with abrupt changes in membrane properties