Equilibrium Positions for UAV Flight by Dynamic Soaring

Abstract

Dynamic soaring is a special flying technique designed to allow UAVs (unmanned aerial vehicles) to extract energy from wind gradient field and enable UAVs to increase the endurance. In order to figure out the energy-extraction mechanisms in dynamic soaring, a noninertial wind relative reference frame of aircraft is built. In the noninertial frame, there is an inertial force which is created by gradient wind field. When the wind gradient ([subscript]GW[/subscript] ) and the components of airspeed ([subscript]vz[/subscript] [subscript]vx[/subscript] ) are positive, inertial force (F) makes positive work to the aircraft. In the meantime, an equilibrium position theory of dynamic soaring is proposed. At the equilibrium positions, the increased potential energy is greater than the wasted kinetic energy when the aircraft is flying upwards. The mechanical energy is increased in this way, and the aircraft can store energy for flight. According to the extreme value theory, contour line figures of the maximum function and the component of airspeed ([subscript]vz[/subscript] ) are obtained to find out the aircraft's lifting balance allowance in dynamic soaring. Moreover, this equilibrium position theory can also help to conduct an aircraft to acquire energy from the environment constantly.

Bingjie Zhu 1 and Zhongxi Hou 1 and Shangqiu Shan 1 and Xinzhu Wang 2

Academic Editor:Hyochoong Bang

1, College of Aerospace Sciences and Engineering, National University of Defense Technology, Changsha, Hunan 410073, China
2, Chongqing Key Laboratory of Heterogeneous Material Mechanics, Department of Engineering Mechanics, Postdoctoral Station of Mechanics, Chongqing University, Chongqing 400040, China

Received 23 June 2015; Revised 22 September 2015; Accepted 30 September 2015; 22 October 2015

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.

1. Introduction

Long endurance flight is a key design criterion for unmanned aerial vehicle [1]. The UAVs' flight endurance can be improved by the progress in power capacity and aerodynamic and engine design [2], but not all such means are equal. Another important alternative for enhancing UAV flight endurance is to seek energy from flight environment. In nature, there are many large-sized birds which make use of environment energy to enhance their flight ability, such as wandering albatross who extracts energy from horizontally moving air and transfers energy to itself. The achieved energy enables wandering albatrosses to fly continuously without flapping their wings [3, 4]. The flight manoeuvre that wandering albatross performs is referred to as dynamic soaring [5]. Early in 1883, Lord Rayleigh published "The Soaring of Birds" in Nature [6], which is regarded as the first paper about dynamic soaring. Dynamic soaring is possible when there is horizontally moving air in the flight environment. Horizontally moving air is also called gradient wind, in which wind speed varies with altitude. In this paper, both the size of the wandering albatrosses with its 3.3 m wing span and its flight performance are taken as a reference model against the dynamic soaring performance of a small UAV.

Dynamic soaring is an important flight technique designed to allow aircraft to extract energy from atmosphere. Based on Lord Rayleigh's observations of albatrosses in the south Atlantic, more recent publications are concerned with wind estimated and dynamic soaring mechanism. In order to estimate wind field for autonomous dynamic soaring, Lawrance et al. [7, 8] provided a method for taking direct observations of the wind during flight and creating a wind to direct future exploration. The wind estimation is the precondition for dynamic soaring. In this paper, a mathematical model is built for the wind field. As for dynamic soaring mechanism, Grenestedt [9] provided insight into how and where energy is extracted from the wind during dynamic soaring. The effects of some aircraft parameters such as mass, drag area, and wing span were discussed. Bower [5] studied the energy transfer mechanisms for a vehicle flying in a spatially and temporally varying wind field. Deittert et al. [10] investigated optimal trajectories for minimal and maximal wind conditions, and the likelihood of favorable winds was predicted based on long-term weather statics and knowledge of minimal and maximal permissible wind strengths. Cone [11] analyzed the dynamic soaring flight of albatross with mathematical models. Peter [12] draw a conclusion that, in an energy neutral cycle of dynamic soaring, whether the vehicle returns to initial velocity and height with no power input depends only upon the maximum lift-to-drag ratio of the vehicle and the wind speed variation.

Scores of research findings reveal the significant meaning of dynamic soaring. Based on the achievements of these findings, some new ideas are studied in this paper. The authors first extend the idea of Grenestedt [9] by investigating the equilibrium position of dynamic soaring and providing insight into energy variation in equilibrium position. The same as other flight styles, they applied forces on the aircraft for dynamic soaring varying continuously. Then there should be some positions of forces balance in the course of flight. At these points, the aircraft stays in a certain speed, with no kinetic energy loss. Of course, if all the points could chain a closed route, the permanent flying would come true when the aircraft flies along the route.

The rest of the paper is organized as follows: the methodology is described in Section 2, including the wind gradient and the mathematical models for the UAV. Section 3 is simulation results and associated discussion. Finally, the concluding remarks are made in Section 4.

2. Mathematical Formulations of the Problem

Inspired by the forces analysis methods in dynamic soaring suggested by Bower [5] and Deittert et al. [10], the problem of equilibrium position for dynamic soaring is analyzed intensively in this paper. The parameters that can influence the equilibrium position include the range of airspeed, the size of wind gradient, and other parameters. Meanwhile, the maximum value of energy acquired in the equilibrium point is analyzed. Once the parameters of an aircraft are known, the wind gradient suited for the equilibrium position of dynamic soaring can be determined. Similarly, in a known wind field, the type of aircraft that has equilibrium point in dynamic soaring can be found out. A set of computational methods for equilibrium position of dynamic soaring is established. The purpose of this section is to provide a theoretical guidance for the optimal dynamic soaring trajectory by means of the analysis of equilibrium position of dynamic soaring. Due to no kinetic energy loss in the position of forces balance, the mechanical energy increases with increasing value of potential energy of the aircraft in the course of upward flying.

2.1. The Wind Field Model

It is assumed that the aircraft flies over a flat surface and the wind field is steady. There is a horizontal shear flow of the air in the three-dimensional space. The flow of the air causes the wind speed to increase with increasing height above the ground. The air density is assumed to be constant. The wind blows along the negative direction of [figure omitted; refer to PDF] -axis, the positive [figure omitted; refer to PDF] -axis points to the sky, and the direction of [figure omitted; refer to PDF] -axis is dependent on right-handed rule. The relation between wind speed [figure omitted; refer to PDF] and height above the surface [figure omitted; refer to PDF] is defined as follows [13, 14]: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the wind speed at height [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is the wind speed at the reference height [figure omitted; refer to PDF] . The variable [figure omitted; refer to PDF] is the aerodynamic roughness length or the roughness factor. It is an experimentally derived constant that accounts for the kind of surface over which the wind is blowing. Typically, a higher value of the roughness length indicates more obstructions on the surface such as trees and buildings. By differentiating (1) with respect to altitude, the following equation can be obtained: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the gradient of the wind field. The variation trends of wind profile and wind gradient profile are depicted in Figure 1.

Figure 1: Wind speed and wind gradient.

[figure omitted; refer to PDF]

Figure 1 shows the wind profile and wind gradient profile for value of the roughness coefficient with [figure omitted; refer to PDF] m/s, [figure omitted; refer to PDF] m, and [figure omitted; refer to PDF] m.

2.2. Flying Model

According to the observed data of dynamic soaring [4], the albatross climbs with headwind, and the wind speed increases with the increasing of flight height. After turning in high altitude, the albatross glides with leeward, and the wind speed decreases with the descending of flight height. An illustration of a typical dynamic soaring trajectory is shown in Figure 2.

Figure 2: The illustration trajectory of dynamic soaring.

[figure omitted; refer to PDF]

Compared with Earth fixed frame, the wind relative reference frame of aircraft is not inertial frame in the wind field. So there is an inertial force [figure omitted; refer to PDF] that works on the aircraft. A noninertial frame is defined, which is the wind relative reference frame of aircraft. In the noninertial frame, there are four applied forces accounted for the aircraft: lift force ( [figure omitted; refer to PDF] ), drag force ( [figure omitted; refer to PDF] ), gravitational force ( [figure omitted; refer to PDF] ), and inertial force ( [figure omitted; refer to PDF] ). The inertial force is created by wind, so the direction of the force is the same as the wind direction. The inertial force can be expressed as [10, 15, 16] [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the unit vector of wind field and [figure omitted; refer to PDF] is the airspeed component in the direction of positive [figure omitted; refer to PDF] -axis. So a point mass flight model with attached forces is adopted, and the definitions of forces, speed, and angles used in the model are shown in Figure 3. The UAV is represented by a three-dimensional point-mass model with 3 DOF. From Figure 3, [figure omitted; refer to PDF] is the resultant force of lift ( [figure omitted; refer to PDF] ) and drag ( [figure omitted; refer to PDF] ), [figure omitted; refer to PDF] is the angle between [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , [figure omitted; refer to PDF] is the resultant force of gravity ( [figure omitted; refer to PDF] ) and inertial force ( [figure omitted; refer to PDF] ), and [figure omitted; refer to PDF] is the reverse force of [figure omitted; refer to PDF] . The wind blows in the direction of position [figure omitted; refer to PDF] .

Figure 3: The applied forces for the UAV.

[figure omitted; refer to PDF]

The barycenter kinetics equation of aircraft in the noninertial frame is represented as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the acceleration vector of airspeed, [figure omitted; refer to PDF] is the vector of lift, [figure omitted; refer to PDF] is the vector of drag, and [figure omitted; refer to PDF] is the vector of gravity. The last expression in (4) is the inertial force.

Using the standard formulations for lift and drag force, [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is air density, [figure omitted; refer to PDF] is wing area, [figure omitted; refer to PDF] is lift coefficient, [figure omitted; refer to PDF] is drag coefficient, [figure omitted; refer to PDF] is the airspeed of the aircraft, [figure omitted; refer to PDF] is the parasitic drag coefficient, and [figure omitted; refer to PDF] is the induced drag coefficient. Note that [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the maximum lift over drag ratio, which, together with [figure omitted; refer to PDF] , is the most important aircraft aerodynamic efficiency parameters and fully defines the aerodynamic drag polar.

In the course of flight, the direction of lift is crosscut to drag all the time, so the resultant force ( [figure omitted; refer to PDF] ) of lift and drag can be indicated as [figure omitted; refer to PDF]

From (9), the airspeed ( [figure omitted; refer to PDF] ) of the aircraft is [figure omitted; refer to PDF]

When the value of [figure omitted; refer to PDF] is equal to [figure omitted; refer to PDF] , the relationship between [figure omitted; refer to PDF] and [figure omitted; refer to PDF] can be attained by (10).

As for the energy variation, first of all, considering the theorem of kinetic energy in the noninertial frame, [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the variation of kinetic energy, [figure omitted; refer to PDF] is the identity element of flight trajectory, and [figure omitted; refer to PDF] . The right-hand side of (11) is infinitesimal energy expressions caused by drag force, inertial force, and gravitational force. The lift force ( [figure omitted; refer to PDF] ) is vertical to the direction of motion, so there is no kinetic energy variation of UAV caused by lift force.

If [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] is the airspeed component in the direction of positive [figure omitted; refer to PDF] , (11) then transforms into [figure omitted; refer to PDF]

The integration of (12) can be expressed as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the kinetic energy at height of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is the kinetic energy at height of [figure omitted; refer to PDF] . After the aircraft fly a period of time, if the kinetic energy of aircraft keeps in a constant value, the energy variation of this process could be indicated by the change of altitude difference, and (13) can now be stated to be [figure omitted; refer to PDF] [figure omitted; refer to PDF] is the variation of flight height. From (14), it can be found that the drag ( [figure omitted; refer to PDF] ) and inertial force ( [figure omitted; refer to PDF] ) bring about the variation of altitude or the variation of potential energy. That is to say, as long as there was wind gradient ( [figure omitted; refer to PDF] ) in the wind field, the aircraft could extract energy from the wind field in the suitable [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . In (14), the drag ( [figure omitted; refer to PDF] ) makes negative work all the time, but inertial force ( [figure omitted; refer to PDF] ) does positive or negative work to the aircraft, depending on the specific cases. If [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , in this case, inertial force ( [figure omitted; refer to PDF] ) does positive work to the aircraft. For the wind field [figure omitted; refer to PDF] , the following relation must hold: [figure omitted; refer to PDF]

That is, the aircraft climbs in the positive direction of [figure omitted; refer to PDF] or dives in the negative direction of [figure omitted; refer to PDF] ; the aircraft extracts energy. If the maneuver is reversed, the aircraft loses energy. The mechanical energy of the aircraft is denoted as [figure omitted; refer to PDF] . The variation of mechanical energy is equal to the energy caused by drag and inertial force. The identity element of mechanical energy can be expressed as [figure omitted; refer to PDF] So the unit mass power of mechanical energy is [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the airspeed. Actually, (17) and (14) mean the same things.

2.3. The Balance Points of Dynamic Soaring

In the point of balance, the composition of driving force vectors is zero, and the resultant force of aircraft is zero. The dynamical equation in noninertial system can be expressed as [figure omitted; refer to PDF] Because the ground speed of the aircraft remains constant in the balance point, that is to say, the kinetic energy of the flight system remains constant. When [figure omitted; refer to PDF] , the potential energy of the flight system increases with increase in flight altitude, and the mechanical energy increases at the same time. On the other hand, when [figure omitted; refer to PDF] , the mechanical energy of the aircraft would not always increase, in which case, some studies need to be done in the future. Here, the balance point is called lifting balance point when [figure omitted; refer to PDF] ; while when [figure omitted; refer to PDF] , the balance point is called descending balance point. When [figure omitted; refer to PDF] , there is no balance point, because there is no inertial force to compensate the drag.

Alternatively, (18) can be expressed in another way: [figure omitted; refer to PDF]

Combining (5)-(6), (10), and (19), there is [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Equation (20) describes the relationship between [figure omitted; refer to PDF] and [figure omitted; refer to PDF] in the balance point. Simultaneously, (20) ensured that the size of aerodynamic force is equal to the resultant force of gravitation and inertial force.

Then, the only restriction to airspeed is the bank angle of air velocity vector, which is defined as [figure omitted; refer to PDF] . There is a unique bank angle to enable forces to be balanced; it is necessary to figure out that the bank angle, which made the included angle between air velocity vector ( [figure omitted; refer to PDF] ) and [figure omitted; refer to PDF] , is [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] is the resultant force of gravity ( [figure omitted; refer to PDF] ) and inertial force ( [figure omitted; refer to PDF] ). In Figure 3, the included angle between airspeed vector ( [figure omitted; refer to PDF] ) and [figure omitted; refer to PDF] can be repressed as [figure omitted; refer to PDF]

When [figure omitted; refer to PDF] , in this case, the aerodynamic force appears to get across the inverse direction of the resultant force of gravitation and inertial force. Hence, [figure omitted; refer to PDF]

With the geometrical relationship of forces balanced, one can acquire the space curve of [figure omitted; refer to PDF] , which is obtained by projecting the four-dimensional space [figure omitted; refer to PDF] into three-dimensional space. The fundamental is to find out a suitable aerodynamic force to compensate the resultant force of gravitation and inertial force.

The resultant force of gravitation and inertial force [figure omitted; refer to PDF] , air velocity vector [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] ; combining (20) and (23), there is [figure omitted; refer to PDF] Let [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the wing loading.

From (24) and (25), [figure omitted; refer to PDF] which is subjected to the constraint [figure omitted; refer to PDF] Obviously, there are two roots for (27), a positive root and a negative root. Moreover, the value of positive root is greater than the absolute value of negative root. Thus, there is [figure omitted; refer to PDF] which is one of constraint conditions for the solution of (26).

3. Simulation Results and Discussions

Because there is no closed form solution for the balance points in dynamic soaring, numerical investigations are carried out instead. As dynamic soaring flight is directly inspired by albatross, wandering albatross' parameters used by Sachs [4] are assumed. Based on the architecture of a wandering albatross, an aircraft demonstrator model is used as a reference. The UAV demonstrator has similarities in dimension as well as the performances with albatross. This will ensure that the performances of the aircraft model are achievable by an engineering-designed vehicle. The parameter values of albatross and UAV demonstrator are given in Table 1, where [figure omitted; refer to PDF] is the maximum lift-to-drag ratio. Compared with the albatross, the UAV demonstrator has a lighter mass. Furthermore, the flying ability close to the ground surface of the UAV demonstrator is limited, while albatross can fly with a wing tip touching the water. For the reason of security, the UAV demonstrator has a smaller wingspan. The aerofoil of UAV demonstrator is E214, which makes the vehicle have a better performance than albatross, reaching a maximum lift-to-drag ratio of 21.

Table 1: Albatross and UAV demonstrator parameter values used.

Parameter	Albatross	UAV demonstrator
Mass [kg]	8.5	5.5
Wing area [figure omitted; refer to PDF] [m2 ]	0.65	0.68088
Aspect ratio	16.81	16.81
Span [m]	3.306	2.61
[figure omitted; refer to PDF] [kg/m3 ]	1.22	1.22
[figure omitted; refer to PDF]	0.033	0.033
[figure omitted; refer to PDF]	20	21

From (26), it can be found that the relationship between variables [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is linear and the value of [figure omitted; refer to PDF] is located in the interval [figure omitted; refer to PDF] , so the range of [figure omitted; refer to PDF] can be defined. The extent of albatross' [figure omitted; refer to PDF] and UAV demonstrator's [figure omitted; refer to PDF] is shown in Figures 4 and 5, respectively.

Figure 4: The extent of [figure omitted; refer to PDF] , three typical values of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] are selected, which can tell the relationship between albatross' [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] .

(a) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Figure 5: The extent of UAV demonstrator's [figure omitted; refer to PDF] , three typical values of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] are selected.

(a) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) When [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , 1.0, and 1.5, the extent of [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

In Figures 4 and 5, the black horizontal dotted line represents [figure omitted; refer to PDF] . The solid curve represents the variation of [figure omitted; refer to PDF] . The dotted curve represents [figure omitted; refer to PDF] . From Figures 4 and 5, it can be found that, no matter what set of parameters ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ) are selected, there always exists a balance point in the course of dynamic soaring. In the range of [figure omitted; refer to PDF] , the maximum value of the solid line is greater than one, and solid line and dotted line go through the original point. That is, [figure omitted; refer to PDF] has corresponding value in the area of [figure omitted; refer to PDF] . Furthermore, in the area of [figure omitted; refer to PDF] , the right-hand limit of (28) includes values which are greater than zero, then the expression of right-hand limit should be in the phase of lifting at zero point, and the derived function of right-hand limit should be greater than zero at zero point; that is, [figure omitted; refer to PDF] Hence, [figure omitted; refer to PDF]

In (30), there is not enough evidence to say that there exists balance point in the area of [figure omitted; refer to PDF] . From Figures 4 and 5, it can be found that the maximum value of the right-hand limit of (26) cannot be less than [figure omitted; refer to PDF] ; that is, [figure omitted; refer to PDF] Equation (31) then is converted to [figure omitted; refer to PDF] where [figure omitted; refer to PDF] In order to obtain the maximum value of (32), another function [figure omitted; refer to PDF] is created: [figure omitted; refer to PDF]

The contour line figure of albatross' [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is expressed in Figures 6 and 7.

Figure 6: The contour line of (albatross) [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]

Figure 7: The contour line of (albatross) [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]

The contour line figure of UAV demonstrator's [figure omitted; refer to PDF] and [figure omitted; refer to PDF] is expressed in Figures 8 and 9.

Figure 8: The contour line of (UAV demonstrator) [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]

Figure 9: The contour line of (UAV demonstrator) [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]

From Figure 6, it is obvious that [figure omitted; refer to PDF] has many values greater than zero. That is to say, (31) has reasonable solutions. From Figure 7, the maximum value of [figure omitted; refer to PDF] for the albatross' lifting balance allowance in dynamic soaring can be found. For example, if the lift coefficient of albatross ( [figure omitted; refer to PDF] ) is 1 and the gradient of the wind field ( [figure omitted; refer to PDF] ) is 1, in Figure 6, the value of [figure omitted; refer to PDF] is 0.4. Compared with Figure 7, [figure omitted; refer to PDF] m/s; combining with (20), it can be found that [figure omitted; refer to PDF] m/s. That is to say, when albatross' airspeed reaches 17.0 m/s, there exists equilibrium point in this wind field. Just getting the value of airspeed is not enough to make conclusions. The restriction to airspeed is the bank angle of air velocity vector, which is defined as [figure omitted; refer to PDF] . With the geometrical relationship of forces balanced, one can acquire the space curve of [figure omitted; refer to PDF] , which is obtained by projecting the four-dimensional space [figure omitted; refer to PDF] into three-dimensional space. The fundamental is to find out a suitable aerodynamic force to compensate the resultant force of gravitation and inertial force. From (24), (25), and (26), the value range of the bank angle is [30°, 60°]. The albatross can reach the equilibrium position in gradient wind under these limiting conditions. Similarly, the equilibrium points for UAV demonstrator in dynamic soaring can be found in Figures 8 and 9.

4. Conclusions

Dynamic soaring gives the albatross an ability to achieve energy from the environment and prolong its flight endurance, which is a guide to consider using dynamic soaring as a propulsive energy source for small UAVs. The paper presents an equilibrium position theory of dynamic soaring which is the forces balance during the flight. At the equilibrium position, the aircraft stays in a certain speed without kinetic energy loss. When the flight direction is upward, the mechanical energy would increase, which can store energy for flight.

In this paper, the wind relative reference frame of aircraft is noninertial, and the inertial force [figure omitted; refer to PDF] created by wind can be used to do work for aircraft. In the noninertial frame, the mechanical energy variation of the aircraft equals the energy caused by drag and inertial force. When the wind gradient [figure omitted; refer to PDF] and [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] are the airspeed component in the direction of positive [figure omitted; refer to PDF] -axis and [figure omitted; refer to PDF] -axis), the inertial force ( [figure omitted; refer to PDF] ) does positive work to the aircraft. The mechanical energy variation for unit mass is influenced by wind gradient, mass ( [figure omitted; refer to PDF] ), and drag coefficient.

Because of the existence of the gradient wind field, there are some balance points during dynamic soaring no matter how parameters are chosen. In the contour line figures of albatross' [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , it is easy to find out that the maximum value of [figure omitted; refer to PDF] , lift coefficient ( [figure omitted; refer to PDF] ), airspeed ( [figure omitted; refer to PDF] ), and bank angle ( [figure omitted; refer to PDF] ) and the gradient of the wind field are required by the dynamic soaring balance points. By means of performances comparison between UAV demonstrator and albatross, it is manageable to apply dynamic soaring in engineering practice.

Finally, in order to make dynamic soaring profit from further research work as soon as possible, dynamic soaring experiments are necessary. Furthermore, an investigation with theory of balances for 6 DOF aircraft model would be helpful for instructing the dynamic soaring experiments.

Acknowledgment

The authors gratefully acknowledge the financial support from the Fundamental Research Funds for the Central Universities (Project no. CDJZR14325501).

Conflict of Interests

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

References

[1] D. Weatherington, U. Deputy Unmanned Aircraft Systems Roadmap, 2005-2030 , Deputy, UAV Planning Task Force, OUSD (AT&L), 2005.

[2] Y. J. Zhao, Y. C. Qi, "Minimum fuel powered dynamic soaring of unmanned aerial vehicles utilizing wind gradients," Optimal Control Applications and Methods , vol. 25, no. 5, pp. 211-233, 2004.

[3] P. L. Richardson, "How do albatrosses fly around the world without flapping their wings?," Progress in Oceanography , vol. 88, no. 1-4, pp. 46-58, 2011.

[4] G. Sachs, "Minimum shear wind strength required for dynamic soaring of albatrosses," Ibis , vol. 147, no. 1, pp. 1-10, 2005.

[5] G. C. Bower Boundary layer dynamic soaring for autonomous aircraft: design and validation [Ph.D. thesis] , Stanford University, Stanford, Calif, USA, 2011.

[6] J. W. S. Rayleigh, "The soaring of birds," Nature , vol. 27, no. 701, pp. 534-535, 1883.

[7] N. R. J. Lawrance, S. Sukkarieh, "Autonomous exploration of a wind field with a gliding aircraft," Journal of Guidance, Control, and Dynamics , vol. 34, no. 3, pp. 719-733, 2011.

[8] N. R. J. Lawrance Autonomous soaring flight for unmanned aerial vehicles [Ph.D. thesis] , University of Sydney, Sydney, Australia, 2011.

[9] J. L. Grenestedt On Dynamic Soaring , of MAE Departmental Seminars, Princeton University Press, Princeton, NJ, USA, 2010.

[10] M. Deittert, C. A. Toomer, A. Richards, A. Pipe, "Engineless unmanned aerial vehicle propulsion by dynamic soaring," Journal of Guidance, Control, and Dynamics , vol. 32, no. 5, pp. 1446-1457, 2009.

[11] J. C. D. Cone A Mathematical Analysis of the Dynamic Soaring Flight of the Albatross with Ecological Interpretations , Virginia Institute of Marine Science, Gloucester Point, Va, USA, 1964.

[12] L. Peter, "Wind energy extraction by birds and flight vehicles," in Proceedings of the 43rd AIAA Aerospace Sciences Meeting and Exhibit, American Institute of Aeronautics and Astronautics, Reno, Nev, USA, January 2005.

[13] P. P. Sukumar, M. S. Selig, "Dynamic soaring of sailplanes over open fields," Journal of Aircraft , vol. 50, no. 5, pp. 1420-1430, 2013.

[14] B. Ricardo, G. Anouck, A. Mariam, S. Joao, "Shear wind estimation," in Proceedings of the AIAA Guidance, Navigation, and Control Conference, American Institute of Aeronautics and Astronautics, Portland, Ore, USA, August 2011.

[15] J. V. Jose, E. J. Saletan Classical Dynamics: A Contemporary Approach , Cambridge University Press, Cambridge, UK, 1998.

[16] Y. J. Zhao, "Optimal patterns of glider dynamic soaring," Optimal Control Applications and Methods , vol. 25, no. 2, pp. 67-89, 2004.