Next Article in Journal
Pd Catalysts Supported on Bamboo-Like Nitrogen-Doped Carbon Nanotubes for Hydrogen Production
Next Article in Special Issue
EMC Component Modeling and System-Level Simulations of Power Converters: AC Motor Drives
Previous Article in Journal
Data-Driven Signal–Noise Classification for Microseismic Data Using Machine Learning
Previous Article in Special Issue
Electromagnetic Field Tests of a 1-MW Wireless Power Transfer System for Light Rail Transit
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Behavioral Modelling Methodology of Switch-Diode Cell for Power Loss Prediction in Electromagnetic Transient Simulation

1
Department of Electrical and Computer Engineering, University of Manitoba, Winnipeg, MB R3T5V6, Canada
2
Manitoba Hydro International, Winnipeg, MB R3P1A3, Canada
*
Author to whom correspondence should be addressed.
Energies 2021, 14(5), 1500; https://doi.org/10.3390/en14051500
Submission received: 14 January 2021 / Revised: 17 February 2021 / Accepted: 4 March 2021 / Published: 9 March 2021
(This article belongs to the Special Issue Electromagnetic Modeling in Power Electronics)

Abstract

:
Modern wide-bandgap (WBG) devices, such as silicon carbide (SiC) or gallium nitride (GaN) based devices, have emerged and been increasingly used in power electronics (PE) applications due to their superior switching feature. The power losses of these devices become the key of system efficiency improvement, especially for high-frequency applications. In this paper, a generalized behavioral model of a switch-diode cell (SDC) is proposed for power loss estimation in the electromagnetic transient simulation. The proposed model is developed based on the circuit level switching process analysis, which considers the effects of parasitics, the operating temperature, and the interaction of diode and switch. In addition, the transient waveforms of the SDC are simulated by the proposed model using dependent voltage and current sources with passive components. Besides, the approaches of obtaining model parameters from the datasheets are given and the modelling method is applicable to various semiconductors such Si insulated-gate bipolar transistor (IGBT), Si/SiC metal–oxide–semiconductor field-effect transistor (MOSFET), and GaN devices. Further, a multi-dimensional power loss table in a wide range of operating conditions can be obtained with fast speed and reasonable accuracy. The proposed approach is implemented in PSCAD/ Electromagnetic Transients including DC, EMTDC, (v4.6, Winnipeg, MB, Canada) and further verified by the hardware setups including different daughter boards for different devices.

1. Introduction

A power electronics (PE) system plays a key role in the process of efficient energy control, conversion, and management. Power semiconductor devices are the core components in a PE system and have a significant impact on system efficiency, reliability, and cost [1]. For decades, silicon-based devices, such as insulated-gate bipolar transistors (IGBTs) [2], metal-oxide-semiconductor field-effect transistors (MOSFETs) [3], are mainly and widely used in various modern PE applications (e.g., Photovoltaics (PV) [4], Power Factor Correction (PFC) [5], power supply [6], and other power converters [7,8]). However, the PE system performance and efficiency are hindered by Si-based devices due to the fundamental material limits. Recently, wide-bandgap (WBG) devices [9,10,11], such as silicon carbide (SiC) MOSFETs [12], enhancement-mode gallium nitride (eGaN) high-electron-mobility transistors (HEMTs) [13,14], have emerged and gained great popularity due to the superior features of fast switching speed and low switching loss. Thereby, the switching frequency can be further increased bringing the merits of size reduction for magnetic components, high power density, and high efficiency. Whereas the increased power losses of semiconductors are typically the main contributor to total loss especially for high-frequency (HF) applications, and the generated heat energy during switching transition may lead to fatigue failure and affect the reliability [15]. Hence, an accurate power loss model, which is applicable for different semiconductors and provides a deep insight into the switching process, is highly desirable for device selection and PE system optimization.
Currently, the ideal switch or two-state resistances model is typically adopted in most electromagnetic transient (EMT) simulators such as PSCAD/EMTDC and MATLAB/Simulink [16]. This simple model is mainly used to evaluate the overall system response and control strategy, and only the conduction loss is roughly considered. The conduction loss can be directly determined by the output curves in the datasheet, while the switching loss is more complicated and can be measured in the double-pulse test (DPT) [17]. Although DPT is widely used and can achieve high accuracy, it typically involves expensive probes and much peripheral bulky equipment such as a high voltage power supply. Designing a testing board with low parasitics is challenging and also significant for WBG devices due to the fast switching. Recently, several physic-based semiconductor models [18,19,20,21], such as simulation program with integrated circuit emphasis (SPICE) models [22,23], have been proposed to accurately describe the transient behaviors of the devices. However, the geometrical parameters for the model are often not available in the datasheet and thus the applicability of the model is limited. Another type of model (i.e., behavioral model [24,25]) has been developed, which focuses more on the external behaviors of the devices instead of the internal physics. As a result, the complexity is reduced and fast simulation speed can be achieved. It is adequate and widely used for system-level study, but more detailed transient concerns are needed to accurately evaluate the switching performance and estimate the power losses.
To have a better description of the switching transients, a lot of analytical loss models have been proposed [26,27]. Piecewise linearizing the switching process of the device is a commonly used method which enables simple and rapid loss estimation [28]. Whereas, the accuracy is still limited due to the ignorance of the parasitics. To improve it, more comprehensive loss models are developed considering various factors, such as temperature-dependent parameters [29,30], interactions between diode and switch [31], cross-talk issue [32], displacement current [33,34], and non-flat miller plateau [35]. Thereby, the switching loss can be obtained by solving the equivalent circuit for each switching sub-stage. Further, the entire switching process of eGaN HEMT in synchronous buck converter application is presented in [36,37] considering the third quadrant operation with the help of the 2-dimensional electron gas (2DEG). However, these methods are complicated involving huge computational burdens, not to mention the convergence issue. The measurement techniques and loss distribution including the capacitive losses for eGaN HEMT are illustrated in [38,39] and the scalable loss estimation method is further proposed based on the measurements. However, the measured data in the datasheet is typically under specific conditions, which limits the applicability and accuracy.
In a PE system, a power switch is typically paired with a diode as a switch-diode cell (SDC) to provide current commutation [40]. This basic commutation cell as shown in Figure 1a is widely used in PE applications and it consists of the active power switch ( S ), diode ( D ), equivalent circuit voltage ( V dd ) and load current ( I L ) [41]. Note that, four configurations of S , namely Si/SiC MOSFET, Si IGBT and eGaN HEMT, are taken into account in this paper and D can be a single positive-intrinsic-negative (PIN) diode, a Schottky barrier diode (SBD), the body diode of MOSFET or the equivalent diode of eGaN HEMT. During switching transition, power loss is resulted mainly from the switching and conduction losses of S as well as the conduction loss and reverse recovery loss of D . In order to estimate these power losses in a PE system, a generalized behavioral modelling method of switch-diode cell in electromagnetic transient simulation (EMT) is proposed and it is an extension of previous work [42,43,44,45]. There are three technical contributions in this paper comparing with the conventional methods.
  • A generalized behavioral model of SDC is proposed which is realized by dependent sources with passive components considering the impacts of parasitics, the temperature-dependent parameters, and the reverse recovery behavior of D . This model is not limited to a specific device and it is applicable to various devices including Si/SiC MOSFET, Si IGBT, and eGaN HEMT. In addition, most of the model parameters can be obtained from the device datasheets by the curve fitting method and no additional measurement is involved. Based on the specific requirement, the model can be modified and integrated into different simulators accordingly.
  • The switching process of the switch-diode cell in a clamped inductive switching circuit is studied analytically based on the equivalent circuits for each switching sub-stage. Accordingly, the semiconductor model is developed and implemented in PSCAD/EMTDC. The switching analysis in this paper is more comprehensive considering the respective features of different devices such as the tail current of IGBT and the third quadrant operation of eGaN HEMT.
  • A DPT setup was designed for experimental verification. To meets the different requirements of different semiconductors, three daughter boards were specifically designed incorporating with the main control board to characterize various devices and measure the corresponding power losses in a wide range of operating conditions. The simulated results are compared with the experimental results and show good agreements within 10% average error range.

2. Model Description

The simulation procedure of the proposed model is demonstrated in Figure 2. Initially, the device requirements for the desired PE application (e.g., voltage, current, temperature, and frequency) are determined. Based on those requirements, a specific semiconductor device is preliminarily selected for modelling and characterization. According to the device’s datasheet, the key model parameters can be extracted by the curve fitting method. Afterwards, the model parameters together with the operating conditions are input to the proposed device model, and a DPT simulation using the proposed model is further carried out. Subsequently, the transient voltage and current waveforms can be obtained, and simultaneously the power loss of the device can be computed. If the simulated results, in terms of switching transient behaviors and power loss, meet the requirements within the acceptable range, then the semiconductor is eventually selected for this application. Otherwise, it is necessary to reselect another device and evaluate the performance until the design is optimized.
In order to understand the switching behaviors of the SDC, a diode-clamped inductive load circuit (i.e., DPT circuit) is taken as an example which is widely used for device characterization. As shown in Figure 1b, the basic commutation unit consists of two complementary switches, one operates as a freewheeling diode D and the other is an active switch S controlled by the gate drive voltage ( v G ) through the external gate resistance R gext . In a typical hard-switching PE system, S is identified by a positive drain current i d (collector current i c for IGBT) direction matching with the direction of I L . Since the commutation time is sufficiently short, I L and V dd are hardly changing during switching transition and thus they are treated as constant current and voltage sources, respectively.
Note that, the crucial circuit parasitic elements are also included as shown in Figure 1b. All the stray inductances in the power loop including the printed circuit board (PCB) trace and device package inductance are lumped and represented by L s , while the common source inductance ( L cs ) of S is considered separately. In addition, the parasitic capacitances of S include gate-drain capacitance ( C gd ), gate-source capacitance ( C gs ) and drain-source capacitance ( C ds ). Besides, the equivalent capacitance of D ( C F ) denotes for junction capacitance of diode. It should be mentioned that when D is configured as the body diode of a switch rather than a single diode,   C F will be the corresponding parasitic capacitance of the switch. During a switching transition, I L commutates between S and D . When a positive v G is given, the gate-source voltage ( v gs ) will increase to turn on S . Subsequently, i d including the channel current ( i ch ), the gate-drain current ( i gd ) and the drain-source current ( i ds ) starts rising, meanwhile the diode forward current ( i F ) declines gradually. When S is fully turn on, the drain-source voltage ( v ds ) decreases to the on-state voltage and the diode forward voltage ( v F ) rises to V dd . The behavioral models of a SDC as illustrated in Figure 1c are proposed to reproduce the switching behaviors of S and D , respectively. The details of the model descriptions including the active switch and diode model are presented as follows.

2.1. Active Switch Model

As shown in Figure 1c, the proposed active switch model consists of two parts, the gate loop and the power loop. It is noted that L cs is shared by both loops and thus each loop includes one L cs in order to decoupling both loops. Additionally, a dependent voltage source ( v Lcs ) is added in the gate loop to reflect the interactive impact of the current source ( i S ) on L cs as expressed by,
v Lcs = L cs · d i S / d t
  • Gate Loop Part
The external v G is typically flipped between V gon (20 or 15 V) and V goff (−5 or 0 V) based on the specific gate drive requirement of the switch. The device internal gate resistance ( R gint ) is merged into R gext as the total gate resistance ( R G ). Furthermore, the gate related junction capacitances (i.e., C gd and C gs ) are represented by the input capacitance ( C iss ) and an additional dependent voltage source ( v mil ). This equivalent v mil becomes valid only when the miller plateau occurs on v gs during switching transition and its value can be computed by
v mil = v th + i ch / g fs
where v th and g fs stand for threshold voltage and transconductance of a switch, respectively. Additionally, i ch during miller plateau period typically equals to I L which will be discussed in Section 2.3.
It is noted that the gate inductance is neglected here for simplicity although it can introduce a slight delay on v gs . In fact, this delay is mainly resulted from L cs and v Lcs due to the fast change of i d . Besides, the gate drive circuit is normally placed close to S in order to minimize the potential oscillation introduced by the gate inductance and thus this inductance is negligible.
  • Power Loop Part
In this model, i d is represented by i S which is the sum of i ch , i gd and i ds . Note that, most of the time, i S is the same as i ch except for the voltage transition period when a displacement current is introduced due to the process of charging and discharging of the parasitic capacitance. In order to reflect the voltage change during switching transition as well as the on-state voltage ( v on ) of S , an equivalent dependent voltage source ( v S ) is adopted here. The value of v on can be determined by R ds ( on ) with I L in (3) or the saturation voltage ( v cesat ) for the case of IGBT.
v on = R ds ( on ) · I L .
Besides, L cs is also included in the power loop part which is associated with L s to influence the transient waveforms. Thereby, the gate loop and power loop parts are decoupled and their interaction is represented by the equivalent dependent sources instead of nonlinear junction capacitances resulting in a reduction of model complexity.

2.2. Diode Model

The static model of a diode typically can be represented by an ideal diode ( D F ) , a forward resistance ( R F ) and a voltage source ( v F 0 ) based on the forward characteristics in the device datasheet. Typically, v F of the diode can be computed by,
v F = R F · i F + v F 0 .
It should be mentioned that, for the case of eGaN HEMT as D , the diode behavior is realized by 2DEG and thus the calculation of v F is based on the reverse conduction characteristic of the GaN device [44] which is highly affected by the gate drive voltage ( v GF ) of GaN device as shown in (5),
v F ( GaN ) = R Fr · i F + v thF v GF
where R Fr and v thF are on-resistance in the third quadrant and threshold voltage of a GaN device, respectively. Since negative v GF is typically provided to avoid the cross-talk issue, higher v F is thus resulted which will increase the conduction loss of D . Notice that, if a positive v GF is provided enabling the channel fully on, the on-state resistance will be the same value in the first quadrant.
Moreover, the dynamic characteristic of D is described by C F in parallel with an equivalent dependent current source ( i re ) for the reverse recovery behavior of D . When D switches from forwarding conduction to off-state, i F cannot be eliminated immediately and it takes a while to extinguish the excess carriers, this time is called reverse recovery time ( t rr ). The reverse recovery process occurs as soon as i F becomes negative, i re can be expressed by [45],
i re = { d i F / d t · t ,   t < t rm I rm · exp ( ( t t rm ) / τ re ) ,   t t rm ,
where τ re denotes decay time constant and i re reaches the peak current ( I rm ) at time t rm . In addition, the slew rate of diode current ( d i F / d t ) typically keeps the same as the turn-off slew rate of S . It needs to notice that, this reverse recovery behavior commonly exists in PIN diode and body diode of S , while the reverse recovery loss is eliminated for the case of SBD or eGaN HEMT and thus i re can be neglected for simplicity. In fact, in these cases, the effect of C F is the main concern which can introduce displacement current resulting in capacitive loss during switching transition.

2.3. Switching Transient Modelling

The switching process of the SDC in the DPT circuit is thoroughly analyzed based on the switching waveforms and the equivalent circuits as follows.
  • Turn-on transition ( t 0 t 3 )
The typical turn-on waveforms along with power loss information are illustrated in Figure 3 considering the case of PIN or body diode as well as the case of SBD or eGaN HEMT. The equivalent circuits during this period are also provided in Figure 4.
Initially, S is in the off-state, and all the I L flows through D . The corresponding v F can be estimated by (5). At t 0 , the gate charging period begins with a positive V gon applying to v G and C iss is charged up through R G . Subsequently, v gs will increase accordingly with the time constant ( τ iss = R G · C iss ). Note that L cs in the gate loop will prolong the turn-on time causing more power losses.
The current rising period begins when v gs goes beyond v th . During this interval, the conductive channel of S is forming and i S starts rising from zero to I L which can be expressed by
i S = g fs · ( v gs v th ) .
The fast change of i S on one hand, will introduce a negative feedback v Lcs from power loop to gate loop due to L cs to further delay the turn-on process. On the other hand, it will result in a total voltage drop ( v L ) on L s and L cs . Simultaneously, v ds decreased by v L as shown in Figure 3.
As I L commutates from D to S , i S reaches I L at t 2 and v gs will be clamped at v mil . At the same time, i F decreases to zero and D enters into the reverse recovery as shown in Figure 3a. This additional i re will add to i S (= I L + i re ) resulting a current spike and thus a bump in v gs according to (2). When i re reaches I rm , it starts declining and the voltage falling period begins. Subsequently, v ds starts decreasing which is controlled by v S and the corresponding slew rate can be determined by,
d v ds / d t = d v S / d t = ( V gon v mil v Lcs ) / ( C gd · R G ) .
As v ds keeps decreasing and v F increases simultaneously, the output capacitance ( C oss ) of S and C F of D is discharged and charged, respectively. Since the voltage of C oss and C F are clamped to V dd , they share the same absolute value of voltage change. The resultant capacitive displacement current for C oss ( i oss ) can be expressed by (9). Additionally, this i oss along with the counterpart for C F ( i CF = C F · d v ds / d t ) will affect i d as can be seen in Figure 3. By applying Kirchhoff’s law, i d can be determined by (10) and i S is modified accordingly to consider these displacement currents.
i oss = i gd + i ds = C oss · d v ds / d t
i d = i S = I L + i re C F · d v ds / d t
Based on Figure 1b along with the above equations, i ch can be further obtained,
i ch = i d i oss = I L + i re ( C oss + C F ) · d v ds / d t
Consequently, during this period, i d includes I L , i re and i CF , while the additional i oss is further added to i ch as shown in Figure 3. As a result, v mil will also change according to (2). This period ends when v ds drop to v on at t 3 . Thereafter, v gs will continue climbing until reaches V gon .
Furthermore, the turn-on waveforms for the case of SBD or eGaN HEMT are presented in Figure 3b. Since the reverse recovery behavior is neglected for these cases as mentioned previously, i re keeps zero and the voltage falling period starts right after i S reaches I L . Apart from that, the turn-on modelling and analysis are the same as the case of the PIN diode.
  • Turn-off transition ( t 4 t 7 )
The turn-off process can be considered as the opposite of turn-on transition and the typical transient waveforms are illustrated in Figure 5. In order to turn off S , V gon is replaced by a negative gate drive signal V goff and thus C iss is discharged through R G resulting the reduction of v gs . As v gs drops to v mil , C gd absorbs nearly all the i g and thus v ds begins to rise which again causes a current decline of i d . When v ds reaches V dd , the miller plateau disappears and i d begins decreasing with v gs which results an additional v L on v ds as shown in Figure 5a.
As v gs drops below v th , i d becomes zero and S turns off completely. However, for the case of Si IGBT, the tail current ( i tail ) is considered due to the recombination of the excess carriers. This i tail will prolong the turn off time and can be modelled by the exponential function [42],
i tail = I tail 0 · exp ( ( t t tail 0 ) / τ tail )
where τ tail stands for carrier transit time and the tailing period starts at t tail 0 with the initial value of I tail 0 . These parameters can be estimated from the turn-off current waveform.
As can be seen from Figure 5b, a notable difference for the case of eGaN HEMT is v gs does not typically plateau due to the much smaller capacitance and it keeps decreasing until reaches V goff . As a result, i ch quickly declines synchronized with v gs based on (7), meanwhile v ds rises slightly and the slew rate is limited by the relatively high C oss at low v ds . In fact, the channel turns off completely before v ds is significantly rising. However, i d does not follow the fast decreasing i ch since its changing rate is limited by the inductances in the power loop. Additionally, C oss is charged by the current difference between i d and i ch resulting a slight increase of v ds which can be expressed as
d v ds / d t = d v S / d t = ( i d i ch ) / C oss
Meanwhile, C F is discharged resulting in a reduction of v F . Additionally, i d thus can be obtained by
i d = i S = I L C F · d v F / d t
Once v gs drops below v th , namely i ch becomes zero, the channel shuts down and I L is shared by C oss and C F . During this period, v ds keeps rising according to (14). When v ds rises to V dd , I L starts commuting to D and the S turn-off transition finishes. Note that if very high R G is used for eGaN HEMT, the turn-off analysis will be the same as the case of MOSFET as shown in Figure 5a.
Based on the above analysis v S can be considered as an open circuit except for voltage rising/falling periods and S on-state. During the voltage transition period, v s is modelled as a dependent voltage source with a voltage slew rate as mentioned previously. In addition, the key expressions of i S for different conditions can be summarized in Table 1. It is noted that when v gs is less than v th , the conduction channel is not established and theoretically, no current is flowing through the device. As a result, i S is modelled with zero ampere under this condition in PSCAD/EMTDC which can be considered as open circuit. In this paper, all analytical equations for v s and i s are implemented and programmed with conditions in the custom programming modules in PSCAD. In this way, it is feasible and convenient to make any modifications as necessary.

3. Power Loss Analysis and Parameter Extraction

In general, the power losses of SDC mainly include conduction loss and switching loss. Typically the conduction losses of S and D can be calculated directly as the product of operating current (i.e., I L ) and the on-state voltage drop based on (3–5). In addition, the reverse recovery loss of D can be estimated based on the reverse recovery charge ( Q rr ) and v F from the device datasheet and the switching loss of S is analyzed as follows.

3.1. Turn on Loss( E o n )

The instantaneous power of S ( p S ) along with E on are presented in Figure 3. Basically, E on consists of the turn-on V-I overlap loss ( E vion ), the reverse recovery related loss ( E rr ) and the capacitive losses ( E oss and E F ) for C oss and C F , respectively. E vion graphically can be divided into two parts, the i d rising period and v ds decreasing period. Hence, E vion can be expressed as
E vion = t 1 t 2 v ds · i d d t + I L · t 2 t 3 v ds d
Since reverse recovery behavior of D and the displacement current of C oss already have been considered in the modelling of i S , therefore the sum of E vion , E rr and E oss , which is actually the measured turn-on loss ( E onm ), can be directly obtained by integrating i d and v ds . This significantly reduces the complexity comparing with the analytical loss model by computing the switching time for each sub-stages. Moreover, according to (11), both the discharging current of C oss and charging current of C F are flowing through the channel of S and thus these capacitive energy losses (i.e., E oss , E F ) are dissipated into the channel. Based on the capacitance curves, E oss can be expressed as
E oss = 0 V dd v ds · C oss d v ds
Since the charging current of C F is provided by V dd and part of the energy is stored in C F , thus the energy loss dissipated in the channel (i.e., E F ) can be derived based on the charge of C F ( Q F ),
E F = V dd · Q F 0 V dd v F · C F d v F = 0 V dd ( V dd v F ) · C F d v F
Consequently, E F theoretically should also be included in E on which can be expressed as,
E on = E onm + E F = t 1 t 3 i d · v ds d t + E F

3.2. Turn off Loss( E o f f )

Generally, the power loss during turn-off transition occurs from t 5 to t 7 which includes the turn-off V-I overlap loss ( E vioff ), E oss and the tailing loss for the case of IGBT ( E tail ). The analysis of E vioff is significantly different for the slow-switching scenario in Figure 5a and the typical fast-switching for eGaN HEMT in Figure 5b. As for the former case, v gs is fixed at v mil and thus i d is relatively constant throughout the voltage rising period. Thereafter, i d decreases significantly meanwhile v ds keeps relatively constant. Hence, E vioff for this case can be graphically calculated as
E vioff = t 5 t 6 v ds · i d d t + t 6 t 7 ( V dd + v L ) · i d d
Similarly, the sum of E vioff and E tail , namely the measured turn-off loss ( E offm ), is typically an integral of i d and v ds . It is noticed that, C F is discharged and the energy is transferred to the inductive load during the voltage rising period resulting in a reduction of i d , while C oss is charged and the corresponding energy (i.e., E oss ) is stored which will be dissipated in the next turn-on transition. Therefore, E oss should be theoretically excluded from E off which can be expressed as
E off = E offm E oss = t 5 t 7 v ds · i d d t E oss
As for the typical eGaN HEMT scenario, v gs skips the plateau and the channel turns off quickly before v ds rises significantly as discussed previously. Afterwards, the energy is commutating between the inductive load and the two capacitances (i.e., C F and C oss ) which is almost lossless. Since the resistive overlap loss only occurs as long as the channel is on, it is significantly reduced for this case due to the relatively low v ds during this time. Nevertheless, E off still can be calculated by (20).

3.3. Parameter Extraction

The key model parameters can be directly extracted from the corresponding curves provided in the device datasheet by the curve fitting method [46,47,48], to avoid the supplementary experiments which are usually not practical. This method is adopted in this paper since it is applicable to different semiconductor devices and provides a relatively balanced tradeoff between accuracy and practicability. As an example, different types of semiconductors from different manufacturers as listed in Table 2 are selected for modelling and validation of the proposed method. The extracting sequence is discussed in detail as follows.
  • Static characteristic
In order to reproduce the switching behavior of S , two key parameters (i.e., v th and g fs ) are considered first. Since v th typically is a temperature-dependent parameter rather than a constant value, it can be fitted by the 2nd order polynomial of junction temperature ( T j ) from the corresponding curve in the datasheet. Likewise, the transfer characteristic of S can be fitted by the quadratic function of v gs and subsequently g fs can be further determined based on v th and (7) as
g fs = i s · k ga / ( i s k gb ) · k gTa · ( T j / T a ) k gTb
where k ga , k gb , k gTa , and k gTb are fitting constants. T a is room temperature which is considered as 25   ° C here. In this way, g fs under the given T j in datasheet can be obtained. The fitted results of transfer characteristics are compared and shown in good agreement with the datasheet in Figure 6. It is also found that there is a positive correlation between T j and g fs for Si IGBT and Si/SiC MOSFET while it shows a negative correlation for the case of eGaN HEMT.
Furthermore, in order to represent the on-state characteristic of S as expressed in (3), the parameter v cesat for IGBT or R ds ( on ) for other cases is needed to be extracted. Typically, both of v cesat and R ds ( on ) are affected by T j and i S according to the curves in the datasheet. Therefore, v cesat can be obtained by the following Equation (22).
v cesat = ( k cea + r cea · i S ) + ( k ceb + r ceb · i S ) · ( T j T a )
where k cea , r cea , k ceb and r ceb are fitting coefficients. Note that the gate voltage is assumed as constant in the parameter extraction for simplicity.
Likewise, R ds ( on ) for the cases of MOSFET and eGaN HEMT can be extracted by
R ds ( on ) = ( k ona + r ona · i S ) + ( k onb + r onb · i S ) · ( T j T a )
where k ona , r ona , k onb and r onb are fitting coefficients. Based on the above equations, the key parameters for different semiconductors in this paper can be extracted and illustrated in Table 3.
  • Parasitic capacitance and inductance
It is a fact that nonlinear capacitances are the key to the dynamic characteristic of the device. Typically, the capacitance curves provided in the datasheet is in the form of C iss ,   C oss and reverse capacitance ( C rss ) which can be mathematically converted to junction capacitances. Generally, these capacitances are voltage-dependent and can be extracted by fitting the curves as
C ( v ) = f ( v ) ,
where f is a general fitting function for extraction of capacitance. In this paper, various f are used for different devices to fit the corresponding curves as shown in Figure 7. Notice that the nonlinear capacitance curves vary from different devices and thus it is reasonable to change f accordingly.
As for the parasitic inductance, only the internal inductance of the device is normally provided in the datasheet, while the stray parasitic inductance is highly related to the specific device package and PCB design. There are two widely used methods for inductance extraction, namely calculation method and experimental method. Based on the PCB and device package specification, the corresponding inductance can be computationally obtained with the help of calculation tools such as the Ansys Q3D Extractor software (v1.0, Canonsburg, PA, USA). According to (1), it also can be extracted from the slew rate of current along with v L during S turn on transition or from the resonant frequency of the power loop in the experimental results. In this paper, the parasitic inductances are initially estimated based on the PCB trace length of the power loop and the gate loop [50] as well as the device package (e.g., 2–5 nH for TO-247 [51]) and further calibrated from the switching waveforms.
  • Diode parameters
According to (4), R F and v F 0 are the key static parameters for D which can be extracted directly from the diode I-V curve in the datasheet and the corresponding values for various temperatures can be estimated by linear interpolation. Furthermore, the third quadrant operation of eGaN HEMT as D is of special concern. Since the corresponding voltage drop is dependent on the gate drive voltage of D , thus it should be fitted by (5) based on the output curves in the third quadrant from the datasheet. The diode I-V curve fitted results for different devices are compared with the datasheet and illustrated in Figure 8. Notice that the conduction performance of the body diode in SiC MOSFET is generally worse than the anti-parallel diode of IGBT and SiC SBD. As for the eGaN HEMT, the reverse voltage drop is highly dependent on the gate drive voltage and the typical value for turn off (i.e., −3 V) will result in considerable conduction loss.
In addition, C F can be obtained using the same method as mentioned above for S from the capacitance curve in the datasheet. According to the previous switching transition analysis, it can be found that the reverse recovery behavior of diode plays a considerable role and the main parameters I rm and Q rr can be extracted from the diode curves as a function of T j and d i F / d t [42],
I rm = k rm 0 + t rm 0 · d i F / d t + ( k rm 1 + t rm 1 · d i F / dt ) · ( T j T a )
Q rr = k rr 0 + t rr 0 · d i F / d t + ( k rr 1 + t rr 1 · d i F / dt ) · ( T j T a )
where k rm 0 ,   t rm 0 , k rm 1 , t rm 1 , k rr 0 , t rr 0 , k rr 1 and t rr 1 are all fitting coefficients. Besides, the reverse recovery time ( t rr ) and τ re can be further determined by
t rr = 2 Q rr / ( d i F / d t )
τ re = 1 / ln 10 · ( t rr I rm / ( d i F / d t ) )

4. Experimental Verification

The objective of the proposed method is to reproduce the switching performance of the SDC and generate the corresponding power loss look-up table (LUT) with reasonable accuracy and fast simulation speed. The models are implemented in PSCAD/EMTDC and validated by comparing with the experimental waveforms and power loss results in the DPT bench for different semiconductor combinations.

4.1. Setup Description

An automatic DPT bench was designed and built for device characterization and loss validation [17]. Since the gate-drive requirements and device package are different for all the devices, three daughter boards were designed accordingly with a power supply (up to 1 kV) and the inductive load (5 mH) as shown in Figure 9. Tektronix High voltage differential probe THDP0200 and current probe TCP0030A were used for Si IGBT/MOSFET voltage and current measurements, respectively, while passive voltage probe (TTP800) and 0.1 Ω current shunt resistor (SDN-414-01) were adopted for SiC MOSFET and eGaN HEMT measurements. The temperature was controlled by a heating block and monitored by a thermal imager (Fluke, TiS40). In the DPT, the desired test conditions (voltage, current, and temperature) were initially set in the personal computer (PC) and all the control signals were given by the microcontroller Arduino on the board. Afterwards, the DC capacitor bank was charged to the desired voltage by the power supply unit, and the device was heated to the desired temperature. When voltage and temperature conditions were ready, two gate pulses were given in sequence to turn on and turn off the device under test (DUT). The switching waveforms and data were obtained by oscilloscope and processed in the PC for transient information and power loss analysis. In order to mitigate the measurement error of power loss due to the asynchrony of voltage and current, it is necessary to calibrate the probes before conducting the DPTs. To further guarantee the accuracy, additional delay time adjustments for the transient waveforms are also needed for the turn-on and turn-off processes. Taking the current as the reference, the calibration time of probe and waveform are provided in Table 4. The calibrating fixture (067-1686-02) from Tektronix was used for calibration of the current and voltage probes. A 10 MHz sinusoidal signal was applied to both probes and the deskew time for voltage probe was adjusted until both measurements were synchronizing. Note that, these calibration times can be different for various probes and DUTs.

4.2. Switching Transient Verification

  • Si IGBT
The daughterboard in Figure 9a was used for Si IGBT and MOSFET tests, and v G was flipped between 15 V and 0 V to control the DUT’s on and off, respectively. The simulated results of switching current and voltage waveforms for IKW40T120 are compared with the DPT measurements in Figure 10. The simulated results demonstrate good agreement with measurements for current and voltage switching waveforms under 25   ° C and 150   ° C . The switching details such as the tail current and the current spike resulting from the reverse recovery of D can be clearly observed. In addition, v ce slightly drops to 500 V as current rising and reaches a peak of 700 V during turn-off transition due to parasitic inductance. Besides, as T j increases from 25   ° C to 150   ° C , the reverse recovery behavior of D becomes more obvious resulting higher current peak (up to 60 A) and the rise of v ce as well as the decline of i c during turn off transition slows down which will increase the switching power loss.
  • SiC MOSFET
A more compact daughterboard as shown in Figure 9b was designed for testing SiC MOSFET. Additionally, the gate drive integrated circuit (IXDN609SI) was adopted as the gate driver onboard to provide 20 V/−5 V drive voltage for SiC MOSFET. Figure 11 shows the simulated switching waveforms of SiC MOSFET under the condition of 600 V and 20 A which match well with measured results. It can be observed that there is only a slight impact of T j on the switching transients in terms of turn-on and turn-off time. Nevertheless, the current still can reach almost 40 A during the turn-on transition due to the reverse recovery behavior of the diode. Besides, it is found that a current ringing occurs during both turn-on and turn-off transitions because of the parasitic resonance. This ringing energy is generally consumed by the HF damping resistance in the circuit. Since either the voltage or current has typically dropped to a low level during the ringing period, thus this ringing loss is neglected in the model for simplicity.
  • Si MOSFET with SiC diode
The DPT results for Si MOSFET with SiC diode using the same daughterboard as for testing Si IGBT are shown in Figure 12. In general, the simulated results match well with the measured results for different R G conditions. As R G increases from 10   Ω to 33   Ω , a half less voltage drop of v ds can be observed during turn-on transition due to the slower current rising speed. Likewise, only a slight increase of v ds can be seen after v ds climbs to V dd . Moreover, it is noted that the current spike is significantly limited comparing with the previous testing using PIN diode due to the merit of zero reverse recovery for SiC SBD. Whereas, there is still a slight current bump causing by the resonance of parasitics as well as the capacitive displacement current as discussed previously.
  • eGaN HEMT
In order to test the eGaN HEMT which is a surface-mount device (SMD), a specific DPT daughterboard was used as shown in Figure 9c. The gate driver provided 6 V/0 V as gate drive voltage to control the lower side GaN switch S , while the upper side SiC SBD served as a freewheeling diode when S turned off. The simulated switching results are compared with the measurements for the two operating conditions as shown in Figure 13. It can be seen that the simulation results are consistent with the experimental results. In the turn-on waveforms, the current rising time is only tens of nanoseconds. After i d reaches I L , it behaves in the similar manner as previous test using Si MOSFET with SiC SBD. However, during the turn-off period, it can be clearly observed that, i d declines significantly and drops to zero almost the same time as v ds reaches steady state while for the other cases of devices, the fast decrease of current typically occurs after v ds climbs to v dd . This is mainly because the channel of eGaN shuts down very fast before v ds increase significantly as discussed in Section 2.3. Hence, when the channel turns off completely, the apparent i d is dominated by the capacitive displacement current which is highly related to the change of v ds .

4.3. Power Loss Verification

With the aim of power loss verification for various devices, the switching losses were measured in the DPT and compared with the simulated results. During the switching transients, p S , which is the product of voltage and current, can be obtained using the math function in oscilloscope and similarly the E onm and E offm can also be obtained by integrating p S . As mentioned previously, the current and voltage probes are calibrated for each test and additional delay time is also added to the waveforms results to keep transient voltage and current synchronous. The captured waveforms and simulated waveforms under the same test conditions are demonstrated in Figure 14 taking Si IGBT as an example. By comparing the measured results with the simulated results, a good agreement can be clearly seen in terms of not only transient voltage and current waveforms but also computed p S , E onm and E offm . Furthermore, the measured power loss results are compared with the simulated loss results for different devices under various operating conditions to validate the proposed method. The average error (ē) is calculated by averaging the absolute value of the error in each case.
Figure 15 shows the power loss results of Si IGBT under different conditions of current, voltage, and temperature. Generally, the total power loss ( E ts ) increases as the operating voltage and current increase, and E off is less than E on except for the high temperature condition. It can be seen that, the average errors of E ts are within 7%, namely 5.1%, 5.5%, and 6.3% for different operating conditions of current, voltage, temperature, respectively.
Likewise, the power loss comparison results for SiC MOSFET are illustrated in Figure 16. Note that the E ts for SiC MOSFET is typically less than 1 mJ which is much less than the counterpart of Si IGBT for similar conditions. It also can be found in Figure 16c that, there is a negative correlation between E on and T j . Since E on is the dominated loss, as T j rises, E ts reduces accordingly though E off increases slightly. It is also noted that ē of E off for various conditions are more than 7% while ē for E on and E ts are still within acceptable range. The reasons for the loss deviation can be the underestimation of parasitics and ringing loss as well as measurement error. Besides, a relatively small amount of loss deviation can still result in a high error percentage when the overall loss is relatively low.
Figure 17 shows E ts results for the combination of Si MOSFET and SiC SBD as the SDC. It can be observed that E ts increases as the operating voltage and current increase. In addition, a slight increase of E ts can be found as the operating temperature rise from 25 °C to 150 °C, while E ts increases significantly when 33 Ω R G is used. Moreover, the power loss results for the case of eGaN HEMT with SiC SBD are illustrated in Figure 18. It should be mentioned that, in order to capture the switching waveforms and the power loss with reasonable accuracy, a 220 Ω gate resistance is used to relatively sacrifice the switching speed and avoid shoot through issue due to the very low v th of GaN. Comparing with the simulated and measured results, a good agreement is achieved at various testing conditions and the error is within an acceptable range, although the ē of E off is slightly higher.

4.4. Discussion

  • Accuracy
Comparing with the original two-state resistance switch model in PSCAD/EMTDC, the proposed model is capable to reproduce the switching transient waveforms considering various impacts of parasitics and interactive behavior of diode. In addition, the thermal effect is also considered to provide reasonably accurate results comparing with the measured results where the temperature is monitored by a thermal imager as shown in Figure 19a. Apart from the switching waveforms, multi-dimensional (i.e., voltage, current, temperature) power loss LUT as shown in Figure 19b also can be obtained simultaneously. The average error is within 10% comparing with measured results for various devices under various conditions. Comparing to the traditional physical model or analytical loss model, no significant advantage is found in the modelling accuracy using the proposed model due to the ignorance of parasitic resonance and some linear assumptions. Nevertheless, the complexity of the proposed model is reduced with no state equations and numerical calculations, and all the model parameters can be extracted from the datasheet.
  • Efficiency
The proposed model uses equivalent dependent voltage and current sources to represent the dynamic characteristics of devices based on the analytical equations for each sub-stage of the switching process. In addition, the gate loop and power loop are decoupled and the complicated numerical calculation as well as solving physical equations are not necessary which can boost the simulation efficiency and avoid convergence issue. In order to obtain an accurate power loss LUT with a wide range of operating conditions, there are numerous permutations to be taken into account and thousands of simulation runs are required instead of repetitive DPTs. For example, it requires around 5400 simulation runs to cover the operating range, namely voltage from 20 V to 600 V with 20 V step, current from 2 A to 60 A with a 2 A step, and temperature from 25 °C to 150 °C with a 25 °C step. Figure 19c demonstrates the time cost of using the SPICE model provided by the manufacturer and using the proposed model in PSCAD. In order to achieve reasonable accuracy, the simulation time step is typically one nanosecond or less. Notice that, it takes thousands of seconds to barely finish around 200 runs in SPICE, while more than 10 times less running time is needed to finish the same number of runs by the proposed model with even less time-step (i.e., 0.1 ns). Furthermore, to finish the whole 5400 simulation runs using the proposed model at 1 ns time-step, it takes less than 300 s which shows the merit of time-saving in generating the power loss LUT.
  • Applicability
The proposed approach can be used to reproduce the switching waveforms and obtain the power loss LUT of SDC configured by various devices such as Si IGBT, Si/SiC MOSFET, eGaN HEMT, and SiC SBD. When it comes to other devices with a new structure such as Cascade GaN, the proposed model cannot be used directly and modifications are needed though the basic modelling method is still applicable. In addition, for each specific device, the curve fitting functions and algorithms for parameter extraction are needed to be adjusted for good fitting results. Apart from PSCAD/EMTDC, the proposed approach can also be applied to other simulators such as MATLAB/Simulink, PLECS, and Saber with respective modifications. Besides, the proposed model provides an insight into the device behavior and most of the elements have clear descriptions. As a result, it is more easily apprehensible than mathematical equations. Conversely, analytical loss models are normally limited to the specific device type or combination and it is difficult to extend the models to various PE applications for loss estimation. Besides, loss measurements are time-consuming, costly, and challenging especially for WBG devices due to the fast switching speed.

5. Conclusions

In this paper, a generalized behavioral modelling approach of the switch-diode cell for power loss prediction is proposed, implemented in PSCAD/EMTDC, and validated by experimental results in double-pulse tests. This proposed model consists of an active switch model and diode model and it can be used for different modern power semiconductors. The modelling approach along with power loss analysis is derived based on the comprehensive switching process analysis in a clamped inductive switching circuit. The static and dynamic characteristics of the switch-diode cell are modelled by dependent voltage and current sources with passive components. In addition, the proposed model is improved by considering the impacts of parasitic elements, interactive behavior of diode, and the temperature-dependent parameters. Besides, the extraction of the model parameters is introduced by curve fitting from the device datasheet. Moreover, the switching transient verification along with power verification is conducted for different devices under a wide range of operating conditions. A good agreement between the simulated results using the proposed model and experimental results can be achieved with less than 10% average error. Consequently, the proposed model provides a good balance in terms of accuracy, efficiency, and applicability.

Author Contributions

Conceptualization, Y.X. and C.N.M.H.; methodology, Y.X., C.N.M.H. and A.G.; software, C.N.M.H. and D.M.; validation, Y.X. and A.G.; formal analysis, Y.X.; investigation, Y.X.; resources, A.G., C.N.M.H. and D.M.; writing—original draft preparation, Y.X.; writing—review and editing, C.N.M.H. and A.G.; visualization, Y.X.; supervision, C.N.M.H.; project administration, C.N.M.H.; funding acquisition, D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the NSERC Collaborative Research and Development Grants, Canada and in part by the Manitoba Hydro International, Canada.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baliga, B.J. Fundamentals of Power Semiconductor Devices; Springer: New York, NY, USA, 2008. [Google Scholar]
  2. Iwamuro, N.; Laska, T. IGBT History, State-of-the-Art, and Future Prospects. IEEE Trans. Electron. Devices 2017, 64, 741–752. [Google Scholar] [CrossRef]
  3. Grant, D.A.; Gowar, J. Power MOSFET Theory and Applications; Wiley: New York, NY, USA, 1989. [Google Scholar]
  4. Islam, M.; Mekhilef, S. Efficient Transformerless MOSFET Inverter for a Grid-Tied Photovoltaic System. IEEE Trans. Power Electron. 2016, 31, 6305–6316. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Yao, C.; Zhang, X.; Chen, H.; Li, H.; Wang, J. Power Loss Model for GaN-Based MHz Critical Conduction Mode Power Factor Correction Circuits. IEEE Trans. Emerg. Sel. Top. Power Electron. 2020, 8, 141–151. [Google Scholar] [CrossRef]
  6. Yang, B.; Wang, J.; Xu, S.; Korec, J.; Shen, Z.J. Advanced Low-Voltage Power MOSFET Technology for Power Supply in Package Applications. IEEE Trans. Power Electron. 2013, 28, 4202–4215. [Google Scholar] [CrossRef]
  7. Marroqui, D.; Garrigós, A.; Torres, C.; Orts, C.; Blanes, J.M.; Gutierrez, R. Interleaved, Switched Inductor and High-Gain Wide Bandgap Based Boost Converter Proposal. Energies 2021, 14, 800. [Google Scholar] [CrossRef]
  8. Waradzyn, Z.; Stala, R.; Mondzik, A.; Skała, A.; Penczek, A. GaN-Based DC-DC Resonant Boost Converter with Very High Eciency and Voltage Gain Control. Energies 2020, 13, 6403. [Google Scholar] [CrossRef]
  9. Millán, J.; Godignon, P.; Perpiñà, X.; Pérez-Tomás, A.; Rebollo, J. A Survey of Wide Bandgap Power Semiconductor Devices. IEEE Trans. Power Electron. 2014, 29, 2155–2163. [Google Scholar] [CrossRef]
  10. Cha, K.H.; Ju, C.T.; Kim, R.Y. Analysis and Evaluation of WBG Power Device in High Frequency Induction Heating Application. Energies 2020, 13, 5351. [Google Scholar] [CrossRef]
  11. Li, X.B.; Ma, H.B.; Yi, J.H.; Lu, S.; Xu, J.P. A Comparative Study of GaN HEMT and Si MOSFET-Based Active Clamp Forward Converters. Energies 2020, 13, 4160. [Google Scholar] [CrossRef]
  12. Cooper, J.A.; Melloch, M.R.; Singh, R.; Agarwal, A.; Palmour, J.W. Status and prospects for SiC power MOSFETs. IEEE Trans. Electron. Devices 2002, 49, 658–664. [Google Scholar] [CrossRef]
  13. Lidow, A.; Rooij, M.d.; Strydom, J.; Reusch, D.; Glaser, J. GaN Transistors for Efficient Power Conversion, 3rd ed.; Wiley: Hoboken, NJ, USA, 2020; pp. 1–23. [Google Scholar]
  14. Jones, E.A.; Wang, F.F.; Costinett, D. Review of Commercial GaN Power Devices and GaN-Based Converter Design Challenges. IEEE Trans. Emerg. Sel. Top. Power Electron. 2016, 4, 707–719. [Google Scholar] [CrossRef]
  15. Bryant, A.T.; Mawby, P.A.; Palmer, P.R.; Santi, E.; Hudgins, J.L. Exploration of power device reliability using compact device models and fast electrothermal simulation. IEEE Trans. Ind. Appl. 2008, 44, 894–903. [Google Scholar] [CrossRef]
  16. Gole, A.M.; Keri, A.; Kwankpa, C.; Gunther, E.W.; Dommel, H.W.; Hassan, I.; Marti, J.R.; Martinez, J.A.; Fehrle, K.G.; Tang, L.; et al. Guidelines for modelling power electronics in electric power engineering applications. IEEE Trans. Power Del. 1997, 12, 505–514. [Google Scholar] [CrossRef]
  17. Ghosh, A.; Ho, C.; Prendergast, J.; Xu, Y. Conceptual Design and Demonstration of an Automatic System for Extracting Switching Loss and Creating Data Library of Power Semiconductors. IEEE Open J. Power Electron. 2020, 1, 431–444. [Google Scholar] [CrossRef]
  18. Mantooth, H.A.; Peng, K.; Santi, E.; Hudgins, J.L. Modelling of Wide Bandgap Power Semiconductor Devices—Part, I. IEEE Trans. Electron. Devices 2015, 62, 423–433. [Google Scholar] [CrossRef]
  19. Kraus, R.; Castellazzi, A. A Physics-Based Compact Model of SiC Power MOSFETs. IEEE Trans. Power Electron. 2016, 31, 5863–5870. [Google Scholar] [CrossRef]
  20. McNutt, T.R.; Hefner, A.R.; Mantooth, H.A.; Berning, D.; Ryu, S. Silicon Carbide Power MOSFET Model and Parameter Extraction Sequence. IEEE Trans. Power Electron. 2007, 22, 353–363. [Google Scholar] [CrossRef]
  21. Mudholkar, M.; Ahmed, S.; Ericson, M.N.; Frank, S.S.; Britton, C.L.; Mantooth, H.A. Datasheet Driven Silicon Carbide Power MOSFET Model. IEEE Trans. Power Electron. 2014, 29, 2220–2228. [Google Scholar] [CrossRef]
  22. Li, H.; Zhao, X.; Su, W.; Sun, K.; You, X.; Zheng, T.Q. Nonsegmented PSpice Circuit Model of GaN HEMT with Simulation Convergence Consideration. IEEE Trans. Ind. Electron. 2017, 64, 8992–9000. [Google Scholar] [CrossRef]
  23. Li, H.; Zhao, X.; Sun, K.; Zhao, Z.; Cao, G.; Zheng, T.Q. A Non-Segmented PSpice Model of SiC mosfet With Temperature-Dependent Parameters. IEEE Trans. Power Electron. 2019, 34, 4603–4612. [Google Scholar] [CrossRef]
  24. Turzynski, M.; Kulesza, W.J. A Simplified Behavioral MOSFET Model Based on Parameters Extraction for Circuit Simulations. IEEE Trans. Power Electron. 2016, 31, 3096–3105. [Google Scholar] [CrossRef]
  25. Wen, Y.; Yang, Y.; Gao, Y. Model parameter calibration method of SiC power MOSFETs behavioural model. IET Power Electron. 2020, 13, 426–435. [Google Scholar] [CrossRef]
  26. Ren, Y.; Xu, M.; Zhou, J.; Lee, F.C. Analytical loss model of power MOSFET. IEEE Trans. Power Electron. 2006, 21, 310–319. [Google Scholar]
  27. Rąbkowski, J.; Skoneczny, H.; Kopacz, R.; Trochimiuk, P.; Wrona, G. A Simple Method to Validate Power Loss in Medium Voltage SiC MOSFETs and Schottky Diodes Operating in a Three-Phase Inverter. Energies 2020, 13, 4773. [Google Scholar] [CrossRef]
  28. Mukunoki, Y.; Konno, K.; Tsubasa, M.; Horiguchi, T.; Nishizawa, A.; Kuzumoto, M.; Hagiwara, M.; Akagi, H. An Improved Compact Model for a Silicon-Carbide MOSFET, and its Application to Accurate Circuit Simulation. IEEE Trans. Power Electron. 2018, 33, 9834–9842. [Google Scholar] [CrossRef]
  29. Jin, M.; Gao, Q.; Wang, Y.; Xu, D. A Temperature-Dependent SiC MOSFET Modelling Method Based on MATLAB/Simulink. IEEE Access 2017, 6, 4497–4505. [Google Scholar] [CrossRef]
  30. Riccio, M.; D’Alessandro, V.; Romano, G.; Maresca, L.; Breglio, G.; Irace, A. A Temperature-Dependent SPICE Model of SiC Power MOSFETs for Within and Out-of-SOA Simulations. IEEE Trans. Power Electron. 2018, 33, 8020–8029. [Google Scholar] [CrossRef]
  31. Shi, B.; Zhao, Z.; Zhu, Y. Piecewise Analytical Transient Model for Power Switching Device Commutation Unit. IEEE Trans. Power Electron. 2019, 34, 5720–5736. [Google Scholar] [CrossRef]
  32. Yuan, D.; Zhang, Y.; Wang, X. An Improved Analytical Model for Crosstalk of SiC MOSFET in a Bridge-Arm Configuration. Energies 2021, 14, 683. [Google Scholar] [CrossRef]
  33. Wang, X.; Zhao, Z.; Li, K.; Zhu, Y.; Chen, K. Analytical Methodology for Loss Calculation of SiC MOSFETs. IEEE Trans. Emerg. Sel. Top. Power Electron. 2019, 7, 71–83. [Google Scholar] [CrossRef]
  34. Xie, R.; Wang, H.; Tang, G.; Yang, X.; Chen, K.J. An Analytical Model for False Turn-On Evaluation of High-Voltage Enhancement-Mode GaN Transistor in Bridge-Leg Configuration. IEEE Trans. Power Electron. 2017, 32, 6416–6433. [Google Scholar] [CrossRef]
  35. Agrawal, B.; Freindl, M.; Bilgin, B.; Emadi, A. Estimating switching losses for SiC MOSFETs with non-flat miller plateau region. In Proceedings of the 2017 IEEE Applied Power Electronics Conference and Exposition (APEC), Tampa, FL, USA, 26–30 March 2017; pp. 2664–2670. [Google Scholar]
  36. Wang, K.; Yang, X.; Li, H.; Ma, H.; Zeng, X.; Chen, W. An Analytical Switching Process Model of Low-Voltage eGaN HEMTs for Loss Calculation. IEEE Trans. Power Electron. 2016, 31, 635–647. [Google Scholar] [CrossRef]
  37. Chen, J.; Luo, Q.; Huang, J.; He, Q.; Du, X. A Complete Switching Analytical Model of Low-Voltage eGaN HEMTs and Its Application in Loss Analysis. IEEE Trans. Ind. Electron. 2020, 67, 1615–1625. [Google Scholar] [CrossRef]
  38. Gareau, J.; Hou, R.; Emadi, A. Review of Loss Distribution, Analysis, and Measurement Techniques for GaN HEMTs. IEEE Trans. Power Electron. 2020, 35, 7405–7418. [Google Scholar] [CrossRef]
  39. Hou, R.; Lu, J.; Chen, D. Parasitic capacitance Eqoss loss mechanism, calculation, and measurement in hard-switching for GaN HEMTs. In Proceedings of the 2018 IEEE Applied Power Electronics Conference and Exposition (APEC), San Antonio, TX, USA, 4–8 March 2018; pp. 919–924. [Google Scholar]
  40. Ho, C.; Breuninger, S.; Pettersson, G.; Escobar, L. Serpa, and A. Coccia. Practical design and implementation procedure of an interleaved boost converter using SiC diodes for PV applications. IEEE Trans. Power Electron. 2012, 27, 2834–2845. [Google Scholar] [CrossRef]
  41. Ho, C.; Breuninger, H.; Pettersson, S.; Escobar, G.; Canales, F. A comparative performance study of an interleaved boost converter using commercial Si and SiC diodes for PV applications. IEEE Trans. Power Electron. 2013, 28, 289–299. [Google Scholar] [CrossRef]
  42. Xu, Y.; Ho, C.N.M.; Ghosh, A.; Muthumuni, D. An Electrical Transient Model of IGBT-Diode Switching Cell for Power Semiconductor Loss Estimation in Electromagnetic Transient Simulation. IEEE Trans. Power Electron. 2020, 35, 2979–2989. [Google Scholar] [CrossRef]
  43. Xu, Y.; Ho, C.N.M.; Ghosh, A.; Muthumuni, D. Design, Implementation and Validation of Electro-Thermal Simulation for SiC MOSFETs in Power Electronic Systems. IEEE Trans. Ind. Appl. 2021. [Google Scholar] [CrossRef]
  44. Reiner, R.; Waltereit, P.; Weiss, B.; Quay, R.; Ambacher, O. Investigation of GaN-HEMTs in Reverse Conduction. In Proceedings of the PCIM Asia, International Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management, Nuremberg, Germany, 16–18 March 2017; pp. 1–8. [Google Scholar]
  45. Lauritzen, P.O.; Ma, C.L. A Simple Diode Model with Reverse Recovery. IEEE Trans. Power Electron. 1991, 6, 188–191. [Google Scholar] [CrossRef]
  46. Jarndal, A.; Essaadali, R.; Kouki, A.B. A Reliable Model Parameter Extraction Method Applied to AlGaN/GaN HEMTs. IEEE Trans. Comput. Aided Design Integr. Circuits Syst. 2016, 35, 211–219. [Google Scholar] [CrossRef]
  47. Jouha, W.; El Oualkadi, A.; Dherbécourt, P.; Joubert, E. Silicon Carbide Power MOSFET Model: An Accurate Parameter Extraction Method Based on Levenberg-Marquardt Algorithm. IEEE Trans. Power Electron. 2018, 33, 9130–9133. [Google Scholar] [CrossRef]
  48. Bryant, A.T.; Kang, X.; Santi, E.; Palmer, P.R.; Hudgins, J.L. Two-step parameter extraction procedure with formal optimization for physics-based circuit simulator IGBT and p-i-n diode models. IEEE Trans. Power Electron. 2006, 21, 295–309. [Google Scholar] [CrossRef]
  49. Ahmed, M.R.; Todd, R.; Forsyth, A.J. Predicting SiC MOSFET Behavior Under Hard-Switching, Soft-Switching, and False Turn-On Conditions. IEEE Trans. Ind. Electron. 2017, 64, 9001–9011. [Google Scholar] [CrossRef]
  50. Howard, J.; Martin, G. High-Speed Digital Design: A Handbook of Black Magic; Pearson Education: London, UK, 1999. [Google Scholar]
  51. Stueckler, F.; Vecino, E. CoolMOS C7 650V Switch in a Kelvin Source Configuration; Infineon Technology. Application Note AN 2013-05; Infineon Technologies: Neubiberg, Germany, 2013. [Google Scholar]
Figure 1. Schematic diagrams of (a) switch-diode cell circuit; (b) DPT circuit; (c) Proposed switch and diode models.
Figure 1. Schematic diagrams of (a) switch-diode cell circuit; (b) DPT circuit; (c) Proposed switch and diode models.
Energies 14 01500 g001
Figure 2. Flow chart of the proposed modelling procedure.
Figure 2. Flow chart of the proposed modelling procedure.
Energies 14 01500 g002
Figure 3. Typical turn on waveforms for S with D (a) PIN or body diode; (b) SBD or eGaN HEMT.
Figure 3. Typical turn on waveforms for S with D (a) PIN or body diode; (b) SBD or eGaN HEMT.
Energies 14 01500 g003
Figure 4. Simplified equivalent circuits during S turn on transition.
Figure 4. Simplified equivalent circuits during S turn on transition.
Energies 14 01500 g004
Figure 5. Typical turn off waveforms for S as (a) Si IGBT or Si/SiC MOSFETs; (b) Typical eGaN HEMT.
Figure 5. Typical turn off waveforms for S as (a) Si IGBT or Si/SiC MOSFETs; (b) Typical eGaN HEMT.
Energies 14 01500 g005
Figure 6. Fitted results of transfer curves for (a) Si IGBT; (b) SiC MOSFET; (c) Si MOSFET; (d) eGaN HEMT.
Figure 6. Fitted results of transfer curves for (a) Si IGBT; (b) SiC MOSFET; (c) Si MOSFET; (d) eGaN HEMT.
Energies 14 01500 g006aEnergies 14 01500 g006b
Figure 7. Parasitic capacitance extraction with the methods of (a) 5th order Gauss function; (b) Analytical equation [49]; (c) Interpolating look-up table (LUT); (d) 3rd order exponential function.
Figure 7. Parasitic capacitance extraction with the methods of (a) 5th order Gauss function; (b) Analytical equation [49]; (c) Interpolating look-up table (LUT); (d) 3rd order exponential function.
Energies 14 01500 g007
Figure 8. Diode static parameter extraction of (a) Si IGBT; (b) SiC MOSFET; (c) Si MOSFET; (d) eGaN HEMT.
Figure 8. Diode static parameter extraction of (a) Si IGBT; (b) SiC MOSFET; (c) Si MOSFET; (d) eGaN HEMT.
Energies 14 01500 g008
Figure 9. Device characterization DPT setups for (a) Si IGBT/MOSFET; (b) SiC MOSFET; (c) eGaN HEMT.
Figure 9. Device characterization DPT setups for (a) Si IGBT/MOSFET; (b) SiC MOSFET; (c) eGaN HEMT.
Energies 14 01500 g009
Figure 10. Switching waveforms of Si IGBTs (a) turn-on; (b) turn-off; @ R G = 15   Ω .
Figure 10. Switching waveforms of Si IGBTs (a) turn-on; (b) turn-off; @ R G = 15   Ω .
Energies 14 01500 g010
Figure 11. Switching waveforms of SiC MOSFETs (a) turn-on; (b) turn-off; @ R G = 5   Ω .
Figure 11. Switching waveforms of SiC MOSFETs (a) turn-on; (b) turn-off; @ R G = 5   Ω .
Energies 14 01500 g011
Figure 12. Switching waveforms of Si MOSFET with SiC SBD (a) turn-on; (b) turn-off; @ 400   V , 20   A , 25 °C.
Figure 12. Switching waveforms of Si MOSFET with SiC SBD (a) turn-on; (b) turn-off; @ 400   V , 20   A , 25 °C.
Energies 14 01500 g012
Figure 13. Switching waveforms of eGaN HEMT with SiC SBD (a) turn-on; (b) turn-off; @ 25   ° C ,   R G = 5   Ω .
Figure 13. Switching waveforms of eGaN HEMT with SiC SBD (a) turn-on; (b) turn-off; @ 25   ° C ,   R G = 5   Ω .
Energies 14 01500 g013
Figure 14. Si IGBT turn on waveforms (a) measured; (b) simulated; and turn off waveforms (c) measured; (d) simulated @ 600   V , 20   A , 25   ° C ,   R G 15   Ω .
Figure 14. Si IGBT turn on waveforms (a) measured; (b) simulated; and turn off waveforms (c) measured; (d) simulated @ 600   V , 20   A , 25   ° C ,   R G 15   Ω .
Energies 14 01500 g014
Figure 15. Si IGBT loss results under different values of (a) current; (b) voltage; (c) Temperature.
Figure 15. Si IGBT loss results under different values of (a) current; (b) voltage; (c) Temperature.
Energies 14 01500 g015
Figure 16. Power losses of SiC MOSFET under different values of (a) current; (b) voltage; (c) temperature.
Figure 16. Power losses of SiC MOSFET under different values of (a) current; (b) voltage; (c) temperature.
Energies 14 01500 g016
Figure 17. Power losses of Si MOSFET under various values of (a) voltage and current, (b) temperature.
Figure 17. Power losses of Si MOSFET under various values of (a) voltage and current, (b) temperature.
Energies 14 01500 g017
Figure 18. Power loss results of eGaN HEMT with SiC SBD (a) turn on, (b) turn off.
Figure 18. Power loss results of eGaN HEMT with SiC SBD (a) turn on, (b) turn off.
Energies 14 01500 g018
Figure 19. DPT results of (a) thermal image; (b) loss LUT for Si IGBT; (c) time cost comparison.
Figure 19. DPT results of (a) thermal image; (b) loss LUT for Si IGBT; (c) time cost comparison.
Energies 14 01500 g019
Table 1. Key expressions for i S in the proposed model.
Table 1. Key expressions for i S in the proposed model.
i S Condition v g s < v t h i d < I L i d > I L Tail Period for IGBTTurn-off for eGaN HEMT
Expression0(7)(10)(12)(14)
Table 2. Semiconductor devices selected for modelling and validation.
Table 2. Semiconductor devices selected for modelling and validation.
Device TypeSi IGBTSiC MOSFETSi MOSFETeGaN HEMTSiC SBD
Part NumberIKW40T120SCT2080KENVHL072N65S3GS66506TSCS220KG
ManufacturerInfineonRohmOn SemiconductorGaN SystemROHM
Table 3. Key fitting coefficients parameters of different semiconductors.
Table 3. Key fitting coefficients parameters of different semiconductors.
Parameter k ga k gb k gTa k gTb k cea , k ona r cea , r ona k ceb , k onb r ceb , r onb
Si IGBT2.05−6.9614.2−0.012−0.001 1 · 10 4 0.9720.0215
SiC MOSFET0.42−0.352.96−0.086 4.1 · 10 4 1.12 · 10 6 0.073 5.1 · 10 4
Si MOSFET3.72−19.213−0.570.048 2.87 · 10 3 0.0340.0012
eGaN HEMT0.23−776.90.51−0.310.067 7.2 · 10 5 8.4 · 10 4 4 · 10 6
Table 4. Calibration time of probes and waveforms.
Table 4. Calibration time of probes and waveforms.
DUT CaseSi IGBTSiC MOSFETSi MOSFETeGaN HEMT
Probes calibration (ns)48264523
Turn-on calibration (ns)35132614
Turn-off calibration (ns)31182512
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, Y.; Ho, C.N.M.; Ghosh, A.; Muthumuni, D. Generalized Behavioral Modelling Methodology of Switch-Diode Cell for Power Loss Prediction in Electromagnetic Transient Simulation. Energies 2021, 14, 1500. https://doi.org/10.3390/en14051500

AMA Style

Xu Y, Ho CNM, Ghosh A, Muthumuni D. Generalized Behavioral Modelling Methodology of Switch-Diode Cell for Power Loss Prediction in Electromagnetic Transient Simulation. Energies. 2021; 14(5):1500. https://doi.org/10.3390/en14051500

Chicago/Turabian Style

Xu, Yanming, Carl Ngai Man Ho, Avishek Ghosh, and Dharshana Muthumuni. 2021. "Generalized Behavioral Modelling Methodology of Switch-Diode Cell for Power Loss Prediction in Electromagnetic Transient Simulation" Energies 14, no. 5: 1500. https://doi.org/10.3390/en14051500

APA Style

Xu, Y., Ho, C. N. M., Ghosh, A., & Muthumuni, D. (2021). Generalized Behavioral Modelling Methodology of Switch-Diode Cell for Power Loss Prediction in Electromagnetic Transient Simulation. Energies, 14(5), 1500. https://doi.org/10.3390/en14051500

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

Article Metrics

Back to TopTop