Abstract

High solar photovoltaic (PV) penetration in the electrical grid can result in undesired effects on the voltage quality, leading to line loss and voltage magnitude increases. One of the main criteria to ensure the safe penetration of high-power solar systems in the main grid is maintaining an acceptable voltage magnitude when a disturbance occurs (e.g., 0.95 and 1.05 per unit) with respect to total installed power generation capacity of PV power plants. This manuscript analyzes the effects of high solar PV penetration per unit of voltage stability using the Fast Voltage Stability Index and total power loss. Moreover, we investigate the flexibility benefits of coordinated voltage control based on a smart inverter of solar PV capacitor banks (SI-CBs) under five cases in a typical microgrid (MG) test model. For the test of the SI-CBs, MG modeling is developed on a modified IEEE 123 test feeder, which includes 11 building load solar PVs with smart inverters and capacitor banks with real-time data from an area in Los Angeles, California, USA. The simulation results are presented to validate the effectiveness of the proposed approach using a real-time MATLAB interface to the Open Distribution System Simulator (OpenDSS).

1. Introduction

Microgrids (MGs) are emerging as critical tools for resilience with distributed energy resources, which have been increasing in the US grid. As one of the world’s largest and most complex human-made systems, MGs are highly vulnerable to three types of external threats: (1) natural disasters, (2) physical attacks, and (3) cyberattacks [13]. Therefore, extensive studies have been conducted on MGs to ensure their stability and reliability in supplying electrical power with many intermittent renewable energy resources [410]. In Ref. [4], the Bronzeville neighborhood of Chicago is considered to present a way to achieve MG integration into the existing distribution management system based on advanced design, control, and test approaches. In Ref. [5], a highly fluctuating load case is considered for modeling an MG power generation cost function, where a real-time monitoring operation is proposed considering practical cases of wind speed, solar irradiance, and battery charge states. Thus, resonance frequency monitoring was the subject of Ref. [7], in which a novel algorithm was proposed based on the wavelet transform, and several cases were considered, such as MG islanded and connected modes.

According to IEEE Standards 1547 and 929 [11, 12], some measures should be considered to ensure the reliability and stability of the grid-tied MG operations and allow integrating more distributed generation (DG)-based renewable energy resource MGs. As a source of renewable energy resources in MGs, solar photovoltaics (PVs) can steadily provide power when the grid is down and reduce energy costs when the grid is available in the US [13]. In addition, the sizing of MG components, such as the PV, battery, and generator capacities, is a paramount issue, as presented in Ref. [14]. The main advantages of solar PV grid integration include the potential of improved power quality and grid reliability and reduced transmission and distribution loss [15, 16], as there is a growing demand for alternative clean energy vs. fossil fuel. However, the integration of solar PV in weak distribution networks can often cause overvoltage problems during higher-level penetration scenarios. Traditionally, overvoltage problems have been mitigated using passive methods, such as network reinforcement, leading to high DG interconnection costs [17]. Furthermore, higher-level solar PV penetration scenarios have been becoming increasingly common in studies because of the potential to affect the life cycle of the operation of the distribution feeder equipment and the increased maintenance costs [18].

Therefore, voltage problems, such as overvoltage or voltage drop, in the distribution network with the penetration of PVs and MGs have been addressed in the literature through impact analyses to develop novel schemes of coordinating voltage control leading to increased PV generation [18, 19]. Table 1 presents a comprehensive review of such impact analyses [2022] or coordinate strategies [2328] of Open Distribution System Simulator (OpenDSS) modeling during solar PV penetration. For example, in Ref. [20], patterns from real measurements were used for a probabilistic load model, and the irradiance, temperature, and load values at each time step were processed in MATLAB and sent to OpenDSS to calculate the three-phase power flow. The model was tested with the comprehensive test feeder, and an error of 2.78% was obtained compared with the original results. Thus, the voltage variation issue due to the solar generation intermittency was the subject of [26], where OpenDSS and MATLAB were used to simulate the proposed feeder over 24 hours. Accordingly, a PI-based voltage control method is suggested to reduce the voltage deviations at the PCC of the PV inverters. In summary, OpenDSS modeling from MATLAB [20, 21, 23, 24, 2629] or CYMDIST [25] through COM server was used in the modified IEEE test feeder 7 [28], 13 [21, 23, 25, 27], 17 [29], 34 [26], 55 [24], and 69 [29] and in American Electric Power (AEP) 12.4 kV [22].

Moreover, to make various cases and a detailed impact analysis during high PV penetration, a chronological Monte Carlo [20] and quasistatic time series [21] were used to analyze the voltage profile that interferes with standard voltage regulation practices [30] and loss [22, 2426], maintaining the voltage limit (voltage profile) [2022, 2426, 28], maximization user profits, and cost minimization [23, 27, 29]. Some approaches maintain the network voltage at an acceptable level using PI control based on Q/P control [26], autonomous coordinate voltage control [24], and conservation voltage reduction algorithms [27] to realize these voltage regulation issues. Some methods maximize the customer benefit based on operation constraints using model predictive control [23] or dynamic programming [25] in OpenDSS modeling. Works in Ref. [31] consider an IEEE 13 nodes test feeder including solar PV generation and distribution static compensator so as to perform harmonic mitigation and power quality analysis. In order to validate the control of the distribution static compensator, several scenarios are considered to investigate disturbances of power quality and harmonics, like insolation of solar radiation, outage of PV plant from grid, and synchronization of PV plant to grid. The purpose of [32] is the amelioration of risks related to solar photovoltaic (PV) systems integration into a commercial electric power grid and the presentation of their mitigation strategies; thereby, two common risk modeling approaches are compared [33]. It is a review about the theoretical works and industrial applications of hybrid AC/DC distribution systems through analyzing their efficiency and energy loss; then, the study can be considered for future research regarding the mitigation of PV energy penetration. In addition, more methods, used to overcome instability issues caused by high penetration of photovoltaic (PV) systems in low-voltage distribution networks (LVDNs), are presented in Ref. [34]. The contributions and shortcomings of the different methods are compared, when dealing with high penetration of PV units in the LVDNs to handle the impact of the overvoltage. These methods are as follows energy storage systems (ESSs) by storing excess energy through batteries during the integration of high penetration of PVs and low load levels, load side management (LSM) by adjusting or controlling power consumption from peak to peak, reactive power support (RPS) by controlling the amount of reactive power by PV inverters through droop control strategy, active power curtailment (APC) in PV inverters to deal with a negative impact of voltage fluctuation that resulted from power intermittency, and grid reinforcement (GR) by changing the network structure to reduce the voltage rise challenge.

The main contributions of the paper, with respect to the existing works described above, are as follows:(i)The paper investigates the issue of high PV penetration effects with respect to the total installed power generation capacity of PV power plants in the MG and presents a coordinating voltage control solution to maintain the voltage limits of Vmin = 0.95 and Vmax = 1.05, which were chosen according to the ANSI standard.(ii)The current existing work with OpenDSS COM interface with MATLAB has only been small size of IEEE model less than 73 such as IEEE 7 [28], 13 [21, 23, 27], 17 [29], 55 [24], and 73 [20] because of the required complex interface skills. However, we have proposed the bigger MG model-based IEEE 123 model compared to the existing work.(iii)In the existing similar works of our topic with OpenDSS-MATLAB [21, 23, 24, 2729], there are three papers, which deal with voltage stability in terms of PU [20, 21, 28]. In addition, there are no cases as extreme cases. However, our paper has more extremely cases such as 45, 60%, 100%, and 200% compared to 10% [20] as well as tried to analyze on the impact of distributed PV in terms of various indicators as power losses, FVSI, and voltage stability compared with other 3 cases.(iv)A new coordinated voltage control, combining two control techniques using smart inverters and capacitor banks, is implemented with the integration of different indexes like mentioned above, leading to a new control algorithm approach different from the existing methods in the open literature as described in [22, 27, 31, 34].(v)The studied MG architecture considers real data from Los Angeles, California, such as the weather conditions and commercial building data. The MG is a modified version of the IEEE test feeder 123 and was simulated using OpenDSS. The results were exported to MATLAB using the COM interface [19].(vi)Different scenarios are considered for testing the robustness of the proposed coordinated voltage control.(vii)A comprehensive approach is proposed, highlighting the steps to achieve a hybrid simulation between OpenDSS and MATLAB and validate the proposed coordinate voltage control.

The rest of the paper is organized as follows. Section 2 presents the configuration of grid-tied MG operations, which include solar PV with different building loads for OpenDSS modeling. Section 3 describes the system framework of the OpenDSS tool COM interface with MATLAB. Section 4 provides the proposed coordinated voltage regulation control strategy of the PV smart inverter (SI) with capacitor banks (CBs). Finally, Section 5 evaluates the simulation results for the proposed system in various cases.

2. Microgrid Configuration

In this paper, the IEEE 123-node test feeder is employed to study the issue of voltage violation due to the high level of solar PV energy penetration. The detailed parameters of the IEEE 123 test model are described in Table 2. The technical details are given in a separated document and can be found at IEEE PES test feeder website [19]. The feeder is a composition of three MGs, separated by different switchers, SW1, …, SW5, as indicated in Figure 1. Several commercial buildings and solar PV power plants were added to the feeder at different locations. The objective is to study the feeder operation using real-time data for the load in different building types and solar irradiance. The data on the electric power consumption of the buildings were provided by the Department of Energy [30]. Thus, solar irradiance can be extracted using the PVGIS website [35].

In total, 10 PV power plants were added to the IEEE 123-node feeder, with a total of power generation capacity of 10 MVA, presented in Figure 1. In this case, the PV power plants supply and support the whole feeder, including the commercial buildings. The main objective is to study and analyze the effect of the PV energy penetration in the proposed MGs.

All data from the Los Angeles area are considered part of the case study for this paper (Figure 2).

The temperature inside the cell module is estimated based on the thermal model provided by [36]

The value of the back-surface module temperature Tm can be calculated using (2). The environmental conditions, namely, solar irradiance, wind speed, and ambient temperature measurements, are employed to calculate the back-surface temperature. These measurements are found in recognized databases, such as NASA and PVGIS [35, 37]:

Two coefficients, a and b, should be determined empirically, and the temperature difference between the cell and module back-surface must be calculated to determine the back-surface temperature. The values were chosen for polymer, thin-film, or steel PV technologies, according to Ref. [25]. Therefore, the PV power output is calculated using

Figure 3 displays the PV power generation under the influence of the weather conditions: temperature, wind speed, and irradiance.

We analyzed the annual load consumption for each type of building used in the MG model. Table 3 lists different types of loads ordered from highest to lowest for the annual load consumption for each building, which is from the load consumption of the interior equipment, fans, cooling, and interior lights. In most cases, the interior equipment and lights contribute to the largest load consumption.

Furthermore, Figure 4 illustrates the combined annual power-sharing between different load types for all buildings. Despite being the primary consumption source for only the HOS, a significant amount of energy is consumed by cooling appliances.

The different building load behaviors are assessed and validated through the OpenDSS, an electric power distribution system simulator supporting distributed energy resource grid integration and grid modernization. Figure 5 presents the hourly consumption profiles for real and reactive power for each building for two days in January 2015.

3. Implementation of OpenDSS COM Interface with MATLAB

In this paper, the simulation framework chosen for the MG model is OpenDSS, developed by the Electric Power Research Institute [38]. It is widely used in the industry to simulate large- and small-scale balanced/unbalanced networks, yielding comprehensive support documentation [39]. It is up-to-date software following the electric grid modeling and analysis needs, specifically those related to the vast amount of renewable energy generation and penetration, smart grid applications in the modern utility grids, and harmonic analyses.

This software was chosen for several reasons. First, it is free and open source, making the model easily transferable, modifiable, and updatable as needed. Furthermore, its COM server has extensive cosimulation capabilities, allowing a driver program to run the simulation from almost any programming language/environment. In this paper, the driver is used in MATLAB, and the OpenDSS executive is described in Figure 6.

In this work, the OpenDSS tool is hosted under MATLAB as the COM interface. Thus, several commands should be used to ensure communication between OpenDSS and MATLAB. Most of these commands are described in Table 4. All commands can be found in Visual Basic for Applications-Distribution System Simulator (VBA-DSS) interface properties and methods under Excel software.

The following program structure can be followed when implementing the COM interface between OpenDSS and MATLAB.

4. Mathematical Modeling (Problem Formulation)

This study considers a radial distribution feeder with n nodes, including the representative load and solar PV with SI for the IEEE 123 test feeder model. We used to Classical DistFlow Equations, which is to express the simplified equation in terms of Power Flow. A well-recognized power solver called DistFlow is used to solve the power flow problems [40].

4.1. Power Flow Analysis

In Figure 7, node 1 represents the feeding substation whose voltage is assumed to be constant, whereas three different load-type models (, ) with combination of Smart Inverter (SI) and Capacitor Bank (CB) are described as follows:(1)Only PV (, ): active and reactive power from solar PV(2)PV-SI ( ): active and reactive power from Smart Inverter of solar PV(3)PV-SI-CB ( ): active and reactive power from the SI of solar PV with CB

The power flow equations with solar PV are given as follows [40]:

The feeder section between nodes n and n–1 could be expressed as . A higher power loss always occurs in the MG when increasing the number of PV installations incorporated into each circuit model. From the model, power loss in section n could be expressed as follows:

4.2. Smart Inverter (SI) of Solar PV

One of the significant effects of the DG is the reverse power flow, while the distribution grid should be designed to move electricity in one direction from large remote generators on the transmission system to customers. Reverse power flow generally occurs when the generation of a distributed electric power plant exceeds the local load demand when >, causing power to flow in the opposite direction than normal, which could make this feeder more vulnerable to voltage instability. The inverter can mitigate the reverse power flow by compensating for the reactive power [41]. The inverter must satisfy the relation S, where, P, Q, and S are inverter’s real power output, var output, and rated capacity, respectively. Thus, the maximum var magnitude that can be injected or absorbed by the inverter of solar PV at any given time can be estimated as  =  . We used Q, which is related to the difference of apparent power and actual real power of the PV inverter as follows:

4.3. Capacitor Bank

Let be the reactive power supplied by CBn where n = 1, 2, 3, …, N, and N is the number of CBs. and denote the weights associated with improving the power factor at the bus and reducing the voltage variation to maintain ANSI (0.95/1.05), respectively. The CBs could be calculated as follows [42]:

4.4. Index of Voltage Stability Using ANSI and Fast Voltage Stability Index

For a couple of decades, voltage stability assessment has received increasing attention due to the complexity of power systems and integration of solar PV in the distribution system. In this work, the voltage stability index in each bus is evaluated using ANSI and the Fast Voltage Stability Index (FVSI). The network’s “PV hosting capacity” is defined as the maximum PV system size that does not exceed the network operating standards when interconnected to any valid bus. Thus, the voltage of all network buses must remain within the ANSI limits of 0.95 ≤  ≤ 1.05 in the steady state. We can estimate the maximum load for weak bus identification using the FVSI [43, 44]. Taking the symbols I as the sending bus and j as the receiving bus, FVSI can be defined as follows [43, 44]:where Z is the line impedance, X represents the line reactance, denotes the reactive power at the receiving end, and indicates the sending end voltage. The general two-bus representation is illustrated in Figure 8.

The value of FVSIij should always be less than 1 to maintain the voltage stability in an MG, and if its value is greater than 1, the MG experiences a voltage drop, which can be controlled by regulating either the distributed energy resources or reactive power compensator [45].

5. Proposed Approach

The steps described by the flowchart are followed to maintain the voltage stability of the proposed MG configuration. After performing the power flow analysis, voltages at all nodes were extracted, and their values were compared to those defined by the ANSI limits (0.95 and 1.05 per unit). The SI control was performed to adjust the voltages if the voltage was not within these limits. Then, the voltage values were verified again if they did not respect the ANSI requirements, and the CBs were used to adjust their values. Finally, the FVSI and power loss were assessed to evaluate the performance of the proposed controls as shown in the Figure 9).

The algorithm of the SI is executed according to the following steps.

The CB control (part b of the flowchart) is based on the power factor adjustment. This power factor should be greater than 0.9 to compensate for the reactive power from the inductive load, maintaining the voltage within ANSI limits. However, the flowchart in Figure 10 (part a) details the SI control.Step 1: Control the SI and extract the values of the active and reactive power of the PV power plant until j exceeds a defined control iteration jmax through the following substeps:(a)Calculate the PV plant output power () based on the actual weather conditions using(b) Check the inverter status: on or off.The objective is to identify the inverter status of the previous iteration (t − 1) and change the status in the actual iteration (t) based on (c)Calculate the inverter’s desired active output power. When the inverter status indicates on, the active output power of the SI is calculated in(d)Check the inverter reactive power limits for the cases of provided and absorbed reactive power.The inverter reactive power limits QL,ac(t) depend on whether the desired reactive power Qd,ac(t) is provided or absorbed (i.e., positive or negative) as illustrated by the following equation:The reactive power limits QLp(t) and QLn(t) are defined by Eqs. (15) and (16) as a function of the active power delivered to the grid Pac(t):(e)Check the inverter capability or the maximum inverter capacity corresponding to its KVA rating and extract the new values of active and reactive power. When the inverter capacity is not exceeded, the expressions for the active and reactive powers, Pac(t) and Qac(t), are given by the following equations:However, Pac(t) and Qac(t) must be reduced to meet the reactive power priority when the inverter capacity is exceeded. The reduced active power under the condition of reactive power priority is given by the following equation:If the desired reactive power is greater than or equal to the inverter capacity , then the new value of the reactive power becomes the capacity of the inverter itself, as given by the following equation:Step 2: Sample the control elements: voltages at all nodes are sampled by the inverter control to verify the need for control action. When necessary, the control action is included in the inverter control list.Step 3: Execute the inverter control action for t and update the reactive power of the PV plant. The desired reactive power is calculated according to the predefined smart inverter volt-var function, calculated for the present monitored voltage using a predefined curve (Figure 11).

6. Results and Discussion

This paper investigated the performance analysis of the proposed SI-CB control strategy during several PV penetrations in 11 buildings in the MG test model for two days from January 1 to 2, 2015. The MG test model was designed using the OpenDSS software platform with MATLAB. Our proposed SI-CB control strategy was verified via the following five cases.

6.1. Case 1

Compare three approaches per unit (V) for 100% PV penetration in 11 buildings in the MG test model (Figure 12).

6.2. Case 2

Compare three approaches per unit (V) for each node in the MG test model during 100% of PV penetration (Figure 13).

6.3. Case 3

Determine the voltage stability to maintain the range of the ANSI limits versus distance per unit (V) and line for the voltage profile for reference and PV with SI-CB control during 100% PV power generation (Figure 14).

6.4. Case 4

Determine the best place for PV installation with SI-CB control in terms of the FVSI for four different levels of PV power generation (Table 5).

6.5. Case 5

Evaluate the total power loss (kW) for comparison between 100% PV with SI control and SI-CB control for two days (Table 6).

Figure 12 illustrates a voltage per unit of all loads in the MG under 100% PV with three different approaches: (a) without control, (b) with SI control, and (c) with SI-CB control. In conjunction with Figure 12(a), Figure 12(b) illustrates that the voltage profiles per unit also maintain acceptable voltages at all points along the distribution feeder between 0.95 and 1.05, except the STM load by applying the proposed controls. Moreover, Figure 12(c) illustrates the improved voltage profiles compared to when the coordinated control of the proposed SI-CB has not been enacted for the STM load.

Most utilities in the US follow the ANSI voltage standards [46] to maintain acceptable service voltages. The service voltage is the point where the electrical systems of the supplier and user are interconnected, which is normally at a smart meter. According to the ANSI standard, the service voltages should be within ±5% of the nominal voltage level (240 V), which translates to 0.95 to 1.05 per unit voltage [47]. Figure 13 indicates that SI control and SI-CB control can improve voltage stability to maintain the ANSI limits for each node during 100% PV power generation in the MG. However, SI control exceeds the 1.05 per unit in the node from 170 to 110 and after 280 nodes during this period. The baseline voltage profile along the feeder with 100% PV with no control reveals that the voltage profile exceeds ANSI (1.05), as depicted in Figure 14. After applying the SI-CB control of PV, the voltage profile is maintained in the range of ANSI limits (0.95–1.05 per unit).

The line with an FVSI close to 1.00 indicates that it is approaching instability; hence, one of the buses to connect this line will experience a sudden voltage drop, leading to voltage collapse. The FVSI of all the lines must be lower than 1 to ensure the stability of the power system. Table 5 indicates the best FVSI for each installation place for each PV power generation. In Table 5, MRA is the best place to install PV with SI-CB control for all cases, whereas HOT is the worst case on January 2, 2019, at 2 : 00 p.m.

When the grid is connected to PV, used in SI control, it results in a power-loss increase. Figure 15 and Table 6 demonstrate that the total power loss in PV with SI control increased slightly by 30 kW compared to without control (100% PV only) at 7 kW. Furthermore, the proposed strategy (100% PV with SI-CB) has slightly (2 kW) less than only SI control (100% PV with SI).

7. Conclusion

In this paper, a coordinate strategy based on SIs and CB control was proposed to mitigate the high PV energy penetration in a modified IEEE 123-node test feeder. The main contribution is to provide a comprehensive approach enhancing the proposed MG configuration in terms of voltage stability and ANSI standards. The proposed approach was achieved using the parallel processing of MATLAB and OpenDSS. Thus, several actual commercial buildings were connected to different buses, and actual weather data were considered to conduct a real feeder operation. Different cases of MG operation were considered for testing the robustness of the proposed algorithm for voltage mitigation. When the grid is connected to PV, the proposed strategy (100% PV with SI-CB) has slightly (2 kW) less than only SI control (100% PV with SI). MRA is the best place to install PV as 0.05 with SI-CB control for all cases, whereas HOT is the worst case as 0.64 on January 2, 2019, at 2:00 p.m. in terms of FVSI. The results prove its effectiveness when dealing with a high penetration level of PV power generation, controlling the SI and CB. The voltage stability was assessed using several metrics, such as power loss and the FVSI, which was less than 1 for the four studied cases. Thus, in all feeder nodes, the pu voltage was maintained in the range of ANSI limits (0.95–1.05 per unit) despite the different sever operation conditions because of the high level of PV penetration. In addition, particular interest was given to the simulation part by presenting essential steps to perform such a study. The presented work can be extended to perform studies for other feeders. The presented work can be extended to study other solutions such as integration of high-power storage solutions or fuel cells. Adding to the power criteria, other physical aspects related to the management of their dynamics, temperatures, and their best economic use are to be investigated. In addition, it is interesting to provide a comprehensive comparison study with other voltage mitigation techniques for the use of the OpenDSS tool to develop benchmarks for microgrid tests.

Abbreviations

Tc:Cell temperature inside the module (°C)
Tm:Measured back-surface module temperature (°C)
Ir:Measured solar irradiance on module (W/m2)
Ta:Ambient air temperature (°C)
Ir0:Reference solar irradiance on module (1000 W/m2)
WS:Wind speed measured at a standard 10 m height (m/s)
ΔT:Temperature difference between the cell and module back surface at an irradiance level of 1000 W/m2
a:Empirically determined coefficient establishing the upper limit for the module temperature at low wind speeds and high solar irradiance
b:Empirically determined coefficient establishing the rate at which the module temperature drops as the wind speed increases
Pmpp:Max power point power at a selected temperature of 25°C and an irradiance of 1.0 kW/m2
f:Correction factor versus temperature variation: f is defined as when the panel temperature increases from 25°C to 75°C, the power drops by about 22%
PPV,n:Real power when only the PV power plant is used
QPV,n:Reactive power when only the PV power plant is used
PPS,n:Real power when a smart inverter of solar PV is installed
QPS,n:Reactive power when a smart inverter of solar PV is installed
PCB,n:Real power from the smart inverter of solar PV with CB
QCB,n:Reactive power from the smart inverter of solar PV with CB
Pn:Complex real power flow
Qn:Complex reactive power flow
Vn:The magnitude of the voltage at node n
Sr:Inverter’s kVA rating
%Sr1:Percentage of inverter’s kVA rating. When the inverter is off, the power from the PV array (DC power) must be greater than this value for the inverter to turn on
%Sr2:Percentage of inverter’s kVA rating. When the inverter is on, it turns off when the power from the array (DC power) drops below this value
Ecurve:Inverter efficiency curve, characterizing the variation of the inverter efficiency as a function of the power from the PV array, Pdc, per unit of inverter’s kVA rating
QLp:Reactive power positive limit when provided
QLn:Reactive power negative limit when absorbed
Qpmax:Maximum value of reactive power that the inverter can provide to the grid
Qamax:Maximum value of reactive power that the inverter can absorb from the grid
%Pmin:Percentage of the Pmpp. The inverter cannot provide/absorb reactive power when its active output power is less than this value
%Pmax:Percentage of the Pmpp. The inverter can provide/absorb reactive power up to its maximum allowed value, Qpmax or Qamax, respectively, when its active output power is greater than this value.

Data Availability

Data are available upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research project was supported in part by PAQ-Collabora (PAR&I-Tk), 2019–2022 program managed by the Tunisian Ministry of High Education and Scientific Research under Grant LSE-ENIT-LR11ES15. This research was also supported by the Basic Study and Interdisciplinary R&D Foundation Fund of University of Seoul (2021).