Optimize an Inverter Liquid Cooling System

This example uses:

In this example you analyze the performance of a liquid cooling system for a three-phase inverter. You run detailed and reduced models (ROM) to find the steady-state temperatures and losses. You compute the optimal size of the heatsink that maximizes the inverter efficiency and minimizes lifetime cost.

Model Overview

Open the model sscv_inverter_liquid_cooling.

model = 'sscv_inverter_liquid_cooling';
open_system(model)
sscv_inverter_liquid_cooling_params;

To drive the Load block, the inverter converts the DC power from the high-voltage battery into three-phase AC power. Conduction losses, switching losses, and reverse recovery losses generate heat in the case. Liquid cooling is effective to dissipate heat in the order of kilowatts. Liquid coolant flows in the inverter case to exchange heat. The Electric Pump System provides the hydraulic power to produce the flow. Then, the hot liquid flows into a heatsink to dissipate the heat. The heatsink is also called radiator. Fans blow air into the radiator to increase heat dissipation via forced convection.

The Three-phase Inverter has a Detailed Inverter variant that consists of six IGBT switching devices and the corresponding body diodes:

The Electric Pump System comprises a 12V DC motor and its H-bridge driver and a small centrifugal pump. The motor regulates the pump flow rate by using PWM switching and a PID controller:

Import Device Parameters

To import datasheet parameters for the IGBTs and Body Diodes in the Detailed inverter model, including turn-on and turn-off switching losses for the IGBT blocks, and reverse recovery losses for the Diode blocks, use the ee_importDeviceParameters utility function.

sscv_inverter_liquid_cooling_import

The plots above show the device characteristics for the imported IGBT and Diode blocks, including switching and reverse recovery losses.

Run a Detailed Simulation

Simulate 4 AC cycles to observe load voltage and current.

set_param([model, '/Three-phase inverter'],'LabelModeActiveChoice', 'Detailed')
stopTime = 4/60; %#ok<NASGU> 
out = sim(model);
figure();
Vabc = out.simlog.Load.N.V.series;
Iabc = out.simlog.Load.wye_impedance.I.series;
subplot(2,1,1);
plot(Vabc.time, Vabc.values('V'));
ylabel('Phase voltages [V]')
grid on
subplot(2,1,2);
plot(Iabc.time, Iabc.values('A'));
xlabel('Time [s]')
ylabel('Phase currents [A]')
grid on

Observe the inverter successfully produces a three-phase current at 60 Hz. It achieves electrical AC steady-state after an initial transient cycle.

Run a Reduced Thermal Simulation

The electrical characteristic times in the inverter subsystem are in the order of milliseconds, but the thermal characteristic times are in the order of minutes. If you are only interested in the thermal steady-state of the system for a particular cooling setup and a fixed value of average conduction and switching losses, you can use the Reduced - Use only for thermal steady-state estimation variant of the Three-phase inverter subsystem.

To run a 1 hour simulation of thermal model with constant conduction and switching losses, but without voltages or currents, at the MATLAB Command Window, enter:

set_param([bdroot, '/Three-phase inverter'],'LabelModeActiveChoice', 'Reduced - Use only for thermal steady-state estimation');
stopTime = 3600;
out = sim(model);

Plot the temperature of the case, heatsink, and coolant parts:

figure();
Tcase = out.simlog.Case_Heat_Exchanger.H.T.series;
Theatsink = out.simlog.Heatsink.T_hs.series;
Tcoolant = out.simlog.Tank.T_I.series;
hold on
plot(Tcase.time, Tcase.values('degC'));
plot(Theatsink.time, Theatsink.values('degC'));
plot(Tcoolant.time, Tcoolant.values('degC'));
hold off
grid on
legend({'Case', 'Heatsink', 'Coolant'})
xlabel('Time [s]')
ylabel('Temperature [degC]')

You observe that it takes almost one hour to reach thermal steady state.

Estimate Thermal Steady-State Iteratively

To find the accurate steady-state losses and temperatures faster, first compute the accurate losses at a given temperature by using the Detailed Inverter option of the Three-Phase Inverter subsystem. Then compute the accurate steady-state temperatures at a given loss by using the Reduced option of the Three-Phase Inverter subsystem.

Finally, iterate between these two processes to estimate the thermal steady-state.

sscv_inverter_liquid_cooling_thermal_ss

Iteration 1 --------
Max Junction temperature difference to previous run = 51.5222 degC
Iteration 2 --------
Max Junction temperature difference to previous run = 8.9189 degC
Iteration 3 --------
Max Junction temperature difference to previous run = 1.4989 degC

                                           LoggingNode                                              Power       SwitchingLosses
    __________________________________________________________________________________________    __________    _______________

    {'sscv_inverter_liquid_cooling.Load.wye_impedance'                                       }         31119             0     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_AH.IGBT'      }        67.212        477.15     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_BL.IGBT'      }        66.877        472.13     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_CL.IGBT'      }        66.591         471.3     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_BH.IGBT'      }        66.176        474.24     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_CH.IGBT'      }            66        459.58     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_AL.IGBT'      }        64.828        471.76     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_AL.Body_Diode'}        60.939        242.43     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_BH.Body_Diode'}        60.686        239.21     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_AH.Body_Diode'}        60.536        238.75     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_CH.Body_Diode'}        60.485        238.77     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_BL.Body_Diode'}        60.326        239.49     
    {'sscv_inverter_liquid_cooling.Three_phase_inverter.Detailed_Inverter.IGBT_CL.Body_Diode'}         59.46        232.96     
    {'sscv_inverter_liquid_cooling.Electric_Pump_System.DC_Motor'                            }        7.4587             0     
    {'sscv_inverter_liquid_cooling.Electric_Pump_System.H_Bridge'                            }        0.2134             0     
    {'sscv_inverter_liquid_cooling.Gmin'                                                     }    0.00035504             0     

System efficiency = 86.0959%

You observe that 3 iterations are enough to find a thermal steady-state with a reasonable accuracy of +-2 degC.

Optimize Heatsink Size for Cost and Efficiency

To find the best option that minimizes the cost and maximizes the efficiency, repeat the steady-state iterations for several heatsink sizes. At the MATLAB Command Window, enter:

displayOutputs = false;
sscv_inverter_liquid_cooling_optimize

************ Case 1: Heatsink with NumFins = 60 ************
System efficiency = 84.8203%
************ Case 2: Heatsink with NumFins = 80 ************
System efficiency = 85.5376%
************ Case 3: Heatsink with NumFins = 120 ************
System efficiency = 86.3134%
************ Case 4: Heatsink with NumFins = 200 ************
System efficiency = 86.8824%
************ Case 5: Heatsink with NumFins = 300 ************
System efficiency = 87.1501%
------------------------------------
Time to run 5 cases = 26.1305 min
************************************

You observe that the option with a good blend of high efficiency and low coolant temperature, with a minimum lifetime cost, is a heatsink with around 120 fins.