Combined electrochemical and thermal modelling with COMSOL ...

29 downloads 0 Views 3MB Size Report
COMSOL multiphysics. Comparison of experimental and simulation results. Summary and outlook. Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando ...
Combination of electrochemical-calorimetric studies on cylindrical lithium ion cells and thermal modelling by COMSOL Multiphysics software C. Ziebert1, A. Ossipova2, M. Rohde1, H.J. Seifert1 1Karlsruhe Institute of Technology (KIT), IAM-AWP 2Ruprecht-Karls University Heidelberg, faculty of physics and astronomy

Institute for Applied Materials – Applied Materials Physics (IAM-AWP)

KIT – University of the State of Baden-Wuerttemberg and National Research Center of the Helmholtz Association

www.kit.edu

Outline

Motivation Introduction Isoperibolic studies Determination of heat dissipation rate and total generated heat Adiabatic studies Combined electrochemical and thermal modelling with COMSOL multiphysics Comparison of experimental and simulation results Summary and outlook

2

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Motivation: Accidents with lithium ion batteries (courtesy Reuters)

(courtesy www.cbsnews.com)

Rechargeable laptop battery – overheating or internal short circuit

(courtesy Auto Bild)

Boeing Dreamliner: – reasons unclear

Chevy Volt: – mechanical impact

Cell Thermodynamics and Safety Studies

- Thermodynamics and phase diagrams strongly influence battery performance - Structural stabilities of active materials need to be taken into account - Battery safety is essential (thermal runaway, abuse, accident) - Thermal management of batteries is important and needs improvement by combination of experimental studies and simulation 3

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Introduction IKEBA Project Integrated Components and Integrated Design of Energy Efficient Battery Systems

Battery Monitoring Automotive uC

Battery Monitoring

5 cooperating partners 4

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

Battery Management Software

Start time: May 2013 Planned duration: 3 years Budget: 7 Million Euro KIT, IAM-AWP

Short institute presentation

Karlsruhe Institute of Technology (KIT) 2006 - Excellence Initiative of the Federal Republic of Germany and the Federal States 1/10/2009 - Foundation of KIT Merger of University of Karlsruhe (TH) and the Forschungszentrum Karlsruhe GmbH Data 2012 Employees: Students: Professors: Budget: Patents: Spin-offs: Campus South

9.254 23.905 359 785 Mio Euro 72 18 Campus North

Research – Teaching – Innovation 5

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Combined electrochemical and thermal characterization at IAM-AWP Enthalpies, heat capacities, heat generation and heat dissipation - Accelerating Rate Calorimeters (ARC, THT / C3 Analysentechnik) - Isothermal Battery Calorimeter (IBC, THT / C3 Analysentechnik) - Differential Scanning Calorimeters (DSC, NETZSCH) - High Temperature Calvet Drop Solution Calorimeter (Alexsys 1000, Setaram) AlexSys 1000

Accelerating Rate (ARC) Battery Calorimeters 6

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

Isothermal Battery Calorimeter (IBC)

Digatron MCT Cell Tester (50 A, 6 V, 24 channels KIT, IAM-AWP

Accelerating Rate Calorimeter (ARC) ARC provides both isoperibolic and adiabatic environments Under isoperibolic conditions the environmental temperature is held constant. Under adiabatic conditions the cell can be studied under conditions of negligible heat loss.

Cell Specifications Cell type:

18650

Cathode material: LiMn2O4 (LMO) Nominal capacity: 1400 mAh

7

Nominal voltage:

3.7 V

Voltage window:

4.2 V to 2.5 V

Core cell weight:

41 g

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Battery Calorimeter

ES-ARC: Ø: 10 cm h: 10 cm

EV-ARC: Ø: 25 cm h: 50 cm

+

EV+ ARC: Ø: 40 cm h: 44 cm

ARC combined with internal or external cycler = Battery Calorimeter

8

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Isoperibolic studies: cycling parameters

Isoperibolic cycling at 30 °C and 1C charge/discharge rate

Current in A Voltage in V

Discharge parameter: - method: constant current (CC) - Umin = 2.5 V - I = 1.4 A → 1C-rate

2

Temperature in °C

Charge parameter: - method: constant current, constant voltage (CCCV) - Umax = 4.2 V - I = 1.4 A → 1C-rate - Imin = 0.14 A

1 0 -1 -2 4.2 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4 34 32 30 0

200

400

600

800 1000 1200

Time in min 9

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Determination of heat dissipation rate and total generated heat Key data for thermal management and safety Conversion of thermal data (temperature, temperature rate) to heat (Joule) and power (Watt) with the aim of understanding of heat release to determine heat removal requirements for thermal management. Energy balance for a cell inside a calorimeter 𝑞𝑞̇ 𝑇𝑇 = 𝑚𝑚 𝑐𝑐𝑝𝑝 Total heat generation rate 𝑞𝑞̇ 𝑇𝑇 Total heat generation

𝑑𝑑𝑑𝑑 − αARC (𝑇𝑇 − 𝑇𝑇𝑒𝑒𝑒𝑒𝑒𝑒 ) 𝑑𝑑𝑑𝑑

Enthalpy accumulation rate 𝑞𝑞̇ 𝑎𝑎𝑐𝑐𝑐𝑐

Heat dissipation rate 𝑞𝑞̇ 𝑑𝑑

𝑡𝑡

𝑞𝑞𝑇𝑇 = 𝑞𝑞𝑑𝑑 +𝑞𝑞𝑎𝑎𝑎𝑎𝑎𝑎 = −αARC ∫0 (𝑇𝑇 − 𝑇𝑇𝑒𝑒𝑒𝑒𝑒𝑒 ) 𝑑𝑑𝑑𝑑 + 𝑚𝑚 𝑐𝑐𝑝𝑝 (𝑇𝑇 𝑡𝑡 − 𝑇𝑇 0 )

S. Hallaj, H. Maleki, J.S. Hong, J.R. Selman, Journal of Power Sources 83, 1-8 (1999)

10

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Effective specific heat capacity U = 3,3V

11

𝑐𝑐𝑝𝑝 = 1.03468 – 0.00718372 ∙ 𝑇𝑇 + 0.00013056 ∙ 𝑇𝑇 2 e.g. at 30 °C cp = 0.94 J/g K Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

𝑐𝑐𝑝𝑝 =

𝑞𝑞 𝑚𝑚∆𝑇𝑇

𝑞𝑞̇ 𝑎𝑎𝑎𝑎𝑎𝑎

𝑞𝑞 = � 𝑈𝑈 ∙ 𝐼𝐼𝐼𝐼𝐼𝐼

𝑑𝑑𝑑𝑑 = 𝑚𝑚 ∙ 𝒄𝒄𝒑𝒑 𝑑𝑑𝑑𝑑

KIT, IAM-AWP

Measurement of calorimeter constant Energy balance for lumped heat transfer system under natural convective heating conditions and 𝐵𝐵𝑖𝑖 ≪ 0.1

60

Solution:

50

2

40

1 Dummy Cell: cylinder of AlMgSi0.5 EN AW-6060 F22 with same dimensions as 18650 cell

30 20

0

100 200 300 400 500 600 Time in min

cp = 0.9 J/g K

m = 46.02 g

at t = 0

ln(𝑇𝑇 − 𝑇𝑇 𝑒𝑒𝑒𝑒𝑒𝑒 )=

3

ln(Tenv-T)

Temperature in °C

Initial condition: 𝑇𝑇 = 𝑇𝑇 𝑒𝑒𝑒𝑒𝑒𝑒

Tenv T

70

𝑑𝑑𝑑𝑑

m 𝑐𝑐𝑝𝑝 𝑑𝑑𝑑𝑑 = αARC(𝑇𝑇 − 𝑇𝑇𝑒𝑒𝑒𝑒𝑒𝑒 )

αARC

𝑚𝑚 𝑐𝑐𝑝𝑝

𝑡𝑡

0 -1 -2 -3

0

A = 3.674 ⋅10 -3 m2

1000

2000 Time in s

3000

4000

Calorimeter constant αARC = (0.04 ± 0.001) W/K

12

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Heat dissipation and total generated heat

Heat dissipation rate in mW

Tenv = 30°C, 1 C charge/1 C discharge 200 150 100 50 0

0

200

Results of integration over one discharge cycle:

cp = 0.94 J/g K

𝑚𝑚 = 43 g

𝑇𝑇(t) = 39.3 °C

𝑇𝑇(0) = 30.13 °C

Total heat generated during 1C discharge at 30 °C: 13

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

400 600 800 1000 1200 Time in min

𝑞𝑞𝑑𝑑 = -322.3 J

𝑞𝑞𝑎𝑎𝑐𝑐𝑐𝑐 = 370.7 J 𝑞𝑞𝑇𝑇 = 48.4 J

KIT, IAM-AWP

C-rate dependence of heat dissipation rate

4.2 0.7A (0.5C) 1.4A (1C) 2.8A (2C) 4.2A (3C)

50 °C

4.0 3.6 3.4 3.2 3.0 2.6

0.4

0.6

0.8

1.0

1.2

1.4

Discharge in Ah

C-Rate dependence of capacity vs. voltage curve during isoperibolic cycling at 50 °C

150

0.25C

100 50 0

0

200

0.5C

{

0.2

1C

200

{

2.4 0.0

250

{

2.8

Heat dissipation rate in mW

Voltage in V

3.8

400 600 Time in min

800

1000

C-rate dependence of heat dissipation rate during isoperibolic cycling at 30 °C

14

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Adiabatic studies: working principle

Why ADI ABATI C? With conditions of NO HEAT LOSS or gain …data showing heat release ie exotherm reaction.... is given as ‘WORST CASE’ DATA Therefore the ARC gives a SIMULATION of what may happen WORST CASE data may then be extrapolated to any LARGER scale As there is no heat loss, the data can be treated to simulate what will happen for any heat loss scenario

Sample Calorimeter

Source: Thermal Hazard Technology (THT)

15

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Adiabatic studies: cycling parameters

Adiabatic cycling at 30 °C and 1C charge/discharge rate

Current in A Voltage in V

Discharge parameter: - method: constant current (CC) - Umin = 2.5 V - I = 1.4 A → 1C-rate

2

Temperature in °C

Charge parameter: - method: constant current, constant voltage (CCCV) - Umax = 4.2 V - I = 1.4 A → 1C-rate - Imin = 0.14 A

1 0 -1 -2 4.2 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4 80 70 60 50 40 30 0

200

400

600

Time in min

16

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Combined electrochemical and thermal modelling with COMSOL multiphysics 1d – electrochemical

Diffusion of Li-ions in the electrodes

2d – thermal (axial-symmetric)

U(t), I(t), SOC

Diffusion of Li-ions in the electrolyte Ohmic Losses

Thermal conductivity Heat capacity

T(t), Q(t)

Heat transport – Temperature distribution

L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 17

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Governing equations 1d – electrochemical model Material balance for Li-ions in the electrodes

2d – thermal model Thermal conductivity (radial)

(Fick‘s second law)

Material balance for Li-ions in the electrolyte

Ohm‘s law for Li-ions

Heat capacity

Energy balance

r radius, c concentration, t + transfer coefficient, ε porosity, a area, J current density, D diffusion constant, φ potential, F Faraday-constant, σ electric conductivity, Li thickness of ithlayer of spirally wounded cell, λ thermal conductivity 0

L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 18

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Comparison of experimental and simulation results Tenv = 40 °C, 1 C Charge/1 C Discharge Current density in A/m2

1 0 -1 4.2 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4

Isoperibolic

Voltage in V

Voltage in V

Current in A

Experiment

∆T in K

∆T in K

2

20 10 0 -10 -20 4.2 4.0 3.8 3.6 3.4 3.2 3.0

2 0

0

200

400

600

800

1000

Time in min 19

Simulation

4

4

0

(1 C ~ 20 A/m2 )

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

0

200

400

600

800

1000

Time in min KIT, IAM-AWP

Comparison of experimental and simulation results Tst = 25 °C, 1 C Charge/1 C Discharge Experiment

Simulation 20 Current in A

1 0 -1

40

Voltage in V

Adiabatic

30 20

0 -10

4.2 4.0 3.8 3.6 3.4 3.2 3.0 50 40 30 20 10

10 0

100 200 300 400 500 600 700 Time in min

20

10

-20

-2 4.2 4.0 3.8 3.6 3.4 3.2 3.0 2.8 2.6 2.4 50

∆T in K

∆T in K

Voltage in V

Current in A

2

0

(1 C ~ 20 A/m2 )

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

0

0

100 200 300 400 500 600 700 Time in min KIT, IAM-AWP

C-rate dependence of maximum temperature increase during discharge isoperibolic

∆T in K

Heat dissipation qd [mW]

Simulation

Experiment 550 500 450 400 350 300 250 200 150 100 50 0

3C 2C

1C 0.5C

0

200

400

∆Tmax in K

Time in s 45 40 35 30 25 20 15 10 5 0

Experiment Simulation

0 21

600 800 1000 1200 Time [min]

1

2

3 C-Rate

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

4

5 KIT, IAM-AWP

Summary and Outlook Combination of electrochemical-calorimetric studies on cylindrical lithium ion cells and thermal modelling by COMSOL Multiphysics software

Results:  Significant temperature rise even for low charge/discharge rates  Analysis of heat generation and dissipation during electrochemical cycling  Good agreement between experimental and simulation results Outlook:  Comparison of different cell chemistries  Assignment of heat effects to electrochemical processes and components

Improvement of Thermal Management and Safety 22

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP

Acknowledgement

Thank you for your attention! This R&D project is partially funded by the Federal Ministry for Education and Research (BMBF) within the framework „IKT 2020 Research for Innovations“ under the grant 16N12515 and is supervised by the Project Management Agency VDI│VDE│IT. supervised by

We gratefully acknowledge the financial support by the Helmholtz program STN and the funding in the Helmholtz Portfolio Project “Electrochemical Storage in the System”

23

Dr. C. Ziebert – 225th ECS Meeting 11.-15.05.2014, Orlando

KIT, IAM-AWP