Source code for pybamm.models.submodels.transport_efficiency.tortuosity_factor

#
# Class for tortuosity factor transport_efficiency
#
import pybamm

from .base_transport_efficiency import BaseModel




class TortuosityFactor(BaseModel):
    """Submodel for user supplied tortuosity factor transport_efficiency

    Parameters
    ----------
    param : parameter class
        The parameters to use for this submodel
    component : str
        The material for the model ('electrolyte' or 'electrode').
    options : dict, optional
        A dictionary of options to be passed to the model.
    """

    def __init__(self, param, component, options=None):
        super().__init__(param, component, options=options)



    def get_coupled_variables(self, variables):
        if self.component == "Electrolyte":
            tor_dict = {}
            for domain in self.options.whole_cell_domains:
                Domain = domain.capitalize()
                eps_k = variables[f"{Domain} porosity"]
                tau_k = self.param.domain_params[domain.split()[0]].tau_e
                tor_k = eps_k / tau_k
                tor_dict[domain] = tor_k
        elif self.component == "Electrode":
            tor_dict = {}
            for domain in self.options.whole_cell_domains:
                if domain == "separator":
                    tor_k = pybamm.FullBroadcast(0, "separator", "current collector")
                else:
                    Domain = domain.capitalize()
                    phi_k = 1 - variables[f"{Domain} porosity"]
                    tau_k = self.param.domain_params[domain.split()[0]].tau_s
                    tor_k = phi_k / tau_k
                tor_dict[domain] = tor_k
        variables.update(self._get_standard_transport_efficiency_variables(tor_dict))

        return variables