International Journal of Power Electronics and Drive Systems (IJPEDS)
Vol. 13, No. 4, December 2022, pp. 2313~2323
ISSN: 2088-8694, DOI: 10.11591/ijpeds.v13.i4.pp2313-2323  2313
Journal homepage: http://ijpeds.iaescore.com
Experimental and simulation approach of cooling system in
3-phase inverter using extended surface
Agus Mukhlisin1,2
, Prisma Megantoro1
, Estiko Rijanto3
, I Nyoman Sutantra2
, Lilik Jamilatul Awalin1
,
Yoga Uta Nugraha1
, Indra Sidharta2
1
Electrical Engineering Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia
2
Mechanical Engineering Department, Institut Teknologi Sepuluh Nopember, Surabaya, Indonesia
3
Research Center for Smart Mechatronics, National Research and Innovation Agency (BRIN), Jakarta Pusat, Indonesia
Article Info ABSTRACT
Article history:
Received Apr 19, 2022
Revised Aug 1, 2022
Accepted Aug 17, 2022
Overheating is a failure mode that significantly affects the reliability of
electronic devices. All electronic devices, including a 3-phase inverter
driving a traction motor, produce heat dissipation. Heat dissipation needs to
be controlled with cooling to prevent overheating. Overheating can be
avoided by increasing cooling or reducing heat dissipation. Heat dissipation
in the 3-phase inverter is caused by the internal resistance of the metal–
oxide–semiconductor field-effect transistor (MOSFET), switching loss, and
other factors. Cooling for the 3-phase inverter can use water coolant or air
coolant. The cooling system is based on the amount of heat dissipation
produced. Cooling of a 3-phase inverter can use air coolant with the addition
of an extended surface area in the heat sink. The heat sink uses aluminum
material, often called pin fin. There are kinds of aluminum available in the
market. We calculated heat generation based on the MOSFET's internal
resistance, switching loss, and other factors. We validated the simulation
results experimentally using a thermal camera. Thus, we could find an
optimal number, dimensions, and aluminum type of fin for the cooling
system in the 3-phase inverter.
Keywords:
3-phase inverter
Cooling system
Extended surface
Heat dissipation
Heat sinks
This is an open access article under the CC BY-SA license.
Corresponding Author:
Agus Mukhlisin
Department of Electrical Engineering, Faculty of Advanced Technology and Multidiscipline
Airlangga University
Gedung Kuliah Bersama, Campus C, Mulyorejo, Surabaya 60115, Indonesia
Email: agus.mukhlisin@ftmm.unair.ac.id
1. INTRODUCTION
Electronic component reliability is influenced by a cooling system [1]. Each electronic component
will produce heat dissipation, including the 3-phase inverter for a motor controller. Heat dissipation needs to
be controlled with a cooling system to prevent overheating. It can also be done by reducing the amount of
heat dissipation produced, such as using the right switching strategy [2]–[4], using field-oriented control
(FOC) [5], [6], using intelligent control [7]–[11], and the others. The most dominant heat source in the motor
controller is losses in the 3-phase inverter. Part of the 3-phase inverter that will produce heat is the metal–
oxide–semiconductor field-effect transistor (MOSFET).
One of the MOSFET failure modes is overheating. The heat generated by MOSFETs must be
transferred quickly to the cooling media with intelligent cooling using a fan [12], thermoelectric cooling [13],
others so that the temperature of the MOSFET does not rise above its ability limit. Heat transfer can occur in
several ways, including the following: conduction, convection, and radiation [14]–[16]. This paper will
 ISSN: 2088-8694
Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323
2314
discuss the MOSFET cooling method using an extended surface (pin fin). The extended surface area of heat
transfer improves the thermal efficiency [17], [18].
Computational fluid dynamic (CFD) simulation can be used to obtain a cooling method that suits your
needs [19], [20]. This method can help to analyze the cooling method without spending time doing
experiments [21]. The numerical method can show the temperature distribution and how much heat must be
cooled [22]. If the final design has been obtained, this design can be used as a reference in making a prototype.
This work presents optimized heat transfer from MOSFET in a 3-phase inverter by using extended surfaces with
some pin fins and dimensions. The calculation of heat loss in the 3-phase inverter is used as the input in the
simulation. The simulation result is validated with the experimental results. The validation is carried out based
on the temperature distribution and the temperature's value at the same point. If the error tolerance that occurs
does not exceed 5%, we consider that an accurate simulation is obtained. Therefore, it can be used to optimize
the number of fins and the dimensions of the fins. Our target is that the 3-phase inverter can operate normally at
a power of 30 kW.
2. METHOD
First, we performed a heat loss calculation in 3-phase inverter for parameter input in the simulation.
Next, the simulation results were validated with the experimental results, where the temperature distribution
was obtained using a thermal camera. After the simulation results are validated, it is possible to optimize the
heat sink according to the heat loss. As a result, the optimal heat sink configuration is obtained with the
controller not overheating.
Heat losses can be seen in the power flow diagram in Figure 1. Based on Figure 1, the 3-phase
inverter losses consist of conduction loss pulse-code modulation PCM, switching loss part submission
warrant PSW, dead time loss PD, and gate charge loss PG, as given in (1) to (4) [23], [24].
𝑃𝐶𝑀 =
1
𝑇𝑆𝑊
∫ 𝑝𝐶𝑀(𝑡)𝑑𝑡
𝑇𝑆𝑊
0
=
1
𝑇𝑆𝑊
∫ 𝑅𝐷𝑆𝑜𝑛. 𝑖𝐷
2(𝑡)𝑑𝑡 = 𝑅𝐷𝑆𝑜𝑛. 𝐼𝐷𝑟𝑚𝑠
2
𝑇𝑆𝑊
0
(1)
𝑃𝑆𝑊 =
1
2
𝑉𝐼𝑁. 𝐼𝑂. (𝑡𝑟 + 𝑡𝑓). 𝑓𝑆𝑊 (2)
𝑃𝐷 = 𝑉𝐷. 𝐼𝑂. (𝑡𝐷𝑟 + 𝑡𝐷𝑓). 𝑓𝑠𝑤 (3)
𝑃𝐺 = (𝑄𝑔−𝐻 + 𝑄𝑔−𝐿). 𝑉
𝑔𝑠. 𝑓𝑠𝑤 (4)
The symbol 𝑅𝐷𝑆𝑜𝑛 is the MOSFET drain-source internal resistance (Ω), TSW is switching time (s), iD is the
current flowing in drain (A), and IDrms is the current flowing from source to drain (A). Vin and Io represent the
battery voltage (V) and the output current (A). tr is the rise time of the MOSFET (s), tf is the fall time of the
MOSFET (s), and fSW is the switching frequency (Hz). VD, tDr, and tDf denote the low-side MOSFET – diode
forward-voltage (Volt), dead-time at rising (s), and dead-time at falling (s). Vgs is the gate drive voltage of the
MOSFET (V), Qg – H is the high side MOSFET gate electric charge (C), and Qg – L is the low side MOSFET
gate electric charge (C). This controller uses the type of MOSFET IRFB4115 with the specifications shown
in Figure 2.
Figure 1. Power flow 3-phase inverter Figure 2. MOSFET IRFB4115 specification [25]
Int J Pow Elec & Dri Syst ISSN: 2088-8694 
Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin)
2315
This controller has the following target specifications; i) maximum voltage operation 150 VDC, ii)
maximum power 30 KW, and iii) using an air-cooling system. We first made a simplified 3D model of the 3-
phase inverter for numerical analysis using CFD software, as shown in Figure 3. The heat source on the 3-
phase inverter is the MOSFET component. So that heat transfer can be analyzed with thermal resistance as
shown in Figure 4.
Figure 3. Simplification of 3D model of the 3-phase inverter
Figure 4. Thermal resistance model controller
The heat source from the MOSFET moves either up to the controller cover or down to the heat sink.
The heat transfer downward will be more dominant than upward. The controller cover traps the air contained
inside the controller container, so it does not allow air circulation. The heat transmitted downward propagates
through the thermal pad.
A thermal pad is a component that serves as an insulator between the base metal of the MOSFET
and the heat sink. The use of the thermal pad must have good thermal conductivity as well as a good
insulator. The motor controller failure can also be caused by improper thermal pad selection. If the heat
generated by the MOSFET is not transferred to the heat sink quickly, the MOSFET will experience
overheating, resulting in failure. Similarly, if the thermal pad cannot isolate the metal base of the MOSFET, a
short circuit will occur, and the controller will fail. Generally, the thermal pad uses silicone material. The
heat transfer then continues to the heat sink fin by conduction. The equation of heat transfer by conduction is
given in [14], [16]:
𝑞𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑜𝑛 = 𝑘. 𝐴. ∆𝑇 (5)
The type of heatsink material affects the quality of heat transfer: the more significant the heat sink material's
thermal conductivity (k), the more remarkable its heat propagation ability. Our heat sink uses 1060 series
 ISSN: 2088-8694
Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323
2316
aluminum material, often found in the local market. Convection sends the heat transfer into the air around the
heat sink fin. The equation of heat transfer by convection is given by (6) to (8) [14], [16]:
𝑞𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛 = ℎ. 𝐴. ∆𝑇 (6)
𝑁𝑢 =
ℎ.𝐿
𝑘𝑓
= 𝑓(𝑅𝑒𝐿, Pr) (7)
𝑅𝑒𝐿 =
𝜌.𝑉.𝐿
𝜇
(8)
The symbols h and A represent the convection coefficient (W/m2
K) and the heat transfer area (m2
). L
is the length of material (meters), kf is the thermal conductivity of the fluid (W/mK), and Nu is the Nusselt
number. Re and Pr denote Reynold’s number and Prandtl’s number. ρ is the density of air (kg/m3
), V is the
velocity of air (m/s), and μ is the friction coefficient.
Convection heat transfer can be optimized by increasing the convection coefficient and cooling area.
Increasing the convection coefficient can be done by increasing the speed of air flowing at the fin heat sink,
for example, by adding a fan, directing the air with a guide fan so that air can be directed precisely at the
heatsink. The cooling area can be enlarged by using variations in dimension and the number of fins on the
heatsink. The quality of the simulation results is determined by the dimensions and meshing used. Our
meshing uses a curvature-based parameter mesh with total nodes of 615,925, using Jacobian points by 4
points. Table 1 lists material types used in the motor controller cooling simulation. Figure 5 shows the
meshing model.
Table 1. Material types for cooling simulation
Component Materials Thermal conductivity (W/mK)
Fin heat sink Aluminum 1060 200
Base plat Aluminum 1060 200
Thermal pad Silicone 3.2
MOSFET plastic body Plastic Silicone 124
MOSFET base metal Aluminum 5052 137
Figure 5. Meshing model
Connection contact set using bounded type. The contact that needs to be defined is between
MOSFET with the thermal pad and the thermal pad with the heat sink. This step ensures that the connections
between different materials are perfectly connected. We conducted numerical simulation using CFD with the
following conditions:
− Surrounding temperature: 28 o
C
− Controlled variable: Convection coefficient (h)
− The source of the heat loss: the MOSFET
First, we calculated heat loss when the motor controller is under a 1.5 kW load as show in Table 2.
Afterward, the simulation can be run, and temperature distribution can be obtained. Simulation results will be
validated using experimental results. Experiments with a maximum load of 30 kW continuously (30 minutes)
cannot be done because the dynamometer will overheat. So, the experiment was carried out with a continues
load of 1.5 kW. Furthermore, the results of these experiments can be used to validate the simulation results.
The validated simulation results can then be used to test the maximum loading, so that the temperature
Int J Pow Elec & Dri Syst ISSN: 2088-8694 
Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin)
2317
distribution is obtained when the maximum loading to obtain the temperature distribution under the
maximum loading. This simulation can also be used to optimize the fin on the heat sink so that cooling is
obtained as needed with minimum fin dimensions.
Figure 6 shows the experimental setup with (a) schema of experiment and (b) installed equipment in
this paper. It is equipped with an eddy current dynamometer paired with a 30 KW motor. The temperature
sensor uses a 10K thermistor recorded in a datalogger (USB to CAN Datalogger). An oscilloscope with four
channels is used to measure voltage and current, both input and output. A thermal camera is used to obtain
the temperature distribution.
Table 2. Calculated heat loss in the motor controller under 1.5 kW load
Type Heat Loss Parameter Value
Conduction Loss RDS(on)= 11 mΩ; ID = 8.15 A; Vin = 80.34 V;
𝑡𝑟=73 ns; 𝑡𝑓=39 ns;𝑓𝑆𝑊=20kHz;
td(on) = 18 ns; td(off) =41ns; Vgs = 23V; Qg =120 nC.
0.617 W
Switching Loss 0.733 W
Deadtime Loss 0.772 W
Gate Charge Loss 0.055 W
Total Heat Loss 2.177 W
(a) (b)
Figure 6. Experiment setup (a) schema and (b) equipment
3. RESULTS AND DISCUSSION
3.1. Simulation validation through experiment
The load was kept constant at 1.5 kW with a motor speed of 950 rpm and variations in wind speed.
We changed wind speed by controlling a fan power to obtain speed 1, speed 2, and speed 3. CAN to USB
datalogger was used to obtain data, among others: temperature of heat sink, phase current, battery current,
battery voltage, and the motor speed. One-time data collection took 35 minutes or until the heat sink
temperature got steady. The result of the experiment is shown in Figure 7.
Figure 7. Motor speed and heat sink temperature
 ISSN: 2088-8694
Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323
2318
The heat sink temperatures in Figure 7 were measured using a thermistor sensor. Yellow, red, blue,
and green line curves are the heat sink temperatures with variations in wind speed 1, speed 2, speed 3, and
natural convection. Each line's steady-state temperature is yellow at 33.37 o
C, red at 33.96 o
C, blue at
34.92 o
C, and green at 39.72 o
C. The greater wind speed yields a lower temperature of the heat sink. This is
in accordance with (6), (7), and (8). The trend line that occurs tends to be steady at a specific time. This
indicates that the cooling is sufficient so that overheating will not occur. The temperature distribution can be
captured in this steady state with a thermal camera.
Figures 8 and 9 show the temperature distribution results obtained from the simulation and
experiment. Figure 8 shows the temperature distribution using h=40 W/m2
K, both experimentally and in
simulation results. Figure 8(a) shows the simulation results where the thermistor sensor laying point shows a
temperature of 33.1 o
C. While Figure 8(b) is a photo taken by the thermal camera where the location of the
thermistor shows the value of 33.3 o
C. These results indicate that the simulation results have a 0.6% error.
The temperature distribution in general also has similarities. Under continuous loading of 1.5 kW, the highest
temperature occurs in the DC-DC converter, yielding 37 o
C.
(a) (b)
Figure 8. Temperature distribution with h=40 W/m2
K (a) simulation result and (b) experiment result
Figure 9 shows the temperature distribution using h=33 W/m2
K for both experimental and
simulation results. Figure 9 (a) shows the simulation results where the thermistor sensor laying point shows a
temperature of 34.3 o
C. In comparison, Figure 9 (b) is a photo taken by the thermal camera where the
location of the thermistor shows the value of 34.9 o
C. The result indicates that the simulation results have a
1.75% error. Based on the results shown in Figures 8 and 9, the simulation results can be accepted with an
error tolerance of ± 1.75%. Furthermore, the simulation method can be used to evaluate the maximum
loading of the 3-phase inverter, i.e., with a load of 30 kW. Table 3 summarizes the data obtained from the
simulation and experimental results.
(a) (b)
Figure 9. Temperature distribution with h=33 W/m2
K (a) simulation result and (b) experiment result
Int J Pow Elec & Dri Syst ISSN: 2088-8694 
Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin)
2319
Table 3. Comparison of simulation with experiment results
Coefficient of Convection Experiment Simulation Percentage of Error
h= 33 W/ m2
K 34.9 o
C 34.3 o
C 1.75 %
h= 40 W/ m2
K 33.1 o
C 33.3 o
C 0.6 %
3.2. Optimization of the heat sink under maximum load
We calculated heat loss produced by each MOSFET when the 3-phase inverter operates under the
maximum load 30 kW. Table 4 shows the parameter values and the heat loss values. Since we have 48
MOSFETs, the total heat loss is 1,048.27 Watts. We used this heat loss value as the input to our CFC model
and conducted a numerical simulation. Figure 10 shows the temperature distribution of the heat sink under a
maximum load of 30 kW. The maximum temperature is 98.9 o
C, so the controller does not overheat. The
temperature where the MOSFET is situated shows a value of 38.5 o
C. Since the temperature is significantly
lower than the critical temperature of the MOSFET, the heat sink fin is not required in this area.
Table 4. Heat loss of each MOSFET in the motor controller under 30 kW Load
Type Heat Loss Parameter Value
Conduction Loss RDS(on)= 11 mΩ; ID = 65.16 A; Vin = 80.34 V;
𝑡𝑟=73 ns; 𝑡𝑓=39 ns;𝑓𝑆𝑊=20kHz;
td(on) = 18 ns; td(off) =41ns; Vgs = 23V; Qg =120 nC.
14.491 W
Switching Loss 3.551 W
Deadtime Loss 3.742 W
Gate Charge Loss 0.055 W
Total Heat Loss 21.839 W
Figure 10. Simulation result under maximum load 30 kW
Furthermore, we optimized the number and length of the heat sink fins with maximum thermal load
and h=60 W/m2
K. The original length of ach heat sink fin is 20 mm. Figure 11 shows the simulation results.
Figure 11(a) shows the temperature distribution using 15 pieces of heat sink fins. The highest temperature of
the MOSFET is 103 o
C. The temperature is lower than the critical operation temperature of the MOSFET,
which is 125 o
C. Figure 11(b) shows the temperature distribution using 12 pieces of heat sink fins. The
highest temperature of the MOSFET is 107 o
C. The result indicates that reducing of fin number to 12 pieces
yielding maximum temperature that is lower than the critical temperature of the MOSFET. Further
optimization is carried out by reducing the fin length to 12 mm and yielding the temperature distribution
shown in Figure 11(c). The maximum temperature of the MOSFET in this condition is 117 o
C, which is
lower than the critical temperature of MOSFET. However, there is an increase in temperature compared to
the result of 12 pieces 20 mm fin. Therefore, the controller under a maximum load of 30 kW is optimally
cooled using the configuration of 12 pieces of fins with 12 mm length, following Figure 11(c).
 ISSN: 2088-8694
Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323
2320
(a)
(b)
(c)
Figure 11. Optimized heat sink in the number and length of the fins: (a) 15 pieces of 20 mm long,
(b) 12 pieces of 20 mm long, and (c) 12 pieces of 12 mm long
Int J Pow Elec & Dri Syst ISSN: 2088-8694 
Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin)
2321
4. CONCLUSION
Thermal analysis of the controller cooling system was carried out using modeling and simulation.
The simulation was validated experimentally with an error of ± 1.75%. Furthermore, the simulation method
can be used to optimize the cooling design for a maximum loading of 30 kW. The results indicate that the
controller does not experience overheating, as demonstrated by the MOSFET temperature. Further
optimization, through simulation was carried out in order to reduce the size and the number of the fin. The
simulation result of the optimized heat sink suggests that the maximum temperature of the MOSFET is lower
than its critical temperature. Hence, the controller did not experience overheating.
ACKNOWLEDGEMENTS
Author thanks to Faculty of Advanced Technology and Multidiscipline, Airlangga University and
Lembaga Pengelola Dana Pendidikan (LPDP) in most cases, sponsor and financial support
acknowledgments. This program support Sustainable Development Goals (SDGs) for affordable and clean
energy also goal 9 for Industry, Innovation and Infrastructure.
REFERENCES
[1] E. Laloya, Ó. Lucía, H. Sarnago and J. M. Burdío, “Heat management in power converters: from state of the art to future ultrahigh
efficiency systems,” IEEE Transactions on Power Electronics, vol. 31, no. 11, pp. 7896-7908, 2016, doi:
10.1109/TPEL.2015.2513433.
[2] N. A. Ramli, A. Jidin, Z. Rasin, and T. Sutikno, “Reduction of total harmonic distortion of three-phase inverter using alternate
switching strategy,” International Journal of Power Electronics and Drive Systems, vol. 12, no. 3, pp. 1598-1608, 2021, doi:
10.11591/ijpeds.v12.i3.pp1598-1608.
[3] D. Christen and J. Biela, “Analytical switching loss modeling based on datasheet parameters for mosfets in a half-bridge,” IEEE
Transactions on Power Electronics, vol. 34, no. 4, pp. 3700-3710, 2019, doi: 10.1109/TPEL.2018.2851068.
[4] A. V. Sant and K. R. Rajagopal, “PM synchronous motor speed control using hybrid fuzzy-pi with novel switching functions,”
IEEE Transactions on Magnetics, vol. 45, no. 10, pp. 4672-4675, 2009, doi: 10.1109/TMAG.2009.2022191.
[5] V. T. Ha, N. T. Lam, V. T. Ha, and V. Q. Vinh, “Advanced control structures for induction motors with ideal current loop response
using field oriented control,” International Journal of Power Electronics and Drive System, vol. 10, no. 4, pp. 1758-1771, 2019,
doi: 10.11591/ijpeds.v10.i4.pp1758-1771.
[6] M. Hasoun, A. el Afia, M. Khafallah, and K. Benkirane, “Field oriented control based on a 24-sector vector space decomposition
for dual three-phase pmsm applied on electric ship propulsion,” International Journal of Power Electronics and Drive Systems,
vol. 11, no. 3, pp. 1175-1187, 2020, doi: 10.11591/ijpeds.v11.i3.pp1175-1187.
[7] D. Brahim, B. Mokhtar, and L. Salah, “Improvement of adaptive fuzzy control to adjust speed for a doubly fed induction motor
drive (DFIM),” International Journal of Power Electronics and Drive Systems, vol. 11, no. 1, pp. 496-504, 2020, doi:
10.11591/ijpeds.v11.i1.pp496-504.
[8] E. Sabouni, B. Merah, and I. K. Bousserhane, “Adaptive backstepping controller design based on neural network for pmsm speed
control,” International Journal of Power Electronics and Drive Systems, vol. 12, no. 3, pp. 1940-1952, 2021, doi:
10.11591/ijpeds.v12.i3.pp1940-1952.
[9] A. Nait Seghir, T. Henni and M. Azira, “Fuzzy and adaptive fuzzy PI controller-based Vector control for permanent magnet
synchronous motor,” 10th Ieee International Conference On Networking, Sensing And Control (ICNSC), 2013, pp. 491-496, doi:
10.1109/ICNSC.2013.6548788.
[10] Jae Sung Yu, Sun Mo Hwang and Chung-Yuen Won, “Performances of fuzzy-logic-based vector control for permanent magnet
synchronous motor used in elevator drive system,” 30th Annual Conference of IEEE Industrial Electronics Society, IECON, 2004,
vol. 3, pp. 2679-2683, doi: 10.1109/IECON.2004.1432229.
[11] N. J. Patil, R. H. Chile and L. Waghmare, “Fuzzy adaptive controllers for speed control of PMSM drive,” International Journal of
Computer Applications, vol. 1, no. 11, pp. 975-8887, 2010, doi: 10.5120/233-387.
[12] A. Deep, J. Singh, Y. Narayan, S. Chatterji and L. Mathew, “Intelligent cooling system for three level inverter,” Communication,
Control and Intelligent Systems (CCIS), pp. 286-289, 2015, doi: 10.1109/CCIntelS.2015.7437925.
[13] U. Sanver, E. Yavuz and C. Eyupoglu, “An electronic control unit for thermoelectric cooling,” IEEE Conference of Russian Young
Researchers in Electrical and Electronic Engineering (EIConRus), 2019, pp. 141-145, doi: 10.1109/EIConRus.2019.8656871.
[14] F. P. Incropera, D. Witt, T. L. Bergman, and A. S. Lavine, Fundamentals of heat and mass transfer, Sixth Edition. New Jersey:
John Wiley & Sons, 2006.
[15] A. Blinov, D. Vinnikov, and T. Lehtla, “Cooling methods for high-power electronic systems,” Scientific Journal of Riga Technical
University, Power and Electrical Engineering, vol. 29, no. 1, pp. 79-86, 2011, doi: 10.2478/v10144-011-0014-x.
[16] Y. A. Çengel, Heat Transfer: A Practical Approach, New York, NY, USA: McGraw-Hill, 2003.
[17] S. S. Kang, “Advanced cooling for power electronics,” 7th International Conference on Integrated Power Electronics Systems
(CIPS), 2012, pp. 1-8.
[18] A. Driss, S. Maalej and M. C. Zaghdoudi, “Thermal modeling of the cooling of a power MOSFET by heat pipes,” International
Conference on Engineering & MIS (ICEMIS), 2017, pp. 1-6, doi: 10.1109/ICEMIS.2017.8273067.
[19] U. SanAndres, G. Almandoz, J. Poza and G. Ugalde, “Design of cooling systems using computational fluid dynamics and
analytical thermal models,” IEEE Transactions on Industrial Electronics, vol. 61, no. 8, pp. 4383-4391, Aug. 2014, doi:
10.1109/TIE.2013.2286081.
[20] Y. Bulut and K. Pandya, “Thermal modeling for power MOSFETs in DC/DC applications,” 5th International Conference on
Thermal and Mechanical Simulation and Experiments in Microelectronics and Microsystems, 2004, pp. 429-433, doi:
10.1109/ESIME.2004.1304074.
[21] O. Spuiber, M. M. De Souza, E. M. S. Narayanan and S. Krishnan, “Analyses of a COOL-MOSFET,” CAS '99 Proceedings. 1999
International Semiconductor Conference (Cat. No.99TH8389), 1999, pp. 131-134 vol. 1, doi: 10.1109/SMICND.1999.810446.
 ISSN: 2088-8694
Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323
2322
[22] M. Langer, Z. Lisik and J. Podgorski, “Numerical simulation of cool MOS transistor,” Experience of Designing and Applications
of CAD Systems in Microelectronics. Proceedings of the VI-th International Conference, CADSM, 2001, pp. 303-304, doi:
10.1109/CADSM.2001.975849.
[23] G. Dusan, M. Purschel, and A. Kiep, “MOSFET power losses calculation using the data-sheet parameters,” Infineon application
note, pp. 1-23, 2006.
[24] ROHM Semiconductor, “Application Note: Calculation of Power Loss (Synchronous),” 2016. [Online]. Available:
http://rohmfs.rohm.com/en/products/databook/applinote/ic/power/switching_regulator/power_loss_appli-e.pdf.
[25] International IOR Rectifier, “IRFB4115PbF HEXFET Power MOSFET Datasheet,” 2014. [Online]. Available:
https://www.infineon.com/dgdl/irfb4115pbf.pdf?fileId=5546d462533600a401535615ba6a1e0f.
BIOGRAPHIES OF AUTHORS
Agus Mukhlisin is a lecturer in Electrical Engineering Department at Faculty of
Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia. He
received a Bachelor’s degree and Master’s degree from Institute Technology of Sepuluh
November, Surabaya, Indonesia in 2015 and 2017. He is currently pursuing Doctoral Program
in Institute Technology of Sepuluh November, Surabaya, Indonesia. His Dissertation Topic
was about Intelligent Vector Control in Motor Controller. His research interests include
System and Control of Electric Vehicle, Battery Management System, Charging Station, and
Dynamic System Simulation. He can be contacted at email: agus.mukhlisin@ftmm.unair.ac.id.
Prisma Megantoro is a lecturer in Electrical & Electronic Engineering
Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University,
Surabaya, Indonesia. He received the bachelor degree and master degree from Universitas
Gadjah Mada, Yogyakarta, Indonesia in 2014 and 2018. His current research is focused on
solar photovoltaic technology, embedded system, and internet of things. He can be contacted
at email: prisma.megantoro@ftmm.unair.ac.id.
Estiko Rijanto has worked for Research Center for Smart Mechatronics - the
National Research and Innovation Agency (BRIN) - Indonesia since 2021. From 2002 to 2021,
he worked at Research Center for Electrical Power and Mechatronics, Indonesian Institute of
Sciences (LIPI). In 1987 he enrolled at Institut Teknologi Bandung (ITB). He completed his
B.Eng. degree at Tokyo University of Agriculture and Technology (TUAT), Tokyo - Japan, in
1993. He received M. Eng. and Dr. Eng. degrees from the same university in 1995 and 1998.
In 2000, he completed his book on “Robust Control: Theory for Application” during his post-
doctoral program at TUAT-VBL. In 2013, he was inaugurated as the Research Professor on
Applied Control Systems by LIPI after presenting his book “Integration of control and
information systems for industrial competitiveness,”. His research interests include control
systems and their applications for mechatronics, electric vehicles, power generations,
renewable energies, and battery management systems. He can be contacted at email:
esti003@brin.go.id.
I Nyoman Sutantra is a lecturer in Mechanical Engineering Department, Institut
Teknologi Sepuluh November, Surabaya, Indonesia. He received a Bachelor’s degree from
Institute Technology of Sepuluh November, Surabaya, Indonesia. He received a Master and
Doctoral degree from The University of Wisconsin – Madison in mechanical department. He
was inaugurated the Professor on mechanical engineering. His research interests include
automotive technology, electric vehicles, hybrid vehicle and steering control. He can be
contacted at email: tantra@me.its.ac.id.
Int J Pow Elec & Dri Syst ISSN: 2088-8694 
Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin)
2323
Lilik Jamilatul Awalin is a lecturer in Electrical & Electronic Engineering
Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University,
Surabaya, Indonesia. She received the bachelor degree from electrical engineering, Widya
Gama Malang, Indonesia in 1977. She received the master degree from Institute Technology
of Sepuluh November, Surabaya, Indonesia in 2004. Finally, she gets the doctoral degree from
University of Malaya, Malaysia in 2014. Her current research is focused on fault localization,
smart grid technology, transmission and distribution network, load shedding, coordination
protection and optimization. She can be contacted at email: lilik.j.a@ftmm.unair.ac.id.
Yoga Uta Nugraha is as a lecturer at Electrical Engineering program in the
Faculty of Advanced Technology and Multidicipline, Universitas Airlangga, Surabaya,
Indonesia. He was born in Sukoharjo, Central of Java. He is a Doctoral candidate in Electrical
Engineering Department of Sepuluh Nopember Institute of Technology Surabaya. He has a
Bachelor (2015) and Master (2020) of Electrical Engineering from Sepuluh Nopember
Institute of Technology Surabaya. His current research is focused on integration on Electric
Vehicle, electric machine, permanent magnet motor, new and renewable energy. He can be
contacted at email: yoga.uta.n@ftmm.unair.ac.id.
Indra Sidharta is a lecturer in Mechanical Engineering Department, Institut
Teknologi Sepuluh November, Surabaya, Indonesia. He received a Bachelor’s degree from
Institute Technology of Sepuluh November, Surabaya, Indonesia in 2004. He received a
Master degree from Hamburg University of Technology in field Materials Science, 2010. His
current research is focused on metal matrix composite, metallurgy, materials of battery. He
can be contacted at email: indra.sidharta@gmail.com.

More Related Content

PDF
Automated Air Cooled Three Level Inverter system using Arduino
PDF
THERMAL ANALYSIS OF HEAT SINK (VARIABLE SHIELD PROFILE) USED IN ELECTRONIC CO...
PDF
30120140503011
PDF
Micro processor based temperature controller on power transistors
PDF
Designing a Heat sink for motor controller unit
PDF
IRJET - Overheat Protection of Induction Motor
PDF
Thermal Characterization and Performance Evaluation of CPU Heat Sink Design
PDF
Design and Analysis of Heat Sink
Automated Air Cooled Three Level Inverter system using Arduino
THERMAL ANALYSIS OF HEAT SINK (VARIABLE SHIELD PROFILE) USED IN ELECTRONIC CO...
30120140503011
Micro processor based temperature controller on power transistors
Designing a Heat sink for motor controller unit
IRJET - Overheat Protection of Induction Motor
Thermal Characterization and Performance Evaluation of CPU Heat Sink Design
Design and Analysis of Heat Sink

Similar to Experimental and simulation approach of cooling system in 3-phase inverter using extended surface (20)

PDF
NUMERICAL ANALYSIS AND SIMULATION OF CONJUGATE HEAT TRANSFER STUDY OF ELECTRO...
PPTX
Simulation of Induction heater using IRFZ44N MOSFET
PDF
Electrically Powered Vehicles Battery Management System: Cooling for Parked V...
PDF
THERMAL ANALYSIS OF AIR FLOW IN A CPU CABINET WITH MOTHERBOARD AND HARD DISK ...
PDF
Three-level modified sine wave inverter equipped with online temperature moni...
PDF
Design, Optimization, and Analysis of Electric vehicle Battery Pack
PDF
IRJET- Thermal Analysis and Management for an Autonomous Underwater Vehicle
DOCX
568 IEEE TRANSACTIONS ON DEVICE AND MATERIALS RELIABILITY, VOL.docx
PDF
How to use CFD for Water Cooled Electronics
PDF
Ic2414251429
PDF
Heat pump design using peltier element For temperature control of the flow cell
PDF
Optimization Potential of Induction Heating Systems by Stefan Schubotz and Ha...
PDF
IRJET - A Review on Design and Optimization of Cooling Plate for Battery Modu...
PDF
Active cooling technique using peltier element for efficient heat mitigation ...
PDF
No 6 g3 solved problems
PDF
Enhancement of Thermo-mechanical Behavior of IGBT Modules through Engineered ...
PDF
THERMAL MODELING AND ANALYSIS OF 3- DIMENSINAL MEMORY INTEGRATION
PDF
How Computer Simulation Helps to Design Induction Heating Systems ASM 2001
PDF
How Computer Simulation Helps to Design Induction Heating Systems ASM 2001
PDF
Selection of Power Semiconductor Switches in M.H.B.R.I. Fitted Induction Heat...
NUMERICAL ANALYSIS AND SIMULATION OF CONJUGATE HEAT TRANSFER STUDY OF ELECTRO...
Simulation of Induction heater using IRFZ44N MOSFET
Electrically Powered Vehicles Battery Management System: Cooling for Parked V...
THERMAL ANALYSIS OF AIR FLOW IN A CPU CABINET WITH MOTHERBOARD AND HARD DISK ...
Three-level modified sine wave inverter equipped with online temperature moni...
Design, Optimization, and Analysis of Electric vehicle Battery Pack
IRJET- Thermal Analysis and Management for an Autonomous Underwater Vehicle
568 IEEE TRANSACTIONS ON DEVICE AND MATERIALS RELIABILITY, VOL.docx
How to use CFD for Water Cooled Electronics
Ic2414251429
Heat pump design using peltier element For temperature control of the flow cell
Optimization Potential of Induction Heating Systems by Stefan Schubotz and Ha...
IRJET - A Review on Design and Optimization of Cooling Plate for Battery Modu...
Active cooling technique using peltier element for efficient heat mitigation ...
No 6 g3 solved problems
Enhancement of Thermo-mechanical Behavior of IGBT Modules through Engineered ...
THERMAL MODELING AND ANALYSIS OF 3- DIMENSINAL MEMORY INTEGRATION
How Computer Simulation Helps to Design Induction Heating Systems ASM 2001
How Computer Simulation Helps to Design Induction Heating Systems ASM 2001
Selection of Power Semiconductor Switches in M.H.B.R.I. Fitted Induction Heat...
Ad

More from International Journal of Power Electronics and Drive Systems (IJPEDS) (20)

PDF
Comparative study by numerical simulation of two methods for automatic flow c...
PDF
Implementation of artificial intelligence for prediction performance of solar...
PDF
Dynamic fuel cell model improvement based on macroscopic energy representation
PDF
Investigation of common mode voltages of single stage boost inverter for five...
PDF
Dynamic modeling of an open cathode PEM fuel cell for automotive energy manag...
PDF
A modified CPS-PWM for capacitor voltage reduction of MMC based variable spee...
PDF
Firefly analytical hierarchy algorithm for optimal allocation and sizing of D...
PDF
Optimal power flow using archimedes optimizer algorithm
PDF
Economic evaluation of induction motor based on motor’s nameplate data and in...
PDF
Assessment of field oriented induction machine control strategy using new gen...
PDF
A current sensor fault diagnosis method based on phase angle shift technique ...
PDF
Optimized speed control with torque ripple reductions of BLDC motor based on ...
PDF
Design and implementation of DC linear actuator and stepper motor for remote ...
PDF
Milestone of the most used maximum power point tracking in solar harvesting s...
PDF
Intelligent torque observer combined with backstepping sliding-mode control f...
PDF
A novel multilevel inverter with reduced components and minimized voltage unb...
PDF
A synchronization technique for single-phase grid applications
PDF
Torque ripple alleviation of a five-phase permanent magnet synchronous motor ...
PDF
A novel interactive technique for load current harmonic reduction for any ran...
PDF
A 1 MHz soft-switching boost DC-DC converter with matching network
Comparative study by numerical simulation of two methods for automatic flow c...
Implementation of artificial intelligence for prediction performance of solar...
Dynamic fuel cell model improvement based on macroscopic energy representation
Investigation of common mode voltages of single stage boost inverter for five...
Dynamic modeling of an open cathode PEM fuel cell for automotive energy manag...
A modified CPS-PWM for capacitor voltage reduction of MMC based variable spee...
Firefly analytical hierarchy algorithm for optimal allocation and sizing of D...
Optimal power flow using archimedes optimizer algorithm
Economic evaluation of induction motor based on motor’s nameplate data and in...
Assessment of field oriented induction machine control strategy using new gen...
A current sensor fault diagnosis method based on phase angle shift technique ...
Optimized speed control with torque ripple reductions of BLDC motor based on ...
Design and implementation of DC linear actuator and stepper motor for remote ...
Milestone of the most used maximum power point tracking in solar harvesting s...
Intelligent torque observer combined with backstepping sliding-mode control f...
A novel multilevel inverter with reduced components and minimized voltage unb...
A synchronization technique for single-phase grid applications
Torque ripple alleviation of a five-phase permanent magnet synchronous motor ...
A novel interactive technique for load current harmonic reduction for any ran...
A 1 MHz soft-switching boost DC-DC converter with matching network
Ad

Recently uploaded (20)

PPTX
Management Information system : MIS-e-Business Systems.pptx
PPTX
AUTOMOTIVE ENGINE MANAGEMENT (MECHATRONICS).pptx
PDF
20250617 - IR - Global Guide for HR - 51 pages.pdf
PDF
Computer organization and architecuture Digital Notes....pdf
PPTX
Module 8- Technological and Communication Skills.pptx
PPTX
Information Storage and Retrieval Techniques Unit III
PPTX
ai_satellite_crop_management_20250815030350.pptx
PDF
Computer System Architecture 3rd Edition-M Morris Mano.pdf
PDF
August 2025 - Top 10 Read Articles in Network Security & Its Applications
PPTX
CONTRACTS IN CONSTRUCTION PROJECTS: TYPES
PPTX
PRASUNET_20240614003_231416_0000[1].pptx
PPTX
CN_Unite_1 AI&DS ENGGERING SPPU PUNE UNIVERSITY
PPTX
Petroleum Refining & Petrochemicals.pptx
PPTX
Amdahl’s law is explained in the above power point presentations
PPTX
ASME PCC-02 TRAINING -DESKTOP-NLE5HNP.pptx
PDF
LOW POWER CLASS AB SI POWER AMPLIFIER FOR WIRELESS MEDICAL SENSOR NETWORK
PDF
Java Basics-Introduction and program control
PPT
Chapter 1 - Introduction to Manufacturing Technology_2.ppt
PDF
Cryptography and Network Security-Module-I.pdf
PDF
UEFA_Carbon_Footprint_Calculator_Methology_2.0.pdf
Management Information system : MIS-e-Business Systems.pptx
AUTOMOTIVE ENGINE MANAGEMENT (MECHATRONICS).pptx
20250617 - IR - Global Guide for HR - 51 pages.pdf
Computer organization and architecuture Digital Notes....pdf
Module 8- Technological and Communication Skills.pptx
Information Storage and Retrieval Techniques Unit III
ai_satellite_crop_management_20250815030350.pptx
Computer System Architecture 3rd Edition-M Morris Mano.pdf
August 2025 - Top 10 Read Articles in Network Security & Its Applications
CONTRACTS IN CONSTRUCTION PROJECTS: TYPES
PRASUNET_20240614003_231416_0000[1].pptx
CN_Unite_1 AI&DS ENGGERING SPPU PUNE UNIVERSITY
Petroleum Refining & Petrochemicals.pptx
Amdahl’s law is explained in the above power point presentations
ASME PCC-02 TRAINING -DESKTOP-NLE5HNP.pptx
LOW POWER CLASS AB SI POWER AMPLIFIER FOR WIRELESS MEDICAL SENSOR NETWORK
Java Basics-Introduction and program control
Chapter 1 - Introduction to Manufacturing Technology_2.ppt
Cryptography and Network Security-Module-I.pdf
UEFA_Carbon_Footprint_Calculator_Methology_2.0.pdf

Experimental and simulation approach of cooling system in 3-phase inverter using extended surface

  • 1. International Journal of Power Electronics and Drive Systems (IJPEDS) Vol. 13, No. 4, December 2022, pp. 2313~2323 ISSN: 2088-8694, DOI: 10.11591/ijpeds.v13.i4.pp2313-2323  2313 Journal homepage: http://ijpeds.iaescore.com Experimental and simulation approach of cooling system in 3-phase inverter using extended surface Agus Mukhlisin1,2 , Prisma Megantoro1 , Estiko Rijanto3 , I Nyoman Sutantra2 , Lilik Jamilatul Awalin1 , Yoga Uta Nugraha1 , Indra Sidharta2 1 Electrical Engineering Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia 2 Mechanical Engineering Department, Institut Teknologi Sepuluh Nopember, Surabaya, Indonesia 3 Research Center for Smart Mechatronics, National Research and Innovation Agency (BRIN), Jakarta Pusat, Indonesia Article Info ABSTRACT Article history: Received Apr 19, 2022 Revised Aug 1, 2022 Accepted Aug 17, 2022 Overheating is a failure mode that significantly affects the reliability of electronic devices. All electronic devices, including a 3-phase inverter driving a traction motor, produce heat dissipation. Heat dissipation needs to be controlled with cooling to prevent overheating. Overheating can be avoided by increasing cooling or reducing heat dissipation. Heat dissipation in the 3-phase inverter is caused by the internal resistance of the metal– oxide–semiconductor field-effect transistor (MOSFET), switching loss, and other factors. Cooling for the 3-phase inverter can use water coolant or air coolant. The cooling system is based on the amount of heat dissipation produced. Cooling of a 3-phase inverter can use air coolant with the addition of an extended surface area in the heat sink. The heat sink uses aluminum material, often called pin fin. There are kinds of aluminum available in the market. We calculated heat generation based on the MOSFET's internal resistance, switching loss, and other factors. We validated the simulation results experimentally using a thermal camera. Thus, we could find an optimal number, dimensions, and aluminum type of fin for the cooling system in the 3-phase inverter. Keywords: 3-phase inverter Cooling system Extended surface Heat dissipation Heat sinks This is an open access article under the CC BY-SA license. Corresponding Author: Agus Mukhlisin Department of Electrical Engineering, Faculty of Advanced Technology and Multidiscipline Airlangga University Gedung Kuliah Bersama, Campus C, Mulyorejo, Surabaya 60115, Indonesia Email: agus.mukhlisin@ftmm.unair.ac.id 1. INTRODUCTION Electronic component reliability is influenced by a cooling system [1]. Each electronic component will produce heat dissipation, including the 3-phase inverter for a motor controller. Heat dissipation needs to be controlled with a cooling system to prevent overheating. It can also be done by reducing the amount of heat dissipation produced, such as using the right switching strategy [2]–[4], using field-oriented control (FOC) [5], [6], using intelligent control [7]–[11], and the others. The most dominant heat source in the motor controller is losses in the 3-phase inverter. Part of the 3-phase inverter that will produce heat is the metal– oxide–semiconductor field-effect transistor (MOSFET). One of the MOSFET failure modes is overheating. The heat generated by MOSFETs must be transferred quickly to the cooling media with intelligent cooling using a fan [12], thermoelectric cooling [13], others so that the temperature of the MOSFET does not rise above its ability limit. Heat transfer can occur in several ways, including the following: conduction, convection, and radiation [14]–[16]. This paper will
  • 2.  ISSN: 2088-8694 Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323 2314 discuss the MOSFET cooling method using an extended surface (pin fin). The extended surface area of heat transfer improves the thermal efficiency [17], [18]. Computational fluid dynamic (CFD) simulation can be used to obtain a cooling method that suits your needs [19], [20]. This method can help to analyze the cooling method without spending time doing experiments [21]. The numerical method can show the temperature distribution and how much heat must be cooled [22]. If the final design has been obtained, this design can be used as a reference in making a prototype. This work presents optimized heat transfer from MOSFET in a 3-phase inverter by using extended surfaces with some pin fins and dimensions. The calculation of heat loss in the 3-phase inverter is used as the input in the simulation. The simulation result is validated with the experimental results. The validation is carried out based on the temperature distribution and the temperature's value at the same point. If the error tolerance that occurs does not exceed 5%, we consider that an accurate simulation is obtained. Therefore, it can be used to optimize the number of fins and the dimensions of the fins. Our target is that the 3-phase inverter can operate normally at a power of 30 kW. 2. METHOD First, we performed a heat loss calculation in 3-phase inverter for parameter input in the simulation. Next, the simulation results were validated with the experimental results, where the temperature distribution was obtained using a thermal camera. After the simulation results are validated, it is possible to optimize the heat sink according to the heat loss. As a result, the optimal heat sink configuration is obtained with the controller not overheating. Heat losses can be seen in the power flow diagram in Figure 1. Based on Figure 1, the 3-phase inverter losses consist of conduction loss pulse-code modulation PCM, switching loss part submission warrant PSW, dead time loss PD, and gate charge loss PG, as given in (1) to (4) [23], [24]. 𝑃𝐶𝑀 = 1 𝑇𝑆𝑊 ∫ 𝑝𝐶𝑀(𝑡)𝑑𝑡 𝑇𝑆𝑊 0 = 1 𝑇𝑆𝑊 ∫ 𝑅𝐷𝑆𝑜𝑛. 𝑖𝐷 2(𝑡)𝑑𝑡 = 𝑅𝐷𝑆𝑜𝑛. 𝐼𝐷𝑟𝑚𝑠 2 𝑇𝑆𝑊 0 (1) 𝑃𝑆𝑊 = 1 2 𝑉𝐼𝑁. 𝐼𝑂. (𝑡𝑟 + 𝑡𝑓). 𝑓𝑆𝑊 (2) 𝑃𝐷 = 𝑉𝐷. 𝐼𝑂. (𝑡𝐷𝑟 + 𝑡𝐷𝑓). 𝑓𝑠𝑤 (3) 𝑃𝐺 = (𝑄𝑔−𝐻 + 𝑄𝑔−𝐿). 𝑉 𝑔𝑠. 𝑓𝑠𝑤 (4) The symbol 𝑅𝐷𝑆𝑜𝑛 is the MOSFET drain-source internal resistance (Ω), TSW is switching time (s), iD is the current flowing in drain (A), and IDrms is the current flowing from source to drain (A). Vin and Io represent the battery voltage (V) and the output current (A). tr is the rise time of the MOSFET (s), tf is the fall time of the MOSFET (s), and fSW is the switching frequency (Hz). VD, tDr, and tDf denote the low-side MOSFET – diode forward-voltage (Volt), dead-time at rising (s), and dead-time at falling (s). Vgs is the gate drive voltage of the MOSFET (V), Qg – H is the high side MOSFET gate electric charge (C), and Qg – L is the low side MOSFET gate electric charge (C). This controller uses the type of MOSFET IRFB4115 with the specifications shown in Figure 2. Figure 1. Power flow 3-phase inverter Figure 2. MOSFET IRFB4115 specification [25]
  • 3. Int J Pow Elec & Dri Syst ISSN: 2088-8694  Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin) 2315 This controller has the following target specifications; i) maximum voltage operation 150 VDC, ii) maximum power 30 KW, and iii) using an air-cooling system. We first made a simplified 3D model of the 3- phase inverter for numerical analysis using CFD software, as shown in Figure 3. The heat source on the 3- phase inverter is the MOSFET component. So that heat transfer can be analyzed with thermal resistance as shown in Figure 4. Figure 3. Simplification of 3D model of the 3-phase inverter Figure 4. Thermal resistance model controller The heat source from the MOSFET moves either up to the controller cover or down to the heat sink. The heat transfer downward will be more dominant than upward. The controller cover traps the air contained inside the controller container, so it does not allow air circulation. The heat transmitted downward propagates through the thermal pad. A thermal pad is a component that serves as an insulator between the base metal of the MOSFET and the heat sink. The use of the thermal pad must have good thermal conductivity as well as a good insulator. The motor controller failure can also be caused by improper thermal pad selection. If the heat generated by the MOSFET is not transferred to the heat sink quickly, the MOSFET will experience overheating, resulting in failure. Similarly, if the thermal pad cannot isolate the metal base of the MOSFET, a short circuit will occur, and the controller will fail. Generally, the thermal pad uses silicone material. The heat transfer then continues to the heat sink fin by conduction. The equation of heat transfer by conduction is given in [14], [16]: 𝑞𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑜𝑛 = 𝑘. 𝐴. ∆𝑇 (5) The type of heatsink material affects the quality of heat transfer: the more significant the heat sink material's thermal conductivity (k), the more remarkable its heat propagation ability. Our heat sink uses 1060 series
  • 4.  ISSN: 2088-8694 Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323 2316 aluminum material, often found in the local market. Convection sends the heat transfer into the air around the heat sink fin. The equation of heat transfer by convection is given by (6) to (8) [14], [16]: 𝑞𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛 = ℎ. 𝐴. ∆𝑇 (6) 𝑁𝑢 = ℎ.𝐿 𝑘𝑓 = 𝑓(𝑅𝑒𝐿, Pr) (7) 𝑅𝑒𝐿 = 𝜌.𝑉.𝐿 𝜇 (8) The symbols h and A represent the convection coefficient (W/m2 K) and the heat transfer area (m2 ). L is the length of material (meters), kf is the thermal conductivity of the fluid (W/mK), and Nu is the Nusselt number. Re and Pr denote Reynold’s number and Prandtl’s number. ρ is the density of air (kg/m3 ), V is the velocity of air (m/s), and μ is the friction coefficient. Convection heat transfer can be optimized by increasing the convection coefficient and cooling area. Increasing the convection coefficient can be done by increasing the speed of air flowing at the fin heat sink, for example, by adding a fan, directing the air with a guide fan so that air can be directed precisely at the heatsink. The cooling area can be enlarged by using variations in dimension and the number of fins on the heatsink. The quality of the simulation results is determined by the dimensions and meshing used. Our meshing uses a curvature-based parameter mesh with total nodes of 615,925, using Jacobian points by 4 points. Table 1 lists material types used in the motor controller cooling simulation. Figure 5 shows the meshing model. Table 1. Material types for cooling simulation Component Materials Thermal conductivity (W/mK) Fin heat sink Aluminum 1060 200 Base plat Aluminum 1060 200 Thermal pad Silicone 3.2 MOSFET plastic body Plastic Silicone 124 MOSFET base metal Aluminum 5052 137 Figure 5. Meshing model Connection contact set using bounded type. The contact that needs to be defined is between MOSFET with the thermal pad and the thermal pad with the heat sink. This step ensures that the connections between different materials are perfectly connected. We conducted numerical simulation using CFD with the following conditions: − Surrounding temperature: 28 o C − Controlled variable: Convection coefficient (h) − The source of the heat loss: the MOSFET First, we calculated heat loss when the motor controller is under a 1.5 kW load as show in Table 2. Afterward, the simulation can be run, and temperature distribution can be obtained. Simulation results will be validated using experimental results. Experiments with a maximum load of 30 kW continuously (30 minutes) cannot be done because the dynamometer will overheat. So, the experiment was carried out with a continues load of 1.5 kW. Furthermore, the results of these experiments can be used to validate the simulation results. The validated simulation results can then be used to test the maximum loading, so that the temperature
  • 5. Int J Pow Elec & Dri Syst ISSN: 2088-8694  Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin) 2317 distribution is obtained when the maximum loading to obtain the temperature distribution under the maximum loading. This simulation can also be used to optimize the fin on the heat sink so that cooling is obtained as needed with minimum fin dimensions. Figure 6 shows the experimental setup with (a) schema of experiment and (b) installed equipment in this paper. It is equipped with an eddy current dynamometer paired with a 30 KW motor. The temperature sensor uses a 10K thermistor recorded in a datalogger (USB to CAN Datalogger). An oscilloscope with four channels is used to measure voltage and current, both input and output. A thermal camera is used to obtain the temperature distribution. Table 2. Calculated heat loss in the motor controller under 1.5 kW load Type Heat Loss Parameter Value Conduction Loss RDS(on)= 11 mΩ; ID = 8.15 A; Vin = 80.34 V; 𝑡𝑟=73 ns; 𝑡𝑓=39 ns;𝑓𝑆𝑊=20kHz; td(on) = 18 ns; td(off) =41ns; Vgs = 23V; Qg =120 nC. 0.617 W Switching Loss 0.733 W Deadtime Loss 0.772 W Gate Charge Loss 0.055 W Total Heat Loss 2.177 W (a) (b) Figure 6. Experiment setup (a) schema and (b) equipment 3. RESULTS AND DISCUSSION 3.1. Simulation validation through experiment The load was kept constant at 1.5 kW with a motor speed of 950 rpm and variations in wind speed. We changed wind speed by controlling a fan power to obtain speed 1, speed 2, and speed 3. CAN to USB datalogger was used to obtain data, among others: temperature of heat sink, phase current, battery current, battery voltage, and the motor speed. One-time data collection took 35 minutes or until the heat sink temperature got steady. The result of the experiment is shown in Figure 7. Figure 7. Motor speed and heat sink temperature
  • 6.  ISSN: 2088-8694 Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323 2318 The heat sink temperatures in Figure 7 were measured using a thermistor sensor. Yellow, red, blue, and green line curves are the heat sink temperatures with variations in wind speed 1, speed 2, speed 3, and natural convection. Each line's steady-state temperature is yellow at 33.37 o C, red at 33.96 o C, blue at 34.92 o C, and green at 39.72 o C. The greater wind speed yields a lower temperature of the heat sink. This is in accordance with (6), (7), and (8). The trend line that occurs tends to be steady at a specific time. This indicates that the cooling is sufficient so that overheating will not occur. The temperature distribution can be captured in this steady state with a thermal camera. Figures 8 and 9 show the temperature distribution results obtained from the simulation and experiment. Figure 8 shows the temperature distribution using h=40 W/m2 K, both experimentally and in simulation results. Figure 8(a) shows the simulation results where the thermistor sensor laying point shows a temperature of 33.1 o C. While Figure 8(b) is a photo taken by the thermal camera where the location of the thermistor shows the value of 33.3 o C. These results indicate that the simulation results have a 0.6% error. The temperature distribution in general also has similarities. Under continuous loading of 1.5 kW, the highest temperature occurs in the DC-DC converter, yielding 37 o C. (a) (b) Figure 8. Temperature distribution with h=40 W/m2 K (a) simulation result and (b) experiment result Figure 9 shows the temperature distribution using h=33 W/m2 K for both experimental and simulation results. Figure 9 (a) shows the simulation results where the thermistor sensor laying point shows a temperature of 34.3 o C. In comparison, Figure 9 (b) is a photo taken by the thermal camera where the location of the thermistor shows the value of 34.9 o C. The result indicates that the simulation results have a 1.75% error. Based on the results shown in Figures 8 and 9, the simulation results can be accepted with an error tolerance of ± 1.75%. Furthermore, the simulation method can be used to evaluate the maximum loading of the 3-phase inverter, i.e., with a load of 30 kW. Table 3 summarizes the data obtained from the simulation and experimental results. (a) (b) Figure 9. Temperature distribution with h=33 W/m2 K (a) simulation result and (b) experiment result
  • 7. Int J Pow Elec & Dri Syst ISSN: 2088-8694  Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin) 2319 Table 3. Comparison of simulation with experiment results Coefficient of Convection Experiment Simulation Percentage of Error h= 33 W/ m2 K 34.9 o C 34.3 o C 1.75 % h= 40 W/ m2 K 33.1 o C 33.3 o C 0.6 % 3.2. Optimization of the heat sink under maximum load We calculated heat loss produced by each MOSFET when the 3-phase inverter operates under the maximum load 30 kW. Table 4 shows the parameter values and the heat loss values. Since we have 48 MOSFETs, the total heat loss is 1,048.27 Watts. We used this heat loss value as the input to our CFC model and conducted a numerical simulation. Figure 10 shows the temperature distribution of the heat sink under a maximum load of 30 kW. The maximum temperature is 98.9 o C, so the controller does not overheat. The temperature where the MOSFET is situated shows a value of 38.5 o C. Since the temperature is significantly lower than the critical temperature of the MOSFET, the heat sink fin is not required in this area. Table 4. Heat loss of each MOSFET in the motor controller under 30 kW Load Type Heat Loss Parameter Value Conduction Loss RDS(on)= 11 mΩ; ID = 65.16 A; Vin = 80.34 V; 𝑡𝑟=73 ns; 𝑡𝑓=39 ns;𝑓𝑆𝑊=20kHz; td(on) = 18 ns; td(off) =41ns; Vgs = 23V; Qg =120 nC. 14.491 W Switching Loss 3.551 W Deadtime Loss 3.742 W Gate Charge Loss 0.055 W Total Heat Loss 21.839 W Figure 10. Simulation result under maximum load 30 kW Furthermore, we optimized the number and length of the heat sink fins with maximum thermal load and h=60 W/m2 K. The original length of ach heat sink fin is 20 mm. Figure 11 shows the simulation results. Figure 11(a) shows the temperature distribution using 15 pieces of heat sink fins. The highest temperature of the MOSFET is 103 o C. The temperature is lower than the critical operation temperature of the MOSFET, which is 125 o C. Figure 11(b) shows the temperature distribution using 12 pieces of heat sink fins. The highest temperature of the MOSFET is 107 o C. The result indicates that reducing of fin number to 12 pieces yielding maximum temperature that is lower than the critical temperature of the MOSFET. Further optimization is carried out by reducing the fin length to 12 mm and yielding the temperature distribution shown in Figure 11(c). The maximum temperature of the MOSFET in this condition is 117 o C, which is lower than the critical temperature of MOSFET. However, there is an increase in temperature compared to the result of 12 pieces 20 mm fin. Therefore, the controller under a maximum load of 30 kW is optimally cooled using the configuration of 12 pieces of fins with 12 mm length, following Figure 11(c).
  • 8.  ISSN: 2088-8694 Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323 2320 (a) (b) (c) Figure 11. Optimized heat sink in the number and length of the fins: (a) 15 pieces of 20 mm long, (b) 12 pieces of 20 mm long, and (c) 12 pieces of 12 mm long
  • 9. Int J Pow Elec & Dri Syst ISSN: 2088-8694  Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin) 2321 4. CONCLUSION Thermal analysis of the controller cooling system was carried out using modeling and simulation. The simulation was validated experimentally with an error of ± 1.75%. Furthermore, the simulation method can be used to optimize the cooling design for a maximum loading of 30 kW. The results indicate that the controller does not experience overheating, as demonstrated by the MOSFET temperature. Further optimization, through simulation was carried out in order to reduce the size and the number of the fin. The simulation result of the optimized heat sink suggests that the maximum temperature of the MOSFET is lower than its critical temperature. Hence, the controller did not experience overheating. ACKNOWLEDGEMENTS Author thanks to Faculty of Advanced Technology and Multidiscipline, Airlangga University and Lembaga Pengelola Dana Pendidikan (LPDP) in most cases, sponsor and financial support acknowledgments. This program support Sustainable Development Goals (SDGs) for affordable and clean energy also goal 9 for Industry, Innovation and Infrastructure. REFERENCES [1] E. Laloya, Ó. Lucía, H. Sarnago and J. M. Burdío, “Heat management in power converters: from state of the art to future ultrahigh efficiency systems,” IEEE Transactions on Power Electronics, vol. 31, no. 11, pp. 7896-7908, 2016, doi: 10.1109/TPEL.2015.2513433. [2] N. A. Ramli, A. Jidin, Z. Rasin, and T. Sutikno, “Reduction of total harmonic distortion of three-phase inverter using alternate switching strategy,” International Journal of Power Electronics and Drive Systems, vol. 12, no. 3, pp. 1598-1608, 2021, doi: 10.11591/ijpeds.v12.i3.pp1598-1608. [3] D. Christen and J. Biela, “Analytical switching loss modeling based on datasheet parameters for mosfets in a half-bridge,” IEEE Transactions on Power Electronics, vol. 34, no. 4, pp. 3700-3710, 2019, doi: 10.1109/TPEL.2018.2851068. [4] A. V. Sant and K. R. Rajagopal, “PM synchronous motor speed control using hybrid fuzzy-pi with novel switching functions,” IEEE Transactions on Magnetics, vol. 45, no. 10, pp. 4672-4675, 2009, doi: 10.1109/TMAG.2009.2022191. [5] V. T. Ha, N. T. Lam, V. T. Ha, and V. Q. Vinh, “Advanced control structures for induction motors with ideal current loop response using field oriented control,” International Journal of Power Electronics and Drive System, vol. 10, no. 4, pp. 1758-1771, 2019, doi: 10.11591/ijpeds.v10.i4.pp1758-1771. [6] M. Hasoun, A. el Afia, M. Khafallah, and K. Benkirane, “Field oriented control based on a 24-sector vector space decomposition for dual three-phase pmsm applied on electric ship propulsion,” International Journal of Power Electronics and Drive Systems, vol. 11, no. 3, pp. 1175-1187, 2020, doi: 10.11591/ijpeds.v11.i3.pp1175-1187. [7] D. Brahim, B. Mokhtar, and L. Salah, “Improvement of adaptive fuzzy control to adjust speed for a doubly fed induction motor drive (DFIM),” International Journal of Power Electronics and Drive Systems, vol. 11, no. 1, pp. 496-504, 2020, doi: 10.11591/ijpeds.v11.i1.pp496-504. [8] E. Sabouni, B. Merah, and I. K. Bousserhane, “Adaptive backstepping controller design based on neural network for pmsm speed control,” International Journal of Power Electronics and Drive Systems, vol. 12, no. 3, pp. 1940-1952, 2021, doi: 10.11591/ijpeds.v12.i3.pp1940-1952. [9] A. Nait Seghir, T. Henni and M. Azira, “Fuzzy and adaptive fuzzy PI controller-based Vector control for permanent magnet synchronous motor,” 10th Ieee International Conference On Networking, Sensing And Control (ICNSC), 2013, pp. 491-496, doi: 10.1109/ICNSC.2013.6548788. [10] Jae Sung Yu, Sun Mo Hwang and Chung-Yuen Won, “Performances of fuzzy-logic-based vector control for permanent magnet synchronous motor used in elevator drive system,” 30th Annual Conference of IEEE Industrial Electronics Society, IECON, 2004, vol. 3, pp. 2679-2683, doi: 10.1109/IECON.2004.1432229. [11] N. J. Patil, R. H. Chile and L. Waghmare, “Fuzzy adaptive controllers for speed control of PMSM drive,” International Journal of Computer Applications, vol. 1, no. 11, pp. 975-8887, 2010, doi: 10.5120/233-387. [12] A. Deep, J. Singh, Y. Narayan, S. Chatterji and L. Mathew, “Intelligent cooling system for three level inverter,” Communication, Control and Intelligent Systems (CCIS), pp. 286-289, 2015, doi: 10.1109/CCIntelS.2015.7437925. [13] U. Sanver, E. Yavuz and C. Eyupoglu, “An electronic control unit for thermoelectric cooling,” IEEE Conference of Russian Young Researchers in Electrical and Electronic Engineering (EIConRus), 2019, pp. 141-145, doi: 10.1109/EIConRus.2019.8656871. [14] F. P. Incropera, D. Witt, T. L. Bergman, and A. S. Lavine, Fundamentals of heat and mass transfer, Sixth Edition. New Jersey: John Wiley & Sons, 2006. [15] A. Blinov, D. Vinnikov, and T. Lehtla, “Cooling methods for high-power electronic systems,” Scientific Journal of Riga Technical University, Power and Electrical Engineering, vol. 29, no. 1, pp. 79-86, 2011, doi: 10.2478/v10144-011-0014-x. [16] Y. A. Çengel, Heat Transfer: A Practical Approach, New York, NY, USA: McGraw-Hill, 2003. [17] S. S. Kang, “Advanced cooling for power electronics,” 7th International Conference on Integrated Power Electronics Systems (CIPS), 2012, pp. 1-8. [18] A. Driss, S. Maalej and M. C. Zaghdoudi, “Thermal modeling of the cooling of a power MOSFET by heat pipes,” International Conference on Engineering & MIS (ICEMIS), 2017, pp. 1-6, doi: 10.1109/ICEMIS.2017.8273067. [19] U. SanAndres, G. Almandoz, J. Poza and G. Ugalde, “Design of cooling systems using computational fluid dynamics and analytical thermal models,” IEEE Transactions on Industrial Electronics, vol. 61, no. 8, pp. 4383-4391, Aug. 2014, doi: 10.1109/TIE.2013.2286081. [20] Y. Bulut and K. Pandya, “Thermal modeling for power MOSFETs in DC/DC applications,” 5th International Conference on Thermal and Mechanical Simulation and Experiments in Microelectronics and Microsystems, 2004, pp. 429-433, doi: 10.1109/ESIME.2004.1304074. [21] O. Spuiber, M. M. De Souza, E. M. S. Narayanan and S. Krishnan, “Analyses of a COOL-MOSFET,” CAS '99 Proceedings. 1999 International Semiconductor Conference (Cat. No.99TH8389), 1999, pp. 131-134 vol. 1, doi: 10.1109/SMICND.1999.810446.
  • 10.  ISSN: 2088-8694 Int J Pow Elec & Dri Syst, Vol. 13, No. 4, December 2022: 2313-2323 2322 [22] M. Langer, Z. Lisik and J. Podgorski, “Numerical simulation of cool MOS transistor,” Experience of Designing and Applications of CAD Systems in Microelectronics. Proceedings of the VI-th International Conference, CADSM, 2001, pp. 303-304, doi: 10.1109/CADSM.2001.975849. [23] G. Dusan, M. Purschel, and A. Kiep, “MOSFET power losses calculation using the data-sheet parameters,” Infineon application note, pp. 1-23, 2006. [24] ROHM Semiconductor, “Application Note: Calculation of Power Loss (Synchronous),” 2016. [Online]. Available: http://rohmfs.rohm.com/en/products/databook/applinote/ic/power/switching_regulator/power_loss_appli-e.pdf. [25] International IOR Rectifier, “IRFB4115PbF HEXFET Power MOSFET Datasheet,” 2014. [Online]. Available: https://www.infineon.com/dgdl/irfb4115pbf.pdf?fileId=5546d462533600a401535615ba6a1e0f. BIOGRAPHIES OF AUTHORS Agus Mukhlisin is a lecturer in Electrical Engineering Department at Faculty of Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia. He received a Bachelor’s degree and Master’s degree from Institute Technology of Sepuluh November, Surabaya, Indonesia in 2015 and 2017. He is currently pursuing Doctoral Program in Institute Technology of Sepuluh November, Surabaya, Indonesia. His Dissertation Topic was about Intelligent Vector Control in Motor Controller. His research interests include System and Control of Electric Vehicle, Battery Management System, Charging Station, and Dynamic System Simulation. He can be contacted at email: agus.mukhlisin@ftmm.unair.ac.id. Prisma Megantoro is a lecturer in Electrical & Electronic Engineering Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia. He received the bachelor degree and master degree from Universitas Gadjah Mada, Yogyakarta, Indonesia in 2014 and 2018. His current research is focused on solar photovoltaic technology, embedded system, and internet of things. He can be contacted at email: prisma.megantoro@ftmm.unair.ac.id. Estiko Rijanto has worked for Research Center for Smart Mechatronics - the National Research and Innovation Agency (BRIN) - Indonesia since 2021. From 2002 to 2021, he worked at Research Center for Electrical Power and Mechatronics, Indonesian Institute of Sciences (LIPI). In 1987 he enrolled at Institut Teknologi Bandung (ITB). He completed his B.Eng. degree at Tokyo University of Agriculture and Technology (TUAT), Tokyo - Japan, in 1993. He received M. Eng. and Dr. Eng. degrees from the same university in 1995 and 1998. In 2000, he completed his book on “Robust Control: Theory for Application” during his post- doctoral program at TUAT-VBL. In 2013, he was inaugurated as the Research Professor on Applied Control Systems by LIPI after presenting his book “Integration of control and information systems for industrial competitiveness,”. His research interests include control systems and their applications for mechatronics, electric vehicles, power generations, renewable energies, and battery management systems. He can be contacted at email: esti003@brin.go.id. I Nyoman Sutantra is a lecturer in Mechanical Engineering Department, Institut Teknologi Sepuluh November, Surabaya, Indonesia. He received a Bachelor’s degree from Institute Technology of Sepuluh November, Surabaya, Indonesia. He received a Master and Doctoral degree from The University of Wisconsin – Madison in mechanical department. He was inaugurated the Professor on mechanical engineering. His research interests include automotive technology, electric vehicles, hybrid vehicle and steering control. He can be contacted at email: tantra@me.its.ac.id.
  • 11. Int J Pow Elec & Dri Syst ISSN: 2088-8694  Experimental and simulation approach of cooling system in 3-phase inverter … (Agus Mukhlisin) 2323 Lilik Jamilatul Awalin is a lecturer in Electrical & Electronic Engineering Department, Faculty of Advanced Technology and Multidiscipline, Airlangga University, Surabaya, Indonesia. She received the bachelor degree from electrical engineering, Widya Gama Malang, Indonesia in 1977. She received the master degree from Institute Technology of Sepuluh November, Surabaya, Indonesia in 2004. Finally, she gets the doctoral degree from University of Malaya, Malaysia in 2014. Her current research is focused on fault localization, smart grid technology, transmission and distribution network, load shedding, coordination protection and optimization. She can be contacted at email: lilik.j.a@ftmm.unair.ac.id. Yoga Uta Nugraha is as a lecturer at Electrical Engineering program in the Faculty of Advanced Technology and Multidicipline, Universitas Airlangga, Surabaya, Indonesia. He was born in Sukoharjo, Central of Java. He is a Doctoral candidate in Electrical Engineering Department of Sepuluh Nopember Institute of Technology Surabaya. He has a Bachelor (2015) and Master (2020) of Electrical Engineering from Sepuluh Nopember Institute of Technology Surabaya. His current research is focused on integration on Electric Vehicle, electric machine, permanent magnet motor, new and renewable energy. He can be contacted at email: yoga.uta.n@ftmm.unair.ac.id. Indra Sidharta is a lecturer in Mechanical Engineering Department, Institut Teknologi Sepuluh November, Surabaya, Indonesia. He received a Bachelor’s degree from Institute Technology of Sepuluh November, Surabaya, Indonesia in 2004. He received a Master degree from Hamburg University of Technology in field Materials Science, 2010. His current research is focused on metal matrix composite, metallurgy, materials of battery. He can be contacted at email: indra.sidharta@gmail.com.