def METHOD_NAME(self):
    for domain in ["negative", "positive"]:
        electrode_type = self.options.electrode_types[domain]
        if electrode_type == "porous":
            intercalation_kinetics = self.get_intercalation_kinetics(domain)
            phases = self.options.phases[domain]
            for phase in phases:
                submod = intercalation_kinetics(
                    self.param, domain, "lithium-ion main", self.options, phase
                )
                self.submodels[f"{domain} {phase} interface"] = submod
            if len(phases) > 1:
                self.submodels[
                    f"total {domain} interface"
                ] = pybamm.kinetics.TotalMainKinetics(
                    self.param, domain, "lithium-ion main", self.options
                )