Analysis of Heat Dissipation and Preheating Module for Vehicle Lithium Iron Phosphate Battery

: The ambient temperature has a great inﬂuence on the discharge and charging performance of a lithium battery, which may cause thermal runaway of the battery pack in extreme cases. In terms of the poor cooling effect caused by only using the cooling bottom plate for liquid cooling and the fact that the battery pack needs to be preheated before it can be used normally, a new cooling structure design was carried out, and a variety of cooling schemes and preheating schemes were proposed for analysis and comparison. The Star ccm+ simulation software was used to analyze and study their liquid cooling performance and preheating performance under different conditions. The best cooling scheme and preheating scheme were obtained by comparing the results of the simulation analysis. The simulation results show that the cooling performance of the cooling scheme using two vertical cooling plates and one cooling bottom plate is the best, and the preheating performance is best when the preheating liquid is used with a certain temperature ﬂow through the preheating pipe of the battery pack for a period of time, and then the battery pack is discharged until the battery pack temperature reaches the working temperature range. The research results have reference value for the control of the ambient temperature of a vehicle lithium iron phosphate battery.


Introduction
New energy vehicles powered by lithium batteries are gaining unprecedented attention and becoming more and more widely used [1]. However, the thermal runaway and damage to the immediate environment associated with this type of high-energy battery have become major safety concerns for electric vehicles [2,3].
Vehicle innovation is an inevitable trend of development under the multiple effects of energy shortage and technological innovation. At present, new energy vehicles are in line with this innovative development trend [4,5]. Its importance is self-evident because the power battery pack is a part of its power. However, a large amount of heat will be generated during the working process of the power battery pack. If the excess heat is not handled properly, it will cause overheating. The vehicle loses part of its efficiency and, in the worst-case scenario, the safety of passengers can be endangered. Therefore, dealing with this part of the heat is the main priority in research on the development of new energy vehicles [6,7]. Meanwhile, at low temperatures, there are restrictions when using lithium batteries: the viscosity of the internal electrolyte will increase, and the activity will decrease. Moreover, this will reduce the capacity of the lithium battery, so they cannot provide the best performance. At this time, thermal runaway or even safety accidents can occur due to the use of high current charging, so the preheating of the battery pack is also very important under low-temperature conditions.
The power battery is generally placed in the battery pack in a certain combination, and the battery pack has anti-collision, waterproof and anti-leakage properties [8]. The main cause of fire in most electric vehicles is thermal runaway caused by external impacts or improper control of the power battery [9]. For example, during spring and summer, electric vehicles are more likely to catch fire. This has caused changes in the performance of the power battery and led to abnormal chemical properties because of environmental factors such as long-term high temperatures or rain during these two seasons [10]. There are loopholes in the battery thermal management system, which causes the heat generated during driving to accumulate, and the temperature may be too high to cause a fire; during a collision with a foreign object, the original battery pack structure is deformed, resulting in leakage of the coolant or tilting of the battery, leading to fire [11]. These existing problems require the power battery pack to be designed to have good sealing performance in order to avoid external liquid infiltration into the battery pack to cause short circuits and avoid electrolyte contact with external oxygen; high hardness performance is also needed in order to avoid external collisions or sharp objects that cause the shell to rupture, as well as good heat insulation performance to avoid affecting the internal temperature of the battery pack when the ambient temperature is too high [12].
Recently, in order to prolong the cycle life and maximize the capability of power cells, a number of battery thermal management (BTM) strategies with external cooling (or heating), such as active cooling with air or liquid and passive cooling with a phase change material or heat pipe, have been developed [13][14][15]. With the recent increase in the battery pack energy density and load, air is unsuitable for use as cooling media in BTMS because of the poor heat capacity and low thermal conductivity [16]. PCM absorbs heat through its own material properties to reduce the temperature of the battery. However, depletion of the available latent heat under extreme conditions-such as high heat density during a high-current charge or discharge, or high ambient temperature-may cause failures of the thermal management systems [17]. Even under mild conditions (ambient temperature < 30 • C and discharge rate < 2 • C), a lack of external cooling may also lead to thermal management system failures [18]. Moreover, liquid-based systems have high heat transfer efficiency with the function of cooling or heating [19]. Therefore, all the research in this paper is based on the liquid cooling of BTMS.
To date, scholars have designed many schemes for heat dissipation in vehicle lithium batteries. However, most of their heat dissipation schemes use a single cooling base plate to dissipate heat, and the preheating model is complicated. There is still little in-depth research on preheating strategies in low-temperature environments. Therefore, the cooling structure presented in this paper adds a vertical cooling plate so that the cooling pipes are buried in the vertical cooling plate and the cooling bottom plate. The structure increases the cooling effect, and, on this basis, the cooling scheme with the best cooling performance is selected. In addition, a new preheating strategy is proposed for low-temperature environmental conditions, and the preheating pipeline is improved, which can effectively shorten the time to reach the battery temperature required for normal operation of the battery pack. At the same time, because the preheating scheme structure is the same as the cooling scheme structure, it can not only ensure that the battery pack can adapt to different working conditions but also saves costs to a certain extent.
In the second section of this paper, the establishment of the simulation model of a single battery module in Star ccm+ and the setting of related parameters are carried out, and then the simulation analysis is carried out under different environments and different discharge rate conditions. After this, in the third section, the design of the new cooling scheme based on liquid cooling is simulated and analyzed, and the conditions of discharge rate, coolant flow rate, number of cooling channels and number of cooling plates are changed to select the best cooling scheme; we then carry out response surface optimization analysis to further optimize the scheme.

Establishment of Single Battery Module Model
In this paper, a single battery module composed of prismatic lithium iron phosphate batteries is used for research and discussion. The size of the square lithium iron phosphate

Establishment of Single Battery Module Model
In this paper, a single battery module composed of prismatic lithium iron ph batteries is used for research and discussion. The size of the square lithium iron ph battery is 17 × 011 × 019 mm 3 , 18 square lithium iron phosphate composed of battery module. The space between individual cells is 1.5 mm. The schematic di shown in Figure 1. (1) The normal working temperature range of lithium iron phosphate batter viously larger than the normal working temperature range of lead-acid batte nickel-metal hydride batteries.
It is more flexible and can adapt to a wide range of actual working conditio (2) Lithium iron phosphate batteries have a larger number of charge and d cycles and a longer cycle life that is much higher than that of lead-acid batte nickel-hydrogen batteries and the service life of lithium iron phosphate batteries (3) A single lithium iron phosphate battery has a higher working voltag means that the use of a lithium iron phosphate battery can reduce the number of in the vehicle power battery pack to a certain extent [22].
(4) A Ni-MH battery has a high self-discharge rate. If the battery is not used f time, a large amount of battery energy will be lost, which will affect the service performance of the battery [23]. In comparison, it is less than 5% of the self-disch of lithium iron phosphate batteries. This means that the lithium iron phosphat can save energy for a long time. Furthermore, when it is not used for a long tim not lose a great deal of battery energy, increasing the additional costs [24].

Simplification of Simulation Model Conditions
In actual production and real-world application, in order to obtain the exact the heat-generating accurate battery total heat Q, the charge and discharge curre battery must be known, as well as polarization resistance and ohmic resistance, culating the total calorific value in this way will bring great calculation difficultie fore, an index of heat generation rate is used to deal with the total heat of lith phosphate batteries during discharge. Moreover, for the convenience of research culation, the heat generation of lithium iron phosphate batteries has been partiall fied. According to Bernardi's research hypothesis, we treat the heated battery as a source of heat. At this time, the heat generation rate of the battery can be expr Equation (1)  This paper selects prismatic lithium iron phosphate batteries for research, and the reasons are as follows: (1) The normal working temperature range of lithium iron phosphate batteries is obviously larger than the normal working temperature range of lead-acid batteries and nickel-metal hydride batteries.
It is more flexible and can adapt to a wide range of actual working conditions [20].
(2) Lithium iron phosphate batteries have a larger number of charge and discharge cycles and a longer cycle life that is much higher than that of lead-acid batteries and nickel-hydrogen batteries and the service life of lithium iron phosphate batteries [21].
(3) A single lithium iron phosphate battery has a higher working voltage, which means that the use of a lithium iron phosphate battery can reduce the number of batteries in the vehicle power battery pack to a certain extent [22].
(4) A Ni-MH battery has a high self-discharge rate. If the battery is not used for a long time, a large amount of battery energy will be lost, which will affect the service life and performance of the battery [23]. In comparison, it is less than 5% of the self-discharge rate of lithium iron phosphate batteries. This means that the lithium iron phosphate battery can save energy for a long time. Furthermore, when it is not used for a long time, it will not lose a great deal of battery energy, increasing the additional costs [24].

Simplification of Simulation Model Conditions
In actual production and real-world application, in order to obtain the exact value of the heat-generating accurate battery total heat Q, the charge and discharge current of the battery must be known, as well as polarization resistance and ohmic resistance, etc. Calculating the total calorific value in this way will bring great calculation difficulties. Therefore, an index of heat generation rate is used to deal with the total heat of lithium iron phosphate batteries during discharge. Moreover, for the convenience of research and calculation, the heat generation of lithium iron phosphate batteries has been partially simplified. According to Bernardi's research hypothesis, we treat the heated battery as a uniform source of heat. At this time, the heat generation rate q v of the battery can be expressed by Equation (1) [27]: where V battery represents the volume of the battery, m 3 ; E represents the electromotive force of the battery, V; U battery represents the terminal voltage across the battery, V; T represents the temperature of the battery, K; dE/dT represents the temperature influence coefficient, which is generally empirically constant in engineering applications.

Single Battery Module Simulation Parameter Setting
The physical characteristics of the prismatic lithium iron phosphate battery used in this paper are as follows.
(1) Density ρ The density of the single prismatic battery used in this paper is 2840 kg/m 3 .
(2) Specific heat capacity c The simplified specific heat capacity calculation equation can be expressed as the weighted average of the specific heat capacity of each part of the battery; the equation is as follows [28]: where c i is the specific heat capacity of the different components of the battery; m i is the corresponding quality; m a is the overall quality of the battery. According to the above equation, the specific heat capacity of the single square battery used in this paper can be derived as 1020 J/kg·K.
(3) Thermal Conductivity The thermal conductivity of a single prismatic battery is anisotropic, and the calculation equation is as follows [29]: where λ 1 , λ 2 , λ 3 represent the thermal conductivity of the positive electrode, negative electrode and separator of the battery, respectively; dx, dy, dz represent the thickness of different battery materials in these three directions, respectively; L, b, h represent the total length, thickness and width of the battery, respectively, and these parameters are provided by the battery manufacturer. According to the above equation, the thermal conductivity of the power battery in the three directions is 1.1, 18.3 and 18.3.

Discharge under Different Ambient Temperature Conditions
In the previous section, the parameters of a single battery module have been set, and the meshing quality also meets the research requirements. We input the above parameters into Star ccm+ to model and simulate a single battery module. Then, we discuss the temperature situation when a single battery module is discharged at the discharge rate of 1 C under different ambient temperatures. The research mainly discusses the maximum temperature and temperature distribution. Then, we plot the temperature of the battery module at 15 • C ambient temperature, 20 • C ambient temperature, 25 • C ambient temperature and 30 • C ambient temperature, respectively, in          Figure 2 shows the temperature of the battery module at the ambient temperature of 15℃. It can be seen from Figure 2 that the maximum temperature of the battery module will first rise at a relatively fast rate. At about 1000s, the rate of the maximum temperature rise will decrease. When at about 2800s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The temperature distribution law of the battery module is that the closer to the center of the module, the higher the temperature. This is because this part of the area is not in direct contact with the air that results in a poorer heat dissipation effect than the surface. The temperature difference of the module is 0.2K. Figure 3 shows the temperature of the battery module at the ambient temperature of 20℃. It can be seen from Figure 3 that the maximum temperature of the battery module will first rise at a relatively fast rate. At about 750s, the maximum temperature rises rate decreases. At about 2900s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.18K.   Figure 4 shows the temperature of the battery module at an ambient temperature of 25℃. The maximum temperature trend of the battery module is similar to that at the ambient temperature of 20℃. It can be seen from Figure 4 that the maximum temperature of the battery module will first rise at a relatively fast rate. At about 650s, the maximum temperature rises rate decreases. At about 2500s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.28K.         Figure 2 shows the temperature of the battery module at the ambient temperature of 15 °C. It can be seen from Figure 2 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 1000 s, the rate of the maximum temperature rise will decrease. When at approximately 2800 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The temperature distribution law of the battery module is that the closer to the center of the module, the higher the temperature. This is because this area is not in direct contact with the air, which results in a poorer heat dissipation effect than at the surface. The temperature difference in the module is 0.2 K. Figure 3 shows the temperature of the battery module at the ambient temperature of 20 °C. It can be seen from Figure 3 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 750 s, the maximum temperature rise rate decreases. At approximately 2900 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.18 K. Figure 4 shows the temperature of the battery module at an ambient temperature of 25 °C. The maximum temperature trend of the battery module is similar to that at the ambient temperature of 20 °C. It can be seen from Figure 4 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 650 s, the maximum temperature rise rate decreases. At approximately 2500 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.28 K. Figure 5 shows the temperature of the battery module at an environmental temperature of 30 °C. Figure 5 illustrates that the greatest temperature change in the battery module is similar to that at the environmental temperature of 20 °C and at the environmental temperature of 25 °C. However, at approximately 2300 s, the highest temperature of the battery module will rise at a relatively fast rate, and this rate is much faster than the rate under the first three environmental temperatures. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.34 K.
By comparing the data in Figures 2-5, it can be found that the maximum temperature variation in the individual battery modules is almost the same at different ambient temperatures, and the temperature distribution of a single module is also substantially similar.   Figure 2 shows the temperature of the battery module at the ambient temperature of 15 • C. It can be seen from Figure 2 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 1000 s, the rate of the maximum temperature rise will decrease. When at approximately 2800 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The temperature distribution law of the battery module is that the closer to the center of the module, the higher the temperature. This is because this area is not in direct contact with the air, which results in a poorer heat dissipation effect than at the surface. The temperature difference in the module is 0.2 K. Figure 3 shows the temperature of the battery module at the ambient temperature of 20 • C. It can be seen from Figure 3 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 750 s, the maximum temperature rise rate decreases. At approximately 2900 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.18 K. Figure 4 shows the temperature of the battery module at an ambient temperature of 25 • C. The maximum temperature trend of the battery module is similar to that at the ambient temperature of 20 • C. It can be seen from Figure 4 that the maximum temperature of the battery module will first rise at a relatively fast rate. At approximately 650 s, the maximum temperature rise rate decreases. At approximately 2500 s, the maximum temperature of the battery module will rise again at a relatively fast rate until the discharge stops. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.28 K. Figure 5 shows the temperature of the battery module at an environmental temperature of 30 • C. Figure 5 illustrates that the greatest temperature change in the battery module is similar to that at the environmental temperature of 20 • C and at the environmental temperature of 25 • C. However, at approximately 2300 s, the highest temperature of the battery module will rise at a relatively fast rate, and this rate is much faster than the rate under the first three environmental temperatures. The battery module temperature distribution law remains unchanged, and the temperature difference is 0.34 K.
By comparing the data in Figures 2-5, it can be found that the maximum temperature variation in the individual battery modules is almost the same at different ambient temperatures, and the temperature distribution of a single module is also substantially similar.

Discharge under Different Discharge Rates
This section studies the temperature of a single battery module under different discharge rates at room temperature. We evaluate a single module discharge at 1 C, 2 C, 3 C. The temperature of the module is shown in Figures 6-8.

Discharge under Different Discharge Rates
This section studies the temperature of a single battery module under different discharge rates at room temperature. We evaluate a single module discharge at 1C, 2C, 3C. The temperature of the module is shown in Figures 6   It can be seen from Figure 6 that the maximum temperature of the battery module first rises at a faster rate when a single battery module is discharged at a discharge rate of 1 C. Then, the rate of ascent slows down, and finally it rises at a faster rate. At 3600 s, the maximum temperature of the battery module can reach approximately 309.5 K. At this time, the temperature difference in the entire battery module is within 0.5 K. Moreover, it can be seen from the temperature distribution diagram of the battery module that the temperature of the battery at the head and tail of the module is relatively low. In comparison, the temperature of the battery in the middle of the module is higher, and the temperature in the area closer to the center is higher, while the temperature in the area farther It can be seen from Figure 6 that the maximum temperature of the battery module first rises at a faster rate when a single battery module is discharged at a discharge rate of 1 C. Then, the rate of ascent slows down, and finally it rises at a faster rate. At 3600 s, the maximum temperature of the battery module can reach approximately 309.5 K. At this time, the temperature difference in the entire battery module is within 0.5 K. Moreover, it can be seen from the temperature distribution diagram of the battery module that the temperature of the battery at the head and tail of the module is relatively low. In comparison, the temperature of the battery in the middle of the module is higher, and the temperature in the area closer to the center is higher, while the temperature in the area farther from the center is lower.
It can be seen from Figure 7 that after a single battery module is discharged at the discharge rate of 2 C for 1800 s, the maximum temperature can reach approximately 319 K, which is around 10 K higher than the maximum temperature when discharged at the rate of 1 C. It can be seen from Figure 8 that after a single battery module is discharged at the discharge rate of 3 C for 1200 s, the maximum temperature can reach approximately 335 K, which is around 25 K higher than the maximum temperature when discharged at the rate of 1 C. It can also be seen that the law of the battery module temperature distribution graph under 2 C and 3 C discharge rates is the same as that of 1 C.

Liquid-Cooled Cooling Model Establishment
Geometric modeling of liquid cooling models was performed using Solidworks. Figure 9 shows the liquid cooling model of the battery pack. The liquid cooling module is composed of two vertical cooling plates and one cooling bottom plate. The battery pack is evenly arranged on the cooling bottom plate and wrapped by two vertical cooling plates. The vertical cooling plate has a thickness of 10 mm and contains a right-angle cooling pipe with a diameter of 8 mm. The thickness of the cooling bottom plate is 30 mm, and it contains a snake-shaped cooling pipe with a diameter of 28 mm.

Discharge under Different Discharge Rates
In this section, Scheme Ⅰ is adopted, i.e., the cooling scheme of the plate containing one snake-shaped cooling channel and two vertical coo containing three cooling channels, to study the temperature of the battery ing by the coolant under the condition of different discharge rates. The coo and outlet of the cooling bottom plate are set to left-in and right-out. Figu temperature distribution results, and Figure 11 shows the maximum temp

Discharge under Different Discharge Rates
In this section, Scheme I is adopted, i.e., the cooling scheme of the cooling bottom plate containing one snake-shaped cooling channel and two vertical cooling plates each containing three cooling channels, to study the temperature of the battery pack after cooling by the coolant under the condition of different discharge rates. The cooling liquid inlet and outlet of the cooling bottom plate are set to left-in and right-out. Figure 10 shows the temperature distribution results, and Figure 11 shows the maximum temperature results.
It can be seen from Figures 10 and 11 that, in Scheme I, when the battery pack is discharged at the rate of 1 C and 2 C, the maximum temperature of the battery pack is approximately 304.7 K and 308 K, respectively, and the temperature difference in the battery pack is around 3.8 K and 5.3 K, respectively. This shows that the battery pack is discharged at these two rates, and the maximum temperature after cooling by the cooling module is ideal; the temperature uniformity is also relatively good. When the battery pack is discharged at the rate of 3 C and 4 C, the maximum temperature of the battery pack is approximately 316.7 K and 329.8 K, respectively, and the temperature difference in the Energies 2021, 14, 6196 9 of 25 battery pack is around 15.2 K and 25.8 K, respectively. In these two cases, the temperature of the battery pack is not ideal and needs further control. By comparing these four cases, the following conclusions can be obtained: as the discharge rate increases, the maximum temperature of the battery pack will also increase, the temperature uniformity will become worse, and the cooling effect of the liquid cooling module will become worse.
Bottom cooling plate Figure 9. Simplified geometric diagram of liquid-cooled battery module.

Discharge under Different Discharge Rates
In this section, Scheme Ⅰ is adopted, i.e., the cooling scheme of the cooling bottom plate containing one snake-shaped cooling channel and two vertical cooling plates each containing three cooling channels, to study the temperature of the battery pack after cooling by the coolant under the condition of different discharge rates. The cooling liquid inlet and outlet of the cooling bottom plate are set to left-in and right-out. Figure 10 shows the temperature distribution results, and Figure 11 shows the maximum temperature results.    It can be seen from Figures 10 and 11 that, in Scheme Ⅰ, when the battery pack is discharged at the rate of 1 C and 2 C, the maximum temperature of the battery pack is approximately 304.7K and 308K, respectively, and the temperature difference in the battery pack is around 3.8 K and 5.3 K, respectively. This shows that the battery pack is discharged at these two rates, and the maximum temperature after cooling by the cooling Figure 11. Maximum temperature of battery modules with different magnifications in Scheme I.

Discharge under Different Coolant Speed Conditions
After the research in the previous section, it is found that when the battery pack is discharged at a 1 C or 2 C rate, the maximum temperature does not exceed 308 K, and the temperature uniformity is also good. Therefore, we do not discuss these two situations further. In order to improve the cooling effect of the cooling module when the battery pack is discharged at a rate of 3 C and 4 C in the first solution, the second solution is adopted. Specifically, the cooling liquid inlet flow rate setting method is modified as follows: when the maximum temperature of the battery pack is lower than 308 K, the coolant inlet flow rate is 0.3 m/s, and when the maximum temperature of the battery pack is higher than 308 K, the coolant inlet flow rate is 0.8 m/s. The order of the inlet and outlet of the coolant remains unchanged. Figure 12 shows the temperature distribution results; Figure 13 shows the maximum temperature results.  In order to compare the effects of the two schemes more intuitively, we compare Scheme Ⅰ with Scheme Ⅱ. Figures 13 and 14 show the comparison results.   In order to compare the effects of the two schemes more intuitively, we compare Scheme Ⅰ with Scheme Ⅱ. Figures 13 and 14 show the comparison results. In order to compare the effects of the two schemes more intuitively, we compare Scheme I with Scheme II. Figures 13 and 14 show the comparison results. Figure 13 shows that, compared to the first solution, increasing the coolant inlet flow rate can effectively improve the cooling effect of the cooling module when the maximum temperature of the battery pack is higher than 308 K. It can also be seen that the maximum temperature of the battery pack in Scheme II is around 0.8 K lower than the maximum temperature of the battery pack in Scheme I when the battery pack is discharged at a Energies 2021, 14, 6196 11 of 25 discharge rate of 3 C. The maximum temperature of the battery pack in Scheme II is around 1.8 K lower than that in Scheme I when the battery pack is discharged at a discharge rate of 4 C. Moreover, the cooling effect in Scheme II is improved. The maximum temperature of the battery pack has been improved. It can be seen from Figure 14 that, compared to Scheme I, the temperature difference in the battery pack in Scheme II is reduced by approximately 0.4 K compared with the temperature difference in the battery pack in Scheme I when the battery pack is discharged at a discharge rate of 3 C. When the battery pack is discharged at 4 C, the temperature difference in the battery pack in Scheme II is reduced by around 2.5 K compared with the temperature difference in the battery pack in Scheme I, and the cooling effect in Scheme II is improved. The temperature distribution of the battery pack under the two different discharge rates has been improved. In summary, Scheme II can improve the cooling effect of the liquid cooling module. In order to achieve a better cooling effect that is in line with actual requirements, the second solution is further studied.
In order to compare the effects of the two schemes more intuitively, we compare Scheme Ⅰ with Scheme Ⅱ. Figures 13 and 14 show the comparison results.

Discharge under Different Numbers of Cooling Channels in the Vertical Cooling Plate
The liquid-cooled cooling modules adopted in Scheme I and Scheme II are composed of two vertical cooling plates and one cooling bottom plate. Two of the vertical cooling plates each contain three right-angle cooling channels, and the cooling bottom plate contains one snake-shaped cooling channel. In order to further study the influence of the number of different cooling channels on the temperature of the battery pack, this section proposes Scheme III and Scheme IV These two new schemes increase the number of right-angle cooling channels contained in two vertical cooling plates. The liquid-cooled cooling module of Scheme III includes two vertical cooling plates that each contain four right-angle cooling channels, and the cooling bottom plate contains one snake-shaped cooling channel. The liquid-cooled cooling module of Scheme IV includes two vertical cooling plates that each contain five right-angle cooling channels, and the cooling bottom plate contains one snakeshaped cooling channel, while the other conditions remain unchanged. The schematic diagram is shown in Figure 15.
From the research in the previous section, it is noted that increasing the coolant inlet flow rate can effectively improve the cooling effect of the cooling module, reduce the maximum temperature of the battery pack and improve the temperature uniformity. Therefore, in this section, an increased coolant inlet flow rate setting and coolant inlet and outlet sequence setting are adopted, and the remaining boundary conditions remain unchanged. The temperature distribution and maximum temperature of the battery pack under the 3 C and 4 C discharge rates of Scheme III and Scheme IV are studied.
In order to compare the maximum temperature of the battery pack more intuitively, this paper plots the maximum temperature of the battery pack for the four scenarios at 3 C and 4 C discharge rates in Figure 16.
tains one snake-shaped cooling channel. In order to further study the influence of the number of different cooling channels on the temperature of the battery pack, this section proposes Scheme Ⅲ and Scheme Ⅳ These two new schemes increase the number of rightangle cooling channels contained in two vertical cooling plates. The liquid-cooled cooling module of Scheme Ⅲ includes two vertical cooling plates that each contain four rightangle cooling channels, and the cooling bottom plate contains one snake-shaped cooling channel. The liquid-cooled cooling module of Scheme Ⅳ includes two vertical cooling plates that each contain five right-angle cooling channels, and the cooling bottom plate contains one snake-shaped cooling channel, while the other conditions remain unchanged. The schematic diagram is shown in Figure 15.
(a) scheme Ⅲ ( b) scheme Ⅳ From the research in the previous section, it is noted that increasing the coolant inlet flow rate can effectively improve the cooling effect of the cooling module, reduce the maximum temperature of the battery pack and improve the temperature uniformity. Therefore, in this section, an increased coolant inlet flow rate setting and coolant inlet and outlet sequence setting are adopted, and the remaining boundary conditions remain unchanged. The temperature distribution and maximum temperature of the battery pack under the 3 C and 4 C discharge rates of Scheme Ⅲ and Scheme Ⅳ are studied. In order to compare the maximum temperature of the battery pack more intuitively, this paper plots the maximum temperature of the battery pack for the four scenarios at 3 C and 4 C discharge rates in Figure 16. It can be seen from Figure 16 that increasing the number of right-angle cooling channels in the vertical cooling plate can effectively reduce the maximum temperature of the battery pack, and the more right-angle cooling channels contained in the vertical cooling plate, the better the maximum temperature control. The maximum temperature of the battery packs in Scheme Ⅰ and Scheme Ⅱ continues to increase and finally stabilizes at a certain value. In comparison, Scheme Ⅲ and Scheme Ⅳ can reduce the maximum temperature of the battery pack after reaching the peak value. This means that the cooling performance of these two schemes is better than that of Scheme Ⅰand Scheme Ⅱ, and the cooling effect of Scheme Ⅳ is better than that of Scheme Ⅲ. However, it can also be seen that the cooling performance of Scheme Ⅳ is only slightly better than that of Scheme Ⅲ. In other words, increasing the number of cooling pipes will no longer have a significant impact on the improvement in cooling performance when the number of cooling pipes in the vertical cooling plate is greater than five. Thus, it will not be discussed further when the number of cooling pipes in the vertical cooling plate is more than five.

Discharge under Different Numbers of Cooling Channels in the Cooling Bottom Plate
Scheme Ⅴ consists of two vertical cooling plates and one cooling bottom plate. Both vertical cooling plates contain three right-angle cooling pipes, and the cooling bottom plate contains two snake-shaped pipes. The cooling module of Scheme Ⅵ includes two vertical cooling plates and one cooling bottom plate. Both vertical cooling plates contain four right-angle pipes, and the cooling bottom plate contains two snake-shaped cooling pipes. The cooling module of Scheme Ⅶ consists of two vertical cooling plates and one cooling bottom plate. The two vertical cooling plates contain five right-angle cooling It can be seen from Figure 16 that increasing the number of right-angle cooling channels in the vertical cooling plate can effectively reduce the maximum temperature of the battery pack, and the more right-angle cooling channels contained in the vertical cooling plate, the better the maximum temperature control. The maximum temperature of the battery packs in Scheme I and Scheme II continues to increase and finally stabilizes at a certain value. In comparison, Scheme III and Scheme IV can reduce the maximum temperature of the battery pack after reaching the peak value. This means that the cooling performance of these two schemes is better than that of Scheme Iand Scheme II, and the cooling effect of Scheme IV is better than that of Scheme III. However, it can also be seen that the cooling performance of Scheme IV is only slightly better than that of Scheme III. In other words, increasing the number of cooling pipes will no longer have a significant impact on the improvement in cooling performance when the number of cooling pipes in the vertical cooling plate is greater than five. Thus, it will not be discussed further when the number of cooling pipes in the vertical cooling plate is more than five.

Discharge under Different Numbers of Cooling Channels in the Cooling Bottom Plate
Scheme V consists of two vertical cooling plates and one cooling bottom plate. Both vertical cooling plates contain three right-angle cooling pipes, and the cooling bottom plate contains two snake-shaped pipes. The cooling module of Scheme VI includes two vertical cooling plates and one cooling bottom plate. Both vertical cooling plates contain four right-angle pipes, and the cooling bottom plate contains two snake-shaped cooling pipes. The cooling module of Scheme VII consists of two vertical cooling plates and one cooling bottom plate. The two vertical cooling plates contain five right-angle cooling pipes, and the cooling bottom plate contains two snake-shaped cooling pipes. The schematic diagram is shown in Figure 17.
The remaining boundary conditions and the order of the coolant inlet and outlet remain unchanged, and a larger coolant inlet flow rate setting is selected. The temperature distribution and maximum temperature of the battery pack under the 3 C and 4 C discharge rates of these three schemes are studied and the temperature distribution and maximum temperature of Scheme IV are compared to those of Scheme VII. Figure 18 shows the maximum temperature comparison of these four schemes.  It can be seen from Figure 18 that there is almost no significant difference in the cooling effect of Scheme Ⅳ in the interval of 0-500 s when the battery pack is discharged at a rate of 3 C. The maximum temperature of the battery pack continues to rise and reaches a peak. The maximum temperature of the battery pack begins to decrease under the action of the cooling module in the interval of 500-1200 s. Moreover, compared to Scheme Ⅴ and  It can be seen from Figure 18 that there is almost no significant difference in the cooling effect of Scheme Ⅳ in the interval of 0-500 s when the battery pack is discharged at a rate of 3 C. The maximum temperature of the battery pack continues to rise and reaches a peak. The maximum temperature of the battery pack begins to decrease under the action of the cooling module in the interval of 500-1200 s. Moreover, compared to Scheme Ⅴ and It can be seen from Figure 18 that there is almost no significant difference in the cooling effect of Scheme IV in the interval of 0-500 s when the battery pack is discharged at a rate of 3 C. The maximum temperature of the battery pack continues to rise and reaches a peak. The maximum temperature of the battery pack begins to decrease under the action of the cooling module in the interval of 500-1200 s. Moreover, compared to Scheme V and Scheme V, the cooling effect of Scheme VI and Scheme VII is better. Moreover, the maximum temperature curves of these two schemes almost overlap, and it can be considered that the cooling effects of these two schemes are almost the same.

Response Surface Optimization
The test data adopt the recorded value of the Starccm+ simulation result to observe the influence of the cooling pipe diameter, the temperature of the coolant and the thermal conductivity of the material on the temperature field of the power battery pack. When using the BBD test method in the Design Expert software for the design, we set the cooling pipe diameter, coolant temperature and material thermal conductivity as three independent variables, which are represented by A, B and C. The specific content is shown in Table 1 below. After determining the independent variables and factor levels of the scheme, we then use the temperature equilibrium degree M of the power battery pack as the dependent variable. We enter the value of the dependent variable under different independent variables in the simulation results obtained by Starccm+ into Design Expert, as shown in Table 2 below. The regression equation model is as follows: In the simulation analysis of the Design Expert software, the regression Equation (4) is used to solve the maximum value, and the following data are obtained: when the diameter of the cooling pipe is 7.5 mm, the coolant temperature is 10.888 • C, and the thermal conductivity of the material is 4.956 W/(m·K); thus, the predicted value of the temperature equilibrium M of the power battery pack is 97.512%, which is the maximum value of these data.
In this section, the diameter of the cooling channel, coolant temperature and thermal conductivity of the material are the three independent variables. Table 3 shows the results of the regression model analysis of variance with the temperature equilibrium degree M as the dependent variable. It can be seen from Table 4 that the p value of the model is less than 0.0001, which means that the model is extremely significant. The correlation coefficient R 2 is 0.9939, and R 2 adj is 0.9878. The values of these two coefficients indicate that the regression model designed by the response surface method is effective. It can also be seen that the order of influence of the three single factors on the temperature equilibrium degree M of the power battery pack is A > B > C. The diameter of the cooling pipe has the greatest influence, the temperature of the coolant has a moderate influence, and the thermal conductivity of the material has the least influence.  Figure 19 shows the normal distribution of residuals and a comparison between residuals and experimental predictions.
influence of the three single factors on the temperature equilibrium degree M of the power battery pack is A > B > C. The diameter of the cooling pipe has the greatest influence, the temperature of the coolant has a moderate influence, and the thermal conductivity of the material has the least influence. Figure 19 shows the normal distribution of residuals and a comparison between residuals and experimental predictions.

Battery Pack Model Preheating Analysis
Electric vehicles not only need to cool down the battery pack while driving, but also need to consider their working conditions in cold areas. When starting directly under low-temperature conditions, due to the reduction in the voltage platform and discharge capacity of the lithium-ion battery, the charging and discharging speed of the battery pack will be reduced, in addition to the problem of battery capacity attenuation. Therefore, the power battery pack must be preheated for electric vehicles under low-temperature working conditions. This section adopts the six geometric structures of the optimal cooling performance among the seven cooling schemes proposed in Section 3 and conducts the preheating research of the battery pack on this basis. By studying the preheating of the battery pack under different conditions of the preheating model, a preheating plan with better preheating performance is obtained. The structure optimization and response surface analysis of the preheating scheme are carried out, and the preheating scheme with the best temperature equilibrium degree is obtained.

Different Preheating Fluid Properties
This paper first adopts the following materials: the preheating liquid in the preheating pipeline is mixed with 50% ethylene glycol and 50% aqueous solution, and the specific material property parameters of 50% water and 50% ethylene glycol solution are shown in Table 4. The temperature of the hydrothermal fluid is 35 • C. The materials of the two vertical plates and the bottom plate remain unchanged; the battery materials are also unchanged. The ambient temperature is set to −15 • C. The heat exchange method between the two vertical plates, the bottom plate and the surface of the battery pack and the external environment are environmental heat exchange. The heat transfer coefficient is 5 W/ m 2 ·K , and the temperature is consistent with the ambient temperature. The inlet of the preheating pipe is set as the speed inlet, the outlet is set as the pressure outlet, and the initial velocity of the preheating liquid is set to 0.5 m/s. Figure 20 shows the preheating performance of the battery packs using the three different concentrations of preheating solution at an ambient temperature of −15 • C. It can be seen that the 50% concentration of mixed solution is used as the preheating solution. It has better preheating performance. Compared with the 70% mixed solution, when the 50% concentration mixed solution is used as the preheating solution, the temperature of the lowest-temperature part of the battery pack can be increased by approximately 3 K.  Figure 21 shows the temperature distribution of the battery pack when a mixed solution of different concentrations is used as the preheating liquid. It can be seen that the temperature distribution trend of the battery pack after using the mixed solution of three different concentrations as the preheating liquid is similar. The temperature of the area near the preheating pipe is higher, and the temperature of the central area and the upper surface of the battery part is lower. This is because when the preheating liquid passes through the preheating pipe, the battery part closer to the pipe absorbs more heat. Figure  22 plots the temperature difference in order to better compare the temperature difference among the three schemes.     Figure 21 shows the temperature distribution of the battery pack when a mixed solution of different concentrations is used as the preheating liquid. It can be seen that the temperature distribution trend of the battery pack after using the mixed solution of three different concentrations as the preheating liquid is similar. The temperature of the area near the preheating pipe is higher, and the temperature of the central area and the upper surface of the battery part is lower. This is because when the preheating liquid passes through the preheating pipe, the battery part closer to the pipe absorbs more heat. Figure 22 plots the temperature difference in order to better compare the temperature difference among the three schemes.  Figure 21 shows the temperature distribution of the battery pack when a mixed solution of different concentrations is used as the preheating liquid. It can be seen that the temperature distribution trend of the battery pack after using the mixed solution of three different concentrations as the preheating liquid is similar. The temperature of the area near the preheating pipe is higher, and the temperature of the central area and the upper surface of the battery part is lower. This is because when the preheating liquid passes through the preheating pipe, the battery part closer to the pipe absorbs more heat. Figure  22 plots the temperature difference in order to better compare the temperature difference among the three schemes.      Figure 21 shows the temperature distribution of the battery pack when a mixed solution of different concentrations is used as the preheating liquid. It can be seen that the temperature distribution trend of the battery pack after using the mixed solution of three different concentrations as the preheating liquid is similar. The temperature of the area near the preheating pipe is higher, and the temperature of the central area and the upper surface of the battery part is lower. This is because when the preheating liquid passes through the preheating pipe, the battery part closer to the pipe absorbs more heat. Figure  22 plots the temperature difference in order to better compare the temperature difference among the three schemes.     It can be clearly seen from the results in Figure 22 that when mixed solutions of different concentrations are used as the preheating solution, the temperature uniformity is slightly different. Among them, the temperature uniformity of the 50% mixed solution is the best, and the temperature uniformity of the 70% mixed solution is the worst. However, the difference between the maximum temperature and the minimum temperature of the battery pack and the difference between the maximum temperature and the minimum temperature of the 50% mixed solution battery pack are not large, and the difference is around 1 K. Figure 23 shows the temperature of the battery pack after preheating when using a 50% mixed solution as the preheating solution under different ambient temperatures. It can be seen that when the ambient temperature is −15 • C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1 K and 2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 • C and −25 • C. When the ambient temperature is −15 • C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 7.54 K; when the ambient temperature is −20 • C, this value is 8.332 K; when the ambient temperature is −25 • C, this value is 9.089 K. It can be clearly seen from the results in Figure 22 that when mixed solutions of different concentrations are used as the preheating solution, the temperature uniformity is slightly different. Among them, the temperature uniformity of the 50% mixed solution is the best, and the temperature uniformity of the 70% mixed solution is the worst. However, the difference between the maximum temperature and the minimum temperature of the battery pack and the difference between the maximum temperature and the minimum temperature of the 50% mixed solution battery pack are not large, and the difference is around 1 K. Figure 23 shows the temperature of the battery pack after preheating when using a 50% mixed solution as the preheating solution under different ambient temperatures. It can be seen that when the ambient temperature is −15 °C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1 K and 2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 °C and −25 °C. When the ambient temperature is −15 °C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 7.54 K; when the ambient temperature is −20 °C, this value is 8.332 K; when the ambient temperature is −25 °C, this value is 9.089 K.   Figure 24 shows the temperature of the battery pack after preheating when 60% mixed solution is used as the preheating liquid under different ambient temperatures. It can be seen that when the ambient temperature is −15 °C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1.2 K and 2.2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 °C and −25 °C. When the ambient temperature is −15 °C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 8.103 K; when the ambient temperature is −20 °C, this value is 8.913 K; when the ambient temperature is −25 °C, this value is 9.724 K.  Figure 24 shows the temperature of the battery pack after preheating when 60% mixed solution is used as the preheating liquid under different ambient temperatures. It can be seen that when the ambient temperature is −15 • C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1.2 K and 2.2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 • C and −25 • C. When the ambient temperature is −15 • C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 8.103 K; when the ambient temperature is −20 • C, this value is 8.913 K; when the ambient temperature is −25 • C, this value is 9.724 K. Figure 25 shows the temperature of the battery pack after preheating when 70% mixed solution is used as the preheating liquid under different ambient temperatures. It can be seen that when the ambient temperature is −15 • C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1.3 K and 2.2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 • C and −25 • C. When the ambient temperature is −15 • C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 8.703 K; when the ambient temperature is −20 • C, this value is 9.573 K; when the ambient temperature is −25 • C, this value is 10.444 K.   Figure 25 shows the temperature of the battery pack after preheating when 70% mixed solution is used as the preheating liquid under different ambient temperatures. It can be seen that when the ambient temperature is −15 °C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1.3 K and 2.2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 °C and −25 °C. When the ambient temperature is −15 °C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 8.703 K; when the ambient temperature is −20 °C, this value is 9.573 K; when the ambient temperature is −25 °C, this value is 10.444 K.   Figures 23-25 show the preheating performance of 50%, 60% and 70% mixed solutions, respectively, at different ambient temperatures. It can be seen that for the mixed solutions of these three concentrations, the following laws are met: the lower the ambient temperature, the lower the temperature that the power battery pack can reach after preheating, and the lower the ambient temperature, the worse the temperature uniformity of the power battery pack. In order to better compare the preheating performance, the preheating performance of mixed solutions of different concentrations at different ambient temperatures is presented as a broken line graph in Figures 26 and 27 for comparison.

Different Ambient Temperature Conditions
It can be seen from Figure 26 that when the ambient temperature is −20 °C, the battery pack with 50% concentration mixed solution as the preheating liquid has the best preheating performance.
Compared with the preheating temperature of the battery pack with 60% concentration mixed solution and the preheating temperature of the battery pack with 70% concentration mixed solution, the battery pack using the 50% concentration mixed solution as    Figure 25 shows the temperature of the battery pack after preheating when 70% mixed solution is used as the preheating liquid under different ambient temperatures. It can be seen that when the ambient temperature is −15 °C, the preheated temperature of the battery pack using the mixed solution of this concentration is approximately 1.3 K and 2.2 K higher than the preheating temperature of the battery pack when the ambient temperature is −20 °C and −25 °C. When the ambient temperature is −15 °C, the difference between the maximum temperature and the minimum temperature of the battery pack after preheating is 8.703 K; when the ambient temperature is −20 °C, this value is 9.573 K; when the ambient temperature is −25 °C, this value is 10.444 K.  show the preheating performance of 50%, 60% and 70% mixed solutions, respectively, at different ambient temperatures. It can be seen that for the mixed solutions of these three concentrations, the following laws are met: the lower the ambient temperature, the lower the temperature that the power battery pack can reach after preheating, and the lower the ambient temperature, the worse the temperature uniformity of the power battery pack. In order to better compare the preheating performance, the preheating performance of mixed solutions of different concentrations at different ambient temperatures is presented as a broken line graph in Figures 26 and 27 for comparison.
It can be seen from Figure 26 that when the ambient temperature is −20 °C, the battery pack with 50% concentration mixed solution as the preheating liquid has the best preheating performance.
Compared with the preheating temperature of the battery pack with 60% concentration mixed solution and the preheating temperature of the battery pack with 70% concentration mixed solution, the battery pack using the 50% concentration mixed solution as  Figures 23-25 show the preheating performance of 50%, 60% and 70% mixed solutions, respectively, at different ambient temperatures. It can be seen that for the mixed solutions of these three concentrations, the following laws are met: the lower the ambient temperature, the lower the temperature that the power battery pack can reach after preheating, and the lower the ambient temperature, the worse the temperature uniformity of the power battery pack. In order to better compare the preheating performance, the preheating performance of mixed solutions of different concentrations at different ambient temperatures is presented as a broken line graph in Figures 26 and 27 for comparison. It can be seen from Figure 27 that when the ambient temperature is −25 °C, the 50% mixed solution has the best preheating performance. Compared with 60% and 70% mixed solutions, the preheating temperature can be increased by approximately 1.3 K and 2.8 K, respectively.  It can be seen from Figure 27 that when the ambient temperature is −25 °C, the 50% mixed solution has the best preheating performance. Compared with 60% and 70% mixed solutions, the preheating temperature can be increased by approximately 1.3 K and 2.8 K respectively.

Different Preheating Strategies
The preheating strategy of preheating the power battery pack with a mixed solution at 35 °C is adopted, and then the power battery pack is discharged for a period of time. In this way, we compare the battery pack's preheating performance at different ambient tem peratures. In order to better compare the two preheating strategies, we take a 50% con centration mixed solution as an example and presented a dotted line graph, shown in Figure 28, for comparison.
It can be seen from Figure 29a that if we only use a preheating method of 50% mixed solution with a certain temperature flowing through the entire power battery pack when the ambient temperature is −15 °C and the power battery pack is preheated with a 50% mixed solution, it will take around 700 s to preheat the temperature of the power battery pack to 285 K. Using the above conditions, it takes around 800 s to preheat the power battery pack to 290 K, and around 1100 s to preheat the power battery pack to 295 K However, the 50% mixed solution with a certain temperature is first used to flow through the entire power battery pack, and then the preheating method of discharging the battery pack takes around 650 s to preheat to 285 K, around 680 s to preheat to 290 K and around 850 s to preheat to 295 K. By comparison, the second preheating strategy has better pre heating performance at an ambient temperature of −15 °C. When the target temperature It can be seen from Figure 26 that when the ambient temperature is −20 • C, the battery pack with 50% concentration mixed solution as the preheating liquid has the best preheating performance.
Compared with the preheating temperature of the battery pack with 60% concentration mixed solution and the preheating temperature of the battery pack with 70% concentration mixed solution, the battery pack using the 50% concentration mixed solution as the preheating liquid can increase the preheating temperature by around 1.5 K and 3.5 K, respectively.
It can be seen from Figure 27 that when the ambient temperature is −25 • C, the 50% mixed solution has the best preheating performance. Compared with 60% and 70% mixed solutions, the preheating temperature can be increased by approximately 1.3 K and 2.8 K, respectively.

Different Preheating Strategies
The preheating strategy of preheating the power battery pack with a mixed solution at 35 • C is adopted, and then the power battery pack is discharged for a period of time.
In this way, we compare the battery pack's preheating performance at different ambient temperatures. In order to better compare the two preheating strategies, we take a 50% concentration mixed solution as an example and presented a dotted line graph, shown in Figure 28, for comparison.
It can be seen from Figure 29a that if we only use a preheating method of 50% mixed solution with a certain temperature flowing through the entire power battery pack when the ambient temperature is −15 • C and the power battery pack is preheated with a 50% mixed solution, it will take around 700 s to preheat the temperature of the power battery pack to 285 K. Using the above conditions, it takes around 800 s to preheat the power battery pack to 290 K, and around 1100 s to preheat the power battery pack to 295 K. However, the 50% mixed solution with a certain temperature is first used to flow through the entire power battery pack, and then the preheating method of discharging the battery pack takes around 650 s to preheat to 285 K, around 680 s to preheat to 290 K and around 850 s to preheat to 295 K. By comparison, the second preheating strategy has better preheating performance at an ambient temperature of −15 • C. When the target temperature is 285 K, it can save approximately 7.1% of the preheating time. When the target temperature is 290 K, it can save approximately 15% of the preheating time. Moreover, when the target temperature is 295 K, it can save approximately 22.7% of the preheating time, which effectively improves the preheating performance.
It can be seen from Figure 28b that when the ambient temperature is −20 • C, the first preheating strategy takes around 730 s to preheat the power battery pack to 285 K, around 900 s to preheat to 290 K and around 1050 s to preheat to 295 K. However, the time required for the second preheating strategy to preheat to the same target temperature is 700 s, 750 and 870 s, respectively. It can be seen that the second preheating strategy can save approximately 4.3% of the preheating time when the target temperature is 285 K, approximately 16.6% of the preheating time when the target temperature is 290 K and approximately 26.7% of the preheating time when the target temperature is 290 K.   It can be seen from Figure 28b that when the ambient temperature is −20 °C, the first preheating strategy takes around 730 s to preheat the power battery pack to 285 K, around 900 s to preheat to 290 K and around 1050 s to preheat to 295 K. However, the time required for the second preheating strategy to preheat to the same target temperature is 700    It can be seen from Figure 28b that when the ambient temperature is −20 °C, the first preheating strategy takes around 730 s to preheat the power battery pack to 285 K, around 900 s to preheat to 290 K and around 1050 s to preheat to 295 K. However, the time required for the second preheating strategy to preheat to the same target temperature is 700 It can be seen from Figure 28c that when the ambient temperature is −25 • C, the first preheating strategy takes around 650 s to preheat the power battery pack to 280 K, around 750 s to preheat to 285 K and around 950 s to preheat to 290 K. However, the time required for the second preheating strategy to preheat to the same target temperature is 600 s, 680 and 800 s, respectively. It can be seen that the second preheating strategy can save approximately 7.6% of the preheating time when the target temperature is 285 K, approximately 9.3% of the preheating time when the target temperature is 290 K and approximately 15.8% of the preheating time when the target temperature is 290 K.
After the above comparison of the two different preheating strategies, the following conclusions can be drawn: when the power battery pack is preheated with a 50% mixed solution, the battery pack temperature reaches a certain value and then discharges it and continues to be heated to a temperature suitable for work. However, this preheating strategy is less effective than the preheating strategy that only uses the mixed solution for preheating, and the lower the ambient temperature, the longer the preheating time required.

Response Surface Optimization
Response surface optimization is mainly used to investigate the influence of the three factors of the diameter of the preheating pipe, the temperature of the preheating fluid and the thermal conductivity of the material on the temperature field of the power battery pack and obtain the optimal value. The test data for optimization adopt the recorded value of the Star ccm+ simulation result. When using the BBD test method in the Design Expert software, we set the diameter of the preheating pipe, the temperature of the preheating fluid and the thermal conductivity of the material as three independent variables, represented by the letters A, B and C, respectively. The details are shown in Table 5 below.
After determining the independent variables and factor levels of the plan, we then use the temperature equilibrium degree M of the power battery pack as the dependent variable and enter the value of the dependent variable under different independent variables in the simulation results obtained by Star ccm+ into Design Expert. The output results are shown in Table 6 below. The regression equation model is as follows: M = 96.23 + 0.0992A − 0.312B − 0.0119C − 0.045AB + 0.0075AC− 0.0153BC + 1.28A 2 − 0.0102B 2 (5) In Design Expert, the regression Equation (5) is solved for the maximum value. When the diameter of the preheating tube is 8.054 mm, the temperature of the preheating liquid is 32.394 • C, and the thermal conductivity of the material is 4.976 W/(m·K). At this time, the predicted value of the temperature equilibrium degree M of the power battery pack is 96.061%, which is the optimal value in these data. This section takes the diameter of the preheating pipe, the temperature of the preheating liquid and the thermal conductivity of the material as the three independent variables, and the regression model variance analysis results with the temperature equilibrium degree M as the dependent variable are shown in Table 7. It can be seen from Table 7 that the p value of the model is less than 0.0001, which means that the model is extremely significant. The correlation coefficient R 2 is 0.9608, and R 2 adj is 0.9213. The values of these two coefficients indicate that the regression model designed by the response surface method is effective. The p value of A 2 is less than 0.001, indicating that the influence of this factor is extremely significant. It can also be seen that the order of influence of the three single factors on the temperature equilibrium degree M of the power battery pack is B > A > C. Therefore, the influence of the temperature of the preheating liquid is the largest, that of the diameter of the preheating pipe is moderate and that of the thermal conductivity of the material is the smallest. Figure 29 shows a normal distribution diagram of the residual and a comparison diagram of the residual values and the experimentally predicted values.

Conclusions
This paper describes the thermal characteristic analysis and simulation discharge analysis of the square lithium iron phosphate battery. The simulation results show that under the same ambient temperature, the greater the discharge rate of the battery, the greater the heat generation, the greater the temperature rise and the more uneven the temperature distribution.
Liquid cooling is adopted to solve the problem wherein air cooling cannot meet the cooling requirements of the battery pack when the battery pack is discharged at a high rate, and the cooling process is simulated. The simulation results show that changing the coolant inlet flow rate can effectively improve the cooling performance of the battery pack under high discharge rate conditions; increasing the number of cooling pipes contained in Changing the order of the inlet and outlet of the coolant can improve the cooling performance. The response surface analysis is performed on the preferred cooling scheme with the best cooling effect, and the best temperature equilibrium is obtained.
A preheating model for a lithium iron phosphate battery is proposed in order to avoid thermal runaway during low-temperature battery charging, and the preheating process is simulated. The simulation results show that a 50% concentration of water and ethylene glycol mixed solution has the best preheating performance. Adopting the strategy of staged preheating and interspersed discharging can shorten the preheating time of the battery pack, and changing the sequence of the inlet and outlet of the preheating liquid can improve the preheating performance of the battery pack. The response surface analysis method is adopted to optimize the parameters of the preheating model, and the best temperature equilibrium is obtained.
The newly designed heat dissipation structure, cooling scheme and preheating scheme presented in this paper can be used in the new energy vehicle lithium battery heat dissipation structure module to keep the vehicle lithium battery temperature within a safe temperature in order to avoid thermal runaway. However, due to limitations of time and the authors' expertise, this paper still has many shortcomings. Given the current research status of this topic at home and abroad, and the shortcomings of this paper, the authors believe that it can be further researched and improved in the following aspects: the temperature of the battery pack under actual vehicle driving conditions should be considered, such as the discharge of the battery pack and the temperature of the battery pack under various common operating conditions. Moreover, when the calculation conditions allow, the heat dissipation of the overall liquid-cooled battery cooling system and the preheating of the preheating battery pack should be considered, and the module model should be simplified or not simplified as much as possible. The overall structure of the battery pack should be designed. In addition to this, the compressive strength of the geometric structure of the battery pack should also be considered.