 Open Access
 Authors : Sherwin Rajkumar , Sanju Nagendran , Princeton Benedict , James Deepak
 Paper ID : IJERTV9IS080323
 Volume & Issue : Volume 09, Issue 08 (August 2020)
 Published (First Online): 02092020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Determination of Thermal Resistances and Optimum Cooling Mode for an Electronic Module used in Electric and Hybrid Vehicles
Sherwin Rajkumara, Sanju Nagendrana,
Princeton Benedicta
a Department of Mechanical Engineering,
LoyolaICAM College of Engineering and Technology (LICET), Chennai, India
Abstract– This paper aims at analyzing a doublesided direct liquid cooling system required for a Silicon (Si) based power electronic converter used in electric/hybrid vehicles. The major parameters involved are flow rate of coolant, power dissipated by the die and the temperatures at various points. These parameters are analyzed for various flow rates, various power sources and the simulated results are verified with the superposition theorem. Comparing the thermal resistance for various combinations of flow rates, an optimum cooling mode is proposed.
Keywords Power electronic converter; Cooling system; Electric/hybrid vehicle; Superposition theorem; Direct liquid cooling; Thermal resistance; Doublesided cooling.

INTRODUCTION
Power Electronic converters are required for converting Alternating Current (AC) to Direct Current (DC) and vice versa in an electric/hybrid vehicle. This means that the power module receives AC from the power socket which should be converted to DC in order to supply it to batteries for charging. In a similar way, the DC power from batteries should be converted to AC by the converter as the traction motors require an AC supply. During the conversion, conduction and commutation losses in the dies are always present, because of which the conversion efficiency of the converter is never 100% [1]. Consequently, this generates heat in the Power Electronic converter. The flux repartition depends on the temperature gradient between the die and the coolant temperature. This temperature gradient depends on several factors such as mode and type of cooling method chosen, the fluid used for cooling, flow rate of the coolant and power dissipated by the module.
Three types of cooling methods namely single high sided cooling, single low sided cooling and doublesided cooling have been reported to be used [1]. In direct liquid cooling method, the coolant is in direct contact with the dissipating surface without an intermediate contact layer. Direct liquid cooling provides promising results as they efficiently carry away the generated heat and so it has been chosen for the study [2]. When forced convection liquid cooling methods are implemented, the power density of the semiconductor package can be greatly improved [3]. This paper aims to determine the thermal resistance during the power conversion process.
James Deepak b
b Assistant Professor, Department of Mechanical Engineering,
LoyolaICAM College of Engineering and Technology (LICET), Chennai, India
A cooling system with embedded Si module having four heat dissipating dies in it is chosen for the analysis. The liquid coolant used to carry away this heat is water. The flow rate of coolant is varied for upper and lower flow of doublesided cooling and analysis is done for steady state. The simulation is done for different boundary conditions and the flux repartition results are obtained. In this way the temperatures of the various dies of the power module are calculated for variation in parameters of cooling system and power sources. The outcomes of various conditions are compared, and the optimum flow rate and cooling method is proposed.
In a direct liquid cooling system, no thermal interface is used. By direct liquid cooling method, nearly 30% reduction in thermal resistance and 40% reduction in size has been achieved in InsulatedGate Bipolar Transistors (IGBTs) [4]. This proves the efficiency of direct liquid cooling systems.

MODEL DESCRIPTION AND SPECIFICATIONS While the major focus in our study is the flux repartition
within the power module, a simplified version of the power module was considered. An actual model of the typical power electronic converter could have as much as 31 million elements, which would require a lot of computation time [1]. Whereas, the simplified model of the module amounted to around 4.6 million elements which significantly reduced the computing time.
This arrangement consists of a case made of plastic which serves as a passage for coolant flow above and/or below the module. The entire packaging is contained within the case. The package consists of copper spreaders, dielectric gel, silver sintering, dies and bumps. The coolant flows through the passage in this arrangement such that the coolant has direct contact with the spreaders on both the sides.
Fig. 1. Entire package crosssection
The die is sintered using silver to the copper spreader on one side. Sintering is preferred in this case as it can offer better electrical and thermal performance.[5] On the other side, a bump is used as an interconnection between the die and the spreader. The bumps serve the purpose of dissipating the heat in a better way, especially for double sided cooling method.[6]
TABLE I. MATERIAL SPECIFICATIONS
Fig. 2. Die Assembly
The silver sintering makes direct contact with the lower spreader for IGBT1 and Diode 1. For IGBT2 and Diode 2, the silver sintering makes direct contact with the upper spreader. Depending on the position on the sintering, the bump takes its position on the opposite side. The die, silver sintering and equivalentbump assemblies are surrounded by dielectric gel to prevent contamination of the electronic components.
Fig. 3 Power sources
The duct where the coolant flows is plain for its entire length. This makes sure that turbulence is negligible. Hence, a laminar flow has been assumed for the analysis.
All dimensions are in mm
Fig. 4 Orthographic and isometric views
Material
Density
[kg/m3]Conductivity
[W/m.K]Specific heat [J/kg.K]
Copper
8900
390
390
Silver sintering
10500
100
232
Dielectric gel
1600
0.7
3000
Die (Si)
3000
125
500
Plastic (ABS)
1000
0.25
1005
Equivalent bump
3888
81
1104
Material
Density
[kg/m3]Conductivity
[W/m.K]Specific heat [J/kg.K]
Copper
8900
390
390
Silver sintering
10500
100
232
Dielectric gel
1600
0.7
3000
Die (Si)
3000
125
500
Plastic (ABS)
1000
0.25
1005
Equivalent bump
3888
81
1104

COMPUTATIONAL MODEL

Initial Setup
All the simulations were performed using ANSYS Fluent* R2 2019. The void conduits were filled so as to account for the space in which the coolant flows. The topology of the entire model was shared. The threedimensional model of the whole package was imported as a STEP file under the geometry of the Fluent* Analysis system in workbench*. Firstly, a Multizone meshing method is introduced for all the bodies.

Meshing
The mapped mesh type is set to be Hexahedral and the free mesh type is set to be Hexahedral Core. The mesh sizes are given in the table below.
TABLE II. MESH SIZE OF COMPONENTS
Body
Mesh Size (mm)
Sintering
5×102
Dies
5×102
Case
5
Fluid Volume
0.5
Equivalent Bump
5×102
Spreader and Dielectric Gel
0.5
As the sintering and die has a lower thickness when compared to the other bodies, a finer mesh has been used for these bodies. The fluid volumes and the region surrounding the die assemblies are also critical in the calculation. So, they have also been assigned with smaller mesh sizes.
* ANSYS Fluent & Workbench are registered trademarks of US based ANSYS Inc. and its subsidiaries.
Fig. 5. Isometric view of fluid volumes, spreaders and die assembly
Fig. 6. Top View of the dielectric gel and die assembly
Fig. 5 and Fig. 6 Clearly indicate the differences in the mesh sizes, which is recognized by the darker and lighter portions.
Fig. 7. Front view of the whole package
Fig. 8. Front View of the fluid volumes, spreader and die assembly
The model was then imported into Fluent. A steady state, pressurebased solver with absolute velocity formulation was used.

Cell zone and Boundary conditions
The flow was considered as Laminar and energy equation has been implemented. The various properties considered were density, specific heat and conductivity. Each component was assigned with the corresponding material property. The power dissipated was represented as a heat flux (100 W) to the corresponding die depending on the case under consideration. The inlets and outlets were mentioned as mass flow type. Firstly, the coolant flow rate was assumed to be 4 lit/min on both the sides of the module. The inlet temperature of the coolant is taken to be 338.15 K at steady state and the reference temperature is taken to be 293 K, which is the ambient temperature. This model was subjected to five different test cases. In the first case, a power source of 100 W was given to IGBT1 present at the immediate left of the module when viewed from the fluid inlet face. Similarly, for case 2, case 3, case 4 and load case, corresponding power sources were assigned. Finally, an analysis was made by varying the flow rate of coolant on both sides of the inlets and thermal resistances were calculated. The thermal resistances for various flow rates [(upper, lower)] of the coolant, (0l/min,8l/min), (2l/min,6l/min), (4l/min,4l/min), (6l/min,2l/min) and (8l/min,0l/min) were analyzed for the power dissipation cases 1 and 4.
The power input to the dies are varied for different cases of analysis as given below.
TABLE III. POWER DISTRIBUTION FOR VARIOUS CASES
IGBT 1
IGBT 2
Diode 1
Diode 2
Case 1
100
0
0
0
Case 2
0
100
0
0
Case 3
0
0
100
0
Case 4
0
0
0
100
Load case
300
300
100
100
All power values are in Watt

PostProcessing
Temperatures of the four dies were extracted from the software using volumeintegral technique of calculation. The outlet temperature of the fluid both in the upper and lower passage was found using mass average method. Reference planes were created exactly at the midsection of each die. The crosssectional views of the temperature profile of the dies were captured using the contours created using these planes.
An isometric view displaying the heat zones for the load case, at various positions of the entire module along with the plastic casing is shown.


RESULTS

Results of simulations
Simulation was performed for steady state condition for a flow rate of 4 lit/min on both the ducts. 200 iterations were set initially, and the solution converged at 163rd iteration.
Fig. 9. Residual graph
A temperature profile was obtained at the cross section of the mid portion of the die which dissipates power. This is represented using contours of temperature. The following results indicate the temperature profiles for various cases.
Fig. 11. Overall Temperature Contour
The mean temperatures of the various dies were calculated for each case using volume integral method. The following table gives the temperatures of the various dies for various cases.
TABLE IV. TEMPERATURE ACROSS DIES
IGBT 1
IGBT 2
Diode 1
Diode 2
Case 1
372.825
347.501
345.278
342.953
Case 2
344.815
373.18
341.51
345.54
Case 3
345.278
342.954
372.832
347.507
Case 4
341.51
345.541
344.819
373.187
Load case
473.759
484.751
412.057
420.316
(a) Case 1
(b) Case 2
(c) Case 3
(d) Case 4
(e) Load case mid sec 1
(f) Load case mid sec 2
(a) Case 1
(b) Case 2
(c) Case 3
(d) Case 4
(e) Load case mid sec 1
(f) Load case mid sec 2
All temperatures are in Kelvin.

Superposition theorem
The difference between the temperature at IGBT1 and inlet temperature of water for case 1 is represented by Tj (1) and
1
1
for case 2 it is represented by Tj (2)
1
respectively. Similarly,
(1) (2)
(1) (2)
for IGBT2, it is Tj2 , Tj2 and so on. P is the power given to the dies during the load case while P1, P2, P3, P4 are the powers given to IGBT1, IGBT2, diode 1 and diode 2 respectively during cases 1, 2, 3 and 4. The general formula for superposition theorem is given as
(1) (1) (2) (2)
(1) (1) (2) (2)
Tj1 =(P/P1 )xTj1 +(P/P2 )xTj1 + …
1
1
+(P/Pn(n))xTj (n)
Tj2 =(P/P (1))xTj (1) +(P/P (2)) xTj (2) + …
1 2 2 2
(n)
(n)
(1) (1) (2) (2)
(1) (1) (2) (2)
+(P/Pn(n))xTj2
Tjn =(P/P1 )xTjn +(P/P2 )xTjn + …
(n)
(n)
+(P/Pn(n))xTjn
In the load case, the power given to the IGBTs is 300W and for diodes it is 100W, the above formula is modified for the IGBT1 as
Fig. 10 MidSection temperature contour
Tj1(IGBT1 load cae by superposition theorem)= (300/100)xTj(1,IGBT1)+(300/100)xTj(2,IGBT1)+ (100/100)xTj(3,IGBT1) + (100/100)xTj(4,IGBT1) (1)
Temperature difference for IGBT1 for superposition theorem was found using eqn. (1). Similar modifications were made to the formula for the other IGBT and the two diodes and calculations were made to find Tj1, Tj2, Tj3, Tj4. Applying the inlet temperature of the coolant (338.15 K) and temperatures of the dies obtained for the four cases in the superposition theorem, results were compared with the temperatures of the load case.
TABLE V. VERIFICATION OF SUPERPOSITION THEOREM
Formulae results
(Temperature difference)
Load case results
(Temperature difference)
Tj1
134.508
135.609
Tj2
145.338
146.601
Tj3
72.815
73.907
Tj4
80.973
82.166
Although the entire module is symmetrical on all sides, a significant difference between Tj1 and Tj2 and between Tj3 and Tj4 (Which represents the temperature difference between the inlet and IGBT1, IGBT2, Diode 1 and Diode 2 respectively) can be observed. This could be attributed to the position of the IGBTs and the diodes with respect to the distance from the coolant inlets. As, IGBT1 is positioned nearer to the inlets than IGBT2, the coolant passes over IGBT1 first. This leads to a slight increase in the coolant temperature before it reaches the spreader area over IGBT2. Therefore, IGBT1 is at a lower temperature compared to IGBT2.

Determination of optimum flow rates
For equal cumulative flow rates of coolant, by varying the flow rates of the upper and lower inlets, thermal resistances were obtained and plotted for two power dissipation cases (Case 1 & Case 4)
Fig.12. Flow Rate vs Thermal Resistance for Case 1
Fig.13. Flow Rate vs Thermal Resistance for Case 4
In both the cases, the thermal resistance has the least value when equal flow rate of the coolant is used on the high side and low side of the dies. In case 1, there is a difference between the thermal resistance values when the flow rate is (2l/min,6l/min) and (6l/min,2l/min) respectively. In this case the die is sintered to the spreader close to the low side. Thermal resistance has a lower value when the flow rate is (2l/min,6l/min). This is because, more heat is dissipated on the lower side than on the higher side for Case 1. So, increasing the flow rate on this side causes better heat transfer. The same phenomenon applies to all the other cases.


CONCLUSIONS
Temperatures calculated using superposition theorem and the temperatures obtained from various diodes of the load case were compared. The results had a maximum deviation of 1.6%. Hence, the superposition theorem was verified. Furthermore, for the same coolant flow rate, the thermal resistance is found to be minimum when there is an equal coolant flow rate on both of the ducts. The thermal resistance decreases by around 17.3% in double sided cooling when compared with single sided cooling. This proves the efficiency of doublesided cooling technique. This concept can be extended to find the transient response of the impedances.
REFERENCES

Sanchez, SÃ©bastien & Nguyen, C. & Cadile, Claudia & Fradin, Jean Pierre & Tounsi, Patrick & Reynes, JeanMichel. (2020). Electrothermal Characterization of DoubleSided Cooling Si Power Module. 10.1007/9783030371616_5.

M. SchneiderRamelow, T. Baumann and E. Hoene, "Design and assembly of power semiconductors with doublesided water cooling," 5th International Conference on Integrated Power Electronics Systems, Nuremberg, Germany, 2008, pp. 17.

B. C. Charboneau et al., "DoubleSided Liquid Cooling for Power Semiconductor Devices Using Embedded Power Packaging," in IEEE Transactions on Industry Applications, vol. 44, no. 5, pp. 16451655, Sept.Oct. 2008.

Hitachi, T. & Hiromichi, G. & Nagaune, F.. (2012). Direct liquid cooling IGBT module for automotive applications. Fuji Electric Rev. 58. 5559.

C. Buttay, B. Allard and R. Riva, "Silver sintering for power electronics integration," 2015 International Conference on Electronics Packaging and iMAPS All Asia Conference (ICEP IAAC), Kyoto, 2015, pp. 554558, doi: 10.1109/ICEP IAAC.2015.7111076.

ChiaChi Tsai, LiLing Liao, YenFu Su, T. Hung and K. Chiang, "Reliability analysis of copper bump interconnection in double sided power module," 2015 16th International Conference on Thermal, Mechanical and MultiPhysics Simulation and Experiments in Microelectronics and Microsystems, Budapest, 2015, pp. 14, doi: 10.1109/EuroSimE.2015.7103093.