pybamm-team
Variable problem of lithium_ionModel
pybamm.lithium_ion.DFN(). Is there a problem with the expression of the two variables related to the heating of positive and negative current collectors? Their heating power is only on the order of 1E-5.
I'm not sure if this detail is a bug, but I think the amount of heat generated by the current collector is a bit small.
The calculations are correct. For example, for a 10 um thick copper collector at a current density of 3 mA/cm2 (one side coat), the ohmic heating is 1.55E-5 W/m3. At the same time, the value completely disagrees with reality. It should be several orders of magnitude higher. This is due to a simple reason. The calculated ohmic heating is for conduction through the (collector) plane; in practice, the heating is related to conduction in the (collector) plane. Getting the correct value would require solving an additional PDE in 2D. I believe that as it stands, ohmic heating for current collectors is a bit misleading, but I guess it is required.
