Characteristics of Ground Effects of an Unmanned Aerial Vehicle with Medium-aspect-ratio Flying Wings during Take-Off

With the continuous development of large and heavy-duty unmanned aerial vehicles (UAVs), take-off and landing technology has become a key link restricting its development. In the analysis of take-off and landing performances of UAVs, ground effect is one of the factors that have to be taken into consideration [1, 2]. Ground effect refers to the nonlinear complex flow field relationships between the aircraft, the atmosphere, and the ground surface during take-off and landing. The large scale of UAVs with medium-aspect-ratio flying wing configurations strengthens the interaction of UAVs with the ground when taking off: this influences the longitudinal aerodynamic characteristics of UAVs. When the take-off height of UAVs with flying wing configurations is different, the ground effects vary correspondingly, so elevator deflection and take-off velocity of UAVs have to be adjusted according to the changes in ground effect. For UAVs with flying wing configurations that are difficult to be controlled in the longitudinal direction, ground effects exert significant influences on manipulation design for safe take-off, therefore, in view of UAVs with medium-aspect-ratio flying wing configurations, it is necessary to study the characteristics of ground effect during take-off and find accurate relationships between longitudinal aerodynamic characteristics with height from the ground, take-off velocity, and angle of attack. Over the past few decades, scholars have studied ground effects by using experimental and numerical simulation methods. By fixing ground boundary conditions and moving ground boundary conditions, the characteristics of ground effects of symmetrical and asymmetrical aerofoil were researched, thus obtaining the influence of height from the ground and angle of attack in take-off on ground characteristics of aerofoil [3-6]. The ground effects of wings with lowand medium-aspect-ratios and a small sweepback angle were investigated [7-11]. The ground effects of largeswept delta wings were also studied under different angles of attack [12-15]. In recent years, the research on ground effects of UAVs has attracted the attention of many scholars. Jones investigated V-shaped tailless UAVs with a sweepback angle of 45° on the leading edge and found that ground effects rise pitching moment and slope of lift curve and reduces induced drag, while the total drag on the UAV increases [16]. Elsewhere In the ground effects of Vand W-shaped UAVs were studied by using the experimental method [17]. Mustachio studied ground effects of unmanned combat aerial vehicles (UCAVs) with a W-shaped configuration and investigated the influences of aileron deflection [18]. Yang researched the impacts of asymmetric ground effects on lateraland heading-directional aerodynamic characteristics of UAVs [19]. Based on the above research, ground effect of UAVs were mainly researched in the experiment by using moving ground to replace the original fixed ground boundary conditions. Although such a research method overcame the deficiency of low accuracy of experimentation when using fixed ground, the experimental costs and complexity were greatly increased. The low-aspect-ratio UAVs were mainly studied in numerical simulation. The take-off process only considered horizontal movement of the ground, while the influences of take-off velocity, take-off angle of attack, and height from the ground were not adequately taken into account. The aim of the study is to research characteristics of ground effects when UAVs with medium-aspect-ratio flying wings take-off to a certain height by using a numerical method. Moreover, the influence of ground effects under different angles of attack, heights from the ground, and take-off velocity on the aerodynamic characteristics of UAVs were investigated, which provides certain aerodynamic data useful for designing control systems governing the take-off process of medium-aspect-ratio UAVs.


Introduction
With the continuous development of large and heavy-duty unmanned aerial vehicles (UAVs), take-off and landing technology has become a key link restricting its development.In the analysis of take-off and landing performances of UAVs, ground effect is one of the factors that have to be taken into consideration [1,2].Ground effect refers to the nonlinear complex flow field relationships between the aircraft, the atmosphere, and the ground surface during take-off and landing.The large scale of UAVs with medium-aspect-ratio flying wing configurations strengthens the interaction of UAVs with the ground when taking off: this influences the longitudinal aerodynamic characteristics of UAVs.When the take-off height of UAVs with flying wing configurations is different, the ground effects vary correspondingly, so elevator deflection and take-off velocity of UAVs have to be adjusted according to the changes in ground effect.For UAVs with flying wing configurations that are difficult to be controlled in the longitudinal direction, ground effects exert significant influences on manipulation design for safe take-off, therefore, in view of UAVs with medium-aspect-ratio flying wing configurations, it is necessary to study the characteristics of ground effect during take-off and find accurate relationships between longitudinal aerodynamic characteristics with height from the ground, take-off velocity, and angle of attack.
Over the past few decades, scholars have studied ground effects by using experimental and numerical simulation methods.By fixing ground boundary conditions and moving ground boundary conditions, the characteristics of ground effects of symmetrical and asymmetrical aerofoil were researched, thus obtaining the influence of height from the ground and angle of attack in take-off on ground characteristics of aerofoil [3][4][5][6].The ground effects of wings with low-and medium-aspect-ratios and a small sweepback angle were investigated [7][8][9][10][11].The ground effects of largeswept delta wings were also studied under different angles of attack [12][13][14][15].
In recent years, the research on ground effects of UAVs has attracted the attention of many scholars.Jones investigated V-shaped tailless UAVs with a sweepback angle of 45° on the leading edge and found that ground effects rise pitching moment and slope of lift curve and reduces induced drag, while the total drag on the UAV increases [16].Elsewhere In the ground effects of V-and W-shaped UAVs were studied by using the experimental method [17].Mustachio studied ground effects of unmanned combat aerial vehicles (UCAVs) with a W-shaped configuration and investigated the influences of aileron deflection [18].Yang researched the impacts of asymmetric ground effects on lateral-and heading-directional aerodynamic characteristics of UAVs [19].
Based on the above research, ground effect of UAVs were mainly researched in the experiment by using moving ground to replace the original fixed ground boundary conditions.Although such a research method overcame the deficiency of low accuracy of experimentation when using fixed ground, the experimental costs and complexity were greatly increased.The low-aspect-ratio UAVs were mainly studied in numerical simulation.The take-off process only considered horizontal movement of the ground, while the influences of take-off velocity, take-off angle of attack, and height from the ground were not adequately taken into account.The aim of the study is to research characteristics of ground effects when UAVs with medium-aspect-ratio flying wings take-off to a certain height by using a numerical method.Moreover, the influence of ground effects under different angles of attack, heights from the ground, and take-off velocity on the aerodynamic characteristics of UAVs were investigated, which provides certain aerodynamic data useful for designing control systems governing the take-off process of medium-aspect-ratio UAVs.

Modelling of ground effects
In the simulation of ground effects in take-off process of an UAV, the calculation model is as shown in Fig. 1.The influences of the ground on flow fields of the UAV need to be taken into account and the direction of free flows velocity should be parallel to the ground in any such calculation.For the above reason, the angle of attack does not show changes in free flows direction, but indicates the angle of the central axis of the UAV to the ground.The simulation of ground effects is based on the principles of relative motion.It is assumed that the UAV is motionless and the velocity of ground movement equals the take-off velocityU,  which is realised by moving-ground boundary conditions.The height h from the ground is defined as the height of the trailing edge of the fuselage above the ground and equals the product of mesh motion Uz and the time since take-off, as realised by way of a dynamic meshing technique.
Fig. 1 The calculation model of ground effect

Dynamic grid technology
In the analysis of ground effects, with the increase of height from the ground of the UAV, the grid between the ground and the UAV extends from the initial small distance to a few meters or more.Therefore, the generation of the dynamic mesh is a key problem.If a dynamic mesh has poor robustness, it may result in negative grid elements, thus causing a failure during calculation.In the meantime, in the simulation of ground effects, the scale of grids near the ground cannot too be large; otherwise the flows in ground boundary layers and disturbance of ground boundary layers to flow fields of the UAV cannot be simulated accurately.In the study, structural grids and a dynamic layering technique were used to update the dynamic mesh used in the simulation.The core idea is to add dynamic layers according to the height of ground grid layers adjacent to the boundary to the motion (Fig. 2).

Fig. 2 The mesh model of dynamic layers
In Fig. 2, the grids in the jth layer are close to the motion boundary-ground, and 0 h indicates the height of grids in the jth layer.When the UAV took off, the height from the ground gradually increased and the grids in the jth layer were initially stretched.When the height of the grids increased to the critical value given by Formula (1), the grid elements in the jth layer were divided by using the constant ratio splitting method.
where: h' and s a indicate the height of ideal grid elements and the splitting factor of each layer.

Governing equation
The governing equation used is the unsteady Navier-Stokes equation and is expressed, in inertial Cartesian coordinates, as follows [20]: where: Q indicates the flow field variable, G , V F , and b q denote the convective term, the dissipation term, and the grid velocity, respectively.The equation is not in closedform.By introducing the thermodynamic relationship for a perfect gas and Stokes assumption, the above equation set forms closed-form equations.In solving the governing equations, the Spalart-Allmaras turbulence model was used and Roe flux-difference splitting scheme was used in space discretization [21].Moreover, the fully implicit dual time stepping method was employed in time-marching, that is, implicit marching and explicit Runge-Kutta marching methods were used in real-and pseudo-time [22].

Verification
In order to verify the validity of the numerical method used in this paper, the wing in the literature [5] was calculated.The cross-sectional airfoil of the wing is NACA4412, and the chord length of the airfoil was 150 mm and the span of the wing was 600 mm, as shown in Fig. 3.   、8 , along with moving ground boundary conditions.The calculation results of ground effects under different angles of attack and heights from the ground were compared with the test results in the literature [5], and comparisons of the results are shown in Figs. 4 and 5.
It can be seen from the above figures that the calculation results coincide with the test results, which indicates that the analysis method for ground effects proposed in the study is reliable.

Results and discussions
A UAV with W-shaped flying wing configurations with zero sweep and torsion angle was adopted for analysis, as shown in Fig. 6.Moreover, it is found that the aspect ratio, span, and chord length of the wing root were  m, respectively.The cross-section of the UAV wing is a thick airfoil with a curved edge at the trailing edge.Owing to it being the lift, drag, and longitudinal characteristics that were studied, a half-model of the UAV was analysed to reduce the computational burden.The coordinate origin was at the head of the UAV and the calculation was carried out in the domain which was 10 times the chord length of the wing root in the X-direction and five times the half-length in the Y-and Z-directions.The structural grid, with its H-H topological structure, was used for grid division and the grids near wall surface and the ground were refined.The number of grid elements was As shown in the Fig. 7, with the increase in takeoff velocity, the suction area and amplitude of the lower surface of the UAV gradually decreased and the negative pressure on the upper surface decreased slightly.

Influence of height from the ground
When the angle of attack was set to ,3 4  the characteristics of the ground effects of the UAV at different heights from the ground were obtained and the calculation results are shown in Fig. 8.A suction area was formed between the lower surface of the UAV and the ground and the suction force reached its peak value at h/cA=0.067.When the height from the ground increased to h/cA=0.153, the effects of air flow between the ground and the surface of the UAV were weakened and the peak suction force on the lower surface decreased.When the height from the ground reached h/cA=0.239, the suction force on the lower surface of the UAV was quasi-stable and remained unchanged when the height from the ground increased further. It can be seen from the change curve of drag coefficient in Fig. 9, a that, with increased height from the ground, the effects of air flow between the ground and the UAV decreased.Therefore, the limitations of the ground on the down wash of wings of the UAV decreased, increasing the induced drag of the UAV, thus increasing its drag coefficient.When the height from the ground reached h/cA=0.5, ground effects exerted only small influences on the aerodynamic characteristics of the UAV and drag coefficient remained unchanged even though the height from the ground increased further.As demonstrated in Fig. 9, b, as the height from the ground decreased, the lift coefficient diminished to some small value due to the influences of the suction area developed on the lower surface of the wings.In particular, when the angle of attack was 2°, the suction amplitude on the lower surface of the wings was large and the lift coefficient was found to be a large negative value.With increasing height from the ground, the suction amplitude on the lower surface gradually decreased, while the lift coefficient gradually increased.When the height from the ground reached h/cA=0.5, ground effects only slightly affected the aerodynamic characteristics of the UAV and the lift coefficient remained unchanged thereafter.As shown in Fig. 9, c, at low altitudes, the pitching stability of the UAV decreased due to the influence of the suction area generated on the lower surface of the wings.

Influence of take-off angle of attack
By using the take-off velocity of =68 U  m/s, the characteristics of ground effects of the UAV under different angles of attack from the ground were obtained and the results are as shown in Fig. 10.
As demonstrated in Fig. 10, when the angle of attack was 2°, a small gap and a large flow velocity were shown between the lower surface and the ground, resulting in a large suction force on the lower surface.When the angle of attack gradually increased, the flow velocity between the lower surface and the ground gradually decreased and the peak suction force on the lower surface decreased.

Conclusions
The characteristics of ground effects on a UAV with medium-aspect-ratio flying wings in the take-off process were investigated from the perspective of the longitudinal aerodynamic characteristics.To study the influences of take-0ff velocity, the height from the ground, and the angle of attack, numerical simulations were conducted under different conditions.The conclusion is as follows: 1.The height from the ground exerted significant influences on the longitudinal aerodynamic characteristics of a medium-aspect-ratio UAV.When the height from the ground was small, a large suction area was formed on the lower surface of the wings.The suction amplitude increased by 37% with a decrease in the angle of attack, which resulted in a negative lift coefficient when the angle of attack was small.As the height of the ground increases, the suction area and amplitude of the lower surface of the wing decrease, and the lift coefficient increases.After the height from the ground reached h/cA=0.5, ground effects exerted little influences on the UAV and the lift coefficient remained unchanged.
2. The angle of attack exerted a certain influence on the longitudinal characteristics of the UAV.When the angle of attack was small, a trend evincing a wide range of down wash was found on the lower surface of the trailing edge of the wings, reducing the lift coefficient.With increase of angle of attack, this trend weakened, thus increasing the lift coefficient.The changes in take-off velocity did not significantly influence the longitudinal aerodynamic characteristics of the UAV.With increasing take-off velocity, the pressures on the upper and lower surfaces of the UAV both reduced, but only slightly.3.In this research on a UAV with medium-aspectratio flying wing configuration, due to the strong suction force generated on the lower surface of the wings, the lift coefficient and the pitching stability decreased with the height of the ground decreased.This differed from the ground effects seen in ground-effect vehicles.In fact, the magnitude of the suction force acting on the lower surface of the wings was affected by the geometric shape of aerofoil to a significant extent, therefore, not all medium-aspect-ratio UAVs have a decreasing lift coefficient at lower height of the ground.

CHARACTERISTICS OF GROUND EFFECTS OF AN UNMANNED AERIAL VEHICLE WITH HIGH-ASPCET-RATIO FLYING WINGS DURING TAKE-OFF S u m m a r y
In view of the take-off process of an unmanned aerial vehicle (UAV) with medium-aspect-ratio flying wing configurations, the influences of ground effects on longitudinal aerodynamic characteristics of the UAV were investigated.By using a numerical simulation method, the effects of different angles of attack, heights from the ground, and take-off velocity on the UAV were studied.Moreover, based on comparisons with test data in previous literature, the numerical calculation method was verified to be accurate.The height from the ground is the primary factor influencing longitudinal aerodynamic characteristics of the UAV.With the increase in height from the ground, the lift and drag coefficients of the UAV, and its pitching stability, increased.When the angle of attack was small, a large suction area was generated on the lower surface of wings with a large suction amplitude therewith.This results in a large range of downwash.As the angle of attack increased, the lower surface of the trailing edge for wings presents weakening downwash.As the take-off velocity rose, the suction area on the lower surface of the UAV, and its amplitude, gradually decreased, and the negative pressure on the upper surface also decreased to a certain extent.The research results can provide a certain theoretical basis for the safe design of take-off operations for medium-aspect-ratio UAVs.

8
Pressure coefficient distribution of model spanwise section at y=2 m: