Research Article | Open Access

Wenru Fan, Bailing Tian, "Adaptive Multivariable Super-Twisting Sliding Mode Controller and Disturbance Observer Design for Hypersonic Vehicle", *Mathematical Problems in Engineering*, vol. 2016, Article ID 5291912, 9 pages, 2016. https://doi.org/10.1155/2016/5291912

# Adaptive Multivariable Super-Twisting Sliding Mode Controller and Disturbance Observer Design for Hypersonic Vehicle

**Academic Editor:**Mohamed Djemai

#### Abstract

A multivariable super-twisting sliding mode controller and disturbance observer with gain adaptation, chattering reduction, and finite time convergence are proposed for a generic hypersonic vehicle where the boundary of aerodynamic uncertainties exists but is unknown. Firstly, an input-output linearization model is constructed for the purpose of controller design. Then, the sliding manifold is designed based on the homogeneity theory. Furthermore, an integrated adaptive multivariable super-twisting sliding mode controller and disturbance observer are designed in order to achieve the tracking for step changes in velocity and altitude. Finally, some simulation results are provided to verify the effectiveness of the proposed method.

#### 1. Introduction

Hypersonic vehicles are intended to be a reliable and cost-effective technology for access to space. During the past decades, a considerable effort has been made by the US Air Force and NASA to further their development. However, the design of control system for hypersonic vehicles is a challenging work due to high speed flight which causes the vehicle to be very sensitive of changes in flight conditions. In addition, the vehicle suffers from severe aerodynamic uncertainties which make the controller design more difficult. Despite the difficulties mentioned above, much effort has been done to develop advanced control technologies for hypersonic vehicle in the past few years.

For the design of control systems for hypersonic vehicles based on linearized dynamical models, several results are available in the literature. For example, Schmidt employed classic and multivariable linear control approach [1] and Groves et al. applied Linear Quadratic Regulator (LQR) technique [2], while Hughes used as well as Linear Parameter Varying (LPV) method [3], and Sigthorsson et al. used implicit model-following control methods [4] to design controller for a linearized hypersonic vehicle model at a specified trim condition. Based on the methods, flight control design is carried out by linearizing the system at a series of operating points and designing separate controllers at each of these points. Finally, the overall flight control system is realized in the philosophy of gain scheduling where the individual gains are interpolated online with respect to some meaningful parameters such as dynamic pressure and Mach number. However, the number of required gains to be designed and scheduled within the controller becomes very large in order to cover different flight missions. In addition, the method involves the lack of guaranteed global robustness, performance, and especially stability [5].

As far as nonlinear control design is concerned, the back-stepping [6], trajectory linearization control [7], robust inversion-based technique [8], sequential loop closure [9], and adaptive control [10] have been developed for hypersonic vehicle in order to improve the robustness and control precision. Although many nonlinear control techniques have been proposed during the last decades, the sliding mode control (SMC) remains a key choice in handing the system with bounded uncertainties and disturbances due to its robustness [11–13]. However, one of the issues of traditional SMC is control chattering caused by high-frequency control switching which restricts its application in real system. Xu et al. [14] combined the adaptive technique with the SMC to design an adaptive sliding mode controller for a generic hypersonic vehicle to track the step commands in velocity and altitude while requiring limited state information. In the method, the behavior of continuous control is achieved via boundary layer technique at the sacrifice of robustness and tracking accuracy to external disturbances and model uncertainties. In addition, it should be noted that the control law designed based on this method is asymptotically stable which means that the convergence rate is at best exponential with infinite settling time. It is obvious that the control law with finite time convergence is more desirable due to the fact that the closed-loop system under finite time control usually demonstrates higher accuracy and better disturbance rejection properties [15]. Preserving the SMC features, the high order sliding mode (HOSM) technique is capable of removing the chattering and improving accuracy and convergent rate [16]. A potential disadvantage of the HOSM to be used in hypersonic vehicle is that the upper bound of uncertainty has to be known exactly in advance. However, it may be difficult to obtain the boundary prior to fly for hypersonic vehicle due to its complex characteristics. Subsequently, the adaptive high order sliding mode control schemes [17, 18] are proposed for hypersonic vehicle where an adaptive law is designed to estimate the uncertainty’s upper bound. In the method, it only requires that the uncertainty is bounded which is a mild assumption for most practical system. However, these methods mentioned above are worst-case-based design and involve the issue of overconservation. Disturbance observer-based control methods provide an effective way to address the issue [19–21]. In the methods, a nonlinear disturbance observer (NDOB) is employed to estimate the uncertainty as well as external disturbance and then a nominal controller is designed for hypersonic vehicle. As a result, the methods obtain not only promising robustness and disturbance rejection performance but also nominal performance recovery. Nevertheless, an assumption of the upper bound of uncertainty has to be known in advance in order to successfully design the NDOB. As noted earlier, this condition may not be available in practice. The first motivation of the research is to propose a practical robust control scheme which is independent on the upper boundary of uncertainty. In addition, it should preserve the advantages of disturbance observer-based methods in [19, 20], such as disturbance rejection, nominal performance recovery, and finite time convergence. Finally, it should be noted that all the methods mentioned above are designed based on single-variable control scheme where only one control input is included. To design multivariable flight control system, the single-variable control scheme is used via decoupling the system into multiple single-variable systems. As pointed in [22], if the multivariable control system is designed directly using the multivariable control scheme, it would have improved chattering reduction property which is the second motivation of the research.

The paper is organized as follows. In Section 2 the vehicle model is introduced and the control objective is stated. The design and stability analysis of adaptive multivariable super-twisting sliding mode controller and disturbance are presented in Section 3. Finally, simulation results are discussed in Section 4, and conclusions are provided in Section 5.

#### 2. Problem Formulation

##### 2.1. Hypersonic Vehicle Model

The longitudinal dynamic model of a generic hypersonic vehicle developed at NASA Langley Research Center is used here depicted by [14]where is control input, fuel equivalence ratio. The state vector for system (1) and (2) is which denotes the velocity, flight path angle (FPA), altitude, angle of attack (AOA), pitch rate, and throttle setting and its time derivative, respectively. Equation (2) is second-order engine dynamics with damping coefficient and natural frequency . The variable denotes the distance between the vehicle and the center of earth. At the trimmed conditions, ft/s, ft, rad, rad, and rad, the lift force , drag force , thrust force , and pitch moment are calculated aswith here is elevator deflection.

##### 2.2. Control Objective

The objective of the research is to determine the control inputs which makes the system output track the desired commands, velocity , and altitude , in finite time in the presence of the following bounded but unknown perturbations:

#### 3. Adaptive Multivariable Super-Twisting Sliding Mode Control Scheme

##### 3.1. Preliminary

In the subsection, several useful lemmas to be used in the design of controller and disturbance observer for hypersonic vehicle are recalled.

Lemma 1 (see [23]). *Consider a non-Lipschitz system in the form of , , and suppose that there exist a Lyapunov function and real numbers , , and , such that is positive for any nonzero and inequality holds. Then, the origin is fast finite time stable and the settling time, depending on the initial state , given by .*

Lemma 2 (see [24]). *Suppose that there exist constant and positive constants such that polynomial is Hurwitz. Then, the integrator chain system is finite time stable under the feedback , where denotes the signum function and satisfies with , .*

##### 3.2. Main Results

Inspired from the research in [22], the single-variable control scheme [25, 26] is extended to the following multivariable super-twisting sliding mode control scheme.

Theorem 3. *Consider the following multivariable system: where and is adaptive gain to be developed soon and suppose that the perturbation satisfies the condition , where the finite boundary exists but is unknown. Then, one has the following claims.**(1) If the gain in (6) is designed aswith positive scalar satisfyingthen vectors and converge to zero in finite time.**(2) For any positive constants , the states of system (6) are globally bounded for any bounded perturbation .**(3) The vectors and converge to zero for any positive parameter if there is no perturbation for system (6), which means in (6).*

*Proof. *The first claim is a direct result of our recent work [27]. Following the conclusion in [27], it follows that there exists a continuous and positive definite Lyapunov functionfor system (6) such thatwith , , , and some bounded constants , , , and provided in [27]. From the definition of adaptive law in (7), it can be observed that the items and in (10) will be positive in finite time for arbitrary constants , , , and and boundary perturbation since . After that, one has implying equivalently in finite time which proves the first claim of Theorem 3.

In order to prove the second claim, we can set , which implies . In this case, (10) is reduced toNext, the item in (11) can be used to dominate for sufficient large . To this end, we rewrite inequality (11) aswhere is a positive scalar satisfying . It is obvious thatIn view of inequality (13) and theorem 4.18 in [28], it can be concluded that the solutions starting in the set with in (13) will remain therein for all future time since is negative on the boundary . On the other hand, if the solutions start outside the set , will decrease monotonically until the solution enters the set . After that, the solution cannot leave the set . Hence, the states of system (6) are globally bounded for any bounded perturbation .

With respect to the third claim, we can set , , and . Then (11) can be reduced to . It follows from Lemma 2 that the vectors and converge to zero for any positive parameter if holds in (6).

Corollary 4. *Consider the following general multivariable uncertainty system:where , are state and control vector and , represent nominal and uncertainty parts, respectively. It is assumed that uncertainty is continuously differentiable with , where exists but is unknown. Furthermore, the following conclusions can be obtained.**(1) If an adaptive multivariable disturbance observer is designed,with and being chosen according to (7) and (8). Then, can be estimated through in finite time.**(2) If multivariable controller is designed, where is disturbance estimation generated from (15) and is arbitrary positive constant. Then, in (14) converges to zero in finite time under the controller (16) and disturbance observer (15).*

*Proof. *Combining (14) and the definitions and , (15) is converted intoBased on Theorem 3, it can be concluded that the estimation error vectors and converge to zero in finite time if the parameters are chosen according to (7) and (8), which implies that can be reconstructed via in finite time. This completes the proof of the first part of Corollary 4. Furthermore, the closed-loop system (14) under controller (16) and disturbance observer (15) has the following form:For convenience, a new vector is introduced so that (18) can be rewritten as , . It follows from the second conclusion in Theorem 3 that this system is input-to-state stable with viewed as input. From estimation error dynamics in (17), it can be observed that is bounded. Therefore, the state in (18) is bounded for arbitrary bounded . Furthermore, it follows from the first conclusion in Theorem 3 that there exists a finite time such that the estimation error for . After that, system (18) is reduced to nominal system: , with bounded initial value . Then, based on the third conclusion in Theorem 3, it can be seen that vector converges to zero in finite time for arbitrary positive scalar . This completes the proof of the second part of Corollary 4.

*Remark 5. *In practice cannot be zero exactly. Therefore, will increase unboundedly. To avoid this drawback, the adaptive law in (7) can be modified as if ; otherwise , where is a positive scalar. The idea of the adaptive law can be interpreted as follows: the gain increases when unacceptably deviates from zero and it stops increasing when is driven into the regime .

*Remark 6. *The introduction of an adaptive disturbance observer in control scheme has two reasons. First, disturbance observer is able to provide estimation for uncertainties and then the estimated value can be used in control input to compensate the actual uncertainties in plant. In this case, the baseline controller with small control gains can be used to stabilize the closed-loop system, which implies the properties of chattering reduction and nominal performance recovery are obtained (see [21]). Furthermore, the gain adaptation is used in disturbance observer in order to avoid dependency on upper boundary of uncertainty which may be difficult to obtain in practice.

##### 3.3. Controller and Disturbance Observer Design

It follows from the results in [14] that the input-output linearized model can be obtained via differentiating velocity and altitude three and four times, respectively, depicted bywhere and are bounded perturbation induced by model parameters uncertainties. All the parameters in (19) can be referred to [14]. For brevity, the details are omitted. Furthermore, the sliding mode manifolds for velocity and altitude are designed as where and . It follows the results in Lemma 2 that and in finite time if the conditions hold and the sliding manifold parameters and are chosen properly. To this end, taking the derivative of and along the trajectories in (20) yields The uncertainties vector can be interpreted as total perturbations which may be caused by model parameter uncertainty and unknown external disturbance. In the simulation, is induced by uncertainties in (5). Furthermore, it is assumed that is continuously differentiable such that holds, where is a finite value but unknown. In fact, is a function of system states. In real applications, the derivatives of system states should be bounded in order to guarantee safe and stable flight for vehicle, which can be confirmed through the results to be provided in Figures 2(d) and 2(e). Consequently, it is a mild assumption. It is obvious that system (21) can be considered as special case for the general multivariable system in (14). Consequently, the adaptive multivariable controller and disturbance observer can be easily designed based on Corollary 4. For brevity, the details are omitted here.

*Remark 7. *It should be noted that the missing derivatives of and in sliding manifolds and control law can be estimated online using an arbitrary order robust finite time convergent differentiator [29].

#### 4. Simulation Results Analysis

##### 4.1. Parameters Setting

The parameters used in the simulation are summarized as follows. The sliding manifold parameters in (20) are taken as , , , , , and . The other parameters and are calculated according to Lemma 2. The parameters used in disturbance observer and controller are chosen as , , and . The tolerance constant mentioned in Remark 5 is 0.001. The step commands for velocity 100 ft/s and altitude 1000 ft are used as the reference commands to be tracked.

##### 4.2. Results Discussion

*Case 1. *The maximum positive uncertainties , , , and are used here. In the case, the simulation results are provided in Figures 1 and 2. Specifically, the response curves for velocity and altitude are shown in Figures 1(a) and 1(b), in which it can be seen that the tracking is achieved in finite time with little overshoot. The results in Figures 1(c)–1(e) show the variations of AOA, FPA, and pitch rate, whereas the behavior of adaptive gain is provided in Figure 1(f). It can be observed from Figure 1(f) that adaptive gain increases when exceeds tolerance constant . In addition, the sliding manifolds for velocity and altitude in (20) are plotted in Figures 2(a) and 2(b). From that, it can be seen that the proposed control method is able to guide and to zero in the presence of bounded but unknown uncertainty. Figures 2(e) and 2(f) show the curves of derivatives in velocity and altitude errors, respectively. The results illustrate that not only the errors for velocity and altitude but also their derivatives converge to zero in finite time. The control variables, fuel equivalence ratio, and elevator deflection are given in Figures 2(c) and 2(f) where the results demonstrate that the control chattering is reduced effectively. In addition, the disturbance estimation values for and are added in Figure 3. Since the actual values for and are difficult to obtain, they are not provided in the simulation. However, the good tracking performance for velocity and altitude in Cases 1 and 2 to be given demonstrates the effectiveness of the control scheme in disturbance rejection.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(a)**

**(b)**

*Case 2. *In order to further verify the effectiveness of the proposed control scheme in disturbance suppression, the Monte Carlo simulation with 200 tests is conducted. The random uncertainty is added in the simulation according to the conditions in (5). All other parameters are the same as that provided in Case 1. In this case, the simulation results are provided in Figure 4. From Figures 4(a) and 4(b), it can be found that the velocity and altitude tracking can be achieved excellently even in the case of random uncertainty. The sliding manifolds for velocity as well as altitude and the corresponding control variables are also plotted in Figure 4. From Figures 4(c) and 4(d), it can be seen that the fuel equivalence ratio and elevator deflection have different steady-state values due to the effect of uncertainties. In fact, the magnitudes of fuel equivalence ratio and elevator deflection have to increase when positive drag force coefficient and negative lift force coefficient uncertainties are included in order to keep the stability of the system. Finally, the simulation results shown in this section demonstrate the effectiveness of the proposed control scheme in achieving the tracking for velocity and altitude.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 5. Conclusion

The finite time tracking for hypersonic vehicle with boundary but unknown perturbations is discussed based on adaptive multivariable super-twisting control scheme. The finite time stability of the integrated controller and disturbance observer is guaranteed under the multivariable control architecture via Lyapunov analysis. The features of the proposed multivariable control scheme are gain adaptation, chattering reduction, nominal performance recovery, and finite time convergence. Finally, the Monte Carlo simulation with 200 tests is provided to demonstrate the effectiveness of the proposed control scheme.

#### Competing Interests

The authors declare that they have no competing interests.

#### Acknowledgments

This work has been supported by National Natural Science Foundation of China (61301244), Tianjin Research Program of Application Foundation and Advanced Technology (12JCZDJC30300), Aeronautical Science Foundation of China Supported by Science and Technology on Aircraft Control Laboratory (20125848004), Independent Innovation Fund of Tianjin University (2013XQ-0022), and Public Fund Supported by Tianjin Key Laboratory of Process Measurement and Control (TKLPMC-201315).

#### References

- D. K. Schmidt, “Optimum mission performance and multivariable flight guidance for airbreathing launch vehicles,”
*Journal of Guidance, Control, and Dynamics*, vol. 20, no. 6, pp. 1157–1164, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. P. Groves, D. O. Sigthorsson, A. Serrani et al., “Reference command tracking for a linearized model of an air-breathing hypersonic vehicle,” in
*Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit*, AIAA-2005-6144, San Francisco, Calif, USA, August 2005. View at: Google Scholar - H. D. Hughes,
*LPV Control for the Longitudinal Dynamics of a Flexible Air-Breathing Hypersonic Vehicle, Mechnaical Engineerin*, North Carolina State University, Raleigh, NC, USA, 2010. - D. O. Sigthorsson, P. Jankovsky, A. Serrani, S. Yurkovich, M. A. Bolender, and D. B. Doman, “Robust linear output feedback control of an airbreathing hypersonic vehicle,”
*Journal of Guidance, Control, and Dynamics*, vol. 31, no. 4, pp. 1052–1066, 2008. View at: Publisher Site | Google Scholar - D. J. Leith and W. E. Leithead, “Survey of gain-scheduling analysis and design,”
*International Journal of Control*, vol. 73, no. 11, pp. 1001–1025, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Q. Zong, Y. H. Ji, F. L. Zeng, and H. Liu, “Output feedback back-stepping control for a generic hypersonic vehicle via small-gain theorem,”
*Aerospace Science and Technology*, vol. 23, no. 1, pp. 409–417, 2012. View at: Publisher Site | Google Scholar - T. A. Adami, J. J. Zhu, M. A. Bolender, D. B. Doman, and M. Oppenheimer, “Flight control of hypersonic scramjet vehicles using a differential algebraic approach,” in
*Proceedings of the AIAA Guidance, Navigation, and Control Conference*, pp. 3680–3698, Keystone, Colo, USA, August 2006. View at: Publisher Site | Google Scholar - Q. Wang and R. F. Stengel, “Robust nonlinear control of a hypersonic aircraft,”
*Journal of Guidance, Control, and Dynamics*, vol. 23, no. 4, pp. 577–585, 2000. View at: Publisher Site | Google Scholar - L. Fiorentini, A. Serrani, M. A. Bolender, and D. B. Doman, “Robust nonlinear sequential loop closure control design for an air-breathing hypersonic vehicle model,” in
*Proceedings of the American Control Conference*, pp. 3458–3463, Seattle, Wash, USA, June 2008. View at: Publisher Site | Google Scholar - L. Fiorentini, A. Serrani, M. A. Bolender, and D. B. Doman, “Nonlinear robust adaptive control of flexible air-breathing hypersonic vehicles,”
*Journal of Guidance, Control, and Dynamics*, vol. 32, no. 2, pp. 401–416, 2009. View at: Publisher Site | Google Scholar - X. H. Yu and O. Kaynak, “Sliding-mode control with soft computing: a survey,”
*IEEE Transactions on Industrial Electronics*, vol. 56, no. 9, pp. 3275–3285, 2009. View at: Publisher Site | Google Scholar - S. D. Gennaro, J. R. Dominguez, and M. A. Meza, “Sensorless high order sliding mode control of induction motors with core loss,”
*IEEE Transactions on Industrial Electronics*, vol. 61, no. 6, pp. 2678–2689, 2014. View at: Publisher Site | Google Scholar - Y. Feng, X. H. Yu, and F. L. Han, “High-order terminal sliding-mode observer for parameter estimation of a permanent-magnet synchronous motor,”
*IEEE Transactions on Industrial Electronics*, vol. 60, no. 10, pp. 4272–4280, 2013. View at: Publisher Site | Google Scholar - H. Xu, M. D. Mirmirani, and P. A. Ioannou, “Adaptive sliding mode control design for a hypersonic flight vehicle,”
*Journal of Guidance, Control, and Dynamics*, vol. 27, no. 5, pp. 829–838, 2004. View at: Publisher Site | Google Scholar - S. H. Li, S. H. Ding, and Q. Li, “Global set stabilisation of the spacecraft attitude using finite-time control technique,”
*International Journal of Control*, vol. 82, no. 5, pp. 822–836, 2009. View at: Publisher Site | Google Scholar | MathSciNet - A. Levant, “Homogeneity approach to high-order sliding mode design,”
*Automatica*, vol. 41, no. 5, pp. 823–830, 2005. View at: Publisher Site | Google Scholar | MathSciNet - Q. Zong, J. Wang, and Y. Tao, “Adaptive high-order dynamic sliding mode control for a flexible air-breathing hypersonic vehicle,”
*International Journal of Robust and Nonlinear Control*, vol. 23, no. 15, pp. 1718–1736, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. L. Tian, W. R. Fan, Q. Zong, J. Wang, and F. Wang, “Adaptive high order sliding mode controller design for hypersonic vehicle with flexible body dynamics,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 357685, 11 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - J. Yang, S. Li, C. Sun, and L. Guo, “Nonlinear-disturbance-observer-based robust flight control for airbreathing hypersonic vehicles,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 49, no. 2, pp. 1263–1275, 2013. View at: Publisher Site | Google Scholar - H. Sun, S. Li, and C. Sun, “Finite time integral sliding mode control of hypersonic vehicles,”
*Nonlinear Dynamics*, vol. 73, no. 1-2, pp. 229–244, 2013. View at: Publisher Site | Google Scholar | MathSciNet - J. Yang, S. Li, and X. Yu, “Sliding-mode control for systems with mismatched uncertainties via a disturbance observer,”
*IEEE Transactions on Industrial Electronics*, vol. 60, no. 1, pp. 160–169, 2013. View at: Publisher Site | Google Scholar - I. Nagesh and C. Edwards, “A multivariable super-twisting sliding mode approach,”
*Automatica*, vol. 50, no. 3, pp. 984–988, 2014. View at: Publisher Site | Google Scholar | MathSciNet - S. Y. Khoo, L. H. Xie, S. K. Zhao, and Z. Man, “Multi-surface sliding control for fast finite-time leader-follower consensus with high order SISO uncertain nonlinear agents,”
*International Journal of Robust and Nonlinear Control*, vol. 24, no. 16, pp. 2388–2404, 2014. View at: Publisher Site | Google Scholar - S. P. Bhat and D. S. Bernstein, “Geometric homogeneity with applications to finite-time stability,”
*Mathematics of Control, Signals, and Systems*, vol. 17, no. 2, pp. 101–127, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. B. Shtessel, M. Taleb, and F. Plestan, “A novel adatpve-gain supertwisting sliding mode controller: methodology and application,”
*Automatica*, vol. 48, no. 5, pp. 759–769, 2012. View at: Publisher Site | Google Scholar | MathSciNet - J. X. Liu, S. Laghrouche, M. Harmouche, and M. Wack, “Adaptive-gain second-order sliding mode observer design for switching power converters,”
*Control Engineering Practice*, vol. 30, pp. 124–131, 2014. View at: Publisher Site | Google Scholar - B. Tian, L. Yin, and H. Wang, “Finite-time reentry attitude control based on adaptive multivariable disturbance compensation,”
*IEEE Transactions on Industrial Electronics*, vol. 62, no. 9, pp. 5889–5898, 2015. View at: Publisher Site | Google Scholar - H. K. Khalil,
*Nonlinear Systems*, Prentice Hall Press, New York, NY, USA, 3rd edition, 2001. - A. Levant, “Higher-order sliding modes, differentiation and output-feedback control,”
*International Journal of Control*, vol. 76, no. 9-10, pp. 924–941, 2003. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2016 Wenru Fan and Bailing Tian. 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.