Merge pull request #1683 from pybamm-team/ecker-thermal-params

Ecker thermal params

Author: rtimms

CHANGELOG.md

@@ -2,6 +2,7 @@
 
 ## Features
 
+-   Added thermal parameters (thermal conductivity, specific heat, etc.) to the `Ecker2015` parameter set from Zhao et al. (2018) and Hales et al. (2019) ([#1683](https://github.com/pybamm-team/PyBaMM/pull/1683))
 -   Added `plot_summary_variables` to plot and compare summary variables ([#1678](https://github.com/pybamm-team/PyBaMM/pull/1678))
 -   The DFN model can now be used directly (instead of `BasicDFNHalfCell`) to simulate a half-cell ([#1600](https://github.com/pybamm-team/PyBaMM/pull/1600))
 
@@ -26,7 +27,7 @@ example notebook ([#1602](https://github.com/pybamm-team/PyBaMM/pull/1602))
 -   Added LG M50 (NMC811 and graphite + SiOx) parameter set from O'Regan 2021 ([#1594](https://github.com/pybamm-team/PyBaMM/pull/1594))
 -   `pybamm.base_solver.solve` function can take a list of input parameters to calculate the sensitivities of the solution with respect to. Alternatively, it can be set  to `True` to calculate the sensitivities for all input parameters ([#1552](https://github.com/pybamm-team/PyBaMM/pull/1552))
 -   Added capability for `quaternary` domains (in addition to `primary`, `secondary` and `tertiary`), increasing the maximum number of domains that a `Symbol` can have to 4. ([#1580](https://github.com/pybamm-team/PyBaMM/pull/1580))
--   Tabs can now be placed at the bottom of the cell in 1+1D thermal models ([#1581](https://github.com/pybamm-team/PyBaMM/pull/1581)) 
+-   Tabs can now be placed at the bottom of the cell in 1+1D thermal models ([#1581](https://github.com/pybamm-team/PyBaMM/pull/1581))
 -   Added temperature dependence on electrode electronic conductivity ([#1570](https://github.com/pybamm-team/PyBaMM/pull/1570))
 -   `pybamm.base_solver.solve` function can take a list of input parameters to calculate the sensitivities of the solution with respect to. Alternatively, it can be set  to `True` to calculate the sensitivities for all input parameters ([#1552](https://github.com/pybamm-team/PyBaMM/pull/1552))
 -   Added a new lithium-ion model `MPM` or Many-Particle Model, with a distribution of particle sizes in each electrode. ([#1529](https://github.com/pybamm-team/PyBaMM/pull/1529))

docs/install/windows-wsl.rst

@@ -13,9 +13,8 @@ Install WSL
 Follow the instructions from Microsoft
 `here <https://docs.microsoft.com/en-us/windows/wsl/install-win10>`__.
 When given the option, choose the Ubuntu 18.04 LTS distribution to
-install. Don’t forget to initialise the Ubuntu installation using the
-instructions given
-`here <https://docs.microsoft.com/en-us/windows/wsl/initialize-distro>`__.
+install. Best practices for setting up a WSL development environment can be found
+`here <https://docs.microsoft.com/en-us/windows/wsl/setup/environment>`__.
 
 Install PyBaMM
 --------------
@@ -50,7 +49,7 @@ directory in bash using the ``cd`` command:
 
 .. code:: bash
 
-   cd PyBaMM 
+   cd PyBaMM
 
 If you are unfamiliar with the linux command line, you might find it
 useful to work through this

pybamm/CITATIONS.txt

@@ -117,6 +117,17 @@
   journal = {Journal of Open Source Software},
 }
 
+@article{Hales2019,
+  title={The cell cooling coefficient: a standard to define heat rejection from lithium-ion batteries},
+  author={Hales, Alastair and Diaz, Laura Bravo and Marzook, Mohamed Waseem and Zhao, Yan and Patel, Yatish and Offer, Gregory},
+  journal={Journal of The Electrochemical Society},
+  volume={166},
+  number={12},
+  pages={A2383},
+  year={2019},
+  publisher={IOP Publishing}
+}
+
 @article{Harris2020,
   title = {{Array programming with NumPy}},
   author = {Harris, Charles R. and Millman, K. Jarrod and van der Walt, St{\'{e}}fan J. and Gommers, Ralf and Virtanen, Pauli and Cournapeau, David and Wieser, Eric and Taylor, Julian and Berg, Sebastian and Smith, Nathaniel J. and others},
@@ -459,4 +470,15 @@
   year = {2017},
   publisher = {Elsevier},
   doi = {10.1016/j.jpowsour.2017.05.110},
-}
+}
+
+@article{Zhao2018,
+  title={Modeling the effects of thermal gradients induced by tab and surface cooling on lithium ion cell performance},
+  author={Zhao, Yan and Patel, Yatish and Zhang, Teng and Offer, Gregory J},
+  journal={Journal of The Electrochemical Society},
+  volume={165},
+  number={13},
+  pages={A3169},
+  year={2018},
+  publisher={IOP Publishing}
+}

pybamm/input/parameters/lithium_ion/cells/kokam_Ecker2015/README.md

@@ -5,3 +5,11 @@ Parameters for a Kokam SLPB 75106100 cell, from the papers
 > Ecker, Madeleine, et al. "Parameterization of a physico-chemical model of a lithium-ion battery I. determination of parameters." Journal of the Electrochemical Society 162.9 (2015): A1836-A1848.
 
 >Ecker, Madeleine, et al. "Parameterization of a physico-chemical model of a lithium-ion battery II. Model validation." Journal of The Electrochemical Society 162.9 (2015): A1849-A1857.
+
+The tab placement parameters are taken from measurements in
+
+> Hales, Alastair, et al. "The cell cooling coefficient: a standard to define heat rejection from lithium-ion batteries." Journal of The Electrochemical Society 166.12 (2019): A2383.
+
+The thermal material properties are for a 5 Ah power pouch cell by Kokam. The data are extracted from
+
+> Zhao, Y., et al. "Modeling the effects of thermal gradients induced by tab and surface cooling on lithium ion cell performance."" Journal of The Electrochemical Society, 165.13 (2018): A3169-A3178.

pybamm/input/parameters/lithium_ion/cells/kokam_Ecker2015/parameters.csv

@@ -6,13 +6,36 @@ Negative current collector thickness [m],1.40E-05,,
 Negative electrode thickness [m],7.4E-05,,
 Separator thickness [m],2E-05,,
 Positive electrode thickness [m],5.4E-05,,
-Positive current collector thickness [m],2.5E-05,,
+Positive current collector thickness [m],1.5E-05,,
 Electrode height [m],1.01E-01,,
 Electrode width [m],8.50E-02,,
-Cell cooling surface area [m2],1.72E-2,,pouch (single layer)
-Cell volume [m3],1.61E-6,,pouch (single layer)
+Negative tab width [m],7E-3,,Hales et al. (2019)
+Negative tab centre y-coordinate [m],4.5E-3,Hales et al. (2019),
+Negative tab centre z-coordinate [m],1.01E-1,Top of cell,
+Positive tab width [m],6.9E-3,Hales et al. (2019),
+Positive tab centre y-coordinate [m],30.9E-3,Hales et al. (2019),
+Positive tab centre z-coordinate [m],1.01E-1,Top of cell,
+Cell cooling surface area [m2],1.72E-2,pouch (single layer),
+Cell volume [m3],1.52E-6,pouch (single layer),
+,,,
+,,,
+# Current collector properties ,,,
+Negative current collector conductivity [S.m-1],58411000,CRC Handbook,copper
+Positive current collector conductivity [S.m-1],36914000,CRC Handbook,aluminium
+,,,
+# Density,,,
+Negative current collector density [kg.m-3],8933,Zhao et al. (2018),
+Positive current collector density [kg.m-3],2702,Zhao et al. (2018),
+,,,
+# Specific heat capacity,,,
+Negative current collector specific heat capacity [J.kg-1.K-1],385,Zhao et al. (2018),
+Positive current collector specific heat capacity [J.kg-1.K-1],903,Zhao et al. (2018),
+,,,
+# Thermal conductivity,,,
+Negative current collector thermal conductivity [W.m-1.K-1],398,Zhao et al. (2018),
+Positive current collector thermal conductivity [W.m-1.K-1],238,Zhao et al. (2018),
 ,,,
 # Electrical,,,
 Nominal cell capacity [A.h], 0.15625, 7.5/48 (parameter set for a single layer cell),
 Typical current [A], 0.15652,,
-Current function [A],0.15652,default current function,
+Current function [A],0.15652,default current function (1C),

pybamm/input/parameters/lithium_ion/experiments/1C_discharge_from_full_Ecker2015/parameters.csv

@@ -3,12 +3,12 @@ Name [units],Value,Reference,Notes
 ,,,
 # Temperature
 Reference temperature [K],296.15,23C,
-Negative current collector surface heat transfer coefficient [W.m-2.K-1],0,Paper does not consider thermal effects
-Positive current collector surface heat transfer coefficient [W.m-2.K-1],0,Paper does not consider thermal effects
-Negative tab heat transfer coefficient [W.m-2.K-1],10,Paper does not consider thermal effects,
-Positive tab heat transfer coefficient [W.m-2.K-1],10,Paper does not consider thermal effects,
-Edge heat transfer coefficient [W.m-2.K-1],0.3,Paper does not consider thermal effects,
-Total heat transfer coefficient [W.m-2.K-1],10,Paper does not consider thermal effects,
+Negative current collector surface heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
+Positive current collector surface heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
+Negative tab heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
+Positive tab heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
+Edge heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
+Total heat transfer coefficient [W.m-2.K-1],10,Assume uniform heat loss,
 Ambient temperature [K],298.15,,
 
 ,,,

pybamm/input/parameters/lithium_ion/negative_electrodes/graphite_Ecker2015/README.md

@@ -10,3 +10,7 @@ The fits to data for the electrode and electrolyte properties are those provided
 by Dr. Simon O’Kane in the paper:
 
 > Richardson, Giles, et. al. "Generalised single particle models for high-rate operation of graded lithium-ion electrodes: Systematic derivation and validation." Electrochemica Acta 339 (2020): 135862
+
+The thermal material properties are for a 5 Ah power pouch cell by Kokam. The data are extracted from
+
+> Zhao, Y., et al. "Modeling the effects of thermal gradients induced by tab and surface cooling on lithium ion cell performance."" Journal of The Electrochemical Society, 165.13 (2018): A3169-A3178.

pybamm/input/parameters/lithium_ion/negative_electrodes/graphite_Ecker2015/parameters.csv

@@ -21,5 +21,10 @@ Negative electrode cation signed stoichiometry,-1,,
 Negative electrode electrons in reaction,1,,
 Negative electrode exchange-current density [A.m-2],[function]graphite_electrolyte_exchange_current_density_Ecker2015,,
 ,,,
+# Density,,,
+Negative electrode density [kg.m-3],1555,Zhao et al. (2018),
+,,,
 # Thermal parameters,,,
+Negative electrode specific heat capacity [J.kg-1.K-1],1437,Zhao et al. (2018),
+Negative electrode thermal conductivity [W.m-1.K-1],1.58,Zhao et al. (2018),
 Negative electrode OCP entropic change [V.K-1],0,,

pybamm/input/parameters/lithium_ion/positive_electrodes/LiNiCoO2_Ecker2015/README.md

@@ -10,3 +10,7 @@ The fits to data for the electrode and electrolyte properties are those provided
 by Dr. Simon O’Kane in the paper:
 
 > Richardson, Giles, et. al. "Generalised single particle models for high-rate operation of graded lithium-ion electrodes: Systematic derivation and validation." Electrochemica Acta 339 (2020): 135862
+
+The thermal material properties are for a 5 Ah power pouch cell by Kokam. The data are extracted from
+
+> Zhao, Y., et al. "Modeling the effects of thermal gradients induced by tab and surface cooling on lithium ion cell performance."" Journal of The Electrochemical Society, 165.13 (2018): A3169-A3178.

pybamm/input/parameters/lithium_ion/positive_electrodes/LiNiCoO2_Ecker2015/parameters.csv

@@ -20,7 +20,11 @@ Positive electrode exchange-current density [A.m-2],[function]nco_electrolyte_ex
 # Interfacial reactions,,,
 Positive electrode cation signed stoichiometry,-1,,
 Positive electrode electrons in reaction,1,,
-
+,,,
+# Density,,,
+Positive electrode density [kg.m-3],2895,Zhao et al. (2018),
 ,,,
 # Thermal parameters,,,
+Positive electrode specific heat capacity [J.kg-1.K-1],1270,Zhao et al. (2018),
+Positive electrode thermal conductivity [W.m-1.K-1],1.04,Zhao et al. (2018),
 Positive electrode OCP entropic change [V.K-1],0,,