Source code for pybamm.models.submodels.thermal.isothermal

#
# Class for isothermal submodel
#
import pybamm

from .base_thermal import BaseThermal




class Isothermal(BaseThermal):
    """
    Class for isothermal submodel.

    Parameters
    ----------
    param : parameter class
        The parameters to use for this submodel
    options : dict, optional
        A dictionary of options to be passed to the model.
    """

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



    def get_fundamental_variables(self):
        # Set the x-averaged temperature to the ambient temperature, which can be
        # specified as a function of space (y, z) only and time
        y = pybamm.standard_spatial_vars.y
        z = pybamm.standard_spatial_vars.z
        # Broadcast t to be the same size as y and z (to catch cases where the ambient
        # temperature is a function of time only)
        t_broadcast = pybamm.PrimaryBroadcast(pybamm.t, "current collector")
        T_x_av = self.param.T_amb(y, z, t_broadcast)
        T_vol_av = self._yz_average(T_x_av)

        T_dict = {
            "negative current collector": T_x_av,
            "positive current collector": T_x_av,
            "x-averaged cell": T_x_av,
            "volume-averaged cell": T_vol_av,
        }
        for domain in ["negative electrode", "separator", "positive electrode"]:
            T_dict[domain] = pybamm.PrimaryBroadcast(T_x_av, domain)

        variables = self._get_standard_fundamental_variables(T_dict)

        return variables




    def get_coupled_variables(self, variables):
        if self.options["calculate heat source for isothermal models"] == "true":
            variables.update(self._get_standard_coupled_variables(variables))
        else:
            zero = pybamm.Scalar(0)
            for var in [
                "Ohmic heating [W.m-3]",
                "X-averaged Ohmic heating [W.m-3]",
                "Volume-averaged Ohmic heating [W.m-3]",
                "Ohmic heating per unit electrode-pair area [W.m-2]",
                "Ohmic heating [W]",
                "Irreversible electrochemical heating [W.m-3]",
                "X-averaged irreversible electrochemical heating [W.m-3]",
                "Volume-averaged irreversible electrochemical heating [W.m-3]",
                "Irreversible electrochemical heating per unit electrode-pair area [W.m-2]",
                "Irreversible electrochemical heating [W]",
                "Reversible heating [W.m-3]",
                "X-averaged reversible heating [W.m-3]",
                "Volume-averaged reversible heating [W.m-3]",
                "Reversible heating per unit electrode-pair area [W.m-2]",
                "Reversible heating [W]",
                "Total heating [W.m-3]",
                "X-averaged total heating [W.m-3]",
                "Volume-averaged total heating [W.m-3]",
                "Total heating per unit electrode-pair area [W.m-2]",
                "Total heating [W]",
                "Negative current collector Ohmic heating [W.m-3]",
                "Positive current collector Ohmic heating [W.m-3]",
                "Surface total cooling [W.m-3]",
                "Surface total cooling [W]",
            ]:
                # All variables are zero
                variables.update({var: zero})

        return variables