1. Introduction
One of the many specialized measuring devices used in photovoltaic systems is a photovoltaic emulator (PVE). It mimics the current–voltage output characteristics of photovoltaic modules. The emulators are particularly useful for evaluating Maximum Power Point Tracking (MPPT) algorithms implemented in DC/AC inverter input circuits. The shape of PV panel characteristics varies depending on the current solar radiation intensity, ambient temperature, wind speed and other environmental variables. The MPPT algorithm in real time ensures maximum power consumption of the inverter [
1,
2]. The power unit in the Maximum Power Point (MPP) is the Watt peak (Wp), and this is the maximum power taken from the PV panel under Standard Test Condition (STC) [
3]. Measurements can be taken regardless of daytime and weather, there is no need to use large-surface expensive photovoltaic panels and energy-intensive light sources illuminating the panel, and complex mounting systems are not required. Solutions of this type are certainly much cheaper and more convenient, and measuring stations are therefore easier to use and take up less space. Such solutions are implemented especially in university laboratories. Laboratory equipment generally represents a compromise between functionality and price. The authors of this paper have designed and built a simple and functional PVE for the testing and laboratory analysis of a single-phase, on-grid inverter.
There are many PVE solutions described in the literature. A significant group represents relatively simple and cheap impulse systems based on boost/buck converters with microprocessor control and a DC source [
4,
5]. Simplicity, high efficiency and low costs are the main advantages of such solutions. However, not all of them can relatively accurately imitate PV characteristics. There are large mismatch errors of emulated curves and long delay time of the dynamic response, approximately 1 s. Better output parameters are achieved by emulators with fast processors and converter drivers [
6], where the relative error of output mismatch reaches a maximum of some 4%. However, these are more sophisticated systems. Some solutions use commercial solar panels powered in parallel with the DC power supply [
7,
8]. Such emulators mimic the PV characteristics quite well and have no limitations as with pulse-switching systems; this means that they are able to reach extreme points on the characteristics and react faster to load resistance changes. The main disadvantage is that the shape of the emulated curve cannot be changed; it is the same as that of the panel used. In order for the emulator to generate the shape of a different characteristic, it is necessary to connect the proper panel. This solution is not practical and additional panels generate additional costs.
More complex PVEs built using pulse systems operate in environments such as Matlab/Simulink, LabView and real-time platforms, e.g., dSPACE, FPGA, FPAA, Arduino [
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19]. Devices of this type can well emulate PV characteristics, including partial shading effects [
10,
11]. Real-time control of PVEs based on both digital and analogue platforms allows for better accuracy of PV characteristics but significantly increases the demand for system resources. This implies high implementation costs for the measuring station. However, it is difficult to compare such emulators taking into account important operating parameters such as the accuracy of the emulated curves, dynamic responses or energy efficiency. This is due to the fact that some authors of such examinations may not report some results. Solutions of this type can also be used to emulate characteristics of several kW but then they have to cooperate with high-power, DC-programmable power supplies [
15,
20,
21]. However, the disadvantage is that the power supplies introduce additional time delays of output reaction to input variables and are expensive laboratory devices.
Another group of emulators is adopted to operate with high power. They can emulate panels connected in series, parallel or serial–parallel. High energy efficiency enables operation in single- and three-phase high-power systems [
22,
23]. Other solutions reproduce I–V characteristic shapes using serial-connected semiconductor diodes. They can relatively well reproduce the shape of PV characteristics; however, there are significant energy losses caused by diode heating, even when emulating low powers [
24]. Other comparable solutions are simple and cheap, but the emulated characteristics are less accurate and contain distortions [
25]. The emulators of smaller powers are usually linear analogue systems. Integrated current–voltage stabilizers or bipolar transistors are used as control elements [
26,
27,
28,
29]. Similar emulators, due to their simple design and control, may have a better output response to changes in input signals but larger mismatch errors in reproduced curves. The disadvantage is also their low energy efficiency, which limits their application to systems with power ranging from several dozen to several hundred watts.
Most of the proposed solutions, except for the diode and linear ones, are switching (pulse) circuits. The control of this type of circuit most often consists of iterative calculations of a single-diode or two-diode PV cell-equivalent circuit model [
30,
31]. Therefore, considerable calculation resources of control systems are used in real time, which introduces significant delay times of the output response with respect to the input variables. There are also difficulties in reproducing the entire characteristic curve of the PV module, especially in extreme areas, i.e. for output voltages and currents close to their maximum values [
32,
33]. Such limitations are related to the PWM pulse control of boost/buck converters. Moreover, pulse devices are a source of electromagnetic radiation; hence, there may be problems with electromagnetic compatibility (EMC).
The aim of the research work presented in this paper was to design, construct and analyze the performance of a simple and low-cost photovoltaic emulator. It will be used for the laboratory testing of low-power PV inverters connected to the LV grid. More details about the emulator and emulation method are given in
Section 3.1,
Section 3.2 and
Section 4.
The authors have preliminarily estimated the approximate cost of the proposed PV emulator: electronic parts—USD 25; mechanical parts, e.g., heat sinks, fans, wires, switches, etc.—USD 25; case—USD 20. In total, the cost is around USD 70. The estimation does not include the cost of a DC power supply, the type of which can be chosen according to current preferences.
2. Single-Diode Model of the PV Cell and the Effect of Model Parameters on the I–V Curve Shape
A photovoltaic cell converts the energy of solar radiation into electricity. The conversion process occurs in the n-type and p-type semiconductors of which the cell is built. In the semiconductor, the flux of solar radiation generates current carriers that are electron–hole pairs. It is an internal photoelectric phenomenon. The p–n junction separates these carriers and a voltage is generated across the junction; therefore, current starts flowing in the external circuit.
A single photovoltaic cell provides a very low amount of current at a low output voltage. A PV panel is constructed using more cells connected in series or combinations in series–parallel. This solution provides the appropriate value of current and voltage depending on the purpose of the panel. Connecting a larger number of cells together does not fundamentally change the shape of the output characteristics, but it shifts the short-circuit current ISC and open-circuit voltage VOC toward larger values.
Various mathematical models of photovoltaic cells/modules have been proposed, including single-diode, double-diode and three-diode models [
34]. For our specific, PVE-oriented application, it is sufficient to employ the simple, single-diode model.
Figure 1 shows a single-diode equivalent circuit model of the photovoltaic cell. The intensity of the photovoltaic current
IPh generated in the photovoltaic generator
IPhgen depends on the intensity and the light spectrum of the flux
φ reaching the surface of the cell. Furthermore,
IPh depends also on the size of this surface, the construction and the methods for semiconductor doping.
Resistances
Rs and
Rsh symbolize parasitic phenomena occurring in the PV cell. The
Rs value is determined by the resistances of connections of metallic electrodes to the semiconductor, as well as the resistances of both metallic fingers and collecting bus-bars. The
Rs resistance is also increased by the semiconductor material resistance of the cathode and anode regions, and external contacts. The
Rsh value depends primarily on the parasitic junction currents, resulting from defects (imperfections) of the crystal structure, as well as the type of semiconductor materials and impurities, cell temperature and irradiation intensity. The dark cell, where
IPh = 0, is described by the basic p–n semiconductor junction theory equation, commonly known as the Shockley diode model equation [
35]:
where
ID is the current of the semiconductor diode,
Isat is the reverse junction saturation current,
q is the electron charge,
V is the conduction direction voltage on the p–n junction,
T is the temperature (K), and
k is the Boltzmann constant.
If the cell is illuminated and the photovoltaic current
IPh starts flowing in the external circuit, the above resistances must be taken into account in the model. The shapes of both the diode current
ID and I–V characteristics depend on the saturation current of the semiconductor junction
Isat, but they do not depend on the photovoltaic current
IPh. The mathematical description of the PV cell model can be presented as [
36,
37]:
where
IPh is the photovoltaic current generated in the cell,
Vcell is the photovoltaic voltage across the cell terminals,
Icell is the photovoltaic current flowing through the cell’s external circuit,
Rs is the series resistance,
Rsh is the shunt resistance, and
A represents the dimensionless diode ideality factor.
As demonstrated in Equation (2), the single-diode model is described using five generally known parameters:
IPh,
Isat,
A,
Rs and
Rsh. It is a simple model, but sufficient to describe the current–voltage characteristics of PV with satisfactory accuracy [
38,
39].
Figure 2 shows a typical I–V characteristic, where the most important points are indicated. Denotations are as follows:
Open-circuit voltage (VOC)—maximum output voltage for open output terminals,
Short-circuit current (ISC)—maximum output current for load resistance equal to 0 (RO = 0 Ω),
Maximum power point (MPP)—the point of the characteristic where the power output on the emulator’s load resistance is maximum and equal to IMPVMP,
Voltage at MPP (VMP),
Current at MPP (IMP).
The most important parameter determining the quality of a PV cell/module is the fill factor (
FF), given as:
where
IMP is the output current value at the MPP, and
VMP is the output voltage value at the MPP.
The FF is a parameter informing about a cell/module’s quality, with higher FF indicating better cell or module quality.
Photovoltaic current IPh is created by converting light radiation into electrical energy. The current is directly proportional to the light stream intensity φ reaching the cell/module surface. The φ increase raises the ISC value and slightly increases the VOC value and vice versa. A serious problem appears if the temperature rises above the STC reference value. An increase in the cell’s operating temperature significantly reduces both the VOC value and the photovoltaic conversion efficiency, while slightly increasing the ISC value. The influence of the temperature on the output voltage and current is determined by the temperature coefficients of the open-circuit voltage and short-circuit current, kV and kI, respectively. They are given as a percentage change in VOC and ISC values per one degree Celsius or Kelvin.
Figure 3a,b shows the influence of parasitic resistances on particular fragments of I–V curves and consequently on its general shape. The slope and shape of the curve between the
VOC and MPP points depend on the resistance
Rs. An increase in
Rs results in a higher
VRs voltage drop, which decreases the output voltage. Consequently, this causes a slope increase for this plot part and reduces the power value at the MPP, as shown in
Figure 3a. For the part determined by the
ISC and MPP points, the resistance
Rsh has a dominant influence on the slope magnitude, as shown in
Figure 3b. The small
Rsh values represent high parasitic leakage currents
Ish, which in effect reduces the flowed photovoltaic cell current
Icell. In this case, unlike for
Rs, higher resistance is more advantageous because it reduces the slope of this curve part and increases the value of maximum power at the MPP.
The coefficient of ideality factor
A corresponds, to some extent, to the curvature (shape) of the I–V characteristics, especially within the knee area, as shown in
Figure 4. According to many researchers, there is no accurate method to determine the exact value of this coefficient. In the model calculation procedures, the initial value is arbitrarily chosen and then adjusted to obtain a relatively good fit for the model curve. Its values are usually taken as 1 <
A < 1.5 [
38].
The actual values of the specific points in I–V curves are published by the manufacturer’s data sheets. In addition, the following data are also provided: temperature coefficients of short-circuit current, open-circuit voltage, power and conversion efficiency. All published values are obtained under STC. Some manufacturers also provide data from measurements under Normal Operation Cell Temperature (NOCT).
3. Proposed PVE Device
3.1. Design Assumptions
During the PVE design, several important assumptions were taken into account. The design assumes a simple construction, low cost and general availability of electronic parts, simple algorithms and control systems and the possibility of easy, simple-to-perform modernization and adaptation to emulating higher-power PV panels.
Ranges of current and voltage operations have been designed. It was established that the set values of the output current would correspond to the
ISC current of the emulated PV panel. However, the set value of the unloaded PVE output voltage will correspond to the
VOC voltage at the terminals of the unloaded panel. These are the extreme points on the I–V characteristic, as shown in
Figure 2. Moreover, the span of the current and voltage ranges of the emulated characteristics was adopted. The output current, corresponding to
ISC, is stepwise adjustable at three ranges—1 A to 2 A, 2 A to 3 A and 3 A to 4 A—and then continuously adjusted within the selected range. The output voltage
VOUT, corresponding to
VOC, is only continuously adjustable within a single range from 28 V to 45 V.
The main assumption of the emulation method is the ability to independently adjust a slope and curvature in specific sections of characteristics, which results in the desired shape of the emulated characteristics, wherein the maximum relative mismatch is assumed to be less than 1%. The values of voltages and currents are set both in steps and continuously, and the minimum and maximum values have been adapted to the operating conditions of the possessed inverter. It was also assumed that the emulated output curves would be continuous, and free from distortions, oscillations, spikes and rapid jump changes, which could have a negative impact on the proper performance of the MPPT algorithms. Relatively high energy efficiency was assumed by limiting dissipated power, especially in the main output control element. Moreover, inexpensive, specialized control and measurement equipment will be used for the proper functioning of the PVE.
3.2. Description of the PVE Circuit
The presented PVE is a modified current/voltage stabilizer. It has a continuously adjustably output voltage and short-circuit current value and co-operates connected to an external DC power supply. It is capable of emulating I–V curves with VOC from 28 V to 45 V and ISC from 1 A to 4 A. Such values correspond to PV panels with powers from around 25 Wp to around 160 Wp. The adjustment and stabilization of the output voltage or current for a given operation mode is realized by using the control element CE, over current protection stages OCP1 and OCP2, potentiometer voltage divider PVD, error amplifier EA and reference source RS. From the control systems point of view, the blocks CE, OCP1, OCP2, PVD, EA and RS constitute the controller.
It is worth emphasizing that this controller, and so the PVE as a whole, is a very complicated, nonlinear, time-varying dynamical system, which is very difficult to tackle in a rigorous mathematical manner. Therefore, our PVE is designed in a graphical way, with support from the simple model in Equation (2).
The CE block contains an n–p–n bipolar power transistor operating as an emitter follower. The output current
IOUT is controlled by the over current protection stage OCP1. The block diagram is shown in
Figure 5 and general output characteristic is shown in
Figure 6.
A low-pass filter (LPF) is installed at the input to reduce higher-frequency interference from external sources that may occur at the DC power terminals. The PVE operates in the automatic regulator mode. The value of the output voltage
VOUT is set using a potentiometer voltage divider (PVD). A portion of this voltage, equal to the
nVOUT value, is fed to one of error amplifier inputs (EA), which is also a voltage comparator. The second input is supplied with the reference voltage
Vref as a voltage drop at either the Zener diode
Dz or resistance
RF, both operating in the RS block. The
Vref and n
VOUT values are compared and their error
Verr controls the power transistor in the CE block. The resistance
Rs is responsible for the slope magnitude between the points
PMax = (
IPmax, VPmax) and
P0V = (0,
Vmax) within the operation mode VM. Under the operation mode VM, current differences Δ
IOUT induce output voltage changes ranging from
VPmax to
VMax. If the output current
IOUT reaches the
IPmax value set in the OCP1, the output voltage achieves the
VPmax value. After exceeding this value, under a further increase in
IOUT, the system shifts to the current operation mode CM, and for output resistance
RO = 0, it reaches the point
PI0 = (
ISC, 0), as shown in
Figure 6. The resistance
RC influences the slope within the CM area, which means the slope between the points
PMax and
PI0. The resistance
R in the OCP1 block sets the value of the output current
IOUT in the assumed range Δ
ISC = 3 A. For ranges 1 A, 2 A, 3 A and 4 A, the values of
R are equal to 1.5 Ω, 0.75 Ω, 0.5 Ω and 0.375 Ω, respectively. In order to better match the emulated curves to the reference ones, an additional OCP2 circuit has been designed, which is coupled together with OCP1. The resistances
RB1 and
RB2 are used to adjust the curvature between the
PMax and
P0V. This solution allows a better adjustment of the characteristic shape for the VM mode and the knee area.
Depending on the operation mode, the PVE works as either a current or voltage stabilizer. It can be assumed that the operating point moves within certain limits around the
PMax. When the operating point moves below the
PMax, the voltage mode VM changes to the current mode CM. In contrast, if the operating point moves above the
PMax, the current mode CM changes to the voltage mode VM, as shown in
Figure 6.
The above output characteristics of the stabilizers are presented in the Cartesian coordinate system as the output voltage versus output current VOUT = f(IOUT). In the following sections of this paper, the curves emulated by PVE will be depicted in the same way as for the PV cells and modules, that means., the output current vs. the output voltage IOUT = f1(VOUT) and for power POUT = f2(VOUT).
4. Emulation Method
This is an original graphical (rather than analytical) method for shaping the emulated output characteristics I–V, consisting of the separate regulation of slopes and curvatures for the VM and CM regions. In this way, it is possible to obtain the shapes of the real PV module curves, which are schematically presented in
Figure 7. To determine the basic shapes of PV characteristics, manufacturer’s data for
VOC,
ISC,
VMP and
IMP are used. Now, the characteristic points
PMax,
PI0 and
P0V of
Figure 6 correspond to the points
MPP = (
VMP,
IMP),
P0I = (0,
ISC) and
PV0 = (
VOC, 0), respectively, in
Figure 7. To precisely match the curvature shapes of particular sections to the I–V curve of the selected PV module, both voltage and current values for a few specific points are required. These values are obtained using the mathematical model [
40] for the selected module.
Particular sections of the emulated characteristics correspond to the operating modes VM and CM of the PVE. The two modes can transfer between each other, so there is a transition region TR, in the vicinity of the MPP. The slopes of the individual sections with their curvatures are proportional to the ΔV and ΔI values, being the maximal voltage and current differences VOC-VMP and ISC-IMP, respectively. These differences are calculated into the relevant values of the resistances RS and RC, on which the slope depends. The curvatures of the sections are adjusted based on auxiliary characteristics’ current and voltage values IA, IA’, VA, VA’. Optionally, more values can be adopted and this depends on the emulator’s adjustment capabilities. Assumed, that the coordinates of these values are fixed and estimated (from the model (2)).
The distances between the MPP and auxiliary values and their positions on the characteristic I–V were selected empirically based on the measurement tests. For this purpose, several modules with different shapes of output curves were selected. The best results were obtained for the following fixed values: VA = 0.2 VMP and VA’ = 0.4 VMP, and IA = 0.75 IMP and IA’ = 0.5 IMP. For such values in each point, the estimated values of the current and voltage are obtained from the mathematical model. Increases in the current differences from IA’-IA and IA-IMP indicate increases in the slope and curvature of the current section CM. Higher voltage differences from VA’-VA and VA-VMP occur for higher slopes in the section VM. The values of calculated differences are converted into appropriate values of resistances RB1 and RB2 in the OCP2 block.
Our original emulated slope shaping method is currently being intensively developed taking into account the influence of variable lighting conditions on the shapes of the emulated curves. In addition, we plan to extend the method with the temperature coefficients kI and kV for the short-circuit current and open-circuit voltage, respectively.
7. Discussion
Based on the simulation and measurement results presented in
Section 6.1, it can be concluded that the ranges of the emulated characteristics Δ
ISC = 3 A and Δ
VOC = 17 V are in accordance with the adopted assumptions. The shapes of the PVE output characteristics obtained from the measurements are very similar to the simulation results, which can also be determined from the obtained values of
PMAX,
IMPP and
VMPP. The PVE correctly emulates the characteristics even far beyond the MPPT area, approaching closely the
ISC point and reaching the
VOC point, which is additionally confirmed by the small values of the output resistance.
Power loss in the PVE circuit is minimized by selecting the appropriate value of the supply voltage. This procedure limits the power loss in the control transistor to around 17 W with the range of emulated power from around 100 W to 155 W. It effectively increases the total energy efficiency of the emulator. The obtained results, presented in
Section 6.2, are satisfactory as for this type of analog device.
An important step in the measurements was to generate the I–V curves of the commercial PV module and compare them to those obtained by means of its mathematical model (2), as presented in
Section 6.3. The values of the absolute differences of the current characteristics and the resulting powers were compared. Very small values were obtained for currents and powers below 30 mA and 820 mW, respectively. Additionally, relative mismatch errors were calculated for both I–V and P–V curves. The maximum error value of ±1% adopted in the assumptions was not exceeded; see
Figure 14. Moreover, the obtained error results were lower than those of 2.5% to 3% reported in the literature [
14]. This proves that the emulator has a very good ability to imitate the desired shapes of the photovoltaic characteristics.
Dynamic response tests, presented in
Section 6.4, showed that the PVE can react quickly to time-varying load conditions of the output and does not generate any interference, harmonic frequencies, current–voltage peaks, etc. The obtained response times—see
Figure 15 and
Figure 16—are shorter compared to similar measurements reported in the analyzed literature. The tests were performed for changes in rising and falling resistances at the PVE output. The oscillograms of
Figure 15 and
Figure 16 recorded minimum and maximum values as well as differences in output voltages and currents.
The emulator is not a pulsed device, so there are no electromagnetic compatibility (EMC) problems. The emulated characteristics, in the entire reproduction range, do not contain ripples, pulsations or harmonic components. Such distortions could interfere with the MPPT algorithms by, e.g., increasing the number of iterations and overloading of the inverter processor.
The disadvantages include the manual input of the desired working conditions and characteristic shapes. This is executed by means of adjustment elements that have been calibrated during testing according to the indications of the measuring equipment. This limits, to a certain extent, the functionality of the proposed PVE; however, from the point of view of the emulator’s intended purpose, which is didactic tasks and research works at the laboratory, this is not so important at present.
The further work plan includes the adaptation of the PVE to emulate the characteristics of photovoltaic panels for output powers of up to 1 kWp. Development works are also being carried out concerning our original method for shaping slopes of emulated photovoltaic curves. Moreover, the automation of settings is under improvement using a combination of analogue and digital circuits and a dedicated application installed on a PC.