A composite electrode model with thermal model

[Ruby8350]

Hello everyone,

I’m trying to run the composite electrode model with the thermal submodel enabled, and I’d like to confirm which additional parameters are required for this setup. I’ve already specified the basic electrochemical properties for the primary/secondary phases (e.g., diffusivity, OCP, and exchange current density).

Here are the model options I used:
options = {
"cell geometry": "pouch",
"thermal": "x-full",
"contact resistance": "true",
"particle phases": ("2", "1"),
"open-circuit potential": (("single", "current sigmoid"), "single"),
}
model = pybamm.lithium_ion.DFN(options=options)

However, I get the following error:
ModelError: Missing variable for submodel 'thermal': 'Negative electrode primary volumetric interfacial current density [A.m-3]'. Check the selected submodels provide all of the required variables.