Next Article in Journal
Economic, Energy and Environmental Efficiency of Road Freight Transportation Sector in the EU
Next Article in Special Issue
Joint Prediction of the State of Charge and the State of Health of Lithium-Ion Batteries Based on the PSO-XGBoost Algorithm
Previous Article in Journal
Energy Efficiency of Heavy-Duty Vehicles in Mexico
Previous Article in Special Issue
Study on the Homogeneity of Large-Size Blade Lithium-Ion Batteries Based on Thermoelectric Coupling Model Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cooling Optimization Strategy for a 6s4p Lithium-Ion Battery Pack Based on Triple-Step Nonlinear Method

School of Energy and Power Engineering, Huazhong University of Science & Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(1), 460; https://doi.org/10.3390/en16010460
Submission received: 30 November 2022 / Revised: 22 December 2022 / Accepted: 28 December 2022 / Published: 31 December 2022
(This article belongs to the Special Issue Advanced Application Technology of Lithium-Ion Batteries)

Abstract

:
In a battery cooling system, by adopting a cooling optimization control strategy, the battery temperature under different external environments and load currents can be adjusted to ensure performance and safety. In this study, two modes of the thermal management system are established for the 6s4p (six serial and four parallel batteries in a stage) battery pack. A single particle model, considering battery aging, is adopted for the battery. Furthermore, a cooling optimization control strategy for the battery is proposed based on the triple-step nonlinear method, and then the optimization effect is validated under two C-rate charge–discharge cycles, NEDC cycles, and US06 cycles. Moreover, an extended PID control strategy is constructed and compared with the triple-step nonlinear method. A comparison of pump power, thermal behavior, and aging performance indicate parallel cooling is more advantageous. This verifies the validity of the triple-step nonlinear method and shows its advantages over the extended PID method. The present study provides a method to investigate the thermal behavior and aging performance of a battery pack in a BTM system, and fills in the research gaps in the cooling optimization control strategy for battery packs.

1. Introduction

The vigorous development and widespread introduction of electric vehicles (EVs) have placed higher requirements on the driving endurance and safety performance of EVs [1]. Lithium-ion batteries (LIBs) are widely used in the automobile industry and in energy storage due to their high energy density [2,3]. However, the size and number of batteries cannot be increased indefinitely due to the limitations in vehicle space. Therefore, increasing the energy density and power density of batteries has become a priority to improve the driving endurance of EVs [4]. As the energy density and power density of the battery increase, the heat generation also increases. This requires an effective cooling method to dissipate the large amount of heat accumulated inside the battery over time [5]. Otherwise, the battery temperature will rise sharply and the battery capacity will continuously decline. In severe cases, it may even lead to thermal runaway or fire and explosion [6,7]. In addition, there are also many studies on the generation of LIBs for systematic research [8,9].
Many investigations have been conducted to better understand and improve battery performance. Some researchers have established battery thermal management (BTM) systems to control the temperature [10,11,12,13]. Studies of BTM system include air-cooling [14,15], liquid-cooling [16,17], phase-change material cooling [18,19], heat pipe-cooling [20,21], and cooling systems combined with the aforementioned methods [22,23,24]. Other researchers have studied the internal aging mechanism and battery aging to lay a foundation for the design of BTM systems [25,26]. Birkl et al. [27] summarized the causes, mechanisms, modes, and effects of battery aging. The causes trigger the mechanisms, and the mechanisms contribute to the modes of aging. The aging modes can be divided into three categories: loss of lithium inventory, loss of active material of the negative electrode (NE), loss of active material of the and positive electrode (PE). These modes contribute to the aging of the battery and the battery capacities decline. A widely studied aging mechanism is the formation of a solid electrolyte interface (SEI) [28,29]. The charge–discharge rate (C-rate, defined as the measurement of the charge and discharge speed, which equals the current required for the battery to charge or discharge its capacity in a specified time) and ambient temperature are generally considered to be the important factors affecting the aging mechanisms [30,31]. Many experimental studies on battery aging and thermal management systems have been conducted. However, experiments on thermal management systems are complicated and time-consuming, especially when multiple charge–discharge cycle tests are required to consider battery aging. Numerical simulations of BTM systems under different EV-related cycle conditions are efficient and low-cost as they can reduce the number of experiments [32].
The difficulty in the simulation of LIBs is related to their multi-domain and multi-physics characteristics. Therefore, various battery models have been used to characterize the performance of LIBs, and these models have been developed to study electrochemical performance [33], aging [34], and even thermal runaway [35]. The single particle (SP) model is a commonly-used reduced-order electrochemical model, which was proposed by Haran et al. [36] based on the pseudo-two-dimensional model of Doyle–Fuller–Newman [37]. In the SP model, it is assumed that both electrodes are composed of multiple spherical particles of uniform size and the current is uniformly distributed on both electrodes. Therefore, each electrode can be approximated as a single spherical particle. The SP model is described by a set of ordinary differential equations, and it is derived directly from a comprehensive electrochemical model, thus explicitly preserving many important battery properties with high computational efficiency. The disadvantage of the SP model is that the electrolyte physics and degradation are not considered, so the accuracy at high C-rates is affected. Many researchers have established advanced SP models to address this deficiency.
Rahimian et al. [38] compared the SP model with a semi-empirical equivalent circuit model to predict battery voltages under constant current charge and discharge for different LIB datasets. The results showed that the SP model outperforms the equivalent circuit model in its predictive power for all datasets, while the computation time of both models was on the same order of magnitude. Li et al. [39] developed an aging model based on the SP model considering SEI layer formation and the stress due to the volume expansion of the particles in the active material. The Li-ion loss model due to SEI layer formation was integrated with an advanced SP model that incorporates electrolytic physics. This model can rapidly predict capacity fade and voltage profile changes as a function of cycle number and temperature with high accuracy. Richardson et al. [40] compared the SP model and the Doyle–Fuller–Newman model for different electrode materials, including lithium nickel manganese cobalt oxide, graphite, and lithium iron phosphate (LFP). Motivated by these differences, correction terms for the SP model were derived. Incorporating this into the SP model resulted in a corrected SP model that was significantly more accurate than the original SP model. Ren et al. [41] proposed an extended SP model with higher accuracy. Based on this model, a new state of charge (SOC) closed-loop estimation algorithm model was constructed. The results showed that the extended SP model can effectively simulate the performance of the battery, and the closed-loop SOC estimation algorithm can correct the initial SOC error without increasing the computational complexity.
In the modeling of the BTM system, in addition to the selection and optimization of a battery model to ensure the accuracy of the simulation, it is also necessary to select an appropriate cooling method. As mentioned above, the current main cooling methods for batteries are air cooling and liquid cooling. However, the research on battery cooling mainly focuses on the structural design and analysis of factors, and research on the optimization and control strategy of battery cooling is still in its infancy. Ji et al. [42] proposed an active temperature control method to equalize the temperature distribution of the batteries by adjusting the coolant flow rate. The optimal cooling flow rate was determined by setting different coolant flow rates to find the lowest battery temperature. This method is complex and costly and is not suitable for complicated operating conditions. Zhu et al. [43] presented a finite-set-based model predictive control strategy for the BTM in hybrid EVs and verified the strategy under different cycle conditions. However, due to a large amount of computational data and high hardware requirements, this optimization process using dynamic programming algorithms was complicated. Ma et al. [44] proposed a cooling optimization strategy for LIBs based on the triple-step nonlinear method and compared the strategy with the PID method under NEDC and US06 cycle conditions. The results showed that the proposed method ensures a stable operating temperature and cooling process of LIBs. Ma et al. [45] further proposed a nonlinear model predictive control to optimize the cooling process of the battery module. The simulation results showed that the nonlinear model predictive control method ensured that the battery works near a target temperature under different working conditions. However, the battery model used in these studies was a lumped thermal model, which cannot reflect the electrochemical performance and aging of the battery. Moreover, the proposed cooling optimization strategies have only been used for the verification of a single battery, and the effect of applying a battery pack has not been studied.
As mentioned above, the advanced SP model is rarely used in the modeling of a BTM system and the cooling optimization strategy has been verified only for a single battery, which limits the application in real scenarios. Therefore, the gaps and deficiencies that are present in the literature are investigated in this study. First, a one-dimensional thermal management model of the 6s4p (six serial and four parallel batteries in a stage) battery pack is established for the first time based on our previous work [32]. The SP model, with a consideration of battery aging, is adopted to the 6s4p battery pack, so the electrochemical performance and degradation of each cell can be obtained more accurately. Secondly, two modes of the thermal management system are developed for the 6s4p battery pack, namely, a series cooling circuit and a parallel cooling circuit. The thermal behavior and aging performance of the battery pack under these two cooling modes are analyzed and compared. Thirdly, a cooling optimization control strategy for the key battery in the 6s4p battery pack is first proposed based on the triple-step nonlinear method, and then the optimization effect of this strategy in the battery pack is validated under two C-rate charge–discharge cycles, NEDC cycles, and US06 cycles. Finally, an extended PID control strategy is constructed and compared with the triple-step nonlinear method to verify the effects of and differences between the two strategies. The present study provides a method to investigate the thermal behavior and aging performance of the battery pack in a BTM system, and fills in the research gaps in cooling optimization for battery packs. The results in the present study demonstrate application of system modeling and the optimization of an automotive battery pack, and are applicable to the design of BTM systems.

2. Materials and Methods

2.1. Charge–Discharge and Nominal Electrochemical Model of the Battery Cell

The characteristics and the geometry of the LFP battery have been discussed in our previous work [29] and are listed in Table 1. The electrochemical reactions can be expressed as in [46].
The reaction at the positive electrode is:
Li 1 i FeP O 4 + i Li + + i e   charge discharge   Li FeP O 4
The reaction at the negative electrode is:
Li j C 6   charge discharge   Li 1 j C 6 + j Li + + j e
Figure 1a shows the charge–discharge model of the battery, including the battery model, thermal boundary condition, electrical connection, and cycle signal of the charge–discharge process. The thermophysical parameters are described by the material parameters module (listed in Table 1) and the part of the battery in contact with the air is described by the ambient conditions module, as shown in Figure 1a. The equation can be expressed as:
Q tra = h conv A b ( T b T amb ) + σ B χ ( T b 4 T amb 4 )
where Q tra is the heat transferred from the battery to the ambient air. h conv is the convective heat transfer coefficient and χ is the surface emissivity, and they are set as 10 W · m 2 · K 1 [32] and 0.8 [47], respectively. A b is the area between the battery and the ambient air. σ B is the Stefan–Boltzmann constant. T b and T amb are the battery temperature and ambient temperature, respectively. The electrochemical model chosen is a simplification of the pseudo-two-dimensional model (proposed by Doyle–Fuller–Newman [37] and described in detail in our previous work [48]) called the single particle (SP) model where the solid phase is modeled as a single particle, as shown in Figure 1b. This model has been widely used in the one-dimensional simulation software, AMESim (Advanced Modeling Environment for performing Simulation of engineering systems). The thickness direction contains five parts: a negative current collector (NCC), a negative electrode (NE), a separator (SEP), a positive electrode (PE), and a positive current collector (PCC), as shown in Figure 1b. In the SP model, the electrode is regarded as the superposition of two continuums, namely the electrolytic phase and a solid matrix, as shown in Figure 1b. The models of solid matrices are microscopic spherical particles in which Li-ions diffuse and react on the surface of the spheres. The main hypotheses [49] of this model are:
(1)
The passive sign convention is used: I > 0 during charge and I < 0 during discharge;
(2)
The particles in both electrodes can be represented as two mean spherical particles;
(3)
The active surface is assumed to be same in all areas: S p = S sep = S n = A ele .
The SP model computes the battery voltage based on its internal state and external solicitations. The voltage is the difference between the potential of its electrodes summed with the electrolyte potential [50]:
U ( t ) = ϕ s ( L ) ϕ s ( 0 ) = U p U n + η ( L ) η ( 0 ) + ϕ e ( L ) ϕ e ( 0 ) ,
where U p and U n are the thermodynamic potential of the PE and NE, respectively.   η ( L ) and η ( 0 ) are the overpotentials at z = L and z = 0, respectively. z is the battery thickness direction, as shown in Figure 1b. ϕ e ( L ) and ϕ e ( 0 ) are the electrolyte potentials at z = L and z = 0, respectively. Thermodynamic potentials are evaluated from the solid concentration of Li in the bulk (R = 0) of the single particle:
U p = U p ref ( x p ) = U p ref ( c s , p b c s , p max ) ,   U n = U n ref ( x n ) = U n ref ( c s , n b c s , n max ) ,
where c s , p b c s , p max and c s , n b c s , n max correspond to the coefficients i and j used in Equations (1) and (2), respectively. c s , p b and c s , n b are the bulk Li concentrations in the solid phase for PE and NE, respectively. c s , p max and c s , n max are the maximum Li concentrations in the PE and NE, respectively. U p ref and U n ref are the reference potentials of the PE and NE, respectively. The nominal physical and chemical phenomena occurring in Li-ion systems can be expressed by the mass conservation of Li+ species, charge conservation, and electrochemical kinetics, as shown in Table 2. All these equations are coupled into the framework of the P2D electrochemical model and have been discussed in previous works [45,48].
The diffusion of the solid phase leads to solid phase overpotential, which can be approximated as the difference between the bulk thermodynamic potential and the thermodynamic potential at the surface of the particles [50]:
U p = U p ref ( x p ) = U p ref ( c s , p b c s , p max ) ,   U n = U n ref ( x n ) = U n ref ( c s , n b c s , n max ) ,
where c s , p R s , p and c s , n R s , n are the Li concentrations at the solid–electrolyte interface phase, where the subscripts p and n denote PE and NE, respectively. The solid ohmic overpotentials can be calculated as [50]:
U p = U p ref ( x p ) = U p ref ( c s , p b c s , p max ) ,   U n = U n ref ( x n ) = U n ref ( c s , n b c s , n max ) ,
where σ p eff and σ n eff are the effective solid phase conductivities of the PE and NE, respectively. L PE and L NE are the thicknesses of the PE and NE, respectively. I is the electric current and A ele is the geometric area of the electrodes. The electrolyte diffusion overpotential due to the Li concentration gradient in the electrolyte can be calculated as [50]:
η e diff = ( 1 t + ) 2 R T F ln c e ( L ) c e ( 0 ) ,
where t + is the Li transference number.   c e ( L ) and c e ( 0 ) are the Li concentrations in the electrolyte phase at z = L and z = 0, respectively. R is the ideal gas constant, T is the temperature, and F is the Faraday constant. The electrolyte ohmic overpotential due to electrolyte conductivity can be calculated as [50]:
η e ohm = I 2 A ele ( L NE κ n eff + 2 L SEP κ SEP eff + L PE κ p eff ) ,
where κ n eff , κ SEP eff , and κ p eff are the electrolyte effective ionic conductivities in the NE, SEP, and PE, respectively.   L SEP is the thickness of SEP. The oxidation and reduction charge transfer coefficients, α ox and α red , are assumed to be equal to 0.5 for both electrodes. The positive and negative kinetic overpotentials can be expressed as [50]:
η p k = 2 R T F ln ( ξ p + ξ p 2 + 1 ) , η n k = 2 R T F ln ( ξ n + ξ n 2 + 1 ) ,
where
ξ p = I 2 i 0 , p S p = R s , p 6 ε s , p i 0 , p A ele L PE I ,   ξ n = I 2 i 0 , p S n = R s , n 6 ε s , n i 0 , n A ele L NE I ,
where S p and S n are the area of PE and NE, respectively. R s , p and R s , n are the particle radii in the PE and NE, respectively. ε s , p and ε s , n are the volume fractions of the active material in the PE and NE, respectively. The exchange current densities can be expressed as [50]:
i 0 , p = k 0 , p F ( c e ) α ox , p ( c s , p R s , p ) α ox , p ( c s , p max c s , p R s , p ) α ox , p i 0 , n = k 0 , n F ( c e ) α ox , n ( c s , n R s , n ) α ox , n ( c s , n max c s , n R s , n ) α ox , n ,
where c e is the Li-ion concentration in the electrolyte phase and α ox , p and α ox , n represent the oxidation charge transfer coefficients of the PE and NE, respectively.   k 0 , p and k 0 , n are the reaction rate constants of the PE and NE, respectively. By combining these equations, the output voltage (Equation (4)) of the battery can be expressed as:
U ( t ) = U p ref ( c s , p R s , p c s , p max ) U n ref ( c s , n R s , n c s , n max ) + η s , p ohm η s , n ohm + η e ohm + η e diff + η p k η n k
It is important to note that Equation (13) does not take into account the battery aging. The expression of battery voltage considering the battery aging is shown in Equation (20).

2.2. Aging Model of the Battery Cell

Battery aging results in capacity and power losses due to several physical mechanisms. The physical mechanism taken into account in this study is the formation of a solid electrolyte interface (SEI) layer on the NE. In addition, its effects on the battery behavior are considered. SEI layer formation is a phenomenon occurring on a negative carbonaceous electrode when its potential is low [49]. Under these conditions, Li starts a reaction with the solvent of the electrolyte which forms the SEI, as shown in Figure 1c. The reaction equation is [51]:
2 EC + 2 e + 2 Li + ( CH 2 OCO 2 Li ) 2 + C 2 H 4
During this reaction, lithium is withdrawn from the NE and trapped in the SEI. This leads to a loss of cycling Li and thus a capacity loss. In addition, the SEI formed is a resistive layer leading to an increase in the resistance. As described in Figure 1c, the electrolyte solvent is transported from the electrolyte bulk to the particle surface of the SEI, it then reacts and forms a new SEI. Equations involved in the aging calculation have been comprehensively discussed in previous work [48] and are listed here in Table 3.
It is worth mentioning that ε SEI * is the instantaneously formed SEI layer porosity. ε SEI is the SEI layer porosity. The ohmic overvoltage can be calculated as [52]:
η SEI = δ SEI κ SEI S n I ,   κ SEI = κ SEI , data ( 1 ε SEI ) γ SEI ,
where κ SEI , data is the index of the SEI layer electronic conductivity data and γ SEI is the SEI Bruggeman exponent. c solv is the concentration of the solvent inside the SEI layer. The bulk solvent concentration, c solv b , is considered to be constant outside the SEI layer, as shown in Figure 1c. Based on these equations, it is possible to evaluate the solvent concentration, c solv , and then the parasitic current of the SEI formation, j SEI , can be calculated. Furthermore, the variation in the SEI porosity, capacity loss, and the SEI layer thickness can be calculated.
For dynamic aging modeling, when the model is initialized, previous aging needs to be taken into account. Firstly, based on the initial capacity loss, the initial insertion rate of the NE can be modified as:
x n , int = x n , 0 + 3600 SOC int 100 ( C nom C loss , int SEI ) A ele L NE ε s , n F c s , n max ,
where x n , 0 is the negative insertion rate when the state of charge is zero (SOC = 0), SOC int is the initial SOC, C nom is the nominal capacity, and C loss , int SEI is the initial cycling Li loss due to SEI formation. Secondly, using the bulk solvent concentration and the battery average temperature, the initial solvent concentrations in the SEI layers can be calculated. With the assumption of a stationary state, it is assumed that the input current is 0 A and all time derivatives are null.
d ε SEI d t = 0 ,   c solv t = 0
There are several different consequences of aging on cell behavior. First, the current resulting from the parasitic reaction needs to be taken into account in the solid phase Li conservation. Therefore, the boundary condition of solid-phase Li+ conservation in Table 2 can be modified as:
D s c s , n r | r = R s = j f a s F + j SEI F
Then, the porosity of the NE is modified due to the clogging of pores by the SEI layer [52]:
ε e , n = ε e , n , int ε s , n 3 δ SEI ( t ) R s , n ,
where ε e , n , int is the initial porosity of the NE. In the event of a SEI layer that is too thick, the porosity is completely clogged. At this point, the battery reaches the end of its life and the model stops. The final effect is the increase in resistance due to the growth of the SEI layer. Another ohmic drop is added to the final voltage to take into account the effect of the SEI on the cell voltage; therefore, Equation (13) can be modified as:
U ( t ) = U p ref ( c s , p R s , p c s , p max ) U n ref ( c s , n R s , n c s , n max ) + η s , p ohm η s , n ohm + η e ohm + η e diff + η p k η n k η SEI
The total current intensity in the NE can be expressed as:
j tot = j f + a s , n j SEI

2.3. Thermal Models of the Battery Cell

The generated heat is described classically as the sum of irreversible and reversible heats [52]:
Q gen = ( ( η s , p ohm η s , n ohm + η e ohm + η e diff + η p k η n k η SEI ) I + T ( d U p d T ( x p ) d U n d T ( x n ) ) I ) ,
the first and second items on the right are the Joule effect and entropic effect in the battery, respectively. The equation of convective and radiative heat transfer is shown in Equation (3). Therefore, the energy balance of the battery cell can be expressed as:
d d t T = 1 m b C p b ( Q gen Q tra ) ,
where m b and C p b are the mass and specific thermal capacity of the battery cell, respectively. The nominal electrochemical model, thermal model, and aging model are coupled to study the battery performance. Several parameters are temperature dependent, and the Arrhenius law [52,53,54] can be used to evaluate these parameters, as shown in Table 4. The parameters of the battery model are summarized in Table 5.

2.4. Thermal Management Model of the 6s4p Battery Pack

The three-dimensional thermal management model of the 6s4p battery pack has been proposed in our previous work [32]. On this basis, two modes of the thermal management system are developed for the 6s4p battery pack using AMESim (Advanced Modeling Environment for performing Simulation of engineering systems) in this study, as shown in Figure 2. The cooling circuits are connected in series (Figure 2a) or in parallel (Figure 2b). Both modes of the thermal management system consist of four cold plates, a centrifugal pump, several thermal-hydraulic pipes, an expansion tank, a radiator, a fan, and a temperature sensor. Different colors and types of lines represent different kinds of connections, including electrical current, heat flow, control signals, and coolant. Every four battery cells in parallel are tightly packed and cooled by cold plates. The convective and radiative heat transfer equations are similar to that in Equation (3). The coolant was selected as a 50% volume concentration of ethylene glycol aqueous solution for the cooling system. The centrifugal pump is used to provide velocity to the coolant. Several thermal-hydraulic pipes connect the entire cooling circuit, and the expansion tank is used to stabilize hydraulic pressure throughout the circuit for safety. The radiator and fan are used to cool the coolant flows from the cold plates, and the temperature sensor is used to monitor the temperature of the coolant before it flows back to the centrifugal pump. Parameters of the cold plate, coolant, and ambient air are shown in Table 6. On this basis, 20 cases (cases A to T) were determined according to the C-rate and coolant mass flow rate, as shown in Table 7.

2.5. Triple-Step Nonlinear Cooling Control Strategy

It can be observed from Figure 2b that the thermal management system of the 6s4p battery pack has a certain symmetry. The input parameters into the coupled model are the current and ambient temperature, and the output parameters are the voltage and temperature of the cell. The temperature of battery 6 is the highest due to its structural design. Therefore, batteries 5 and 6 were selected as the key batteries for the study of the cooling strategy. The thermal models of batteries 5 and 6 are shown in Figure 3a. It can be observed that the thermal models are extremely complicated due to the consideration of various thermal terms. However, the thermal model of a single battery is suitable for linearization under certain conditions, as shown in Figure 3b. Moreover, the heat generation of each cell is the same due to the use of the same coupled model and initial conditions, as shown in Figure 3c. Only the heat generation of the cell and the heat transfer with the coolant are considered; a similar simplification can be found in the literature [44].
Based on the above model description and Figure 2b, battery 5 was selected for analysis of the cooling optimization strategy (see Figure 3b) and the thermal model can be expressed as:
Q b 5 = m b C p , b d d t T b 5 = Q gen ,   b 5 Q b 5 co ,
where Q b 5 , Q gen ,   b 5 , and Q b 5 co are the heat change, heat generation, and heat transferred from battery 5 to the coolant, respectively. According to the thermal model description and Equation (22), Q gen ,   b 5 can be simplified to:
Q gen ,   b 5 = Q gen ,   b 5 = U ir I + U r I T b 5 ,
where subscripts ir and r represent irreversible and reversible, respectively. Q b 5 co can be expressed as:
Q b 5 co = h b 5 co A b 5 co ( T b 5 T co ) ,
where h b 5 co and A b 5 co are the heat transfer coefficient and the contact area between battery 5 and the coolant, respectively. T b 5 and T co are the temperature of battery 5 and the coolant, respectively. The heat transfer coefficient, h b 5 co , can be calculated as [55]:
h b 5 co = N u co · k co D ,
where N u co and k co are the Nusselt number and the thermal conductivity of coolant, respectively. D is the equivalent diameter of the pipe, which is 35 mm. For the forced convection between the coolant and battery 5, turbulent flow is considered when the length of the pipe is much longer. Thus, N u co can be obtained by the experimental correlation of Sieder–Tate [56]:
N u co = 0.027 ( R e co ) 0.8 ( P r co ) 1 3 ( μ co μ w ) 0.14 ,
where R e co and P r co are the Reynolds number and Prandtl number of the coolant, respectively. μ co is the dynamic viscosity of the coolant and μ w is the dynamic viscosity near the pipe wall. R e co can be calculated as:
R e co = ρ co v co D μ co ,
where ρ co and v co are the density and flow velocity of the coolant, respectively. From Equations (27)–(29), the heat transfer coefficient, h b 5 co , can be expressed as:
h b 5 co = 0.027 ( R e co ) 0.8 ( P r co ) 1 3 ( μ co μ w ) 0.14 k co D = a 1 v co 0.8 ,
where a 1 = 0.027 ( ρ co ) 0.8 ( P r co ) 1 3 k co ( μ co ) 0.66 ( μ w ) 0.14 D 0.2 . Let u = v co 0.8 , from Equations (24)–(30), then the temperature of battery 5 ( T b 5 ) can be expressed as:
d d t T b 5 = U ir I + U r I T b 5 m b C p , b a 1 A b 5 co ( T b 5 T co ) m b C p , b u
The input is the current load, I , and the coolant flow rate, v , and the output is the temperature of battery 5.
As mentioned above, a new cooling and control strategy should be proposed to avoid too much cooling or an excessive increase in battery temperature. Fortunately, a thermal model of a single battery is available for the design of triple-step nonlinear cooling control. After the triple-step nonlinear cooling control strategy was solved for battery 5, this strategy was applied to the 6s4p battery pack to validate its applicability. Moreover, an extended PID cooling control strategy was also adopted to compare the results with those using a triple-step nonlinear cooling control strategy. Both the triple-step nonlinear method and extended PID method adjust the temperature of battery 5 by changing the coolant flow rate. The detailed design process is shown in Figure 4.
The design of steady state control, the reference variable feed-forward control, and the error feedback control are three main steps of the triple-step nonlinear method [44,57]. The design processes for battery 5 are illustrated in Figure 3d,e, respectively. T tv is the target temperature of battery 5, T b 5 is the actual temperature of battery 5, and u is a control law. Let A ( T b 5 ) = U ir I + U r I T b 5 m b C p , b , B ( T b 5 ) = a 1 A b 5 co ( T b 5 T co ) m b C p , b . Then, Equation (31) can be rewritten as:
d d t T b 5 = A ( T b 5 ) + B ( T b 5 ) u
Equation (32) can be adapted to different methods [58].
Step 1: Steady-state-like Control.
The control law, u s , changes dynamically with the output. The purpose of u s is to gradually adjust the temperature of battery 5 to a stable state. When a stable state is reached, d d t T b 5 = 0 , so combining Equations (31) and (32), we then obtain:
u s = A ( T b 5 ) B ( T b 5 ) = U ir I + U r I T b 5 a 1 A b 5 co ( T b 5 T co ) , B ( T b 5 ) 0
Step 2: Reference Variable Feed-forward Control.
The target temperature, T tv , is pre-set and has a prediction effect on its differential, d d t T tv . Therefore, introducing the tracking target temperature control law can effectively reduce the adjustment time. The control law, u f , is added based on Step 1. By letting d d t T b 5 = d d t T tv , we obtain:
u = u s + u f u f = 1 B ( T b 5 ) · d d t T tv = m b C p , b a 1 A b 5 co ( T b 5 T co ) · d d t T tv
Step 3: Error Feedback Control.
Based on the controller design method of the error system, the control law, u e , is obtained. This control law can further reduce the impact caused by the external disturbance and modeling deviation of the system and guarantee the steady-state of the system.
u = u s + u f + u e u e = 1 B ( T b 5 ) ( d d t T tv d d t T b 5 )
The error is the deviation between the target temperature and actual temperature of battery 5:
e = T tv T b 5
Based on Equations (35) and (36), the derivative of the input deviation can be expressed as:
d d t e = B ( T b 5 ) u e
The Lyapunov function is selected to solve the error feedback control law [44]:
V = 1 2 e 2 + 1 2 K 0 ( e d t ) 2 d d t V = e · d d t e + K 0 e e d t
Let d d t V = K 1 e 2 , and combine Equations (37) and (38), then the control law can be expressed as:
B ( T b 5 ) · u e = K 1 e K 0 e d t
where K 1 > 0 and K 0 > 0 . When d d t V 0 , a stable state is achieved. By simplifying Equation (39), the control law can be expressed as:
u e = K 1 e B ( T b 5 ) + K 0 B ( T b 5 ) e d t u e = m b C p , b K 1 a 1 A b 5 co ( T b 5 T co ) e m b C p , b K 0 a 1 A b 5 co ( T b 5 T co ) e d t
Equations (33), (34), and (40) are summarized as the final control law:
u = A ( T b 5 ) B ( T b 5 ) + 1 B ( T b 5 ) · d d t T tv + K 1 e B ( T b 5 ) + K 0 B ( T b 5 ) e d t
The final control law is:
u = U ir I + U r I T b 5 a 1 A b 5 co ( T b 5 T co ) m b C p , b a 1 A b 5 co ( T b 5 T co ) · d d t T tv m b C p , b K 1 a 1 A b 5 co ( T b 5 T co ) e m b C p , b K 0 a 1 A b 5 co ( T b 5 T co ) e d t
We now consider the system modeling error, d i. Under the framework of state stability theory, the robust stability of the system needs to be analyzed. Let ζ = e d t , then combining Equations (37) and (39), we obtain:
e = ζ ˙ d e d t = K 1 e K 0 ζ + di ,
and Equation (43) can be expressed in the matrix form as:
d e ¯ d t = H e ¯ + G di ,
where e ¯ = [ ξ   e ] T , H = [ 0 K 0   1 K 1 ] , and G = [ 0   1 ] T . Equation (44) is a linear system. As mentioned above, in the description of the triple-step nonlinear method, the asymptotic stability of the system is demonstrated without considering the disturbance, di . According to the conclusion of input state stability of a linear system, the closed-loop system has robust stability [59].

3. Model Validation

The whole system is established in the one-dimensional software package AMESim, including the electrochemical model, thermal model, aging model, and cooling control model. Each module provided by AMESim is modeled based on physical mechanisms. The connection relationship between the modules in the system is stated in the above description, and the practicability of each module has been proven through extensive research [44,45,58]. Therefore, AMESim is an experimental platform that can be used to replace the real experimental platform to explore the battery performance in the thermal management system. Moreover, the model of a 5 A · h LFP/C battery at 25 °C was adopted for the validation of the battery model used in this study, as shown in Figure 5. The equipment used in the experiment has been introduced in previous work [51]. In this work, a cycle consists of a constant-current (CC) charging process and a constant-current discharging process. The experimental charge–discharge data were then compared with the simulated results in this study, as shown in Figure 6.
Figure 6a shows the voltage curves of the charge–discharge process at different C-rates. It can be seen that the simulated curves agree well with the experimental results. Then, the experimental temperature data at different C-rates were then compared with the simulated results, as shown in Figure 6b. The cell temperature can reach 47.5 °C at five C-rates after one cycle of the charge–discharge process. The temperature curves obtained by simulation at different C-rates were also in good agreement with the experimental values. Finally, the decrease in capacity was compared and validated. The aging mechanism of the SEI formation was considered in the simulations for all conditions. The initial capacity loss, C loss , int SEI , was 0.01 A · h according to the battery factory parameters, which was set as the initial condition in the modeling, as shown in Table 6. The loss of capacity (LOC) (the ratio of the lost capacity to nominal capacity, 5 A · h ) was calculated by Equation (44) at different C-rates, and then compared with the experimental data, as shown in Figure 6c. The comparison presents good agreement, which indicates the model accuracy. The simulated results were slightly different from the experimental data, which was possibly due to some model parameters not being totally applicable to the experimental tests. In addition, the initial conditions and boundary conditions of the battery in the experiment and simulation are not exactly the same. The simulated battery temperature was also compared with that in the literature [32,60] during the charge–discharge process, showing a good agreement. For example, the battery temperature in Figure 6b is verified in the literature [32].

4. Results and Discussion

4.1. Performance of Battery Cell and 6s4p Battery Pack under Natural Cooling

The thermal behavior and aging performance of a battery cell and the 6s4p battery pack were first investigated under natural cooling. The charge–discharge process of a single cell was first charged and then discharged and cycled at different C-rates and ambient temperatures ( T amb ). The temperature and LOC evolutions of a single cell are shown in Figure 7. Taking the time required for one charge–discharge cycle at one C-rate as the standard, Figure 7a,b shows the battery temperature and LOC at C-rates from one to five. As the C-rate increases, the time for a single charge–discharge cycle becomes shorter, so the number of cycles increases. It can be observed that both the battery temperature and LOC increase with the increase in the C-rate. The maximum battery temperature can even exceed 50 °C at a C-rate of five, which is essentially consistent with the results of our previous research [32]. The LOC continues to increase as the charge–discharge cycles proceed, and also correspondingly increases with the increasing C-rate. The initial LOC is 0.01 A · h , accounting for 0.2%. The final LOC exceeds 0.21% at a C-rate of five, which means that the newly added LOC reaches 0.01% after only about five charge–discharge cycles. For a C-rate of five, the evolutions of battery temperature and LOC at different T amb are presented, as shown in Figure 7c,d. Both the battery temperature and LOC increase when the T amb increases from 25 °C to 65 °C. The maximum battery temperature and LOC can even reach 80 °C and 1% at T amb = 65   , respectively. A thermal runaway can occur in a battery cell when the temperature exceeds 80 °C, which is beyond the scope of this study. Therefore, both the C-rate and T amb are critical for the thermal behavior and aging performance of a battery cell.
Before studying the thermal behavior and aging performance of the 6s4p battery pack, the temperature and aging performance of a single cell are further summarized and analyzed. The temperature, LOC, SEI layer thickness, and SEI layer porosity of a single battery at different C-rates and ambient temperatures are shown in Figure 8. As shown in the figures, the C-rate plays an important role in the thermal behavior during the charge–discharge process of the battery cell. Moreover, the temperature rise is relatively uniform with the T amb increasing from 25 °C to 65 °C, regardless of the C-rate. This indicates that although the battery temperature is greatly affected by the T amb , the difference in the effect under different C-rates is small, as shown in Figure 8a. However, this situation is not the same for the aging performance of the battery cell, as shown in Figure 8b–d. As the C-rate and T amb increase, the uptrend of LOC, SEI layer thickness, and SEI layer porosity is slow at first and then sharp. It reveals that the high C-rate and T amb extensively accelerate the battery aging, while normal C-rate and T amb have less effect. Since the initial LOC is considered to be 0.2% (0.01 A · h ), the initial SEI layer thickness and porosity are, 4.03848 nm and 0.01, respectively.
Based on the above analysis, the thermal behavior and aging performance of the 6s4p battery under natural cooling are then studied. The heat conduction between battery cells and heat convection and radiation between the cell and ambient air are considered. The evolutions of maximum battery temperature, LOC, SEI layer thickness, and SEI layer porosity in the 6s4p battery pack at T amb = 25   are shown in Figure 9. The maximum battery temperature in the 6s4p battery pack can reach 80 °C at 5 C-rate, even if the T amb is only 25 °C, as shown in Figure 9a. This means that the rapid heat generation of the battery pack at a high C-rate leads to a sharp increase in temperature, which requires the intervention of a thermal management system for timely cooling and heat dissipation. In addition, the SEI formation continuously occurs, resulting in the continuous increase of SEI layer thickness and porosity, and thus the LOC of the battery. Specifically, the maximum and minimum temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack at different C-rates are shown in Figure 10. As shown in the figures, squares of different styles and colors represent different meanings. Where green represents the initial value, red represents the maximum value, and blue represents the minimum value. The increase of the C-rate has a linear effect on the battery temperature, while the effect on the aging performance has a sudden change. The aging at low C-rates is small or even negligible (such as 1 C and 2 C), but the aging at high C-rates is not negligible (such as 4 C and 5 C).

4.2. Performance of 6s4p Battery Pack under Two Modes Cooling System

After exploring the performance of the 6s4p battery pack under natural cooling, the thermal behavior and aging performance of the 6s4p battery pack in two cooling modes are further investigated and compared. For both parallel and series cooling, the pump rotary speeds of 1000 and 3000 r · min 1 are considered to provide the kinetic energy for coolant. According to the pump performance in the AMESim, the corresponding mechanical power provided by the pump to the coolant is 7.168 W and 218.873 W, and 2.620 W and 81.552 W under parallel and series cooling, respectively, as shown in Table 7. This indicates that parallel cooling is more energy efficient while series cooling requires a higher pump rotary speed to obtain the power. Further, the volumetric flow rate and velocity of the coolant can be obtained. It is worth mentioning that the flow in the four cold plates is the same under parallel cooling, which means a quarter of the total flow rate in each cold plate, as shown in Table 7. Therefore, for parallel cooling, the mass flow rates of each cold plate are, respectively, 0.081 and 0.276 kg · s 1 when the pump rotary speeds are 1000 and 3000 r · min 1 . For series cooling, the mass flow rates of each cold plate are always, respectively, 0.118 and 0.408 kg · s 1 . Combining the C-rates from 1 to 5, 20 cases (cases A to T) can be divided, as shown in Table 7.
Figure 11a,b show the evolutions of the maximum battery temperature in the 6s4p battery pack for the 20 cases. For both parallel and series cooling, the ambient temperature is set to 25 °C. Overall, under the same C-rate and coolant flow rate, the maximum battery temperature of the series and parallel cooling systems is close. From the specific values obtained from the simulation, the maximum battery temperature in the series connection is slightly higher. In addition, it is noteworthy that the maximum battery temperature in the 6s4p battery pack does not always decrease with the coolant flow rate increases. This phenomenon is particularly evident in parallel cooling, as shown in Figure 11a. When the C-rate is high, such as 4 C and 5 C, the heat generation of the battery is higher. A Larger coolant flow rate allows for faster heat exchange and more heat removal. Therefore, the maximum battery temperature at a higher coolant flow rate (Case F and G) is lower than that at a lower coolant flow rate (Case A and B). The opposite situation occurs when the C-rate is low, such as 1 C and 2 C. The battery does not generate much heat at this time, and a smaller coolant flow rate can carry out heat exchange for a longer time. Therefore, the maximum battery temperature at a lower coolant flow rate (Case D and E) is lower than that at a higher coolant flow rate (Case I and J). Further, the evolutions of the maximum temperature difference in the 6s4p battery pack for the 20 cases, as shown in Figure 11c,d. It can be observed that the maximum temperature difference in the 6s4p battery pack decreases with the coolant flow rate increasing, regardless of the C-rate. As the coolant flow rate increases, the maximum temperature difference decreases more with parallel cooling than that with series cooling. This indicates that increasing the coolant flow rate can better reduce the temperature difference between the cells in the battery pack when parallel cooling is used. These results also demonstrate that at a high C-rate, a further increase in the coolant flow rate is required to control the temperature and temperature difference of the battery pack. Furthermore, parallel cooling is relatively more advantageous.
Based on the above analysis, the temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack under parallel cooling are further analyzed. After the charge–discharge process, the maximum and minimum battery temperatures at high and low coolant flow rates are calculated at different C-rates, as shown in Figure 12a. For high C-rates (4 C and 5 C), the maximum and minimum temperature at a higher coolant flow rate (0.276 kg · s 1 ) is lower than those at a lower coolant flow rate (0.081 kg · s 1 ). For low C-rates (1 C and 2 C), the maximum and minimum temperatures at a higher coolant flow rate are higher than those at a lower coolant flow rate. These results also validate the results in Figure 11. Further, the characteristics of this change caused by the difference in C-rate and coolant flow rate are also reflected in the parameters of battery aging, as shown in Figure 12b–d. The use of a higher coolant flow rate at high C-rates can not only reduce the maximum temperature of the battery pack but also reduce the formation of SEI, thereby reducing capacity loss. Moreover, applying a higher coolant flow rate at low C-rates, not only is the optimization of battery temperature and aging performance not as good as those at a lower coolant flow rate, but a higher coolant flow rate also consumes more pump power. Therefore, high coolant flow is required at high C-rates, and low coolant flow is more suitable at low C-rates. All these show that it is necessary to select the appropriate coolant flow to ensure the performance of the battery. This requires timely adjustment of the coolant flow rate according to the working conditions of the battery.

4.3. Cooling Optimization of 6s4p Battery Pack under Different Cycle Conditions

As mentioned above, the triple-step nonlinear cooling control strategy is first solved for a single cell, and then this strategy is applied to the 6s4p battery pack to validate its applicability. Moreover, an extended PID cooling control strategy is also adopted to compare the results with those using a triple-step nonlinear cooling control strategy. Both the triple-step nonlinear method and the extended PID method adjust the temperature by changing the coolant flow rate. The effect of two methods is verified under constant current (2 C-rate), New European Driving Cycle (NEDC), and Urban Severe Driving Cycle (US06). Each working condition circulates two times to test the effect of control and optimization. NEDC consists of 4 urban cycles and 1 suburban cycle, and US06 represents an aggressive driving cycle. The initial temperature of the 6s4p battery pack is 40 °C, and the temperature of coolant and ambient is 25 °C.
Under 2 C-rate charge–discharge cycles, the heat generation of each battery cell is shown in Figure 13a. It can be seen that the heat generation in the discharge process is higher than that in the charging process. At the end of the charge and discharge process, the heat generation increases sharply. The well-tuned four parameters of the extended PID method are: K i p = 2 × 10 4 (inner proportional element), K p = 6 × 10 2 (proportional element), K i = 2 × 10 6 (integral element), K d = 1 × 10 4 (derivative element). The comparison of coolant velocity under extended PID and triple-step nonlinear control is shown in Figure 13b. The coolant velocity increases continuously as the charge–discharge proceeds under the extended PID control, while under the triple-step nonlinear control, the coolant velocity increases smoothly. This indicates less pump power under the triple-step nonlinear control. The temperature comparison of battery 6, battery 5, and cold plate under extended PID and triple-step nonlinear control is shown in Figure 13c. As can be seen from the figure, under the triple-step nonlinear control, the temperature of the cold plate can be controlled at about 25 °C in a short time, while the extended PID control takes a little more time, which verifies the validity of the triple-step nonlinear method. Further, through the temperature evolution of batteries 5 and 6, the application effect of the two control methods in the 6s4p battery pack can be analyzed. The temperature of batteries 6 and 5 can be controlled to be below 35 °C and 32 °C in about 10 min, respectively. Although the subsequent charge–discharge process causes some temperature fluctuations, the temperature of batteries 6 and 5 can still be controlled below 34 °C and 31 °C, respectively. In addition, whether it is the temperature of the battery or the cold plate, the difference between the two methods is small, and the triple-step nonlinear method has a slight advantage. This slight advantage can be inferred from the aging performance of the batteries, as shown in Figure 13d. The LOC of batteries 6 and 5 under triple-step nonlinear control is smaller than that under the extended PID control.
Under NEDC cycles, the heat generation of each battery cell is shown in Figure 14a. The heat generation also increases dramatically at the end of each urban and suburban cycle. The well-tuned four parameters of the extended PID method are: K i p = 3 × 10 3 , K p = 1 × 10 2 , K i = 2 × 10 8 , K d = 1 × 10 4 . The comparison of coolant velocity under extended PID and triple-step nonlinear control is shown in Figure 14b. The coolant velocity increases continuously under the extended PID control, while under the triple-step nonlinear control, the coolant velocity correspondingly changes with the change of heat generation. This indicates a more efficient control of coolant flow rate and less pump power consumption under the triple-step nonlinear control. The temperature comparison of battery 6, battery 5, and the cold plate is also shown in Figure 14c. The triple-step nonlinear method still shows a good control effect on the temperature of the cold plate. The temperature of batteries 6 and 5 can be controlled to be below 34 °C and 32 °C in about 6 min, respectively. During the subsequent cycle, the temperature of batteries 6 and 5 can further be controlled to be below 29 °C and 28 °C, respectively. The triple-step nonlinear method still has a slight advantage, which can be inferred from the aging performance of the batteries, as shown in Figure 14d.
Under US06 cycles, the heat generation of each battery cell is shown in Figure 15a. The aggressive driving cycles lead to dramatic and fluctuating heat generation. The well-tuned four parameters of the extended PID method are: K i p = 3 × 10 3 , K p = 1.5 × 10 2 , K i = 2 × 10 7 , K d = 1 × 10 4 . The comparison of coolant velocity under extended PID and triple-step nonlinear control is shown in Figure 15b. The coolant velocity still increases continuously under the extended PID control, while under the triple-step nonlinear control, the coolant velocity correspondingly fluctuates with the change of the heat generation. This also indicates a more efficient control of coolant flow rate and less pump power consumption under the triple-step nonlinear control. The temperature comparison of battery 6, battery 5, and the cold plate is also shown in Figure 15c. The triple-step nonlinear method still shows a good control effect on the temperature of the cold plate. The temperature of batteries 6 and 5 can be controlled to be below 38 °C and 34 °C in about 9 min, respectively. During the subsequent cycle, the temperature of batteries 6 and 5 can further be controlled to be below 36 °C and 32 °C, respectively. The triple-step nonlinear method still has a slight advantage, which can be inferred from the aging performance of the batteries, as shown in Figure 15d.

5. Conclusions

Two modes of thermal management system for 6s4p battery pack are established based on one-dimensional modeling. The single particle model with the consideration of aging is adopted for the battery. Both the battery temperature and loss of capacity (LOC) increase when the C-rate and ambient temperature ( T amb ) increase. The maximum battery temperature and LOC can even reach 80 °C and 1% at T amb = 65   , respectively. High C-rate and T amb can extensively accelerate the battery aging, while normal C-rate and T amb have weaker effect. The maximum battery temperature in the 6s4p battery pack can reach 80 °C at 5 C-rate, even if the T amb is only 25 °C. The comparison of pump power, thermal behavior, and aging performance indicate parallel cooling is relatively more advantageous. When the C-rate is high (4 C and 5 C), the maximum battery temperature at a higher coolant flow rate (0.276 kg · s 1 ) is lower than that at a lower coolant flow rate (0.081 kg · s 1 ). When the C-rate is low (1 C and 2 C), the maximum battery temperature at a lower coolant flow rate is lower than that at a higher coolant flow rate. The characteristics of this change caused by the difference in C-rate and coolant flow rate are also reflected in the battery aging. Further, a cooling control strategy for the key battery is firstly proposed based on the triple-step nonlinear method, and then the optimization effect of this strategy in the battery pack is validated under 2 C charge–discharge cycles, NEDC cycles, and US06 cycles. Moreover, the extended PID control strategy is constructed and compared with the triple-step nonlinear method. The results show that the temperature and LOC of key batteries and the temperature of the cold plate are lower than that of the extended PID method. This verifies the validity of the triple-step nonlinear method and shows an advantage compared with the extended PID method.

Author Contributions

Conceptualization, H.Z.; methodology, H.Z.; software, H.Z. and H.C.; validation, H.Z. and H.C.; formal analysis, H.Z., H.C. and H.F.; investigation, H.Z.; resources, H.F.; data curation, H.C.; writing—original draft preparation, H.Z.; writing—review and editing, H.Z., H.C. and H.F.; funding acquisition, H.F.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Q tra Heat transferred from the battery to the ambient (W) u e Error feedback control law
hconvConvective heat transfer coefficient ( W · m 2 · K 1 ) e Input error
A b Contact area between the battery and the ambient air ( m 2 )Greek symbols
T b Battery temperature (K) χ Surface emissivity
T amb Ambient temperature (K) σ B Stefan-Boltzmann constant
S p Active area of positive electrode ( m 2 ) ρ b Density of battery ( kg · m 3 )
S sep Active area of separator ( m 2 ) η ( L ) Overpotential at z = L
S n Active area of negative electrode ( m 2 ) η ( 0 ) Overpotential at z = 0
A ele Active area of electrode ( m 2 ) ϕ e ( L ) Electrolyte potential at z = L
I Electric current (A) ϕ e ( 0 ) Electrolyte potential at z = 0
k b Thermal conductivity of battery ( W · m 1 · K 1 ) ψ Hysteresis factor
C p , b Specific thermal capacity of battery ( J · kg 1 · K 1 ) ε e Volume fraction of electrolyte
U p Thermodynamic potential of the positive electrode (V) σ eff Effective solid phase conductivity ( S · m 1 )
U n Thermodynamic potential of the negative electrode (V) ϕ s Potential of solid phase (V)
zBattery thickness direction κ eff Electrolyte effective ionic conductivity ( S · m 1 )
c s , p b Bulk Li concentration in the solid phase for positive electrode ( mol · m 3 ) κ D eff Electrolyte effective ionic conductivity ( S · m 1 )
c s , n b Bulk Li concentration in the solid phase for negative electrode ( mol · m 3 ) α ox Charge transfer coefficient of the anodic reaction
c s , p max Maximum Li concentration in the positive electrode ( mol · m 3 ) α red Charge transfer coefficient of the reduction reaction
c s , n max Maximum Li concentration in the negative electrode ( mol · m 3 ) η Overpotential (V)
U p ref Reference potential of the positive electrode (V) γ e Electrolyte Bruggeman exponent
U n ref Reference potential of the negative electrode (V) σ Solid phase conductivity
f ts Transition factor ε s Volume fraction of the active material
xInsertion rate η p diff Solid phase overpotential of positive electrode (V)
U p ch , ref Charge reference potentials of the positive electrode (V) η n diff Solid phase overpotential of negative electrode (V)
U p dch , ref Discharge reference potentials of the positive electrode (V) η s , p ohm Solid ohmic overpotential of positive electrode (V)
U n ch , ref Charge reference potentials of the negative electrode (V) η s , n ohm Solid ohmic overpotential of negative electrode (V)
U n dch , ref Discharge reference potentials of the negative electrode (V) η e diff Electrolyte diffusion overpotential (V)
T ref Reference temperature η e ohm Electrolyte ohmic overpotential (V)
c s Li concentration in the active material particles ( mol · m 3 ) β SEI Reaction charge transfer coefficient of the SEI layer formation
r Radial coordinate inside a spherical particle (m) δ SEI SEI layer thickness (m)
D s Solid phase diffusion coefficient ( m 2 · s 1 ) κ SEI SEI layer conductivity ( S · m 1 )
j f Current per volume unit ( A · m 3 ) ε SEI SEI layer porosity
a s Specific interfacial surface area ( m 1 ) ε SEI * Instantaneous formed SEI layer porosity
F Faraday constant ( C · mol 1 ) ρ SEI SEI layer density ( kg · m 3 )
c e Concentration of the electrolyte ( mol · m 3 ) η SEI SEI ohmic overvoltage (V)
D e eff Effective electrolyte phase diffusion coefficient ( m 2 · s 1 ) κ SEI , data Index of the SEI layer electronic conductivity data ( S · m 1 )
t + Lithium transference number γ SEI SEI Bruggeman exponent
L NE Thickness of the positive electrode ε e , n Porosity of the negative electrode
L PE Thickness of the negative electrode ρ co Density of the coolant ( kg · m 3 )
i 0 Exchange current density (A) μ co Dynamic viscosity of the coolant ( Pa · s )
R Ideal gas constant μ w Dynamic viscosity near the wall in the pipe ( Pa · s )
T Temperature
D e Electrolyte phase diffusion coefficient ( m 2 · s 1 )Subscripts and superscripts
k 0 Reaction rate constant ( m 2.5 · mol 0.5 · s 1 )convConvection
f ± Electrolyte activity coefficientbBattery
R s Radius of spherical particle (m)sSolid phase
c s , p R s , p Li concentration at the positive electrode-electrolyte interface phase ( mol · m 3 )pPositive electrode
c s , n R s , n Li concentration at the negative electrode-electrolyte interface phase ( mol · m 3 )nNegative electrode
F Faraday constantmaxMaximum
j SEI Parasitic current density for SEI layer formation ( A · m 2 )eElectrolyte
k SEI SEI formation reaction constant ( m 7 · mol 2 · K 1 )effEffective
c solv partedg Concentration of the solvent at the particle edge ( mol · m 3 )diffDifference
U SEI Solvent reduction potential of the SEI layer (V)intInitial
c solv Concentration of the solvent ( mol · m 3 )irIrreversible
D solv SEI layer solvent diffusivity ( m 2 · s 1 )rReversible
c solv b Bulk concentration of the solvent ( mol · m 3 )coCoolant
C loss SEI Cycling lithium loss due to SEI formation ( A · h )tvTarget value
LOC SEI Loss of capacity due to SEI formationb5Battery 5
C tot Total capacity ( A · h )
M SEI SEI molar mass ( kg · mol 3 )Acronyms
x n , 0 Negative insertion rate when the state of charge is zeroEVsElectric vehicles
C nom Nominal capacity ( A · h )LIBsLithium-ion batteries
m b Mass of the battery (kg)BTMBattery thermal management
j tot Total current intensity ( A · m 3 )NENegative electrode
Q gen Generated heat of battery (W)PEPositive electrode
Q b 5 co Heat transferred from cell 5 to coolant (W)SEIsolid electrolyte interface
h b 5 co Heat transfer coefficient between cell 5 and coolant ( W · m 2 · K 1 )SPSingle particle
A b 5 co Contact area between cell 5 and coolant ( m 2 )LFPLithium iron phosphate
N u co Nusselt number of the coolantSOCstate of charge
R e co Reynolds number of the coolantAMESimAdvanced Modeling Environment for performing Simulation of engineering systems
P r co Prandtl number of the coolantNCCNegative current collector
v co Flow velocity of the coolant ( m · s 1 )SEPSeparator
u s Steady-state-like control lawPCCPositive current collector
u f Reference variable feed-forward control lawLOCLoss of capacity

References

  1. Al-Zareer, M.; Dincer, I.; Rosen, M.A. A Review of Novel Thermal Management Systems for Batteries. Int. J. Energy Res. 2018, 42, 3182–3205. [Google Scholar] [CrossRef]
  2. Patil, M.S.; Seo, J.H.; Panchal, S.; Jee, S.W.; Lee, M.Y. Investigation on Thermal Performance of Water-Cooled Li-Ion Pouch Cell and Pack at High Discharge Rate with U-Turn Type Microchannel Cold Plate. Int. J. Heat Mass Transf. 2020, 155, 119728. [Google Scholar] [CrossRef]
  3. Panchal, S.; Gudlanarva, K.; Tran, M.K.; Fraser, R.; Fowler, M. High Reynold’s Number Turbulent Model for Micro-Channel Cold Plate Using Reverse Engineering Approach for Water-Cooled Battery in Electric Vehicles. Energies 2020, 13, 1638. [Google Scholar] [CrossRef] [Green Version]
  4. Ye, M.; Guo, H.; Xiong, R.; Yu, Q. A Double-Scale and Adaptive Particle Filter-Based Online Parameter and State of Charge Estimation Method for Lithium-Ion Batteries. Energy 2018, 144, 789–799. [Google Scholar] [CrossRef]
  5. Tran, M.K.; Panchal, S.; Chauhan, V.; Brahmbhatt, N.; Mevawalla, A.; Fraser, R.; Fowler, M. Python-Based Scikit-Learn Machine Learning Models for Thermal and Electrical Performance Prediction of High-Capacity Lithium-Ion Battery. Int. J. Energy Res. 2022, 46, 786–794. [Google Scholar] [CrossRef]
  6. Hekmat, S.; Molaeimanesh, G.R. Hybrid Thermal Management of a Li-Ion Battery Module with Phase Change Material and Cooling Water Pipes: An Experimental Investigation. Appl. Therm. Eng. 2020, 166, 114759. [Google Scholar] [CrossRef]
  7. Li, L.; Chen, X.; Yuan, Q.; Wang, T.; Ji, H.; Papović, S.; Raleva, K.; Pan, F.; Yang, T.; Li, J. Effects of Minor Mechanical Deformation on the Lifetime and Performance of Commercial 21700 Lithium-Ion Battery. J. Electrochem. Soc. 2022, 169, 060544. [Google Scholar] [CrossRef]
  8. Latini, D.; Vaccari, M.; Lagnoni, M.; Orefice, M.; Mathieux, F.; Huisman, J.; Tognotti, L.; Bertei, A. A Comprehensive Review and Classification of Unit Operations with Assessment of Outputs Quality in Lithium-Ion Battery Recycling. J. Power Sources 2022, 546, 231979. [Google Scholar] [CrossRef]
  9. Ghaeminezhad, N.; Wang, Z.; Ouyang, Q. A Review on Lithium-Ion Battery Thermal Management System Techniques: A Control-Oriented Analysis. Appl. Therm. Eng. 2023, 219, 119497. [Google Scholar] [CrossRef]
  10. Wu, W.; Wang, S.; Wu, W.; Chen, K.; Hong, S.; Lai, Y. A Critical Review of Battery Thermal Performance and Liquid Based Battery Thermal Management. Energy Convers. Manag. 2019, 182, 262–281. [Google Scholar] [CrossRef]
  11. Kim, J.; Oh, J.; Lee, H. Review on Battery Thermal Management System for Electric Vehicles. Appl. Therm. Eng. 2019, 149, 192–212. [Google Scholar] [CrossRef]
  12. Aswin Karthik, C.; Kalita, P.; Cui, X.; Peng, X. Thermal Management for Prevention of Failures of Lithium Ion Battery Packs in Electric Vehicles: A Review and Critical Future Aspects. Energy Storage 2020, 2, e137. [Google Scholar] [CrossRef] [Green Version]
  13. Lin, J.; Liu, X.; Li, S.; Zhang, C.; Yang, S. A Review on Recent Progress, Challenges and Perspective of Battery Thermal Management System. Int. J. Heat Mass Transf. 2021, 167, 120834. [Google Scholar] [CrossRef]
  14. Mohammadian, S.K.; Zhang, Y. Thermal Management Optimization of an Air-Cooled Li-Ion Battery Module Using Pin-Fin Heat Sinks for Hybrid Electric Vehicles. J. Power Sources 2015, 273, 431–439. [Google Scholar] [CrossRef]
  15. Mohammadian, S.K.; Rassoulinejad-Mousavi, S.M.; Zhang, Y. Thermal Management Improvement of an Air-Cooled High-Power Lithium-Ion Battery by Embedding Metal Foam. J. Power Sources 2015, 296, 305–313. [Google Scholar] [CrossRef]
  16. An, Z.; Shah, K.; Jia, L.; Ma, Y. A Parametric Study for Optimization of Minichannel Based Battery Thermal Management System. Appl. Therm. Eng. 2019, 154, 593–601. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, J.; Li, H.; Li, W.; Shi, J.; Wang, H.; Chen, J. Thermal Characteristics of Power Battery Pack with Liquid-Based Thermal Management. Appl. Therm. Eng. 2020, 164, 114421. [Google Scholar] [CrossRef]
  18. Mortazavi, B.; Yang, H.; Mohebbi, F.; Cuniberti, G.; Rabczuk, T. Graphene or H-BN Paraffin Composite Structures for the Thermal Management of Li-Ion Batteries: A Multiscale Investigation. Appl. Energy 2017, 202, 323–334. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, W.; Ling, G.; Zhuang, L.; Liang, Z. The Effect of Reducing the Thermal Contact Resistance on the Performance of Battery Thermal Management System. Int. J. Energy Res. 2021, 45, 9970–9982. [Google Scholar] [CrossRef]
  20. Liu, F.F.; Lan, F.C.; Chen, J.Q.; Li, Y.G. Experimental Investigation on Cooling/Heating Characteristics of Ultra-Thin Micro Heat Pipe for Electric Vehicle Battery Thermal Management. Chinese J. Mech. Eng. 2018, 31, 53. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Wei, K. Experimental and Numerical Study of a Passive Thermal Management System Using Flat Heat Pipes for Lithium-Ion Batteries. Appl. Therm. Eng. 2020, 166, 114660. [Google Scholar] [CrossRef]
  22. Li, X.; He, F.; Zhang, G.; Huang, Q.; Zhou, D. Experiment and Simulation for Pouch Battery with Silica Cooling Plates and Copper Mesh Based Air Cooling Thermal Management System. Appl. Therm. Eng. 2019, 146, 866–880. [Google Scholar] [CrossRef]
  23. Dan, D.; Yao, C.; Zhang, Y.; Zhang, H.; Zeng, Z.; Xu, X. Dynamic Thermal Behavior of Micro Heat Pipe Array-Air Cooling Battery Thermal Management System Based on Thermal Network Model. Appl. Therm. Eng. 2019, 162, 114183. [Google Scholar] [CrossRef]
  24. Zhang, L.; Yuan, Q.; Hu, S.; Xu, X. Research on Performance of Thermal Management System Integrated with Multiple Heat Exchange Methods. Ionics 2022, 28, 789–799. [Google Scholar] [CrossRef]
  25. Kabir, M.M.; Demirocak, D.E. Degradation Mechanisms in Li-Ion Batteries: A State-of-the-Art Review. Int. J. Energy Res. 2017, 41, 1963–1986. [Google Scholar] [CrossRef]
  26. Xiong, R.; Pan, Y.; Shen, W.; Li, H.; Sun, F. Lithium-Ion Battery Aging Mechanisms and Diagnosis Method for Automotive Applications: Recent Advances and Perspectives. Renew. Sustain. Energy Rev. 2020, 131, 110048. [Google Scholar] [CrossRef]
  27. Birkl, C.R.; Roberts, M.R.; McTurk, E.; Bruce, P.G.; Howey, D.A. Degradation Diagnostics for Lithium Ion Cells. J. Power Sources 2017, 341, 373–386. [Google Scholar] [CrossRef]
  28. Kamyab, N.; Weidner, J.W.; White, R.E. Mixed Mode Growth Model for the Solid Electrolyte Interface (SEI). J. Electrochem. Soc. 2019, 166, A334–A341. [Google Scholar] [CrossRef]
  29. Yang, S.-c.; Hua, Y.; Qiao, D.; Lian, Y.-b.; Pan, Y.-w.; He, Y.-l. A Coupled Electrochemical-Thermal-Mechanical Degradation Modelling Approach for Lifetime Assessment of Lithium-Ion Batteries. Electrochim. Acta 2019, 326, 134928. [Google Scholar] [CrossRef]
  30. Sun, S.; Guan, T.; Cheng, X.; Zuo, P.; Gao, Y.; Du, C.; Yin, G. Accelerated Aging and Degradation Mechanism of LiFePO4/Graphite Batteries Cycled at High Discharge Rates. RSC Adv. 2018, 8, 25695–25703. [Google Scholar] [CrossRef] [PubMed]
  31. Liang, J.; Gan, Y.; Yao, M.; Li, Y. Numerical Analysis of Capacity Fading for a LiFePO4 Battery under Different Current Rates and Ambient Temperatures. Int. J. Heat Mass Transf. 2021, 165, 120615. [Google Scholar] [CrossRef]
  32. Zhang, H.; Li, C.; Zhang, R.; Lin, Y.; Fang, H. Thermal Analysis of a 6s4p Lithium-Ion Battery Pack Cooled by Cold Plates Based on a Multi-Domain Modeling Framework. Appl. Therm. Eng. 2020, 173, 115216. [Google Scholar] [CrossRef]
  33. An, Z.; Jia, L.; Wei, L.; Yang, C. Numerical Modeling and Analysis of Thermal Behavior and Li+ Transport Characteristic in Lithium-Ion Battery. Int. J. Heat Mass Transf. 2018, 127, 1351–1366. [Google Scholar] [CrossRef]
  34. Zhao, P.; Li, M.; Kang, J.; Rizzoni, G. Analysis of Fading Characteristics of a Lithium Ion Battery Based on an Integration Model. Int. J. Heat Mass Transf. 2017, 104, 1317–1324. [Google Scholar] [CrossRef]
  35. Dong, T.; Peng, P.; Jiang, F. Numerical Modeling and Analysis of the Thermal Behavior of NCM Lithium-Ion Batteries Subjected to Very High C-Rate Discharge/Charge Operations. Int. J. Heat Mass Transf. 2018, 117, 261–272. [Google Scholar] [CrossRef]
  36. Haran, B.S.; Popov, B.N.; White, R.E. Determination of the Hydrogen Diffusion Coefficient in Metal Hydrides by Impedance Spectroscopy. J. Power Sources 1998, 75, 56–63. [Google Scholar] [CrossRef]
  37. Fuller, T.F.; Doyle, M.; Newman, J. Simulation and Optimization of the Dual Lithium Ion Insertion Cell. J. Electrochem. Soc. 1994, 141, 1–10. [Google Scholar] [CrossRef] [Green Version]
  38. Rahimian, S.K.; Rayman, S.; White, R.E. Comparison of Single Particle and Equivalent Circuit Analog Models for a Lithium-Ion Cell. J. Power Sources 2011, 196, 8450–8462. [Google Scholar] [CrossRef]
  39. Li, J.; Adewuyi, K.; Lotfi, N.; Landers, R.G.; Park, J. A Single Particle Model with Chemical/Mechanical Degradation Physics for Lithium Ion Battery State of Health (SOH) Estimation. Appl. Energy 2018, 212, 1178–1190. [Google Scholar] [CrossRef]
  40. Richardson, G.; Korotkin, I.; Ranom, R.; Castle, M.; Foster, J.M. Generalised Single Particle Models for High-Rate Operation of Graded Lithium-Ion Electrodes: Systematic Derivation and Validation. Electrochim. Acta 2020, 339, 135862. [Google Scholar] [CrossRef]
  41. Ren, L.; Zhu, G.; Kang, J.; Wang, J.V.; Luo, B.; Chen, C.; Xiang, K. An Algorithm for State of Charge Estimation Based on a Single-Particle Model. J. Energy Storage 2021, 39, 102644. [Google Scholar] [CrossRef]
  42. Ji, B.; Song, X.G.; Cao, W.P.; Pickert, V. Active Temperature Control of Li-Ion Batteries in Electric Vehicles. IET Conf. Publ. 2013, 2013, 3–7. [Google Scholar] [CrossRef]
  43. Zhu, C.; Lu, F.; Zhang, H.; Mi, C.C. Robust Predictive Battery Thermal Management Strategy for Connected and Automated Hybrid Electric Vehicles Based on Thermoelectric Parameter Uncertainty. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 6, 1796–1805. [Google Scholar] [CrossRef]
  44. Ma, Y.; Mou, H.; Zhao, H. Cooling Optimization Strategy for Lithium-Ion Batteries Based on Triple-Step Nonlinear Method. Energy 2020, 201, 117678. [Google Scholar] [CrossRef]
  45. Ma, Y.; Ding, H.; Mou, H.; Gao, J. Battery Thermal Management Strategy for Electric Vehicles Based on Nonlinear Model Predictive Control. Meas. J. Int. Meas. Confed. 2021, 186, 110115. [Google Scholar] [CrossRef]
  46. Xu, M.; Zhang, Z.; Wang, X.; Jia, L.; Yang, L. Two-Dimensional Electrochemical-Thermal Coupled Modeling of Cylindrical LiFePO4 Batteries. J. Power Sources 2014, 256, 233–243. [Google Scholar] [CrossRef]
  47. Kim, G.H.; Pesaran, A.; Spotnitz, R. A Three-Dimensional Thermal Abuse Model for Lithium-Ion Cells. J. Power Sources 2007, 170, 476–489. [Google Scholar] [CrossRef]
  48. Li, C.; Zhang, H.; Zhang, R.; Lin, Y.; Fang, H. On the Characteristics Analysis and Tab Design of an 18650 Type Cylindrical LiFePO4 Battery. Appl. Therm. Eng. 2021, 182, 116144. [Google Scholar] [CrossRef]
  49. Santhanagopalan, S.; Guo, Q.; Ramadass, P.; White, R.E. Review of Models for Predicting the Cycling Performance of Lithium Ion Batteries. J. Power Sources 2006, 156, 620–628. [Google Scholar] [CrossRef]
  50. Prada, E.; Di Domenico, D.; Creff, Y.; Bernard, J.; Sauvant-Moynot, V.; Huet, F. Simplified Electrochemical and Thermal Model of LiFePO4-Graphite Li-Ion Batteries for Fast Charge Applications. J. Electrochem. Soc. 2012, 159, A1508–A1519. [Google Scholar] [CrossRef]
  51. Zhang, H.; Li, C.; Chen, H. Analysis of Prismatic Lithium-Ion Battery Degradation Based on an Electrochemical-Thermal-Degradation Model. Int. J. Energy Res. 2022, 46, 23658–23681. [Google Scholar] [CrossRef]
  52. Prada, E.; Di Domenico, D.; Creff, Y.; Bernard, J.; Sauvant-Moynot, V.; Huet, F. A Simplified Electrochemical and Thermal Aging Model of LiFePO4-Graphite Li-Ion Batteries: Power and Capacity Fade Simulations. J. Electrochem. Soc. 2013, 160, A616–A628. [Google Scholar] [CrossRef] [Green Version]
  53. Ploehn, H.J.; Ramadass, P.; White, R.E. Solvent Diffusion Model for Aging of Lithium-Ion Battery Cells. J. Electrochem. Soc. 2004, 151, A456. [Google Scholar] [CrossRef] [Green Version]
  54. Jin, X.; Vora, A.; Hoshing, V.; Saha, T.; Shaver, G.; García, R.E.; Wasynczuk, O.; Varigonda, S. Physically-Based Reduced-Order Capacity Loss Model for Graphite Anodes in Li-Ion Battery Cells. J. Power Sources 2017, 342, 750–761. [Google Scholar] [CrossRef] [Green Version]
  55. Choi, Y.S.; Kang, D.M. Prediction of Thermal Behaviors of an Air-Cooled Lithium-Ion Battery System for Hybrid Electric Vehicles. J. Power Sources 2014, 270, 273–280. [Google Scholar] [CrossRef]
  56. Sahel, D.; Ameur, H.; Benzeguir, R.; Kamla, Y. Prediction of Heat Transfer Development in a Smooth Tube. J. Eng. Phys. Thermophys. 2018, 91, 682–687. [Google Scholar] [CrossRef]
  57. Zhao, H.; Gao, B.; Ren, B.; Chen, H. Integrated Control of In-Wheel Motor Electric Vehicles Using a Triple-Step Nonlinear Method. J. Frankl. Inst. 2015, 352, 519–540. [Google Scholar] [CrossRef]
  58. Wang, F.; Hao, N.; Song, L.; Chen, H. Triple-Step Nonlinear Control Design for Road Vehicles after a Tire Blow-out on the Highway. In Proceedings of the 2016 12th World Congress on Intelligent Control and Automation (WCICA), Guilin, China, 12–15 June 2016; pp. 1414–1419. [Google Scholar] [CrossRef]
  59. Li, X.; Zhang, X.; Song, S. Effect of Delayed Impulses on Input-to-State Stability of Nonlinear Systems. Automatica 2017, 76, 378–382. [Google Scholar] [CrossRef]
  60. Xu, M.; Wang, R.; Zhao, P.; Wang, X. Fast Charging Optimization for Lithium-Ion Batteries Based on Dynamic Programming Algorithm and Electrochemical-Thermal-Capacity Fade Coupled Model. J. Power Sources 2019, 438, 227015. [Google Scholar] [CrossRef]
Figure 1. Model diagram of the battery: (a) charge and discharge model; (b) electrochemical model (SP model); and (c) aging model (SEI layer formation on the negative electrode).
Figure 1. Model diagram of the battery: (a) charge and discharge model; (b) electrochemical model (SP model); and (c) aging model (SEI layer formation on the negative electrode).
Energies 16 00460 g001
Figure 2. Two thermal management modes of the 6s4p battery pack: (a) series cooling circuit and (b) parallel cooling circuit.
Figure 2. Two thermal management modes of the 6s4p battery pack: (a) series cooling circuit and (b) parallel cooling circuit.
Energies 16 00460 g002
Figure 3. Battery thermal management model and control strategies: (a) schematic of battery thermal management model used in the study; (b) simplified battery thermal management model for designing a triple-step nonlinear control; (c) coupled model of the battery cell; (d) triple-step nonlinear cooling control strategy; and (e) extended PID cooling control strategy.
Figure 3. Battery thermal management model and control strategies: (a) schematic of battery thermal management model used in the study; (b) simplified battery thermal management model for designing a triple-step nonlinear control; (c) coupled model of the battery cell; (d) triple-step nonlinear cooling control strategy; and (e) extended PID cooling control strategy.
Energies 16 00460 g003
Figure 4. Optimization process for the key battery.
Figure 4. Optimization process for the key battery.
Energies 16 00460 g004
Figure 5. Schematic diagram of the experimental setup.
Figure 5. Schematic diagram of the experimental setup.
Energies 16 00460 g005
Figure 6. Comparison of voltage, temperature, and capacity fading in simulations and experiments: (a) voltage comparison at different C-rates; (b) temperature comparison at different C-rates; and (c) comparison of decrease in capacity at different C-rates.
Figure 6. Comparison of voltage, temperature, and capacity fading in simulations and experiments: (a) voltage comparison at different C-rates; (b) temperature comparison at different C-rates; and (c) comparison of decrease in capacity at different C-rates.
Energies 16 00460 g006
Figure 7. Evolutions of the temperature and LOC of a single battery: (a) battery temperature at different C-rates; (b) LOC at different C-rates; (c) battery temperature at different T amb when C-rate is 5; (d) LOC at different T amb when C-rate is 5.
Figure 7. Evolutions of the temperature and LOC of a single battery: (a) battery temperature at different C-rates; (b) LOC at different C-rates; (c) battery temperature at different T amb when C-rate is 5; (d) LOC at different T amb when C-rate is 5.
Energies 16 00460 g007
Figure 8. The temperature, LOC, SEI layer thickness, and SEI layer porosity of a single battery: (a) battery temperature at different C-rates and T amb ; (b) LOC at different C-rates and T amb ; (c) SEI layer thickness at different C-rates and T amb ; (d) SEI layer porosity at different C-rates and T amb .
Figure 8. The temperature, LOC, SEI layer thickness, and SEI layer porosity of a single battery: (a) battery temperature at different C-rates and T amb ; (b) LOC at different C-rates and T amb ; (c) SEI layer thickness at different C-rates and T amb ; (d) SEI layer porosity at different C-rates and T amb .
Energies 16 00460 g008
Figure 9. Evolutions of the maximum temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack at different C-rates: (a) the maximum temperature; (b) the maximum LOC; (c) the maximum SEI layer thickness; (d) the maximum SEI layer porosity.
Figure 9. Evolutions of the maximum temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack at different C-rates: (a) the maximum temperature; (b) the maximum LOC; (c) the maximum SEI layer thickness; (d) the maximum SEI layer porosity.
Energies 16 00460 g009
Figure 10. The maximum and minimum temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack at different C-rates: (a) the maximum and minimum temperature; (b) the maximum and minimum LOC; (c) the maximum and minimum SEI layer thickness; (d) the maximum and minimum SEI layer porosity.
Figure 10. The maximum and minimum temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack at different C-rates: (a) the maximum and minimum temperature; (b) the maximum and minimum LOC; (c) the maximum and minimum SEI layer thickness; (d) the maximum and minimum SEI layer porosity.
Energies 16 00460 g010
Figure 11. Evolutions of the maximum temperature and maximum temperature difference of the cells in the 6s4p battery pack at different conditions: (a) the maximum temperature under parallel cooling; (b) the maximum temperature under series cooling; (c) the maximum temperature difference under parallel cooling; (d) the maximum temperature difference under series cooling.
Figure 11. Evolutions of the maximum temperature and maximum temperature difference of the cells in the 6s4p battery pack at different conditions: (a) the maximum temperature under parallel cooling; (b) the maximum temperature under series cooling; (c) the maximum temperature difference under parallel cooling; (d) the maximum temperature difference under series cooling.
Energies 16 00460 g011
Figure 12. The temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack under parallel cooling: (a) battery temperature; (b) LOC; (c) SEI layer thickness; (d) SEI layer porosity.
Figure 12. The temperature, LOC, SEI layer thickness, and SEI layer porosity of the cells in the 6s4p battery pack under parallel cooling: (a) battery temperature; (b) LOC; (c) SEI layer thickness; (d) SEI layer porosity.
Energies 16 00460 g012
Figure 13. Parameters evolutions of key batteries at 2 C-rate: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Figure 13. Parameters evolutions of key batteries at 2 C-rate: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Energies 16 00460 g013
Figure 14. Parameters evolutions of key batteries under NEDC: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Figure 14. Parameters evolutions of key batteries under NEDC: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Energies 16 00460 g014
Figure 15. Parameters evolutions of key batteries under US06: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Figure 15. Parameters evolutions of key batteries under US06: (a) heat generation of battery; (b) comparison of coolant velocity under extended PID and triple-step nonlinear control; (c) comparison of battery temperature under extended PID and triple-step nonlinear control; (d) comparison of LOC under extended PID and triple-step nonlinear control.
Energies 16 00460 g015
Table 1. The technical specifications and thermophysical parameters of the battery cell.
Table 1. The technical specifications and thermophysical parameters of the battery cell.
ParameterValue
Cell length140 mm
Cell width70 mm
Cell thickness7 mm
Nominal voltage3.2 V
Charge cut-off voltage3.65 V
Discharge cut-off voltage2 V
Nominal capacity5 A · h
Energy density245 Wh · L 1
Positive electrode materialLiFePO4
Negative electrode materialGraphite-based carbon
Electrolyte materialCarbonate based
Thermal conductivity, k b ( W · m 1 · K 1 )55.66
Density, ρ b ( kg · m 3 )2318.07
Specific thermal capacity, C p , b ( J · kg 1 · K 1 )1056.08
Table 2. Equations in the electrochemical model.
Table 2. Equations in the electrochemical model.
Physical and Chemical MechanismsEquationsBoundary Conditions
Solid phase: conservation of Li+ species t c s 1 r 2 r ( r 2 D s r c s ) = 0 D s r c s | r = 0 = 0 ,   D s r c s | r = R s = j f a s F
Electrolyte phase: conservation of Li+ species t ε e c e z ( D e eff z c e ) ( 1 t + ) j f F = 0 z c e | z = 0 = z c e | z = L = 0
Solid phase: charge conservation z ( σ eff z ϕ s ) j f = 0 σ e eff z ϕ s | z = 0 = σ eff z ϕ s | z = L = I A ele ,
z ϕ s | z = L NE = z ϕ s | z = L L PE = 0
Electrolyte phase: charge conservation z ( κ eff z ϕ e ) + z ( κ D eff z ln c e ) + j f = 0 z ϕ e | z = 0 = z ϕ e | z = L = 0
Electrochemical kinetics j f = a s i 0 [ exp ( α ox F R T η ) exp ( α red F R T η ) ]
Electrolyte ionic diffusivity D e eff = D e ε e γ e
Electrolyte ionic conductivity κ eff = κ ε e γ e
Electrolyte ionic diffusional conductivity κ D eff = 2 R T F κ eff ( t + 1 ) ( 1 + d ln f ± d ln c e )
Solid phase electronic conductivity σ eff = σ ε s γ s
Specific interfacial surface area a s = 3 ε s R s
Table 3. Equations solved for the aging model.
Table 3. Equations solved for the aging model.
Physical and Chemical MechanismsEquationsBoundary Conditions
SEI formation reaction kinetics j SEI = 2 F k SEI c solv partedg c s , n ( r = R s , n ) 2
exp ( β SEI F R T · ( δ SEI I κ SEI S n ) )
exp ( β SEI F R T ( ϕ s , n U SEI ) )
Solvent diffusion in the SEI c solv t = D solv 2 c solv y 2 d δ SEI d t c solv y D solv c solv y | y = R s , n d δ SEI d t c solv partedg = j SEI 2 F ,
c solv ( r = R s , n + δ SEI ) = ε SEI c solv b
Loss of capacity variation d C loss SEI d t = j SEI S n 3600
LOC SEI = C loss SEI C tot × 100
SEI porosity variation d ε SEI d t = j SEI S n 3600 C loss SEI ( 1 ε SEI ) ( 1 1 ε SEI 1 ε SEI * )
SEI layer thickness δ SEI = 3600 C loss SEI M SEI 2 F ρ SEI S n ( 1 ε SEI )
Table 4. Equations for the temperature-dependent parameters.
Table 4. Equations for the temperature-dependent parameters.
ParametersEquationsRefs.
Active material in solid phase
Solid phase diffusivity ( D s ) D s , p = 5.9 × 10 20 exp ( 50,000 8.314 ( 1 T 1 298 ) ) [50,52]
D s , n = 3 × 10 15 exp ( 35,000 8.314 ( 1 T 1 298 ) ) [50,52]
Solid phase conductivity ( σ ) σ p = 0.7 exp ( 10,000 8.314 ( 1 T 1 298 ) ) [50,52]
σ n = 100 exp ( 11,000 8.314 ( 1 T 1 298 ) ) [50,52]
Electrode reaction rate constant ( k 0 ) [50,52]
k 0 , n = 1.5 × 10 10 exp ( 20,000 8.314 ( 1 T 1 298 ) ) [50,52]
Electrolyte
Electrolyte ionic diffusivity ( D e ) D e = 1 × 10 4 × 10 4.43 54 T 229 exp ( c e 0.26795 T 2 + 169.99 T 25,342 ) [50,52]
Electrolyte ionic conductivity ( κ ) κ = ( 1.5069 × 10 2 + 8.8078 × 10 5 T 9.4814 × 10 8 T 2 ) c e exp ( ( 1.5799 × 10 3 + 7.8183 × 10 6 T 1.0513 × 10 8 T 2 ) c e 1.2167 ) [50,52]
SEI formation
SEI formation reaction rate constant ( k SEI ) k SEI = 2.9465 × 10 26 exp ( 149,150 8.314 ( 1 T 1 298 ) ) [50,52]
SEI solvent diffusivity ( D solv ) D solv = 5.5098 × 10 20 exp ( 79,781 8.314 ( 1 T 1 298 ) ) [50,52]
SEI conductivity ( κ SEI ) κ SEI = 1.7 × 10 4 exp ( 10,578 8.314 ( 1 T 1 298 ) ) [50,52]
Table 5. Parameters of the battery model.
Table 5. Parameters of the battery model.
ParametersValueRef.
SOC int (%)100Estimated
A ele ( m 2 )0.395[51]
C nom ( A · h )5[51]
α ox 0.5[51]
α red 0.5[51]
R ( J · mol 1 · K 1 )8.314
F ( C · mol 1 )96,487
Positive electrode
L PE ( μ m )70[51]
x p at SOC 0%0.74[51]
x p at SOC 100%0.035[51]
ε e , p 0.332[51]
ε s , p 0.42[51]
R s , p ( μ m )0.11[51]
γ s , p 2.1[51]
c s , p max ( mol · m 3 )22,806[51]
Negative electrode
S n ( m 2 )0.395[51]
L NE ( μ m )34[51]
x n at SOC 0%0.0132[51]
x n at SOC 100%0.811[51]
ε e , n 0.33[51]
ε s , n 0.555[51]
R s , n ( μ m )5[51]
γ s , n 2.3[51]
c s , n max ( mol · m 3 )31,370[51]
Aging due to SEI layer formation
C loss , int SEI ( A · h )0.01[52]
ε SEI , int 0.01[52]
U SEI (V)0.04[52]
ρ SEI ( kg · m 3 )1690[52]
c solv b ( mol · m 3 )4541[52]
β SEI 0.5[52]
γ SEI 1.5[52]
M SEI ( kg · mol 1 )0.162[52]
ε SEI * { ε SEI * = 0.03 I S n + 0.01   ( 0 I S n 6 ) ε SEI * = 0.1525 I S n 0.725   ( 6 < I S n 10 ) ε SEI * = 0.8   ( 10 < I S n 20 ) [52]
Electrolyte
c e ( mol · m 3 )1200[51]
t + 0.363[51]
Separator
L SEP ( μ m )25[51]
ε e , sep 0.54[51]
γ sep 1.5[51]
Table 6. Thermal physical parameters for the cold plate, coolant, and ambient air.
Table 6. Thermal physical parameters for the cold plate, coolant, and ambient air.
ParameterCold PlateCoolantAmbient Air
Initial temperature (K)298.15298.15298.15
Thermal conductivity,
k ( W · m 1 · K 1 )
2370.4156-
Density, ρ ( kg · m 3 )270010691.1691
Specific thermal capacity,
C p ( J · kg 1 · K 1 )
8973310-
Convective heat transfer coefficient,
h ( W · m 2 · K 1 )
--10
Dynamic viscosity, μ co ( kg · m 1 · s 1 )-0.004563
Table 7. Parameters for the two modes of thermal management system under different working conditions.
Table 7. Parameters for the two modes of thermal management system under different working conditions.
Cooling TypeParallelSeries
Pump rotary speed (r·min−1)1000300010003000
Mechanical power provided by the pump to the coolant (W)7.168218.8732.62081.552
Volumetric flow rate (L·min−1)18.325 (4.581)62.3352 (15.584)6.6547623.0198
Velocity (m·s−1)3.868 (0.967)13.153 (3.288)1.4044.858
Mass flow rate (kg·s−1)0.325 (0.081)1.104 (0.276)0.1180.408
C-rates
5 CCase ACase FCase KCase P
4 CCase BCase GCase LCase Q
3 CCase CCase HCase MCase R
2 CCase DCase ICase NCase S
1 CCase ECase JCase OCase T
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, H.; Chen, H.; Fang, H. Cooling Optimization Strategy for a 6s4p Lithium-Ion Battery Pack Based on Triple-Step Nonlinear Method. Energies 2023, 16, 460. https://doi.org/10.3390/en16010460

AMA Style

Zhang H, Chen H, Fang H. Cooling Optimization Strategy for a 6s4p Lithium-Ion Battery Pack Based on Triple-Step Nonlinear Method. Energies. 2023; 16(1):460. https://doi.org/10.3390/en16010460

Chicago/Turabian Style

Zhang, Hongya, Hao Chen, and Haisheng Fang. 2023. "Cooling Optimization Strategy for a 6s4p Lithium-Ion Battery Pack Based on Triple-Step Nonlinear Method" Energies 16, no. 1: 460. https://doi.org/10.3390/en16010460

APA Style

Zhang, H., Chen, H., & Fang, H. (2023). Cooling Optimization Strategy for a 6s4p Lithium-Ion Battery Pack Based on Triple-Step Nonlinear Method. Energies, 16(1), 460. https://doi.org/10.3390/en16010460

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop