Research Article  Open Access
Yongheng Yang, Frede Blaabjerg, "LowVoltage RideThrough Capability of a SingleStage SinglePhase Photovoltaic System Connected to the LowVoltage Grid", International Journal of Photoenergy, vol. 2013, Article ID 257487, 9 pages, 2013. https://doi.org/10.1155/2013/257487
LowVoltage RideThrough Capability of a SingleStage SinglePhase Photovoltaic System Connected to the LowVoltage Grid
Abstract
The progressive growing of singlephase photovoltaic (PV) systems makes the Distribution System Operators (DSOs) update or revise the existing grid codes in order to guarantee the availability, quality, and reliability of the electrical system. It is expected that the future PV systems connected to the lowvoltage grid will be more active with functionalities of lowvoltage ridethrough (LVRT) and the grid support capability, which is not the case today. In this paper, the operation principle is demonstrated for a singlephase gridconnected PV system in a lowvoltage ridethrough operation in order to map future challenges. The system is verified by simulations and experiments. Test results show that the proposed power control method is effective and the singlephase PV inverters connected to lowvoltage networks are ready to provide grid support and ridethrough voltage fault capability with a satisfactory performance based on the grid requirements for threephase renewable energy systems.
1. Introduction
More and more singlephase photovoltaic (PV) systems are connected to the public grid mainly because of the matured PV technology and the declined price of the PV module cell [1]. As it is reported by PHOTON International, there was 27.7 GW of global PV generation systems installed in 2011, nearly 21 GW in Europe, which makes accumulated installed global PV capacity be around 70 GW at the end of 2011. It can be foreseen that the future gridconnected PV systems will play a very important role in the regulation of the conventional power systems.
However, behind the thriving penetration of singlephase gridconnected PV systems, the phenomenon of largescale Distributed Generation (DG) systems connected to the public grid also raises some concern on the availability, quality, and reliability of the grid. On this basis, many grid requirements are published in recent years by the Transmission System Operators, Distribution System Operators, or international standard committees (e.g., IEC and IEEE) in order to regulate the interaction between DGs and the public grid [2–9]. In spite of this fact, most of these grid requirements are focused on threephase wind power systems or DGs connected to medium or highlevel grid. There are few grid codes and papers describing the demands for PV systems connected to lowlevel grid, like singlephase PV systems [6, 10–15].
Due to the fast growing installations of gridconnected PV systems, it is expected that the future PV systems will be more active and “smart” with the functionalities of lowvoltage ridethrough capability and grid support, which are required for threephase wind power systems. In the future, the gridconnected PV system needs to fulfill the basic grid requirements, like power quality demand, and also provides some ancillary services likewhatthe conventional power plant does. For instance, the PV inverters should inject some reactive current when the grid presents a voltage sag fault. In that case, the control strategies should be developed in order to help the PV system ridethrough during the grid faults.
For singlephase applications, the calculation of active and reactive power is a big challenge, which leads to the difficulty of direct instantaneous power control. Unlike threephase systems, it is difficult to employ directly a dqrotating synchronous reference frame in singlephase systems to implement positive sequence control or constant reactive power control [7, 8, 16]. Possible solutions for singlephase applications are based on the singlephase PQ theory, instantaneous power control, or the droop control concept [13, 17–20]. Moreover, the grid fault detection and synchronization techniques are also challenges for the singlephase system when operating in the grid fault mode.
In this paper, the overall control strategy for singlephase PV systems in the grid faulty mode operation is demonstrated based on the grid requirements defined for threephase wind power systems as it is expected that they will be the basic requirements for singlephase PV systems. A transport delay grid fault detection method is adopted in this study. A power control method based on the singlephase PQ theory and the ProportionalResonant (PR) current control scheme is proposed. Finally, the system is tested in simulation and verified by experiments. The results show that the singlephase PV systems are ready for providing ancillary service in the future and the grid detection method has a significant impact on the performance of the singlephase PV systems under grid faults.
2. Grid Requirements
The grid requirements for DGs are updated and revised based on the development and penetration level of DGs and issued in the form of national grid code [2–6], such as E.ON grid code and CEI grid code published by Comitato Elettrotecnico Italiano. In these grid codes, the basic demands are imposed like the power quality (the THD level of the injected current) and the antiislanding requirement. Considering the penetration level of DGs connected to the grid, combined requirements are expected to be added into the grid codes in the future.
For example, in the German grid code, the DGs should handle ridethrough grid fault, known as lowvoltage ridethrough (LVRT) capability, and support the grid during the voltage recovery. Similar requirements can be found in other countries, like Denmark, Italy, Spain, and USA in which the DGs share a high percentage of power generation. The LVRT capability is defined so that the DG should stay connected to the grid within a short time and inject some reactive power to support the grid when the grid presents a voltage fault. Different LVRT curves with defined stayconnected time are shown in Figure 1.
As it is shown in Figure 1, the DGs should stay connected to the grid when the voltage drops to 0 V in the German grid code in 0.15 seconds. At the same time, the system should inject some reactive current in order to support the grid recovery according to Figure 2.
These grid requirements are set to ensure the safety of utility maintenance personnel and also other devices connectedtothe public grid and to avoid the grid collapse due to voltage faults. It is expected that these grid requirements will be recommended to be used to regulate singlephase PV systems, since the impact of largescale singlephase PV systems on the lowlevel public grid cannot be ignored. The discussion in this paper is based on these requirements.
3. SinglePhase GridConnected Photovoltaic System in the Grid Faulty Mode Operation
3.1. System Configuration
In consideration of the whole conversion efficiency, a singlestage PV system is connected to the grid through a transformer, which is shown in Figure 3. In some cases, a voltageboost stage can be used between the DC capacitor and the inverter, depending on the output voltage level of the PV panels. The boost stage makes the system flexible track the Maximum Power Point (MPP) of the PV panel and also makes it possible to handle the power variation across the DC capacitor at twice the grid fundamental frequency [9].
The control system shown in Figure 3 consists of a voltage sag detection unit, an Orthogonal Signal Generator (OSG) for power calculation, the power profiles including MPPT control in normal operation, active and reactive power controllers for current reference generation, a grid current controller, and a PWM generation block.
As it can also be seen in Figure 3, in order to implement grid fault operation mode for this system, a voltage sag generator is used by switching the series resistor and the load resistance . During the LVRT mode operation, the switch is open in order to limit the shortcircuit effect on the upstream grid and the switch is closed, while in normal operation, is closed and is open.
An filter () is used as the output filter in order to limit the highorder harmonics coming from the inverter switching behavior. Thus, neglecting all the filter losses, the system at the AC side can be described by in the time domain and in the Laplace domain, where is the injected grid current, is the inverter output voltage, and is the grid voltage. Hence, the equivalent model of the output filter can be given as it is shown in Figure 4.
3.2. Voltage Sag Detection
Grid condition information is very important for the control system to perform special functionalities. Similarly, as a part of the grid condition detection, the voltage sag detection is the way to identify the voltage fault and it determines the dynamic performance of the voltage sag compensator [21, 22] and the behavior of the whole control system. It is easy to understand that the system should switch from normal operation mode to grid fault operation mode as soon as possible, once the voltage sag is detected. Hence, it is essential to have a precise and fast fault detection unit.
Recently, many voltage sag detection methods have been proposed, such as Root Mean Square (RMS) method, peak value method (OSG based sag detection techniques), the missing voltage technique, and wavelet transform method [10, 21–25]. The RMS method computes the RMS value of the input signal by using a running average window of one cycle. Therefore, theoretically, this method will introduce one cycle delay in order to give the correct value. The wavelet transform based voltage sag detection is becoming of high interest [24], because of its fast detection speed. However, the implementation of this method is unfortunately not easy in real time.
In this paper, the peak value method is used to detect the voltage fault. The implementation of this method is easy, and by means of the OSG system, it is possible and flexible to calculate the average active and reactive power instantaneously. Since the T/4 Delay structure is used to build up the OSG system as it is shown in Figure 5, the system will introduce a quarter period time delay of the input fundamental period. The peak value of the input voltage can be given by where is the orthogonal signal of the input voltage generated by shifting a quarter period of the input signal.
3.3. Power Calculation and Power Profiles
In different operation modes, the power profiles are different. Firstly, the average active power and reactive power of the system are computed using Discrete Fourier Transform (DFT) structure, as it is shown in Figure 6.
Under normal grid conditions, the system is operating in Maximum Power Point Tracking (MPPT) mode in order to deliver as much energy as possible to the grid, and the power factor is normally near to one. According to the grid requirements, in grid faulty mode operation, the singlephase PV inverter should inject some reactive power into the grid to support the grid voltage during recovery. The amount of reactive power should be injected dependent on the voltage sag depth and the inverter rating current. At the same time, the PV panels should switch to nonMPPT operation mode in order to avoid tripping the inverter overcurrent protection [10, 26, 27]. The operation modes are described in Figure 7.
By employing Park transform () to the grid voltage and current in the OSG systems, the active power and reactive power delivered to the grid can also be calculated and expressed as where the subscripts “” and “” denote the “” and “” components of the grid voltage/current in a rotating synchronous reference frame.
As it is stated above, in normal MPPT operation mode, the system is operating at unity power factor. Then, neglecting the switching power losses, the active power and reactive power can be obtained: in which and are the amplitude of the grid voltage and current in MPPT operation mode, respectively, is the tracked maximum output power of the PV panels.
In faulty grid operation mode, according to Figure 2, the required reactive current under voltage sags can be given as in which , , and are the amplitude values of the grid instantaneous voltage, initial voltage before grid faults, and the nominal grid voltage, , are the nominal current and the initial reactive current before a grid failure, and p.u.. By substituting (6) into (4), the required active power and reactive power reference in the grid faulty mode operation can be given as where is the output power of the PV panel in the grid faulty mode operation.
Normally, the nominal current is limited by the inverter maximum rating current, ≤ , where is the inverter maximum rating current and the initial reactive current is zero, = 0 A.
According to (7) and Figure 7, the PV panels should “derate” its output power when the grid presents a voltage fault, as it can be illustrated in Figure 8. The derated power level ( or ) is dependent on the voltage sag depth and the nominal current. If the voltage drop is small, either by decreasing the PV voltage reference or increasing it can limit the PV output power. However, it should be noted that the PV voltage must be in the inverter input voltage range.
3.4. Power Controllers and Current Reference Generation
Similar to the power calculation in the rotating reference frame, the active power and reactive power for a singlephase application can be computed with the help of OSG systems: in which the subscripts “” and “” denote the “” and “” components of the grid voltage/current in a stationary reference frame.
In this way, the currents in a stationary reference frame can be expressed as Since the active power and reactive power are DC quantities in steady state, it is possible to use PI controllers for power regulation, and thus the current reference can be given by where “*” denotes the reference signal, and are the transfer functions of the PI controllers for active power and reactive power, respectively, and which can be given by
3.5. Current Controller with Harmonic Compensation
The existing current control methods, such as Proportional Resonant (PR) [28, 29], Deadbeat control (DB), Repetitive controller (RC) [30], and PI control which needs a Park transform () [29], can be adopted in different operation modes. In this case, the PR controller is used. For the purpose of the suppression of the injected current harmonics, resonant harmonic compensation method is employed, which is shown in Figure 9.
Hence, the current closedloop transfer function () can be expressed as in which is the transfer function of PR controller with harmonics compensation and is the processing and PWM delay [28]. and have the following expressions: where is the proportional gain, and are the resonant and harmonic compensator gains, is the harmonic order, is the fundamental frequency, and is the sampling period.
Hereafter, according to the above illustration, the whole control diagram can be structured, as it is shown in Figure 10. It can be seen from Figures 9 and 10 that the harmonics in the grid voltage will propagate to the grid current reference . However, the cascaded harmonic compensators will suppress the harmonics if the gains around the resonance frequencies are high enough. The Bode diagram of the PR controller with selective harmonics compensation is shown in Figure 11, which can be used to tune the parameters, , , and . It is shown in Figure 11 that the selection of a large leads to a wide frequency band around the resonance point [8].
(a)
(b)
It is also worth to point out that the voltage variation of the DC capacitor (DClink voltage) at twice the grid fundamental frequency will introduce even harmonics in the injected current. Moreover, in faulty grid operation mode, this distortion is much more severe. In this case, the harmonic compensation method by cascading several resonant controllers will cause heavy computation burden and become less practical, when the higher order harmonics need to be compensated [30].
4. Simulation and Experimental Results
Referring to Figures 3 and 10, the whole system is simulated in MATLAB in order to verify the performance of the proposed power control method. In this system, the PV panel is modeled by means of a 3D lookup table [31]. The nominal maximum power is 1030 W, the voltage and current at the maximum power point are 400 V and 2.56 A, respectively. The Incremental Conductance (INC) MPPT method is adopted to capture the maximum output power of the PV system. This INC method is modified according to Figures 7 and 8 in order to limit the output power when the grid voltage drops, presenting a grid fault. The other simulation parameters are listed in Table 1. The PI controller parameters for this simulation are 6.2, 1.5, 1, and 50. Simulation results are provided in Figures 12 and 13.

(a)
(b)
(c)
(a)
(b)
As it is shown in Figures 12 and 13, the singlestage singlephase PV system is subjected to a grid fault (about 0.35 p.u.) at s. This voltage sag lasts for about 320 ms. Since the detection unit needs a few milliseconds () to give the correct operation condition of the grid, instead of an immediate response, the MPPT and the control units continue to operate as in a normal condition and thus the grid current rises in this time interval. It happens similarly after the fault is cleared.
Once the grid fault is detected, the MPPT unit is disabled and the control system is switched to grid fault operation mode as it is shown in Figure 7. Soon after, the system starts to inject reactive current as required and limit the active power output to prevent the inverter from overcurrent protection. After the fault is cleared (voltage rises to 0.9 p.u.), the system goes back to its normal operation, and again, it is trying to track the maximum output power of the PV panels. However, as it is shown in Figures 12 and 13, it may take some time because of the MPPT perturbing process.
In order to validate the effectiveness of the proposed control method, the system is also tested experimentally. Instead of some real PV panels, two DC sources are used as the input power. The whole system consists mainly of(i)two Delta Elektronika SM 3005 DC power supply, the nominal DC voltage, 400 V,(ii)two DC capacitors connected in parallel, the equivalent capacitance, μF,(iii)a Danfoss VLT FC302 3phase inverter which is configured as a singlephase inverter,(iv)a sag generator as shown in Figure 3, = 19.2 Ω and = 20 Ω,(v)a 5 kVA 3phase transformer with leakage inductance 2 mH and resistance 0.01 Ω per phase,(vi)a DS 1103 dSPACE control system.The PI controller parameters for active power regulation are 1.2 and 50, while the parameters for reactive power control are 1 and 50. The other parameters of the singlephase system are shown in Table 1. A voltage drop is generated by switching the manual switches and , which are shown in Figure 14. The experimental results are given in Figure 15.
(a)
(b)
(c)
(d)
As it is shown in the experimental results, when the grid is subjected to a voltage sag (0.45 p.u.), the detection and control units of the singlephase system cannot react to this sudden change immediately, and thus the grid current is going to enter into a zone where the current reference is not appropriate until the detection unit outputs the sag signal. During the LVRT operation mode, the system is injecting some reactive power into the power grid and at the same time the active power is limited in consideration of the current limitation of the PV inverter, which is shown in Figure 15(b). When the voltage sag is cleared (the voltage amplitude goes to 90% of the nominal value), the system returns to its normal operation mode and is injecting current at unity power factor.
It is also shown in Figure 15 that both the start of the voltage sag and the end of that present a transient time delay ( and ) which is approximately 5 ms (T/4). Firstly, it is because of the fault detection time delay (T/4 Delay sag detection method) as illustrated in the previous paragraph. The power profile might be another reason since it is a step change when the voltage goes back to 90% of the nominal value. The response of the current controller may deteriorate the dynamics as well. Therefore, fast and accurate grid detection and reasonable power profile are necessary for a singlephase system in the grid fault mode operation.
However, the simulation and experimental results demonstrate the flexibility of a singlephase system to provide grid support functionality. The power control method in the test is effective to do so and it can do it fast.
5. Conclusion
Based on the existing grid requirements, this paper discussed the potential of a singlestage singlephase system operating in grid fault condition. The power control method proposed in this paper is effective when the system is switched to grid fault operation mode. It can be concluded that the future singlephase gridconnected PV systems are ready to be more active and more “smart” in the regulation of power grid. Like the conventional generation systems, the singlephase PV systems can also offer some ancillary service in the presence of an abnormal grid condition. However, the performance is dependent on the grid detection and control methods.
Acknowledgment
The authors would like to thank China Scholarship Council (CSC) for supporting this Ph.D. project in the Department of Energy Technology at Aalborg University, Denmark.
References
 EPIA, “EPIAPhotovoltaic Market Report, 2011,” http://www.epia.org/. View at: Google Scholar
 F. Iov, A. D. Hansen, P. E. Sørensen, and N. A. Cutululis, “Mapping of grid faults and grid codes,” Tech. Rep., Risø National Laboratory, Technical University of Denmark, 2007. View at: Google Scholar
 E. ON GmbH, “Grid Code—High and extra high voltage,” http://www.eonnetz.com/. View at: Google Scholar
 M. Braun, G. Arnold, and H. Laukamp, “Plugging into the Zeitgeist,” IEEE Power and Energy Magazine, vol. 7, no. 3, pp. 63–76, 2009. View at: Publisher Site  Google Scholar
 Comitato Elettrotecnico Italiano, “CEI 021: Reference technical rules for connecting users to the active and passive LV distribution companies of electricity,” http://www.ceiweb.it/. View at: Google Scholar
 H. Kobayashi, “Fault ride through requirements and measures of distributed PV systems in Japan,” in Proceedings of the IEEEPES General Meeting, pp. 1–6, July 2012. View at: Google Scholar
 R. Teodorescu, M. Liserre, and P. Rodriguez, Grid Converters For Photovoltaic and Wind Power Systems, WileyIEEE Press, 2011.
 F. Blaabjerg, R. Teodorescu, M. Liserre, and A. V. Timbus, “Overview of control and grid synchronization for distributed power generation systems,” IEEE Transactions on Industrial Electronics, vol. 53, no. 5, pp. 1398–1409, 2006. View at: Publisher Site  Google Scholar
 S. B. Kjaer, J. K. Pedersen, and F. Blaabjerg, “A review of singlephase gridconnected inverters for photovoltaic modules,” IEEE Transactions on Industry Applications, vol. 41, no. 5, pp. 1292–1306, 2005. View at: Publisher Site  Google Scholar
 Y. Yang, F. Blaabjerg, and Z. Zou, “Benchmarking of grid fault modes in singlephase gridconnected photovoltaic systems,” in Proceedings of the IEEE Energy Conversion Congress and Exposition (ECCE '12), pp. 4370–4377, September 2012. View at: Google Scholar
 K. Fujii, N. Kanao, T. Yamada, and Y. Okuma, “Fault ride through capability for solar inverters,” in Proceedings of the European Power Electronics (EPE '11), pp. 1–9, 2011. View at: Google Scholar
 M. Z. C. Wanik, A. Mohamed, A. F. A. Kadir, and I. Erlich, “Low voltage ride through capability of fuel cell generation system connected to low voltage network,” in Proceedings of the IEEE Conference on Clean Energy and Technology (CET '11), pp. 369–373, June 2011. View at: Google Scholar
 R. A. Mastromauro, M. Liserre, T. Kerekes, and A. Dell'Aquila, “A singlephase voltagecontrolled gridconnected photovoltaic system with power quality conditioner functionality,” IEEE Transactions on Industrial Electronics, vol. 56, no. 11, pp. 4436–4444, 2009. View at: Publisher Site  Google Scholar
 C. H. Benz, W. T. Franke, and F. W. Fuchs, “Low voltage ride through capability of a 5 kW gridtied solar inverter,” in Proceedings of the 14th International Power Electronics and Motion Control Conference (EPEPEMC '10), pp. T1213–T1220, September 2010. View at: Publisher Site  Google Scholar
 G. M. S. Azevedo, G. Vazquez, A. Luna, D. Aguilar, and A. Rolan, “Photovoltaic inverters with fault ridethrough capability,” in IEEE International Symposium on Industrial Electronics (IEEE ISIE '09), pp. 549–553, July 2009. View at: Publisher Site  Google Scholar
 P. Rodriguez, A. Luna, R. MunozAguilar, F. Corcoles, R. Teodorescu, and F. Blaabjerg, “Control of power converters in distributed generation applications under grid fault conditions,” in Proceedings of the European Conference on Cognitive Ergonomics (ECCE '11), pp. 2649–2656, September 2011. View at: Google Scholar
 M. Saitou and T. Shimizu, “Generalized theory of instantaneous active and reactive powers in singlephase circuits based on Hilbert transform,” in Proceedings of the IEEE 33rd Annual Power Electronics Specialists Conference (PESC '02), pp. 1419–1424, June 2002. View at: Google Scholar
 S. Dasgupta, S. K. Sahoo, and S. K. Panda, “Singlephase inverter control techniques for interfacing renewable energy sources with microgridPart I: parallelconnected inverter topology with active and reactive power flow control along with grid current shaping,” IEEE Transactions on Power Electronics, vol. 26, no. 3, pp. 717–731, 2011. View at: Publisher Site  Google Scholar
 S. A. Khajehoddin, M. K. Ghartemani, A. Bakhshai, and P. Jain, “A power control method with simple structure and fast dynamic response for singlephase gridconnected DG systems,” IEEE Transactions on Power Electronics, vol. 28, no. 1, pp. 221–233, 2013. View at: Google Scholar
 J. C. Vasquez, R. A. Mastromauro, J. M. Guerrero, and M. Liserre, “Voltage support provided by a droopcontrolled multifunctional inverter,” IEEE Transactions on Industrial Electronics, vol. 56, no. 11, pp. 4510–4519, 2009. View at: Publisher Site  Google Scholar
 D. Kai, K. W. E. Cheng, X. D. Xue et al., “A novel detection method for voltage sags,” in Proceedings of the International Conference on Power Electronics Systems and Applications (ICPESA '06), pp. 250–255, November 2006. View at: Google Scholar
 D.M. Lee, T. G. Habetler, R. G. Harley, T. L. Keister, and J. R. Rostron, “A voltage sag supporter utilizing a PWMswitched autotransformer,” IEEE Transactions on Power Electronics, vol. 22, no. 2, pp. 626–635, 2007. View at: Publisher Site  Google Scholar
 B. Bae, J. Lee, J. Jeong, and B. Han, “Lineinteractive singlephase dynamic voltage restorer with novel sag detection algorithm,” IEEE Transactions on Power Delivery, vol. 25, no. 4, pp. 2702–2709, 2010. View at: Publisher Site  Google Scholar
 C. Fitzer, M. Barnes, and P. Green, “Voltage sag detection technique for a dynamic voltage restorer,” IEEE Transactions on Industry Applications, vol. 40, no. 1, pp. 203–212, 2004. View at: Publisher Site  Google Scholar
 O. C. MonteroHernández and P. N. Enjeti, “A fast detection algorithm suitable for mitigation of numerous power quality disturbances,” IEEE Transactions on Industry Applications, vol. 41, no. 6, pp. 1684–1690, 2005. View at: Publisher Site  Google Scholar
 S.F. Chou, C.T. Lee, P.T. Cheng, and F. Blaabjerg, “A reactive current injection technique for renewable energy converters in low voltage ridethrough operations,” in Proceedings of the IEEEPES General Meeting, pp. 1–7, July 2011. View at: Google Scholar
 C.T. Lee, C.W. Hsu, and P.T. Cheng, “A lowvoltage ridethrough technique for gridconnected converters of distributed energy resources,” IEEE Transactions on Industry Applications, vol. 47, no. 4, pp. 1821–1832, 2011. View at: Publisher Site  Google Scholar
 M. Ciobotaru, R. Teodorescu, and F. Blaabjerg, “Control of singlestage singlephase PV inverter,” in Proceedings of the European Conference on Power Electronics and Applications (EPE ’05), pp. 1–10, September 2005. View at: Google Scholar
 D. Dong, T. Thacker, R. Burgos, F. Wang, and D. Boroyevich, “On zero steadystate error voltage control of singlephase PWM inverters with different load types,” IEEE Transactions on Power Electronics, vol. 26, no. 11, pp. 3285–3297, 2011. View at: Google Scholar
 S. Jiang, D. Cao, Y. Li, J. Liu, and F. Z. Peng, “LowTHD, fasttransient, and costeffective synchronousframe repetitive controller for threephase UPS inverters,” IEEE Transactions on Power Electronics, vol. 27, no. 6, pp. 2994–3005, 2012. View at: Google Scholar
 Plexim GmbH, “PLECS 3. 3. 1 User Manual,” 2012, http://www.plexim.com. View at: Google Scholar
Copyright
Copyright © 2013 Yongheng Yang and Frede Blaabjerg. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.