Research Article | Open Access

# Synergetic Control of Grid-Connected Photovoltaic Systems

**Academic Editor:**Md. Rabiul Islam

#### Abstract

It is important to improve the dynamic performance and the low-voltage ride-through (LVRT) capability of a grid-connected photovoltaic (PV) system. This paper presents synergetic control for the control of a grid-connected PV system. Modeling of a grid-connected PV system is described, and differential-algebra equations are obtained. Two control strategies are used in normal operation and during LVRT of a PV system. Practical synergetic controllers with two control strategies are synthesized. The mathematical expressions are derived for computing control variables. The design of the synergetic controllers does not require the linearization of the grid-connected PV system. A grid-connected PV system with synergetic controllers is simulated in Simulink surroundings. The control performance is studied in normal operation and during LVRT. Simulation results show that the synergetic controllers are robust and have good dynamic characteristics under different operation states.

#### 1. Introduction

The world is faced with serious problems of energy depletion and environmental pollution. The research and development of photovoltaic (PV) technologies have become a hot topic in the world [1]. Solar PV is now used around the world as an important technology for the conversion of solar energy because of its cleanliness and security. The solar PV capacity increased 25% over 2014 to a record 50 GW, lifting the global total to 227 GW. The solar PV industry is one of the fastest growing high-tech industries [2].

The control of a PV system is an important and difficult task. A grid-connected PV system mainly includes maximum power point tracking (MPPT) and the control of the DC-AC converter. Much research on MPPT has been conducted, and various MPPT algorithms have been proposed. The maximum power point (MPP) may be located by the perturbation and observation (P&O) algorithms [3], the incremental conductance (InC) algorithm [4], and the artificial neural network algorithm [5].

Reference [6] proposes a modified InC algorithm. The algorithm eliminates the division calculations involved in its structure and improves the variable step size, which only depends on the PV power change. Reference [7] improves P&O algorithm. This approach combines ant colony optimization with the traditional P&O method to yield faster and efficient convergence. This improved P&O algorithm can recognize global MPP under partially shaded conditions. A new MPPT algorithm is proposed in [8]. This scheme uses a gray wolf optimization technique to track the global peak of a PV array under partial shading conditions. It can solve the problems such as lower tracking efficiency, steady-state oscillations, and transients as encountered in P&O.

A two-stage three-phase grid-connected PV system in [9] contains a DC-DC boost converter and a DC-AC VSC converter. Pulse width modulator (PWM) signals fire the converters. The DC-AC VSC control system has an external voltage control loop and an internal current control loop. The external control loop regulates the DC link voltage, and the internal control loop regulates grid currents. The grid currents are transformed into *d*-axis and *q*-axis currents to accomplish vector control. Two control loops adopt the proportional-integral (PI) algorithm. Reference [10] presents a novel sliding-mode (SM) control for grid-connected PV systems. A systematic adaptive procedure to calculate the band of the hysteresis comparators is developed to improve the performance of the SM controller. A vector controller can keep the maximum power delivery of the PV system. Reference [11] uses a probabilistic wavelet fuzzy neural network (PWFNN) to structure the reactive power controller for a grid-connected PV system. The balance of the active power between the PV array and the DC-AC converter during grid faults is controlled by the DC link voltage. The controller can improve the operation of the grid-connected PV system during LVRT. Reference [12] suggests injecting the maximum rated current to maximize the inverter power capability during LVRT. The strategy combines a proper balance between positive- and negative-current sequences. High- and low-power production scenarios limit the inverter output current to the maximum rated value and avoid active power oscillations. Reference [13] presents a new adaptive PI controller using the continuous mixed p-norm (CMPN) algorithm for enhancing the LVRT capability of grid-connected PV systems. The adaptive PI controller is used to control the DC-AC converter. The gains of the PI controller are changed by the CMPN algorithm online without the need to fine-tune or optimize.

Reference [14] reviews the general synergetic control design procedure. Synergetic control is applied to a DC-DC boost converter, deriving a basic control law. An adaptive control strategy gives better trade-off between large-signal stability and load step response time. Reference [15] introduces a practical synergetic controller to regulate the buck converters that coordinate pulse current charging of batteries. Simulation and experiment results show that the synergetic controller is robust for such nonlinear dynamic systems and achieves better performance than the standard PI controller.

This paper presents a design for the synergetic controllers for an internal current control loop of a three-phase grid-connected PV system. The theory of synergetic control is first described. Modeling of a grid-connected PV system is introduced, and differential-algebra equations are obtained. Synergetic controllers are derived in detail for normal operation and during LVRT of a PV system. Two control strategies are used during normal and LVRT operations of PV systems. The design of a synergetic controller does not require linearization of the PV system. The parameters and of the two synergetic controllers are the same, and equals . The parameters of the synergetic controllers are easy to determine, and the proposed control schemes are easy to achieve. An example given in Matlab is adapted according to the synergetic controllers for assessing the performance of synergetic control of the grid-connected PV system. The synergetic control is chatter-free, and the simulation results demonstrate the effectiveness of the proposed control schemes.

#### 2. Synergetic Control Theory

Synergetic control is a state space control method based on modern mathematics. Synergetic control is applicable to the control of nonlinear, dynamic, and high-dimensional systems. It can be perfectly analyzed by mathematical expressions.

The nonlinear state equation of a controlled system is where is the state vector of the controlled system , is a continuous nonlinear function, and is the control vector of .

The macrovariables are defined for each input channel as a function of the state variables. The synergetic controller directs the system to move into the manifold from any initial motion point where is the macrovector of .

The dynamic evolution of the macrovariable towards the manifolds is defined as follows [14]: where defines the rate of convergence of the system.

Substituting from (2) into (3) yields

The control vector can be acquired by substituting (1) into (4). The system can be controlled to stay in the desired manifold.

#### 3. Modeling of the Grid-Connected PV System

The main circuit of the two-stage grid-connected PV system is described in Figure 1. A PV array is connected to a power grid via a DC boost converter and a three-phase voltage source converter (VSC) [9].

The differential equations (5), (6), (7), and (8) can be written according to Figure 1. where is the grid voltage, is the output voltage of the DC-AC VSC converter, and is the alternating current.

Equation (8) is based on a three-phase () reference frame. , , and reference frames are shown in Figure 2 [16].

Equation (8) is transformed into (10) from the reference frame to the rotating reference frame using the sinus-based Park transformation. where , and is the Park transformation matrix given in (11). Consider

Equation (10) may be written in (12) and (13). Consider where and in (12) and (13) are control variables that control the DC-AC VSC converter.

The MPPT algorithm computes the duty cycle to control the DC boost converter. The relationship of inputs and outputs of the DC boost converter is depicted by (16) and (17). Consider where is the duty cycle of the DC boost converter.

The instantaneous active and reactive powers are defined by (18) and (19). Consider

Selecting , (18) and (19) change into the two following equations:

The model above can be used to design the synergetic control of the grid-connected PV system.

#### 4. Synergetic Control of a Grid-Connected PV System in Normal Operation

Control strategies of the DC boost converter and the DC-AC VSC converter in the grid-connected PV system must be made. MPPT is implemented in the DC boost converter in normal operation. There are a number of MPPT algorithms; however, this paper does not analyze them.

The DC-AC VSC converter is controlled using the synergetic control presented in this paper. The control variables and in (12) and (13) are derived by synergetic control theory.

The reference value of may be obtained from the external voltage control loop. where that is, where is the gain of the proportional term, is the gain of the integral term, is the reference value of the DC voltage , and subscript denotes a reference value.

There are 2 control variables, and therefore, 2 macrovariables must be selected. The first macrovariable is

Substituting (24) into (4), (25) is obtained:

Substituting (12), (23), and (24) into (25), the first control variable is computed by (26) as follows:

The second macrovariable is selected where is a constant.

Substituting (27) into (4) and considering (13), the second control variable is as follows:

The control variables and are computed by (26) and (28), respectively, to guarantee system stability in normal operation. The synergetic control scheme of the grid-connected PV system is shown in Figure 3.

#### 5. Synergetic Control during LVRT

The PV system should stay connected and support the grid with reactive power during the voltage dip. Therefore, the reference value of is , depending on the voltage magnitude. Figure 4 depicts the principle of voltage support in the event of grid faults [17].

The abscissa in Figure 4 stands for , and the ordinate is . and , where is the present voltage during the fault, is the rated voltage, is voltage before the fault, is the reactive current, is the rated current, and is the reactive current before the fault. If a voltage dip is more than 10% of the rated voltage, the generator should provide a reactive current amounting to at least 2% of the rated current for each percent of the voltage dip within 20 ms after fault recognition [17]. If a voltage dip is more than 50% of the rated voltage, the generator must inject the grid with a reactive power of 100% of the rated current.

The limitation of *I*_{d} is pu, so the current will not be greater than the rated current.

If the maximum power of the PV array at MPP is less than the power , the maximum power can be injected into the grid with . A MPPT is used during LVRT to obtain the maximum power and economic benefits, and the control strategies are the same as those in normal operation.

If the maximum power of the PV array at MPP is more than the power , is set, and the power output of the PV array equals for the power balance. The following control strategies are used:

The first macrovariable is selected as

Substituting (29) into (4) and considering (12), the first control variable is computed as follows:

Equation (28) is also used for the second control variable during LVRT.

A MPPT is not used, and the duty cycle *D* is determined by a PI controller to regulate the DC link voltage.
where is the initial value of .

The duty cycle of the DC boost converter is determined by (31) and (32) during LVRT. The DC-AC VSC converter is controlled on the basis of (28) and (30).

The control strategies described above are shown in Figure 5.

Equations (26) and (30) are derived according to the synergetic control algorithm. is variable and is computed by (22) in normal operation. Consequently, (26) is obtained for the control of the PV system in normal operation. is a constant, and (30) is derived for the second control strategy during LVRT. Equations (26) and (30) are used under different operating conditions of the grid-connected PV system.

The DC link voltage should remain stable to maintain good operation of the PV system. This requires a power balance in the PV system. If the power injected into the grid by the DC-AC converter is less than the output power of the PV array, will increase. If the output power of the converter is the same as the output power of the PV array, will not change. The two control strategies can satisfy the power balance in the PV system in normal operation and during a LVRT.

The PV array works at the MPP because of the MPPT in normal operation. If the output power of the converter is less than the output power of the PV array and increases, the DC voltage regulator will generate greater on the basis of (22). The synergetic control will result in more output power of the converter and achieve a power balance in the PV system after the regulation process.

The synergetic control maintains the output power of the converter for the second control strategy during LVRT. If the output power of the PV array is greater than and increases, the DC voltage regulator will give a smaller *D* according to (31) and (32). The voltage of the PV array will increase according to (17), and its output power will decrease due to the power-voltage characteristics of the PV array. The power balance in the PV system will be achieved once more, and will return to its reference value.

The DC-AC converter may be damaged due to the large current that passes through it. The two control strategies can prevent the converter from overcurrent in normal operation and during LVRT.

#### 6. Case Studies

Matlab software provides an example titled *Detailed Model of a 100-kW Grid-Connected PV Array* [9]. The controllers of the example are adapted to assess the performance of synergetic control of the grid-connected PV system.

The example includes a PV array with an open-circuit voltage of 321 V. is 500 V, and the rated AC voltage is 260 V. A distribution transformer has a voltage ratio of 25 kV/260 V.

##### 6.1. Normal Operation Simulation

MPPT used in the case study is based on incremental conductance with an integral controller that can ensure that the system operates in MPP when the radiation intensity and temperature change rapidly.

The DC-AC converter is controlled using a synergetic control scheme. Equations (26) and (28) are rewritten according to the symbols in the Matlab example [9]:

Figure 6 shows the VSC main controller containing synergetic control.

The VDC regulator in Figure 6 is constructed on the basis of (22). The inputs of the PLL and measurements block are primary voltages and currents of the distribution transformer. This block tracks the frequency and phase of a sinusoidal three-phase signal and performs Park transformation from a three-phase (*abc*) reference frame to a *dq*0 reference frame. The synergetic control block in Figure 6 is a Matlab function block. It includes Matlab code. The code computes 2 control variables, and , with (33) and (34). The control variables are transformed into the reference values of three-phase voltages to generate PWM and control the VSC converters.

The parameters and are selected and applied in (22), (33), and (34). The settings are and for the simulation in normal operation. The standard test conditions are 1000 W/m^{2} irradiance and 25°C temperature. The change of sun irradiance is shown in Figure 7. Figure 8 describes the change of temperature for the simulation.

The MPPT regulator changes the duty cycle to regulate the PV voltage for tracking maximum power. At sec, the MPPT is enabled, and Figure 9 shows the PV array output power. Maximum power is 100.4 kW at the standard test conditions. VSC converters deliver the power to the grid by synergetic control.

The root mean square (RMS) of the primary voltage of the distribution transformer is constant. The RMS of the current is directly proportional to the power. Figure 10 shows the primary current of the distribution transformer.

Figures 9 and 10 indicate that the grid-connected PV system can track maximum power and deliver the power to the grid when the radiation and temperature change rapidly. The system operates stably.

The bus voltage may change in normal operation. Simulation of voltage fluctuation is performed at the standard test conditions. Figure 11 depicts a 25 kV bus voltage . The RMS of decreases by 10% at sec.

The current is shown in Figure 12. increases when the bus voltage decreases.

The maximum power does not change because sun irradiance and temperature are constant at the standard test conditions. Figure 13 demonstrates that the power injected into the grid is almost constant during the grid voltage fluctuation. This means that the synergetic control performs well in normal operation.

##### 6.2. LVRT Simulation

LVRT is the capability of electric generators to stay connected to the grid during short periods of voltage dip. LVRT is an important feature of the generator control system. There are several standards for LVRT requirements [17, 18]. Figure 14 shows the LVRT requirement in [19]. These requirements may be used for large solar power installations.

Three cases are used to test the performances of the synergetic control during a LVRT. The three-phase fault block in Simulink is connected in utility grid to simulate a three-phase short circuit with arc resistances.

The feature of the PV system is tested first when the bus voltage decreases by 30% due to a grid fault. The reference value of the reactive current is pu on the basis of Figure 4. The limitation of is , and the power is pu. The maximum power of the PV array at MPP is 0.48 pu when it is simulated at an irradiance of 500 W/m^{2} and a temperature of 25°C. The first control strategy is used because the maximum power of the PV array at MPP is less than . is set; the other settings are equivalent to those used in the normal operation simulation. Equations (22), (33), and (34) and MPPT are used. The RMS of the voltage decreases by 30% at s, and the simulation results are shown in Figures 15–17.

Figure 15 depicts the reactive power injected into the grid. The PV system provides the grid with a reactive power of 42 kVar during the LVRT.

The active power injected into the grid is shown in Figure 16. The injected active power is 47.2 kW when the reactive power injected into the grid is 42 kVar at a 30% decrease of the grid voltage. The maximum power of the PV array at MPP is delivered to the grid, and the solar energy is fully utilized.

The rated primary current is 3.3 A. Figure 17 shows that the primary current is less than its rated value. Therefore, the PV system may stay connected to the grid and provide the grid with 47.2 kW and 42 kVar at a 30% decrease of the grid voltage due to grid fault.

Then, simulations will test the performances of the synergetic control when the grid voltage decreases to 0.2 pu. The second control strategy is used because according to Figure 4 and . The duty cycle *D* is computed according to (31) and (32). Equation (30) is rewritten using the symbols in [9]:

Equations (31), (32), (34), and (35) are used to compute the control variables for controlling the DC boost converter and the VSC converters during LVRT. The parameters are and . The settings are and pu.

A three-phase short circuit occurs in the power distribution system at sec for the LVRT simulation. The bus voltage decreases to 0.2 pu during short circuit. and are shown in Figures 18 and 19, respectively. The 25 kV bus voltage in Figure 18 decreases to 0.2 pu after sec.

Small fluctuations of primary current arise after the voltage dip. The peak value of during LVRT is less than double of that in normal operation.

Currents and are depicted in Figure 20. and are regulated to 0 and −1 pu, respectively, after the dip. The PV system stays connected and supports the grid with reactive power.

The DC link voltage is shown in Figure 21. It fluctuates slightly and then remains at 500 V after the dip.

Another simulation is made for testing the performance of the synergetic control at a bus voltage of 0.05 pu. The second control strategy is used. A three-phase short circuit is applied in the bus at sec, and the bus voltage decreases to 0.05 pu to simulate the voltage of the fault arc. The current settings are and pu, to avoid a large current. The reference value of current is pu, to provide the grid with reactive power after sec. Figure 22 shows three-phase currents flowing in the choke. The maximum value of the current is less than 2.0 pu. Synergetic control can limit the current peak during LVRT.

The example in [9] has the current regulator with PI control and simulates only the normal operation of the PV system. The regulator from (31) and (32) is added to the example to simulate the LVRT of the PV system. and pu are fixed. Two parameters of the current regulator with PI control are changed from 0.3 and 20 (for normal operation) to 0.03 and 0.002, respectively, for LVRT. Therefore, parameter tuning of the current regulator with PI control is difficult. Chatter has been the main obstacle for sliding-mode control systems [10]. Probabilistic wavelet fuzzy neural networks (PWFNNs) include a membership layer, probabilistic layer, wavelet layer, and rule layer [11]. PWFNNs are very complicated, such that it is difficult to determine the PWFNN parameters. The gains of the PI controller in [13] must be changed online using the CMPN algorithm for improving the LVRT capability of grid-connected PV systems.

The parameters of the synergetic controller described in this paper are invariant during normal operations and during LVRT, and parameter equals . There is no need to change the parameters of the synergetic controller with a complicated algorithm. Therefore, parameter tuning of the synergetic controller is easy, and the controller is sufficiently robust for use in the PV system. The design of the synergetic controller is simple and does not require linearization of the PV system. In addition, the proposed control schemes are easy to realize. The synergetic control is chatter-free and displays good static and dynamic performance.

Asymmetric short circuit faults may occur in power systems. Unbalanced grid voltages are comprised of positive, negative, and zero sequence voltage components. The zero sequence is not considered here because of three-wire systems. Unbalanced grid voltage sags will cause performance deterioration of the converter. Negative sequences result in DC-link voltage ripples and harmonic power. The positive and negative sequence currents should be controlled simultaneously to improve control performance. Reference [20] used separate current controllers for positive and negative sequences. Synergetic control may be used in two current controllers to control positive and negative sequence currents separately. It is possible to use synergetic control for improving the performance of PV systems during unbalanced voltage sags.

#### 7. Conclusions

Solar PV is an important renewable energy technology and does not generate pollution. PV systems are developing rapidly, and most PV systems are grid-connected. Research on the control of grid-connected PV systems contributes to the improvement of the operation of the distribution network and the PV system.

Synergetic control can be used for the control of a grid-connected PV system. The design of a synergetic controller does not require the linearization of the PV system. The mathematical expressions for computing control variables can be derived according to the synergetic control algorithm and the mathematical model of a grid-connected PV system. Two control strategies are used in normal operation and during LVRT. The parameters and of the two synergetic controllers are the same. The parameters of the synergetic controllers are easy to determine, and results indicate that the synergetic controllers are robust. The grid-connected PV system can obtain a maximum power point and inject the power into the grid when the radiation and temperature change rapidly. The DC link voltage and AC currents are limited, and the DC-AC VSC converters are not damaged during LVRT. The PV system can run with pu continuously and supply reactive power to the grid when the grid voltage decreases to 0.05 pu.

Synergetic control has good dynamic characteristics in normal operation and during LVRT and is the alternative solution for grid-connected PV systems.

#### Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (51177074, 51407095) and the Jiangsu Province Natural Science Foundation (BK20151548).

#### References

- A. Rosa,
*Fundamentals of Renewable Energy Processes*, Academic Press, Boston, MA, USA, Third ed edition, 2012. View at: Publisher Site - REN21,
*Renewables 2016 Global Status Report*, REN21 Secretariat, Paris, France, 2016, http://www.ren21.net/wp-content/uploads/2016/06/GSR_2016_Full_Report_REN21.pdf. - N. Femia, G. Petrone, G. Spagnuolo, and M. Vitelli, “Optimization of perturb and observe maximum power point tracking method,”
*IEEE Transactions on Power Electronics*, vol. 20, no. 4, pp. 963–973, 2005. View at: Publisher Site | Google Scholar - K. H. Hussein, I. Muta, T. Hoshino, and M. Osakada, “Maximum photovoltaic power tracking: an algorithm for rapidly changing atmospheric conditions,”
*IEE Proceedings-Generation, Transmission and Distribution*, vol. 142, no. 1, pp. 59–64, 1995. View at: Publisher Site | Google Scholar - P. Q. Dzung, L. D. Khoa, H. H. Lee, L. M. Phuong, and N. T. D. Vu, “The new MPPT algorithm using ANN based PV,” in
*Proceedings of the International Forum on Strategic Technology*, pp. 402–407, Ulsan, South Korea, 13–15 October 2010. View at: Publisher Site | Google Scholar - N. Zakzouk, M. Elsaharty, A. Abdelsalam, A. Helal, and B. Williams, “Improved performance low-cost incremental conductance PV MPPT technique,”
*IET Renewable Power Generation*, vol. 10, no. 4, pp. 561–574, 2016. View at: Publisher Site | Google Scholar - K. Sundareswaran, V. Vigneshkumar, P. Sankar, S. Simon, P. Nayak, and S. Palani, “Development of an improved P&O algorithm assisted through a colony of foraging ants for MPPT in PV system,”
*IEEE Transactions on Industrial Informatics*, vol. 12, no. 1, pp. 187–200, 2016. View at: Publisher Site | Google Scholar - S. Mohanty, B. Subudhi, and P. Ray, “A new MPPT design using grey wolf optimization technique for photovoltaic system under partial shading conditions,”
*IEEE Transactions on Sustainable Energy*, vol. 7, no. 1, pp. 181–188, 2016. View at: Publisher Site | Google Scholar - MathWorks,
*Detailed Model of a 100-kW Grid-Connected PV Array*, MathWorks, Natick, MA 01760-2098, USA, 2015, http://www.mathworks.com/examples/simpower/mw/sps_product-power_PVarray_grid_det-detailed-model-of-a-100-kw-grid-connected-pv-array. - N. Kumar, T. Saha, and J. Dey, “Sliding-mode control of PWM dual inverter-based grid-connected PV system: modeling and performance analysis,”
*IEEE Journal of Emerging Selected Topics in Power Electronics*, vol. 4, no. 2, pp. 435–444, 2016. View at: Publisher Site | Google Scholar - F. Lin, K. Lu, and T. Ke, “Probabilistic wavelet fuzzy neural network based reactive power control for grid-connected three-phase PV system during grid faults,”
*Renewable Energy*, vol. 92, pp. 437–449, 2016. View at: Publisher Site | Google Scholar - J. Sosa, M. Castilla, J. Miret, J. Matas, and Y. Al-Turki, “Control strategy to maximize the power capability of PV three-phase inverters during voltage sags,”
*IEEE Transactions on Power Electronics*, vol. 31, no. 4, pp. 3314–3323, 2016. View at: Publisher Site | Google Scholar - H. M. Hasanien, “An adaptive control strategy for low voltage ride through capability enhancement of grid-connected photovoltaic power plants,”
*IEEE Transactions on Power Apparatus and Systems*, vol. 31, no. 4, pp. 3230–3237, 2016. View at: Publisher Site | Google Scholar - E. Santi, A. Monti, D. Li, K. Proddutur, and R. Dougal, “Synergetic control for dc-dc boost converter: implementation options,”
*IEEE Transactions on Industry Applications*, vol. 39, no. 6, pp. 1803–1813, 2003. View at: Publisher Site | Google Scholar - Z. Jiang and R. Dougal, “Synergetic control of power converters for pulse current charging of advanced batteries from a fuel cell power source,”
*IEEE Transactions on Power Electronics*, vol. 19, no. 4, pp. 1140–1150, 2004. View at: Publisher Site | Google Scholar - MathWorks,
*abc to dq0, dq0 to abc*, MathWorks, Natick, MA 01760-2098, USA, 2015, http://www.mathworks.com/help/physmod/sps/powersys/ref/abctodq0dq0toabc.html. *Grid Code High and Extra High Voltage*, E.ON Netz GmbH, Bayreuth, Germany, 2006.- BDEW,
*Technical Guideline: Generating Plants Connected to the Medium-Voltage Network*, BDEW, Berlin, Germany, 2008, http://www.bdew.de. - DEIF,
*LVRT Capability-Test Results*, http://www.deifwindpower.com/wind-turbine-solutions/control-systems/lvrt-test-results - H. Chong, R. Li, and B. Jim, “Unbalanced-grid-fault ride-through control for a wind turbine inverter,”
*IEEE Transactions on Industry Applications*, vol. 44, no. 3, pp. 845–856, 2008. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2017 Junjie Qian et al. 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.