Verma 2020
Verma 2020
ISA Transactions
journal homepage: www.elsevier.com/locate/isatrans
Research article
article info a b s t r a c t
Article history: The conventional maximum power point tracking (MPPT) algorithm shows best performance under
Received 3 May 2019 uniform insolation but when photovoltaic (PV) array is partially irradiated, the Power vs Voltage (P–
Received in revised form 21 November 2019 V) plot consists of multiple local maxima power point (LMPP) and one global maxima power point
Accepted 3 January 2020
(GMPP). The conventional MPPT algorithm may track local peak and fluctuate around it, resulting
Available online 11 January 2020
in lower power tracking. To eradicate this drawback of conventional algorithm, the solar PV system
Keywords: requires the synthesis of modified controller which is able to discriminate between local and global
Photovoltaic system peak point. Along with implementing modified MPPT controller, to minimise the adverse effect of
Partial shading condition partial shading on PV system, different PV array arrangements like series-parallel (SP), honey comb
Perturb & observe (HC), total cross tied (TCT) etc. may be used. Author(s) in the present study, has proposed asymmetrical
Type-1 FLC interval type-2 fuzzy logic control (IT-2 AFLC) based MPP algorithm for tracking global peak in partial
IT-2 FLC shading condition (PSC) with different PV array arrangements. The presented algorithm has been
compared with other approaches viz. perturb & observe (P&O) and type-1(T-1) FLC for GMPP tracking,
fill factor, shading losses, mismatch loss and efficiency to establish its superiority. For evaluating the
efficiency of different algorithms, the EN50530 MPPT efficiency test has been performed under dynamic
condition. The proposed algorithm has been developed under MATLAB/Simulink environment.
© 2020 ISA. Published by Elsevier Ltd. All rights reserved.
https://doi.org/10.1016/j.isatra.2020.01.009
0019-0578/© 2020 ISA. Published by Elsevier Ltd. All rights reserved.
252 P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263
electrical connection. Mismatch loss caused by partial shading are does not require professional expertise in comparison with above
minimised and generated power is enhanced. Braun et al. [14] intelligent control techniques. T-2 FLC was introduced which is
optimised the PV array rearrangement for different weather con- capable to cover all phases of its design. Unlike type-1 fuzzy, in-
ditions. Overall, 4%–5% gain in power was achieved using op- terval type 2 fuzzy logic controller is represented by two bounded
timisation technique under partial shading. This technique suf- membership functions, a lower and upper membership function.
fers from poor performance when modules exceed the size of Footprint of uncertainty (FOU) is the new dimension of IT-2
adaptive bank under huge and moderate shading. FLC. Interval type-2 FLC has output processing stage which has
Various literature is available on MPPT algorithms but, most type-reducer followed by defuzzifier. In [40,41] Altin proposes
of the presented algorithm works well under uniform weather interval T-2 FLC based MPPT algorithm for PV system but the
condition, as PV system exhibits only one peak point under such presented techniques were not tested for SPV system under par-
condition. But under partial shading condition or changing en- tially shaded conditions. In [42] authors have developed adaptive
vironmental condition, PV array output exhibits multiple peaks estimator for unstable output error system. The measurement of
i.e. LMPP and GMPP. MPPT algorithm may lurch around LMPP. The noisy system has been used to approximate the optimal control
conventional MPPT algorithms like P&O, Fractional Open Voltage policy. In [43] authors have designed model free adaptive control
(FOCV), Incremental Conductance (I&C) etc. were tested under which can take into accounts internal-external parametric and
uniform shading condition in [15–19] although these techniques non-parametric uncertainties.
were easy to understand and implement but the conventional In this paper, the authors have proposed a novel, intelli-
algorithms may lags behind in tracking global peak under PSC’s gent asymmetrical interval type-2 FLC based MPPT algorithm
and also show ripple. Advancement in control strategy leads to for stand-alone PV system under PSC. The conventional MPPT
the better tracking of GMPP [20–25]. Many MPPT techniques algorithm works on P&O steps and uses hill-incrementing method
based on intelligent, adaptive, genetic algorithms were intro- for successive iterations. These techniques continuously compare
duced in research articles. Authors in [26] presented new control present and previous power and voltage values. The conventional
and optimisation techniques such as sliding mode control, par- algorithm iteration step is based on pre-determined voltage step
ticle swarm optimisation, predictive control, cuckoo search algo- size ∆V. So, when conventional algorithms reach first peak of PV
rithm, etc. for control and better performance of grid connected array i.e. mostly local peak, it lurches around and stops searching
PV system. The presented strategies not only improve the effi- in the forward direction. Conventional methods tracks global
ciency but also reduces the power generated cost. A. Menadi [27] peak power under uniform irradiation as PV array depicts only
implemented fuzzy logic controller based MPPT algorithm on one peak [44]. But under other conditions it track local peak point
PV system connected to grid. The MATLAB/Simulink results are with less power. This also reduces the efficiency of PV system.
conformed with hardware results that shows MPPT algorithm This drawback of conventional algorithms can be overcome by the
tracks maximum power but the MPPT performance has been intelligent asymmetrical IT-2 FLC based MPPT algorithm. Mem-
tested under STC, weather change has not been considered in bership functions of the proposed asymmetrical IT-2 FLC based
the presented studies. Fuzzy Logic and P&O based MPPT track- algorithm are optimised to track global peak under transient
ing algorithm was presented in [28]. The presented technique conditions. Performance assessment of proposed asymmetrical
takes advantages of both the algorithms, utilising slow and fast interval type-2 FLC with other approaches viz. P&O, T-1 FLC in
change in insolation of P&O and fast-tracking ability of fuzzy terms of GMPP tracking, shading losses, mismatch loss, fill factor,
but it suffers from drift problem under rapidly changing irradi- EN50530 MPPT efficiency test has been performed under dynamic
ance level. In [29] applications of enhanced versions of FLC viz. condition establish its superiority.
modified hill climbing and adaptive P&O fuzzy logic based maxi- Although, the interval type-2 fuzzy based MPPT for PV sys-
mum power point tracking has been presented by the authors. tem under uniform insolation has already been implemented
All the improved version of FLC’s were compared for simula- but asymmetrical interval type-2 FLC under PSC with various
tion and hardware results. Asymmetrical type-1 FLC based MPPT PV array arrangements is still not available in literature. The
algorithm was implemented in [30] under both uniform and non- technique shall be helpful for the system designers for designing
uniform insolation condition . The implemented asymmetrical stand-alone solar PV system.
FLC improves the steady state as well as transient response of A 1.28 kW stand-alone solar PV system has been designed and
the system. In [31] authors presented novel Beta-Parameter based simulated for the proposed studies. Fig. 1 shows the schematic
type-1 FLC where a third input, a variable beta parameter is diagram of the proposed stand-alone PV system under PSC. The
introduced. The presented MPPT technique has been used for system is composed of PV array, boost converter, MPPT con-
rapidly changing environmental situations i.e. transient condition troller to locate highest peak and load. The methodology involves
using T-1 FLC. It has been found in the literature [29–31] that Simulink modelling of stand-alone PV system with various PV
the authors have used type 1 FLC for GMPP tracking under PSC array arrangements. SimPower/Simulink software is used for the
but the different PV array arrangement were not considered simulation studies. The mathematical modelling and design of the
and MPPT efficiency test under transient condition has not been components of proposed system are explained in the following
performed. The algorithm in [32] presented adaptive neuro fuzzy sections.
inference system based MPPT algorithm. The ANFIS model for
tracking MPP takes temperature and insolation as input and duty 2. PV array
cycle as the output. Further, artificial neural network [33], genetic
algorithm [34], grey wolf optimisation [35], ant-colony optimi- PV cell, is a basic unit of PV array, acts as an electrical trans-
sation [36], simulated annealing [37] were introduced for MPPT. ducer which converts light (photon) energy into electrical energy.
But these methods of MPPT need expert knowledge and huge The ideal solar cell circuit have current source with an anti-
data [31,38,39]. Conventional methods, called hard computing parallel diode. Shunt and series resistance are incorporated due
methods, need readjustment of controller parameters frequently to ohmic resistances and leakage current [45]. Fig. 2 presents
and cannot take into account the uncertainty and vagueness in circuit diagram of one diode model of PV cell with shunt and
measurands. If the data have uncertainties as in case of partial series resistance to make practical PV cell. PV cells are connected
shading condition, interval type-2 fuzzy method provides bet- in series/parallel combination to make PV modules, which are
ter solutions. Interval type-2 fuzzy logic based MPPT algorithm further connected to make PV array of desired output.
P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263 253
Table 1
Maximum Power obtained from 4 × 4 PV array under various arrangements.
PV array arrangement Pnom (W) Vm (V) Im (A)
S 1280 283.2 4.52
SP 1280 70.8 18.08
TCT 1280 70.8 18.08
where, Io/p and Vo/p is the output current and voltage of PV array,
Isat indicates reverse saturation current of diode, s and p are the
number of cell in series and parallel, q is the elementary charge
(1.602*10−19 C), A indicates ideal factor, K stands for Boltzmann
Fig. 2. Equivalent circuit of one diode model of PV cell.
constant (1.38*10−23 J), Tac indicates actual temperature.
In the present study, various 4 × 4 PV array arrangement viz.
S, SP and TCT of 1.28 kW has been considered. The specification
The output current of solar PV cell Io/p(cell) is given by using of the PV array is given in Appendix 1. Fig. 3 shows the various
Kirchhoff’s law: arrangements of PV array.
In Table 1, the maximum power obtained for the various PV
Io/p(PVCell) = Iphc − Id − Ish (1) array arrangements are listed.
[ ( ( )) ]
q Vo/p(cell) + Io/p(cell) ∗ Rs 3. Boost converter
Io/p(PVcell) = Iphc − Isat ∗ exp −1
AkTac
Boost converter converts variable dc from PV array into stable
Vo/p(cell) + Io/p(cell) Rs fixed dc. The boost converter circuit has IGBT switch, inductor,
− (2)
Rsh shunt capacitor and diode. Fig. 4 shows the circuit diagram of dc–
dc boost converter. The values of inductance, capacitance, duty
The output current for the PV array Io/p is given as ratio and resistive load can be obtained as:
( )
[ ( ( )) ] Vi/p ∗ Vo/p − Vi/p
q Vo/p + Io/p ∗ Rs (s/p) L= (4)
Io/p = p.Iphc − p.Isat ∗ exp −1 (∆I fsw Vo/p )
sAkTac
(3)
Vo/p + Io/p Rs (s/p)
( )
Vi/p
− α =1− (5)
Rsh (s/p) Vo/p
Fig. 5. (a) Type-1 FLC and (b) Blurring of type-1 membership function.
Fig. 6. Proposed Interval type-2 fuzzy logic based MPPT algorithm under PSC.
P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263 255
Fig. 7. Input variable (a) Error membership functions (b) Change in Error membership functions.
Fig. 10. Power vs voltage plot under various shading scenarios. 5. Simulation results
Fig. 11. (a)–(c) Power vs. Time plot for various MPPT techniques with different PV array arrangements for case 1 (STC).
P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263 257
Fig. 12. Comparative Analysis of Maximum power tracked using various MPPT techniques with different PV array arrangements for case 1(STC).
Table 3 Table 6
Shading loss (W) in various MPPT techniques under study with different PV Efficiency (%) of various MPPT techniques under study with different PV array
array arrangements. arrangements.
MPPT techniques Series Series-Parallel Total cross tied MPPT techniques Series Series-Parallel Total cross tied
P&O 302 297 253 P&O 77.90 78.32 81.83
Case 1 T-1 287 254.5 252 Case 1 T-1 79.12 81.71 81.91
IT-2 282 252 252 IT-2 79.52 81.91 81.91
P&O 342 345 321.5 P&O 88.99 88.70 90.93
Case 2 T-1 339 328 310 Case 2 T-1 89.21 90.23 92.03
IT-2 330 326 305 IT-2 90.13 90.51 92.50
P&O 766 762 763.5 P&O 86.82 87.50 87.24
Case 3 T-1 765 761.5 748.5 Case 3 T-1 86.94 87.58 89.78
IT-2 752 747 716.5 IT-2 89.18 90.33 95.18
P&O 574 561 453 P&O 82.38 83.89 96.49
Case 4 T-1 574 559 439 Case 4 T-1 82.38 84.13 98.13
IT-2 555 555 435 IT-2 84.39 84.59 98.59
Table 4
type-2 fuzzy algorithm under TCT arrangement. On the other side,
Fill Factor in various MPPT techniques under study with different PV array
arrangements.
series PV array arrangement gives the lowest maximum power.
MPPT techniques Series Series-Parallel Total cross tied
To express these results more clearly, the same results has been
presented on bar chart in Fig. 12 for case 1.
P&O 0.558 0.561 0.586
Case 1 T-1 0.566 0.585 0.586 Case 2 (Left corner shading): In this case, the system has been
IT-2 0.569 0.586 0586 tested with left corner shading with different PV array arrange-
P&O 0.535 0.533 0.547 ments. The PV panels were irradiated with 650 W/m2 at left
Case 2 T-1 0537 0.543 0.553
corner. In Fig. 13, the IT-2 FLC algorithm has been compared with
IT-2 0.542 0.544 0.556
other algorithms for different PV array arrangements under left
P&O 0.301 0.310 0.309
corner shading. In this case, the global maximum peak power
Case 3 T-1 0.301 0.311 0.319
IT-2 0.309 0.319 0.338 i.e. 975 W has been tracked by interval type-2 fuzzy algorithm
under TCT arrangement. To express these results more clearly,
P&O 0.412 0.426 0.491
Case 4 T-1 0.412 0.428 0.499 the same results has been presented on bar chart in Fig. 14 for
IT-2 0.423 0.430 0.501 case 2.
Case 3 (Top shading): In this case, the system has been tested
with top PV panels shaded with different PV array arrangements.
Table 5 The PV panels were irradiated with 400 W/m2 at top. In Fig. 15,
Mismatch losses (%) in various MPPT techniques under study with different PV
the IT-2 FLC has been compared with other algorithms for dif-
array arrangements.
ferent PV array arrangements under top shading. In this case,
MPPT techniques Series Series-Parallel Total cross tied
the global maximum peak power i.e. 563.5 W has been tracked
P&O 22.07 21.6 18.16
by interval type-2 fuzzy algorithm under TCT arrangement. On
Case 1 T-1 20.8 18.2 18.0
IT-2 20.4 18.0 18.0 the other side, series PV array arrangement gives the lowest
maximum power. To express these results more clearly, the same
P&O 11.0 11.2 9.06
Case 2 T-1 10.72 9.67 7.7 results has been presented on bar chart in Fig. 16 for case 3.
IT-2 9.86 9.48 7.4
Case 4 (Patch Shading): In this case, the system has been tested
P&O 13.1 12.5 12.7 with patch shading with different PV array arrangements. The PV
Case 3 T-1 13.0 12.0 10.2
panels were irradiated with 300 W/m2 at centre. In Figs. 15 and
IT-2 10.8 9.9 4.8
17, the IT-2 FLC has been compared with other algorithms for
P&O 17.6 16.1 3.5
different PV array arrangements under patch shading. In this case,
Case 4 T-1 17.6 15.8 1.8
IT-2 15.4 15.4 1.4 the global maximum peak power i.e. 845 W has been tracked
by interval type-2 fuzzy algorithm under TCT arrangement. On
the other side, series PV array arrangement gives the lowest
maximum power tracked by P&O. To express these results more
clearly, the same results has been presented on bar chart in Fig. 18
for case 4.
258 P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263
Fig. 13. (a)–(c) Power vs. Time plot for various MPPT techniques with different PV array arrangements for case 2 (Left corner shading).
Fig. 14. Comparative Analysis of Maximum power tracked using various MPPT techniques with different PV array arrangements for case 2 (Left corner shading).
Further, simulation results of the proposed MPPT algorithm the given PV array. It can be given by:
has been compared with other algorithms viz. conventional P&O
PMP = VMP ∗ IMP (14)
and Intelligent T-1 FLC under various shading scenarios along
with different PV arrangements for GMPP tracking, shading (b) Shading Losses: The power loss due to shading is called
losses, mismatch loss, fill factor and the EN 50530 MPPT efficiency shading loss. Shading loss is the difference in power be-
test. tween the maximum power obtained from an array un-
der STC (PMP ,withoutshading ) and the total maximum available
(a) GMPP Tracking: Under partially shaded condition, the non- power under PSC (PMP ,shading ) [48]. It can be represented as:
linear P–V characteristic exhibits multiple peaks. These
multiple peaks contain many LMPPs and single GMPP. The
point at which voltage and current value is maximum for PMP ,shadinglosses = PMP ,withoutshading − PMP ,shading (15)
P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263 259
Fig. 15. (a)–(c) Power vs. Time plot for various MPPT techniques with different PV array arrangements for case 3 (Top shading).
Fig. 16. Comparative Analysis of Maximum power tracked using various MPPT techniques with different PV array arrangements for case 3 (Top shading).
(c) Mismatch Losses: It can be defined as, the difference of value of fill factor is equivalent to unity, then the system is
individual maximum power (Pmpi ) of PV array and global more efficient.
maximum power (PGMPP ). Table 2 shows the individual PGMPP
maximum powers under different shading scenarios. Mis- FF = (17)
Voc ∗ Isc
match loss (Pmml ) is given by:
(e) MPPT Efficiency test EN50530: Under dynamic weather
Pmml = Pmpi − PGMPP (16)
condition, the assessment of the proposed algorithm is
(d) Fill Factor (FF): It can be defined as the ratio of maximum done with EN 50530 MPPT efficiency test. EN 50530 effi-
global power (PGMPP ) to the product of open circuit voltage ciency test is to supply the solar insolation in trapezoidal
(Voc )and short circuit current (Isc ) of the PV system. If the waveform with different ramp inclination. MPPT average
260 P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263
Fig. 17. (a)–(c) Power vs. Time plot for various MPPT techniques with different PV array arrangements for case 4 (Patch Shading).
Fig. 18. Comparative Analysis of Maximum power tracked using various MPPT techniques with different PV array arrangements for case 4 (Patch Shading).
tracking efficiency of the proposed controller is evaluated be seen that the performance of TCT PV array arrangement is
using the given formula [49]: better than any other configuration. Also, the proposed asymmet-
rical interval type-2 FLC algorithm has minimum shading losses,
Po/pPV improved fill factor and minimum mismatch losses as compared
ηMPPT = (18)
to the other algorithm. Table 6 shows the efficiency of proposed
Ptheo PV
MPPT algorithm compared with other algorithms. From Table 6, it
Tables 3–5 represents shading loss, fill factor and mismatch can be seen that IT-2 based AFLC with TCT PV array arrangement
loss (%) respectively in various MPPT algorithms under study with has highest efficiency.
Dynamic Results: The proposed asymmetrical interval type-2
different PV array arrangements. The proposed algorithm has
FLC MPPT algorithm has been tested under dynamic condition.
been compared with other control algorithm in terms of tracked Authors(s) have tested the proposed algorithm with EN 50530
GMPP, shading losses, fill factor and mismatch losses. Also, the MPPT efficiency standard. Fig. 19 shows two trapezoidal wave-
efficiency of the proposed algorithm has been evaluated. It can form representing slow and fast change at insolation level of
P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263 261
Fig. 19. Dynamic change in insolation level for EN 50530 standard MPPT efficiency test.
Fig. 20. Power vs. Time plot for various MPPT techniques with different PV array arrangements under dynamic condition.
400 W/m2 /s and 800 W/m2 /s respectively. The insolation level Table 7 shows the dynamic response of various MPPT tech-
changes slowly from 1–6 s and the insolation level change is niques. The proposed MPPT algorithm tracks the maximum global
faster from 7–10 s. power i.e. 1028 W (1–6 s) when PV array is configured in total
Fig. 20 shows the Power Vs. Time plot of the proposed algo- cross tied manner which is highest in comparison with other
rithm compared with other algorithms under dynamic condition. PV array arrangements. Also, the settling time is minimum for
Fig. 20(a), (b) and (c) shows the comparative analysis when PV the proposed algorithm. Hence, it can be concluded that TCT
array are connected in S, SP and TCT respectively. configuration gives better results.
262 P. Verma, R. Garg and P. Mahajan / ISA Transactions 100 (2020) 251–263
Table 7
Dynamic response of various MPPT techniques under study.
0–1 (s) 1–06 (s)
PV array arrangement MPPT techniques
Power (W) Settling time (s) Power (W) Settling time (s)
P&O 84 0.62 975 3.33
SERIES T-1 86 0.60 993 3.33
IT-2 86 0.60 998 3.33
P&O 85 0.62 983 3.4
SERIES-PARALLEL
T-1 96 0.20 1025 3.15
IT-2 96 0.20 1028 3.15
P&O 85 0.55 1022 3.26
TOTAL CROSS TIED
T-1 96 0.20 1024 3.10
IT-2 96 0.20 1028 3.10
[31] Li X, Wen H, Hu Y, Jiang L. A novel beta parameter based fuzzy- [40] Altin N. Interval type-2 fuzzy logic controller based maximum power point
logic controller for photovoltaic MPPT application. Renew Energy tracking in photovoltaic systems. Adv Electr Comput Eng 2013;13:65–70.
2019;130:416–27. [41] Altin N. The type-2 fuzzy logic controller-based maximum power point
[32] Kharb RK, Shimi SL, Chatterji S, Ansari MF. Modeling of solar PV module tracking algorithm and the quadratic boost converter for Pv system. J
and maximum power point tracking using ANFIS. Renew Sustain Energy Electron Mater 2018;47:4475–85.
Rev 2014;33:602–12. [42] Tutsoy O, Colak S. Adaptive estimator design for unstable output error sys-
[33] Rizzo SA, Scelba G. ANN based MPPT method for rapidly variable shading tems: A test problem and traditional system identification based analysis
conditions. Appl Energy 2015;145:124–32. 2015;229:902–16.
[34] Daraban S, Petreus D, Morel C. A novel MPPT (maximum power point [43] Tutsoy O, Erol D, Tugal H. Design of a completely model free adaptive
tracking) algorithm based on a modified genetic algorithm specialized on control in the presence of parametric, non-parametric uncertainties and
tracking the global maximum power point in photovoltaic systems affected random control signal delay. ISA Trans 2018;76:67–77.
by partial shading. Energy 2014;74:374–88. [44] Miyatake M, Veerachary M, Toriumi F, Fujii N, Ko H. Maximum power
[35] Mohanty S, Subudhi B, Member S, Ray PK. A new MPPT design using grey point tracking of multiple photovoltaic arrays: A PSO approach. IEEE Trans
wolf optimization technique for photovoltaic system under partial shading Aerosp Electron Syst 2011;47:367–80.
conditions 2015;7:1–8. [45] Bellia H, Youcef R, Fatima M. A detailed modeling of photovoltaic module
[36] Sundareswaran K, Vigneshkumar V, Sankar P, Simon SP, Srinivasa Rao using MATLAB. NRIAG J Astron Geophys 2014;3:53–61.
Nayak P, Palani S. Development of an improved P & O algorithm assisted [46] Zadeh LA. Fuzzy sets. Inf Control 1965;8:338–53.
through a colony of foraging ants for MPPT in PV system. IEEE Trans Ind [47] Zadeh LA. The concept of a linguistic variable and its application to
Informatics 2016;12:187–200. approximate reasoning-I. Inform Sci 1975;8:199–249.
[37] Lyden S, Haque ME. A simulated annealing global maximum power point [48] Bingöl O, Özkaya B. Analysis and comparison of different PV array config-
tracking approach for PV modules under partial shading conditions. IEEE urations under partial shading conditions. Sol Energy 2018;160:336–43.
Trans Power Electron 2016;31:4171–81. [49] Mara W, Piotrowicz M. Calculation of dynamic MPP-tracking efficiency of
[38] Seyedmahmoudian M, Horan B, Rahmani R, Oo AMT, Stojcevski A. Effi- PV-inverter using recorded irradiance. In: Proc 20th int conf mix des integr
cient photovoltaic system maximum power point tracking using a new circuits syst - mix 2013. 2013. p. 431–4.
technique. Energies 2016;9:1–18.
[39] Gosumbonggot J, Fujita G. Partial shading detection and global maximum
power point tracking algorithm for photovoltaic with the variation of
irradiation and temperature. Energies 2019;12:202.