Linear and Nonlinear System Modeling
Linear and Nonlinear System Modeling
See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Linear and Nonlinear
System Modeling
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Phillip Carmical (pcarmical@scrivenerpublishing.com)
Martin Scrivener (martin@scrivenerpublishing.com)
100 Cummings Center, Suite 541J
Beverly, MA 01915-6106
Publishers at Scrivener
Scrivener Publishing
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Linear and Nonlinear
Souvik Ganguli
Tamal Roy
Edited by
and
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
This edition first published 2024 by John Wiley & Sons, Inc., 111 River Street, Hoboken, NJ 07030, USA
and Scrivener Publishing LLC, 100 Cummings Center, Suite 541J, Beverly, MA 01915, USA
© 2024 Scrivener Publishing LLC
For more information about Scrivener publications please visit www.scrivenerpublishing.com.
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or
transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or other-
wise, except as permitted by law. Advice on how to obtain permission to reuse material from this title
is available at http://www.wiley.com/go/permissions.
For details of our global editorial offices, customer services, and more information about Wiley prod-
ucts visit us at www.wiley.com.
ISBN 978-111984-742-7
Set in size of 11pt and Minion Pro by Manila Typesetting Company, Makati, Philippines
10 9 8 7 6 5 4 3 2 1
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Contents
Preface xi
1 Assessment of Faults in Hybrid System Connected
with Main Grid 1
Aveek Chattopadhyaya, Niladri Mukherjee
and Surajit Chattopadhyay
1.1 Introduction 1
1.2 Hybrid System Connected with Main Grid 3
1.3 FFT Results in Different Conditions, Respective Bar Diagram,
and Observations 4
1.4 Inter-Harmonic Group Analysis, Results, and Observations 4
1.5 Statistical Parameter Analysis Based on Discrete Wavelet
Transform, Results, and Observations 7
1.6 Algorithm to Determine Non-Identical Conditions 9
1.7 Specific Outcome of This Chapter 11
1.8 Conclusions 12
References 12
2 Diversified Harmonics Modeling for Power System
Stability Analysis 15
Tamal Roy, Debopoma Kar Ray and Surajit Chattopadhyay
2.1 Introduction 15
2.2 Classification 16
2.2.1 Steady-State Stability 16
2.2.2 Transient Stability 16
2.2.3 Dynamic Stability 17
2.3 Power Equation 18
2.4 Maximum Power 19
2.5 Nonlinearity and Harmonics 19
2.6 Active Power, Load Angle, and Reactance 20
2.7 Effects of Harmonics on Stability Model 21
v
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
vi Contents
xi
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
welcome. We are grateful to all who encouraged and helped us directly or
Tamal Roy
Suman Lata Tripathi
Souvik Ganguli
Kolkata, India
June 2024
indirectly during the preparation of the manuscript.
xii Preface
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
1
Assessment of Faults in Hybrid
System Connected with Main Grid
Aveek Chattopadhyaya1*, Niladri Mukherjee2 and Surajit Chattopadhyay 3
3
Department of Electrical Engineering, GaniKhan. C. Institute of Engineering and
Technology, Malda, India
Abstract
In this chapter, fault assessment in grid-connected microgrid (solar and wind-
based hybrid system) was performed by signal processing-based waveform anal-
ysis. For this analysis, wind and solar system-based grid-connected microgrid
(hybrid system) is considered. System current waveforms were captured for fault
identification. Total harmonic distortion (THD) based on fast Fourier transform
(FFT) and inter-harmonic group analysis were considered as tools to identify the
presence of different harmonics in non-identical environments. In Park plane,
one space vector was also calculated from the captured three-phase current sig-
nals. Different parameters based on discrete wavelet transform (DWT) were also
computed from space vector for assessing the symmetrical and unsymmetrical
faults in hybrid system, which is connected to the main grid.
1.1 Introduction
The demand for renewable energy resources has escalated significantly
over the past few years, due to the fact that the fossil fuels will run out in
the near future and the harmful environmental effects of the fossil fuels. In
hybrid system, wind & solar systems are used for generating power where
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(1–14) © 2024 Scrivener Publishing LLC
1
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
2 Linear and Nonlinear System Modeling
wind turbines generators are used for wind generation and the solar pan-
els are used for solar generation. The power generation of solar and wind
system strongly depends on the weather conditions. Distributed generation
plays a vital role in today’s huge energy demand scenario. Hybrid genera-
tion system is a system, which is having renewable/non-renewable sources.
Hybrid systems are interconnected with renewable/non-renewable sources
to supply power to local load and grid/micro grid. Different issues like con-
trol, optimal sizing, grid connection of hybrid systems having solar photo-
voltaic (PV), wind power, fuel cells etc. have been vividly ventilated by the
researchers [1]. Different issues regarding benefits and drivers of ‘hybrid
renewable energy systems (HRES)’ for off grid power generation has been
discussed in [2]. A simulink model is developed of grid coupled hybrid sys-
tem in MATLAB where the hybrid system consists of wind and photovoltaic
systems [3]. To simulate a hybrid system some basic systems are required.
In PV/wind, hybrid systems three important parts are needed [4]. Firstly,
PV/wind system, wind energy conversion (WEC) system, power elec-
tronics devices and MPPT algorithm [3]. The hybrid system encompasses
PV-wind system, DC/DC converter and energy (wind) conversion system
[4]. E. Natsheh et al. proposed a control strategy for the optimization of
smart grid performance where a hybrid system is connected with a grid [5].
In [6], a modular PV and fuel cell based grid connected hybrid system has
been developed in MATLAB and ‘160 Wp solar module’ and ‘5 kW solid
oxide fuel cell (SOFC)’ was used for simulation. N. T. Pathan et al. discussed
steady state performance of a hybrid system connected with main grid [7].
In this study permanent - magnet - synchronous – generator has been used
for wind turbine. Two different system (wind and solar) are synchronized
to supply a common grid and a local load. The output of the system fed to
a common grid and a local load while solar system generates dc power and
wind system produces pulsating AC. The main difficulty is synchronization.
The output power from wind system is converting into DC then connected
with solar system. Afterword it is fetched to a PWM inverter and finally syn-
chronized [9]. For the reliable operation of grid connected microgrid system,
it is very much obvious that, different faults should be detected properly as
early as possible. In [8], authors designed and simulated and studied hybrid
system consists of wind and solar in isolated mode [8]. A. González et al.
suggested sensitivity analysis based optimally sizes of ‘grid-connected photo-
voltaic–wind power systems’ in [10]. Using this system, cost of electricity has
been lowered as compared to purchase of electricity. S. Swain proposed an
Active Crowbar Protection (ACB_P) system to improve the system's power
quality by enhancing the DFIG's fault-ride through (FRT) capacity. The sug-
gested scheme's functionality is confirmed by running simulations on a 1.7
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Assessment of Faults in Hybrid System 3
THREE
IRRADIANCE PHASE
FAULT
Ir
PV GRID
BOOST
TRANSMISSION LINE REACTOR
MODULE CONVERTER INVERTER
Temp
V,I
TEMPERATURE
MPPT
Wind Speed
Wind Universal
Generator Bridge
rectifier
Pitch Angle
Figure 1.1 Single line representation of hybrid system connected with main grid.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
4 Linear and Nonlinear System Modeling
then directly connected to the output of the boost converter of the solar
system for proper synchronization because both the output is DC. In
DC, the synchronizing process is easier than the ac synchronizing pro-
cess because in the ac system, the main concerns are phase sequence and
frequency mismatch. To overcome these problems, synchronization was
done through DC. The combined output of the hybrid system was then
fed to the inverter. Inverter output voltage was then stepped up to 25 kV
through a step-up transformer. The output of the transformer was trans-
mitted through 15-km line and fed to a 25-kV grid. Here, non-identical
faults were contemplated in the grid side. In fault and normal situations, all
three-phase currents of bus (after transformer) was considered at proper
sampling frequency for analysis. The three-phase current was captured and
normalized then using Park transformation matrix, one space vector was
found. After finding out the current space vector (in Park plane), it was
processed by DWT. In DWT, the mother wavelet was considered as ‘db4’
as it is well suited for this analysis. From the DWT coefficients, different
parameters like kurtosis, mean, RMS, and skewness were calculated in
different conditions to detect and identify the grid side disturbances. Fast
Fourier transform (FFT) and inter-harmonics groups were also computed
on captured system current signals to assess the abnormal conditions.
100
80
60
40
20
0
Normal L-G L-L-G L-L-L L-L
DC 0.23 23.67 13.9 136.01 38.49
300
250
200
150
100
50
0
Normal L-G L-L-G L-L-L L-L
THD 1.04 9.26 6.02 287.38 99.2
Magnitude
20
15
10
5
0
Ig1 Ig2 Ig3
Normal 12 12 12
LG 13.32 12.47 12.243
LLG 13.32 12.4715 12.25
LL 14.724 15.0045 15.8348
LLL 20.86 19.3341 21.5753
Figure 1.4 Group analysis (inter-harmonics) of ‘C’ phases current in non-identical situations.
40
LLG FAULT
30 LLL FAULT
LG FAULT
20
10
0
1 2 3 4 5 6 7 8 9
DWT decomposition level----->
(a)
0.78
0.77
LL FAULT
LLG FAULT LG FAULT
0.76
NORMAL CONDITION
0.75
LLL FAULT
0.74
0.73
1 2 3 4 5 6 7 8 9
DWT decomposition level----->
(b)
Figure 1.5 (a) Kurtosis values of approximation coefficients. (b) Mean values of
approximation coefficients.
difference pattern between normal and different fault states. Figures 1.8(a)
and (b) shows the RMS and skewness of detail coefficients results. Up
to level two, a very small difference between normal and other fault cir-
cumstances are shown in Figure 1.8(a). The difference between level two
and four has remained constant for L-G problems, but it has increased in
character for other types of conditions. The highest and minimum differ-
ences between conventional and L-G and L-L-L fault situations have been
reported. Figure 1.8(b) shows some zig-zag changes of skewness values in
different conditions.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Assessment of Faults in Hybrid System 9
2.5
2
LG FAULT
1.5 LL FAULT
1
0.5
NORMAL CONDITION LG FAULT
0
1 2 3 4 5 6 7 8 9
DWT decomposition level---->
(a)
6
LLL FAULT
5 LL FAULT LLG FAULT
1
LG FAULT
0 NORMAL CONDITION
–1
–2
–3
1 2 3 4 5 6 7 8 9
DWT decomposition level----->
(b)
Figure 1.6 (a) RMS values of approximation coefficients. (b) Skewness values of
approximation coefficients.
800
600
LLG FAULT
NORMAL CONDITION
400 LL FAULT
LG FAULT
LLL FAULT
200
0
1 2 3 4 5 6 7 8 9
DWT decomposition level---->
(a)
5 LL FAULT
4
LG FAULT
3
LLG FAULT
2 NORMAL CONDITION
1 LLL FAULT
0
–1
–2
–3
1 2 3 4 5 6 7 8 9
DWT decomposition level----->
(b)
Figure 1.7 (a) Detail coefficients kurtosis values. (b) Mean values of detail coefficients.
LLG FAULT
0.15 LG FAULT
0.1
0.05
0
1 2 3 4 5 6 7 8 9
DWT decomposition level---->
(a)
Skewness values
6
Skewness values of detail coefficients--->
4
NORMAL CONDITION
LL FAULT
3
–1
1 2 3 4 5 6 7 8 9
DWT decomposition level----->
(b)
Figure 1.8 (a) RMS from DWT coefficients (detail). (b) Skewness values of detail
coefficients.
1.8 Conclusions
In this chapter, grid connected one standard hybrid PV and wind power
system were considered. The mentioned system consisted of four import-
ant segments, which are: wind system, solar system, wind energy conver-
sion (WEC) system, and double-input DC/DC converter. Fast Fourier
transform (FFT)-based THD, inter-harmonics group, Park plane, and
DWT-based different parameters were used for fault analysis in wind and
solar-based hybrid system. The minimal THD was noticed at conventional
situation and the highest THD value was observed at L-L-L symmetrical
fault situation. In this chapter, Park plane was used to extract the space
vector and then mean, RMS, skewness, and kurtosis were computed based
on DWT to identify and discriminate the grid abnormalities. Discussed
techniques can be used to develop fast and reliable protection scheme in
any grid-connected microgrid system. After studying this chapter, readers
can gather knowledge about the protection of microgrid system.
References
1. Krishna, S, K., & Kumar, S., K. (2015). A review on hybrid renewable energy
systems. Renew. Sustain. Energy Rev., Science Direct, 52, 907–916.https://doi.
org/10.1016/j.rser.2015.07.187.
2. Mohammed, S., Y., Mustafa, W., M., & Bashir, N. (2014). Hybrid renew-
able energy systems for off-grid electric power: Review of substantial
issues. Renew. Sustain. Energy Rev., Science Direct, 35, 527–539. https://doi.
org/10.1016/j.rser.2014.04.022.
3. Rana, A., & Ilyas, Md. (2015). Implementation of a wind/pv hybrid system
using matlab/simulink. Int. J. Adv. Res. Electr. Electron. Instrum. Eng., 4(7),
5941-5948. https://doi.org/10.15662/ijareeie.2015.0407019.
4. Harini, M., Ramaprabha, R., & Mathur, L., B. (2012). Modeling of grid con-
nected hybrid wind/pv generation system using matlab. ARPN J. Eng. Appl.
Sci., 7(9), 1157-1161.ISSN 1819-6608.
5. Natsheh, M., E., Albarbar, A., & Yazdani, J. (2011). Modeling and control
for smart grid integration of solar/wind energy conversion system. 2nd
IEEE PES International Conference and Exhibition on Innovative Smart Grid
Technologies (ISGT Europe), Manchester, UK, 1-8. https://doi./org/10.1109/
ISGTEurope.2011.6162643.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Assessment of Faults in Hybrid System 13
6. Bayrak, G., & Cebeci, M. (2013). Grid connected fuel cell and PV hybrid
power generating system design with Matlab Simulink. Int. J. Hydrog., Science
Direct, 39(16), 8803-8812.https://doi.org/10.1016/j.ijhydene.2013.12.029.
7. Pathan, T., N., & Adhau, P. S., & Sable, M. (2017). MPPT for grid connected
hybrid wind driven PMSG - solar PV power generation system with single
stage converter. J. Electr. Power Syst. Eng., 3(1), 41-49.
8. Jyothi, A., S., P., Lakshmanrao, Kurian, P., C., & Singh, K., B. (2014). Design
and simulation of wind and solar hybrid system under isolated mode of
operation. Int. J. Ind. Electron. Electr. Eng., 2(3), 23-27.
9. Ren, H., Wu, Q., Gao, W., & Zhou, W. (2016). Optimal operation of a
grid-connected hybrid PV/fuel cell/battery energy system for residential
applications. Energy, Science Direct. 113, 702-712. https://doi.org/10.1016/j.
energy.2016.07.091.
10. González, A., Riba, R., J., Rius, A., & Puig, R. (2015). Optimal sizing of a hybrid
grid-connected photovoltaic and wind power system. Appl. Energy, Science
Direct. 154, 752–762. https://doi.org/10.1016/j.apenergy.2015.04.105.
11. Swain, S., & Ray, K., P. (2017). Short circuit fault analysis in a grid con-
nected DFIG based wind energy system with active crowbar protec-
tion circuit for ride through capability and power quality improvement.
Electrical Power and Energy Systems. 84, 64 –75. https://doi.org/10.1016/J.
IJEPES.2016.05.006.
12. Mukherjee, N., Chattopadhyaya, A., Chattopadhyay, S., &Sengupta, S.
(2020). Discrete Wavelet Transform and Stockwell Transform based
Statistical Parameters Estimation for Fault Analysis in Grid Connected Wind
Power System. IEEE Syst. J. 14(3), 4320 – 4328. https://doi.org/10.1109/
JSYST.2020.2984132.
13. Datta, S., Chattopadhyaya, A., Chattopadhyay, S., & Das, A. (2021). Harmonic
Distortion, inter-harmonic Group Magnitude and Discrete Wavelet Transform
based Statistical Parameter Estimation for Line to Ground Fault Analysis in
Microgrid System. MFIIS-2020, Published in IET Digital Library. 177 – 184.
https://doi.org/10.1049/icp.2021.1087.
14. Datta, S., Chattopadhyaya, A., Chattopadhyay, S., & Das, A. (2020). Line to
ground and line to line fault analysis in IEEE standard 9 bus system. AMSE-
IIETA Journal 19(504(2A)), 10-18. https://doi.org/10.18280/mmc_a.931-402.
15. Chattopadhyay, S., Chattopadhyaya, A., & Sengupta, S. (2014). Measurement
of harmonic distortion and skewness of stator current of induction motor at
crawling in clarke plane. IET Sci. Meas. Technol. 8(6), 528-536. https://doi.
org/10.1049/iet-smt.2013.0082.
16. Chattopadhyay, S., Chattopadhyaya, A., & Sengupta, S. (2014). Analysis of
stator current of induction motor used in transport system at single phas-
ing by measuring phase angle, symmetrical components, skewness, kurtosis
and harmonic distortion in Park plane. IET Electr. Syst. Transp., 4(1), 1 – 8.
https://doi.org/10.1049/iet-est.2012.0048.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
2
Diversified Harmonics Modeling for
Power System Stability Analysis
Tamal Roy1, Debopoma Kar Ray1 and Surajit Chattopadhyay2*
1
Dept. of Electrical Engineering, MCKVIE, West Bengal, India
2
Dept. of Electrical Engineering, GKCIET, West Bengal, India
Abstract
This chapter presents diversified harmonics modeling for power system networks.
The chapter starts with a description of a simplified steady-state stability model.
Network stability is classified based on load change and disturbance. Then the
presence of harmonics is introduced step by step. Then a brief discussion is pre-
sented on how the model will look in the presence of harmonics. The change of
operating points due to harmonics is discussed in a case study.
2.1 Introduction
The term ‘stability’ in the power system indicates the state of ‘synchronism’
maintained during operation. Power system networks involve synchro-
nous generators along with other types of generation units. Moreover, they
feed power to a variety of loads of diverse and dynamic characteristics.
Also, power system networks face different types of faults and other types
of disturbances. All these make the operation in synchronism a very diffi-
cult and challenging task to be maintained. The maximum power that can
be delivered to maintain stability is known as the stability limit.
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(15–32) © 2024 Scrivener Publishing LLC
15
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
16 Linear and Nonlinear System Modeling
2.2 Classification
System characteristics are greatly influenced by load condition, disturbance,
method, and capacity to maintain stable operation. Concerning these, power
system stability (PSS) is classified into three categories, as follows:
a. Steady-state stability
b. Transient stability
c. Dynamic stability
operating points (HOPs) due to fault. Also, the change of harmonic load
angle due to fault at different frequencies is observed. Based on the obser-
vation logic can be made to judge whether the system is healthy or faulty.
If there is any fault, then fault bus along with the type of fault is identified.
The specific outcome of this method is that, here, fault bus and type of fault
can also be assessed from other buses and feeders which may not be con-
nected directly with fault bus by using harmonic voltage angle and power
flow monitoring, and logic is developed accordingly. This method can be
extended for other multi-bus systems and for identification of other types
of faults. Also, steady state stability analysis can be done with the aid of this
approach [10].
EV
P sin (2.1)
X
Generator Load
Pm
O δ
EV
Pmax = (2.2)
X
The above expression shows that maximum power increase with the
increase of generated and load voltage and decreases with the increase of
system reactance.
installed in the system. It may be noted that there are other major sources
of harmonics in the network like converters, faults, transients, and nonlin-
ear loads. Thus, all these harmonics produced in the system also introduce
nonlinearity in the network.
Vj = voltage at bus j
Vk = voltage at bus k
δjk = difference between voltage of bus j and bus k
X = reactance between bus j and bus k
V jVk
P sin jk (2.3)
X
Thus, the power flow of the network and its stability will depend on the
reactance of the network. Network reactance can be expressed in matrix
form. Let,
The reactance matrix can be derived from the admittance matrix of the
network. Let,
[ X ] [Y ] 1 (2.4)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Diversified Harmonics Modeling for PSS Analysis 21
Y11 Y1n
[Y ]n n Y21 Y2n (2.6)
Y31 Ynn
X R j L (2.7)
X j L (2.8)
|X| |j L| L (2.9)
1 1
Y (2.10)
X j L
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
22 Linear and Nonlinear System Modeling
1 1
|Y |
j L L
(2.11)
|X| |j L| L 2 f L (2.12)
where,
2 f (2.13)
The system is expected to operate with 5% (it may vary based on regula-
tion guidelines) deviation at the following frequency:
f 50 Hz or 60 Hz (2.14)
i 2 fi (2.15)
where,
i 1, 2, 3, 4, 5.n
i 1, 3, 5.n (2.16)
i 1, 2, 3, 5.n (2.17)
i = fractional (2.18)
|Xi | |j i L| i L (2.19)
V j i Vk i
Pjk i sin jk i (2.20)
X jk i
V j i Vk i
Pjk i max (2.21)
X jk i
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
24 Linear and Nonlinear System Modeling
Pi
Pim
Pci
O δi
Figure 2.3 Harmonic power, operating point, and corresponding maximum power.
i =1 (2.22)
V j 1Vk 1
Pjk 1 max (2.23)
X jk 1
P1
Pm1
P1c
O δ1
Figure 2.4 Fundamental power, operating point, and corresponding maximum power.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Diversified Harmonics Modeling for PSS Analysis 25
Line
Figure 2.5 Two-bus network model used for harmonic stability study.
100
90
80
Bus voltages (V)
70
60
50
40
30
20
10
0
(10Hz) (20Hz) (30Hz) (40Hz) (70Hz) (80Hz) (90Hz) (150Hz)
Frequency components (Hz)
Bus 1 Bus 2
Figure 2.6 Harmonic components of bus voltages excluding fundamental component for
300-kW load.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
28 Linear and Nonlinear System Modeling
90
85
80
75
70
(50Hz) (10Hz) (20Hz) (30Hz) (40Hz) (70Hz) (80Hz) (90Hz) (150Hz)
90
80
70
Bus voltages (V)
60
50
40
30
20
10
0
(10Hz) (20Hz) (30Hz) (40Hz) (50Hz) (70Hz) (80Hz) (90Hz) (150Hz)(250Hz)(350Hz)(450Hz)(550Hz)
Frequency components (Hz)
Bus 1 Bus 2
Figure 2.8 Harmonic components of bus voltages excluding fundamental component for
150-kW load.
Observation:
The major observations from the case study are as follows:
2.10 Conclusions
Different issues of PSS concerning harmonic power flow in the network
were presented in this chapter. At first, system stability was defined.
Considering gradual load and large disturbances, power system stability
was classified. The power equation needed for stability study was expressed.
Different useful parameters were pointed out that are important for stabil-
ity consideration. Then, a detailed discussion was made on how harmonics
influence system reactance, operating points that are important for stabil-
ity analysis. Then the case study results were provided, which show that
operating points of different frequency components are different. As the
magnitudes of harmonic components are significantly less than that of the
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Diversified Harmonics Modeling for PSS Analysis 31
References
1. C. H. Joe, J. Juan, G. Sanchez, “Flexible AC Transmission Systems,” in Power
System Modeling, Computation, and Control, IEEE, 2020, pp.421-485, DOI:
10.1002/9781119546924.ch14.
2. J. L. R. Torres, F. G. Longatt, “Implementation of a Real Phasor Based Vulnerability
Assessment and Control Scheme: The Ecuadorian WAMPAC System,” in
Dynamic Vulnerability Assessment and Intelligent Control: For Sustainable Power
Systems, IEEE, 2018, pp.389-411, DOI: 10.1002/9781119214984.ch18.
3. D. Xu, F. Blaabjerg, C. Wenjie, N. Zhu, “Grid Codes for Wind Power
Generation Systems,” in Advanced Control of Doubly Fed Induction Generator
for Wind Power Systems, IEEE, 2018, pp.43-63, DOI: 10.1002/9781119172093.
ch3.
4. D. Xu, F. Blaabjerg, W. Chen, N. Zhu, “Resonant Control of DFIG Under
Grid Voltage Harmonics Distortion,” in Advanced Control of Doubly Fed
Induction Generator for Wind Power Systems, IEEE, 2018, pp.195-235, DOI:
10.1002/9781119172093.ch8.
5. H. Farhangi, G. Joos, “Analysis and Studies Using Recommended Models,” in
Microgrid Planning and Design: A Concise Guide, IEEE, 2019, pp.57-63, DOI:
10.1002/9781119453550.ch4.
6. A. A. Sallam, O. P. Malik, “Distributed Generation,” in Electric Distribution
Systems, IEEE, 2019, pp.487-533, DOI: 10.1002/9781119509332.ch18.
7. S Chattopadhyay, “Nanogrids and Picogrods and their integration with electric
vehicles”, IET, London, ISBN: 978-1-83953-482-9, 2022.
8. S Chattopadhyay, A Das, “Overhead Electric Power Lines: Theory and prac-
tice”, IET, London, ISBN: 9781839533112, 2021.
9. S Chattopadhyay, M Mitra, S Sengupta, “Electric Power Quality”, Springer,
Netherland, ISBN: 978-94-007-0635-4, 2011.
10. D. Kar Ray (2019). “Harmonic Power Flow Analysis and its Applications in
Electrical Systems”, Ph. D. Thesis, UCSTA, CU, Kolkata, India.
11. T. Zabaiou, L. A. Dessaint, I. Kamwa, “Preventive Control Approach for
Voltage Stability Improvement using Voltage Stability Constrained Optimal
Power Flow based on Static Line Voltage Stability Indices”, IET Gen. Trans.
Distr., Vol. 8 (5), pp. 924-934, doi:10.1049/iet-gtd. 2013. 0724, 2014.
12. P. M. Anderson, A. A. Fouad, “Power System Control and Stability”, IEEE
Press, ISBN (Print): 978-0-471-23862-1, 2003.
13. P. W. Sauer, M. A. Pai, “Power System Dynamics and Stability”, Stipes
Publishing Co., ISBN-10:1588746739, ISBN-13:978-1588746733, 2016.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
32 Linear and Nonlinear System Modeling
2
Electrical Engineering Department, GNIT, Panihati, Kolkata, India
3
Electrical Engineering Department, G.K.C. Institute of Engineering and Technology,
Malda, India
4
Electrical Engineering Department, Jadavpur University, Kolkata, India
Abstract
Reliability has become one of the most important properties of a power distribu-
tion network in modern era. Nowadays, a distribution network faces threats of
cyberattacks on it along with weather-related events, such as storms and floods.
Microgrids gained extreme popularity for enhancing overall reliability and resil-
ience of conventional power networks. It consists of a set of distributed energy
sources and loads present in a small region, which acts as a separate controllable
set-up in relation to the main grid. However, widespread introduction of microg-
rid introduces some new technical challenges. Unlike traditional power networks,
generators, energy storage facilities, and loads are present at all the levels. Thus,
generation, distribution, and consumption can occur at any level, making the
direction of power flow very much uncertain. Deep understanding of a microgrid
operation is still at its infancy. In order to develop a better understanding, stan-
dardized forms of microgrid networks are of paramount importance. This chapter
showcases three most widely used standardized networks: IEEE 14 bus microgrid
network, IEEE 9 bus microgrid network, and IEC 61850-7-420 standard microg-
rid network.
Keywords: DC microgrid, AC microgrid, hybrid microgrid, energy storage,
energy resource, power distribution, protection schemes
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(33–48) © 2024 Scrivener Publishing LLC
33
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
34 Linear and Nonlinear System Modeling
3.1 Introduction
In recent times, reliability, i.e., the ability to stay ‘operational in the
event of adversity’ is becoming one of the most desirable properties of
a power distribution network. The word ‘adversity’ refers to weather-
related events, such as storms and floods, and also physical/cyberattacks
on power network. Evolutions in the area of distributed energy resource
(DER) and energy storage system are making microgrids extremely pop-
ular for enhancing overall reliability and resilience of conventional power
networks. Microgrid consists of a set of distributed energy sources and
loads present in a small area, which acts as a separate controllable set-up
with respect to the main grid, having a distinct electrical borderline. It
is capable of working either connected with the traditional power grid
or on its own as an island. However, widespread introduction of micro-
grid will bring forth its own sets of technical challenges. The main chal-
lenges deal with the significant impact on existing networks created due
to the deployment of a number of DERs and design of proper protection
schemes which can operate satisfactorily in the much lower voltage and
current levels present in the microgrid environment. Unlike traditional
power networks, generators, energy storage facilities, and loads are pres-
ent at every level of a microgrid network. Generator, distributor, and also
consumer can stay at any level of a microgrid network and the direction
of power flow is also uncertain. Complete understanding of a micro-
grid operation is still at a very early phase. Thus, to facilitate a secure
operation in such a multifaceted environment, where so many different
parameters and variables are to be monitored, controlled and communi-
cated, extensive research work is needed. In order to do so, standardized
forms of microgrid networks are of paramount importance. This chapter
showcases three most widely used standardized networks: IEEE 14 bus
microgrid network, IEEE 9 bus microgrid network, and IEC 61850-7-420
standard microgrid network.
PV Wind Battery
Grid Farm Farm Storage
DC Bus
PV Wind
Grid Farm Farm
AC Bus
Transformer Transformer
AC Loads AC Loads
from it and the main conventional grid can also get the excess power of
the microgrid. Different voltage levels can be achieved through the use
of transformers. Perfect synchronization of phase angle and frequency is
needed for the operation of AC microgrids. High power applications need
this type of microgrids.
AC Bus
Rectifier
DC Bus
Battery
Storage PV Farm DC Loads HEV
T1
SS 1 SS 2
T3
T2 Yes T4
T2
No Blackout
Success?
Bulk
PV Cell
Generator
Bus B2 Bus B1
Neutral
Grounding
Furnace
Bus B4 Bus B3
Non-linear
Load
Diesel
Generator 1
Battery
Charging Bus B5
System
Linear
Bus B6 Load 1
Bus B7
Bus B8
Diesel
Linear Generator 2
Load 2
Bus B9
has an arrangement for battery charging, B2 has linear load, and B3 has
linear load 1 connected to it. The length of transmission lines connecting
one bus to its adjacent bus is 1 km. Operating voltage and frequency is
kept at 400 V (line to line) and 50 Hz, respectively. Bulk generator bus is
considered as the swing bus. Bulk generator is rated at 100 MW, PV system
has a capacity of 40 MW, and both the diesel generators are rated at 50 MW.
Bus 8 Bus 9
Load C T3 G2
G2 T2 Line 7-8 Line 8-9
Wind
farm A Bus 3
Bus 2
Line 5-7
Line 6-9
Bus 7 PV
farm A
Bus 5 Bus 6
Line 4-5
Line 4-6
Bus 4
BESS T1
Bus 1
G1
Utility Grid
CB1
CB2
CB3 CB4
Load2 DG3
Load3 DG4 DG5
Microgrid
3.8 Conclusion
The presence of sources like small PV panels, wind turbines, and diesel
generators in a microgrid network grid brings forth marked differences
between the structure of a traditional power grid and a microgrid. These
radical changes in structure create challenges in operation and protection
of a microgrid network which were not known by the power engineers.
Traditional load flow control methods have also been rendered ineffective
in a microgrid environment along with the popular protection schemes.
Hence, rigorous study and research work is needed to address the vari-
ous issues faced in a microgrid environment. The standardized microgrid
networks play an important role in studying the nature of microgrids in
different modes of operations as well as devising new operation and pro-
tection strategies.
References
1. S. Datta, A. Chattopadhyaya, S. Chattopadhyay & A. Das (2020), “S-Transform
Based Kurtosis Analysis for Detection of LG and LL Faults in 14 Bus Microgrid
System,” IETE J. Res., Vol. 2, no. 1, pp. 10_20.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
46 Linear and Nonlinear System Modeling
2. U. Datta, A. Kalam, & J. Shi (2020), “Battery Energy Storage System for
Aggregated Inertia-Droop Control and a Novel Frequency Dependent State-
of-Charge Recovery”, Energies, Vol. 13, no. 8, pp. 1_18.
3. T. S. Ustun, C. Ozansoy, & A. Zayegh (2011), “Distributed Energy Resources
(DER) object modeling with IEC 61850–7–420,” Australasian Universities
Power Engineering Conference (AUPEC 2011), Brisbane, Australia, pp. 1_6.
4. M. Mishra, M. Sahani, & P.K.Rout (2017), “An islanding detection algorithm
for distributed generation based on Hilbert–Huang transform and extreme
learning machine”, Sustain. Energy Grids Netw., vol. 9, no. 1, pp. 13_ 26.
5. M. Manohar, E. Koley, & S. Ghosh (2017), “Reliable protection scheme for
PV integrated microgrid using an ensemble classifier approach with real-
time validation”, IET Sci. Meas. Technol., vol. 12, no. 2, pp. 200_208.
6. S. Datta, S. Chattopadhyay and A. Das (2017), “FFT based Harmonic
Assessment of Line to Ground Fault in 14 Bus Microgrid System,” 2017
International Conference on Modelling and Simulation in Science, Technology
and Engineering Mathematics (MS-17), Kolkata, India, pp. 73_86.
7. N. Mukherjee, A. Chattopadhyaya, S. Chattopadhyay, & S. Sengupta (2020),
“Discrete-Wavelet-Transform and Stockwell-Transform-Based Statistical
Parameters Estimation for Fault Analysis in Grid-Connected Wind Power
System”, IEEE Syst. J., vol. 14, no. 3, pp. 4320_4328.
8. T. S. Ustun, C. Ozansoy, & A. Zayegh (2012), “Modeling of a Centralized
Microgrid Protection System and Distributed Energy Resources According
to IEC 61850-7-420”, IEEE Trans. Power Syst., Vol. 27, no. 3, pp. 1560_1567.
9. S. Dhar, R. K. Patnaik, & P. K. Dash (2017), “Fault Detection and Location
of Photovoltaic Based DC Microgrid Using Differential Protection Strategy”,
IEEE Trans. Smart Grid, vol. 9, no. 5, pp. 4303_4312.
10. T. Ghanbari, & E. Farjah (2013), “Unidirectional Fault Current Limiter: An
Efficient Interface between the Microgrid and Main Network”, IEEE Trans.
Power Syst., vol. 28, no. 2, pp. 1591_1598.
11. T. S. Ustun, C. Ozansoy, & A. Zayegh (2011), “Extending IEC 61850-7-420
for distributed generators with fault current limiters”, IEEE PES Innov. Smart
Grid Technol., Perth, Australia, pp. 1_8.
12. B. M. Eid, N. A. Rahim, J. Selvaraj, & A. H. El. Khateb (2016), “Control
Methods and Objectives for Electronically Coupled Distributed Energy
Resources in Microgrids”, IEEE Syst. J., Vol. 10, no. 2, pp. 446_452.
13. D. Rebollal, M. Carpintero-Rentería, D. Santos-Martín, & M. Chinchilla (2021),
“Microgrid and Distributed Energy Resources Standards and Guidelines
Review: Grid Connection and Operation Technical Requirements,” Energies,
Vol. 14, no. 3, pp. 1_24.
14. Y. Usama, & X. Lu (2014), ‘Design and implementation of a wavelet analy-
sis-based shunt fault detection and identification module for transmission
lines application’, IET Gener. Transm. Distrib., Vol. 8, no. 3, pp. 431_441.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Comparative Study of Microgrid Networks 47
15. D. Guillen, M. Paternina, & J. Ortiz-Bejar (2018), ‘Fault detection and clas-
sification in transmission lines based on a PSD index’, IET Gener. Transm.
Distrib., Vol. 12, no. 18, pp. 4070_4078.
16. A. Ali, Q. Khan, & A. Hussain (2015), ‘Fault modelling and detection
in power generation, transmission and distribution systems’, IET Gener.
Transm. Distrib., Vol. 16, no. 3, pp. 2782_2791.
17. D. Salomonsson, L. Soder, & A. Sannino, ‘Protection of Low-Voltage DC
Microgrids’, IEEE Trans. Power Deliv., Vol. 24, pp. 1045_1053.
18. S. Wen, Y. Wang, & Y. Tang (2019), ‘Proactive frequency Control based
on Ultra-Short-Term Power Fluctuation Forecasting for High Renewable
Penetrated Power Systems’, IET Renew. Power Gener., Vol. 13, no. 12, pp.
2166_2173.
19. J. Chen, C. Kuo, & S. Chen (2019), ‘DC-side fault detection for photovoltaic
energy conversion system using fractional-order dynamic-error-based fuzzy
Petri net integrated with intelligent meters’, IET Renew. Power Gener., Vol. 10,
no. 9, pp. 1318_1327.
20. Y. M. Makwana, & R. B. Bhalja (2017), ‘Islanding detection technique based
on superimposed components of voltage’, IET Renew. Power Gener., Vol. 11,
no.11, pp. 1371_1381.
21. M. Manohar, E. Koley, & S. Ghosh (2019), ‘Enhancing the reliability of pro-
tection scheme for PV integrated microgrid by discriminating between array
faults and symmetrical line faults using sparse auto encoder’, IET Renew.
Power Gener., Vol. 13, no. 2, pp. 308_317.
22. R. Mohanty, U. S. M. Balaji, & A. K. Pradhan (2016), ‘An Accurate Non-
iterative Fault-Location Technique for Low-Voltage DC Microgrid’, IEEE
Trans. Power Deliv., Vol. 31, no. 2, pp. 475_481.
23. W. Qiu, Q. Tang, J. Liu, Z. Teng, & W. Yao (2019), ‘Power Quality Disturbances
Recognition Using Modified S – Transform and Parallel Stack Sparse Auto-
encoder’, Electr. Power Syst. Res., Vol. 174, no. 1, pp. 1_10.
24. C. Lin, W. Gao, & M. F. Guo (2019), ‘Discrete Wavelet Transform Based
Triggering Method for Single Phase Earth Fault in Power Distribution
Systems’, IEEE Trans. Power Deliv., Vol. 34, no. 5, pp. 2058_2068.
25. R. Kumar, B.Singh, D. T. Shahani, A. Chandra, & K. Al-Haddad (2015),
‘Recognition of Power Quality Disturbances Using S – Transform Based
ANN Classifier and Rule Based Decision Tree’, IEEE Trans. Ind. Appl., Vol.
51, no. 2, pp. 1249_1258.
26. M. Ikhide, S. B. Tennakoon, H. Ha, A. L. Griffiths, S. Subramanian, & A. J.
Adamczyk (2019), ‘A transient and non-unit-based protection technique for
DC grids based on the rate-of-change (R-o-C) of the fault induced travelling
wave components’, Sustain. Energy Grids Netw., Vol. 17, no. 1, pp. 100_11.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
4
Application of Active Power Filter
in the Hybrid Power System to
Regulate the Grid Voltage
Sarita Samal1*, Rudranarayan Dash1, Arjyadhara Pradhan1
and Prasanta Kumar Barik2
2
Dept. of Mechanical &Electrical Engineering, CAET, OUAT,
Bhubaneshwar, Odisha, India
Abstract
This work outlines a strategy for integrating solar photovoltaic (SPV) and wind
energy system (WES) connections. Both rely on storage batteries for energy deliv-
ery, which are one of the distribution system’s key advantages. This technique is
used to control the voltage at the point of common coupling (PCC) while the SPV
system is not producing any power. Using the SPV and batteries as storage energy,
the overall system was able to achieve bidirectional grid voltage regulation. During
voltage disruptions like sag, swell, and interruption, the bidirectional inverter
functions as a normal active power filter, regulating the grid voltage. When the
SPV system is using a battery charger, particularly at night time to charge the bat-
teries, it will behave like a conventional switch, and energy will be stored. The pro-
posed technique established the idea of indirect feeder link pathways to regulate
voltage regulation, in which as the voltage level rises, a large amount of reverse
power flow from WES is automatically regulated by employing the SPV system
inverter to charge and store the power. Simulink blocks in MATLAB can be used
to model the suggested system. The proposed strategy to adjust the feeder voltage
by exchanging real active power through the energy storage batteries could be
verified by simulation results. Further, the proposed system may be tested with
different recent control approaches of active filters for better performance.
Keywords: Active power filter, solar PV, wind energy, voltage sag, swell
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(49–66) © 2024 Scrivener Publishing LLC
49
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
50 Linear and Nonlinear System Modeling
4.1 Introduction
The performance of the electrical power generation has greatly deterio-
rated recently as a consequence of increased power utilization. The use
of nonlinear loads in the system results in harmonics, which affect the
current and voltage waveforms and reduce the power factor of the supply
system [1]. Consequently, nonlinear loads in the system must be regulated
for harmonics. Because passive filters are inexpensive and reliable, they
have been employed for harmonic reduction in the past, but they have
several problems, including resonance and mistuning, which cause their
interruption in the network [2, 3]. The passive filters, have the drawbacks
of being costly and having large inverter losses [4, 5]. As a consequence,
this chapter presents a series active power filter (APF) incorporated with
a hybrid power system in conjunction with multiple control approaches to
eliminate all of these shortcomings and produce superior harmonic mit-
igation at a lower overall cost [6, 7]. The authors of References [8] and
[9] examined the dynamics of a combined series APF and battery energy
storage (BES) combined arrangement to reduce resonance issues during
transient conditions. A PI (proportional integration) control was used to
drive the voltage source control (VSC). The authors of References [10]
and [11] projected a PI-controlled series APF with a combined storage of
a BES and a super capacitor to increase voltage stability. In References [12]
and [13], an enhanced modular multilevel converter method for dynamic
power allocation based on BES units in a PI-controlled with combined
shunt active power filter (SAPF) and BES system have been implemented.
In the majority of investigations, the combination of series APF and BES
assistance with power systems uses conventional PI control scheme for the
VSC [14, 15]. Hence, in this chapter, a combined effect of solar PV, wind
energy, and BS integrated series APF is presented and its effects are veri-
fied concerning conventional APF. Using the SPV and batteries as storage
energy, the overall system was able to achieve bidirectional grid voltage
regulation. During voltage disruptions like sag, swell, and interruption,
the bidirectional inverter functions as a normal APF, regulating the grid
voltage. When the SPV system is using a battery charger, particularly at
night time to charge the batteries, it will behave like a conventional switch,
and energy will be stored. The proposed technique established the idea of
indirect feeder link pathways to regulate voltage regulation, in which as
the voltage level rises, a large amount of reverse power flow from WES is
automatically regulated by employing the SPV system inverter to charge
and store the power.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Application of APF in the Hybrid Power System 51
AC PCC
Series Transformer
AC source 3-phase Transformer
Nonlinear load
Buck-Boost
DC-link
Converter
capacitor
DC BUS
Series APF
AC/DC DC/DC
converter converter
Hybrid energy system
+ –
Battery
Wind-system Solar PV
q Vpv I pv Rse
I pv N p I ph N p I 0 exp 1 (4.1)
N s AkT
Rs
Id
Rp
IP
point tracking. Perturb and observe (P&O) provides a simple and quick
solution to changes in solar radiation. Because of this, peak energy is pro-
duced via P&O regardless of the amount of irradiation. Figure 4.3 displays
a flowchart for the P&O MPPT. The SPV system’s parameters are listed in
Table 4.1 [18].
Start
P(n)–P(n–1) YES
=0
NO
NO P(n)–P(n–1) YES
>0
NO YES NO YES
V(n)–V(n–1) V(n)–V(n–1)
<0 >0
D = D + ΔD D = D – ΔD D = D + ΔD D = D – ΔD
RETURN
L1
D1
Vin D Sw V0
C R
L
R
Tr
Diode
DC Output
25 C1
SPV
DC
400 C2
Insolation
P&O
MPPT
Wind Blade
Connected
Storage system DC converter load
Discrete,
Wp [po]
Goto5
s = 2e-05 s
Te powergui
Tm
[Idc]
Speed3 Ws Wm T Duty Cycle
Goto1
+ –i
12 Manual Switch A A V+
+ + Idc Vdc
B B Ld Goto3 v+– Cd
+
wind speed [Vdc] V–
–
C C Boost Converter
with Load
wind generation Universal Bridge
Vdc Vdc
Duty cycle
Idc Idc
MPPT controller
1
P0 CP , R 2V 3 (4.2)
2
Where P0 = mechanical power, β = pitch angle, λ = tip-speed ratio, λ
= ΩR/V, R = blade radius, and ρ = air density. These parameters are used
to model and develop the WEC system. Figure 4.6 shows the WES model
and Figure 4.7 displays the simulation model of WES. The symbols used in
WES are listed in Table 4.2.
Internal
Resistance IBattery
Controlled VBattery
Voltage Source
VB Integrator
+ VL – I
IS S IO
IC +
R
+
VS VO
– C
D –
Operating Principle
The input inductor, which results in an unanticipated shift in the input cur-
rent, is what allows the DC to DC converter to work. The inductor absorbs
energy from the input and stores it when the switch is switched on. The
energy is released when the switch is turned off.
Figure 4.9 shows how a buck converter works. The buck converter’s
metal oxide silicon field effect transistor (MOSFET) is switched on and the
subsequent semiconductor is disabled as a result of the increased square
wave frequency. The current passes through the magnetic field, charging
C, and powers the load if the magnetic field’s current is larger than the first
semiconductor’s gate terminal’s current. The positive voltage applied to the
cathode disables the Schottky diode D1.
L1 D1
+ +
d
DC + DC
Input Output
VIN C1 Load VOUT > VIN
g
s
Square wave MOSFET
driving switch Switch
800
700
Voltage in (V) 600
500
400
300
200
100
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
a ia vca
vsa ib va
Electrical b vcb
Power Grid vsb vb Load
ic vcc
c
vsc vc
N
vsa vca*
vsb
vsc Controller vcb* Inverter
ia
ib vcc*
ic
Vdc + –
in the dq0 reference plane after being translated from abc coordinates to
dq0 coordinates. After comparison, both are transformed back into ABC
reference frames. The phase locked loop (PLL) required for Park’s trans-
formation and inverse Park’s transformation can be produced. The series
converter’s six switches are linked to a pulse width modulation (PWM)
generator, which provides synchronized switching pulses. By providing
proper switching pulses, the series converter’s goals can be met. The gen-
eration of the switching pulses is obtained from the three-phase reference
signal that is acquired from two controllers. The pulses required for VSI
operation are produced by a hysteresis voltage controller by comparing the
expected voltage of the output (Vc*) to the real series APF output voltage
(Vc). The modeling of series APF control scheme is shown in Figure 4.13.
The series APF detects a voltage drop in a sag condition and injects
the necessary voltage through the series coupling transformer. It ensures
Lowpass
Lowpass Filter
abc
0
dq0 dq0
Constant2
sin_cos Terminator1 abc 1
Freq sin_cos Vabc1
Terminator2 abc_to_dq0 0
1 Vabc (pu) wt dq0_to_abc
Transformation
Vabc Terminator Constant1 Transformation
Sin_Cos
3-phase PLL
Terminator3
Vdc
PLL
Vdc*
PI
controller
S1
Isa
dq abc Vsa
K Vsa* S2
VLa
LPF controller
Vsb S3
Isb Vsb* S4
K controller
S5
VLb Vsc
Vsc* S6
Isc abc LPF dq controller
K
VLc
that the load terminal is kept at the rated voltage. The series APF recog-
nizes a voltage swell as it happens and injects the compensating voltage to
cancel it out, maintaining the evaluated voltage across the load terminal.
This method is equivalent to maintaining the rated voltage under uneven
source voltage. After the series inverter had injected the compensation
voltage in series with the source voltage, it was discovered that the load
voltage was similar to the combination of the source and replacement volt-
ages. When the source voltage is non-sinusoidal and unbalanced, the load
bus’s rated sinusoidal voltage is controlled by injecting a series voltage at
the necessary magnitude and phase angle. The system is therefore safe-
guarded against voltage-related PQ problems including voltage sag, swell,
imbalance, flicker, and so forth. Figure 4.14 demonstrates the function of
the PI controller.
Sag Compensation
To validate the series APF operation, voltage sag is injected into the system.
The series APF is employed at the load terminal to address voltage-related
power quality (PQ) abnormalities such voltage sag, swell, inequalities, flicker
notches, and so on. The load voltage profile exhibited in Figure 4.15(a)
matches the voltage sag observed on the time axis from 0.1 s to 0.3 s. A volt-
age decrease is detected by the series APF during a sag state, and it injects the
necessary voltage through the series coupling transformer. It ensures that the
load terminal is kept at the rated voltage. As seen in Figure 4.15(b), when
the APF is turned on, it injects compensating voltage into the PCC to make
up for load voltage sag. The voltages at the source and the load are there-
fore equivalent. The voltage of the load following the adjustment is shown
in Figure 4.15(c). In order to keep the rated voltage across the load terminal
constant, the series APF operates by quickly detecting voltage changes at the
source and injecting a compensation voltage. The voltage on the DC-link
capacitor is revealed in Figure 4.15(d). The capacitor voltage is practically
constant up to 700 V with a settling period of 0.25 s, as shown in the result.
Sag compensation
To test the effectiveness of the series APF, voltage sag is injected into the
system for 0.1 s to 0.3 s of the voltage waveform, shown in Figure 4.16(a).
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
62 Linear and Nonlinear System Modeling
600
400
Voltage in (V)
200
–200
–400
–600
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
(a)
100
Voltage in (V)
50
–50
–100
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
(b)
600
400
Voltage in (V)
200
0
–200
–400
–600
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
(c)
800
700
600
Voltage in (V)
500
400
300
200
100
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
(d)
Figure 4.15 Profile obtained Case 1. (a) Load voltage before compensation, (b) compensating
voltage injected by APF, (c) load voltage after compensation, and (d) DC-link capacitor
voltage.
600
400
Voltage in (V) 200
–200
–400
–600
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
(a)
100
Voltage in (V)
50
0
–50
–100
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
(b)
600
400
Voltage in (V)
200
0
–200
–400
–600
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Time in (S)
(c)
800
700
600
Voltage in (V)
500
400
300
200
100
0
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Time in (S)
(d)
Figure 4.16 Profile obtained Case 2. (a) Load voltage before compensation, (b) compensating
voltage injected by APF, (c) load voltage after compensation, and (d) DC-link capacitor
voltage.
capacitor. The result demonstrates that the capacitor voltage is nearly con-
stant up to 700 V with a settling time of 0.025 s.
The suggested hybrid energy system with series APF can adjust for
current disturbances at the load side while also maintaining a consistent
DC-link capacitor voltage. Nonlinear load is considered when evaluating
the suggested APF’s performance. The load voltage sag can be efficiently
corrected, and the settling time of the DC-link capacitor voltage is reduced
to 0.025 s, as shown by the numerical comparison of results.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
64 Linear and Nonlinear System Modeling
4.5 Conclusion
This work focuses on power quality improvement of power system with
nonlinear load using series APF. The performance of the APF is studied
with and without integrating a hybrid energy system across its DC-link
side. Moreover, the performance of APF is improved by regulating the
DC-link capacitor voltage by integrating hybrid renewable energy-based
sources based across it. From the obtained results, it is found that load
voltage sag can be compensated effectively and also, the settling time of
DC-link capacitor voltage is reduced to 0.025 s. During voltage disruptions
like sag, swell and interruption, the bidirectional inverter functions as a
normal APF, regulating the grid voltage. When the SPV system is using
a battery charger, particularly at night time to charge the batteries, it will
behave like a conventional switch, and energy will be stored.
References
1. Z. Zeng, H. Yang, S. Tang, and R. Zhao, “Objective-oriented power quality
compensation of multifunctional grid-tied inverters and its application in
microgrids,” IEEE Trans. Power Electron., vol. 30, no. 3, pp. 1255–1265, 2015.
2. H. Akagi and H. Fujita, “A new power line conditioner for harmonic com-
pensation in power systems,” IEEE Trans. Power Deliv., vol. 10, no. 3, pp.
1570–1575, 1995.
3. S. Rahmani, K. Al-Haddad, and F. Fnaiech, “A series hybrid power filter to
compensate harmonic currents and voltages,” in Proceedings of the IEEE 2002
28th Annual Conference IECON 02 Industrial Electronics Society, pp. 644–
649, Sevilla, Spain, November 2002
4. E. R. Ribeiro and I. Barbi, “Harmonic voltage reduction using a series active
filter under difffferent load conditions” IEEE Trans. Power Electron., vol. 21,
no. 5, pp. 1394–1402, 2006.
5. F. Z. Peng, H. Akagi, and A. Nabae, “A new approach to harmonic compensa-
tion in power systems-a combined system of shunt passive and series active
fifilters,” IEEE Trans. Ind. Appl., vol. 26, no. 6, pp. 983–990, 1990.
6. S. Janpong, K. L. Areerak, and K. N. Areerak, “A literaturesurvey of neural
network applications for shunt active power filters,” Int. J. Electr. Comput.
Eng., vol. 5, no. 12, 2011.
7. B. Singh, G. Bhuvaneswari, and S. R. Arya, “Review on powerquality solution
technology,” Asian Power Electron. J., vol. 6, no. 2, 2012.
8. C. Kumar and M. K. Mishra, “An improved hybrid DSTATCOM topology to
compensate reactive and nonlinear loads,” IEEE Trans. Ind. Electron., vol. 61,
no. 12, pp. 6517–6527, 2014
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Application of APF in the Hybrid Power System 65
Abstract
This chapter focuses on dynamic modeling of drone using MATLAB. The control
of diverse drone movements is vital for understanding and managing drone flight
dynamics. The main objective is to replicate movements like ascending, descend-
ing, forward-backward, and lateral translations through a comprehensive simula-
tion model within MATLAB. This simulation will accurately represent the drone’s
flight dynamics and interactions within its environment. The understanding of
drone maneuverability and control is serving as a foundational step towards devel-
oping more sophisticated control systems for these aerial platforms by facilitating
advancements in drone technology across various industries.
5.1 Introduction
Drones, also known as unmanned aerial vehicles (UAVs), have become
increasingly popular in recent years due to their versatility and usefulness
in a wide range of industries (Shelare et al., 2021) [1]. These aerial vehi-
cles are equipped with advanced technology and can be used for a variety
of purposes, including surveillance, reconnaissance, agriculture, disaster
management, and filmmaking. Drones are particularly valuable because
Email: tri.suman78@gmail.com
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(67–88) © 2024 Scrivener Publishing LLC
67
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
68 Linear and Nonlinear System Modeling
they can access remote or dangerous areas, collect real-time data, and per-
form tasks with precision and efficiency.
Understanding the complex dynamics of drone flight is crucial for
optimizing their performance and ensuring safe and reliable operation.
Drones can move in six directions (six degrees of freedom or 6-DOF),
which means they can travel up, down, forward, back, left, and right, and
they can also rotate in all three dimensions (roll, pitch, and yaw) (Najm
and Ibraheem, 2019) [1]. Comprehending these 6-DOF is essential for
effective drone control because it directly affects how we can manage the
drone’s position, orientation, and overall flight dynamics. This under-
standing is the foundation for developing precise control strategies,
enhancing operational capabilities, and navigating the challenges inher-
ent in drone flight.
MATLAB’s Simulink is a powerful platform for modeling, simulating,
and analyzing dynamic systems, making it an ideal tool for understanding
the complexities of drone flight. This chapter aims to provide a compre-
hensive understanding of the parameters governing drone flight, enabling
proficient navigation of the MATLAB and Simulink environment.
5.3 Methodology
In the process of designing a drone control system using MATLAB, the
understanding of the 6-DOF in drone flight (thrust or up/down, pitch or
forward/back, roll or left/right, and yaw or twist left and right) is essential
for applying effective engineering. This knowledge serves as the founda-
tion for developing a control system that enables the drone to hover and
move seamlessly.
The design of a control system for a drone involves a thorough grasp
of hardware components and sensors, such as ultrasound, camera, pres-
sure sensor, and inertial measurement unit (IMU). Addressing the under-
actuated nature of the system, with four actuators for 6-DOF, presents a
challenge. The configuration and spin directions of motors are pivotal,
requiring a motor mixing algorithm to convert intuitive commands into
precise motor speeds. The final control system, built and simulated in
Simulink, undergoes tuning and testing before generating flight code for
the onboard microcontroller. This methodology leverages sensor data for
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
70 Linear and Nonlinear System Modeling
state estimation and deploys algorithms for precise motor control, ensur-
ing desired drone movements in three-dimensional (3D) space.
actuators, precisely governing the motor actions. Meanwhile, the flag acts
as a safety measure, ensuring the drone’s flight remains stable and secure,
particularly in unforeseen circumstances.
of pitch, roll, and yaw on its orientation, we aim to gain insights into the
effectiveness and responsiveness of our drone control system.
Simulation
Values
–0.3
–2
–2
0
0
Time
100 s
40 s
70 s
0s
5s
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Dynamic Modeling of Drone Control 75
Table 5.2 Varying the values of pitch parameter of the drone flight control.
Time Values Simulation
0s 0
20 s 0.002
40 s 0.002
60 s –0.002
(Continued)
in the table with respect to the values of pitch parameters. The observed
drone movement in the X-axis coordinate, corresponding to changes in
the pitch parameter values, provides valuable insights into the impact of
pitch adjustments on the drone’s spatial dynamics. Initially set at 0, the
drone remains stationary at the origin. However, as the pitch parameter
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Dynamic Modeling of Drone Control 77
Table 5.2 Varying the values of pitch parameter of the drone flight control.
(Continued)
Time Values Simulation
80 s –0.002
100 s 0
and avoid obstacles. Without roll control, the drone would fly in a straight
line, severely limiting its ability to navigate and perform complex maneu-
vers. Figure 5.4 shows the drone flight control model for roll parameter.
The observed variations in the roll parameter values of the drone’s flight
control, as outlined in Table 5.3, indicate significant effects on the drone’s
lateral movement along the Y-axis coordinate. Initially, when the roll
parameter is set to 0, the drone remains stationary at the origin. However,
as the roll parameter increases to 0.002, a positive value, the drone exhib-
its a consistent lateral movement to the right over the course of 40 s. This
behavior is in line with the fundamental principle of drone flight control,
where adjustments in the roll parameter contribute to the drone’s angular
rotation around its longitudinal axis, consequently influencing its lateral
displacement.
Conversely, when the roll parameter transitions from 0.002 to –0.002, a
negative value, the drone’s motion reverses, leading to a lateral movement
to the left. This behavior aligns with the expectation that a negative roll
parameter induces a rotation in the opposite direction, effectively steer-
ing the drone to the left along the Y-axis. The subsequent return to a roll
parameter value of 0 results in the drone retracing its path back to the
starting point. These observed movements underscore the critical role of
the roll parameter in controlling the drone’s lateral displacement, show-
casing its direct influence on the direction and magnitude of the drone’s
movement in the Y-axis coordinate.
(Continued)
Dynamic Modeling of Drone Control 79
Table 5.3 Varying the values of roll parameter of the drone flight control.
Simulation
Values
–0.002
–0.002
0.002
0.002
0
Time
20 s
40 s
60 s
80 s
0s
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
80 Linear and Nonlinear System Modeling
Table 5.3 Varying the values of roll parameter of the drone flight control. (Continued)
Time Values Simulation
100 s 0
(Continued)
Table 5.4 Varying the values of yaw parameter of the drone flight control.
Simulation
Values
–1
0
1
Time
20 s
40 s
60 s
0s
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
82 Linear and Nonlinear System Modeling
Table 5.4 Varying the values of yaw parameter of the drone flight control.
(Continued)
Time Values Simulation
80 s –1
100 s 0
representing the rotation around the vertical axis, plays a crucial role in
determining the drone’s orientation and heading. As observed in the sim-
ulation, when the yaw parameter is set to 0, the drone maintains its initial
heading. When the yaw value transitions from 0 to 1, the drone exhibits
a clockwise twist or rotation to the right. Conversely, when the yaw value
changes from 1 to –1, the drone initiates an anti-clockwise twist or rota-
tion to the left. This dynamic response is indicative of the drone’s ability
to change its orientation based on yaw parameter variations, highlighting
the direct influence of this control parameter on the drone’s rotational
behavior.
Furthermore, the periodic oscillation observed in the yaw values (0, 1,
–1) at regular time intervals (20, 40, 60 s) reflects a deliberate and con-
trolled alteration in the drone’s heading. This pattern repeats, and as the
yaw parameter values oscillate between 1 and –1 at subsequent intervals
(80 s), the drone continues to twist from right to left. The simulation
demonstrates a clear correlation between yaw parameter adjustments and
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Dynamic Modeling of Drone Control 83
Table 5.5 Varying the values of all parameters of the drone flight control.
Time Thrust value Pitch value Roll value Yaw value
0s 0 0 0 0
20 s –1 0.002 0.002 0.5
40 s –0.8 0.002 0.002 –0.5
60 s –1.2 –0.002 –0.002 0.5
80 s –1 –0.002 –0.002 –0.5
100 s 0.2 0 0 0
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
(Continued)
Table 5.6 Varying the values of all parameters of the drone flight control.
84 Linear and Nonlinear System Modeling
Simulation
Time
20 s
40 s
60 s
0s
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Dynamic Modeling of Drone Control 85
Table 5.6 Varying the values of all parameters of the drone flight control.
(Continued)
Time Simulation
80 s
100 s
5.6 Applications
Drones, also known as UAVs, have risen from futuristic fantasy to com-
monplace reality. These versatile machines, equipped with rotors and
remote control, are buzzing into action across diverse fields, transforming
industries and redefining possibilities. Surveillance, environmental moni-
toring and photography and videography industries are some applications
in drone control [3, 4].
Surveillance. Drones have become the eyes in the sky, transforming
surveillance and monitoring. These UAVs equipped with powerful cam-
eras and sensors can gather real-time data from inaccessible or dangerous
areas, providing a comprehensive view that traditional methods struggle
to match. Imagine soaring above a crowd, capturing faces and activities
with unparalleled detail, or swiftly surveying vast landscapes, pinpoint-
ing anomalies hidden from ground-based observers. This is the power of
drone surveillance, offering unparalleled awareness for security, disaster
response, environmental protection, and even wildlife conservation. From
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
86 Linear and Nonlinear System Modeling
5.7 Conclusion
In conclusion, this chapter provided an overview of the design of drone
control and simulation using MATLAB and its Simulink environment. The
chapter provided a comprehensive understanding of the parameters gov-
erning drone flight, enabling proficient navigation of the MATLAB and
Simulink environment. An effective control system capable of precisely
maneuvering the drone was designed and implemented. Through the utili-
zation of Simulink, realistic drone simulations were conducted to evaluate
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Dynamic Modeling of Drone Control 87
References
1. Shelare SD, Aglawe KR, Waghmare SN, Belkhode P (2021b) Advances
in water sample collections with a drone – A review. Mater. Today Proc.
47:4490–4494. https://doi.org/10.1016/j.matpr.2021.05.327
2. Najm AA, Ibraheem IK (2019) Nonlinear PID controller design for a
6-DOF UAV quadrotor system. Eng. Sci. Technol., an International Journal
22(4):1087–1097. https://doi.org/10.1016/j.jestch.2019.02.005
3. MATLAB: https://in.mathworks.com/discovery/drone-simulation.html
4. Goyal, N., Sharma, S., Rana, A. K., & Tripathi, S. L. (Eds.). (2022). Internet of
Things: Robotic and Drone Technology. CRC Press, Boca Raton, UK.
5. Thillaiarasu, N., Tripathi, S. L., & Dhinakaran, V. (Eds.). (2022). Artificial
Intelligence for Internet of Things: Design Principle, Modernization, and
Techniques. CRC Press, Boca Raton, UK.
6. Tripathi, S. L., Ganguli, S., Kumar, A., & Magradze, T. (Eds.). (2022).
Intelligent green technologies for sustainable smart cities. John Wiley & Sons,
USA, doi: 10.1002/9781119816096
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
6
Development of New Bioinspired
Hybrid Algorithms for Parameter
Modeling of Photovoltaic Panels
Souvik Ganguli*, Shilpy Goyal and Parag Nijhawan
Abstract
Three hybrid bioinspired techniques are developed in this chapter. Two popular
nature-dependent methods, viz., the firefly method, abbreviated as firefly algo-
rithms (FA) and the bacterial foraging algorithm (BFA) were integrated utilizing
the benefits of each. Firefly algorithms is employed for exploring the entire search
domain while the three improved BFA schemes proposed by Supriyono have been
utilized to search locally in the global optimization methods. Mainly, the step size
of chemotaxis in bacterial foraging technique is modified with linear, quadratic,
and exponential adaptations to improve upon the basic BFA technique. The bioin-
spired techniques not only outperform their parent methods, but they also surpass
some of the most cutting-edge metaheuristic strategies. A bunch of unconstrained
standard functions are taken up in this work to test the efficiency of the suggested
approaches. The statistical analysis of the optimal values as well their non-para-
metric test also validate this fact. Further, these hybrid algorithms are also suc-
cessfully used for the evaluation of two-diode model parameters using a novel
objective function coined as weighted sum of square error (WSSE). The test results
of the two different photovoltaic systems truly prove that the proposed method is
effective in obtaining the unknown model parameters. The variation of the cost
function with the iteration counts also verify the faster convergence and accuracy
of the results. In future, multi-dimensional diode modeling can be considered to
identify their undetermined parameters with the proposed technique. Instead of
arbitrary weights, optimized weights can be assigned to bring about even better
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(89–108) © 2024 Scrivener Publishing LLC
89
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
90 Linear and Nonlinear System Modeling
results. Parameter modeling of various photovoltaic panels can also benefit from
a multi-objective strategy.
6.1 Introduction
In bioinspired techniques, hybrid algorithms have been developed by inte-
grating two or more algorithms to improve or enhance overall search effi-
ciency. Many researchers have quite often tried to utilize the merits of each
algorithm for the common good [1]. In this chapter, three new bioinspired
hybrid propositions are presented by integrating the firefly algorithm (FA)
[2] with some improved forms of bacterial foraging [3]. The bacterial for-
aging technique, suggested by Passino, is evolved on the basis of the for-
aging phenomenon of E. Coli bacteria that stays within the intestine of the
human body. But, bacterial foraging algorithm (BFA) uses a lot of user-de-
fined parameters, suffers from the drawback of premature convergence,
and its efficacy heavily decreases with dimensions and problem complexity
[4], for which an integrated approach of FA and BFA is conceptualized.
A bacterial foraging system comprises four salient components, namely,
chemotaxis followed by swarming, reproduction, and finally, the elimina-
tion-dispersal. Out of these four key steps of the algorithm, adaptive varia-
tion of the chemotactic step, is widely popular in the literature.
Supriyono [5] proposed one such scheme to bring about the linear,
quadratic, and exponential formulation of the chemotactic step which
improved the standard BFA to a considerable extent. This strategy is used
to yield a modified bacterial foraging algorithm, coined as MBFA in this
chapter. Thus, three global optimization methods are presented combining
FA with MBFA. Firefly algorithm is used for exploration purpose while
three MBFA methods are employed to execute the exploitation feature. The
hybrid algorithms are referred to in this chapter as FABFA-1, FABFA-2,
and FABFA-3, respectively. Both unimodal and multimodal benchmark
functions of appreciable dimensions are considered to test the efficacy of
these new hybrid algorithms. A sufficient number of recently developed
metaheuristic algorithms are set up for a comparative study. Statistical
tests, both parametric and non-parametric, are conducted to validate the
experimental outcomes. The convergence curves also prove the accuracy
of results and the quality of solutions.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 91
q V IRS
I I PV I 01 exp 1
1 KT
q V IRS V IRS
I 02 exp 1 (6.1)
2 KT Rp
where ‘IPV’ represents the photocurrent, while ‘I01’ and ‘I02’ typify the
reverse saturation currents of the respective diodes, ‘Rp’ and ‘RS’ denote,
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
92 Linear and Nonlinear System Modeling
RS
I
IPV I01 I02 RP
V
On solving,
Under the shorted conditions, the voltage ‘V’ will be zero while I = ISC,
denoting the short-circuit current. In that case, Equation 6.1 can be written
as
qI SC RS qI SC RS I SC RS
I SC I PV I 01 exp 1 I 02 exp 1
1KT 2 KT Rp
(6.4)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 93
Thus,
qI SC RS qI SC RS I SC RS
I PV I SC I 01 exp 1 I 02 exp 1
1KT 2 KT Rp
(6.5)
q Vm I m RS
Im I PV I 01 exp 1
1 KT
q Vm I m RS Vm I m RS
I 02 exp 1 (6.6)
2 KT Rp
Similarly, the error function for the short-circuit condition can be for-
mulated with Equation 6.5, as given below.
qI SC RS
errSC I SC I 01 exp 1
1KT
qI SC RS I SC RS
I 02 exp 1 I PV (6.8)
2 KT Rp
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
94 Linear and Nonlinear System Modeling
q Vm I m RS
errm I PV I 01 exp 1
1 KT
q Vm I m RS Vm I m RS
I 02 exp 1 Im (6.9)
2 KT Rp
Thus, the objective function, namely the weighted sum of square error
(WSSE) is framed as
2 2
err w1 errOC w2 errSC w3 errm2 (6.10)
functions. Moreover, these algorithms are further utilized to find out the
different parameters of the two-diode models. Quite a substantial number
of methods are employed to compare with the proposed approaches. The
statistical measures of the results are also reported for the test functions.
In order to verify the significance of the findings, the Wilcoxon’s rank-sum
test [22] is also performed. The plot of the fitness function against the total
number of iterations also proves better convergence speed and accuracy of
the suggested techniques over the compared methods.
proves to be the preferred one based on the results obtained. Even the per-
formance of FABFA-1 and FABFA-2 are better on numerous occasions as
compared to the other techniques compared. To assess the implication of
the outcomes, the Wilcoxon test measuring the rank-sum of the datasets
[22] was performed. The p-values are calculated in this evaluation. The
approximate orders of the results are reported in Table 6.3. The p-values
having the numerical figure greater than 0.05 is considered to be insig-
nificant and hence underlined. Mostly, the p-values are similar on most
occasions, and hence, represented with the order only. Few distinct values
near 0.05 are also indicated in this table.
It is clear from Table 6.3 that almost 95% of the results are significant.
Thus, the results of the proposed bioinspired techniques are validated by
the non-parametric test. To test an engineering application, the proposed
methods are applied for the parameter modeling of the two-diode models
of some commercial solar cells. The datasheets of the solar cells used are
reported in Table 6.4 for the benefit of curious readers.
The various terms used in the Table 6.4 carry their usual significance.
The population size considered for optimizing the PV cell parameters is
50, while the iteration count for this problem is taken up as 1000. Various
methods are employed to find the seven unidentified parameters (IPV,
α1, α 2, RS, Rp, I01, I02) of the two-diode model. The parameter limits are
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 97
(Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
98 Linear and Nonlinear System Modeling
(Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 99
(Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
100 Linear and Nonlinear System Modeling
(Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
102 Linear and Nonlinear System Modeling
Table 6.4 Information sheet for the selected photovoltaic modules [9].
Name of manufacturer Schutten Solar Solarworld
Model number STM6 40-36 Pro. SW255
Cell type Mono-crystalline Poly-crystalline
Vm (V) 16.98 30.90
Im (A) 1.50 8.32
Voc (V) 21.02 38.00
Isc (A) 1.663 8.88
Ns (number of cells) 36 60
T (ºC) 51 25
Table 6.6 Optimal two-diode model parameters using several methods for STM6 40-36 and Pro. SW255.
Model Algorithms IPV α1 α2 RS Rp I01 I02 WSSEmin
STM6 40-36 FABFA-1 1.64 1.07 2.15 0.50 240.72 7.05E–09 4.41E–07 1.23E–05
FABFA-2 1.65 1.08 2.11 0.50 239.87 7.07E–09 4.43E–07 1.21E–05
FABFA-3 1.65 1.09 2.13 0.49 240.34 7.10E–09 4.44E–07 1.15E–05
FA [2] 1.78 1.11 2.78 0.48 235.89 4.29E–07 2.33E–07 0.0056
BFA [3] 1.23 1.56 2.12 0.50 234.18 7.01E–07 3.52E–07 0.0367
GWO [10] 1.68 0.50 1.32 0.50 234.68 0.00000 2.35E–07 0.0010
DA [11] 1.65 1.08 2.14 0.50 241.73 7.12E–09 4.46E–07 0.0001
SCA [12] 1.79 0.50 1.25 0.50 138.37 0.00000 9.02E–07 0.0255
SSA [13] 1.67 4.57 1.26 0.21 159.74 7.29E–07 1.02E–07 0.0003
(Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
104 Linear and Nonlinear System Modeling
Table 6.6 Optimal two-diode model parameters using several methods for STM6 40-36 and Pro. SW255. (Continued)
Model Algorithms IPV α1 α2 RS Rp I01 I02 WSSEmin
Pro. SW255 FABFA-1 8.88 1.50 5.81 0.034 260.81 7.12E–07 8.15E–07 2.57E–05
FABFA-2 8.89 1.49 5.81 0.033 262.67 7.09E–07 8.18E–07 2.05E–05
FABFA-3 8.90 1.51 5.83 0.032 263.23 7.08E–07 8.21E–07 1.38E–05
FA [2] 9.29 0.50 1.12 0.581 178.51 6.87E–07 4.45E–07 0.6428
BFA [3] 9.35 0.52 1.14 0.672 189.21 0.0000 2.39E–07 0.7826
GWO [10] 9.33 0.50 1.09 0.467 218.24 0.00000 1.65E–07 0.5651
DA [11] 8.89 1.51 5.84 0.033 263.26 7.08E–07 8.20E–07 0.0002
SCA [12] 9.70 0.50 1.11 0.500 500.00 0.00000 2.44E–07 1.0399
SSA [13] 9.40 0.50 1.11 0.452 190.02 0.00000 2.13E–07 0.6491
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 105
102
FABFA-3
SSA
GWO
Objective space SCA
Best score obtained so far
DA
100
10–2
10–4
100 200 300 400 500 600 700 800 900 1000
Iteration
Figure 6.2 Convergence plots for DDM (Schutten Solar STM6 40-36).
Objective space
2 FABFA-3
10
SSA
GWO
SCA
DA
Best score obtained so far
100
10–2
10–4
100 200 300 400 500 600 700 800 900 1000
Iteration
below can also be tested in the future to acquire the diode model parame-
ters for photovoltaic systems. Some typical forms of the BFA approach can
be adapted for integration with FA for improved outcomes. A multi-di-
mensional diode model may also be constructed to produce more precise
conclusions. A multi-objective strategy incorporating conflicting aims can
be tried in the near future to find the parameters of the two-diode model.
Since the bioinspired approaches used are stochastic methods, the results
are somewhat subject to the parameter choice and the selection of termi-
nation conditions.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
106 Linear and Nonlinear System Modeling
6.5 Conclusions
Three hybrid bioinspired approaches are discussed in this chapter. The FA
and the BFA are two well-known algorithms with natural-world inspira-
tions that can be combined to maximize their individual strengths. In the
global optimization method, FA is used to probe the entire search space,
while Supriyono’s technique of enhanced BFA schemes is employed for the
local search. Linear, quadratic, and exponential adjustments to the che-
motactic size of the BFA are primarily used to enhance the fundamental
BFA technique. Bioinspired methods not only outperform their ances-
tor methods but also some recently developed metaheuristic techniques.
Using roughly ten widely used test functions, we are able to demonstrate
the efficacy of our method. Moreover, the statistical methods and the
non-parametric test verify the experiment results. In addition, the WSSE
has proven useful with these hybrid algorithms for testing parameters in
the two-diode model. Results from both photovoltaic models indicate that
the proposed method works very well. Convergence plots guarantee accu-
rate and rapid convergence of results. In the future, the proposed method
can be used to calculate the parameters of a more complex diode model.
Assigning optimized weights, as opposed to arbitrary weights, can lead
to even better results. Multi-objective approaches involving two or more
competing goals are another option for determining the parameters of var-
ious photovoltaic panels.
References
1. Hatamlou, A. (2017). A hybrid bio-inspired algorithm and its application.
Appl. Intell., 47(4), 1059-1067.
2. Yang, X. S. (2010). Firefly algorithm, stochastic test functions and design
optimisation. Int. J. Bio-Inspired Comput., 2(2), 78-84.
3. Passino, K. M. (2012). Bacterial foraging optimization. In Innovations and
Developments of Swarm Intelligence Applications (pp. 219-234). IGI Global.
4. Das, S., Biswas, A., Dasgupta, S., & Abraham, A. (2009). Bacterial foraging
optimization algorithm: theoretical foundations, analysis, and applications.
In Foundations of computational intelligence volume 3 (pp. 23-55). Springer,
Berlin, Heidelberg.
5. Supriyono, H. (2012). Novel bacterial foraging optimisation algorithms with
application to modelling and control of flexible manipulator systems (Doctoral
dissertation, University of Sheffield).
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
New Bioinspired Hybrid Algorithms 107
6. Biswas, P. P., Suganthan, P. N., Wu, G., & Amaratunga, G. A. (2019). Parameter
estimation of solar cells using datasheet information with the application of
an adaptive differential evolution algorithm. Renew. Energy, 132, 425-438.
7. Abd Elaziz, M., & Oliva, D. (2018). Parameter estimation of solar cells diode
models by an improved opposition-based whale optimization algorithm.
Energy Convers. Manage., 171, 1843-1859.
8. Jordehi, A. R. (2016). Parameter estimation of solar photovoltaic (PV) cells:
A review. Renew. Sustain. Energy Rev., 61, 354-371.
9. Jadli, U., Thakur, P., & Shukla, R. D. (2017). A new parameter estimation
method of solar photovoltaic. IEEE J. Photovoltaics, 8(1), 239-247.
10. Mirjalili, S., Mirjalili, S. M., & Lewis, A. (2014). Grey wolf optimizer. Adv.
Eng. Software, 69, 46-61.
11. Mirjalili, S. (2016). Dragonfly algorithm: a new meta-heuristic optimization
technique for solving single-objective, discrete, and multi-objective prob-
lems. Neural Comput. Appl., 27(4), 1053-1073.
12. Mirjalili, S. (2016). SCA: a sine cosine algorithm for solving optimization
problems. Knowledge-Based Syst., 96, 120-133.
13. Mirjalili, S., Gandomi, A. H., Mirjalili, S. Z., Saremi, S., Faris, H., & Mirjalili,
S. M. (2017). Salp Swarm Algorithm: A bio-inspired optimizer for engineer-
ing design problems. Adv. Eng. Software, 114, 163-191.
14. Gupta, P., & Sharma, A. K. (2019). Clustering-based Optimized HEED pro-
tocols for WSNs using bacterial foraging optimization and fuzzy logic sys-
tem. Soft Comput., 23(2), 507-526.
15. Chen, H., Zhang, Q., Luo, J., Xu, Y., & Zhang, X. (2020). An enhanced
Bacterial Foraging Optimization and its application for training kernel
extreme learning machine. Appl. Soft Comput., 86, 105884.
16. Selvakumar, B., & Muneeswaran, K. (2019). Firefly algorithm based feature
selection for network intrusion detection. Comput. Secur., 81, 148-155.
17. Jain, L., & Katarya, R. (2019). Discover opinion leader in online social net-
work using firefly algorithm. Expert Syst. Appl., 122, 1-15.
18. Tang, K., Xiao, X., Wu, J., Yang, J., & Luo, L. (2017). An improved multilevel
thresholding approach based modified bacterial foraging optimization. Appl.
Intell., 46(1), 214-226.
19. Chen, Y. P., Li, Y., Wang, G., Zheng, Y. F., Xu, Q., Fan, J. H., & Cui, X. T.
(2017). A novel bacterial foraging optimization algorithm for feature selec-
tion. Expert Syst. Appl., 83, 1-17.
20. Pang, B., Song, Y., Zhang, C., Wang, H., & Yang, R. (2019). Bacterial foraging
optimization based on improved chemotaxis process and novel swarming
strategy. Appl. Intell., 49(4), 1283-1305.
21. Yong, Y., Ma, L., Zhao, J., & Shen, X. (2020, October). A Modified Bacterial
Foraging Optimizer with Adaptive Chemotactic Step in Dynamic Search
Region. In International Conference on Intelligent Computing (pp. 41-52).
Springer, Cham.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
108 Linear and Nonlinear System Modeling
22. Wilcoxon, F., Katti, S. K., & Wilcox, R. A. (1970). Critical values and prob-
ability levels for the Wilcoxon rank sum test and the Wilcoxon signed rank
test. Selected tables in mathematical statistics, 1, 171-259.
23. Jamil, M., & Yang, X. S. (2013). A literature survey of benchmark functions
for global optimization problems. arXiv preprint arXiv:1308.4008.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
7
Power Quality Improvement by
Using PV-Integrated DSTATCOM
Pushpanjali Shadangi1, Sushree Diptimayee Swain1* and Pravat Kumar Ray2
1
Department of Electrical Engineering, O.P. Jindal University, Raigarh (C.G.), India
2
Department of Electrical Engineering, NIT, Rourkela, India
Abstract
This chapter presents a new photovoltaic system (PV)-based compensator for
the elimination of power quality-related issues. In this chapter, PV-based distrib-
uted static compensator (PV-DSTATCOM) with hybrid synchronous reference
frame theory (HSRF) current generation technique is proposed. The motive of
this research is to integrate renewable energy source to the power grid, as the
world is facing huge electric power demand worldwide. Photovoltaic system is
chosen due to its popularity and decaying cost among the other renewable energy
sources. Integration of PV with the conventional power grid system causes huge
power quality issues and these issues can be reduced by using proper compensa-
tion techniques. To meet power quality standards under IEEE, we are using very
effective control strategies based on soft computing techniques. The proposed
PV-DSTATCOM is designed to compensate the harmonics and to maintain the
reactive power transfer in the presence of unbalanced loads. The efficacy of this
system is validated by using MATLAB simulations as well as experimental data.
Keywords: Power quality issues, PV-DSTATCOM, hybrid synchronous
reference frame theory, nonlinear loads, grid stability
7.1 Introduction
As we all know, the use of non-conventional sources of energy is very
much popular and also becoming necessary to meet the required energy
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(109–126) © 2024 Scrivener Publishing LLC
109
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
110 Linear and Nonlinear System Modeling
synchronous reference frame (HSRF) (Swain, Ray, & Mohanty, 2017) control
technique, which is designed to correct supply power factor, eliminate har-
monics, provides load balancing, and to improve the load terminal voltage.
Photovoltaic-based distributed static compensator (PV-DSTATCOM)
(Agarwal, Hussain, & Singh, 2017) (Singh, Kandpal, & Hussain, 2018) pro-
vides shunt compensation (Zhai et al., 2020) (de Almeida Carlos, Jacobina,
& dos Santos, 2018) and injects the current at the point of common cou-
pling (PCC). The PV-DSTATCOM, with some improvements, is able to
provide fast corrections for the dynamically changing loads.
Three algorithms are stated in this chapter, namely, instantaneous reac-
tive power theory (IRPT) (Chen, Li, Zeng, & Wang, 2018) (Wu et al.,
2021), modified instantaneous reactive power theory (M-IRPT) (Patnaik &
Panda, 2013), and HSRF (Swain et al., 2017). Two current control methods
are also considered here with the performances for DSTATCOM (Kumar,
Mishra, & Mekhilef, 2020) systems. These current control methods are
namely, indirect current control (ICC) (Patjoshi & Mahapatra, 2013)
(Meng et al., 2017) (Park, Kwon, & Choi, 2017) and direct current control
(DCC) (Raab, Kramer, & Ackva, 2021) (Mutschler & Marcks, 2002).
This chapter is organized as: Section 7.1 is the introduction. Section 7.2
describes the PV-based DSTATCOM model. Section 7.3 demonstrates the
controller design and different control algorithms along with the ICC and
DCC methods. Section 7.4 discusses the simulation results obtained. Section
7.5 illustrates the implementation of the proposed controller in laborato-
ry-developed prototype, and the last, section 7.6, concludes the chapter.
the controller at MPP. The MPPT is associated there to maintain the pulse
width of converter/inverters. There are various methods to track the MPP
but among all, the perturb and observe technique is more popular.
In this chapter, PV-DSTATCOM is proposed and designed for PQ
improvement. The proposed system is designed by using nonlinear model-
ing. The PV panel model is integrated with DSTATCOM with parallel DC
capacitor along with the PI controller to get the harmonic-free power as
per our requirement. It is known to us that by the use of various inductive,
power electronics-based loads, we cannot completely eliminate the har-
monics from the source supply but we can reduce it up to a certain extent by
using appropriate control algorithm. The PV-DSTATCOM is realized using
three-phase, three-wire, IGBT-based voltage source inverter with parallel
DC capacitor. Control algorithm is mainly designed for supply power factor
improvements, to reduce harmonics, and to maintain voltage stability at the
grid. The PV-DSTATCOM provides parallel compensation and inoculate
the current at common coupling point. The PV-DSTATCOM with some
modifications is able to provide rapid corrections for the rapidly changing
loads such as furnaces and traction. Phase lock loop (PLL) is introduced
to generate firing pulse for insulated gate bipolar transistor (IGBT). The
PV-integrated DSTATCOM is shown in Figure 7.1. The experimental hard-
ware set up of the proposed method is also shown in Figure 7.2.
Rs Ls ISa
Va
Ls Three Phase
Rs ISb
Vb Three Wire
Non-Linear
Rs Ls ISc Load
Vc
Lf Lf Lf
Rf Rf Rf
cd
PV Panel
IS IL
3-Phase RL
Voltage Source
(100 V,50 Hz) XL
If
Lf
Current Sensor
Cdc
Measurement Unit
ADC
DAC PWM
Interfacing Inverter
Board
Oscilloscope
dSPACE
PC
Gate signal
Hysteresis control
PI Limiter Imax
Circuit Reference current generation
Vdc-ref Controller
– Vdc
Vdc PI
Vdc
Vdc* Controller
Ploss
Vsa
Clarke p
Vsb Transformation p ica*
Vsc LPF
a-b-c
Power Reference Reference icb*
ILa Calculation q –q Current
Clarke Current
ILb –1 Calculation Calculation icc*
Transformation
ILc
1 1 v sa
1
v 2 2 2
v sb (7.1)
v 3 3 3
0 v sc
2 2
1 1 iLa
1
i 2 2 2
iLb (7.2)
i 3 3 3
0 iLc
2 2
Equations 7.1 and 7.2 can be used for the calculation of active and reac-
tive powers. The power is shown in Equation 7.3:
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Power Quality Improvement Using PV-DSTATCOM 115
p v v i
(7.3)
q v v i
The real and reactive powers, average, and reactive components are
given by Equations 7.4 and 7.5:
p p p (7.4)
q q q (7.5)
this power is then fed to the low pass filter (LPF) and hence, the power can
be obtained as in Equations 7.6 and 7.7:
p p ploss (7.6)
q 0 (7.7)
i *c 1 v v p ploss
(7.8)
i *c v 2
v 2
v v q
1 0
i *Ca
2 1 3 i *C
i *Cb (7.9)
*
3 2 2 i *C
i Cc
1 3
2 2
Vdc
Vdc
–
PI Controller
Vdc*
PLoss
Vsa Vd Vd
Park's p ica*
Vsb a-b-c to d-q p
Conversion Transformation a-b-c
Vsc p-q LPF
Reference Reference icb*
Vq Vq Powers
ILa Calculation
Current Current
i
ILb Park's -1 Calculation Calculation cc*
Transformation q -q
ILc
1 1 1 v sa
vd 2 cos t sin t 2 2
vsb
vq 3 sin t cos t 0 3 3
2 2 vsc
(7.10)
v cos t sin t vd
v sin t cos t vq (7.11)
load voltage and source currents are detected simultaneously to get the
harmonics. Among various control schemes available, the HSRF method
is superior, it is very efficient and a reliable scheme to avoid the voltage
unbalance, reduce harmonic content, and avoid voltage instability. Here,
three-phase load voltages (VLa, VLb, VLc) and source currents (Isa, Isb,
Isc) are transformed from a-b-c to rotating reference frame axis.
The voltages received can be further separated into two parts, DC and
dynamic components, as given below in Equation 7.12:
Vd Vd Vd
Vd Vq Vq (7.12)
The load voltages are VLa, VLb, VLc, and these are distorted and unbal-
anced types. If d-q-components are sent to the filter (LPF), it is subtracted
from the fundamental signal. In this way, the DC components vanish and
we get the harmonic signals, refer to Equation 7.13:
Vd Vd Vd
Vq Vq Vq (7.13)
Harmonic contained load voltages are given by the below Equation 7.14.
cos sin
VLah
2 2 Vd
VLbh cos( ) sin( ) (7.14)
3 3 Vq
VLch 2 2
cos( ) sin( )
3 3
Id Id I
d
Id
Iq I q (7.15)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
118 Linear and Nonlinear System Modeling
cos sin
I sah
2 2 I d
I sbh cos( ) sin( ) (7.16)
3 3 I q
I sch 2 2
cos( ) sin( )
3 3
isa
K
LPF Vca*
VLa –
isb a-b-c to d-q d-q to a-b-c – Vcb*
K
Transformation Transformation
VLb LPF Vcc*
isc –
K
VLc
Vdc –
PLL Vs
Vdc* Vdc
Figure 7.6 Hybrid control approach-based synchronous reference frame (HSRF) method.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Power Quality Improvement Using PV-DSTATCOM 119
isa
*
usa isa
Vsa
* isb
Unit Vector usb isb Gate Pulses
Vsb Controller
Generator u * isc
sc
Vsc isc
Vdc Imax
–
Vdc
*
Vdc
I * sa I max * usa
I * sb I max * usb
I * sc I max * usc (7.20)
I *ca I * sa I La
I *cb I * sb I Lb
I *cc I * sc I Lc (7.21)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
120 Linear and Nonlinear System Modeling
– i * ica –
usa isa* ca
Vsa
u isb* – i * icb – Gate Pulses
Vsb Unit Vector sb cb
Controller
Generator u
sc isc * – icc* icc –
Vsc
Vdc Imax
Vdc –
Vdc *
Both the methods ICC and DCC have both been analyzed to show good
good total harmonic distortion (THD) spectrum but ICC shows superior-
ity over DCC. The Table 7.1 shows the THD after compensation for both
the control methods ICC and DCC for the range of order of harmonics.
200 10
Vs (V)
IL (V)
0 0
−200 −10
0 0.02 0.04 0.06 0.08 0.1 0 0.02 0.04 0.06 0.08 0.1
Time (s) Time (s)
(a) (b)
50 1
Is (A)
Is (A)
0 0
−50 −1
0 0.02 0.04 0.06 0.08 0.1 0 0.02 0.04 0.06 0.08 0.1
Time (s) Time (s)
(c) (d)
1
Is (A)
−1
0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Time (s)
(e)
Figure 7.9 (a) Power supply voltage during steady-state response, (b) current at load
without filter, (c) source current for resistive load post compensation, (d) source current
for resistive-inductive load post compensation, and (e) source current for resistive-
inductive-EMF load post filtration.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
122 Linear and Nonlinear System Modeling
Mag (% of Fund)
Mag (% of Fund)
20 20
10 10
0 0
0 5 10 15 20 0 5 10 15 20
Harmonic order Harmonic order
(a) (b)
FFT analysis FFT analysis
THD = 34.29% THD = 3.73%
Mag (% of Fund)
Mag (% of Fund)
20 0.8
0.6
10 0.4
0.2
0 0
0 5 10 15 20 0 5 10 15 20
Harmonic order Harmonic order
(c) (d)
Mag (% of Fund)
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 5 10 15 20 0 5 10 15 20
Harmonic order Harmonic order
(e) (f)
Figure 7.10 Source current pre-compensation THD response. (a) For Case 1, (b) Case 2,
(c) Case 3. (d) For Case 1, (e) Case 2, and (f) Case 3 post-compensation THD response.
Case 3: (diode rectifier along with R-L-E load during steady-state condi-
tion): Figure 7.9(e), Figure 7.10(c) and (f): simulation results of current at
source, THD of current at source for pre- and post-compensation during
steady-state condition.
After doing all the analysis, it was concluded that before applying any
compensation in the system, we are getting more THD values for all the
cases, as follows: 24.29%, 24.70%, and 34.29%, respectively, but when we
introduced the compensation, the THD was drastically reduced to 3.73%,
4.23%, and 4.53%, respectively, for all the three cases stated above. Hence,
it is concluded that a suitable algorithm has to be chosen for the operation
of DSTATCOM to get the satisfactory output response.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Power Quality Improvement Using PV-DSTATCOM 123
IL 1A/div.
Is
10ms/div. 10ms/div.
(a) (b)
1A/div.
10ms/div.
(c)
Figure 7.11 Laboratory experimental results for (a) current at load end pre compensation,
current at source end post compensation for Case 1 (b) source current post compensation
for Case 2, and (c) source current post filtration for Case 3.
7.6 Conclusion
Implemented Photo-Voltaic Distribution Static Compensator
(PV-DSTATCOM) with HSRF control algorithm as a reference voltage gen-
eration approach was proposed. The conventional control technique based
PV-DSTATCOM was found to be very simple in realization and also achieve
good performance when considering various load perturbations. But there is
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
124 Linear and Nonlinear System Modeling
References
Agarwal, R. K., Hussain, I., & Singh, B. (2017). Implementation of LLMF control
algorithm for three-phase grid-tied spv-dstatcom system. IEEE Trans. Ind.
Electron., 64 (9), 7414-7424. doi:10.1109/TIE.2016.2630659
Alam, S. J., & Arya, S. R. (2020). Control of UPQC based on steady state linear
Kalman filter for compensation of power quality probelms. Chin. J. Electr.
Eng., 6 (2), 52-65. doi:10.23919/CJEE.2020.000011
Bernstein, A., Wang, C., Dall’Anese, E., Le Boudec, J.-Y., & Zhao, C. (2018). Load
flow in multiphase distribution networks: Existence, uniqueness, non-
singularity and linear models. IEEE Trans. Power Syst., 33 (6), 5832-5843.
doi: 10.1109/TPWRS.2018.2823277
Chen, J., Li, H., Zeng, D., & Wang, G. (2018). High-impedance fault detec-
tion based on single-phase instantaneous reactive power theory. In 2018
IEEE power energy society general meeting (pesgm) (p. 1-5). doi:10.1109/
PESGM.2018.8586198
De Almeida Carlos, G. A., Jacobina, C. B., & dos Santos, E. C. (2018). Alternative
breed of three-phase four-wire shunt compensators based on the cascaded
transformer with single dc link. IEEE Trans. Ind. Appl., 54 (3), 2492-2505.
doi: 10.1109/TIA.2018.2798594
Dongre, A. A., Mishra, Y. P., & Kumar Dubey, A. (2020). Carrier PWM based
capacitor supported dynamic voltage restorer for voltage sag and swell mit-
igation in distribution system. In Michael faraday iet international summit
2020 (mfiis 2020) (Vol. 2020, p. 273-278). doi:10.1049/icp.2021.1066
Jo, S.-H., Son, S., & Park, J.W. (2013). On improving distortion power quality
index in distributed power grids. IEEE Trans. Smart Grid, 4 (1), 586-595.
doi: 10.1109/TSG.2012.2222943
Kalla, U. K., Singh, B., Murthy, S. S., Jain, C., & Kant, K. (2018). Adaptive sliding
mode control of standalone single-phase micro grid using hydro, wind, and
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Power Quality Improvement Using PV-DSTATCOM 125
Swain, S. D., Ray, P. K., & Mohanty, K. B. (2017). Improvement of power quality
using a robust hybrid series active power filter. IEEE Trans. Power Electron.,
32 (5), 3490-3498. doi: 10.1109/TPEL.2016.2586525
Tan, Y., & Wang, Z. (2019). Incorporating unbalanced operation constraints of
three-phase distributed generation. IEEE Trans. Power Syst., 34 (3), 2449-
2452. doi: 10.1109/TPWRS.2019.2895559
Verma, A., & Singh, B. (2021). CAPSA based control for power quality correc-
tion in PV array integrated EVCS operating in standalone and grid con-
nected modes. IEEE Trans. Ind. Appl., 57 (2)1789-1800. doi: 10.1109/
TIA.2020.3041812
Wu, F., Yang, B., Hu, A., Zhang, Y., Ge, W., Ni, L., Zha, Y. (2021).Inertia and
damping analysis of grid-tied photovoltaic power generation system
with dc voltage droop control. IEEE Access, 9, 38411-38418. doi: 10.1109/
ACCESS.2021.3059687
Yang, H., Shen, W., Yu, Q., Liu, J., Jiang, Y., Ackom, E., & Dong, Z. Y. (2020).
Coordinated demand response of rail transit load and energy storage system
considering driving comfort. CSEE J. Power Energy Syst., 6 (4), 749-759. doi:
10.17775/CSEEJPES.2020.02590
Zhai, H., Zhou, F., Zhu, C., Yi, H., Wang, Z., Tao, R., & Wei, T. (2020). An optimal
compensation method of shunt active power filters for system- wide volt-
age quality improvement. IEEE Trans. Ind. Electron., 67 (2), 1270-1281. doi:
10.1109/TIE.2019.2899561
Zhang, J., Liu, Y., Ding, K., Feng, L., Hamelmann, F. U., & Chen, X. (2020). Model
parameter analysis of cracked photovoltaic module under outdoor condi-
tions. In 2020 47th IEEE photovoltaic specialists conference (PVSC) (p. 2).
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
8
Modeling and Simulation of Current
Transformer to Study Its Behaviors
in Different Conditions
Aveek Chattopadhyaya1* and Surajit Chattopadhyay2
2
Department of Electrical Engineering, GaniKhan C. Institute of Engineering and
Technology, Malda, W.B., India
Abstract
Current transformer (CT) is a pertinent device in electrical systems. This chap-
ter aims to discuss modeling and simulation of current transformer to study its
behaviors. For this purpose, CT secondary signal is considered. Applying wave-
let transform (WT) and fast Fourier transform (FFT), different behaviors of CT
are analyzed. Wavelet-transform-based multi-resolution analysis (MRA) and root
mean square (RMS) are also used to ascertain the saturation of CT and to calculate
the exact primary current during saturation.
8.1 Introduction
For protection and measuring purposes, Current transformer (CT) is used.
Current transformer is widely used in electrical systems for smooth and
reliable operation of the electrical system. For this purpose, proper mod-
eling, simulation, and behaviors of CT in different conditions are utmost
required things. Otherwise, maloperation of CT can affect the performance
of the system. In Reference [1], computer simulation is discussed to show
the saturation and as a repercussion the deformation of the CT secondary
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(127–142) © 2024 Scrivener Publishing LLC
127
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
128 Linear and Nonlinear System Modeling
due to direct current (DC) component in the CT primary side. In this con-
dition, core hysteresis and the effect of residual flux on the CT performance
were also investigated. In Reference [2], saturation of CT was detected in
transient conditions using sampled value of CT secondary current [2]. Y.
C. Kang et al. proposed first, second, and third difference-based functions
of CT secondary waveform to detect CT saturation [3]. In this analysis, the
first difference of the CT secondary signal was used to identify the initia-
tion of saturation. The above discussed method was tested for 100:5 CT.
Negative and positive sequence impedance-based protection scheme is
developed throughout saturation and ratio discrepancy [4]. The developed
algorithm was verified in digital platform. N. Villamagna et al. proposed
component (symmetrical)-based saturation detection technique in differ-
ent fault conditions, where faults were created in EMTP/ATP simulation
platform [5]. A digital signal processor (TMS3206701) was used to test the
CT saturation detection algorithm, where the difference (third) of second-
ary waveform was considered for this purpose [6]. A. Hooshyar et al. pro-
posed a method to calculate primary component of fault signal in existence
of DC component and saturation [7]. The method is based on the look up
table of fault currents and least square errors. Waveform processing and
artificial intelligence (AI)-based method proposed by researchers to iden-
tify and rectify saturation for smooth operation about electrical systems
[8–9]. In this analysis, signal processing such as wave transform (WT) was
used to identify saturation and AI-like artificial neural network was also
used to correct saturation. A. Chattopadhyaya et al. proposed different sig-
nal processing and signal processing-based statistical parameter analysis
to assess the distorted signal in power systems [10–17]. In Reference [18],
authors proposed a statistical criterion-based technique to detect CT satu-
ration and reconstruction of CT distorted signal during saturation. Here,
authors also discussed the statistical methods to discriminate fault current
from inrush current. Monitoring the flux density, a method was developed
to identify saturation [19].
The chapter aims to present the behaviors of CT in non-identical envi-
ronments and its signal estimation on those conditions by modeling and
simulation of CT. Knowledge of behaviors of CT is very essential for appli-
cation and design of a CT. In this chapter, CT behaviors will also be dis-
cussed considering CT secondary signal and CT core flux. This assessment
is achieved by FFT and statistical parameters based on discrete wavelet
transform (DWT). This chapter will also discuss to calculate CT primary
current under CT saturated conditions. Current transformer (CT) satura-
tion for different reasons is also considered here. Here, modeling and cir-
cuit simulation of CT is discussed, which is followed by behaviors analysis,
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Modeling and Simulation of CT 129
CT Primary Current
Breaker
RL Load
AC Source Secondary
Current
SCOPE
Primary Current Secondary Voltage
Source Measuring Unit
CT Measuring Unit
Multimeter
Case Study I:
Time constant (L/R) of the circuit is 3
In Case study I, the CT performance was observed when time constant
(L/R) of the circuit is 3. Fast Fourier transform-based spectrum analysis
was done on the primary and secondary signal of CT and CT core flux
under this condition (time constant is 3). In spectrum analysis, each fre-
quency component consists of first positive peak (FPP), second positive
peak (SPP), and negative peak (NP). Figure 8.3(a) and 8.3(b) depicts the
CT primary signal and secondary signal, and core flux, respectively. In
Figure 8.3(a), the primary and secondary signal of CT is overlapping each
other, which indicates normal condition. Spectrum analysis of CT primary
signal, secondary signal, and core flux are shown in Figure 8.3(c), 8.3(d),
and 8.3(e), respectively where no harmonics are observed in this condi-
tion. Results of the spectrum analysis dictate that for time constant 3, the
CT operates in linear zone.
1
I1 & V2
–1
–2
–3
–4
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
time in second
(a)
flux Vs time
1
0.8
0.6
0.4
flux
0.2
–0.2
–0.4
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
time in second
(b)
Peak Magnitude Spectrum called by Simulink
4000
Magnitude based on “Base Peak” - Parameter
3000
2000
1000
–1000
–2000
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(c)
Figure 8.3 (a) CT primary signal and secondary signal. (b) Core flux of CT. (Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
132 Linear and Nonlinear System Modeling
10
–2
–4
–6
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(d)
Peak Magnitude Spectrum called by Simulink
1.5
Magnitude based on “Base Peak” - Parameter
0.5
–0.5
–1
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(e)
Figure 8.3 (Continued) (c) Spectrum analysis of CT primary current. (d) Spectrum
analysis of CT secondary voltage. (e) Spectrum analysis of CT core flux.
100
50
I1 & V2---->
–50
–100
–150
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time in second---->
(a)
flux Vs time
20
15
10
5
flux---->
–5
–10
–15
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time in second---->
(b)
× 104 Peak Magnitude Spectrum called by Simulink
12
Magnitude based on “Base Peak” - Parameter
10
8
6
4
2
0
–2
–4
–6
–8
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(c)
Figure 8.4 (a) CT primary signal and secondary signal. (b) Core flux of CT.
(c) Spectrum analysis of CT primary current. (Continued)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
134 Linear and Nonlinear System Modeling
250
150
100
50
–50
–100
–150
–200
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(d)
30
20
10
–10
–20
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
(e)
Figure 8.4 (Continued) (d) Spectrum analysis of CT secondary voltage. (e) Spectrum
analysis of CT core flux.
–1
–2
–3
–4
0 2 4 6 8 10 12 14 16 18
Order of Harmonic
Figure 8.5 Spectrum anqalysis of CT secondary voltage (at breaker closing time is
0.02 second).
higher frequency components and output of low pass filter contains lower
order frequencies. Detail coefficients are the output of high pass filter and
approximation or approximate coefficients are the output of low pass filter.
In this chapter, at different switching times, the CT secondary signal
was decomposed up to DWT level nine and at each and every level, skew-
ness values were calculated to detect CT saturation. Figure 8.6 to Figure
8.9 delineate the detail and approximate or approximation coefficients
skewness values (up to DWT decomposition level nine). In Figure 8.9,
–1
Skewness
Approximate
–2
–3
Detail
–4
–5
–6
1 2 3 4 5 6 7 8 9
DWT decomposition level
–1
Skewness
–2
Detail
–3
–4
–5
–6
1 2 3 4 5 6 7 8 9
DWT decomposition level
–1
Skewness
–2
Detail
–3
–4
–5
–6
1 2 3 4 5 6 7 8 9
DWT decomposition level
–1
Approximate
–2
Skewness
–3 Detail
–4
–5
–6
–7
–8
1 2 3 4 5 6 7 8 9
DWT decomposition level
the skewness values of detail coefficients are totally different from the above
limned figures. Another important observation is that in normal condition
(switching time 1.25/60), the intersecting point between approximation
and detail coefficients skewness values is between level four (4) and five
(5) (in DWT level), but for otherdelineated figures, this point is in between
DWT decomposition level two (2) and three (3). From these skewness val-
ues, the CT saturation or distortion of CT secondary signal can be assessed
easily.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
138 Linear and Nonlinear System Modeling
2
1
0
–1 First
–4 notch
–3
–4
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (sec)--->
Figure 8.10 Primary and secondary signal, DWT (detail level1) coefficient.
0
–20
–40 First and six
–60 notch
–80
–100
–4
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Time (sec)--->
For 1000-Amp CT primary current, at DWT detail level one (1) coef-
ficients, no notches were observed but the maximum value of 0.00002351
(which is almost zero) is considered for analysis. Due to flow of high pri-
mary current (25 kAmp), the starting notches positive peaks are inflated
with respect to normal condition. For analysis of saturation, these values
were used. From Table 8.1 results, it is clear that CT saturation in each
cycle can be detected by measuring the peak values of starting and ending
notches. Table 8.1 is used to delineate the magnitude of starting and ending
notches based on which CT saturation was detected.
Table 8.1 Starting and ending notches for different primary current.
Positive
Known Starting peak value No. of Positive peak
primary notch (starting ending value ( end
current number notch) notch notch)
1000 Amp 1 0.00002351 1 0.00002351
2 0.00002351 2 0.00002351
3 0.00002351 3 0.00002351
4 0.00002351 4 0.00002351
5 0.00002351 5 0.00002351
6 0.00002351 6 0.00002351
7 0.00002351 7 0.00002351
8 0.00002351 8 0.00002351
25,000 Amp 1 9.505 1 4.281
2 11.74 2 1.95
3 12.56 3 3.723
4 12.61 4 3.335
5 8.281 5 3.595
6 12.95 6 3.786
7 13.33 7 4.518
8 6.307 8 4.783
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
140 Linear and Nonlinear System Modeling
20
15
10
5
0
0 2000 4000 6000 8000 10000 12000 14000 16000
Primary Current in amps
Figure 8.12 Different primary currents versus RMS value of approximate coefficients.
8.8 Conclusion
In this chapter, circuit simulation and behaviors of CT were discussed by
signal processing-based techniques. Saturation of CT was discussed in
different conditions by signal assessment of CT secondary side. The CT
saturation due to switching transient and high primary current was also
discussed. The presence of different harmonics in CT saturated condi-
tions was analyzed by FFT and DWT-based statistical parameter analysis.
The MRA-based RMS technique was also discussed to measure the actual
primary current (RMS) during CT saturation. All these things are very
important for proper protection of electrical systems.
References
1. O’Kelly, D. (1992). Calculation of the transient performance of protective
current transformers including core hysteresis. IEE Proceedings-c, 139(5),
455-460.https://10.1049/ip-c.1992.0063.
2. Evans, F. J., Wells, G. (1970). Use of a sampling scheme to detect transient
saturation in protective current transformers. IEEE Trans. Instrum. Meas.,
19(3), 144-147. https://10.1109/TIM.1970.4313884.
3. Kang, Y. C., Ok, H., Kang, S. H., Crossley, P. A. (2004). Design and evalua-
tion of an algorithm for detecting current transformersaturation. IEE Proc.-
Gener. Transm. Distrib., 151(1), 27-35. https://10.1049/ip-gtd:20030982.
4. Sachdev,M. S., Sidhu,T. S., Gill, H. S. (2000). A busbar protection technique
and its performance during CT saturation and CT ratio-mismatch. IEEE
Trans. Power Delivery, 15(3), 895-901. https:// 10.1109/61.871349.
5. Villamagna, N., Crossley, P. A. (2006). A CT saturation detection algorithm
using symmetrical components for current differential protection. IEEE
Trans. Power Deliv., 21(1), 38-45. https://10.1109/TPWRD.2005.848654.
6. Kang, Y. C., Ok, S. H., Kang, S. H. (2004). A CT saturation detection algorithm.
IEEE Trans. Power Deliv., 19(1), 78-85. https://10.1109/TPWRD.2003.820200.
7. Hooshyar, A., Pasand, M. S. (2012). Accurate measurement of fault currents
contaminated with decaying DC offset and CT saturation. IEEE Trans. Power
Deliv., 27(2), 773-783. https://10.1109/TPWRD.2011.2176965.
8. Hong, Y. Y., Chang-Chian, P. C. (2008). Detection and correction of distorted
current transformer current using wavelet transform and artificial intelligence.
IET Gener. Transm. Distrib., 2(4), 566–575. https://10.1049/iet-gtd:20070383.
9. Yu, D. C., Cummins, J. C., Wang, Z., Yoon, H. Z., Kojovic, L. A. (2001). Correction
of Current Transformer sistorted Secondary currents due to saturation using
Artificial Neural Networks. IEEE Trans. Power Deliv., 16(2), 189-194. https://
10.1109/61.915481.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
142 Linear and Nonlinear System Modeling
Abstract
Multilevel inverters have discovered more extensive extent of utilization in mod-
erate and also in high-power adjustable-speed drives. This chapter introduces a
vector control scheme of induction motor drive which includes pulse width mod-
ulations for reducing harmonics and total harmonic distortion (THD). For better
control of induction motor, indirect vector control has been applied which offers
advantages such as removal of flux sensor, more dynamic responses, and mini-
mum torque pulses is applied. The inverter named neutral point clamped inverter
is applied for observing dynamic control of the motor drive along with efficiency.
The main attention of this chapter is to improve the performance of indirect vector
controller. The THD analysis proves the better operation of induction motor as
compared to conventional voltage source inverter-fed induction motor drive. By
the help of MATLAB simulation, the dynamic performance as well as steady-state
of multilevel inverter-based drive are analyzed.
Keywords: Vector control, neutral point clamped inverter (NPCI),
synchronously rotating reference frame (SRRF)
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(143–158) © 2024 Scrivener Publishing LLC
143
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
144 Linear and Nonlinear System Modeling
9.1 Introduction
The acknowledgement engine has just been freed from fixed speed machine
to variable speed with the assistance of various control techniques. If a suit-
able control technique is not implemented properly, the benefits of variable
frequency operation cannot be totally recognized. The assumption of method
assumes a significant job in deciding the general feasibility, its attributes,
and execution of the engine drive. During ordinary activity of the induction
drive, the chosen control procedure is needed to guarantee that the activity
is restricted to the zone of high working torque. It must be considered that
it synchronizes the ratings of inverter which helps to reduce the switching
losses as well as the system loss. By the help of control strategies rather than
design of the system, we are able to manage the overload conditions or fault
conditions. One of the most common motor is induction motor, which is
used in most industries, signal control systems, and in the main AC power-
driven applications. And the abovementioned are significant benefits of
induction motor (IM) drives. It is simple and possesses rugged construction,
lower cost, and less maintenance. There are varieties of induction motors
available for study. Efficient control of three-phase induction motor requires
proper concept of properties of motors. In this work of closed loop induc-
tion motor drive, a suitable converter named tandem converter was applied,
which proved to minimize the THD of the whole system and display better
performance. The induction motor performance encourages the implemen-
tation in AC motors. For the better control of induction motor drive several
typologies have been developed in the last many years. Decoupling mecha-
nism that works fine in case of separately excited DC motor was applied to
AC asynchronous machine. Here, the field oriented control was identified
to be suitable for obtaining the outputs. The outputs being the mechanical
and electrical as rotor current and speed of the machine. Along with that,
the chapter includes the THD performance of the machine using a suitable
inverter topology. The approach is to synthesize the machine currents in
order to generate better decoupling mechanism. As field oriented control is
one of the most suitable approach, the same is used for closed loop control.
The neutral point clamped multilevel inverters (NPCMLI) are also known
as diode clamped inverter technically. Here in the voltage level, the middle
point is referred to as the inverter’s neutral point. The increased inverter’s
potential has doubled the level of voltage of the device without precise volt-
age regulation. In various inverter topologies, voltage stress is a central point
influencing the operation of multilevel inverters (MLI). In diode clamped
MLI, diodes are being implemented in the circuit to restrict the stress on
power devices. If N = number of levels, then the voltage sources are counted
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
MLI-FED Closed Loop Control and Analysis of Induction Motor 145
as [N-1], switching conductor devices are assumed as 2[N-1], and the num-
ber of diodes required are [N-1][N-2]. It provides low THD at the output, it
operated at both high-level and low-level frequencies, and low rating devices
can be used; these are some of the advantages of NPCMLIs [1].
The objective of the chapter is to demonstrate the study and research
outcome of NPCMLI of five-level with loads operating at different or wide
modulation index. Here, the switches are selected in such a way that it
minimizes the common mode voltage (CMV) of the inverter and the influ-
ence of the harmonics on the current at the load side is well taken care.
After proper analysis of the model and with the result obtained, five-phase,
three-level neutral point clamped (NPC) inverter proved dynamic for
small changes in load parameters. In this chapter, the study of the model
along with the outputs is done with MATLAB simulation [2–4].
Various other researchers worked in quite similar areas which may also be
corelated for the future scope of the study. Here the improvement in speed
control technique of induction motor drive has offered a wider scope of appli-
cations. This study analyzed the reducing harmonic distortions and CMV by
using appropriate switching techniques in case of five-level NPC inverter feed-
ing three-phase asynchronous motor drive [5–7]. Mostly, NPC is applied for
aircrafts, motor drives, and also in HVDC and AC transmission lines. In this
method, transformation of three phases into two phase speed and time depen-
dent system abc to dq is considered. One of the components defines the mag-
netic flux and other corelates the torque of the motor. Vector control topology is
one of the unique concepts which involves decoupling of flux and torque com-
ponents in order to obtain dynamic responses. Proportional-integral (PI) con-
trollers are typically used in vector control to maintain the measured currents
at their reference levels. It improves torque response, dynamic speed accuracy,
and low power consumption; these are some advantages of the vector control
method. Here, the operation of indirect vector control is discussed in detail.
The fundamental idea behind vector control is to decouple the machine flux
and torque which lets us operate on any of the parameters independently with-
out hampering the performance of the other. The instantaneous component
of the stator currents here get transformed into a frame rotating at synchro-
nous speed that may be aligned with the stator flux, rotor flux , and airgap
flux vectors. In the process, two components, i.e., direct (d)-axis and quadra-
ture(q)-axis components are produced. Both the d and q axes components con-
trol the torque as well as machine flux, respectively. The control becomes easier
and simpler and the performance of the asynchronous machines becomes
dynamic in nature with faster change in load. In our case, the machine used is
a three-phase induction motor. The machine drive is considered to be in vector
control only if the stator current vector component and the decoupled com-
ponents in synchronously rotating reference frames (SRRF) match each other.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
146 Linear and Nonlinear System Modeling
+ S1 + S1 + S1
VC1 C1 VC1 C1 VC1 C1
S2 S2 S2
S3 S3 S3
VC2 C2 VC2 C2 VC2 C2
S4 S4 S4
Vdc
Input S1’ S1’ S1’
VC3 C3 S2’ VC3 C3 S2’ VC3 C3 S2’
S3’ S3’ S3’
VC4 C4 VC4 C4 VC4 C4
S4’ S4’ S4’
– – –
Vdc INVERTER IM
Ibes*
Iaes*
Ices*
Current Error Block
Ibs*
Ias*
Ics*
Field Orientation and Reference
Current Generation
Field weakening
Estimator Block for ids* iqs* ω2 *
Block
Ref Torque
T*
Limiter
Torque
T
Speed Controller
ωre
ωr* ωr
acurrent
Discrete
Ts = 5e-05 s.
DC +
A
Tm is_abc
DC Voltage Source
A
m Wm
m
B B
C Te
DC -
Demax
C
Reference current generation
torque
field orientation and
is_abc*
pulses
Te*
speeda
is_abc
W+
1
–
2
Clock
time
Figure 9.3 Simulation of vector control of IM drive using neutral point diode clamped
inverter (NPC).
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
148 Linear and Nonlinear System Modeling
dimr *
ids (n) imr * (n) r (9.1)
dt
T * (n)
iqs (n) (9.2)
Kimr * (n)
* iqs * (n)
2 (n) (9.3)
r imr (n)
Lr
r
Rr
3 p M.
K
2 2 1 r
Lr (1 r )M . (9.5)
Lr 3
And r 1 where M Lm , Lm is the inductance of magnetism.
M 2
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
MLI-FED Closed Loop Control and Analysis of Induction Motor 149
1 +
w*
+
2 – +
Kp
w
1
– Gain 1
z
+
Unit Delay
1
Ki + Te*
+
Gain
Unit Delay 1
1
z
2 + 1
imr + ids
-K- du/dt
Tr Derivative
1
T
0.1209 u[1]/(u[2]*u[3]) 2
iqs
k Fcn
0.17296 u[1]/(u[2]*u[3]) 3
w2
Tr 1 Fcn 1
regarding the motor’s feedback speed. The system’s output is the excitation
current. The magnetization current remains the same or constant when it is
less than the minimum speed. For speed higher the minimum rated speed,
the magnetization current changes inversely proportional to the speed [14–
19]. The mathematical relation of the mentioned detail is demonstrated,
flux.The current iqs is the q-axis component of the stator and is actually
controlled for smooth operation. The decoupling or the independence in
the operation between these two components is possible only if the iden-
tification of the rotor flux is done correctly or measured accurately. Direct
flux sensors or flux estimators work fine for the abovementioned [20, 21].
The speed controller produces an output torque and using limiter, (T*) is
considered to be the reference torque. The field weakening controller leads
to an output as magnetizing current, which in turn is proportional to flux.
These input signals are fed to the control system or field-oriented control-
ler for estimating the components, torque (iqs*), and the flux (ids*).
Here, the slip frequency is also estimated along with the other two. The
basic objective remains to operate the AC asynchronous motor similarly
like separately excited DC motor. This phenomenon is achieved as the
torque and flux components remain decoupled in the SRRF [22–25].
9.3 Results
Here, the simulation results of the neutral point diode clamped MLIs
with resistive-inductive load are presented. The Figures 9.5(a) and 9.5(b)
includes the current waveform and harmonic study of five-level MLI
inverter simulated in the MATLAB software. The results helped in imple-
menting the inverter feeding the induction motor with specific parameters.
The inverter was used for performing closed loop simulation and analyz-
ing the output characteristics. The outputs of the induction motor were
also demonstrated in the figures. The outputs can be found as dynamics
of speed which are classified as starting dynamics, reversal dynamics, and
load perturbations.
Figure 9.5 load current waveform of five-level MLI connected to RL
load. The THD analysis is also done on the five-level inverter, as shown in
the figure below.
Simulation was performed on a 1-HP induction motor for the control
logic and the results are explained for starting and reversal below. The
specifications of squirrel cage induction motor are mentioned in the Tables
9.1 and 9.2.
2.5
1.5
1
current(amp)
0.5
–0.5
–1
–1.5
–2
–2.5
2.5 2.55 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 3
time(s)
(a)
FFT analysis
Fundamental (50Hz), THD = 4.60%
2.5
Mag (% of Fundamental)
1.5
0.5
0
0 2 4 6 8 10 12 14 16 18 20
Harmonic order
(b)
Figure 9.5 (a) Lload current waveform of 5 level MLI connected to RL load. (b) THD
analysis of MLI with load current connected to RL load.
torque is limited with a value twice the rated value of motor and the result
of which the initial current was limited proportional to the value of torque
limit while the motor gains the speed at required torque in the beginning.
At the speed error reaching nearly equal to zero value in rad/s, the current
in the winding decreases to a value equal to no load and the torque devel-
oped raises to a value balancing the load torque as found in the responses
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
MLI-FED Closed Loop Control and Analysis of Induction Motor 153
in the figure shown below. The speed reaches the desired value at steady
state condition and developed torque becomes minimum, as shown in
Figure 9.6.
200
Ref Iabcs
0
–200
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
200
Act Iabcs
0
–200
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
300
200
Wact
100
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
400
Torque
200
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
time(sec)
200
Ref Iabcs
–200
0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4
200
Act Iabcs
0
–200
0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4
200
Wact
0
–200
0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4
0
Torque
–200
–400
0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4
time(sec)
the currents in the stator showing the responses of the braking is followed
by the reverse of phase current while rotating the rotor in the opposite
motion. As the asynchronous motor runs through the dynamic situation at
no load before and after the reversal of the rotor rotation, the current of the
converter at steady state condition are very much identified to be similar
in both directions of the rotor rotation, i.e., starting and reversal of motor
rotations. In phase sequences, the current magnitudes are observed to be
different in different directions.
controller by keeping the speed fall or rise in the safe limits during constant
load variations. High overshoot in speed may lead to instability in main-
taining safety in motor application. The motor, as a result, settles at the
value of the reference speed. The load perturbations are demonstrated in
Figure 9.8. In the figure, it can be observed that the load has been applied
and removed at instants and the electrical and mechanical changes in the
output characteristics are responding immediately.
The load perturbations are demonstrated in Figure 9.9. The dramatic
nature of the load current and speed can be seen in the figure.
5
Ref Iabcs
0
–5
2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3
5
Act Iabcs
0
–5
2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3
200
Wact
0
–200
2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3
10
Torque
0
–10
2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3
time(sec)
80
Magnitude
60
40
20
0
0 2 4 6 8 10 12 14 16 18 20
Harmonic order
9.4 Conclusion
The NPC topology was adopted for higher level of power applications as it
can achieve better harmonic reduction and the prior existing control strate-
gies help to reduce switching losses. Using this configuration, the stress level
of voltage on high frequency power electronics devices is reduced dynami-
cally compared to conventional low level inverters. Initial uses involved the
railroad traction areas and steel industries in Japan and Europe. There are
several advantages of this structure, which are, major reduction of switch-
ing losses in particular and provides low THD at the output. Here, the study
was applied to induction motor for the realization of the potential of the
methodology. The concept is applicable for smaller and medium applica-
tions. But its advantages can be tested on high rated applications. Analysis
was done using MATLAB simulation tool and the signals can be generated
to realize the same on real time test bed. The converter topology and the
loop controller are well synchronized for obtaining results and the potential
can still be extended during further studies in future.
This scheme can be used in DFIG wind turbine to develop and regulate
the speed of doubly fed induction generator (DFIG) for optimized power
generation especially when the wind energy conversion system is grid tied.
The above results show that the system is dynamic and is reliable. More
reliable algorithm can be implemented for enhancing the results and can
be implemented when hardware for the same is considered. Field program-
mable gate arrays (FPGA) can be one solution for furthering the hardware
and find effective results for practical applications. The model reference
adaptive systems can be one of several other approaches here. The above
results and interface may also have interdisciplinary extensions as the soft
computing techniques are extensively used for smoothening or enhancing
the applicability. Machine learning algorithms also have a large presence in
the recent scenario of technology development.
References
1. Anwer Hossain, Md., Tariquzzaman, Md., Kumar Podder, Amit., Shahadath
Hossain Sabbir, Md., Habibullah, Md. (2019). Predictive Current Control of
a Simplified Three-Level Neutral-Point Clamped Inverter.3rd International
Conference on Electrical, Computer & Telecommunication Engineering (ICECTE).
2. Das Adhikary, Rabindranath., Halder, Tapas., Balow, Writwik. (2019). Diode
& Neutral Point Clamped Five-Level Inverter For the Power Quality Issues.
978-1-5386-6722-4/19/$31.00 © IEEE.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
MLI-FED Closed Loop Control and Analysis of Induction Motor 157
3. Pillay, T.L., Saha, A.K. (2018). A comparative analysis of three phase npc
multi-level inverter IEEEPES/IAS.
4. Bheemaiah, Ch., RanjanMuduli, Utkal., Kumar Behera, Ranjan. (2018).
Performance Comparison of Five-Phase Three-Level NPC to Five-Phase
Two-Level Voltage Source Inverter. 978-1-5386-9316-2/18/$31.00 IEEE.
5. Novotny, D.W. and Lipo, T.A. (1996).Vector Control and Dynamics of AC
Drives. Oxford University Press.
6. N. Mistry, Jigar., Vakharia, Krishna., Tailor, Dhaval., Patel, Brijesh. (2016).
Software and Hardware Implementation of Five level Neutral Point Clamped
Inverter 978-1-4673-9939-5/16/$31.00 © IEEE.
7. Rodriguez, Jose, Bernet, Steffen, K. Steimer, Peter. (2010) A Survey on
Neutral-Point-Clamped Inverters. IEEE Trans. Ind. Electron., VOL. 57, NO.
7, JULY.
8. Rodriguez, J., Bernet, S., Wu, B., Pontt, J.O., and Kouro, S. (2007). Multilevel
voltage-source-converter topologies for industrial medium-voltage drives.
IEEE Trans. Ind. Electron. vol. 54, no. 6, pp. 2930–2945, Dec.
9. Tekwani, P. N., Kanchan, R. S., Gopakumar, K., and Vezzini, A. Fivelevel
inverter topology for induction motor drives with common-mode voltage
elimination in complete modulation range. (2007). EPE J, vol. 17, no. 1,
pp. 11-23, January March.
10. IMECS, Maria., Trzynadlowski, Andrzej M. (2005). Vector Control Schemes
for Tandem-Converter Fed Induction Motor Drives. IEEE Trans. Ind.
Electron., VOL. 20, NO. 2, MARCH.
11. Shiri, A., Vahedi, A., Shoulaie, A. (2006).The effect of parameter variations
on the performance of indirect vector controlled induction motor drive.Power
Electronicsand Motion Control 5thInternationalConference Vol- 2, Page(s):1–5.
14-16 Aug.
12. Texas instruments. (2003). Ac induction motor control (Three phase), sys-
tem document Digital control system group, 30 April.
13. Leite, V., Araujo, R., Freitas, D. (2003) A new online identification meth-
odology for flux and parameters estimation of vector controlled induction
motors. Electric Machines and Drives Conference, Page(s):449-455 vol. 1. 1-4
June.
14. Barrero, F., Perez, J., Millan, R., Franquelo, L.G. (1999) Self-commissioning
for voltage-referenced voltage-fed vector controlled induction motor drives.
Industrial Electronics Society, The 25th Annual Conference of the IEEE,
Page(s):1033-1038 vol.3. 29 Nov. - 3 Dec.
15. Blasco-Gimenez, R., Asher, G.M., Sumner, M., Cilia, J., Bradley, K.J. (1996)
Field weakening at high and low speed for sensor less vector controlled
induction motor drives.Power Electronics and Variable Speed Drives, Sixth
International Conference, Page(s):258–261, 23-25 Sept.
16. Krishnan, R., Bharadwaj, A.S. (1991) A review of parameter sensitivity and
adaptation in indirect vector controlled induction motor drive systems.
Power Electronics, IEEE Transactions on Vol-6, Page(s):695–703. Oct.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
158 Linear and Nonlinear System Modeling
17. Dr. Tez, E.S. (1995) A Simple Understanding of Field-orientation for Ac Motor
Control, Loughborough University of Technology, pp: (3/1-3/4),
18. Prof. Zinger, Don S. (1990). Field Oriented Control of Induction Motors,
University of Akron, Vol-1, pp: 391-396.
19. Krishnan R. (1996). Review of flux-weakening in high performance vector
controlled induction motor drives.Industrial Electronics, Proceedings of the
IEEE International Symposium Page(s):917-922 vol-2, 17-20 June.
20. Ohdachi, Y., Kawase, Y., Hirako, M. (1995) Dynamic analysis of vector con-
trolled induction motor using finite element method. IEEE Trans., Vol-31,
Page(s) 1904–1907. May.
21. Singh, Bhim., Singh, B.N., Singh, B.P., A. and Chandra, K., Haddad, Al,
(1995) Unity Power Factor Converter-Inverter Fed Vector Controlled Cage
Motor Drive Without Mechanical Speed Sensor. IEEE Trans.
22. Singh, Bhim., Ghatak, Choudhuri, Sumit. (2005) Simplified Model for
Simulation of Vector Controlled Induction Motor Drive (VCIMD) Using
Matlab. International Conference On Computer Applications in Electrical
Engineering Recent Advances, Sep. 29-Oct. 1.
23. Tahami, F., Shojaei. (2006) A Novel Fault Tolerant Reconfigurable Concept
for Vector Control of Induction Motors. Power Electronics and Motion
Control, 12th International Conference. Page(s):1199-1204 Aug.
24. Shiri, A., Vahedi, A., Shoulaie, A. (2006) The effect of parameter variations
on the performance of indirect vector controlled induction motor drive.
Power Electronics and Motion Control 5thInternational Conference. Vol-2,
Page(s):1–5. 14-16 Aug.
25. Liu, Xu., Ruan, Yi., Zhang, Chaoyi., Sheng, Huanqing., Yang, Yong. (2008).
On Speed Sensorless Vector Control System for Induction Motor Based on
Estimating Speed by Torque Current Differential. Proceedings of the 27th
Chinese Control Conference. July 16-18.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
10
Hybrid Grey Wolf Optimizer for
Modeling and Control of Electric Drives
Souvik Ganguli1* and Prasanta Sarkar2
1
Department of Electrical and Instrumentation Engineering, Thapar Institute of
Engineering and Technology, Patiala, Punjab, India
2
Department of Electrical Engineering, National Institute of Technical Teachers’
Training & Research, Kolkata, West Bengal, India
Abstract
Permanent magnet synchronous motor drives, or PMSM drives for short, use
integrated nature-dependent metaheuristic algorithms in this chapter to perform
reduced order modeling and controller design in a unified domain. The firefly
algorithm (FA) and the grey wolf optimizer (GWO) are two of the most important
algorithms in computer science. Together, they form a new hybrid architecture
called the hybrid grey wolf optimizer (HGWO). First, a lower-order model of a
permanent magnet synchronous motor drive, comprising of speed and current
controllers, was obtained through a reduction using a signal-processing-based
identification technique. The anonymous proportional-integral (PI) controller
gains can be roughly estimated by matching the reduced system developed in cas-
cade with the controller to a desired system. The gains of the delta operator are
nearly replicated by the gains of the continuous-time controller. This allows for the
development of a unified controller for the drive. Thus, the proposed algorithm is
useful not only for order diminution but also for determining the control param-
eters of PMSM drives. In addition, many commonplace and industrial uses of ac
motor drives of this type can benefit in future research from this unified approach
study.
Keywords: Firefly algorithm (FA), grey wolf optimizer (GWO),
model reduction, controller design, reference model matching
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(159–172) © 2024 Scrivener Publishing LLC
159
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
160 Linear and Nonlinear System Modeling
adjustments in the exploration phase, while the grey wolf optimization meth-
ods is utilized as a local search mechanism to verify exploitation capabilities.
The method is activated with the aid of FA to kick off the quest. The calcula-
tion is repeated for a fixed number of times to find the optimal global posi-
tion in the necessary solution space. Therefore, the GWO algorithm begins
with the result that FA generated. Then, GWO is used as the search method
to accelerate the convergence to the optimal outcome. The combined solu-
tion of FA and GWO methods can successfully find an optimal value in this
fashion.
The particle swarm optimization (PSO) [25], the Harris hawks optimiza-
tion (HHO) [26], the arithmetic optimization algorithm (AOA) [27], the war
strategy algorithm (WSA) [28], and the wild horse optimization (WSO) [29]
are some of the methods used to compare and contrast with the proposed
method. For the sake of comparison, we also use a method that combines the
chaotic grey wolf approach with PSO. The prime motive behind model order
reduction approach is to carry out the control synthesis for a practical low
order PI controller in this chapter. The PID parameters are typically obtained
through auto-tuning or by using traditional old PID tuning rules [30]. In the
Truxal approach [19], which is used in precise model matching control [20],
the coefficients of the reference system are calculated to achieve the desired
outcome, and the controller gains are assessed in a way that is very similar
to the complete closed-loop system. The most significant drawback of this
exact model matching (EMM) method is that the controller’s implementa-
tion in hardware is not assured. To get around this limitation of the EMM
procedure, AMM [21] will be used instead of the current design schemes.
Therefore, the HGWO technique is used in conjunction with AMM tech-
nology to retrieve the controller’s coefficients. Similarly, a set of alternative
algorithms is used to compare with the HGWO method.
7.94 s 13.23
GHOS (s) 7 3
(10.1)
9.6 10 s 3.99 10 3 s 2 6.97 s 34.63
Ultimately, we want to create a PMSM drive controller that works well. Due
to the higher-order nature of the input-output connection produced from the
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Hybrid Grey Wolf Optimizer 165
available data, the first step is to simplify the model so that it can be used with a
normal PI controller. Due to its ability to reduce drive-related noise and inter-
ruption, PI controllers are frequently chosen over traditional PID controllers.
For a sampling rate of about 1 kHz, the discrete-delta domain transfer function
of the aforementioned model is denoted by Equation 10.1 as
When creating the reduced models in the delta domain, we assume pop-
ulation sizes of 20 and 100 and use a maximum of 100 iterations. All pub-
lished algorithms are considered, with their default parameters (proposed,
parent, and standard heuristics) taken into account. The unified domain’s
reduced systems are listed in Table 10.1. The minimum values for the fit-
ness function are also listed for your convenience. The optimal outcome is
also highlighted in italics in Table 10.1. The ancestor procedures and a few
other methods are used to draw comparisons.
Some well-known strategies in the literature have been compared with
the proposed HGWO method. This approach, as demonstrated in Table
10.1, has the lowest fitness value. Table 10.2 details the estimated transient
parameters for the simplified models. Results in the table with bolded
results are the most similar to the model’s parameters. In addition to the
methods used for comparing, the newly developed chaotic grey wolf opti-
mizer-particle swarm optimization (CGWOPSO) techniques are also uti-
lized in this table. Ten standard one-dimensional chaotic maps [33] are
used to develop CGWOPSO-1 to CGWOPSO-10 algorithms.
Once again, the integrated method outperforms the alternative algo-
rithms discussed in this section. In addition, the reduced-order models are
evaluated using the standard error metrics from systems and control theory.
The integral of absolute error (IAE) is calculated. It is first of all compared
with the ten variants of the hybrid technique CGWOPSO. Figures 10.1 and
10.2, respectively, show the comparison of HGWO with the different cha-
otic versions used with the hybrid CGWOPSO methods in two parts.
From Figures 10.1 and 10.2, it is clear that the IAE value for the pro-
posed HGWO method is quite less than as compared to CGWOPSO meth-
ods. Further, the parent techniques as well as standard methods are also
compared with HGWO in terms of IAE. The graph is shown in Figure 10.3.
The error corresponding to IAE with the suggested HGWO method
is not only less than the parent FA and GWO approaches but also found
less than some of the latest methods with which comparison is carried out
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
166 Linear and Nonlinear System Modeling
Table 10.1 Model reduction and the minimum cost function of the test system
in the unified domain of analysis.
Second-order models in the Minimum cost
Methods unified domain function
HGWO (proposed) 990.5 1436.0144 2.9563e–04
2
843.3 3759.2
in Figure 10.3. Other metrics of performance like ITAE, ISE, and ITSE
can also be tried out in future to provide the assessment of the proposed
method over the techniques available in the literature.
The final step now is to use the integrated delta domain to determine
the PI controller’s settings. For controller synthesis, the AMM method is
employed. The proposed algorithm allows for a heuristic determination of
the controller’s parameters; Table 10.3 displays the obtained results. Along
with determining the optimal values for the regulatory parameters, the
heuristic scheme also determines the respective fitness value as well. Once
again, bold letters denote the bare minimum fitness level.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Hybrid Grey Wolf Optimizer 167
Comparison of IAE
0.0059
Error values 0.0058
0.0057
0.0056
0.0055
0.0054
0.0053 O
-1
-2
-3
-4
-5
W
SO
SO
SO
SO
SO
HG
OP
OP
OP
OP
OP
W
W
CG
CG
CG
CG
CG
Methods
Figure 10.1 Comparison of IAE between HGWO and CGWOPSO 1–5 methods.
Comparison of IAE
0.0058
0.00575
0.0057
Error values
0.00565
0.0056
0.00555
0.0055
0.00545
0.0054
0.00535
O
-6
-7
-8
-9
0
-1
W
SO
SO
SO
SO
SO
HG
OP
OP
OP
OP
OP
W
W
CG
CG
CG
CG
CG
Methods
Figure 10.2 Comparison of IAE between HGWO and CGWOPSO 6–10 methods.
Comparison of IAE
0.006
0.0059
0.0058
Error values
0.0057
0.0056
0.0055
0.0054
0.0053
0.0052
HGWO FA [23] GWO PSO HHO AOA WSA WHO
[24] [25] [26] [27] [28] [29]
Methods
Figure 10.3 Comparison of IAE between HGWO and some standard heuristics.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Hybrid Grey Wolf Optimizer 169
Table 10.3 Parameters of the controller and the cost function values of the
reduced-order PMSM drive system in the combined domain.
Fitness function
Algorithms KP KI value
HGWO 1.8679 1.4967e–05 7.3715e–05
GWOPSO-1 0.0010 4.3879 7.6759e–05
GWOPSO-2 0.0041 2.1390 7.9023e–05
GWOPSO-3 0.0187 2.8646 7.3853e–05
GWOPSO-4 0.0009 4.3969 7.6776e–06
GWOPSO-5 0.0011 3.2887 7.7516e–05
GWOPSO-6 0.0019 4.3868 7.7428e–05
GWOPSO-7 0.0011 4.3923 7.6763e–05
GWOPSO-8 0.0027 4.3908 7.8690e–05
GWOPSO-9 0.0010 4.3888 7.6759e–05
GWOPSO-10 0.0014 4.3861 7.6850e–05
FA [23] 0.7172 3.3084e–04 7.5926e–05
GWO [24] 1.8012 1.3171e–04 7.388e–05
PSO [25] 0.8281 4.4387 7.7782e–05
HHO [26] 0.0013 4.3772 7.6833e–05
AOA [27] 0.0043 2.2451 7.5692e–05
WSA [28] 0.0025 4.4372 7.6717e–05
WSO [29] 0.0021 4.3328 7.5622e–05
10.4 Conclusions
The PMSM drive’s discrete delta domain controller architecture and order
reduction are discussed here. The PMSM drive’s reduced-order models are
enhanced using a constrained global optimization method. For this pur-
pose, a hybrid approach involving the GWO and FA methods was used
effectively. In addition, its PI controller’s parameters are derived using
ideas from approximate model matching. This technique is advantageous
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
170 Linear and Nonlinear System Modeling
because it relies solely on the feedback from the output. The control mecha-
nism was developed in such a way that it closely tracks the target model. In
this chapter, we also discussed one practical use of the combined domain
approach for modeling and control of the PMSM drive system. Similarly,
induction motor drives, which are a common type of drive, can have their
size minimized and a controller designed for them using the same kinds of
smart methods. Similar modeling and control methods can be applied to
other common drives, like switched reluctance motor (SRM) and brushless
direct current (BLDC) motor drive systems.
References
1. Patil Rahul, D., Saravana Prakash, P., Sandeep, N., & Verma, A. K. (2022).
Performance Evaluation of PMSM Drive Using ANFIS based Speed
Controller Under SVPWM Technique. In 2022 IEEE 10th Power India
International Conference (PIICON) (pp. 1-6).
2. Fatemimoghadam, A., Yan, Y., Iyer, L. V., & Kar, N. C. (2022). Permanent
Magnet Synchronous Motor Drive Using Deep-Neural-Network-Based
Vector Control for Electric Vehicle Applications. In 2022 International
Conference on Electrical Machines (ICEM) (pp. 2358-2364).
3. Kumar, A., Chandekar, A., Deshmukh, P. W., & Ugale, R. T. (2023).
Development of electric vehicle with permanent magnet synchronous motor
and its analysis with drive cycles in MATLAB/Simulink. Mater. Today Proc.,
72, 643-651.
4. Su, G., Wang, P., Guo, Y., Cheng, G., Wang, S., & Zhao, D. (2022).
Multiparameter identification of permanent magnet synchronous motor
based on model reference adaptive system—Simulated annealing particle
swarm optimization algorithm. Electronics, 11(1), 159.
5. Zhen, S., Li, Y., Liu, X., Wang, J., Chen, F., & Chen, X. (2023). A Lyapunov-
based robust control for permanent magnet synchronous motor in the mod-
ular joint of collaborative robot. Robotica, 41(5), 1389-1406.
6. Zhou, W., Hua, Q., Yao, Y., Kong, L., Xie, A., & Zhang, D. (2022). On-line
parameter identification and self-tuning PI controller for permanent magnet
synchronous motor. In 2022 2nd International Conference on Robotics and
Control Engineering (pp. 26-31).
7. Nguyen, T. P., Hsieh, M. F., Huynh, T. A., & Doan, Q. V. (2022). Direct
Flux-Vector Oriented Control for Interior Permanent Magnet Synchronous
Motor. In 2022 IEEE Ninth International Conference on Communications and
Electronics (ICCE) (pp. 501-506).
8. Huo, D., Yang, Z., Wang, Y., Wang, B., & Gong, C. (2023, March). Novel
Design Technique of Fuzzy Adaptive PI Regulator for Permanent Magnet
Synchronous Motor. In Proceedings of the 3rd International Symposium on
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Hybrid Grey Wolf Optimizer 171
22. Ganguli, S., Kaur, G., & Sarkar, P. (2020). A new hybrid algorithm for identi-
fication in the unified delta framework. In AIP Conference Proceedings (Vol.
2207, No. 1, p. 040002). AIP Publishing LLC.
23. Yang, X. S. (2010). Firefly algorithm, stochastic test functions and design
optimisation. Int. J. Bio-Inspired Comput., 2(2), 78-84.
24. Mirjalili, S., Mirjalili, S. M., & Lewis, A. (2014). Grey wolf optimizer. Adv.
Eng. Software, 69, 46-61.
25. Kennedy, J., & Eberhart, R. (1995, November). Particle swarm optimization.
In Proceedings of ICNN’95-International Conference on Neural Networks (Vol.
4, pp. 1942-1948).
26. Heidari, A. A., Mirjalili, S., Faris, H., Aljarah, I., Mafarja, M., & Chen, H.
(2019). Harris hawks optimization: Algorithm and applications. Future
Gener. Comput. Syst., 97, 849-872.
27. Abualigah, L., Diabat, A., Mirjalili, S., Abd Elaziz, M., & Gandomi, A. H.
(2021). The arithmetic optimization algorithm. Comput. Methods Appl.
Mech. Eng., 376, 113609.
28. Ayyarao, T. S., RamaKrishna, N. S. S., Elavarasan, R. M., Polumahanthi, N.,
Rambabu, M., Saini, G., ... & Alatas, B. (2022). War strategy optimization
algorithm: a new effective metaheuristic algorithm for global optimization.
IEEE Access, 10, 25073-25105.
29. Naruei, I., & Keynia, F. (2021). Wild horse optimizer: A new meta-heuris-
tic algorithm for solving engineering optimization problems. Eng. Comput.,
1-32.
30. O’Dwyer, A. (2009). Handbook of PI and PID controller tuning rules. Imperial
college press.
31. Ayyar, K. (2013). Investigations on model order reduction and design of speed
controller for electric drive systems.
32. Krishnan, R. (2001). Electric motor drives: modeling, analysis, and control
(Vol. 626). New Jersey: Prentice Hall.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
11
Parameter Estimation of First-Order RC
Model of Lithium-Ion Batteries in Electric
Vehicles Using Slime Mould Algorithm
Ramdutt Arya1, Shatrughan Modi2 and Souvik Ganguli3*
1
Centre for Energy & Environment, Malaviya National Institute of Technology,
Jaipur, Rajasthan, India
2
School of Computing, Indian Institute of Information Technology, Una,
Himachal Pradesh, India
3
Department of Electrical and Instrumentation Engineering,
Thapar Institute of Engineering and Technology, Patiala, Punjab, India
Abstract
Electric vehicles (EVs) are being considered as the promising solution for car-
bon-free transportation and the future of mobility. Battery electric vehicles (BEV)
are the true EVs that completely run from electricity supplied by their batteries
(Li-ion). For dynamic simulation of BEV and use of the battery effectively, it is
essential to model the battery and estimate its parameters accurately. Generally,
estimation of the battery parameters requires complex, time-consuming, and
expensive methods. This chapter focuses on a simple equivalent circuit methods
(ECM) is used for battery modeling. Estimation of battery parameters is done by
comparing the proposed model output to the known catalogue output with the
help of a newly formed heuristic optimizer-slime mould algorithm (SMA). Slime
mould algorithm is based on the morphological transformations, oscillation, and
foraging of slime mould found in nature, and have shown remarkable exploration
and exploitation capabilities. The evaluation of 31 parameters for first-order ECM
is done by taking the fitness function as the sum of squared-error of model voltage
and catalogue voltage, for charging as well as discharging scenarios. Additionally,
four well-known algorithms, ant lion optimizer (ALO), dragonfly algorithm (DA),
grey wolf optimizer (GWO), and sine-cosine algorithm (SCA), are examined for
comparison of performance and efficiency in terms of the parameters, accuracy
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(173–196) © 2024 Scrivener Publishing LLC
173
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
174 Linear and Nonlinear System Modeling
11.1 Introduction
The global environment is facing huge challenges. Air pollution and GHG
(Green House Gases) emissions due to the burning of fossil fuels are among
the major challenges. This burning of fossil fuel is done mainly in the
transportation, power, and industrial sectors. The transport sector obtains
90% of its energy need from burning of fossil fuels [1]. Invalid source spec-
ified. In order to decrease dependency on fossil fuels and replacing them
with clean and environmentally friendly solutions, the electrification of the
transport sector is the main focus of every nation. In this perspective, elec-
tric vehicles (EVs) have become the best choice to fulfil the need for a clean
source of transportation that produces very low or no carbon emission and
give enhanced fuel economy. This has increased the demand for EVs as
these are now being considered as the favorable transportation mode [2]
Invalid source specified.
Electric vehicles are those vehicles in which there is an electric drive or
vehicles in which propulsion power (either fully or partially) is received
from the electricity. There are mainly three types of EVs: hybrid electric
vehicle (HEV), plug-in hybrid electric vehicle (PHEV), and battery electric
vehicle (BEV). Among these, BEVs are considered as the true EV and fully
green mode of transport if the battery is charged from the green sources
only [3].
Electric vehicles’ success mainly depends on their battery/battery
pack. The battery serves as the energy source for the EV. This is a crucial
part of the EV, as most of the important aspects related to EV like range,
charging time, cost, weight, size, and safety are mainly dependent on it [4].
The battery pack consists of various cells packed and connected in a spe-
cific manner to achieve a certain voltage/power requirement. These cells
are made of materials that possess electrochemical properties. There are
various types of batteries available in the market based on materials like
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Parameter Estimation of First-Order RC Model 175
4y
R0 ( 1 2 x 3 x 2 )e ( 5 6 x 7 x2 ) (11.1)
11 y
R1 ( 8 9 x 10 x 2 )e ( 12 13 x 14 x2 ) (11.2)
18 y
C1 ( 15 16 x 17 x 2 )e ( 19 20 x 21 x2 ) (11.3)
25 y
V0 ( 22 23 x 24 x 2 )e
( 26 27 y 28 y2 29 y3) 30 x 31 x2 (11.4)
where α1, α2, α3, ……, α31 are coefficients of the equations formulated as a
mix of polynomials and exponential and V0 represents the terminal voltage
when no-load is connected (open circuit voltage).
Equations 11.1–11.5 are used to estimate the parameters and terminal
voltage of the battery when the load is connected to its terminals for charging
as well as discharging scenario for the first-order RC equivalent model.
During operations when current is kept constant, the time-dependent,
nonlinear expression for the battery terminal voltage is represented by the
following polynomial equation [14],
Q0
VM V0 IR1 e t / R1C1
I (R0 R1 ) (11.5)
C1
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
180 Linear and Nonlinear System Modeling
Here, Q0(0) is the nominal rated capacity of the battery at the time of
start and ‘t’ is the time duration since charging and discharging starts, both
these values are provided by the manufacturer.
In Equations 11.1–11.4, by replacing x with Cr and y with SoCcr for
charging scenario, and replacing with Dr and (1-DoDdr) for discharging
scenario, respectively, we can obtain the parameters for the battery model.
For discharging scenario, (1-DoDdr) is chosen, as with the increase in
DoDdr , V0 decreases [15].
Charge rate (Cr) and Dr can be controlled as per requirements of the
application. The Cr and Dr are kept constant for a specific BEV operation
and can be determined from the following equations,
I
Cr = (11.6)
Q0
I
Dr = (11.7)
Q0
The SOCcr and DODdr can be determined from the following equations,
I t
SoCcr SoCini (11.8)
Q0 .3600
I t
DoDdr DoDini (11.9)
Q0 .3600
Using the above equations, the model which is formed can be deployed
as a battery pack in any BEV simulation after estimation of polynomial
coefficients and battery parameters along with the battery charging and
discharging characteristics identification [15].
best fits the voltage-SoC characteristics curve taken from the manufactur-
er’s battery.
Here, the fitness is formulated in the form of sum of the squared error
between voltage values obtained from the model presented and catalogue
values at various sample data points for both chargings as well as discharg-
ing cases. This fitness function can be mathematically expressed as
n
F (x ) i 1 (Vi M ViC )2 (11.10)
with constraint as
0 1 ,0 2 ,0 3 ,. . . ..0 31
where Vi C denotes the ith voltage value taken from the catalogue and Vi M
is the ith value of model output voltage obtained from Equation 11.5 for the
first-order RC equivalent. The optimization is done with the help of SMAa
newly formed metaheuristic optimization technique at different SoC and
DoD values.
When a vein gets close to a food supply, the bio-oscillator sends out
a wave that spreads outward. This enhances the flow of cytoplasm inside
the vein. This rate of flow of cytoplasm determines the thickness of the
vein, i.e., fast flow will lead to a thick vein and vice versa. With the help of
this positive and negative feedback combination, the optimal path can be
established by the slime to connect the food source in better and efficient
way [22].
It is a three-stage process to discover an optimal path to the food. The
first one is searching and approaching the food that depicts the exploration,
oscillation, and wrapping of the food, which demonstrates the exploitation
characteristics of slime mould [18].
Xb (t ) vb. W . X A (t ) X B (t ) , r p
X (t 1) (11.11)
vc. X (t ), r p
vb [ a, a] (11.13)
t
a arctanh 1 (11.14)
max _ t
bF S(i)
1 r . log 1 , condition
bF wF
W Smellindex(i)
bF S(i)
1 r . log 1 , other
bF wF
(11.15)
here, condition tells that half of the population is ranked by S(i), r rep-
resents the random value in the [0,1] interval. max_t shows the maximum
number of iterations, bF and wF represents the best fitness and the worst
fitness value, respectively, obtained in the current iteration. SmellIndex
represents the sorted fitness values.
➢ Oscillation
The propagation wave, generated by the biological oscillator, redirects the
slime mould’s cytoplasm in its veins to a more favorable location in terms
of food density. Variations of the venous width of slime mould are being
simulated by W, vb, and vc. Synergistic characteristics between vb, and vc
resemble the selective nature of slime mould. Once it finds a good source
of food, it keeps trying to search for a better source of food. To achieve
a higher quality of food from the previous one, it separates some of its
organic matter to explore other locations, instead of investing all of it only
a single source of food [18].
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
184 Linear and Nonlinear System Modeling
11.4 Methodology
In this work, SMA is employed. This metaheuristic algorithm is used to
estimate the 31 polynomial coefficients of the first-order RC equivalent
circuit model of the battery, for the charging scenario as well as discharg-
ing scenario. It further helps to determine the modeled battery parameters
R0, R1, C1, and V0, which are nonlinear multivariable functions of SOCcr,
DODcr, Cr, and Dr. The values of Cr and Dr are taken as 1 and remains con-
stant for a particular BEV. The values of SOCcr and DODcr were taken from
the manufacturer datasheet. These modeled battery parameters are esti-
mated with four other well-known metaheuristic algorithms, viz., ALO,
GWO, DA, and SCA using Equations 11.1–11.4, which in turn calculate
the terminal voltage from Equation 11.5 and compared with catalogue val-
ues to determine the best fitness using the fitness function, Equation 11.19.
For the discharging scenario, we performed the abovementioned process
except that we used the value of x as Dr and y as (1-DODcr) because the
value of terminal voltage is decreasing with the increase in discharge.
The optimal parameters are estimated from the SMA and the other
four standard algorithms taken as reference. The SMA-estimated optimal
parameters are then evaluated with those obtained from alternative opti-
mization methods. The rate of convergence in estimating parameters and
fitness values from each algorithm is also determined. The fitness value
is defined as the sum of squared error. For the statistical analysis, we also
determined the best and worst fitness values, as well as the mean and stan-
dard deviation of these values.
Apart from the detailed statistical analysis, the obtained results are
were tested for significance. To check the significance, there are various
parametric and non-parametric tests. Parametric tests are used when
the data is uniformly distributed or follow some pattern, while non-
parametric tests are those in which there is no uniformity in the data
or the data is minimal [23]. In our case, the obtained results data do not
have any pattern.
We used the non-parametric Wilcoxon rank-sum test to obtain the
results of our first-order ECM for charging as well as discharging case. This
test provides the p-values for every algorithm compared with the proposed
algorithm [24]. Generally, the p-values, which are less than 0.05, confirm
that the results obtained by the proposed algorithm are significant in com-
parison with other algorithms [25]. This p-value can be decided by our
requirements to check the significance. Here, we have taken it as a standard
choice, i.e., p < 0.05.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Parameter Estimation of First-Order RC Model 185
Convergence curve
25
DA
GWO
ALO
SCA
20 SMA
Best Fitness value
15
10
0
0 50 100 150 200 250 300 350 400 450 500
Iteration
(i)
Convergence curve
45
DA
GWO
40 ALO
SCA
SMA
35
30
Best Fitness value
25
20
15
10
0
0 50 100 150 200 250 300 350 400 450 500
Iteration
(ii)
Figure 11.1 (i) Convergence curve for the first-order ECM charging. (ii) Convergence
curve for the first-order ECM discharging.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
190 Linear and Nonlinear System Modeling
Table 11.3 Fitness value analysis for the results of first-order RC ECM.
Charging Discharging
Standard Standard
Algorithm Best Worst Mean deviation Best Worst Mean deviation
SMA 1.25876 1.47741 1.34339 0.06224 0.08789 0.13355 0.10736 0.01183
DA 1.34443 1.75473 1.48143 0.09209 0.11789 2.57200 0.88574 0.65209
GWO 1.32021 1.53535 1.39020 0.05435 0.10643 0.19519 0.14297 0.02716
ALO 1.25498 1.41495 1.30158 0.04459 0.10831 0.79773 0.17523 0.14888
SCA 1.41466 2.40946 1.67595 0.27255 0.15419 6.62060 1.24497 1.67679
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Parameter Estimation of First-Order RC Model 191
Table 11.5 Modified p-values for Wilcoxon test for first-order RC ECM with
addition of Holm-Bonferroni corrections.
Corrected p – values
Proposed
Battery ECM method ALO DA GWO SCA
First-order SMA 2.5639E– 1.5497E– 3.3338E– 6.6283E–
charging 02 05 02 07
First-order SMA 6.9989E– 3.7039E– 2.3025E– 2.7182E–
discharging 06 07 05 07
0 10 20 30 40 50 60 70 80 90 100
2 groups have mean ranks significantly different from Group 1
(i)
Click on the group you want to test
0 10 20 30 40 50 60 70 80 90 100
4 groups have mean ranks significantly different from Group 1
(ii)
Figure 11.2 (i) Kruskal-Wallis test results for the First-order ECM charging. (ii) Kruskal-
Wallis test results for the first-order ECM discharging.
11.6 Conclusions
This chapter demonstrates the mathematical modeling of the first-order
RC ECM of battery which constitutes 31 parameters for charging as well
as discharging scenarios. The battery chosen is of Li-ion type, which is the
most suitable for EV applications. Then, the estimation of the parameters
of this model is done with the help of a newly formed metaheuristic tech-
nique named SMA. The results obtained from SMA are then contrasted
with four other standard metaheuristics techniques like DA, GWO, ALO,
and SCA. Simple statistical analysis like mean, standard deviation, and best
and worst fitness values are done. The compelling and statistical results
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Parameter Estimation of First-Order RC Model 193
demonstrate that SMA has the advantage over other algorithms in the
sense that the sum of square error is the least as compared to standard
algorithms. Slime mould algorithm (SMA) is useful as it gives the least
error and very good standard deviation. This tells the high reliability and
efficiency of SMA over other algorithms. From the convergence curves, we
see that there is not much significant advantage of SMA in terms of con-
vergence. The non-parametric assessments, viz., Kruskal-Wallis test and
Wilcoxon rank-sum test with added Holm-Bonferroni corrections, also
validate the better significance of SMA. We see that the p-values are quite
low as compared to the threshold value (0.05) taken. Overall, the results
obtained from SMA are good and acceptable in comparison to other algo-
rithms. The SMA is useful for the modeling of the Li-ion battery, especially
for EV applications where the range of computation time and accuracy is
large. Although the first-order ECM captures the battery behavior consid-
erably well, it is not able to furnish any useful information regarding the
lifetime of the battery. The performance is somewhat below average when
the number of nonlinearity increases. Since the convergence of SMA is
quite fast, we can include other critical battery parameters one by one and
assess the performance of this algorithm. Also, we can hybridize the SMA
with other good metaheuristic algorithms like ALO and GWO to yield bet-
ter results. The SMA is a newly formed heuristic algorithm that can further
be modified/improved for better results and suitability to various applica-
tions and superior performance.
References
1. F. Perera, “Pollution from Fossil-Fuel Combustion is the Leading
Environmental Threat to Global Pediatric Health and Equity: Solutions
Exist.,” Int. J. Environ. Res. Public Health, vol. 15, no. 1, December 2017.
2. S. V. Anand, “Global Environmental Issues,” Open Access Sci. Rep., January
2013.
3. J. Larminie and J. Lowry, Electric Vehicle Technology Explained, 2nd ed., John
Wiley & Sons Ltd, 2012.
4. K. V. Vidyanandan, “Batteries for Electric Vehicles,” November 2019.
5. C. Iclodean, B. Varga, N. Burnete, D. Cimerdean and B. Jurchiş, “Comparison
of Different Battery Types for Electric Vehicles,” IOP Conf. Ser.: Mater. Sci.
Eng., vol. 252, p. 012058, October 2017.
6. L. Gao, S. Liu and R. A. Dougal, “Dynamic lithium-ion Battery model for
system simulation,” Components and Packaging Technologies, IEEE Trans.
on, vol. 25, pp. 495-505, October 2002.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
194 Linear and Nonlinear System Modeling
22. T. Nakagaki, H. Yamada and T. Ueda, “Interaction between cell shape and
contraction pattern in the Physarum plasmodium.,” Biophys. Chem., vol. 84,
no. 3, pp. 195-204, May 2000.
23. S. S. MJ Campbell, “Parametric and Non-parametric tests for comparing
two or more groups,” [Online] https://www.healthknowledge.org.uk/public-
health-textbook/research-methods/1b-statistical-methods/parametric-non-
parametric-tests, 2020.
24. F. Wilcoxon, S. K. Katti and R. A. Wilcox, Critical values and probability levels
for the Wilcoxon rank sum test and the Wilcoxon signed rank test / [by], [Pearl
River, N.Y.]: American Cyanamid, 1963.
25. S. A. McLeod, “What a p-value tells you about statistical significance,” Simply
psychology, pp. 1-4, 2019.
26. Y. O. S. E. F. HOCHBERG, “A sharper Bonferroni procedure for multiple
tests of significance,” Biometrika, vol. 75, pp. 800-802, December 1988.
27. N. Breslow, “A Generalized Kruskal-Wallis Test for Comparing K Samples
Subject to Unequal Patterns of Censorship,” Biometrika, vol. 57, p. 579–594,
1970.
28. “EIG (2000),” [Online] http://www.eigbattery.com, 2000.
29. D. I. MacKenzie, J. D. Nichols, J. A. Royle, K. H. Pollock, L. L. Bailey and J.
E. Hines, “Chapter 3 - Fundamental Principals of Statistical Inference,” in
Occupancy Estimation and Modeling (Second Edition), Second Edition ed.,
D. I. MacKenzie, J. D. Nichols, J. A. Royle, K. H. Pollock, L. L. Bailey and J. E.
Hines, Eds., Boston, Academic Press, 2018, pp. 71-111.
30. J. M. Cimbala, “Hypothesis Testing,” Hypothesis Testing, pp. 1-5, 2014.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
12
Harmonic Distortion-Based Performance
Analysis and Fault Diagnosis of
Inverter Connected with BLDC
Motor Using Starting Transients
Surajit Chattopadhyay*, Chiranjit Sain, Purnendu Burui, Sk Rased Ali
and Soumya Saha
Abstract
This work deals with “harmonic distortion-based performance analysis and fault
diagnosis of inverter connected with brushless DC (BLDC) motor using starting
transients.” In the field of electrical engineering, advancement has been observed
almost in all fields like machinery, inverter technology, grid technology, and con-
trol methodology. Simultaneously, changes have been observed in material science,
power electronic devices, sensor technology, communications, etc. In this era, any
utility system depends on multi-disciplinary units and technologies. Therefore,
the performance of a device or unit, or system depends on not only itself but also
on various aspects. In this respect, the performance of the inverter also depends on
various aspects. In addition, faults may occur outside the system as well as inside
the system. It has a become a great challenge to monitor the performance of a sys-
tem and to detect faults effectively and accurately. Here, an attempt has been taken
to find some solutions for performance analysis and fault detection conveniently
and accurately. Inverters of various configurations consisting of different power
electronics devices connected with solar photovoltaic resources and BLDC motors
have been considered as a system of study. The output of the inverter unit has been
captured by the data acquisition system (DAS). Current signals are then analyzed
by the fast Fourier transform (FFT). The performance analysis and fault diagnosis
were carried out. A comparative study was performed that shows different results
Tamal Roy, Suman Lata Tripathi and Souvik Ganguli (eds.) Linear and Nonlinear System Modeling,
(197–214) © 2024 Scrivener Publishing LLC
197
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
198 Linear and Nonlinear System Modeling
for different Inverter configurations. The modeling and simulation performed for
various types of inverters like insulated-gate bipolar transistor (IGBT), and metal-
oxide-semiconductor field-effect transistor (MOSFET)-based inverter units have
shown satisfactory outcomes both in healthy and faulty conditions. Starting tran-
sient-based total harmonic distortion (THD) shows significant variations in fault
condition from the healthy condition in those inverters. Also, intermediate fault
conditions can be discriminated against THD values. In the end, major achieve-
ments and scopes for extension of this work have been presented.
Keywords: Brushless DC motor, fault diagnosis, harmonic distortion, inverter
switch, performance analysis, semiconductor switch, starting transients, total
harmonic distortion
12.1 Introduction
a. Need
In the last few decades, technology has advanced in several directions. In
the field of electrical engineering, advancement has been observed almost
in all fields like machinery, inverter technology, grid technology, and con-
trol methodology. Simultaneously, changes have been observed in material
science, power electronic devices, sensor technology, communications, etc.
In this era, any utility system depends on multi-disciplinary units and tech-
nologies. Therefore, the performance of a device or unit or system depends
on not only itself but also on various aspects. In this respect, the perfor-
mance of the inverter also depends on various aspects. In addition, faults
may occur outside the system as well as inside the system. It has become a
great challenge to monitor the performance of a system and to detect faults
effectively and accurately.
Along with the development of engineering and technology, various
mathematical tools have been introduced by different researchers in the
last few decades. Some of them have been observed very effective in ana-
lyzing images, signals, etc.
Therefore, it has become necessary to find some solutions to monitor
performance, detect faults, and cope with other various challenges.
b. Recent Research
Inverter technology and signal assessment have gone through tremendous
advancement. At the beginning of this work, a detailed literature survey
was carried out. Major observations are presented below.
M. Daniel Pradeep (2015) explained the voltage control method to con-
trol the speed of the brushless DC (BLDC) motor [1]. Abhishek Jain (2015)
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 199
C. Chhlonh et al. (2019) used the fuzzy controller to control the speed
of the rear wheel using BLDC [28]. A. Azarudeen et al. (2017) analyzed the
performance of the pulse width modulation (PWM) controller used for
the BLDC motor [29]. M. Ebadpour et al. (2017) introduced a model for
dual BLDC motor drives connected in parallel [30]. M. Nicola et al. (2019)
utilized a back electromotive force (EMF) observer to achieve sensorless
control of the motor [31]. A. Usman et al. (2019) made a detailed analysis
of faults that occurred due to the demagnetization effect [32].
N. Seth et al. (2017) made a detailed review on inverter topologies that
are useful for industries [33]. Tomoya Sugimoto et al. (2019) applied T
type boost inverter that gives better performance with the help of switch-
ing capacitors [34]. Tomoya Sugimoto et al. (2018) reduced switching loss
using multilevel topology with multiple numbers of inverters [35]. Pablo
Montero Robina et al. (2018) utilized a converter with a five-level diode
combination to interface with the solar grid [36].
S. Bernet et al. (1999) compared the performance of different types of
inverters [37]. H. Miyazaki et al. (2000) used IGBT-based clamped inverter
for industry applications [38]. H. A. Mantooth et al. (1997) made a sim-
ulation of a PWM inverter considering electro-thermal features [39]. A.
Ramamurthy et al. (1999) modeled voltage change during switching off the
condition for IGBT applications [40].
B. Chen et al. (2015) modeled MOSFET-based microinverter with high
efficiency [41]. O. Sivkov et al. (2018) compared silicon-based inverters
[42]. H. Kernstock et al. (2012) used soft switching methods to improve
efficiency [43]. S. Mangkalajarn et al. (2019) used class-E inverter to
achieve optimum performance [44].
Literature survey shows that a lot of works have been carried out to model,
design, and to analyze the performance of different converters used in drives’
applications. However, very few works have been noticed to utilize harmonics
patterns for performance analysis and fault diagnosis. This has motivated the
authors to find some solutions for performance analysis and fault detection
conveniently and accurately. The inverter of various configurations consisting
of different power electronics devices connected with solar photovoltaic (PV)
resource and BLDC motor have been considered as a system of study.
c. Challenges
The challenges and needs mentioned in the previous sections have moti-
vated the work team to find some solutions for performance analysis and
fault detection conveniently and accurately. The inverter of various config-
urations consisting of different power electronic devices connected with
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 201
12.2 Modeling
With the advancement of various types of inverter technology and the
advancement of electrical machinery, the use of inverter-connected motor
applications has increased to a great extent. Therefore, it has become a nec-
essary and urgent job to monitor the performance of the inverter-driven
motor unit.
In this work, an attempt has been taken to monitor the performance of
a solar PV-driven inverter that feeds power to a BLDC motor and to detect
short circuit faults inside the inverter unit.
Modeling was done by using MATLAB (R2018a) with different types
of inverter configurations. This analysis was done in two conditions, as
follows:
PowerGUI tool was used to carry out fast Fourier transform (FFT) anal-
ysis of the phase currents.
1. Energy unit
2. Converter unit
3. Motor unit
4. Data capture unit
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
202 Linear and Nonlinear System Modeling
I
N
V
E BLDC
DC Source R
T Motor
E
R
Capture Phase
Currents, RPM
Torque and
Analysis through
MATLAB
Figure 12.1 Block diagram of BLDC motor with hysteresis current control.
BATTERY
DATA
ACQUISITION
SYSTEM
ANALYSIS THE
SIGNAL USING
MATLAB
b. Energy Unit
The energy unit consists of a battery unit. The unit is charged by solar PV
resources through the charge controller.
c. Converter Unit
The converter unit consists of inverters to provide a three-phase supply to
the motor unit. Mainly, there are three types of inverter configurations that
were considered for modeling and simulation.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 203
Control
1-D T(u) iReq
iabc iabc
Scopes
S
R TL w
iabc R trq C
TL = f(w^2)
BLDC
Inverter
Figure 12.3 Modeling of BLDC motor with hysteresis current control (MATLAB 18).
a. IGBT Inverter
Total harmonic distortion (THD) for IGBT inverter output currents
drawn by the BLDC motor during starting conditions are presented in
Table 12.1. The THD values for different phases are compared, as shown
in Figure 12.4.
5000
4000
3000
2000
1000
0
t=0.0 to 0.5 t=0.5 to 1 t=0.0 to 1
THD of ia 518.7 6086.57 1756.69
THD of ib 538.24 6016.06 1841.55
THD of ic 531.7 6413.79 1803.47
Time(s)
THD of ia THD of ib THD of ic
8000
6000
4000
2000
0
t=0.0 to 0.5 t=0.5 to 1 t=0.0 to 1
THD of ia 553.09 11977.67 1690.01
THD of ib 549.47 14690.87 1530.75
THD of ic 525.41 7118.72 1682.7
Time
b. MOSFET Inverter
The THD for MOSFET inverter output currents drawn by the BLDC motor
during starting conditions is presented in Table 12.2. The THD values for
different phases are compared, as shown in Figure 12.5.
Table 12.3 THD for ideal semiconductor switch inverter in the healthy condition.
Duration t = 0.0 to 0.5 (s) t = 0.5 to 1 (s) t = 0.0 to 1 (s)
THD of ia 533.09 11977.01 1690.01
THD of ib 549.47 14,690.87 1530.75
THD of ic 525.41 7118.72 1682.7
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
206 Linear and Nonlinear System Modeling
8000
6000
4000
2000
0
t=0.0 to 0.5 t=0.5 to 1 t=0.0 to 1
THD of ia 533.09 11977.01 1690.01
THD of ib 549.47 14690.87 1530.75
THD of ic 525.41 7118.72 1682.7
Time
THD of ia THD of ib THD of ic
Figure 12.6 THD for ideal semiconductor switch inverter in the healthy condition.
530
520
510
500
Ideal Semiconductor
MOSFET IGBT
Switch
THD of ia 553.09 518.7 533.09
THD of ib 549.47 538.24 549.47
THD of ic 525.41 531.7 525.41
Time(s)
THD of ia THD of ib THD of ic
Figure 12.7 THD comparison of different types of inverters for t = 0.0 s to 0.5 s.
10000
8000
6000
4000
2000
0
Ideal
MOSFET IGBT Semiconductor
Switch
THD of ia 11977.67 6086.57 11977.01
THD of ib 14690.87 6016.06 14690.87
THD of ic 7118.72 6413.79 7118.72
Time(s)
THD of ia THD of ib THD of ic
1200
1000
800
600
400
200
0
Ideal
MOSFET IGBT Semiconductor
Switch
THD of ia 1690.01 1756.69 1690.01
THD of ib 1530.75 1841.55 1530.75
THD of ic 1682.7 1803.47 1682.7
Time(s)
THD of ia THD of ib THD of ic
1500
1000
500
0
Case1 Case2 Case3 Case4 Case5 Case6 Case7 Case8 Case9
THD of Ia 1594.95 1729.75 552.74 1875.4 2429.62 983.78 984.49 984.57 984.57
THD of Ib 1502.18 1680.88 633.55 1399.9 605.37 983.78 984.49 984.57 0
THD of Ic 1582.22 1897.72 614.55 1304.9 528.72 983.77 984.49 984.57 984.44
Figure 12.10 Variation of THD in the output current of inverter made by IGBT switches.
extreme faults are noted and presented in Figure 12.10. The figure clearly
shows that the fault in the power electronic switch is reflected in the varia-
tion of equivalent resistance that results in a significant change in THD val-
ues. The THDs are found slightly different in different phases at incipient
fault conditions. The THD becomes constant when switches get completely
damaged in cases 6, 7, 8, and 9.
1500
1000
500
0
Case 1 Case 2 Case 3 Case 4 Case 5 Case 6 Case 7 Case 8 Case 9
THD of Ia 1776.7 1656.6 565.59 2844.5 604.83 983.74 984.49 984.57 984.57
THD of Ib 1776.7 1656.6 565.59 2844.5 604.83 983.78 984.49 984.57 984.57
THD of Ic 1776.7 1656.6 565.59 2844.5 604.83 983.78 984.49 984.57 984.57
Figure 12.11 Variation of THD in the output current of inverter made by MOSFET switches.
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 209
1500
1000
500
0
Case 1 Case 2 Case 3 Case 4 Case 5 Case 6 Case 7 Case 8 Case 9
THD of Ia 1776.72 1656.6 565.59 2844.46 604.83 983.78 984.49 984.57 984.57
THD of Ib 1687.57 1879.42 592.65 2493.63 861.63 983.78 984.49 984.57 984.46
THD of Ic 1554.36 1703.49 530.91 1366.54 326.78 983.77 984.49 984.57 984.46
Figure 12.12 Variation of THD in the output current of Inverter made by ISS switches.
presented here. From these outcomes, the THD for different intermediate
and extreme faults are noted and presented in Figure 12.11. The figure clearly
shows that the fault in the power electronic switch is reflected in the varia-
tion of equivalent resistance that results in a significant change in THD val-
ues. Unlike inverters with IGBT, here, THD values are the same in all phases.
During completely damaged condition (equivalent to opening or breaking of
the switching path), THD becomes constant in the last four cases.
12.7 Conclusion
A comparative study performed in the previous sections shows different
results. The modeling and simulation performed for IGBT and MOSFET-
based inverter units showed satisfactory outcomes capable of differenti-
ating the healthy and faulty conditions. Starting transient-based THD
shows significant variations in fault conditions from the healthy condition
in those inverters. Also, intermediate fault conditions can be discrimi-
nated against using THD values. However, ideal switch-based inverter
units have shown satisfactory outcomes in healthy conditions and failed
to run the motor to provide necessary torque and speed to the mechan-
ical load during fault conditions. Hence, in this work, THD-based fault
analysis of the inverter with ideal switch could not be completed using
starting transient. Major achievements of this work are as follows:
(a) modeling of battery-inverter-BLDC motor system, (b) effective utili-
zation of various inverter configurations with different power electronic
switches, (c) capture of starting transient, (d) determination of harmonic
distortion and THD values of starting transients at healthy condition, (e)
determination of harmonic distortion and THD values of starting tran-
sients at faulty condition, (f) comparison of THD values, and (g) discrimi-
nation of healthy and faulty conditions. The work may further be extended
to carry out the following: (a) the work may be extended for performance
analysis and fault diagnosis of motor units, (b) motor other than BLDC
can be used in further study, (c) different energy resources can be used in
further study, (d) grid-connected inverter unit can be tested, (e) steady-
state and switching off transients can be analyzed, (f) other mathemati-
cal tools can be used to analyze current signatures, and (g) the work may
be extended for analysis and diagnosis from remote end using Internet of
Things and the Cloud.
References
1. M. Daniel Pradeep, “A Novel Method of Speed and Voltage Control of BLDC
Motor” IJCSMC, Vol. 4, Issue. 3, March 2015.
2. Abhishek Jain, “Controlling of Permanent Magnet Brushless Dc Motor
Using Instrumentation Technique” Int. J. Adv. Eng. Res. Dev. Volume 2, Issue
1, Jan-2015.
3. Yasser Ali Almatheel, “Speed Control of Dc Motor Using Fuzzy Logic
Controller” 2017 International Conference on Communication, Control,
Computing And
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 211
Four-Quadrant Operation,” in IEEE Trans. Ind. Inf., vol. 14, no. 1, pp. 146-
155, Jan. 2018, DOI: 10.1109/TII.2017.2723962.
17. G. Liu, S. Chen, S. Zheng and X. Song, “Sensorless Low-Current Start-Up
Strategy of 100-kW BLDC Motor with Small Inductance,” in IEEE Trans. Ind.
Inf., vol. 13, no. 3, pp. 1131-1140, June 2017, DOI: 10.1109/TII.2016.2607158.
18. S. Chen, W. Sun, K. Wang, G. Liu, and L. Zhu, “Sensorless High-Precision
Position Correction Strategy for a 100 kW@20 000 r/min BLDC Motor with
Low Stator Inductance,” in IEEE Trans. Ind. Inf., vol. 14, no. 10, pp. 4288-
4299, Oct. 2018, DOI: 10.1109/TII.2018.2793947.
19. R. Horvat, K. Jezernik and M. Čurkovič, “An Event-Driven Approach to the
Current Control of a BLDC Motor Using FPGA,” in IEEE Trans. Ind. Electron.,
vol. 61, no. 7, pp. 3719-3726, July 2014, DOI: 10.1109/TIE.2013.2276776.
20. X. Zhou, X. Chen, M. Lu and F. Zeng, “Rapid Self-Compensation Method
of Commutation Phase Error for Low- Inductance BLDC Motor,” in IEEE
Trans. Ind. Inf., vol. 13, no. 4, pp. 1833-1842, Aug. 2017, DOI: 10.1109/
TII.2017.2653812.
21. S. Neethu and V. Sreelekha, “Hysteresis controller based open phase fault
tolerant control of BLDC motor drives,” 2014 International Conference on
Power Signals Control and Computations (EPSCICON), 2014, pp. 1-6, DOI:
10.1109/EPSCICON.2014.6887514.
22. S. K. Chari, R. Dhiman and R. Saxena, “Novel and Robust Hysteresis Current
Control Strategies for a BLDC Motor: A Simulation Study and Inverter
Design,” 2018 2nd IEEE International Conference on Power Electronics,
Intelligent Control and Energy Systems (ICPEICES), 2018, pp. 841-846, DOI:
10.1109/ICPEICES.2018.8897424.
23. H. Suryoatmojo, F. G. Cladella, V. Lystianingrum, D. C. Riawan, R.
Mardiyanto, and D. Ariana, “Performance of BLDC Motor Speed Control
Based on Hysteresis Current Control Mechanism,” 2018 International
Seminar on Intelligent Technology and Its Applications (ISITIA), 2018, pp.
147-152, DOI: 10.1109/ISITIA.2018.8710910.
24. Y. K. Lee, “Torque Ripple and Switching Power Loss Minimization with
Constant Band Hysteresis Current Controller for BLDC Motor,” 2019 IEEE
PES Asia-Pacific Power and Energy Engineering Conference (APPEEC), 2019,
pp. 1-4, DOI: 10.1109/APPEEC45492.2019.8994737.
25. Ebin Joseph T and Sreethumol M V, “Speed control of BLDC motor drive
under Direct Torque Control scheme with modified integrator,” 2015 IEEE
International Conference on Electrical, Computer and Communication
Technologies (ICECCT), 2015, pp. 1-5, DOI: 10.1109/ICECCT.2015.7225931.
26. J. E. Muralidhar and P. V. Aranasi, “Torque ripple minimization & closed-
loop speed control of BLDC motor with hysteresis current controller,” 2014
2nd International Conference on Devices, Circuits and Systems (ICDCS), 2014,
pp. 1-7, DOI: 10.1109/ICDCSyst.2014.6926122.
27. M. S. Trivedi, R. K. Keshri and V. Ravikiran, “Predictive hysteresis current
control of PM BLDC motor in stationary plane,” 2017 IEEE Transportation
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Harmonic Distortion-Based Performance Analysis 213
Tamal Roy, PhD, received his PhD from Jadavpur University in 2016.
In 2008, he joined the Department of Electrical Engineering at Hooghly
Engineering and Technology College as a Lecturer with 15 years of aca-
demic experience. Since 2011, he has been working as an assistant profes-
sor in the Electrical Engineering Department of the MCKV Institute of
Engineering and presently is Head of the Department. His current research
interests include adaptive control, uncertainty modelling, and robust con-
trol of nonlinear systems.
215
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Index
217
Downloaded from https://onlinelibrary.wiley.com/doi/ by ibrahim ragab - Universita Di Firenze Sistema , Wiley Online Library on [09/09/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
218 Index
Green Energy: Solar Energy, Photovoltaics, and Smart Cities, Edited by Suman
Lata Tripathi and Sanjeevikumar Padmanaban, ISBN: 9781119760764.
Covering the concepts and fundamentals of green energy, this volume,
written and edited by a global team of experts, also goes into the practical
applications that can be utilized across multiple industries, for both the
engineer and the student.