Battery Pack Temperature Distribution Simulation with COMSOL and MATLAB

Transcription

1 Battery Pack Temperature Distribution Simulation with COMSOL and COMSOL CONFERENCE EUROPE M. Masomtob 1,2,4, C. Schaeper 1,2, W. Zhou 1 and D. U. Sauer 1,2,3 1 Electrochemical Energy Conversion and Storage Systems Group, Institute for Power Electronics and Electrical Drives (ISEA), RWTH Aachen University, Germany 2 Juelich Aachen Research Alliance, JARA-Energy, Germany 3 Institute for Power Generation and Storage Systems (PGS), E.ON ERC, RWTH Aachen University, Germany 4 National Metal and Materials Technology Center (MTEC), National Science and Technology Development Agency (NSTDA), Thailand batteries@isea.rwth-aachen.de Excerpt from the Proceedings of the 2012 COMSOL Conference in Milan

2 Why COMSOL and COMSOL MATLAB 3D CAD import FEM and CFD Monitoring Large System Integration Objective of this Work Combination between COMSOL and for a liquid cooling system of a battery pack 2

3 Electric Car Electric Car 3

4 Battery Pack 80 s Battery Pack z y x Inlet Outlet 26 Battery Cells 4

5 Diagram of the Cooling System First Floor SIDE A Fourth Floor Inlet Outlet Second Floor SIDE B Third Floor 5

6 Simplifying the Cooling System First Floor SIDE A Fourth Floor Inlet Outlet Second Floor SIDE B Third Floor Coolant inlet T=constant Series Cooling Coolant outlet 1 st BM 2 nd BM 3 r d BM 9 th BM 10 th BM Flow rate in each channel, 12.5 % 6

7 Simulation problem Coolant inlet T=constant Series Cooling Coolant outlet 1 st BM 2 nd BM 3 r d BM 9 th BM 10 th BM Q= f(t) T = constant value Heat generation in each battery module depends on electric load. 7

8 Techniques for Simulation 10 Battery Modules in COMSOL Only One at a time COMSOL Only COMSOL with MATLAB Time of mesh Generation Data Computer Performance - Medium Medium Not Required Medium Automatic High Medium Low Total runtime - Normal Faster 8

9 Soildworks 3D CAD Making 3D CAD from Simpify Parasolid.file 9

10 COMSOL Setting 3D CAD from Parasolid.file Material Property Boundary Condition Mesh Generation Module Laminar Flow Transient Study To test model in COMSOL m.file 10

11 MATLAB of modify 3D CAD from Manual Modification So that the from COMSOL can exchange parameters between m.files. - Normal Values - Values in form of txt.file, in case the values as a function of time m.file from COMSOL 11

12 MATLAB of modify How to import txt.file? model.func.create('int1', 'Interpolation'); model.func('int1').model('mod1'); model.func('int1').set('source', 'file'); model.func('int1').set('filename', 'D:\running model\version1\heat_loss.txt'); model.func('int1').setindex('funcs', 'Heat', 0, 0); model.func('int1').importdata; Time (s) How to use txt.file? Heat_Loss.txt : : : : model.physics('nitf').feature('hs1').set('ptot', 1, 'Heat(t)'); How to use normal parameter (Constant Value)? Heat Generation(w) model.physics('nitf').feature('inl1').set('u0in', 1, 'Flow/Area'); 12

13 MATLAB of modify 3D CAD from Manual Modification So that the from COMSOL can exchange parameters between m.files. - Normal Values - Values in form of txt.file, in case the values as a function of time m.file from COMSOL of 13

14 MATLAB of modify 3D CAD from of **One all s** 14

15 MATLAB for running management 3D CAD from Functions Start running each m.file of Battery Module in serial Set parameters both txt.file and normal values Collect and save the results to txt.file Calibrate the results of previous battery module before setting the inputs of next battery module Show runtime of each m.file of Battery Module 15

16 3D CAD from MATLAB Process BM-1 Input of First battery module Inlet temperature is constant, 298 K. Flow rate of coolant is 1 L/min Heat generation as function of time(txt.file). T = Ti 16

17 Temperature (K) Heat Gerneration per (W) 3D CAD from MATLAB Process BM T_Inlet_BM1 Q_Generation Input Flow Rate 1L/min Time (s) T = Ti 17

18 3D CAD from MATLAB Process BM-1 of First battery module Outlet temperature depends on heat generation of BM-1(txt.file) Temperatures of each battery cells (txt.file) Runtime ranges between 34 to 37 hours. T = Ti 18

19 Temperature (K) Heat Gerneration per (W) 3D CAD from MATLAB Finished Process BM T_Inlet_BM1 T_Outlet_BM1 Q_Generation Flow Rate 1L/min Time (s) T = Ti 19

20 3D CAD from MATLAB Process BM-2 Input of Second battery module Inlet temperature is outlet temperature of first battery module (txt.file). Flow rate remains constant at 1 L/min Heat generation as a function of time (txt.file). of Second battery module Outlet temperature depends on heat generation of BM-1 and BM-2 (txt.file) Temperatures of each battery cell (txt.file) T = Ti 20

21 Temperature (K) Heat Gerneration per (W) 3D CAD from MATLAB Finished Process BM T_Inlet_BM1 T_Outlet_BM1 T_Outlet_BM2 Q_Generation Flow Rate 1L/min Time (s) T = Ti 21

22 3D CAD from Module MATLAB Finished Process 10 BMs 299 T_Outlet_BM1 T_Outlet_BM5 T_Outlet_BM9 T_Outlet_BM2 T_Outlet_BM6 T_Outlet_BM10 T_Outlet_BM3 T_Outlet_BM7 Q_Generation Temperature (K) T_Inlet_BM1 T_Outlet_BM4 T_Outlet_BM Flow Rate 1L/min Heat Gerneration per (W) Time (s) T = Ti Manop Masomtob 22

23 Conclusions The m.file of the s from COMSOL is able to run with the to investigate the temperature distribution of the cooling system in the battery pack time of each battery module in COMSOL with MATLAB is faster than running time in COMSOL, 2 times ( Intel(R) Core (TM) i7 CPU 2.80 GHz 2.80 GHz and 6 GB of RAM) Furture work The m.file of from COMSOL will be connected with Simulink/MATLAB for Large system 23

24 Thank you for your attention 24

25 Battery Pack Temperature Distribution Simulation with COMSOL and COMSOL CONFERENCE EUROPE M. Masomtob 1,2,4, C. Schaeper 1,2, W. Zhou 1 and D. U. Sauer 1,2,3 1 Electrochemical Energy Conversion and Storage Systems Group, Institute for Power Electronics and Electrical Drives (ISEA), RWTH Aachen University, Germany 2 Juelich Aachen Research Alliance, JARA-Energy, Germany 3 Institute for Power Generation and Storage Systems (PGS), E.ON ERC, RWTH Aachen University, Germany 4 National Metal and Materials Technology Center (MTEC), National Science and Technology Development Agency (NSTDA), Thailand batteries@isea.rwth-aachen.de