Control of a Variable-Impedance Fault Current Limiter to Assist Low-Voltage Ride-Through of Doubly Fed Induction Generators

A fault current limiter (FCL) may be applied to assist the low-voltage ride-through (LVRT) of a doubly fed induction generator (DFIG). FCLs with fixed impedance, lack the flexibility to adjust their impedance to adapt to different LVRT scenarios. The direct switch-in and -out of the fixed-impedance FCL yields transient electromagnetic oscillations in the DFIG, which need to be addressed. In this paper, a variable-impedance FCL is implemented at the stator side of the DFIG to assist its LVRT, and a novel methodology is proposed to control the impedance of the FCL, with which the stator current oscillation is effectively constrained and the smooth switch-out of the FCL is realized to avoid continued active power consumption of the FCL and to restore the DFIG to its pre-fault working condition. Analysis of the LVRT transient is carried out, which lays the foundation for the control methodology to determine the impedance of the FCL based on calculation of the optimization goal. The feasibility and effectiveness of the control to the variable-impedance FCL are verified by the numerical analysis results, which compare the LVRT simulation results with the application of the fixed-impedance and the variable-impedance FCLs.


Introduction
Due to their the merits, which include variable-speed constant-frequency operation capability and low investment requirement, DFIGs have become one of the most widely applied types of wind turbine generators (WTGs) [1][2][3]. Stability issues with the integration of WTGs have become a hot research topic [4] and are becoming increasingly complicated with the development of the offshore wind farms. A comprehensive review on the power transmission technology involved in offshore wind farms was carried out in [5], where the capabilities of WTGs to realize LVRT and support the system frequency through the active power control were introduced in detail. The LVRT capability is a major concern in the fault scenarios of wind power systems [6,7]. To assist the LVRT of the DFIG, both software control schemes and hardware protection schemes have been developed [8]. The former usually involves improvement to the converter control of the DFIG, e.g., adjusting the power/current control references according to the LVRT scenario [9], and eliminating the impact of stator flux oscillation yielded by the voltage dip on the converter control through feed-forward compensation [10,11]. The latter assists the LVRT with the application of hardware protections, e.g., the crowbar [12,13], FCL [14,15], dynamic voltage restorer (DVR) [16,17], STATCOM [18], and so on. Software control schemes do not require extra investment in protection equipment; however, their ability to deal with severe voltage dips switch it out. However, the selection of the coefficient in the ramp function has not been studied, which will affect the LVRT performance. To control the variable-impedance FCL for the desired LVRT effect, a control strategy needs to be developed on the basis of the LVRT analysis.
In this paper, the control of a variable-impedance FCL placed at the stator side of a DFIG is studied to improve LVRT performance. First, the control scheme to RSC of the DFIG considering the switch-in of the FCL is developed. Then, analysis on the LVRT transient is carried out, based on which an index that considers both the ability to suppress current oscillation and the active power consumption of the FCL is calculated. During the LVRT process, the impedance of the FCL is controlled to obtain the optimal LVRT index. Numerical analysis results show that, with this control scheme, the variable impedance FCL is able to suppress the electromagnetic oscillations during the LVRT. Also, the impedance is gradually reduced with constrained oscillations during the switch-out process. Simulation results show that the variable-impedance FCL has advantages over the fixed-impedance FCL in the following aspects: superior LVRT performance, flexibility in control target setting, and adaptivity to different voltage dips. Improvements over the ramp function to control the impedance of the variable-impedance FCL are also verified based on comparison of LVRT performances.
This paper is organized as follows. Section 2 introduces the implementation of the variable-impedance FCL to assist the LVRT of the DFIG, and how the RSC control is adjusted in coordination with the switch-in of the FCL. Section 3 describes the analysis on the LVRT transient of the DFIG with the application of the variable-impedance FCL and optimizes the control to the impedance of the FCL. The effectiveness of the proposed control scheme to improve the LVRT performance with the variable-impedance FCL is verified by the numerical analysis results in Section 4. Conclusions are given in Section 5.

Control Principal of Bridge-Type FCL for LVRT of DFIG
A single-machine infinite-bus (SMIB) model of DFIG is illustrated in Figure 1, where V and I are sinusoidal voltage and sinusoidal current, subscripts s, r, GSC, and PCC denote the stator, rotor, grid-side converter, and point of common coupling. To assist the LVRT of the DFIG, the FCL is placed on the stator side of the DFIG. A bridge-type FCL with variable impedance is applied to the LVRT of the DFIG. Its configuration is shown in Figure 2. The bridge-type FCL contains a bridge path and a shunt path. Four diodes D1, D2, D3, and D4 form the bridge configuration. A DC reactor with a small value L dc is connected in series with the IGBT switch, and the inherited resistance of the DC reactor is represented by R dc . The free-wheeling diode D5 is equipped to the dc reactor to guarantee secure operation. The resistor R FCL and inductor L FCL form a shunt path, which provides impedance to limit the fault current when the IGBT switches out to drive the flowing current to the shunt path. When a fault is detected, the FCL is switched in by shutting down the IGBT path. The impedance of the FCL is controlled by adjusting the duty ratio of the control signal applied to the IGBT switch. The equivalent impedance of the FCL is given by Equation (1), where d is the duty ratio of the switch signal, ω 0 is the synchronous rotating speed, and Z eq FCL is the equivalent impedance of the FCL. When voltage dips at the generator terminal are detected, the variable-impedance FCL is switched in to suppress the fault current oscillation. The resistance of the FCL helps to increase damping of the amplitude of current oscillation. Once the oscillations have been sufficiently damped, the continued switch-in of the FCL consumes large quantities of active power. This constrains the active power output capability of the DFIG during the LVRT, resulting in the degraded frequency stability of the wind power system.
After the clearance of the fault, the voltage recovery process also yields intense electromagnetic oscillations, which requires the FCL to be switched in to suppress the current oscillations. After damping of the oscillations is achieved, the FCL needs to be switched out to avoid high stator voltage compared to normal operation value. Abrupt switch-out of the FCL needs to be avoided as it results in abrupt stator voltage loss, leading to the undesired voltage dip ride-through transient. With the variable-impedance FCL, the FCL is switched out in such a way that the impedance of the FCL is gradually decreased to relieve the transient oscillations.

Rotor-Side Converter Control Coordinated with Switch-In of the FCL
When the FCL is switched in at the stator side, the configuration of the DFIG changes. The control scheme to the RSC for the active and reactive power control of the DFIG needs to be adjusted accordingly. The adjustment to the RSC control is derived below.
The flux and voltage equations of the DFIG are given by V r = R r I r + jsω 0 ψ r + pψ r (3) where ψ is the flux, L m is the mutual inductance, s is the slip, p is the differential operator that calculates the derivative with respect to time, e.g., pψ s = dψ s /dt, and ω 0 is the synchronous rotating speed. Equations (2)-(5) are the vector model of the DFIG, which is widely applied in the LVRT analysis, as can be referred in [10].
For the vector control with direct-axis (d-axis) voltage orientation, i.e. V sd = V s , V sq = 0, decoupled control of the active power and reactive power of the stator is realized by regulating the d-axis stator current and the quadrature-axis (q-axis) stator current respectively. With the constant stator voltage, steady-state stator flux is obtained. The corresponding rotor current reference is calculated based on Equation (4).
As for the RSC control, the rotor current is regulated by adjusting the rotor voltage. Combing Equations (2)-(5), the rotor voltage is written as Rotor current related terms Cross-coupling terms The PI control to regulate the rotor current by adjusting the rotor voltage is designed based on the above rotor voltage equation. As for the rotor current related terms, they are replaced by the PI control link. The input of the PI control is the difference between the rotor current and its reference. The feed-forward terms are added to compensate for the cross-coupling terms and the stator flux-related terms. The change rate of the stator flux may be obtained based on the stator voltage equation (Equation (2)). However, when the FCL placed at the stator side is switched-in for the LVRT, the stator voltage equation is rewritten as Combining Equations (6) and (7), the d-axis and q-axis PI rotor current control of the RSC, considering the switch-in of the FCL, is given by where k p and k i are the proportional and integral gains of the PI control, respectively, and superscript ref denotes the reference value. The difference between the rotor current and its reference goes through the PI control link to adjust the rotor voltage to obtain the desired rotor current output. The integral link outputs the integral of the rotor current difference as a part of the rotor voltage output.
The corresponding diagram of RSC control is shown in Figure 3. The rotor current control is realized in the d-q coordinate, as show in Figure 3 where the direct axis and the quadrature axis rotor current references are directly prescribed. To implement this control scheme, a, b, and c phase quantities are firstly transformed to an α-β coordinate through Clark transformation and then transformed to a d-q coordinate through Park transformation. Through the above-mentioned transformation process, the d-q axis rotor current reference may be transformed back to phase a, b, and c curves.

Analysis of LVRT
To analyze the impact of FCL impedance on the LVRT of DFIGs, analysis of the LVRT transient is carried out.
The dynamic the rotor current may be obtained by solving the second-order differential equation combining the rotor voltage equation (Equation (6)) and the rotor voltage control equation (Equations (8) and (9)), as given by Solving Equation (10), the analytical expression of rotor current is derived as Equation (11): where superscript ini denotes the initial value.
As mentioned before, the impact of stator flux oscillation on the rotor current control is eliminated by the feed-forward compensation. In this case, the rotor current can be regarded as an external condition to analyze the transient stator flux and current. The main factor to consider is the voltage dip depth. With the switch-in of the FCL, its impact on the response of the stator winding to the voltage dip needs to be considered. To obtain the analytical solution to the stator flux and current, the procedure is illustrated in Figure 4. Given the fast response of the inner-loop current control, the rotor current is considered to instantly reach the reference. Based on this assumption, the stator flux and current are derived based on the differential equation of the stator voltage (Equation (7)).
Assuming the rotor current is controlled to its reference, Equation (7) is rewritten as a first-order differential equation of the stator flux, as given by Equation (12): Solving Equation (12), the analytical expressions of the stator flux and the stator current are given by Equation (13): where τ is the damping coefficient of the stator flux, k is the voltage dip depth. Superscript ∞ denotes the steady-state value.
To illustrate dynamics of the stator flux during the LVRT, the vector diagram of the stator flux is shown in Figure 5.
During the LVRT, the stator flux transits from the initial value ψ ini s to the post-fault steady-state value, ψ ∞ s .The dynamic stator flux has two components: the post-fault steadystate component and the difference between ψ ini s and ψ ∞ s . The latter decays with the exponential damping e −τt and rotates in the complex plane in a clockwise direction at the electrical angular speed (dL FCL /L s + 1)ω 0 . As time goes on, this component gradually decreases as the stator flux approaches its post-fault steady-state value.
The analysis of the LVRT transient during voltage recovery may be derived similarly.

Control of the Variable-Impedance FCL for LVRT of DFIGs
When applying a variable-impedance FCL to the LVRT, two aspects need to be compromised in the impedance control. At the beginning period of the voltage drop, a larger impedance is preferred to reduce the oscillation of the stator current. When the current oscillation has been damped to a small aptitude, the FCL impedance needs to gradually decrease to reduce the active power consumption of the FCL. To achieve a compromise between the above two effects, the impedance of the FCL is controlled based on a weighted sum of two terms. The corresponding optimization goal is written as Equation (14), which determines the control target of the FCL impedance. For each time interval, the optimal FCL impedance is updated.
where W 1 and W 2 are the weight coefficients.
In Equation (14), the first term is the deviation of the stator current, which characterizes the effect of the damping to the current oscillation. The second term is the resistance of the FCL, which reflects the active power consumption. The weighted sum of the two terms works as a comprehensive criterion to control the FCL impedance. The calculation of the first term relies on the analytical calculation of the stator current, as previously derived in Equation (13). Considering that the impedance varies throughout the LVRT transient, once the optimal FCL impedance is determined for each time interval, the calculation results of the steady-state stator flux and the damping coefficient of the stator flux are updated to obtain the stator current dynamics.
As for the voltage recovery process, the impedance of the FCL is controlled in a similar way that inserts a large impedance to suppress transient oscillation yielded by the voltage rise, then gradually reduces the impedance to restore the DFIG to the pre-fault operation.
With the damping to the stator current oscillation, the second term in the optimization criterion given by Equation (14) ensures that the optimal impedance will eventually come to zero, i.e., the FCL will eventually be switched out. Thus, it is an appropriate optimization target for the control of the variable-impedance FCL. Still the weight coefficients need to be tuned to improve the LVRT performance. Setting W 1 to a high value results in slow switch-out of the FCL, while setting W 2 to a high value yields degraded performance in stator current oscillation suppression.

Numerical Analysis
A case study was conducted to analyze effect on the LVRT of the variable-impedance FCL. A test system with a DFIG, as shown in Figure 1, was adopted for the numerical analysis. Related parameters of the test system are listed in Table 1 [32]. Parameters of the C p function are given by [33]. The introduced LVRT scenario is as follows: the stator voltage drops to from 1 p.u. to 0.2 p.u. at t = 0.1 s and lasts for 0.625 s (the LVRT requirement of the Chinese grid code) [34]. At t = 0.725 s, the stator voltage recovers to 1 p.u.

Verification of LVRT Analysis
In this section, the accuracy of the analysis of the LVRT transient is validated. Under normal operation, the rotor current reference is 0.4891-j0.3239 p.u. During the LVRT, the rotor current reference is set to 1.05-j0.3 p.u. The resistance and reactance of the FCL are set to 0.1 p.u. and 0.15 p.u., respectively. As for the rotor current and stator current, the calculation results obtained from the analytical model are compared with the time-domain simulation results, as shown in Figures 6 and 7.  The analytical calculation and simulation results were obtained based on the same modelling of the DFIG. The simulation results were directly obtained by solving the differential-algebra equation set of the DFIG model using the numerical algorithm. Simulations are easily carried out, but they are incapable of providing direct quantification of the studied parameters (rotor and stator current in this case) and require great computational effort. The analytical calculation derives the expression of the studied parameter based on the same system model. Thus, the analytical calculations are capable of obtaining the same results as the simulations do.
As seen from the above figures, the analytical calculation results are consistent with the dynamic simulation results, which verifies the accuracy of the analytical LVRT analysis. It can also be seen that the rotor current quickly reaches the reference value, which verifies the assumption adopted to solve the stator current.
As seen in Figure 6, the rotor current increases during the transient, which results from the change in the rotor current reference. During the LVRT of the DFIG, the rotor current reference is adjusted to 1.05-j0.3 p.u. from its pre-fault value, 0.4891-j0.3239 p.u. The increase in the direct axis rotor current reference (from 0.4891 to 1.05) leads to the increase in the rotor current. The direct axis rotor current is related to the active power control of the DFIG, and it is increased to enhance the active power output capability, which helps to support the frequency stability of the wind power system during the LVRT.
Based on the analytical expression of the stator current (Equation (13)), the properties of the transient stator current oscillation shown in Figure 7 were analyzed. The decay of the stator current oscillation is related to the time constant, i.e., (R s + dR FCL )/L s , while the period of the stator current oscillation is obtained based on the electrical angular speed, i.e., (dL FCL /L s + 1)ω 0 .

Comparison between LVRT Effects with Variable-and Fixed-Impedance FCL
In this section, the LVRT effects of the variable-and fixed-impedance FCL are compared to show the improvement brought by the variable impedance. The impedance of the FCL is adjusted to strike a balance between the damping of the current oscillation and the active power consumption of the FCL. The maximum resistance and reactance of the FCL are set to 0.5 p.u. and 0.75 p.u., respectively. The resistance of the FCL, based on the optimization target (Equation (14)), is shown in Figure 8. The reactance of the FCL is in proportion to the resistance, thus has a similar profile. It is not given here to avoid redundancy.
From the above figure, it can be seen that the resistance of the FCL is set to a high value to suppress the transient oscillation of the stator current, and then it gradually decreases to reduce the consumption of the stator active power. A similar switch-out process of the FCL with its impedance gradually approaching zero is observed at the voltage recovery stage.
The LVRT effect of the variable-impedance FCL is compared to those with the fixedimpedance FCLs. Three impedance values, 0 + j0 p.u. (i.e., No FCL), 0.2 + j0.3 p.u., and 0.4 + j0.6 p.u., are set in the fixed-impedance FCL group. The comparison of the stator current, stator voltage, rotor voltage, active power output of the stator, and the active power consumption of the FCL are shown in Figures 9-13, respectively.      As seen from Figures 9 and 11, with the higher impedance of the FCL, the current oscillation is suppressed and damped at a faster rate. Considering that the change rate of the stator flux is introduced to the RSC control as the feed-forward compensation term, the rotor voltage rise is also suppressed with the higher impedance of the FCL. These are the advantages of adopting an FCL with high impedance to assist the LVRT.
In the meantime, the high impedance FCL reduces the active power output capability of the DFIG during the LVRT. As for the simulation scenario adopted in this section, the rotor current references are set to the same values, which means that the same active power output of the DFIG is delivered to the integrated grid. As seen from Figures 12 and 13, apparently the active power of the stator is increased with the high impedance FCL to support the stator voltage (seen in Figure 10). However, the increased stator active power does not contribute to the active power output delivered to the integrated grid but dissipates in the FCL. As shown in Figure 12, with the fixed impedance 0.4+j0.6 p.u., the active power of the stator already exceeds its pre-fault value. In this case, if the active current reference is further increased to provide more active power for the grid frequency support, the stator active power will rise to a higher value with the high impedance FCL, leading to serious deacceleration of the DFIG.
With control of the variable-impedance FCL, the advantages of the high impedance FCL scheme and the no FCL scheme are retained while their disadvantages are avoided. With the adjustable impedance value ranging from 0+j0 to 0.5+0.75 p.u., the large impedance is applied at the early stage of the LVRT to effectively constrain the stator current oscillation. Then the impedance is gradually reduced to switch out the FCL. Its active power consumption is reduced while the stator current oscillation is still constrained to a low amplitude.
As for the above simulations of the fixed-impedance FCLs, the transient process of their switch-out are not considered, which is another important aspect to consider for the operation of the FCL, as they eventually need to be switched out to let the DFIG return to its pre-fault condition during the voltage recovery. In Figure 10, the fixed-impedance FCL stays switched in during the whole LVRT process. Thus, during the voltage recovery stage (after 0.725 s), the grid voltage recovers to the reference while the generator terminal voltage (stator voltage) is further boosted due to the switched-in FCL, i.e., the stator voltage rises to a value higher than the pre-fault reference. This is likely to cause violations of the security constraint of the stator voltage. Thus, the fixed-impedance FCL needs to be switched out during the voltage recovery stage.
As for the voltage dip ride-through process, the fixed-impedance FCL needs to be switched out to restore the DFIG to normal operation condition, i.e., the stator voltage is restored to the pre-fault reference. Additionally, switching the FCL out early helps to avoid the impact of accumulated heat on the security properties of the FCL. The smooth switch-out is another advantage of the control to the variable-impedance FCL. With the switch-out process of the FCL considered, the LVRT effects are further compared, as shown in Figures 14-16. The fixed-impedance FCLs are switched in at t = 0.1 s and t = 0.725 s to deal with the voltage dip and recovery. They stay switched in for 0.1 s, then they are switched out.   As can be seen from the above figures, for both the voltage dip and voltage recovery processes, abrupt changes in the stator voltage occur immediately following the direct switch-out of the fixed-impedance FCL. Once again, oscillations of the stator current are the result, especially with the large impedance. For example, with a fixed-impedance of 0.4 + j0.6 p.u., the switch-out of the FCL yields severe voltage dips, leading to high amplitude stator current oscillation. The rotor voltage is increased to compensate for the stator flux oscillation. In contrast, with control of the variable-impedance FCL, its impedance is gradually decreased to switch out the FCL. The amplitude of the stator current oscillation is low during the voltage dip ride-through, which also helps to quickly restore the DFIG to its pre-fault working condition at the voltage recovery stage.

Analysis of the Impact of Weight Coefficient Settings
Control of the variable-impedance FCL is affected by the settings of the weight coefficients. With the weight coefficient W 1 set to 1, the LVRT effects of different values of the weight coefficient, W 2 , are compared in Figures 17-19.
As can be seen from the above figures, control of the variable-impedance FCL may be alternated by adjusting the weight coefficients in the optimization target. With a larger value of W 2 , the control of the variable-impedance FCL puts more emphasis on reducing the active power consumption of the FCL, in which case the impedance value is decreased at a faster rate (seen from Figure 17) to reduce the consumed active power (seen from Figure 19). This is realized at the cost of reduced damping of the stator current oscillation (seen from Figure 18). The flexibility of the control over the variable-impedance FCL lies in adjusting the weight coefficients, whereby the FCL may be controlled to deal with intense oscillations yielded by severe voltage dips (increasing W 1 ) or to quickly switch out without consuming large amounts of active power (increasing W 2 ).

Adaptivity to Different Levels of Voltage Dips
Compared with the fixed-impedance FCL, another advantage of the variable-impedance FCL is that its impedance may be flexibly adjusted to adapt to different levels of voltage dips. In this section, a different LVRT scenario is introduced with the voltage dropping to 0.6 p.u. The impedance of the fixed-impedance FCL is set to 0.5 + j0.75 p.u., which equals the maximum impedance value of the variable-impedance FCL. The LVRT effects with fixed-and variable-impedance FCLs, as well as the no FCL condition are compared in Figures 20 and 21.  As can be seen from Figure 20, when a minor voltage dip occurs, the switch-in of large impedance FCL (0.5 + j0.75 p.u.) yields a high stator voltage that exceeds its nominal value. When the stator voltage constraint is considered, the fixed-impedance FCL may not be allowed to switch in, which will result in the undamped stator current oscillation as shown by the simulation result of the no FCL scenario in Figure 21. With the variableimpedance used to assist the LVRT, the high impedance is only applied at the early stage of the LVRT to suppress the inrush fault current. The stator overvoltage is relieved by gradually reducing the impedance value. Still, during this process, the oscillation of the stator current is constrained to comparatively small amplitudes (seen from Figure 21). More security constraints like the stator voltage constraint may be added as constraints to control the impedance of the FCL for the improved LVRT performance.

Verification of Improvement to the Control Strategy for Variable-Impedance the FCL
In this section, the effects of the control strategies for the variable-impedance FCL are compared to show improvements due to the proposed control scheme. With the same LVRT scenario applied in Section 4.2, the impedance of the FCL is controlled by the ramp function as well as the proposed control scheme. With the ramp function, the impedance of the FCL is controlled with different slope coefficients. The resistance of the FCL with the ramp function control is illustrated in Figure 22. As shown in Figure 22, the impedance of the FCL is controlled by ramp functions with large and small slope coefficients. In the case of the large slope coefficient, the impedance of the FCL is decreased at a faster rate, i.e., the FCL is more quickly switched out after its switch-in to deal with abrupt stator voltage changes.
With different control strategies for the variable-impedance FCL, the stator current, stator voltage, and active power consumption of the FCL are compared in Figures 23-25, respectively.   From the above comparison results it can be seen that the effect of the variableimpedance FCL controlled by the ramp function to assist the LVRT is affected by the slope coefficient. With a large slope coefficient, the impedance of the FCL is quickly reduced, thus it fails to provide sufficient damping to the stator current oscillation (seen from Figure 23). With a small slope coefficient, the impedance of the FCL remains effective longer, which not only yields a long period of stator overvoltage (seen from Figure 24), but also consumes large amount of active power during the LVRT (seen from Figure 25). With the proposed strategy, the need for effective damping of current oscillation and timely switch-out of the FCL are compromised, thus the proposed control strategy achieves improved overall LVRT performance in comparison with the ramp function control.

Conclusions
In this paper, a novel control methodology for variable-impedance FCLs to assist the LVRT of DFIGs was proposed. First, the RSC control of the DFIG was coordinated with the FCL implemented on the stator side. Then the analysis of the LVRT transient of the DFIG was carried out to quantify the LVRT effects with the different impedances of the FCL. During the LVRT process, including the voltage dip and the voltage recovery stages, the impedance of the FCL is controlled based on optimization to a comprehensive index considering both the ability to suppress the stator current oscillation and the active power consumption of the FCL. Based on the numerical analysis results, the following conclusions were drawn: (i) The LVRT analysis of the DFIG shows satisfactory accuracy compared with timedomain simulation results.
(ii) With the proposed control scheme of the variable-impedance FCL, its impedance is gradually decreased to ensure timely switch-out of the FCL while constraining the stator current oscillation. The control scheme works like a soft switch for inserted FCL impedance, which helps to improve the LVRT performance compared with fixed-impedance FCLs as well as variable-impedance FCLs with the comparatively simple control (e.g., the ramp function control).
(iii) The control scheme of the variable-impedance FCL provides flexibility in control, e.g., the control target may be flexibly adjusted by changing the weight coefficients of the optimization index. Moreover, new optimization targets (e.g., reactive current injection capability) and constraints (e.g., stator voltage constraint) may be added as extensions to the existing scheme to realize the desired LVRT performance through control of the variable-impedance FCL.

Conflicts of Interest:
The authors declare no conflict of interest.