electronics

Article
Feasible Evaluation and Implementation of Shunt Active Filter
for Harmonic Mitigation in Induction Heating System
Rahul Raman 1 , Pradip Kumar Sadhu 1 , Ritesh Kumar 2 , Shriram Srinivasarangan Rangarajan 2,3, *,
Umashankar Subramaniam 4 , Edward Randolph Collins 3,5 and Tomonobu Senjyu 6, *

1 Department of Electrical Engineering, Indian Institute of Technology (ISM), Dhanbad 826004, India
2 Department of Electrical and Electronics Engineering, S R University, Hasanparthy, Warangal 506371, India
3 Department of Electrical and Computer Engineering, Clemson University, Clemson, SC 29634, USA
4 Renewable Energy Lab, College of Engineering, Prince Sultan University, Riyadh 11586, Saudi Arabia
5 College of Engineering & Technology, Western Carolina University, Cullowhee, NC 28723, USA
6 Department of Electrical and Electronics Engineering, Faculty of Engineering, University of the Ryukyus,
Okinawa 903-0213, Japan
* Correspondence: [email protected] (S.S.R.); [email protected] (T.S.)

Abstract: This paper proposes a shunt active filter (SAF) for harmonic mitigation and reactive power
compensation in an induction heating (IH) system. The high-frequency switching in the resonant
inverter of the IH system generates a considerable number of high-frequency harmonics. The latter
flow back to the supply side and causes a wide variety of problems. The predominant ones are the
deterioration of the power quality, distortion in the grid voltage, and malfunctioning of the protective
equipment. These harmonics need to be attenuated as per the IEEE 519-1992 and IEEE 519-2014
standards. To overcome these problems, an SAF based on the modified version of instantaneous
power theory was placed between the power supply and the IH equipment. Moreover, the proposed
model could work in unbalanced and non-sinusoidal input voltage conditions, as well as dynamic
Citation: Raman, R.; Sadhu, P.K.;
Kumar, R.; Rangarajan, S.S.; conditions with changing reference currents. The feasibility of the proposed SAF-based IH system
Subramaniam, U.; Collins, E.R.; was verified by a series of simulation results and a comparative analysis of the THD of the input
Senjyu, T. Feasible Evaluation and currents. The power quality issues were successfully addressed, which signifies the ability and
Implementation of Shunt Active effectiveness of the proposed model to mitigate the problems caused by harmonics and enhance the
Filter for Harmonic Mitigation in power factors.
Induction Heating System. Electronics
2022, 11, 3464. https://doi.org/ Keywords: EMI; harmonics; induction heating; PSIM; RFI; shunt active filter
10.3390/electronics11213464

Academic Editor: Hany M. Hasanien

Received: 4 October 2022 1. Introduction

Accepted: 22 October 2022
In recent decades, shunt active filters (SAFs) have greatly attracted researchers inves-
Published: 26 October 2022
tigating harmonics elimination in induction heating (IH) systems [1]. The latter plays an
Publisher’s Note: MDPI stays neutral extremely important role in various industrial applications. With the advancements in
with regard to jurisdictional claims in power electronics technology, they have proved to be a great boon for many industries and
published maps and institutional affil- are expected to become the heart and soul of many more soon [2]. Various new topologies
iations.
of resonant inverters have been devised for use in high-frequency induction heating appli-
cations [3]. However, their operation at high frequency is associated with some inherent
problems. Electromagnetic interference (EMI) and radio frequency interference (RFI) are
quite common problems in IH equipment [4].
Copyright: © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
When power electronic switches, such as an IGBT or a MOSFET, are operated at
This article is an open access article
high frequency, they produce a considerable number of high-frequency harmonics. These
distributed under the terms and have an immanent propensity to flow back towards the input power supply, which has
conditions of the Creative Commons a very adverse effect on the power quality [5]. Some of the most common problems are
Attribution (CC BY) license (https:// the malfunctioning of the protective equipment and the non-sinusoidal nature of the grid
creativecommons.org/licenses/by/ voltage that gets distorted due to the backflow of the high-frequency harmonics [6].
4.0/).

Electronics 2022, 11, 3464. https://doi.org/10.3390/electronics11213464 https://www.mdpi.com/journal/electronics

Electronics 2022, 11, 3464 2 of 15

Thus, their generation must be kept below a certain specified value, beyond which it
can become dangerous, and for this purpose, various regulations have been formulated
in different parts of the globe. IEEE 519 is a very commonly used harmonic restriction
standard in many countries, and it puts up a restriction on the total harmonic distortion
(THD) in the supply current [7]. Good-quality filters are required to attenuate the harmonics
and reduce the problems they cause in the input power supply. The resonant inverter
generates two types of harmonic currents, namely switching frequency harmonics and
high-frequency harmonics. Both need to be attenuated by a considerable amount to follow
the stringent harmonic standards [8].
The current harmonics can prove to be more dangerous than the voltage harmonics.
The former has a tendency to enter the grid where they are responsible for the malfunction
of sensitive apparatus connected to the bus (IEEE Working Group on Non-sinusoidal
Situations) [9]. In the past few decades, the control of harmonic currents has been a
major problem for design engineers. However, to overcome these problems, many new
technologies have come into existence. Some of the more famous are passive, active, and
hybrid filters, the compensation of magnetic flux, and the injection of DC ripple voltage.
These all have some unique qualities of harmonic attenuation, but they suffer from some
limitations as well [10]. The voltage harmonics and current distortion were considerably
reduced even with the application of passive harmonic filters, and these were the most
traditional and common methods for harmonic attenuation. However, they suffered from
the problem of poor dynamic performance, and they were able to attenuate only those
harmonic components for which they had been previously tuned. This serious drawback
restricts their application to only a particular type of load. They also have a tendency to
interact with other loads present in the system, which in turn leads to resonance conditions,
and thus the entire system can become prone to malfunctioning [11].
The attenuation qualities of active harmonic filters are much better than those of their
classical counterparts. The compensation principle of the former was proposed around
1970, after which they brought some revolutionary changes to the field of power quality
improvement. The shunt active filter (SAF) provides a great remedy to the problems
related to reactive and harmonic currents [12]. The non-linear loads produce reactive
and harmonic currents which the SAF eliminates by producing canceling currents of the
same magnitude but in a 180◦ phase shift. The repetitive controllers were quite useful for
eliminating the harmonics in periodic signals with steady-state tracking [13]. However,
problems arise when the reference current changes. This is because of the inability of the
repetitive controllers to work in the first cycle in the conditions. Thus, they are capable of
working only in steady-state conditions and fail to perform in dynamic conditions [14].
Different researchers have proposed different kinds of techniques to overcome the
problems caused by harmonics in IH equipment. Pal et al., in 2015, proposed new equip-
ment consisting of a modified half-bridge inverter topology to overcome the problems
caused by harmonics [15]. The input current waveform was highly distorted without the
use of a filter, and the THD was 44.99% while that with a filter was found to be 17.90%.
However, the presence of even this amount of harmonic current is not acceptable, as per
the stringent harmonic standards, and might lead to a wide variety of harmonic-related
problems along with the deterioration of the power quality [16].
Bojoi et al., in 2008, proposed an IH system comprising a shunt active power filter
(SAPF) having a DSP controller [17]. Proportional-sinusoidal signal integrators (P-SSIs)
were used for the purpose of controlling the current. These controllers are tuned for
different harmonic components, and they are based on the principle of selective harmonic
compensation. They generally comprise two or more current-controlled loops which
create complexity in the controller design as the proper current control with multiple
loops and with a resonant controller is quite difficult. Moreover, the parameters of the
resonant controllers need to be tuned properly, which further increases the complexity of the
controller design [18]. The slow transient response also needs to be improved, and for this
purpose, the SSI constant for the fundamental frequency controller must be increased, which
Electronics 2022, 11, 3464 3 of 15

decreases the selectivity. In most of the DSP-based current controller schemes, such as the
P-SSI (proportional-sinusoidal signal integrator), the PI-MRI (proportional integral-multi
reference frame integrator), and the P-SSI-SRF (proportional-sinusoidal signal integrator
with synchronous reference frame), a delay compensation method, starting with the 17th
harmonic itself, Is essential. Some advanced versions, such as the PI-RES proportional
integral-resonant controller) remain stable in this range of harmonic compensation [18].
However, these also suffer from stability problems when attempts are made to compensate
for 29th or higher-order harmonics [19]. Furthermore, in the case of IH equipment where
the resonance frequency is quite high, the delay imposed by the sampling time can cause
an adverse effect on the performance of the system as well as its frequency.
Herrera et al. performed a comparative analysis of the various approaches of instan-
taneous power theory and found that all of them gave similar results under balanced
and sinusoidal voltage [20]. However, none of them performed satisfactorily under un-
balanced and non-sinusoidal voltage conditions. Li et al. proposed a method to reduce
the harmonic distortion on the source side, and also to reduce the inherent impact on the
rectifier transformer [21]. Wada et al. proposed a control technique to avoid the ‘whack-
a-mole’ phenomenon caused by the PFC correction capacitors in distribution feeders of
long length [22]. Mehrasa et al. proposed a controller based on multi-level converters
employing the direct Lyapunov method to attenuate the current harmonics and compensate
for the reactive power [23]. Chaoui et al. proposed a PI-controlled SAF for power quality
improvement [24]. None of the aforementioned methods reported in the literature worked
in unbalanced and non-sinusoidal input voltage conditions or in dynamic conditions with
a changing reference current. The method discussed in this paper addresses all of these
issues, and this is one of the key novelties of the proposed scheme.
To overcome these problems, an IH system comprising a shunt active filter (SAF) based
on instantaneous power theory will be proposed in the present work. The latter, commonly
known as PQ theory, has a large number of benefits over its classical counterparts. It
offers dynamic compensation of harmonic currents as well as dynamic correction of power
factors [25]. Moreover, in systems where a neutral wire is present, the zero-sequence current
is also compensated instantaneously and dynamically. The currents supplied by the source
to the load are also reduced as well as being balanced. The instantaneous three-phase
power delivered by the source to the load is turned to a constant value [26].

2. DC Bus Capacitor vs. Active Filtering Methods

Some researchers in the past have tried to eliminate the harmonics generated by the
resonant converter by simply putting a DC bus capacitor between the upper and lower
rails of the resonant converter. However, further research has emphasized the use of active
filtering techniques over DC bus capacitors, despite of the extra cost added by the former
in many applications.
Li et al., in 2009, proposed the use of an active filtering method as a substitute for
the DC bus capacitor [27]. They also evaluated the practical feasibility of active filtering
methods. Their study concluded that a system that contains only capacitors to address
harmonic issues has poor reliability. The capacitor has a poor shelf life, and the temperature
rating is also quite limited. Their study also concluded that the passive capacitor filter can
ring with the filter inductance with injected current harmonics. On the other hand, the
spurious resonance effect in the DC-link has easily been removed using the active filtering
technique. Moreover, an elastic performance is ensured due to the programmable nature of
the filtering function. It was also observed that the size of the capacitor in microfarads is
dictated by the rise in temperature that occurs due to the heating caused by the harmonic
currents, and not by the bandwidth of the filter.
Bose and Kastha, in 1991, also proposed the replacement of capacitors in a power
electronic system by a high-frequency current-fed type active filter [28]. If the capacitor is
used as a filter, it will be a one-pole type filter that offers attenuation of only 20 db/decade
above its cutoff frequency. When the harmonic spectrum is large, the capacitor is unable to
Electronics 2022, 11, 3464 4 of 15

absorb significant lower-order current harmonics if it is performing satisfactorily with the

higher-order harmonics. The problem of conducted EMI requires the use of filters designed
for specific cases. Its behavior changes in different frequency ranges and diverse solutions
are required to address such complicated problems.
The aforementioned research also concluded that the cost becomes extravagant when
a single filter is used over a broad range of frequencies and at different power levels.
Moreover, the counterpart active filtering technique effectively compensates for both power
factors and harmonic distortions and provides benefits such as the use of a single unit to
compensate for numerous non-linear loads [29].

3. Proposed Structure of IH Equipment with SAF

Figure 1 shows the basic structure of the IH equipment with the application of a SAF at
the input side. It also shows a controller based on instantaneous power theory that is used
to control the compensation currents produced by the SAF. A high-frequency alternating
magnetic field is produced by the series resonant inverter that cuts the workpiece and
produces a heating effect due to the eddy current loss and hysteresis loss. The latter occurs
mainly because of the continuous magnetization and demagnetization of the workpiece
by the alternating magnetic field. One more factor that greatly governs the heating phe-
nomenon is the skin effect, due to which the current is mainly confined in the outer surface
of the workpiece and thus produces more heat there.

Figure 1. Proposed design of a shunt active filter for induction heating equipment.

The instantaneous power theory that was introduced by Akagi was used for determin-
ing the harmonic distortions [30,31].
According to the notations used in Figure 1:
√
(i I H ) a = ( I I H1 ) a
2 Sinωt + (ei I H ) a (1)
√  ◦

(i I H )b = ( I I H1 )b 2 Sin ωt − 120 + (ei I H )b (2)
√  ◦

(i I H )c = ( I I H1 )c 2 Sin ωt + 120 + (ei I H )c (3)
Electronics 2022, 11, 3464 5 of 15

where (IIH1 )a , (IIH1 )b , and (IIH1 )c represent, respectively, the root mean square values of
the fundamental component of currents (iIH )a , (iIH )b , and (iIH )c , that are flowing in the IH
System. The high-frequency harmonics generated because of the power electronic switches
are given by the quantities
(ei I H ) a , (ei I H )b (ei I H )c .
The various quantities may be represented as follows:
       
(e I H ) a (i In ) a (i I H ) a (iSAF ) a
eI H = (e I H )b , i In = (i In )b , i I H = (i I H )b , iSAF = (iSAF )b  (4)
(e I H )c (i In )c (i I H )c (iSAF )c

Now, using the Clarke transformation, the three-phase quantities that are in the a-b-c
coordinate system can be converted to the α-β-0 coordinate system:
1
 r  1 −√21 (e I H ) 


(e I H )α −
√2 a
(e I H ) β  = 2 
 0
3
− 23 (e I H )b 

(5)
3 √1 2
( e I H )0 √1 √1 (e I H )c
2 2 2
   
(e I H )α (e I H ) a
(e I H ) β  = [( Mat) ](e I H )  (6)
IH b
( e I H )0 (e I H )c
where:
1
−√12
 
r 1 −
√2
2 3
[( Mat) I H ] =  0 − 23  (7)

3 √1 2
1
√ √1
2 2 2

Similarly, the three phase currents that are flowing in the induction heating system can
be converted from the a-b-c coordinate system to the α-β-0 coordinate system as follows:
   
(i I H )α (i I H ) a
(i I H ) β  = [( Mat) ](i I H )  (8)
IH b
( i I H )0 (i I H )c

The three-phase instantaneous real power can be expressed in the a-b-c coordinate
system as well as in the α-β-0 coordinate system as follows:

p abc = p a + pb + pc (9)

p abc = (e I H ) a (i I H ) a + (e I H )b (i I H )b + (e I H )c (i I H )c (10)
pαβ0 = (e I H )α (i I H )α + (e I H ) β (i I H ) β + (e I H )0 (i I H )0 (11)
pαβ0 = p I H + p0 (12)
The two components of pαβ0 are the instantaneous real power (pIH ) and the instanta-
neous zero sequence power (p0 ).

p I H = ( e I H )α (i I H )α + ( e I H ) β (i I H ) β (13)

p0 = ( e I H )0 ( i I H )0 (14)
Because of the absence of the neutral wire in the three-phase system, the zero sequence
components (eIH , iIH , and p0 ) can be assumed to be missing.
Thus, r " #
1
2 1 −√ 2 −√21
[( Mat) I H ] = 3 (15)
3 0 2 − 23
Electronics 2022, 11, 3464 6 of 15

r " #(e ) 
1
− −√21  I H a 
 
(e I H )α 2 1 √2
= (e I H ) (16)
(e I H ) β 3 0 2
3
− 23 (e )b
IH c
r " # (i ) 

1
− −√21  I H a 
 
(i I H )α 2 1 √2
= (i I H ) (17)
(i I H ) β 3 0 2
3
− 23 (i )b
IH c

pαβ = p I H = (e I H )α (i I H )α + (e I H ) β (i I H ) β (18)
The instantaneous imaginary power can also be expressed in the α-β coordinate
system as follows:
→ → → → →
q I H k = ( e I H )α i × (i I H ) β j + ( e I H ) β j × (i I H )α i (19)

The above expressions may be written in matrix form as follows:

 " #
( e I H )α ( e I H ) β (i I H )α
 
pI H
= (20)
qI H −(e I H ) β (e I H )α (i I H ) β

 " # −1 
(e I H )α (e I H ) β
 
(i I H )α pI H
= (21)
(i I H ) β −(e I H ) β (e I H )α qI H
" #
(e I H )α −(e I H ) β p I H
  
(i I H )α 1
= (22)
(i I H ) β (e I H )2α + (e I H )2β (e I H ) β (e I H )α qI H
" #
(e I H )α −(e I H ) β
  
(i I H )α 1 pI H
=
(i I H ) β (e I H )2 +(e I H )2β (e I H ) β (e I H )α 0
" α # (23)
(e I H )α −(e I H ) β

1 0
+ 2 2
(e I H )α +(e I H ) β (e I H ) β (e I H )α qI H

Since the instantaneous active and reactive powers constitute components that are
associated with both the fundamental and harmonic components of voltages and currents,
they can be represented as follows:

p I H = p I H + peI H (24)

q I H = q I H + qeI H (25)
where:
p I H represents the fundamental component of instantaneous active power. It corre-
sponds to the real power that is absorbed by the induction heating equipment, and it is
associated with the fundamental components of current iIH and voltage eIH ;
peI H represents the harmonic component of instantaneous active power. It corresponds
to the real power that is absorbed by the induction heating equipment, and it is associated
with the harmonic components of current iIH and voltage eIH ;
q I H represents the fundamental component of instantaneous reactive power, and it is
associated with the fundamental components of current iIH and voltage eIH ;
qeI H represents the harmonic component of instantaneous reactive power, and it is
associated with the harmonic components of current iIH and voltage eIH .
We will now consider a new quantity, iIHF , that represents the fundamental active
component of the current absorbed by the induction heating equipment. It can be expressed
in the α-β coordinate system in the form of a matrix as follows:

 " # −1 
(e I H )α (e I H ) β
 
(i I HF )α pI H
= (26)
(i I HF ) β −(e I H ) β (e I H )α 0
Electronics 2022, 11, 3464 7 of 15

" #
(e I H )α −(e I H ) β

pI H

(i I HF )α
 (e I H ) β (e I H )α 0
= (27)
(i I HF ) β (e I H )2α + (e I H )2β
Using the above expressions, the reference current (iSAF ∗ ) can be obtained in the α-β

coordinate system as follows:

" # 
∗ )
(iSAF
  
(i I H )α (i I HF )α
∗ )
α = − (28)
(iSAF β (i I H ) β (i I HF ) β

Putting the values of Equations (22) and (27) in Equation (29), we get:
" # " # −1   " # −1 
∗ )
(iSAF (e I H )α (e I H ) β (e I H )α (e I H ) β

pI H pI H
∗
α = − (29)
(iSAF ) β −(e I H ) β (e I H )α qI H −(e I H ) β (e I H )α 0
" #
(e I H )α −(e I H ) β

peI H
" #
∗ )
(iSAF (e I H ) β (e I H )α q I H + qeI H
α = (30)
∗
(iSAF ) β (e I H )2α + (e I H )2β
One of the most important components of the shunt active filter (SAF) is the inverter.
It has a capacitor whose voltage corresponds to the DC link voltage (Edc ). This voltage
needs to be regulated to keep it within appropriate levels. This process is compulsory for
the proper functioning of the SAF. The voltage (Edc ) of the DC link increases when the SAF
receives the energy that is stored in the capacitor. If it does not receive energy, the capacitor
will discharge, and thus there is a reduction in Edc . A regulating power (Pregulating ) which
is a function of the reference voltage (Eref ) and the DC link voltage (Edc ) is used to regulate
the value of the instantaneous active power.
( peI H )regulated = peI H − pregulating (31)

pregulating = k ( Ere f − Edc ) (32)

where k is a constant.
( peI H )regulated = peI H − k( Ere f − Edc ) (33)
Thus, considering the regulation of the voltage (Edc ) of the DC link and using the
above expressions, the reference currents can be expressed as follows:
" #" #
(e I H )α −(e I H ) β ( peI H )regulated
" #
∗ )
(iSAF (e I H ) β (e I H )α q I H + qeI H
α = (34)
∗
(iSAF ) β ( e ) + ( e )2
2
IH α IH β

The reference currents that are in the α-β coordinate system can be converted back to
the a-b-c coordinate system as follows:
 
 ∗ ) 
(iSAF 1 0 " #
√ ∗ )
r
∗
a 2 − 3 (iSAF
 (i
SAF )b
=  1 α (35)
 
∗
3 −21 2√ (iSAF ) β

∗
(iSAF )c 3
2 − 2

4. The Case of Unbalanced and Non-Sinusoidal Input Voltage

The phase voltages (eIH )a , (eIH )b , and (eIH )c at the input of the IH equipment might be
balanced or unbalanced. The former mainly contains positive sequence component (eIH )+1
while the latter also contains negative and zero sequence components at the fundamental
frequency. Moreover, the harmonics from any sequence components are also present in
the latter case (IEEE Working Group on Non-sinusoidal Situations, 1996). Thus, a positive
Electronics 2022, 11, 3464 8 of 15

sequence voltage detector (PSVD) needs to be employed to detect the fundamental positive
sequence components of the input phase voltages (eIH )a , (eIH )b , and (eIH )c . The controller for
the IH equipment based on instantaneous power theory shown in Figure 2 was upgraded by
adding a PSVD block to enable it work in unbalanced and non-sinusoidal voltage conditions.
The PSVD block instantaneously extracts the fundamental positive sequence voltages
(eIH )a , (eIH )b , and (eIH )c . Thus, the SAF compensates for the input current of the IH
equipment in such a manner that only the active part of the fundamental positive sequence
current (iIH )+1 is supplied by the power supply. The working of the PSVD, which consists
of a phase-locked loop (PLL)circuit, is also very closely related to the PQ theory. The
basic structure of the PSVD is shown in Figure 3. The input phase voltages that are in the
a-b-c coordinate system can be transformed into the α-β coordinate system to get (eIH )α
and (eIH )β . The PLL circuit produces auxiliary currents (iIH )0 α and (iIH )0 β , which have
arbitrary magnitudes and are derived only from the positive sequence current (IIH )+1 at
the fundamental frequency. These voltages and currents are used to compute the real and
imaginary auxiliary powers p0 IH and q0 IH , respectively.

Figure 2. Proposed structure of the controller based on instantaneous power theory.

Figure 3. Proposed structure of a positive sequence voltage detector.

The IH system currents in the α-β coordinate system can be expressed in terms of
symmetrical components as follows:
∞ √ ∞ √
(i I H )α = ∑ 3( I I H )+m Sin(ωm t + θ+m ) + ∑ 3( I I H )−m Sin(ωm t + θ−m ) (36)
m =1 m =1

∞ √ ∞ √
(i I H ) β = ∑ − 3( I I H )+m Cos(ωm t + θ+m ) + ∑ 3( I I H )−m Cos(ωm t + θ−m ) (37)
m =1 m =1
Electronics 2022, 11, 3464 9 of 15

Considering the case of the fundamental (m = 1) positive sequence component,

Equations (37) and (38) can be expressed as follows:
√
(i I H )α = 3( I I H )+1 Sin(ω1 t + θ+1 ) (38)
√
(i I H ) β = − 3( I I H )+1 Cos(ω1 t + θ+1 ) (39)

Since, the amplitudes of auxiliary currents (iIH )0 α and (iIH )0 β are randomly selected,
they may be considered as a unity to reduce the complexity. Moreover, the phase angles
(θ+1 ) of the aforesaid currents are also randomly selected, and thus they are taken as zero
for simplicity. The only crucial aspect is the accurate determination of the fundamental
frequency (ω1 ) by the PLL circuit.

(i I H )0 α = Sin(ω1 t) (40)

(i I H )0 β = −Cos(ω1 t) (41)

The average parts of auxiliary powers p0 IH and q0 IH are obtained only from the funda-
mental component of the positive sequence voltage (eIH )+1 . The auxiliary currents (iIH )0 α
and (iIH )0 β are also obtained only from (iIH )+1 . Moreover, because of the fundamental compo-
nent of the negative sequence voltage and the other voltage components, the harmonics are
present only in the oscillating component of the auxiliary power that is not included in the
inverse voltage calculation. The instantaneous voltages (eIH )0 α and (eIH )0 β , which represent
the fundamental components of the positive sequence voltage, reobtained as follows:
" #
(e I H )0 α (i I H )0 α −(i I H )0 β p0 I H
  
1
= (42)
(e I H )0 β 0
(i I H )0 2 + (i I H )0 2 (i I H ) β (i I H ) α
0 q0 I H
α β

Using Equations (41) and (42), the sum of the square of the auxiliary currents (iIH )0 α
and (iIH )0 β can be calculated.
2 2
(i I H )0 α + (i I H )0 β = Sin2 (ω1 t) + Cos2 (ω1 t) (43)

2 2
(i I H )0 α + (i I H )0 β = 1 (44)

The positive sequence phase voltages (eIH )0 a , (eIH )0 b , and (eIH )0 c can be obtained from
(eIH )0 0
α and (eIH ) β by using the inverse Clarke transformation as follows:

(e I H )0 a
  r " #T 
1
1 − −√12 (e I H )0 α

(e I H )0  = 2 √2 (45)
b
3 0 3
− 23 (e I H )0 β
(e I H )0 c 2

5. Phase-Locked Loop Circuit

The PLL circuit present in the PSVD continuously tracks the fundamental frequency of
the voltage fed to the IH equipment and helps it to work effectively in highly unbalanced
and non-sinusoidal voltage conditions. Figure 4 shows the functional block diagram
of a PLL Circuit which automatically finds out the phase angle and frequency of the
fundamental positive sequence component of the input phase voltage. The algorithm of
the PLL circuit is based on the expression of instantaneous power, which is hypothetical.
However, it has been assumed as a variable that has the same dimensions as that of power.
It can be expressed as follows:

p0 abc = (e I H ) a (i I H )0 a + (e I H )b (i I H )0 b + (e I H )c (i I H )0 c (46)

p0 abc = (e I H ) ab (i I H )0 a + (e I H )cb (i I H )0 c (47)

Electronics 2022, 11, 3464 10 of 15

The feedback signals (iIH )0 c (ωt) and (iIH )0 a (ωt) are generated by the PLL circuit by
simply finding the time integral of the angular frequency at the output of the PI controller.
The average value of P0 abc can be expressed as follows:
p0abc = 3(e I H )0 +1 (i I H )0 +1 CosΦ (48)

To operate it in a stable condition, the average value of P0 abc should be equal to zero.
Moreover, the feedback current (iIH )0 a (ωt) should lead the phase voltage by 90◦ .

p0abc = 0 (49)

Figure 4. Block diagram of a PLL circuit.

6. Simulation Diagram and Results

In the present work, the simulation of a high-frequency full-bridge series resonant
inverter has been carried outby considering the equivalent circuit parameters. PSIM 2022.1,
which is quite famous in the field of power electronics and motor drives, was used for the
simulation of various equivalent circuits. First, an IH system with an SAFwas proposed as
shown in Figure 1, and the mathematical model of the controller was developed. Figure 2
gives the proposed structure of the controller based on instantaneous power theory. Figure 3
shows the proposed structure of a PSVD that was incorporated to help deal with the
unbalanced input voltage conditions. Figure 4 shows the block diagram of the PLL circuit
which was used in Figure 3. Finally, all these modifications were added to the structure
outlined in Figure 2, which was based on the original instantaneous power theory. Figure 5
shows the proposed structure of a controller based on the extended instantaneous power
theory which is capable of working even in unbalanced and non-sinusoidal input voltage
conditions. Next, the simulations were performed to verify the necessity and importance of
the SAF in the IH equipment. Initially, simulations were carried out without the application
of a low pass filter (LPF), but subsequently they were carried out with a shunt active filter
(SAF)installed at the input of the IH equipment.

Figure 5. Proposed structure of the controller based on extended instantaneous power theory.
Electronics 2022, 11, 3464 11 of 15

These were finally compared to see the effect of the SAF on the harmonics present in
the input current waveform. The equivalent simulation circuit diagrams of the induction
heating equipment without a filter and with an SAF are shown in Figures 6 and 7, respec-
tively. The corresponding waveforms of the input current obtained after a simulation in
PSIM are also shown in Figures 8 and 9, respectively. The FFT analysis of these waveforms is
shown in Figures 10 and 11, respectively. The equivalent simulation circuit diagram, which
represents the basic structure of the IH equipment, consists of a three-phase sinusoidal
AC source followed by a three-phase star-delta transformer that is fed to a three-phase
uncontrolled bridge rectifier.

Figure 6. Simulation circuit diagram of the IH equipment without an LPF.

Figure 7. Simulation circuit diagram of the IH equipment with an SAF.

Figure 8. Waveform of the input phase current of the IH equipment without an LPF.
Electronics 2022, 11, 3464 12 of 15

Figure 9. FFT of the input phase current of the IH equipment without an LPF.

Figure 10. Waveform of the input phase current of the IH equipment with an SAF.

Figure 11. FFT of the input phase current of the IH equipment with an SAF.

The output of the rectifier acts as an input to a high-frequency series resonant full-
bridge inverter. The modeling of the IH coil was carried out with a combination of an
inductor and a resistor in series connection with an external capacitor. The series combina-
tion of the RLC acts as a load to the high-frequency inverter, which is resonant in nature.
Along with these essential components, a few more components, such as voltmeters and
ammeters, were used for measurement purposes.

7. Calculation of THD from the FFT Analysis of the Input Phase Current Waveform
A fast Fourier transform, commonly known as FFT, is an algorithm that is meant for
sampling a signal over a specified time period and then dividing it into various frequency
Electronics 2022, 11, 3464 13 of 15

components. The waveforms of the input phase current in the time domain were converted
to the frequency domain with the FFT option available in PSIM.

7.1. Without Filter

The FFT analysis of the input phase current waveform showed the presence of
10 predominant harmonic components that are responsible for harmonic distortions. The
measure option available in PSIM was used to find the RMS values of various components
present in the input current waveform. The fundamental component of the current was
found to be 50.403 A without the application of a filter. The highest magnitude of harmonic
current was found to be 13.372 A, while the lowest measured value was 1.049 A. The other
harmonics having magnitudes lower than this were neglected as they hardly contributed
to any current distortions. The THD is expressed as follows:
s
∞
∑ In,r.m.s
2
n =2
THD =
I1,r.m.s

The values of the magnitudes of various harmonic components are expressed as follows:

THD
= (13.3722 + 5.8172 + 2.9142 + 1.4722 + 1.3112 + 1.2892


+1.2192 + 1.1492 + 1.0572 + 1.0492 )0.5 /50.403 × 100

= 30.92%

7.2. With a Shunt Active Filter

The FFT analysis of the input phase current waveform obtained after the application
of an SAF shows that most of the predominant harmonic components were attenuated. The
fundamental component was found to be 59.134 A while the harmonic components were
0.729 A and 0.312 A. s
∞
2
∑ In,r.m.s
n =2
THD = I1,r.m.s

THD
n o
0.5
= (0.7292 + 0.3122 ) /59.134 × 100
= 1.34%

8. Conclusions and Future Work

In the present paper, an IPT-based shunt active filter for harmonic mitigation in an
induction heating system has been proposed. Moreover, an exhaustive analysis of its
performance under unbalanced and non-sinusoidal input voltage conditions, as well as
dynamic conditions with changing reference currents, has been carried out. The modified
version of the IPT caters to the specific needs of IH systems for industrial use. The control
strategy followed in the proposed method compensated for the harmonic currents and
reactive power of the IH system, and the DC bus voltage of the SAF was almost maintained
at the reference value under all conditions. A detailed mathematical analysis of the pro-
posed model was carried out, followed by its validation via simulations in PSIM. The FFT
analysis of the input phase current waveform shows that almost all the high-frequency
harmonic components were attenuated with the application of an SAF. Its THD was found
to be 30.92% without a filter, while the THD with an SAF was just 1.34%. The results
show that the developed model could reduce the THD to a level of less than 5%, which
is in accordance with IEEE 519-1992 and IEEE 519-2014 standards. This would also be
considerably less than the previous models proposed in the literature.
In order to validate the simulation results, it will be necessary to set up the shunt
active filter based IH system hardware in future work. Its practical implementation and its
Electronics 2022, 11, 3464 14 of 15

measurements will be realized using a dSPACE DS1104 R&D controller with TMS320F240
as a slave DSP. Moreover, its ability to provide harmonic mitigation, reactive power compen-
sation, and the balancing of unbalanced and non-linear loads will be under consideration.
Experimental results will be compared with the simulation results to determine the efficacy
and viability of the proposed method.

Author Contributions: Conceptualization, R.R. and P.K.S.; methodology, R.R., R.K. and P.K.S.;
software, R.R., R.K. and P.K.S.; validation, R.R. and P.K.S.; formal analysis, R.R., R.K. and P.K.S.; inves-
tigation, R.R. and P.K.S.; resources, R.R. and P.K.S.; data curation, R.R. and P.K.S.; writing—original
draft preparation, R.R. and P.K.S.; writing—review and editing, R.R.; visualization, R.R., P.K.S., S.S.R.,
U.S., E.R.C. and T.S.; supervision, P.K.S., S.S.R., U.S., E.R.C. and T.S.; project administration, R.R.,
P.K.S., S.S.R., U.S., E.R.C. and T.S. All authors have read and agreed to the published version of
the manuscript.
Funding: This research received no external funding.
Data Availability Statement: Not applicable.
Acknowledgments: The authors would like to thank the Indian Institute of Technology (ISM),
Dhanbad for providing the necessary facilities to carry out this research.
Conflicts of Interest: The authors declare no conflict of interests.

References
1. Kumar, A.; Kumar Sadhu, P.; Kumar Mohanta, D.; Bharata Reddy, M.J. An Effective Switching Algorithm for Single Phase Matrix
Converter in Induction Heating Applications. Electronics 2018, 7, 149. [CrossRef]
2. Sarnago, H.; Lucía, O.; Mediano, A.; Burdío, J.M. Class-D/DE dual-mode-operation resonant converter for improved-efficiency
domestic induction heating system. IEEE Trans. Power Electron. 2013, 28, 1274–1285. [CrossRef]
3. Acero, J.; Burdio, J.M.; Barragan, L.A.; Navarro, D.; Alonso, R.; Ramon, J.; Monterde, F.; Hernandez, P.; Llorente, S.; Garde, I.
Domestic induction appliances. IEEE Ind. Appl. Mag. 2010, 16, 39–47. [CrossRef]
4. Khlifi, K.; Haddouk, A.; Hlaili, M.; Mechergui, H. Harmonic Pollution Caused by Non-Linear Load: Analysis and Identification.
Int. J. Energy Environ. Eng. 2018, 12, 510–517.
5. Ko, A.; Swe, W.; Zeya, A. Analysis of harmonic distortion in non-linear loads. Int. J. Comput. Internet Manag. 2011, 19, 661–666.
6. Singh, M.; Mahapatra, S. Implementation of Passive Filters for Harmonics Reduction. Int. J. Adv. Sci. Technol. 2015, 78, 1–12.
[CrossRef]
7. IEEE. Standard 519. IEEE recommended practice and requirement for harmonic control in electric power systems. In IEEE Std
519-2014 (Revision of IEEE Std 519-1992); IEEE: Piscataway, NJ, USA, 2014; pp. 1–29.
8. Das, J.C. Passive filters—Potentialities and limitations. IEEE Trans. Ind. Appl. 2004, 40, 232–241. [CrossRef]
9. Abaza, A.; El-Sehiemy, R.A.; Said, M.; Ghoniem, R.M.; Barakat, A.F. Implementation of an Electronically Based Active Power
Filter Associated with a Digital Controller for Harmonics Elimination and Power Factor Correction. Electronics 2022, 11, 2205.
[CrossRef]
10. Sao, J.K.; Naayagi, R.T.; Panda, G.; Patidar, R.D.; Swain, S.D. SAPF Parameter Optimization with the Application of Taguchi SNR
Method. Electronics 2022, 11, 348. [CrossRef]
11. Panati, P.U.; Ramasamy, S.; Ahsan, M.; Haider, J.; Rodrigues, E.M.G. Indirect Effective Controlled Split Source Inverter-Based
Parallel Active Power Filter for Enhancing Power Quality. Electronics 2021, 10, 892. [CrossRef]
12. Hoon, Y.; MohdRadzi, M.A.; MohdZainuri, M.A.A.; Zawawi, M.A.M. Shunt Active Power Filter: A Review on Phase Synchro-
nization Control Techniques. Electronics 2019, 8, 791. [CrossRef]
13. Das, S.; Ray, P.; Sahoo, A.; Ramasubbareddy, S.; Babu, T.; Kumar, N.; Elavarasan, R.; Mihet-Popa, L. A Comprehensive Survey on
Different Control Strategies and Applications of Active Power Filters for Power Quality Improvement. Energies 2021, 14, 4589.
[CrossRef]
14. De Araujo Ribeiro, R.L.; de Azevedo, C.C.; de Sousa, R.M. A robust adaptive control strategy of active power filters for
powerfactorcorrection, harmonic compensation, and balancing of nonlinear loads. IEEE Trans. Power Electron. 2011, 27, 718–730.
[CrossRef]
15. Pal, P.; Sadhu, P.K.; Pal, N.; Sanyal, S. An exclusive design of EMI-RFI suppressor for modified half bridge inverter fitted induction
heating equipment. Int. J. Mechatron. Electr. Comput. Technol. IJMEC 2015, 5, 2084–2100.
16. Saad, S.; Zellouma, L. Fuzzy logic controller for three-level shunt active filter compensating harmonics and reactive power. Electr.
Power Syst. Res. 2009, 79, 1337–1341. [CrossRef]
17. Bojoi, R.; Griva, G.; Profumo, F.; Cesano, M.; Natale, L. Shunt active power filter implementation for induction heating applications.
COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2008, 27, 445–459. [CrossRef]
18. Green, T.C.; Marks, J.H. Control techniques for active power filters. IEEE Proc. Electr. Power Appl. 2005, 152, 369–381. [CrossRef]
Electronics 2022, 11, 3464 15 of 15

19. Zainuri, M.A.A.M.; Radzi, M.A.M.; Soh, A.C.; Mariun, N.; Rahim, N.A. Dc-link capacitor voltage control for single-phase shunt
active power filter with step size error cancellation in self-charging algorithm. IET Power Electron. 2016, 9, 323–335. [CrossRef]
20. Herrera, R.S.; Salmeron, P.; Kim, H. Instantaneous reactive power theory applied to active power filter compensation: Different
approaches assessment and experimental results. IEEE Trans. Ind. Electron. 2008, 55, 184–196. [CrossRef]
21. Li, Z.; Hu, T.; Abu-Siada, A. A minimum side-lobe optimization window function and its application in harmonic detection of an
electricity gird. Energies 2019, 12, 2619. [CrossRef]
22. Wada, K.; Fujita, H.; Akagi, H. Considerations of a shunt active filter based on voltage detection for installation on a long
distribution feeder. IEEE Trans. Ind. Appl. 2002, 38, 1123–1130. [CrossRef]
23. Mehrasa, M.; Pouresmaeil, E.; Akorede, M.F.; Jørgensen, B.N.; Catalão, J.P. Multilevel converter control approach of active power
filter for harmonics elimination in electric grids. Energy 2015, 84, 722–731. [CrossRef]
24. Chaoui, A.; Gaubert, J.P.; Krim, F.; Champenois, G. PI controlled three-phase shunt active power filter for power quality
improvement. Electr. Power Compon. Syst. 2007, 35, 1331–1344. [CrossRef]
25. Waware, M.; Agarwal, P. A Review of Multilevel Inverter Based Active Power Filter. Int. J. Comput. Electr. Eng. 2011, 3, 196–205.
[CrossRef]
26. Shyu, K.-K.; Yang, M.-J.; Chen, Y.-M.; Lin, Y.-F. Model Reference Adaptive Control Design for a Shunt Active-Power-Filter System.
IEEE Trans. Ind. Electron. 2008, 55, 97–106. [CrossRef]
27. Li, S.; Ozpineci, B.; Tolbert, L.M. Evaluation of a current source active power filter to reduce the dc bus capacitor in a hybrid
electric vehicle traction drive. In Proceedings of the 2009 IEEE Energy Conversion Congress and Exposition, San Jose, CA, USA,
20–24 September 2009; pp. 1185–1190. [CrossRef]
28. Bose, B.K.; Kastha, D. Electrolytic capacitor elimination in power electronic system by high frequency active filter. In Proceedings
of the Conference Record of the 1991 IEEE Industry Applications Society Annual Meeting, Dearborn, MI, USA, 28 September–
4 October 1991; Volume 1, pp. 869–878. [CrossRef]
29. Abbas, A.S.; El-Sehiemy, R.A.; Abou El-Ela, A.; Ali, E.S.; Mahmoud, K.; Lehtonen, M.; Darwish, M.M. Optimal harmonic
mitigation in distribution systems with inverter based distributed generation. Appl. Sci. 2021, 11, 774. [CrossRef]
30. Akagi, H. Active harmonic filters. Proc. IEEE 2005, 93, 2128–2141. [CrossRef]
31. Akagi, H.; Watanabe, E.H.; Aredes, M. Instantaneous Power Theory and Applications to Power Conditioning; IEEE Press: Piscataway,
NJ, USA, 2007.