{ "package": "batt", "modules": { "pathsim_batt.cells.pybamm_cell": { "name": "pathsim_batt.cells.pybamm_cell", "description": "", "docstring_html": "", "classes": [ { "name": "_CellBase", "description": "Shared base for PyBaMM cell blocks.", "docstring_html": "

Shared base for PyBaMM cell blocks.

Discretises the PyBaMM model at construction time and exposes its ODE\nright-hand side to PathSim's numerical integrator via the DynamicalSystem\ninterface. The differential state vector of the discretised model becomes\nthe PathSim state; PathSim's chosen solver advances it in time.

Only PyBaMM models that produce a pure ODE after discretisation are\nsupported (i.e. models with no algebraic variables, such as SPMe and SPM).\nModels that result in a DAE system (e.g. DFN) are not supported and will\nraise NotImplementedError at construction time.

Because the SPMe/SPM family of models is stiff, users should prefer an\nimplicit solver (e.g. ESDIRK43, BDF) when constructing the\nPathSim Simulation.

Subclasses set _thermal_option and _pybamm_output_vars to select the\nthermal sub-model and define which PyBaMM variables map to the block's\noutput ports (SOC is always appended last). _pybamm_output_vars must\ncontain _TERMINAL_VOLTAGE_VAR; its position in the list is found\ndynamically.

\n", "source": "class _CellBase(DynamicalSystem):\n \"\"\"Shared base for PyBaMM cell blocks.\n\n Discretises the PyBaMM model at construction time and exposes its ODE\n right-hand side to PathSim's numerical integrator via the ``DynamicalSystem``\n interface. The differential state vector of the discretised model becomes\n the PathSim state; PathSim's chosen solver advances it in time.\n\n Only PyBaMM models that produce a pure ODE after discretisation are\n supported (i.e. models with no algebraic variables, such as SPMe and SPM).\n Models that result in a DAE system (e.g. DFN) are not supported and will\n raise ``NotImplementedError`` at construction time.\n\n Because the SPMe/SPM family of models is stiff, users should prefer an\n implicit solver (e.g. ``ESDIRK43``, ``BDF``) when constructing the\n PathSim ``Simulation``.\n\n Subclasses set ``_thermal_option`` and ``_pybamm_output_vars`` to select the\n thermal sub-model and define which PyBaMM variables map to the block's\n output ports (SOC is always appended last). ``_pybamm_output_vars`` must\n contain ``_TERMINAL_VOLTAGE_VAR``; its position in the list is found\n dynamically.\n \"\"\"\n\n _thermal_option: str = \"\"\n _thermal_extra_options: dict[str, str] = {}\n _pybamm_output_vars: list[str] = []\n\n def __init__(\n self,\n model: pybamm.BaseBatteryModel | None = None,\n parameter_values: pybamm.ParameterValues | None = None,\n initial_soc: float = 1.0,\n pybamm_solver: pybamm.BaseSolver | None = None,\n ) -> None:\n self._initial_soc = float(initial_soc)\n\n try:\n self._v_idx = self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)\n except ValueError:\n raise TypeError(\n f\"{type(self).__name__}._pybamm_output_vars must contain \"\n f\"'{_TERMINAL_VOLTAGE_VAR}'.\"\n ) from None\n\n if model is None:\n model = pybamm.lithium_ion.SPMe(\n options={\"thermal\": self._thermal_option, **self._thermal_extra_options}\n )\n\n self._parameter_values = _prepare_parameter_values(parameter_values)\n try:\n self._v_lower = float(self._parameter_values[\"Lower voltage cut-off [V]\"])\n self._v_upper = float(self._parameter_values[\"Upper voltage cut-off [V]\"])\n except KeyError as exc:\n raise ValueError(\n f\"parameter_values is missing a voltage cut-off entry: {exc}. \"\n \"Ensure your parameter set defines both 'Lower voltage cut-off [V]' \"\n \"and 'Upper voltage cut-off [V]'.\"\n ) from exc\n\n pybamm_solver = pybamm_solver or pybamm.CasadiSolver(mode=\"safe\")\n\n sim = pybamm.Simulation(\n model,\n parameter_values=self._parameter_values,\n solver=pybamm_solver,\n )\n sim.build(initial_soc=self._initial_soc, inputs=_DEFAULT_INPUTS)\n\n all_out_vars = self._pybamm_output_vars + [\"Discharge capacity [A.h]\"]\n casadi_objs = sim.built_model.export_casadi_objects(\n all_out_vars,\n input_parameter_order=list(_DEFAULT_INPUTS.keys()),\n )\n\n t_sym = casadi_objs[\"t\"]\n x_sym = casadi_objs[\"x\"]\n z_sym = casadi_objs[\"z\"]\n p_sym = casadi_objs[\"inputs\"]\n\n if z_sym.numel() > 0:\n raise NotImplementedError(\n f\"{type(self).__name__}: the supplied PyBaMM model has \"\n f\"{z_sym.numel()} algebraic variable(s) after discretisation \"\n \"(DAE system). Only pure ODE models are supported. \"\n \"Use SPMe or SPM instead of DFN.\"\n )\n\n rhs_fn = casadi.Function(\"rhs\", [t_sym, x_sym, p_sym], [casadi_objs[\"rhs\"]])\n jac_fn = casadi.Function(\n \"jac_rhs\", [t_sym, x_sym, p_sym], [casadi_objs[\"jac_rhs\"]]\n )\n\n out_var_fns = {}\n for idx, var_name in enumerate(all_out_vars):\n var_expr = casadi_objs[\"variables\"][var_name]\n out_var_fns[var_name] = casadi.Function(\n f\"outvar_{idx}\", [t_sym, x_sym, p_sym], [var_expr]\n )\n\n self._casadi_rhs = rhs_fn\n self._jac_rhs_eval = jac_fn\n self._out_var_fcns = out_var_fns\n self._q_nominal = float(self._parameter_values[\"Nominal cell capacity [A.h]\"])\n\n q_nominal = self._q_nominal\n initial_soc_val = float(initial_soc)\n pybamm_output_vars = list(self._pybamm_output_vars)\n\n def _pack(u):\n return casadi.DM([float(u[0]), float(u[1])])\n\n def func_dyn(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n return np.array(rhs_fn(t, xv, p)).flatten()\n\n def jac_dyn(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n return np.array(jac_fn(t, xv, p))\n\n v_lower = self._v_lower\n v_upper = self._v_upper\n v_idx = self._v_idx\n\n def func_alg(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n outputs = [float(out_var_fns[n](t, xv, p)) for n in pybamm_output_vars]\n q_dis = float(out_var_fns[\"Discharge capacity [A.h]\"](t, xv, p))\n soc = max(0.0, min(1.0, initial_soc_val - q_dis / q_nominal))\n outputs.append(soc)\n V = outputs[v_idx]\n if V <= v_lower:\n raise StopSimulation(f\"undervoltage: V={V:.4f} V <= {v_lower} V\")\n if V >= v_upper:\n raise StopSimulation(f\"overvoltage: V={V:.4f} V >= {v_upper} V\")\n return np.array(outputs)\n\n x0_fn = casadi.Function(\"x0\", [p_sym], [casadi_objs[\"x0\"]])\n\n y0 = np.array(x0_fn(casadi.DM(list(_DEFAULT_INPUTS.values())))).flatten()\n\n super().__init__(\n func_dyn=func_dyn,\n func_alg=func_alg,\n initial_value=y0,\n jac_dyn=jac_dyn,\n )\n\n def __len__(self) -> int:\n return len(self._pybamm_output_vars) + 1\n\n def reset(self) -> None:\n super().reset()", "bases": [ "pathsim.blocks.DynamicalSystem" ], "methods": [ { "name": "__init__", "description": "", "docstring_html": "", "source": "def __init__(\n self,\n model: pybamm.BaseBatteryModel | None = None,\n parameter_values: pybamm.ParameterValues | None = None,\n initial_soc: float = 1.0,\n pybamm_solver: pybamm.BaseSolver | None = None,\n) -> None:\n self._initial_soc = float(initial_soc)\n\n try:\n self._v_idx = self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)\n except ValueError:\n raise TypeError(\n f\"{type(self).__name__}._pybamm_output_vars must contain \"\n f\"'{_TERMINAL_VOLTAGE_VAR}'.\"\n ) from None\n\n if model is None:\n model = pybamm.lithium_ion.SPMe(\n options={\"thermal\": self._thermal_option, **self._thermal_extra_options}\n )\n\n self._parameter_values = _prepare_parameter_values(parameter_values)\n try:\n self._v_lower = float(self._parameter_values[\"Lower voltage cut-off [V]\"])\n self._v_upper = float(self._parameter_values[\"Upper voltage cut-off [V]\"])\n except KeyError as exc:\n raise ValueError(\n f\"parameter_values is missing a voltage cut-off entry: {exc}. \"\n \"Ensure your parameter set defines both 'Lower voltage cut-off [V]' \"\n \"and 'Upper voltage cut-off [V]'.\"\n ) from exc\n\n pybamm_solver = pybamm_solver or pybamm.CasadiSolver(mode=\"safe\")\n\n sim = pybamm.Simulation(\n model,\n parameter_values=self._parameter_values,\n solver=pybamm_solver,\n )\n sim.build(initial_soc=self._initial_soc, inputs=_DEFAULT_INPUTS)\n\n all_out_vars = self._pybamm_output_vars + [\"Discharge capacity [A.h]\"]\n casadi_objs = sim.built_model.export_casadi_objects(\n all_out_vars,\n input_parameter_order=list(_DEFAULT_INPUTS.keys()),\n )\n\n t_sym = casadi_objs[\"t\"]\n x_sym = casadi_objs[\"x\"]\n z_sym = casadi_objs[\"z\"]\n p_sym = casadi_objs[\"inputs\"]\n\n if z_sym.numel() > 0:\n raise NotImplementedError(\n f\"{type(self).__name__}: the supplied PyBaMM model has \"\n f\"{z_sym.numel()} algebraic variable(s) after discretisation \"\n \"(DAE system). Only pure ODE models are supported. \"\n \"Use SPMe or SPM instead of DFN.\"\n )\n\n rhs_fn = casadi.Function(\"rhs\", [t_sym, x_sym, p_sym], [casadi_objs[\"rhs\"]])\n jac_fn = casadi.Function(\n \"jac_rhs\", [t_sym, x_sym, p_sym], [casadi_objs[\"jac_rhs\"]]\n )\n\n out_var_fns = {}\n for idx, var_name in enumerate(all_out_vars):\n var_expr = casadi_objs[\"variables\"][var_name]\n out_var_fns[var_name] = casadi.Function(\n f\"outvar_{idx}\", [t_sym, x_sym, p_sym], [var_expr]\n )\n\n self._casadi_rhs = rhs_fn\n self._jac_rhs_eval = jac_fn\n self._out_var_fcns = out_var_fns\n self._q_nominal = float(self._parameter_values[\"Nominal cell capacity [A.h]\"])\n\n q_nominal = self._q_nominal\n initial_soc_val = float(initial_soc)\n pybamm_output_vars = list(self._pybamm_output_vars)\n\n def _pack(u):\n return casadi.DM([float(u[0]), float(u[1])])\n\n def func_dyn(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n return np.array(rhs_fn(t, xv, p)).flatten()\n\n def jac_dyn(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n return np.array(jac_fn(t, xv, p))\n\n v_lower = self._v_lower\n v_upper = self._v_upper\n v_idx = self._v_idx\n\n def func_alg(x, u, t):\n xv = casadi.DM(x.reshape(-1, 1))\n p = _pack(u)\n outputs = [float(out_var_fns[n](t, xv, p)) for n in pybamm_output_vars]\n q_dis = float(out_var_fns[\"Discharge capacity [A.h]\"](t, xv, p))\n soc = max(0.0, min(1.0, initial_soc_val - q_dis / q_nominal))\n outputs.append(soc)\n V = outputs[v_idx]\n if V <= v_lower:\n raise StopSimulation(f\"undervoltage: V={V:.4f} V <= {v_lower} V\")\n if V >= v_upper:\n raise StopSimulation(f\"overvoltage: V={V:.4f} V >= {v_upper} V\")\n return np.array(outputs)\n\n x0_fn = casadi.Function(\"x0\", [p_sym], [casadi_objs[\"x0\"]])\n\n y0 = np.array(x0_fn(casadi.DM(list(_DEFAULT_INPUTS.values())))).flatten()\n\n super().__init__(\n func_dyn=func_dyn,\n func_alg=func_alg,\n initial_value=y0,\n jac_dyn=jac_dyn,\n )", "signature": "(model: pybamm.BaseBatteryModel | None = None, parameter_values: pybamm.ParameterValues | None = None, initial_soc: float = 1.0, pybamm_solver: pybamm.BaseSolver | None = None)", "parameters": [ { "name": "model", "type": "pybamm.BaseBatteryModel | None", "default": "None", "description": "" }, { "name": "parameter_values", "type": "pybamm.ParameterValues | None", "default": "None", "description": "" }, { "name": "initial_soc", "type": "float", "default": "1.0", "description": "" }, { "name": "pybamm_solver", "type": "pybamm.BaseSolver | None", "default": "None", "description": "" } ], "returns": "None", "method_type": "method" }, { "name": "reset", "description": "", "docstring_html": "", "source": "def reset(self) -> None:\n super().reset()", "signature": "()", "parameters": [], "returns": "None", "method_type": "method" } ], "attributes": [ { "name": "_thermal_option", "description": "", "type": "str", "value": "''" }, { "name": "_thermal_extra_options", "description": "", "type": "dict[str, str]", "value": "{}" }, { "name": "_pybamm_output_vars", "description": "", "type": "list[str]", "value": "[]" }, { "name": "_initial_soc", "description": "", "type": null, "value": "float(initial_soc)" }, { "name": "_v_idx", "description": "", "type": null, "value": "self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)" }, { "name": "_parameter_values", "description": "", "type": null, "value": "_prepare_parameter_values(parameter_values)" }, { "name": "_v_lower", "description": "", "type": null, "value": "float(self._parameter_values['Lower voltage cut-off [V]'])" }, { "name": "_v_upper", "description": "", "type": null, "value": "float(self._parameter_values['Upper voltage cut-off [V]'])" }, { "name": "_casadi_rhs", "description": "", "type": null, "value": "rhs_fn" }, { "name": "_jac_rhs_eval", "description": "", "type": null, "value": "jac_fn" }, { "name": "_out_var_fcns", "description": "", "type": null, "value": "out_var_fns" }, { "name": "_q_nominal", "description": "", "type": null, "value": "float(self._parameter_values['Nominal cell capacity [A.h]'])" } ], "parameters": [ { "name": "model", "type": "pybamm.BaseBatteryModel | None", "default": "None", "description": "" }, { "name": "parameter_values", "type": "pybamm.ParameterValues | None", "default": "None", "description": "" }, { "name": "initial_soc", "type": "float", "default": "1.0", "description": "" }, { "name": "pybamm_solver", "type": "pybamm.BaseSolver | None", "default": "None", "description": "" } ] }, { "name": "_CoSimCellBase", "description": "Shared base for co-simulation PyBaMM cell blocks.", "docstring_html": "

Shared base for co-simulation PyBaMM cell blocks.

Wraps pybamm.Simulation.step() in a periodic Wrapper event, so\nPyBaMM advances on discrete macro-steps while PathSim sees a zero-order-held\noutput signal between events. This allows using PyBaMM models that produce\nDAE systems after discretisation (e.g. DFN), because PyBaMM owns the\ndifferential-algebraic solve internally.

Subclasses set _thermal_option, _pybamm_output_vars and port labels.\n_pybamm_output_vars must contain _TERMINAL_VOLTAGE_VAR; its position\nin the list is found dynamically.

\n", "source": "class _CoSimCellBase(Wrapper):\n \"\"\"Shared base for co-simulation PyBaMM cell blocks.\n\n Wraps ``pybamm.Simulation.step()`` in a periodic ``Wrapper`` event, so\n PyBaMM advances on discrete macro-steps while PathSim sees a zero-order-held\n output signal between events. This allows using PyBaMM models that produce\n DAE systems after discretisation (e.g. DFN), because PyBaMM owns the\n differential-algebraic solve internally.\n\n Subclasses set ``_thermal_option``, ``_pybamm_output_vars`` and port labels.\n ``_pybamm_output_vars`` must contain ``_TERMINAL_VOLTAGE_VAR``; its position\n in the list is found dynamically.\n \"\"\"\n\n _thermal_option: str = \"\"\n _thermal_extra_options: dict[str, str] = {}\n _pybamm_output_vars: list[str] = []\n\n def __init__(\n self,\n model: pybamm.BaseBatteryModel | None = None,\n parameter_values: pybamm.ParameterValues | None = None,\n initial_soc: float = 1.0,\n pybamm_solver: pybamm.BaseSolver | None = None,\n dt: float = 1.0,\n ) -> None:\n self._initial_soc = float(initial_soc)\n\n try:\n self._v_idx = self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)\n except ValueError:\n raise TypeError(\n f\"{type(self).__name__}._pybamm_output_vars must contain \"\n f\"'{_TERMINAL_VOLTAGE_VAR}'.\"\n ) from None\n self._dt = float(dt)\n if self._dt <= 0.0:\n raise ValueError(\"dt must be positive\")\n\n if model is None:\n model = pybamm.lithium_ion.SPMe(\n options={\"thermal\": self._thermal_option, **self._thermal_extra_options}\n )\n\n self._model = model\n self._parameter_values = _prepare_parameter_values(parameter_values)\n try:\n self._v_lower = float(self._parameter_values[\"Lower voltage cut-off [V]\"])\n self._v_upper = float(self._parameter_values[\"Upper voltage cut-off [V]\"])\n except KeyError as exc:\n raise ValueError(\n f\"parameter_values is missing a voltage cut-off entry: {exc}. \"\n \"Ensure your parameter set defines both 'Lower voltage cut-off [V]' \"\n \"and 'Upper voltage cut-off [V]'.\"\n ) from exc\n self._pybamm_solver = pybamm_solver or pybamm.IDAKLUSolver()\n self._q_nominal = float(self._parameter_values[\"Nominal cell capacity [A.h]\"])\n\n self._sim = self._build_sim()\n\n self._last_outputs: npt.NDArray[np.float64] = self._initial_outputs()\n\n super().__init__(func=self._discrete_step, T=self._dt, tau=self._dt)\n\n # ensure outputs are valid before first scheduled sample\n self.outputs.update_from_array(self._last_outputs)\n\n def _build_sim(self) -> pybamm.Simulation:\n \"\"\"Create and build a fresh ``pybamm.Simulation`` with default inputs.\"\"\"\n sim = pybamm.Simulation(\n self._model,\n parameter_values=self._parameter_values,\n solver=self._pybamm_solver,\n )\n sim.build(initial_soc=self._initial_soc, inputs=_DEFAULT_INPUTS)\n return sim\n\n def _initial_outputs(self) -> npt.NDArray[np.float64]:\n \"\"\"Compute outputs at t=0 from the built PyBaMM model using default inputs.\n\n Uses the same CasADi export approach as ``_CellBase`` to evaluate each\n output variable at the initial state vector. The evaluation uses\n ``_DEFAULT_INPUTS`` (0 A current, 298.15 K ambient temperature) because\n the wired input ports are not yet available at construction time. The\n resulting open-circuit voltage is therefore physically meaningful but does\n not account for a non-default initial temperature or a non-zero current at\n t=0.\n \"\"\"\n all_out_vars = self._pybamm_output_vars + [\"Discharge capacity [A.h]\"]\n casadi_objs = self._sim.built_model.export_casadi_objects(\n all_out_vars,\n input_parameter_order=list(_DEFAULT_INPUTS.keys()),\n )\n t_sym = casadi_objs[\"t\"]\n x_sym = casadi_objs[\"x\"]\n z_sym = casadi_objs[\"z\"]\n p_sym = casadi_objs[\"inputs\"]\n p0 = casadi.DM(list(_DEFAULT_INPUTS.values()))\n x0 = casadi.Function(\"x0\", [p_sym], [casadi_objs[\"x0\"]])(p0)\n # Algebraic initial conditions (empty for ODE models such as SPMe;\n # non-empty for DAE models such as DFN).\n z0 = casadi.Function(\"z0\", [p_sym], [casadi_objs[\"z0\"]])(p0)\n\n outputs: list[float] = []\n for name in self._pybamm_output_vars:\n fn = casadi.Function(\n \"v\", [t_sym, x_sym, z_sym, p_sym], [casadi_objs[\"variables\"][name]]\n )\n outputs.append(float(fn(0.0, x0, z0, p0)))\n\n q_dis_fn = casadi.Function(\n \"q\",\n [t_sym, x_sym, z_sym, p_sym],\n [casadi_objs[\"variables\"][\"Discharge capacity [A.h]\"]],\n )\n q_dis = float(q_dis_fn(0.0, x0, z0, p0))\n soc = max(0.0, min(1.0, self._initial_soc - q_dis / self._q_nominal))\n outputs.append(soc)\n\n return np.array(outputs, dtype=np.float64)\n\n def _discrete_step(self, current: float, t_amb: float) -> npt.NDArray[np.float64]:\n inputs = {\n \"Current function [A]\": float(current),\n \"Ambient temperature [K]\": float(t_amb),\n }\n self._sim.step(dt=self._dt, inputs=inputs, save=False)\n\n sol = self._sim.solution\n outputs = [float(sol[n].entries[-1]) for n in self._pybamm_output_vars]\n q_dis = float(sol[\"Discharge capacity [A.h]\"].entries[-1])\n soc = max(0.0, min(1.0, self._initial_soc - q_dis / self._q_nominal))\n outputs.append(soc)\n\n self._last_outputs = np.array(outputs, dtype=np.float64)\n self.outputs.update_from_array(self._last_outputs)\n V = outputs[self._v_idx]\n if V <= self._v_lower:\n raise StopSimulation(f\"undervoltage: V={V:.4f} V <= {self._v_lower} V\")\n if V >= self._v_upper:\n raise StopSimulation(f\"overvoltage: V={V:.4f} V >= {self._v_upper} V\")\n return self._last_outputs\n\n def __len__(self) -> int:\n return len(self._pybamm_output_vars) + 1\n\n def reset(self) -> None:\n super().reset()\n self._sim = self._build_sim()\n self._last_outputs = self._initial_outputs()\n self.outputs.update_from_array(self._last_outputs)", "bases": [ "pathsim.blocks.Wrapper" ], "methods": [ { "name": "__init__", "description": "", "docstring_html": "", "source": "def __init__(\n self,\n model: pybamm.BaseBatteryModel | None = None,\n parameter_values: pybamm.ParameterValues | None = None,\n initial_soc: float = 1.0,\n pybamm_solver: pybamm.BaseSolver | None = None,\n dt: float = 1.0,\n) -> None:\n self._initial_soc = float(initial_soc)\n\n try:\n self._v_idx = self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)\n except ValueError:\n raise TypeError(\n f\"{type(self).__name__}._pybamm_output_vars must contain \"\n f\"'{_TERMINAL_VOLTAGE_VAR}'.\"\n ) from None\n self._dt = float(dt)\n if self._dt <= 0.0:\n raise ValueError(\"dt must be positive\")\n\n if model is None:\n model = pybamm.lithium_ion.SPMe(\n options={\"thermal\": self._thermal_option, **self._thermal_extra_options}\n )\n\n self._model = model\n self._parameter_values = _prepare_parameter_values(parameter_values)\n try:\n self._v_lower = float(self._parameter_values[\"Lower voltage cut-off [V]\"])\n self._v_upper = float(self._parameter_values[\"Upper voltage cut-off [V]\"])\n except KeyError as exc:\n raise ValueError(\n f\"parameter_values is missing a voltage cut-off entry: {exc}. \"\n \"Ensure your parameter set defines both 'Lower voltage cut-off [V]' \"\n \"and 'Upper voltage cut-off [V]'.\"\n ) from exc\n self._pybamm_solver = pybamm_solver or pybamm.IDAKLUSolver()\n self._q_nominal = float(self._parameter_values[\"Nominal cell capacity [A.h]\"])\n\n self._sim = self._build_sim()\n\n self._last_outputs: npt.NDArray[np.float64] = self._initial_outputs()\n\n super().__init__(func=self._discrete_step, T=self._dt, tau=self._dt)\n\n # ensure outputs are valid before first scheduled sample\n self.outputs.update_from_array(self._last_outputs)", "signature": "(model: pybamm.BaseBatteryModel | None = None, parameter_values: pybamm.ParameterValues | None = None, initial_soc: float = 1.0, pybamm_solver: pybamm.BaseSolver | None = None, dt: float = 1.0)", "parameters": [ { "name": "model", "type": "pybamm.BaseBatteryModel | None", "default": "None", "description": "" }, { "name": "parameter_values", "type": "pybamm.ParameterValues | None", "default": "None", "description": "" }, { "name": "initial_soc", "type": "float", "default": "1.0", "description": "" }, { "name": "pybamm_solver", "type": "pybamm.BaseSolver | None", "default": "None", "description": "" }, { "name": "dt", "type": "float", "default": "1.0", "description": "" } ], "returns": "None", "method_type": "method" }, { "name": "_build_sim", "description": "Create and build a fresh ``pybamm.Simulation`` with default inputs.", "docstring_html": "

Create and build a fresh pybamm.Simulation with default inputs.

\n", "source": "def _build_sim(self) -> pybamm.Simulation:\n \"\"\"Create and build a fresh ``pybamm.Simulation`` with default inputs.\"\"\"\n sim = pybamm.Simulation(\n self._model,\n parameter_values=self._parameter_values,\n solver=self._pybamm_solver,\n )\n sim.build(initial_soc=self._initial_soc, inputs=_DEFAULT_INPUTS)\n return sim", "signature": "()", "parameters": [], "returns": "pybamm.Simulation", "method_type": "method" }, { "name": "_initial_outputs", "description": "Compute outputs at t=0 from the built PyBaMM model using default inputs.", "docstring_html": "

Compute outputs at t=0 from the built PyBaMM model using default inputs.

Uses the same CasADi export approach as _CellBase to evaluate each\noutput variable at the initial state vector. The evaluation uses\n_DEFAULT_INPUTS (0 A current, 298.15 K ambient temperature) because\nthe wired input ports are not yet available at construction time. The\nresulting open-circuit voltage is therefore physically meaningful but does\nnot account for a non-default initial temperature or a non-zero current at\nt=0.

\n", "source": "def _initial_outputs(self) -> npt.NDArray[np.float64]:\n \"\"\"Compute outputs at t=0 from the built PyBaMM model using default inputs.\n\n Uses the same CasADi export approach as ``_CellBase`` to evaluate each\n output variable at the initial state vector. The evaluation uses\n ``_DEFAULT_INPUTS`` (0 A current, 298.15 K ambient temperature) because\n the wired input ports are not yet available at construction time. The\n resulting open-circuit voltage is therefore physically meaningful but does\n not account for a non-default initial temperature or a non-zero current at\n t=0.\n \"\"\"\n all_out_vars = self._pybamm_output_vars + [\"Discharge capacity [A.h]\"]\n casadi_objs = self._sim.built_model.export_casadi_objects(\n all_out_vars,\n input_parameter_order=list(_DEFAULT_INPUTS.keys()),\n )\n t_sym = casadi_objs[\"t\"]\n x_sym = casadi_objs[\"x\"]\n z_sym = casadi_objs[\"z\"]\n p_sym = casadi_objs[\"inputs\"]\n p0 = casadi.DM(list(_DEFAULT_INPUTS.values()))\n x0 = casadi.Function(\"x0\", [p_sym], [casadi_objs[\"x0\"]])(p0)\n # Algebraic initial conditions (empty for ODE models such as SPMe;\n # non-empty for DAE models such as DFN).\n z0 = casadi.Function(\"z0\", [p_sym], [casadi_objs[\"z0\"]])(p0)\n\n outputs: list[float] = []\n for name in self._pybamm_output_vars:\n fn = casadi.Function(\n \"v\", [t_sym, x_sym, z_sym, p_sym], [casadi_objs[\"variables\"][name]]\n )\n outputs.append(float(fn(0.0, x0, z0, p0)))\n\n q_dis_fn = casadi.Function(\n \"q\",\n [t_sym, x_sym, z_sym, p_sym],\n [casadi_objs[\"variables\"][\"Discharge capacity [A.h]\"]],\n )\n q_dis = float(q_dis_fn(0.0, x0, z0, p0))\n soc = max(0.0, min(1.0, self._initial_soc - q_dis / self._q_nominal))\n outputs.append(soc)\n\n return np.array(outputs, dtype=np.float64)", "signature": "()", "parameters": [], "returns": "npt.NDArray[np.float64]", "method_type": "method" }, { "name": "_discrete_step", "description": "", "docstring_html": "", "source": "def _discrete_step(self, current: float, t_amb: float) -> npt.NDArray[np.float64]:\n inputs = {\n \"Current function [A]\": float(current),\n \"Ambient temperature [K]\": float(t_amb),\n }\n self._sim.step(dt=self._dt, inputs=inputs, save=False)\n\n sol = self._sim.solution\n outputs = [float(sol[n].entries[-1]) for n in self._pybamm_output_vars]\n q_dis = float(sol[\"Discharge capacity [A.h]\"].entries[-1])\n soc = max(0.0, min(1.0, self._initial_soc - q_dis / self._q_nominal))\n outputs.append(soc)\n\n self._last_outputs = np.array(outputs, dtype=np.float64)\n self.outputs.update_from_array(self._last_outputs)\n V = outputs[self._v_idx]\n if V <= self._v_lower:\n raise StopSimulation(f\"undervoltage: V={V:.4f} V <= {self._v_lower} V\")\n if V >= self._v_upper:\n raise StopSimulation(f\"overvoltage: V={V:.4f} V >= {self._v_upper} V\")\n return self._last_outputs", "signature": "(current: float, t_amb: float)", "parameters": [ { "name": "current", "type": "float", "default": null, "description": "" }, { "name": "t_amb", "type": "float", "default": null, "description": "" } ], "returns": "npt.NDArray[np.float64]", "method_type": "method" }, { "name": "reset", "description": "", "docstring_html": "", "source": "def reset(self) -> None:\n super().reset()\n self._sim = self._build_sim()\n self._last_outputs = self._initial_outputs()\n self.outputs.update_from_array(self._last_outputs)", "signature": "()", "parameters": [], "returns": "None", "method_type": "method" } ], "attributes": [ { "name": "_thermal_option", "description": "", "type": "str", "value": "''" }, { "name": "_thermal_extra_options", "description": "", "type": "dict[str, str]", "value": "{}" }, { "name": "_pybamm_output_vars", "description": "", "type": "list[str]", "value": "[]" }, { "name": "_initial_soc", "description": "", "type": null, "value": "float(initial_soc)" }, { "name": "_v_idx", "description": "", "type": null, "value": "self._pybamm_output_vars.index(_TERMINAL_VOLTAGE_VAR)" }, { "name": "_dt", "description": "", "type": null, "value": "float(dt)" }, { "name": "_model", "description": "", "type": null, "value": "model" }, { "name": "_parameter_values", "description": "", "type": null, "value": "_prepare_parameter_values(parameter_values)" }, { "name": "_v_lower", "description": "", "type": null, "value": "float(self._parameter_values['Lower voltage cut-off [V]'])" }, { "name": "_v_upper", "description": "", "type": null, "value": "float(self._parameter_values['Upper voltage cut-off [V]'])" }, { "name": "_pybamm_solver", "description": "", "type": null, "value": "pybamm_solver or pybamm.IDAKLUSolver()" }, { "name": "_q_nominal", "description": "", "type": null, "value": "float(self._parameter_values['Nominal cell capacity [A.h]'])" }, { "name": "_sim", "description": "", "type": null, "value": "self._build_sim()" }, { "name": "_last_outputs", "description": "", "type": "npt.NDArray[np.float64]", "value": "self._initial_outputs()" } ], "parameters": [ { "name": "model", "type": "pybamm.BaseBatteryModel | None", "default": "None", "description": "" }, { "name": "parameter_values", "type": "pybamm.ParameterValues | None", "default": "None", "description": "" }, { "name": "initial_soc", "type": "float", "default": "1.0", "description": "" }, { "name": "pybamm_solver", "type": "pybamm.BaseSolver | None", "default": "None", "description": "" }, { "name": "dt", "type": "float", "default": "1.0", "description": "" } ] }, { "name": "CellElectrical", "description": "Cell block — electrical outputs only, external thermal coupling.", "docstring_html": "

Cell block — electrical outputs only, external thermal coupling.

PathSim integrates the electrochemical state via the discretised PyBaMM\nODE. Temperature dynamics live outside this block: wire Q_dot to a\nLumpedThermal (or similar) block and feed its temperature output back\nto T_cell.

Note

The SPMe/SPM ODE is stiff. Use an implicit solver (e.g.\nESDIRK43, BDF) when constructing the PathSim\nSimulation to avoid prohibitively small step sizes.

Parameters

model : pybamm.BaseBatteryModel or None: PyBaMM lithium-ion model. Defaults to isothermal SPMe with heat\nsource calculation enabled.
parameter_values : pybamm.ParameterValues or None: PyBaMM parameter set. Defaults to Chen2020.
initial_soc : float: Initial state of charge (0–1). Default 1.0.
pybamm_solver : pybamm.BaseSolver or None: PyBaMM solver used only during model build / discretisation.\nDefaults to CasadiSolver(mode="safe").

Inputs

I (0) : current [A], positive = discharge\nT_cell (1) : cell temperature [K] from external PathSim thermal block

Outputs

V (0) : terminal voltage [V]\nQ_dot (1) : total heat generation [W]\nSOC (2) : state of charge (0–1)

\n", "source": "class CellElectrical(_CellBase):\n \"\"\"Cell block — electrical outputs only, external thermal coupling.\n\n PathSim integrates the electrochemical state via the discretised PyBaMM\n ODE. Temperature dynamics live outside this block: wire ``Q_dot`` to a\n ``LumpedThermal`` (or similar) block and feed its temperature output back\n to ``T_cell``.\n\n .. note::\n The SPMe/SPM ODE is stiff. Use an implicit solver (e.g.\n ``ESDIRK43``, ``BDF``) when constructing the PathSim\n ``Simulation`` to avoid prohibitively small step sizes.\n\n Parameters\n ----------\n model : pybamm.BaseBatteryModel or None\n PyBaMM lithium-ion model. Defaults to isothermal SPMe with heat\n source calculation enabled.\n parameter_values : pybamm.ParameterValues or None\n PyBaMM parameter set. Defaults to ``Chen2020``.\n initial_soc : float\n Initial state of charge (0–1). Default 1.0.\n pybamm_solver : pybamm.BaseSolver or None\n PyBaMM solver used only during model build / discretisation.\n Defaults to ``CasadiSolver(mode=\"safe\")``.\n\n Inputs\n ------\n I (0) : current [A], positive = discharge\n T_cell (1) : cell temperature [K] from external PathSim thermal block\n\n Outputs\n -------\n V (0) : terminal voltage [V]\n Q_dot (1) : total heat generation [W]\n SOC (2) : state of charge (0–1)\n \"\"\"\n\n _thermal_option = \"isothermal\"\n _thermal_extra_options = {\"calculate heat source for isothermal models\": \"true\"}\n _pybamm_output_vars = [\n \"Terminal voltage [V]\",\n \"Total heating [W]\",\n ]\n\n input_port_labels = {\"I\": 0, \"T_cell\": 1}\n output_port_labels = {\"V\": 0, \"Q_dot\": 1, \"SOC\": 2}", "bases": [ "pathsim_batt.cells.pybamm_cell._CellBase" ], "methods": [], "attributes": [ { "name": "_thermal_option", "description": "", "type": null, "value": "'isothermal'" }, { "name": "_thermal_extra_options", "description": "", "type": null, "value": "{'calculate heat source for isothermal models': 'true'}" }, { "name": "_pybamm_output_vars", "description": "", "type": null, "value": "['Terminal voltage [V]', 'Total heating [W]']" }, { "name": "input_port_labels", "description": "", "type": null, "value": "{'I': 0, 'T_cell': 1}" }, { "name": "output_port_labels", "description": "", "type": null, "value": "{'V': 0, 'Q_dot': 1, 'SOC': 2}" } ], "parameters": [] }, { "name": "CellElectrothermal", "description": "Cell block — coupled electrical and thermal model.", "docstring_html": "

Cell block — coupled electrical and thermal model.

PathSim integrates the full electrochemical + thermal state (via the\ndiscretised PyBaMM ODE). The cell temperature is part of the PyBaMM\nstate vector and is read back as output port T. Supply a\ntime-varying ambient / coolant temperature via T_amb to couple to a\npack-level thermal model.

Note

The SPMe/SPM ODE is stiff. Use an implicit solver (e.g.\nESDIRK43, BDF) when constructing the PathSim\nSimulation to avoid prohibitively small step sizes.

Parameters

model : pybamm.BaseBatteryModel or None: PyBaMM lithium-ion model. Defaults to SPMe(thermal="lumped").
parameter_values : pybamm.ParameterValues or None: PyBaMM parameter set. Defaults to Chen2020.
initial_soc : float: Initial state of charge (0–1). Default 1.0.
pybamm_solver : pybamm.BaseSolver or None: PyBaMM solver used only during model build / discretisation.\nDefaults to CasadiSolver(mode="safe").

Inputs

I (0) : current [A], positive = discharge\nT_amb (1) : ambient / coolant temperature [K]

Outputs

V (0) : terminal voltage [V]\nT (1) : cell temperature [K] (part of PyBaMM state)\nQ_dot (2) : total heat generation [W]\nSOC (3) : state of charge (0–1)

\n", "source": "class CellElectrothermal(_CellBase):\n \"\"\"Cell block — coupled electrical and thermal model.\n\n PathSim integrates the full electrochemical + thermal state (via the\n discretised PyBaMM ODE). The cell temperature is part of the PyBaMM\n state vector and is read back as output port ``T``. Supply a\n time-varying ambient / coolant temperature via ``T_amb`` to couple to a\n pack-level thermal model.\n\n .. note::\n The SPMe/SPM ODE is stiff. Use an implicit solver (e.g.\n ``ESDIRK43``, ``BDF``) when constructing the PathSim\n ``Simulation`` to avoid prohibitively small step sizes.\n\n Parameters\n ----------\n model : pybamm.BaseBatteryModel or None\n PyBaMM lithium-ion model. Defaults to ``SPMe(thermal=\"lumped\")``.\n parameter_values : pybamm.ParameterValues or None\n PyBaMM parameter set. Defaults to ``Chen2020``.\n initial_soc : float\n Initial state of charge (0–1). Default 1.0.\n pybamm_solver : pybamm.BaseSolver or None\n PyBaMM solver used only during model build / discretisation.\n Defaults to ``CasadiSolver(mode=\"safe\")``.\n\n Inputs\n ------\n I (0) : current [A], positive = discharge\n T_amb (1) : ambient / coolant temperature [K]\n\n Outputs\n -------\n V (0) : terminal voltage [V]\n T (1) : cell temperature [K] (part of PyBaMM state)\n Q_dot (2) : total heat generation [W]\n SOC (3) : state of charge (0–1)\n \"\"\"\n\n _thermal_option = \"lumped\"\n _pybamm_output_vars = [\n \"Terminal voltage [V]\",\n \"X-averaged cell temperature [K]\",\n \"Total heating [W]\",\n ]\n\n input_port_labels = {\"I\": 0, \"T_amb\": 1}\n output_port_labels = {\"V\": 0, \"T\": 1, \"Q_dot\": 2, \"SOC\": 3}", "bases": [ "pathsim_batt.cells.pybamm_cell._CellBase" ], "methods": [], "attributes": [ { "name": "_thermal_option", "description": "", "type": null, "value": "'lumped'" }, { "name": "_pybamm_output_vars", "description": "", "type": null, "value": "['Terminal voltage [V]', 'X-averaged cell temperature [K]', 'Total heating [W]']" }, { "name": "input_port_labels", "description": "", "type": null, "value": "{'I': 0, 'T_amb': 1}" }, { "name": "output_port_labels", "description": "", "type": null, "value": "{'V': 0, 'T': 1, 'Q_dot': 2, 'SOC': 3}" } ], "parameters": [] }, { "name": "CellCoSimElectrical", "description": "Cell block (co-simulation) — electrical outputs only, external thermal coupling.", "docstring_html": "

Cell block (co-simulation) — electrical outputs only, external thermal coupling.

PyBaMM advances internally on discrete macro-steps of dt via\npybamm.Simulation.step(). PathSim receives zero-order-held outputs\nbetween macro-steps.

This mode supports PyBaMM models that result in DAE systems (e.g. DFN).

Parameters

model : pybamm.BaseBatteryModel or None: PyBaMM lithium-ion model. Defaults to isothermal SPMe with heat\nsource calculation enabled.
parameter_values : pybamm.ParameterValues or None: PyBaMM parameter set. Defaults to Chen2020.
initial_soc : float: Initial state of charge (0–1). Default 1.0.
pybamm_solver : pybamm.BaseSolver or None: Solver used by PyBaMM for the internal time stepping.\nDefaults to IDAKLUSolver().
dt : float: Co-simulation macro-step size [s]. Must be > 0.

\n", "source": "class CellCoSimElectrical(_CoSimCellBase):\n \"\"\"Cell block (co-simulation) — electrical outputs only, external thermal coupling.\n\n PyBaMM advances internally on discrete macro-steps of ``dt`` via\n ``pybamm.Simulation.step()``. PathSim receives zero-order-held outputs\n between macro-steps.\n\n This mode supports PyBaMM models that result in DAE systems (e.g. DFN).\n\n Parameters\n ----------\n model : pybamm.BaseBatteryModel or None\n PyBaMM lithium-ion model. Defaults to isothermal SPMe with heat\n source calculation enabled.\n parameter_values : pybamm.ParameterValues or None\n PyBaMM parameter set. Defaults to ``Chen2020``.\n initial_soc : float\n Initial state of charge (0–1). Default 1.0.\n pybamm_solver : pybamm.BaseSolver or None\n Solver used by PyBaMM for the internal time stepping.\n Defaults to ``IDAKLUSolver()``.\n dt : float\n Co-simulation macro-step size [s]. Must be > 0.\n \"\"\"\n\n _thermal_option = \"isothermal\"\n _thermal_extra_options = {\"calculate heat source for isothermal models\": \"true\"}\n _pybamm_output_vars = [\n \"Terminal voltage [V]\",\n \"Total heating [W]\",\n ]\n\n input_port_labels = {\"I\": 0, \"T_cell\": 1}\n output_port_labels = {\"V\": 0, \"Q_dot\": 1, \"SOC\": 2}", "bases": [ "pathsim_batt.cells.pybamm_cell._CoSimCellBase" ], "methods": [], "attributes": [ { "name": "_thermal_option", "description": "", "type": null, "value": "'isothermal'" }, { "name": "_thermal_extra_options", "description": "", "type": null, "value": "{'calculate heat source for isothermal models': 'true'}" }, { "name": "_pybamm_output_vars", "description": "", "type": null, "value": "['Terminal voltage [V]', 'Total heating [W]']" }, { "name": "input_port_labels", "description": "", "type": null, "value": "{'I': 0, 'T_cell': 1}" }, { "name": "output_port_labels", "description": "", "type": null, "value": "{'V': 0, 'Q_dot': 1, 'SOC': 2}" } ], "parameters": [] }, { "name": "CellCoSimElectrothermal", "description": "Cell block (co-simulation) — coupled electrical and thermal model.", "docstring_html": "

Cell block (co-simulation) — coupled electrical and thermal model.

PyBaMM advances internally on discrete macro-steps of dt via\npybamm.Simulation.step(). PathSim receives zero-order-held outputs\nbetween macro-steps.

This mode supports PyBaMM models that result in DAE systems (e.g. DFN).

Parameters

model : pybamm.BaseBatteryModel or None: PyBaMM lithium-ion model. Defaults to SPMe(thermal="lumped").
parameter_values : pybamm.ParameterValues or None: PyBaMM parameter set. Defaults to Chen2020.
initial_soc : float: Initial state of charge (0–1). Default 1.0.
pybamm_solver : pybamm.BaseSolver or None: Solver used by PyBaMM for the internal time stepping.\nDefaults to IDAKLUSolver().
dt : float: Co-simulation macro-step size [s]. Must be > 0.

\n", "source": "class CellCoSimElectrothermal(_CoSimCellBase):\n \"\"\"Cell block (co-simulation) — coupled electrical and thermal model.\n\n PyBaMM advances internally on discrete macro-steps of ``dt`` via\n ``pybamm.Simulation.step()``. PathSim receives zero-order-held outputs\n between macro-steps.\n\n This mode supports PyBaMM models that result in DAE systems (e.g. DFN).\n\n Parameters\n ----------\n model : pybamm.BaseBatteryModel or None\n PyBaMM lithium-ion model. Defaults to ``SPMe(thermal=\"lumped\")``.\n parameter_values : pybamm.ParameterValues or None\n PyBaMM parameter set. Defaults to ``Chen2020``.\n initial_soc : float\n Initial state of charge (0–1). Default 1.0.\n pybamm_solver : pybamm.BaseSolver or None\n Solver used by PyBaMM for the internal time stepping.\n Defaults to ``IDAKLUSolver()``.\n dt : float\n Co-simulation macro-step size [s]. Must be > 0.\n \"\"\"\n\n _thermal_option = \"lumped\"\n _pybamm_output_vars = [\n \"Terminal voltage [V]\",\n \"X-averaged cell temperature [K]\",\n \"Total heating [W]\",\n ]\n\n input_port_labels = {\"I\": 0, \"T_amb\": 1}\n output_port_labels = {\"V\": 0, \"T\": 1, \"Q_dot\": 2, \"SOC\": 3}", "bases": [ "pathsim_batt.cells.pybamm_cell._CoSimCellBase" ], "methods": [], "attributes": [ { "name": "_thermal_option", "description": "", "type": null, "value": "'lumped'" }, { "name": "_pybamm_output_vars", "description": "", "type": null, "value": "['Terminal voltage [V]', 'X-averaged cell temperature [K]', 'Total heating [W]']" }, { "name": "input_port_labels", "description": "", "type": null, "value": "{'I': 0, 'T_amb': 1}" }, { "name": "output_port_labels", "description": "", "type": null, "value": "{'V': 0, 'T': 1, 'Q_dot': 2, 'SOC': 3}" } ], "parameters": [] } ], "functions": [ { "name": "_prepare_parameter_values", "description": "Copy *parameter_values* (defaulting to Chen2020) and mark both", "docstring_html": "

Copy parameter_values (defaulting to Chen2020) and mark both\ndriving inputs as PyBaMM "[input]" placeholders.

\n", "source": "def _prepare_parameter_values(\n parameter_values: pybamm.ParameterValues | None,\n) -> pybamm.ParameterValues:\n \"\"\"Copy *parameter_values* (defaulting to Chen2020) and mark both\n driving inputs as PyBaMM ``\"[input]\"`` placeholders.\"\"\"\n if parameter_values is None:\n parameter_values = pybamm.ParameterValues(\"Chen2020\")\n parameter_values = parameter_values.copy()\n parameter_values[\"Current function [A]\"] = \"[input]\"\n parameter_values[\"Ambient temperature [K]\"] = \"[input]\"\n return parameter_values", "signature": "(parameter_values: pybamm.ParameterValues | None)", "parameters": [ { "name": "parameter_values", "type": "pybamm.ParameterValues | None", "default": null, "description": "" } ], "returns": "pybamm.ParameterValues" } ] }, "pathsim_batt.thermal.lumped": { "name": "pathsim_batt.thermal.lumped", "description": "", "docstring_html": "", "classes": [ { "name": "LumpedThermal", "description": "Single-node lumped thermal model.", "docstring_html": "

Single-node lumped thermal model.

\nm C_p \\frac{dT}{dt} = \\dot{Q} - U A (T - T_{\\mathrm{amb}})\n

Parameters

mass : float: Thermal mass [kg].
Cp : float: Specific heat capacity [J kg⁻¹ K⁻¹].
UA : float: Overall heat transfer conductance [W K⁻¹].
T0 : float: Initial temperature [K]. Default 298.15 K.

Inputs

Q_dot (0) : heat generation rate [W]\nT_amb (1) : ambient temperature [K]

Outputs

T (0) : cell temperature [K]

\n", "source": "class LumpedThermal(DynamicalSystem):\n \"\"\"Single-node lumped thermal model.\n\n .. math::\n\n m C_p \\\\frac{dT}{dt} = \\\\dot{Q} - U A (T - T_{\\\\mathrm{amb}})\n\n Parameters\n ----------\n mass : float\n Thermal mass [kg].\n Cp : float\n Specific heat capacity [J kg⁻¹ K⁻¹].\n UA : float\n Overall heat transfer conductance [W K⁻¹].\n T0 : float\n Initial temperature [K]. Default 298.15 K.\n\n Inputs\n ------\n Q_dot (0) : heat generation rate [W]\n T_amb (1) : ambient temperature [K]\n\n Outputs\n -------\n T (0) : cell temperature [K]\n \"\"\"\n\n input_port_labels = {\"Q_dot\": 0, \"T_amb\": 1}\n output_port_labels = {\"T\": 0}\n\n def __init__(self, mass=0.065, Cp=750.0, UA=0.5, T0=298.15):\n # input validation\n if mass <= 0:\n raise ValueError(f\"'mass' must be positive but is {mass}\")\n if Cp <= 0:\n raise ValueError(f\"'Cp' must be positive but is {Cp}\")\n if UA < 0:\n raise ValueError(f\"'UA' must be non-negative but is {UA}\")\n\n # store parameters\n self.mass = float(mass)\n self.Cp = float(Cp)\n self.UA = float(UA)\n\n def _fn_d(x, u, t):\n (T,) = x\n Q_dot, T_amb = u\n return np.array([(Q_dot - self.UA * (T - T_amb)) / (self.mass * self.Cp)])\n\n def _fn_a(x, u, t):\n return x.copy()\n\n super().__init__(\n func_dyn=_fn_d,\n func_alg=_fn_a,\n initial_value=np.array([float(T0)]),\n )", "bases": [ "pathsim.blocks.DynamicalSystem" ], "methods": [ { "name": "__init__", "description": "", "docstring_html": "", "source": "def __init__(self, mass=0.065, Cp=750.0, UA=0.5, T0=298.15):\n # input validation\n if mass <= 0:\n raise ValueError(f\"'mass' must be positive but is {mass}\")\n if Cp <= 0:\n raise ValueError(f\"'Cp' must be positive but is {Cp}\")\n if UA < 0:\n raise ValueError(f\"'UA' must be non-negative but is {UA}\")\n\n # store parameters\n self.mass = float(mass)\n self.Cp = float(Cp)\n self.UA = float(UA)\n\n def _fn_d(x, u, t):\n (T,) = x\n Q_dot, T_amb = u\n return np.array([(Q_dot - self.UA * (T - T_amb)) / (self.mass * self.Cp)])\n\n def _fn_a(x, u, t):\n return x.copy()\n\n super().__init__(\n func_dyn=_fn_d,\n func_alg=_fn_a,\n initial_value=np.array([float(T0)]),\n )", "signature": "(mass = 0.065, Cp = 750.0, UA = 0.5, T0 = 298.15)", "parameters": [ { "name": "mass", "type": null, "default": "0.065", "description": "" }, { "name": "Cp", "type": null, "default": "750.0", "description": "" }, { "name": "UA", "type": null, "default": "0.5", "description": "" }, { "name": "T0", "type": null, "default": "298.15", "description": "" } ], "returns": null, "method_type": "method" } ], "attributes": [ { "name": "input_port_labels", "description": "", "type": null, "value": "{'Q_dot': 0, 'T_amb': 1}" }, { "name": "output_port_labels", "description": "", "type": null, "value": "{'T': 0}" }, { "name": "mass", "description": "", "type": null, "value": "float(mass)" }, { "name": "Cp", "description": "", "type": null, "value": "float(Cp)" }, { "name": "UA", "description": "", "type": null, "value": "float(UA)" } ], "parameters": [ { "name": "mass", "type": null, "default": "0.065", "description": "Thermal mass [kg]." }, { "name": "Cp", "type": null, "default": "750.0", "description": "Specific heat capacity [J kg⁻¹ K⁻¹]." }, { "name": "UA", "type": null, "default": "0.5", "description": "Overall heat transfer conductance [W K⁻¹]." }, { "name": "T0", "type": null, "default": "298.15", "description": "Initial temperature [K]. Default 298.15 K." } ] } ], "functions": [] } } }