Simulation of Thermal Behaviour of a Lithium Titanate Oxide Battery

One of the reasonable possibilities to investigate the battery behaviour under various temperature and current conditions is the development of a model of the lithium-ion batteries and then by employing the simulation technique to anticipate their behaviour. This method not only can save time but also they can predict the behaviour of the batteries through simulation. In this investigation, a three-dimensional model is developed to simulate thermal and electrochemical behaviour of a 13Ah lithium-ion battery. In addition, the temperature dependency of the battery cell parameters was considered in the model in order to investigate the influence of temperature on various parameters such as heat generation during battery cell operation. Maccor automated test system and isothermal battery calorimeter were used as experimental setup to validate the thermal model, which was able to predict the heat generation rate and temperature at different positions of the battery. The three-dimensional temperature distributions which were achieved from the modelling and experiment were in well agreement with each other throughout the entire of discharge cycling at different environmental temperatures and discharge rates.


Introduction
Lithium-ion batteries are one of the most developing categories of batteries on the market these days because of their high energy density and capacity. A large amount of energy is stored inside them and they have great sensitivity to the operating conditions. Therefore, safety is an important issue in lithium-ion batteries. In addition, demands on safety of these batteries is increasing with their utilization in more applications.
With the intention of reaching out to safety requirements of the lithium-ion batteries on electronic device applications, researchers are resuming to do supplementary investigations on the essential issues in relation to the lithium-ion batteries.
System safety, cycle life, and cell performance are influenced by temperature distribution in the cell. Consecutively, it depends on heat dissipation rate at surface of the cell and heat generation rate within the cell.
Although lithium-ion batteries are susceptible to extreme heat load under severe or abnormal functional conditions, thermal management has been one of the considerable issues in developing lithium-ion batteries in hybrid electric vehicle and battery system applications.
A pseudo 2D electrochemical model for modelling electrochemical systems subject to realistic automotive operation situations was proposed [1]. The model was developed for a lithium ion battery. It consists of complicated electrochemical phenomena, which were generally eliminated in online battery performance forecasters such as over potentials owing to mass transport restrictions and the A procedure was designed for dependency modelling of the discharge behaviour of a lithium-ion battery cell on the environmental temperature [12]. The two-dimensional modelling of the potential was validated by the modelling outcomes. The heat generation rates as a function of the position on the electrodes and the discharge time were determined in order to anticipate the temperature distributions of the lithium-ion battery [12]. This was according to the modelling outcomes of current density and the potential distributions. The temperature distributions, which were achieved from experimental measurements, were in good agreement with the modelling [12].
A battery cell model, which is flexible to investigate the thermal, electrochemical, physical phenomena, and advance over extensive length scales in battery cell systems of different assemblies, is necessary.
Unfortunately, thermal parameter measurement explanations and electrical parameters determination for lithium-ion batteries were not conveniently found in the literature. Many researchers commonly address the thermal parameters without reporting measurement procedures. A detailed description of thermal parameter measurement is reported in this investigation. Notwithstanding, to the author's best knowledge, only very few publications [13][14][15] are available in the literature that discuss the thermal simulation of lithium-ion batteries by considering all of the influential parameters such as thermal, electrical, and chemical processes on the thermal behaviour of the lithium-ion batteries. In addition, most of the previous studies did not take into account all of the electrochemical phenomena. In this investigation, the Multi-Scale Multi-Dimensional (MSMD) battery module was used for a lithium titanate oxide battery, which to the author´s best knowledge, it has not been done yet. The investigated model is able to determine the surface temperature distribution of the battery cell at various operating conditions with high accuracy.

The Battery Modelling
A 13 Ah pouch type commercial lithium-ion battery cell with dimensions of 204 mm width, 129 mm length, and 7.7 mm thickness and a lithium titanate oxide based anode was modelled for all simulations. The picture of the battery cell inside fixture, which was chosen for this investigation, is illustrated in Figure 1. temperature distributions of the lithium-ion battery [12]. This was according to the modelling outcomes of current density and the potential distributions. The temperature distributions, which were achieved from experimental measurements, were in good agreement with the modelling [12]. A battery cell model, which is flexible to investigate the thermal, electrochemical, physical phenomena, and advance over extensive length scales in battery cell systems of different assemblies, is necessary.
Unfortunately, thermal parameter measurement explanations and electrical parameters determination for lithium-ion batteries were not conveniently found in the literature. Many researchers commonly address the thermal parameters without reporting measurement procedures. A detailed description of thermal parameter measurement is reported in this investigation. Notwithstanding, to the author's best knowledge, only very few publications [13][14][15] are available in the literature that discuss the thermal simulation of lithium-ion batteries by considering all of the influential parameters such as thermal, electrical, and chemical processes on the thermal behaviour of the lithium-ion batteries. In addition, most of the previous studies did not take into account all of the electrochemical phenomena. In this investigation, the Multi-Scale Multi-Dimensional (MSMD) battery module was used for a lithium titanate oxide battery, which to the author´s best knowledge, it has not been done yet. The investigated model is able to determine the surface temperature distribution of the battery cell at various operating conditions with high accuracy.

The Battery modelling
A 13 Ah pouch type commercial lithium-ion battery cell with dimensions of 204 mm width, 129 mm length, and 7.7 mm thickness and a lithium titanate oxide based anode was modelled for all simulations. The picture of the battery cell inside fixture, which was chosen for this investigation, is illustrated in Figure 1.
In accordance with the construction and geometry of the battery cell, a three-dimensional model was constructed in ANSYS (2018) and the battery geometry was generated and analysed in an appropriate manner for additional analysis. Different components of positive and negative current tabs are illustrated in Figure 2. The negative tab and positive tab are used to accumulate the current flow via the battery cell. As shown in the figure five parts comprise the model: (1) Positive current tab meshing structure; (2) Interior part of positive current tab; (3) Positive current tab; (4) Skin of positive current tab; (5) Contact region of positive current tab.  In accordance with the construction and geometry of the battery cell, a three-dimensional model was constructed in ANSYS (2018) and the battery geometry was generated and analysed in an appropriate manner for additional analysis. Different components of positive and negative current tabs are illustrated in Figure 2. The negative tab and positive tab are used to accumulate the current flow via the battery cell. As shown in the figure five parts comprise the model: The active volume, contact region, and skin of the battery cell are illustrated in Figure 3. The active volume demonstrates the stacked construction, comprising separator layers, negative and positive active materials, and aluminium foils. The thin skin enclosing the active volume. Geometrical structured meshing of lithium titanate oxide battery cell is shown in Figure 4. As shown in Figure 3 four parts comprise the model:   The active volume, contact region, and skin of the battery cell are illustrated in Figure 3. The active volume demonstrates the stacked construction, comprising separator layers, negative and positive active materials, and aluminium foils. The thin skin enclosing the active volume. Geometrical structured meshing of lithium titanate oxide battery cell is shown in Figure 4. As shown in Figure 3 four parts comprise the model:  The active volume, contact region, and skin of the battery cell are illustrated in Figure 3. The active volume demonstrates the stacked construction, comprising separator layers, negative and positive active materials, and aluminium foils. The thin skin enclosing the active volume. Geometrical structured meshing of lithium titanate oxide battery cell is shown in Figure 4. As shown in Figure 3 four parts comprise the model:     The active volume, contact region, and skin of the battery cell are illustrated in Figure 3. The active volume demonstrates the stacked construction, comprising separator layers, negative and positive active materials, and aluminium foils. The thin skin enclosing the active volume. Geometrical structured meshing of lithium titanate oxide battery cell is shown in Figure 4. As shown in Figure 3 four parts comprise the model:

Determination of Thermal Parameters
The general energy balance differential equation, which explains the distribution of generated local, conducted, and accumulated heats and the variation of temperature within a battery cell or pack, can be written as follows: where C p : Specific heat capacity λ: Volumetric heat generation ρ: Physical mass k i : Thermal conduction in direction i T: Temperature Free convection is the main heat transfer process from the surfaces of the battery cell. The dissipated heat flux from the battery surface to the surrounding can be considered by both the convection and the radiation heat contributions: where σ: Stefan-Boltzmann constant ε: Emissivity of the battery cell surface T a : Ambient temperature h: Convective heat transfer T s : Battery surface temperature Natural convection on a surface is a function of the orientation besides the geometry of the surface. In addition, it depends on the thermos physical properties of the fluid and the variation of temperature on the surface [16]. The complicatedness of the fluid flow makes it hard to achieve straightforward analytical relations for natural convection. The Rayleigh number, which controls the flow regime in natural convection, is defined as the product of the Prandtl and Grashof numbers [16]: where v: Kinematics viscosity of the fluid β: Coefficient of volume expansion g: Gravitational acceleration δ: Characteristic length of the geometry The Nusselt number for natural convection could be determined in the following form [16]: where the constants n and C depend on the flow and the geometry of the surface. In this investigation natural convection process was considered for heat transfer from the surfaces of the battery cell. A procedure was used for the determination of the thermal parameters of the lithium-ion batteries. The method is able to determine a single thermal parameter such as specific heat capacity or thermal conductivity.
The heat capacity is quantifiable physical parameters of a substance, which characterizes the amount of heat, which is required to alter the temperature of the substance by one degree. In other words, heat capacity or specific heat is the amount of heat, which is needed to heat or cool 1 kg of a material by 1 • C. It determines how quickly a battery heats or cools down from its primary temperature in a given surrounding condition. The outcome of the heat capacity characterization process depends on the amount of absorbed heat, which consecutively depends on the ambient circumstances of the process.
In this investigation, isothermal battery calorimeter was employed for determination of specific heat. In order to measure the heat capacity of the battery cell with known mass (m) a thermal procedure was selected. At first, the battery was placed in the Maccor chamber for 3 h to be equilibrated at initial temperature T 1 . The chamber temperature was set to T 1 . Then it was placed rapidly in the isothermal battery calorimeter. The isothermal battery calorimeter temperature was set to T 2 . After a period, the battery cell temperature reached to the chamber temperature. This procedure was repeated several times and the average amount was considered. The heat (Q) which was transferred between the isothermal battery calorimeter and battery cell was measured by using heat flux sensors inside chamber. Consequently, by having the temperature difference (T 1 − T 2 ) and the mass (m) of the battery cell the heat capacity was calculated by using the following equation: Transient and steady state methods are the main methods, which could be used for the determination of the thermal conductivity of batteries. Guarded hot plate is a steady state method, which can approximately estimate the thermal conductivity. In this method, the battery is placed between a heat sink and heat source. The thermal conductivity could be determined by knowing the battery cell thickness, temperature difference across the battery, and heat flux.
where q: Heat flux ∆T: Temperature gradient ∆x: Thickness To assure of the accuracy of the results different ways were used for determining the battery cell density. The battery density was calculated based on the battery material. In addition, the density of the lithium-ion battery was calculated by measuring the mass and volume of the battery cell. In order to determine precisely the physical and thermal parameters of the lithium titanate oxide battery cell and to assure of the accuracy of the previous estimations a procedure was selected. In this method, the battery cell was divided to different parts such as negative current collector, negative electrode, separator, positive electrode and positive current collector. The battery cell cross-section along with the thickness of different layers is illustrated in Figure 5. The material properties of the battery cell were estimated by employing the following formulations [17]: Energies 2019, 12, 679 where α: Thermal diffusivity T N c : Thickness of negative current collector T N e : Thickness of negative electrode T S : Thickness of separator T P e : Thickness of positive electrode T P c : Thickness of positive current collector K P c : Thermal conductivity of positive current collector K P e : Thermal conductivity of positive electrode K S : Thermal conductivity of separator K N e : Thermal conductivity of negative electrode K N c : Thermal conductivity of negative current collector ρ P c : Density of positive current collector ρ P e : Density of positive electrode ρ S : Density of separator ρ N e : Density of negative electrode ρ N c : Density of negative current collector C P P c : Heat capacity of positive current collector C P P e : Heat capacity of positive electrode C PS : Heat capacity of separator C P N e : Heat capacity of negative electrode C P N c : Heat capacity of negative current collector σ P c : Electric conductivity of positive current collector σ P e : Electric conductivity of positive electrode σ N c : Electric conductivity of negative current collector σ N e : Electric conductivity of negative electrode

Determination of Electrical Parameters
A 2RC equivalent circuit model was used in this investigation. The model is presented in Figure  6. To achieve the parameters of the equivalent circuit model different loading profile were applied to the battery cell. The loads consist of charge and discharge cycles with different C-rates. The voltage variation of load profiles at different temperature and C-rates is illustrated in Figure 7. An example of the voltage response when a 4 C-rate (52 Ah) current is applied to the battery cell at 27 °C is

Determination of Electrical Parameters
A 2RC equivalent circuit model was used in this investigation. The model is presented in Figure 6. To achieve the parameters of the equivalent circuit model different loading profile were applied to the battery cell. The loads consist of charge and discharge cycles with different C-rates. The voltage variation of load profiles at different temperature and C-rates is illustrated in Figure 7. An example of the voltage response when a 4 C-rate (52 Ah) current is applied to the battery cell at 27 • C is illustrated in Figure 8.

Determination of Electrical Parameters
A 2RC equivalent circuit model was used in this investigation. The model is presented in Figure  6. To achieve the parameters of the equivalent circuit model different loading profile were applied to the battery cell. The loads consist of charge and discharge cycles with different C-rates. The voltage variation of load profiles at different temperature and C-rates is illustrated in Figure 7. An example of the voltage response when a 4 C-rate (52 Ah) current is applied to the battery cell at 27 °C is illustrated in Figure 8.

Determination of Electrical Parameters
A 2RC equivalent circuit model was used in this investigation. The model is presented in Figure  6. To achieve the parameters of the equivalent circuit model different loading profile were applied to the battery cell. The loads consist of charge and discharge cycles with different C-rates. The voltage variation of load profiles at different temperature and C-rates is illustrated in Figure 7. An example of the voltage response when a 4 C-rate (52 Ah) current is applied to the battery cell at 27 °C is illustrated in Figure 8.  The values of open circuit voltage, resistances, and capacitances of the 2RC equivalent circuit model are illustrated in Figure 9. The parameters of the 2RC equivalent circuit model were determined by using the following equations [18,19]: The values of open circuit voltage, resistances, and capacitances of the 2RC equivalent circuit model are illustrated in Figure 9. The parameters of the 2RC equivalent circuit model were determined by using the following equations [18,19]: Energies 2019, 12, 679 9 of 15 where V 0 : Battery voltage before the discharge current pulse is applied. R 0 : Ohmic resistance I: Amplitude of the current pulse. V 1 : Battery voltage one seconds (t 1 ) after the discharge current is applied. R 1 , C 1 : Resistance and capacitance of the first RC network V 2 : Battery voltage ten seconds (t 2 ) after the discharge current is applied. V 3 : Battery voltage eighteen seconds (t 3 ) after the discharge current is applied. R 2 , C 2 : Resistance and capacitance of the second RC network determined by using the following equations [18,19]: where V0: Battery voltage before the discharge current pulse is applied. R0: Ohmic resistance I: Amplitude of the current pulse. V1: Battery voltage one seconds (t1) after the discharge current is applied. R1, C1: Resistance and capacitance of the first RC network V2: Battery voltage ten seconds (t2) after the discharge current is applied. V3: Battery voltage eighteen seconds (t3) after the discharge current is applied.  Table 1. These functions were expressed in fifth order polynomial form:   The coefficients of fifth order polynomial are illustrated in Table 1. These functions were expressed in fifth order polynomial form:

Modelling Method
The model was solved in ANSYS by employing the Multi-Scale Multi-Dimensional (MSMD) battery module. The model combines the principal design parameters of the battery cell such as corresponding physical parameters, materials, and dimensions to computational fluid dynamics and heat transfer. In addition, the battery model is able to simulate a single battery cell or a battery pack to investigate their electrochemical and thermal behaviour. At the solution phase of the model, unsteady state problem and the thermal time interdependent were solved numerically assuming the heat generation in the battery cell as a dynamic source. The amount of heat generation inside a lithium-ion battery cell, which was proportionate to the temperature and current rate, was measured by an isothermal battery calorimeter and was considered as an input to the thermal model. The thermal and electrical fields were solved by using the following equations [17]:

Measurement of Heat Generation Rates
Heat generation inside batteries is a complicated process, which could be divided into reversible and irreversible parts. Experimental setup, which was used in this investigation, is illustrated in Figure 10. Isothermal battery calorimeter was employed for heat loss measurement. The detailed experimental procedures and setup, as well as corresponding equipment and materials, could be found in [20]. Maccor automated system was used as a battery cycler for the whole experiments. The battery cycler charged and discharged the battery with different current rates. Heat flux determination of lithium titanate oxide battery cell by using isothermal calorimeter is shown in Figure 11. determination of lithium titanate oxide battery cell by using isothermal calorimeter is shown in Figure 11.

Experimental Validation
Several sets of experiments were accomplished to validate model performance at different temperature by using a 13 Ah battery, which was fabricated by Altairnano and is shown in Figure 4. A 45A discharge process was used to validate the temperature, which was anticipated by the model. FLIR thermal camera and contact thermocouples were used to monitor the surface temperature of the battery cell at different positions. Experimental temperature results and temperature simulation for 45 A discharge process is illustrated in Figure 12. It could be seen from the figure that the temperature increase of the simulation is in good agreement with the experimental data. This demonstrate that the model is capable of simulating the real battery cell. The experimental temperature increase data, which was used in this investigation, was the average amount of several thermocouples and FLIR thermal camera measurements. In addition, the simulation temperature increase was the average value of entire battery cell temperature. determination of lithium titanate oxide battery cell by using isothermal calorimeter is shown in Figure 11.

Experimental Validation
Several sets of experiments were accomplished to validate model performance at different temperature by using a 13 Ah battery, which was fabricated by Altairnano and is shown in Figure 4. A 45A discharge process was used to validate the temperature, which was anticipated by the model. FLIR thermal camera and contact thermocouples were used to monitor the surface temperature of the battery cell at different positions. Experimental temperature results and temperature simulation for 45 A discharge process is illustrated in Figure 12. It could be seen from the figure that the temperature increase of the simulation is in good agreement with the experimental data. This demonstrate that the model is capable of simulating the real battery cell. The experimental temperature increase data, which was used in this investigation, was the average amount of several thermocouples and FLIR thermal camera measurements. In addition, the simulation temperature increase was the average value of entire battery cell temperature.  Figure 11. A heat generation analysis of lithium titanate oxide battery cell for 300 K.

Experimental Validation
Several sets of experiments were accomplished to validate model performance at different temperature by using a 13 Ah battery, which was fabricated by Altairnano and is shown in Figure 4. A 45A discharge process was used to validate the temperature, which was anticipated by the model. FLIR thermal camera and contact thermocouples were used to monitor the surface temperature of the battery cell at different positions. Experimental temperature results and temperature simulation for 45 A discharge process is illustrated in Figure 12. It could be seen from the figure that the temperature increase of the simulation is in good agreement with the experimental data. This demonstrate that the model is capable of simulating the real battery cell. The experimental temperature increase data, which was used in this investigation, was the average amount of several thermocouples and FLIR thermal camera measurements. In addition, the simulation temperature increase was the average value of entire battery cell temperature. During the charge and discharge cycling, the battery cell was made to be isothermal by situating it in an isothermal battery calorimeter. The calorimetric measurements were used as a heat generation source in the battery cell. A Maccor automated test system was employed as the discharge apparatus with the intention of monitoring the current and voltage. The simulated and experimental outcomes for 52 A discharge process and 300 K environmental temperature are illustrated in Figure 13. The experimental temperature data, which were used for model validation, are the value of surface temperature of the battery at four different locations, which were measured by four contact thermocouples. Temperature value of the experimental data compare good to the simulation, demonstrating that the model could simulate the real battery cell.

Simulation Results and Discussion
In this investigation, several quantities for the volumetric heat generation were considered. The values were measured by an isothermal battery calorimeter for various load profiles. The model is simulated in both time dependent and steady state environment to determine the temperature spatial distribution over the battery surface. In addition, the modelling is able to show the maximum value of surface temperature of the battery as a function of time for different environment temperature and discharge current rates. During the charge and discharge cycling, the battery cell was made to be isothermal by situating it in an isothermal battery calorimeter. The calorimetric measurements were used as a heat generation source in the battery cell. A Maccor automated test system was employed as the discharge apparatus with the intention of monitoring the current and voltage. The simulated and experimental outcomes for 52 A discharge process and 300 K environmental temperature are illustrated in Figure 13. The experimental temperature data, which were used for model validation, are the value of surface temperature of the battery at four different locations, which were measured by four contact thermocouples. Temperature value of the experimental data compare good to the simulation, demonstrating that the model could simulate the real battery cell. During the charge and discharge cycling, the battery cell was made to be isothermal by situating it in an isothermal battery calorimeter. The calorimetric measurements were used as a heat generation source in the battery cell. A Maccor automated test system was employed as the discharge apparatus with the intention of monitoring the current and voltage. The simulated and experimental outcomes for 52 A discharge process and 300 K environmental temperature are illustrated in Figure 13. The experimental temperature data, which were used for model validation, are the value of surface temperature of the battery at four different locations, which were measured by four contact thermocouples. Temperature value of the experimental data compare good to the simulation, demonstrating that the model could simulate the real battery cell.

Simulation Results and Discussion
In this investigation, several quantities for the volumetric heat generation were considered. The values were measured by an isothermal battery calorimeter for various load profiles. The model is simulated in both time dependent and steady state environment to determine the temperature spatial distribution over the battery surface. In addition, the modelling is able to show the maximum value of surface temperature of the battery as a function of time for different environment temperature and discharge current rates.

Simulation Results and Discussion
In this investigation, several quantities for the volumetric heat generation were considered. The values were measured by an isothermal battery calorimeter for various load profiles. The model is simulated in both time dependent and steady state environment to determine the temperature spatial distribution over the battery surface. In addition, the modelling is able to show the maximum value of surface temperature of the battery as a function of time for different environment temperature and discharge current rates.
The temperature distributions of the battery cell were determined as a function of time at different discharge rates. As could be seen in Figures 12 and 13, the overall temperature distributions, which were achieved from the model and experiment, are in good agreement with each other.
The heat, which is dissipated from the battery, and the heat generated inside it are approximately equal during low current rates. Therefore, fast equilibrium could be reached. In another word, most of generated heat will be transferred to the surrounding through free or forced convection. The evolution of the uttermost temperature of the battery cell is confined at low current rates. The phenomenon demonstrates the minor rises in the surface temperature. Notwithstanding, the difference among the maximum temperatures, which were attained, from the modelling and experiment was lower than that between the corresponding minimum ones.
The simulation was accomplished at different discharge current rates ranging from 0.25 C to 9 C with 0.25 C interval. The modelling discharge profiles agree good with those, which were gained from experimental. The corresponding heat loss from the battery cell was shown in Figure 11. In accordance with the findings, the temperature increase sharply to a specific point. As anticipated, the position of the hottest area is seen near the negative tab of the battery cell throughout the discharge process. In addition, non-uniform temperature propagation was observed.
Temperature distribution of the battery cell at different discharge rates ranging from 0.25 C to 9 C is illustrated in Figure 14. As could be seen from the distributions of temperature over the volume of the battery cell, the temperature contours showed moderate slopes at low current rates. On the contrary, the temperature contours demonstrate sharper slopes at higher current rates. It demonstrates quicker temperature increase during the discharge of the battery cell, which is due to the higher heat generation inside the battery cell. At high current rates, achieving equilibrium occurs in a longer time. This phenomenon demonstrates that the modelling discharge curves agree well with those which were achieved from the experiments. The temperature distributions of the battery cell were determined as a function of time at different discharge rates. As could be seen in Figures 12 and 13, the overall temperature distributions, which were achieved from the model and experiment, are in good agreement with each other.
The heat, which is dissipated from the battery, and the heat generated inside it are approximately equal during low current rates. Therefore, fast equilibrium could be reached. In another word, most of generated heat will be transferred to the surrounding through free or forced convection. The evolution of the uttermost temperature of the battery cell is confined at low current rates. The phenomenon demonstrates the minor rises in the surface temperature. Notwithstanding, the difference among the maximum temperatures, which were attained, from the modelling and experiment was lower than that between the corresponding minimum ones.
The simulation was accomplished at different discharge current rates ranging from 0.25 C to 9 C with 0.25 C interval. The modelling discharge profiles agree good with those, which were gained from experimental. The corresponding heat loss from the battery cell was shown in Figure 11. In accordance with the findings, the temperature increase sharply to a specific point. As anticipated, the position of the hottest area is seen near the negative tab of the battery cell throughout the discharge process. In addition, non-uniform temperature propagation was observed.
Temperature distribution of the battery cell at different discharge rates ranging from 0.25 C to 9 C is illustrated in Figure 14. As could be seen from the distributions of temperature over the volume of the battery cell, the temperature contours showed moderate slopes at low current rates. On the contrary, the temperature contours demonstrate sharper slopes at higher current rates. It demonstrates quicker temperature increase during the discharge of the battery cell, which is due to the higher heat generation inside the battery cell. At high current rates, achieving equilibrium occurs in a longer time. This phenomenon demonstrates that the modelling discharge curves agree well with those which were achieved from the experiments.  electrode is much higher than that of the positive electrode. This phenomenon leads to lower temperatures in the neighbourhood of the current collecting tab of the positive electrode compared to the negative electrode [21][22][23].
Electrochemical reaction rate was increased due to higher temperature gradient. This phenomenon could be described by higher current rates in some parts of the cell owing to high temperature gradient. The minimum and maximum temperatures, which, were collected from the modelling and experiment, are in good agreement with each other over the entire scope of battery cell surface at the different discharge rates. Notwithstanding, some difference was seen among the discharge curves, which were achieved from the experiment and model in close proximity to the end period of discharge. The highest discrepancy was seen for high discharge rates.

Conclusions
The principal objective of this investigation was to develop a precise, computationally efficient and simplified Dual Potential Multi-Scale Multi-Dimensional (MSMD) Battery model. A procedure was used to simulate the thermal behaviour of a lithium-ion battery at different current rates and environmental temperature. The three-dimensional temperature distribution of the battery was anticipated as a function of the discharge time by using the model. By comparing the modelling discharge curves with the experimental outcomes at different environmental temperatures and discharge rates the modelling was validated. The parameters of equivalent circuit model were determined from multi-pulse charge and discharge data. An average specific heat capacity was considered for the battery cell in the time dependent and unsteady state simulation. To assist the suggested modelling method, calorimetric experiments were accomplished. By using the heat generation, which was measured by the isothermal calorimeter, the model was simulated to demonstrate the temperature distribution. A great temperature discrepancy was seen in battery surface at high current rates. This phenomenon could be described by high amount of heat generation due to higher temperature gradient. In addition, the value of surface temperature of the battery was determined by the model, was compared to the experimental data, and was in good agreement with the model data. Greater temperature gradients were seen at the battery cell surfaces owing to higher current rate. Subsequently, design of an appropriate thermal management system specifically during high current rates charging and discharging could play a fundamental role in preventing great temperature growth of the battery cell. The simulation methodology, which was demonstrated in this investigation, might contribute to the development of a battery cell thermal management system, which enables the temperature evolution of lithium-ion batteries as a subordinate of time to be more precisely anticipated. In addition, it can assist to anticipate the evolution of the thermal, electrical, and chemical processes.