def METHOD_NAME(self):
    np.random.seed(10)
    spacing = 1e-1
    # 2d network
    self.net = op.network.Cubic(shape=[10, 10, 1], spacing=spacing)
    self.net.add_model_collection(
        op.models.collections.geometry.spheres_and_cylinders)
    self.net.regenerate_models()
    # phase and physics
    self.air = op.phase.Air(network=self.net)
    # add diffusive conductance model
    self.air['pore.temperature'] = 300
    self.air.add_model(propname='pore.diffusivity',
                       model=op.models.misc.linear,
                       m=1.860793056e-06,
                       b=-0.0005375624384,
                       prop='pore.temperature')
    self.air.add_model(propname='throat.diffusive_conductance',
                       model=pm.diffusive_conductance.generic_diffusive)
    # add thermal conductance model
    del self.air.models['pore.thermal_conductivity']
    self.air.add_model(propname='pore.thermal_conductivity',
                       model=op.models.misc.constant,
                       value=0.0262,
                       regen_mode='constant')
    self.air.add_model(propname='throat.thermal_conductance',
                       model=pm.thermal_conductance.generic_thermal)
    # settings
    self.tfd_settings = {
        "conductance": "throat.diffusive_conductance",
        "quantity": "pore.concentration",
        "cache_A": False,
        "cache_b": False
    }
    self.tfc_settings = {
        "conductance": "throat.thermal_conductance",
        "quantity": "pore.temperature",
        "pore_volume": "pore.heat_capacity",
        "cache_A": False,
        "cache_b": False
    }
    self.pardiso = op.solvers.PardisoSpsolve()
    self.rk45 = op.integrators.ScipyRK45(verbose=True)
    # First algorithm, transient fourier conduction
    self.tfc = op.algorithms.TransientReactiveTransport(network=self.net,
                                                        phase=self.air)
    self.net['pore.heat_capacity'] = self.net['pore.volume']*1229.2875
    self.tfc.settings._update(self.tfc_settings)
    self.tfc.set_value_BC(self.net.pores("left"), 400)
    self.T0 = np.ones(self.tfc.Np) * 300
    self.T0[self.net.pores('left')] = 400
    # Second algorithm, transient fickian diffusion
    self.tfd = op.algorithms.TransientReactiveTransport(network=self.net,
                                                        phase=self.air)
    self.tfd.settings._update(self.tfd_settings)
    self.tfd.set_value_BC(self.net.pores("left"), 100)
    self.c0 = np.ones(self.tfd.Np) * 50
    self.c0[self.net.pores('left')] = 100
    # add variable props to algs
    self.tfd.settings.variable_props.add('pore.temperature')
    self.tfc.settings.variable_props.add('pore.concentration')
    # multiphysics algorithm
    algs = [self.tfc, self.tfd]
    self.tmp = op.contrib.TransientMultiPhysics(algs, network=self.net)