def METHOD_NAME(self):
    np.random.seed(0)
    self.net = op.network.Cubic(shape=[4, 3, 1], spacing=1.0)
    self.geo = op.geometry.GenericGeometry(network=self.net,
                                           pores=self.net.Ps,
                                           throats=self.net.Ts)
    self.geo['throat.conduit_lengths.pore1'] = 0.1
    self.geo['throat.conduit_lengths.throat'] = 0.6
    self.geo['throat.conduit_lengths.pore2'] = 0.1
    self.phase = op.phase.Phase(network=self.net)
    self.phase['pore.diffusivity.X'] = 1e-9
    self.phase['throat.valence.X'] = 1
    self.phys = op.physics.GenericPhysics(
        network=self.net, phase=self.phase, geometry=self.geo
    )
    self.phys['throat.diffusive_conductance.X'] = 1e-16
    self.phys['throat.hydraulic_conductance'] = 1e-15
    self.phys['throat.ionic_conductance'] = 1e-15
    self.sf = op.algorithms.StokesFlow(network=self.net, phase=self.phase)
    self.sf.set_value_BC(pores=self.net.pores('right'), values=1)
    self.sf.set_value_BC(pores=self.net.pores('left'), values=0)
    self.sf.run()
    self.phase.update(self.sf.results())
    self.oc = op.algorithms.OhmicConduction(network=self.net, phase=self.phase)
    self.oc.settings['conductance'] = 'throat.ionic_conductance'
    self.oc.settings['quantity'] = 'pore.potential'
    self.oc.set_value_BC(pores=self.net.pores('front'), values=0.05)
    self.oc.set_value_BC(pores=self.net.pores('back'), values=0.)
    self.oc.run()
    self.phase.update(self.oc.results())
    self.adm = op.algorithms.NernstPlanck(
        network=self.net, phase=self.phase, ion='X'
    )
    self.adm.settings._update({"cache": False})
    self.adm.set_value_BC(pores=self.net.pores('right'), values=2)
    self.adm.set_value_BC(pores=self.net.pores('left'), values=0)