Research Article  Open Access
ChangJian Ru, RuiXuan Wei, "Safe Control for Spiral Recovery of Unmanned Aerial Vehicle", Abstract and Applied Analysis, vol. 2014, Article ID 983624, 10 pages, 2014. https://doi.org/10.1155/2014/983624
Safe Control for Spiral Recovery of Unmanned Aerial Vehicle
Abstract
With unmanned aerial vehicles (UAVs) widely used in both military and civilian fields, many events affecting their safe flying have emerged. That UAV’s entering into the spiral is such a typical safety issue. To solve this safety problem, a novel recovery control approach is proposed. First, the factors of spiral are analyzed. Then, based on control scheduling of state variables and nonlinear dynamic inversion control laws, the spiral recovery controller is designed to accomplish guidance and control of spiral recovery. Finally, the simulation results have illustrated that the proposed control method can ensure the UAV autonomous recovery from spiral effectively.
1. Introduction
Unmanned aerial vehicles (UAVs), with many unique advantages, have been widely operated in both military and civilian fields [1–5]. They can carry on all kinds of complex mission, such as aerial mapping [6], disaster monitoring [7], search and rescue [8], and reconnaissance and attacking [9]. With the UAVs widely used in the entire airspace and the increasing flight frequency, they have suffered from safety issues such as the crash, out of control, and collision with other planes [10, 11]. Particularly, two typical safety events happened on UAVs in recent years. One is that, in June 2012, a U.S. military UAV named as Global Hawk was stepping into the spiral during the routine training, eventually leading to spin and crash. The other is that, in 2011, a U.S. military UAV was intercepted and captured successfully by Iran, using the electronic jamming. So, how to effectively solve these safety problems and ensure UAVs fly safely and reliably has become an important problem to be addressed. The work in this paper is carried out just under such a background, but we mainly focus on one safety issue, that is, UAV’s spiral.
Currently, the relevant researches on the spiral of aircraft are mainly in the field of manned aircraft. In foreign countries, since the 1950s of last century, the USA has carried on the research of spiral recovery for manned aircraft. The representative outcome was released by Langley Research Center, where a kind of idea for spiral recovery was proposed and discussed in theory. In the subsequent air shows, American pilots successfully used this method to recover from the stall and the spiral. In China, aimed at the flight characteristics of manned fighters, a series of methods for spiral judgment and recovery were developed and formed into flight specifications and institutions for different types of aircraft. For example, for J6 fighter, the approach named as “flat, middle and push” is presented to the spiral recovery. In addition, other scholars have discussed the spiral technology for the aerobatics [12].
But at the same time, there is little research on the spiral recovery for the UAVs. For example, as is reported by a Russian website in 2008, a kind of professional technology is developed to avoid the UAV entering into the spiral. In China, individual scholar focuses on the problem of UAV’s spiral recovery under the condition of rudder locked and designs a control law which seems to be effective for spiral recovery [13]. But this study is only suitable for the moment of spiral mode, and it does not consider the unstable spiral caused by the large attack angle. From the above analysis, it can be found that there are very few relevant literatures and researches on the UAV’s spiral recovery, but for UAV the problem of the spiral is certainly a very important problem affecting UAV’s flying safely.
As is known, the UAV is a very complex nonlinear system. Particularly, when UAV is flying in the state of spiral, the nonlinear characteristics will be particularly evident. In order to improve the performance of nonlinear systems, several researches have recently proposed some typical stability and vibration control methods, such as faulttolerant control [14], slidingmode control [15, 16], feedback control [17–19], sampleddata control [20], optimization control [21, 22], and energytopeak control [23, 24]. The effectiveness of these approaches has been verified, respectively, by different simulations or experiments, but they cannot be directly introduced to solve the control problem of UAV’s spiral recovery. According to an expert from the US Air Force Research Laboratory, named Bruce T. Clough, autonomous safe control would be the most effective measures to improve the safety of the UAV [25]. Based on this idea, we try to propose the autonomous control method for the spiral recovery, ensuring the UAV recovery from the spiral effectively and making it fly safely.
In this paper, we contribute to the design of the controller to make the UAV recover from spiral autonomously. Firstly, we analyze the reasons why UAV may enter into the state of spiral. Secondly, the structure of the controller for spiral recovery is given and the control sequence of state variables is discussed. Then, based on the nonlinear dynamic inversion [26], the control law for spiral recovery is designed.
The rest of this paper is organized as follows. Section 2 is the cause formulation and analysis of UAV’s spiral. In Section 3, we focus on the design process of the controller for spiral recovery. Simulation experiments are provided to illustrate the effectiveness of the approach presented in this paper in Section 4. Finally, the conclusion and the further work are depicted in Section 5.
2. The Cause Formulation and Analysis of UAV’s Spiral
Spiral is an abnormal motion for an UAV. It is very likely to happen when an UAV is in the stalling state along with rotating around three axes. At that time, the UAV will descend sharply along the spiral trajectory with a small radius. Based on some relevant literature, the cause of UAV’s spiral can be divided into three categories. They are wing selfrotation, direction divergence, and rolling divergence, respectively. Next, we will analyze these causes.
(1) Wing SelfRotation. The coefficient of normal force, which is perpendicular to the wings of UAV, can be written as For the lowspeed UAV, it has a larger lift coefficient, but a smaller drag coefficient. So, the lift coefficient can be used to analyze the change of rolling damping coefficient at a large angle of attack.
When an UAV is flying at a large angle of attack, if the rolling motion happens, for the sinking wings, the angle of attack will increase, while the lift coefficient will decrease. On the contrary, for the rising wings, the angle of attack will decrease, while lift coefficient will increase. At that time, damping in roll will change the symbol, and the damping effect will help the rolling, leading to selfrotation of the wings.
(2) Direction Divergence. When an UAV is in the stall condition, the nose will deflect automatically. That is because the UAV has lost the static stability. When the UAV emerges the direction divergence, different lifts values between two wings caused by the sideslip may make the UAV both deflect and roll. The pitching inertia moment generated during above process may make the angle of attack further increase.
(3) Rolling Divergence. Rolling divergence usually occurs when the UAV is flying at a negative angle of attack. At this time, if the UAV’s lateral static stability , it will have the lateral static stability. Otherwise, if the absolute value of negative angle of attack exceeds that of the stall, that is, , UAV will lose the lateral static stability. At this time, it will be very likely to cause the rolling divergence.
3. Controller Design for Spiral Recovery
In this section, we will design the controller design for spiral recovery. This controller should be switched automatically when an UAV enters into the spiral. Meanwhile, it can carry on the recovery control of spiral autonomously according to the actual dynamic of spiral. The design of controller should take into account the following two parts: control sequence and control law.
3.1. Control Sequence of State Variables
When an UAV is flying in the stall, there will be coupling interaction and location constraints existing among each control surface, and then the effect of control surfaces may reduce [27, 28]. Consequently, it is difficult to control each state variable of UAV system simultaneously during the process of spiral recovery. According to the practical experience of spiral recovery used for manned fighters, the effective approach is to prevent the UAV from rotation firstly. If the angle of attack is controlled with a large angular rate of rotation, it will be insufficient for manipulation torque generated by pushing the stick forward to overcome the upward inertia torque. At that time, the angle of attack cannot reduce down. Even the rotation will become more quickly due to the effect of rolling and partial torque generated by pushing the stick forward. Meanwhile, the control of the attitude angle mainly depends on the ailerons. However, the efficiency of manipulating ailerons is lower at a high angle of attack. Therefore, the prerequisite for controlling the attitude angles is UAV at a small angle of attack.
In summary, one must grasp the control sequence of state variables accurately during the process of spiral recovery and give the outputs corresponding to control these relevant variables. By this way, it can prevent the spiral from becoming more divergent and make the most effect out of controlling surfaces. So, the control sequence of spiral recovery can be concluded in Table 1.

3.2. Control Law for Spiral Recovery Based on Nonlinear Dynamic Inversion
If UAV is in the spiral, the nonlinear characteristics will become particularly evident. Therefore, based on linear smalldisturbance motion model, control laws of three channels cannot be independently designed. It is well known that, for flight vehicle design, the most widely used method is nonlinear dynamic inversion control [29]. The fundamental principle of this method is that nonlinear inversion and nonlinear function are used to eliminate the nonlinearity of the controlled objective, thus forming global linearization. In order to illustrate this principle better, take the following affine nonlinear system as an example:
If is invertible, then let , where pseudocontrol input can be selected as . Because it is impossible to obtain the accurate invertible model of the system, , are used to express approximate compensation and inversion of this system. The structure of dynamic inversion control system is shown in Figure 1.
The dynamic inversion control approach can automatically adapt to the large change of flight conditions and configuration at the fixed gain. So, it will be more effective to control the UAV in the state of maneuvering flight. In this section, we will combine singular perturbation theory with nonlinear dynamic inversion to design control law of UAV’s spiral recovery.
As is analyzed above, the UAV’s nonlinearity is very evident when it is in the spiral. Consequently, the sixdegreeoffreedom motion model [30] should be introduced here, which is shown as follows:
According to (3), it can be clearly seen that this motion consists of 12 state variables. Among these variables, , , and are the rate of roll angle, pitch rate, and yaw rate, respectively; , , and are the angle of attack, sideslip, and roll, respectively; , , and are velocity, the angle of yaw, and trajectory inclination, respectively; , , and are the projection of UAV’s centroid on the ground coordinate system. Define a vector which consists of these above state variables as follows:
Based on the movement rule of UAV and singular perturbation theory, the state variables can be divided into three parts, which are fastspeed variables , the middlespeed variables , and the slowspeed variables or . By this way, nonlinear equation (3) can be divided into three different timescale subsystems. Based on the state feedback control, the structures of these three firstorder linear decoupling control structures can be obtained. At this time, these three different speed circuits can be designed independently. For the design of UAV’s spiral recovery, we are mainly concerned about the fastspeed and midspeed variables. The structure of nonlinear dynamic inversion control law designed is shown in Figure 2. Next, we will illustrate the design of these two loops.
(1) The Control Law of FastSpeed Loop. The fastspeed loop is used to control these variables , , . According to the given inputs of angular, , , , the deviation of control surface can be solved. The equation of fastspeed loop can be written as the following form: where , is a vector that consisted of eight state variables, and is the output of inner loop, and it is also the input control of UAV objective. The first item in the right of (5) represents nonlinear coupling torque, and the second item of that represents manipulating torque generated by control surfaces. These two items can be, respectively, shown as follows:
Let the expected angular acceleration rate of fastspeed loop meet the ideal closedloop dynamic response, which is as follows: where represents the expected angular acceleration and represents the command signal of angular acceleration generated by the lowspeed loop. The variables , , and are PID feedback gain of fastspeed loop, and their values are 20~30 rad/s. According to the dynamic inversion principle, the control input can be obtained as
Equations (7) and (8) are combined to form an integrated fastspeed loop, achieving the control law design of fastspeed loop for UAV. The control structure corresponding to the fastspeed loop is shown in Figure 3.
(2) The Control Law of MiddleSpeed Loop. The middlespeed loop is used to control these variables , , and . According to nonlinear dynamic inversion theory, the inputs , , can be solved. If deflection of control surfaces is ignored, this loop can be written as the following affine form: where and each element of is shown as follows:
Let the closedloop dynamic response of expected inputs of angular rate , , of fastspeed loop be with the following form: where , , are the bandwidth of middlespeed loop, which is selected in 2~4 rad/s.
Based on the idea of dynamic inversion, the output of control command can be obtained as follows:
Equations (11) and (12) are combined to form the integrated middlespeed loop, and at last control law design of UAV’s middlespeed loop is achieved. The control structure corresponding to the middlespeed loop is shown in Figure 4.
4. Simulation Analysis and Verification
To verify the effectiveness of the proposed method, three different comparative simulation scenarios are carried on. All the computations and experiments are operated on a computer with Core i3 CPU, 3.30 GHz, and Windows XP operating systems. These experiments are depicted as follows. Section 4.1 depicts how to make the UAV enter into the spiral. From Section 4.2 to Section 4.4, the effects of spiral recovery under three different conditions are compared.
4.1. Entering into the Spiral
First, we should try to make an UAU enter into the spiral. Here, the UAV used in the simulation is the fixedwing with the angle critical of attack 30°. At beginning, the UAV is supposed that it is flying in level, with height 2000 m and velocity 42 m/s. we cut off the control module of spiral recovery, and assume UAV is controlled by the autopilot. Carry on the operation of BHU, and the UAV will come into the spiral after a short time. The change of attitude angle and angular rate are shown in Figures 5(a), 5(b), and 5(c). As is shown in Figures 5(a), 5(b), and 5(c), the attack angle is changing periodically and is greater than the critical attack angle. Threedimensional trajectory of UAV flight is shown in Figure 5(d). According to Figure 5, it can be seen that UAV is entering in the spiral after the maneuvering operation.
(a) The response curve of attack angle
(b) The response curve of roll angle and yaw angle
(c) The response curve of attitude angular rate
(d) Flight trajectory of UAV’s entering into the spiral
4.2. Autonomous Control of Spiral Recovery
In this section, the work is to verify whether the proposed method can make the UAV recovery from spiral autonomously. Assume that UAV’s initial state and the operation of its coming into spiral are the same as those in Section 4.1. Once the UAV enters into the state of spiral, the control module of spiral recovery can be switched automatically. Under the action of spiral recovery controller, the change curves of UAV’s attitude angle and angular rate are shown in Figures 6(a)–6(c).
(a) The response curve of attack angle
(b) The response curve of roll angle and yaw angle
(c) The response curve of attitude angular rate
(d) Flight trajectory of UAV’s spiral recovery
From these figures, it can be seen that when the UAV has entered into the spiral, the control module of spiral recovery will be switched automatically. Then, given , , and , control the angular velocity of rotation. At s, UAV will cognize that the rotation rate is slower, and then, given , , and , control the attack angle and attitude angle. At s, UAV will cognize the attack angle is lower, and then, given , , . At s, UAV will stop rolling and achieve the desired attitude angle, ultimately recovering from the spiral effectively. The threedimensional flight trajectory of spiral recovery is shown in Figure 6(d). According to this figure, it can be clearly seen that the UAVs have succeeded to recover from the spiral. That means autonomous control of spiral recovery is effective.
4.3. The Method Named as “Flat, Middle and Push”
As is shown in Introduction, an effective way to spiral recovery in mannedfighter field is named as “flat, middle and push.” Therefore, in order to further verify the effectiveness of proposed method, this method is introduced to the spiral recovery of the UAV.
When UAV is in the spiral, imitate the ground operator to carry on recovery control of the spiral, using the method “flat, middle and push.” The process of this method is shown as follows. First, make the rudder flat. Second, put the control handle to the middle position and then push the control handle to the end in the direction of the spiral. Third, when the vehicle stops rotating, put the control handle to the middle position.
The change curves of UAV’s attitude angle and angular rate are shown from Figure 7(a) to Figure 7(c), and the threedimensional trajectory of spiral recovery is shown in Figure 7(d). From these figures, it can be seen that although this method can realize recovery control of the spiral, the control method is weak, and it takes a long time to recover from the spiral. Thus, it may lead to large loss of the height, making it possible for UAV to crash.
(a) The response curve of attack angle
(b) The response curve of roll angle and yaw angle
(c) The response curve of attitude angular rate
(d) Flight trajectory of UAV’s spiral recovery
4.4. Recover from the Spiral without considering the Control Sequence
In the previous section, the proposed method is verified. However, if we do not consider the control sequence of the controller, what will the effect of spiral recovery be? Regardless of the control sequence, give desired outputs directly, and then, with nonlinear dynamic inversion control laws, the change curves of UAV’s attitude angle and angular rate are shown in Figures 8(a), 8(b), and 8(c). Threedimensional trajectory of spiral recovery is shown in Figure 8(d). From these figures, it can be clearly seen that the change of the angle of attack and the angular rates of attitude are more complex than those in Figures 6 and 7. These figures have shown that if spiral recovery control is carried on regardless of control sequence, the control effect will weaken, thus making the spiral more complex to some degree.
(a) The response curve of attack angle
(b) The response curve of roll angle and yaw angle
(c) The response curve of attitude angular rate
(d) Flight trajectory of UAV’s spiral recovery
To verify the control effect of the three simulations above, the time of spiral recovery of the three comparative simulations above is shown in Table 2. According to the above table, it can be clearly seen that the proposed method in this paper to recover from the spiral is the shortest one. This means this method can make the UAV recover from the spiral autonomously. Finally, the effectiveness of the proposed control method has been verified.

5. Conclusions
In this paper, an autonomous control method for UAV’s spiral recovery is proposed, which is based on nonlinear dynamical inversion. The simulation results have verified that the proposed control method can make the UAV recover from the spiral effectively. Compared with the other two methods, mentioned in the paper, the time of recovery from the spiral is the shortest for this method. Our further study will focus on the robustness of this control method and carry on the practical flight experiment on a small UAV of our research team to verify its usage value.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the National Science Foundation of China (NSFC) under Grant no. 61105012 and the National Aviation Science Foundation of China (NASFC) under Grant no. 20135896027.
References
 D. M. Stipanovic, G. Inalhan, and C. J. Tomlin, “Decentralized overlapping control of a formation of unmanned aerial vehicles,” Automatica, vol. 40, no. 8, pp. 1285–1296, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 D. Gu, “A differential game approach to formation control,” IEEE Transactions on Control Systems Technology, vol. 16, no. 1, pp. 85–93, 2008. View at: Publisher Site  Google Scholar
 Y. Xu, “Nonlinear robust stochastic control for unmanned aerial vehicles,” Journal of Guidance, Control, and Dynamics, vol. 32, no. 4, pp. 1308–1319, 2009. View at: Publisher Site  Google Scholar
 C. J. Ru, R. X. Wei, J. Dai, D. Shen, and L. P. Zhang, “Autonomous reconfiguration control method for UAV's formation based on Nash bargain,” Acta Automatica Sinica, vol. 39, no. 8, pp. 1349–1359, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 L. Jieun, S. K. Hyeong, and K. Youdan, “Formation geometry center based formation controller design using Lyapunov stability theory,” KSAS International Journal, no. 2, pp. 71–76, 2008. View at: Google Scholar
 M. Bryson and S. Sukkarieh, “Cooperative localisation and mapping for multiple UAVs in unknown environments,” in Proceedings of the IEEE Aerospace Conference, pp. 1–12, March 2007. View at: Publisher Site  Google Scholar
 P. Doherty and P. Rudol, “A UAV search and rescue scenario with human body detection and geolocalization,” in Advances in Artificial Intelligence, vol. 4830 of Lecture Notes in Computer Science, pp. 1–13, Springer, Berlin, Germany, 2007. View at: Google Scholar
 S. Sasa, Y. Matsuda, M. Nakadate, and K. Ishikawa, “Ongoing research on disaster monitoring uav at jaxas aviation program group,” in Proceedings of the SICE Annual Conference, pp. 978–981, 2008. View at: Google Scholar
 Q.J. Fan and Z. Yang, “Research status of coordinated formation flight control for multiUAVs,” Acta Aeronautica et Astronautica Sinica, vol. 30, no. 4, pp. 683–691, 2009. View at: Google Scholar
 R. Loh, Y. Bian, and T. Roe, “UAVs in civil Airspace: safety requirements,” IEEE Aerospace and Electronic Systems Magazine, vol. 24, no. 1, pp. 5–17, 2009. View at: Publisher Site  Google Scholar
 L. Pallottino, E. M. Feron, and A. Bicchi, “Conflict resolution problems for air traffic management system solved with mixed integer programming,” IEEE Transactions on Intelligent Transportation Systems, vol. 3, no. 1, pp. 1756–1766, 2002. View at: Publisher Site  Google Scholar
 M. TingLiang, “The application discussion of spiral technology on aerial acrobatics,” Technology & Market, vol. 18, no. 7, pp. 91–92, 2011. View at: Google Scholar
 C. Shuqin, Z. Wencheng, and X. Jianhui, “Control law designed for UAV to step out spiral motion when rudder deadlock happens,” Electronics Optics & Control, vol. 17, no. 9, pp. 57–59, 2010. View at: Google Scholar
 J. Yang, “Faulttolerant control of a class of asynchronous sequential machines with permanent faults,” Automatica, vol. 50, no. 3, pp. 989–993, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 H. Li, J. Yu, C. Hilton, and H. Liu, “Adaptive slidingmode control for nonlinear active suspension vehicle systems using TS fuzzy approach,” IEEE Transactions on Industrial Electronics, vol. 60, no. 8, pp. 3328–3338, 2013. View at: Publisher Site  Google Scholar
 H. Zhang, J. Wang, and Y. Shi, “Robust ${H}_{\infty}$ slidingmode control for Markova jump systems subject to intermittent observations and partially known transition probabilities,” System & Control letters, vol. 62, no. 12, pp. 1114–1124, 2013. View at: Google Scholar
 H. Li, X. Jing, and H. R. Karimi, “Outputfeedbackbased H∞ control for vehicle suspension systems with control delay,” IEEE Transactions on Industrial Electronics, vol. 61, no. 1, pp. 436–446, 2014. View at: Publisher Site  Google Scholar
 H. Zhang, R. Wang, J. Wang, and Y. Shi, “Robust finite frequency ${H}_{\infty}$ staticoutputfeedback control with application to vibration active control of structural systems,” Mechatronics, vol. 24, no. 4, pp. 354–366, 2014. View at: Google Scholar
 H. Zhang and J. Wang, “Combined feedbackfeedforward tracking control for networked control systems with probabilistic delays,” Journal of the Franklin Institute, vol. 351, no. 6, pp. 3477–3489, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Liu, F. Ding, and Y. Shi, “An efficient hierarchical identification method for general dualrate sampleddata systems,” Automatica, vol. 50, no. 3, pp. 962–970, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Karimi, “Optimal vibration control of vehicle enginebody system using Haar functions,” International Journal of Control, Automation and Systems, vol. 4, no. 6, pp. 714–724, 2006. View at: Google Scholar
 H. R. Karimi, “Robust ${H}_{\infty}$ filter design for uncertain linear systems over network with networkinduced delays and output quantization,” Modeling, Identification and Control, vol. 30, no. 1, pp. 27–37, 2009. View at: Publisher Site  Google Scholar
 H. Zhang, X. Zhang, and J. Wang, “Robust gainscheduling energytopeak control of vehicle lateral dynamics stabilization,” Vehicle System Dynamics, vol. 52, no. 3, pp. 309–340, 2014. View at: Google Scholar
 H. Zhang, S. Yang, and W. Junmin, “On energytopeak filtering for nonuniformly sampled nonlinear systems: a Markovian jump system approach,” IEEE Transactions on Fuzzy Systems, vol. 22, no. 1, pp. 1063–6706, 2014. View at: Google Scholar
 B. T. Clough, “Unmanned aerial vehicles: autonomous control challenges, a researcher’s perspective,” in Proceedings of the 2nd AIAA “Unmanned Unlimited” Systems, Technologies, and Operations—Aerospace, San Diego, Calif, USA, September 2003. View at: Google Scholar
 S. Sieberling, Q. P. Chu, and J. A. Mulder, “Robust flight control using incremental nonlinear dynamic inversion and angular acceleration prediction,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 6, pp. 1732–1742, 2010. View at: Publisher Site  Google Scholar
 Y. Chen, X. Dong, and J. Xue, “Constrained adaptive control allocation for multieffector flight control system,” Systems Engineering and Electronics, vol. 33, no. 5, pp. 1118–1123, 2011. View at: Publisher Site  Google Scholar
 J. P. Dyhr, K. A. Morgansen, T. L. Daniel, and N. J. Cowan, “Flexible strategies for flight control: an active role for the abdomen,” The Journal of Experimental Biology, vol. 216, no. 9, pp. 1523–1536, 2013. View at: Publisher Site  Google Scholar
 Y. Shu and S. Tang, “Integrated robust dynamic inversion design of missile guidance and control based on nonlinear disturbance observer,” in Proceeding of the 4th International Conference on Intelligent HumanMachine Systems and Cybernetics (IHMSC '12), pp. 42–45, Nanchang, China, August 2012. View at: Publisher Site  Google Scholar
 Y. Han and Y. Fan, “Control system of automatic landing of UAV based on nonlinear dynamic inversion,” Acta Aeronautica et Astronautica Sinica, vol. 29, pp. 66–72, 2008. View at: Google Scholar
Copyright
Copyright © 2014 ChangJian Ru and RuiXuan Wei. 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.