Prerpints.org logo
Preprint
Article

Experimental Characterization of Propeller Induced Flow (PIF) below a Multi-Rotor UAV

Altmetrics

Downloads

135

Views

43

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

09 January 2024

Posted:

10 January 2024

You are already at the latest version

Alerts
Abstract
The availability of multi-rotor UAVs with lifting capacities of several kilograms, allows for a new paradigm in atmospheric measurement techniques, i.e. the integration of research-grade sonic anemometers for airborne turbulence measurements. With their ability to hover and move very slowly, this approach yields an unrevealed flexibility compared to mast-based sonic anemometers, for a wide range of boundary layer investigations that require an accurate characterization of the turbulent flow. For an optimized sensor placement, potential disturbances by the propeller-induced flow (PIF), as well as the potential impact of the sensor weight and angular momentum on the flight performance in the case of a boom-mounted sensor, have to be considered.
Keywords: 
Subject: Environmental and Earth Sciences  -   Atmospheric Science and Meteorology

1. Introduction

The proper characterization of atmospheric turbulence is a critical factor for the understanding of the structure and dynamics of the atmospheric boundary layer (ABL) [1,2]. Accurate wind and temperature measurements with high spatial and temporal resolution are thus crucial for a wide range of scientific applications in basic and applied ABL research. Mast-, tower- and bridge-based measurements with ultrasonic anemometers [3,4,5,6] have over more than five decades developed into the golden standard for turbulence measurements in experimental ABL research [7].
Ultrasonic anemometers sample the three-dimensional flow field and the sonic temperature with a high temporal resolution of typically 10 to 100 Hz. The latest generation of research-grade ultrasonic anemometers from various providers has a rated accuracy in the order of 0.05   m / s and 0.10   m / s for the vertical and horizontal wind velocity components, respectively. The combined accuracy and robustness (no moving parts) of these instruments establish them as state-of-the-art sensors for high-resolution in-situ observation of turbulent velocity and temperature fluctuations. However, masts or towers as sensor carriers for sonic anemometers limit considerably the measurement flexibility. Recent studies in basic ABL meteorology [8,9] and wind energy meteorology [10,11] highlight the need for an enhanced comprehension of key ABL processes. This necessitates the advancement of our measurement techniques beyond the traditional mast-based approach.
Lidar remote sensing offers one pathway in this direction. While scanning Doppler wind lidars, such as short-range or long-range WindScanner systems, offer valuable wind measurements [12,13], their flexibility and applicability are constrained by several factors. Unlike ultrasonic anemometers, they do not capture the virtual temperature, which is a crucial parameter for comprehensive ABL studies. Their high cost and lower effective sampling frequency limit widespread usage. In addition, the inherent spatial averaging over the probe volume, typically in the order of 20 m for pulsed lidar systems, hampers the analysis of higher-frequency turbulence characteristics.
A first step towards increased measurement flexibility using sonic anemometers was their deployment with the help of tethered balloons [14,15,16] and kites or blimps [17,18]. These studies have demonstrated that these systems can provide reliable turbulence data, when the sensor is mounted correctly, i.e. far enough from the carrier platform to avoid flow distortion, and when the sensor’s motion is recorded and corrected. Lifting a sonic anemometer with a battery power supply for an appropriate measurement time of at least 30 minutes demands kites, balloons, or blimps of considerable size. The deployment of such systems brings additional infrastructural and logistical requirements, concerning, e.g., winch systems and gas supply, limiting the flexibility of deployment. For many of the tethered systems, there is also an upper operational limit of wind speed in the order of 10 m / s .
As a consequence of the rapid development in the field of uncrewed aerial vehicles (UAVs) over the last two decades, drones have also found their way as flexible, mobile, and cost-efficient sensor carriers in atmospheric research [19,20,21]. The commercial availability of corresponding airframes with sufficient payload capacities and the accessibility of freely programmable open-source autopilot solutions make UAVs now also well-suited as sensor platforms for atmospheric turbulence measurements. Turbulence measurements on fixed-wing systems, with typical cruising speeds of 15 m / s –25 m / s , usually rely on multi-hole probes [22,23,24,25,26,27,28] and require complex correction and compensation algorithms for the attitude and, in particular, the relatively high horizontal speed of the aircraft [29,30,31,32]. With typical flight times ranging from 30 minutes to several hours, those systems can measure turbulence along the flight path over larger areas. For applications that require stationary measurements, e.g., for the determination of coherence of turbulence for structural design [13,33], for missions over highly heterogeneous surfaces, or for atmospheric profiling with high spatial resolution, i.e., slow ascent rates over a fixed point, e.g., for the investigation of the stable ABL [8,9], rotary-wing UAVs are the obvious choice. The reduced endurance compared to fixed-wing systems, in the order of half an hour to an hour, is for such applications compensated by the ability to hover or move very slowly. This makes rotary-wing UAVs also suitable for operating very accurately close to the ground or in the vicinity of buildings or other structures, such as wind turbines.
Multi-rotor drones bear consequently a large potential as suitable sensor-carrier for sonic anemometers, and corresponding approaches are reported in the literature [34,35,36,37,38,39]. So far, those approaches focus primarily on the measurement of the mean horizontal wind speed, often carrying miniaturized sonic anemometers with measurement geometries not fully capable of providing reliable measurements of the vertical wind component. The few studies applying research-grade sonic anemometers [35,36,37] on the UAVs have not yet proven the ability to measure the full spectrum of undisturbed ambient turbulence.
To reach that goal it is required to place the sensor well outside the propeller-induced flow (PIF). This can be realized by either mounting the sonic anemometer on a sufficiently long extension arm beside or above the drone or by flying it as a sling load far below the UAV. Both strategies require information on the PIF created by the drone in operation to select positions with undisturbed conditions or at least minimized flow distortion of an acceptable level. In general, the PIF decreases with distance from the propellers, placing the wind sensor far away from the rotors is thus a simple and effective strategy to mitigate the PIF influence on the wind measurements [36,40]. A rigidly attached mass, positioned away from an airframe’s center of gravity, e.g., by a fixed boom, introduces, however, angular momentum and inertia that complicate in-flight stabilization and negatively impact flight dynamics. Locating optimal positions near the drone’s center of mass that minimize flow distortion from the PIF can enhance the design of drone-based systems for accurate ambient turbulence measurement with sonic anemometers.
Computational fluid dynamics (CFD) studies offer a practical alternative to wind tunnel tests [41] for in-flight measurements of drones capable of carrying research-grade sonic anemometers, which would otherwise necessitate very large wind tunnel facilities. Corresponding CFD simulations are well established in the field of UAVs [42], and their application spans from the investigation of propeller efficiency, performances, and workloads [43,44] to the characterization of the PIF. The latter studies investigate, however, the PIF features mainly concerning its effect on flight stability, and are thus focusing on the near-field flow around the drone, rather than for distances relevant for ultrasonic sensor placement [45,46,47,48]. To the authors’ knowledge, the first attempt to use CFD for sensor placement considerations on a large multi-rotor drone has been described in Ghirardelli et al. [49]. The simulations are performed for an airframe of the size and properties identical to the one used in the experimental study presented hereinafter. In this study, we carry out a first evaluation of the corresponding numerical simulations. For this, we have designed and performed a low-cost experiment to directly measure the PIF under controlled indoor conditions.
The manuscript is organized as follows: Section 2 details the experimental setup, covering the selected UAV, its positioning, and the anemometer-equipped measurement-rack design. Section 3 outlines the measurement strategy, the experiment execution, and the data processing techniques. Section 4 discusses the PIF measurement results, and compares them with CFD simulations and environmental visualizations. Section 5 presents a conclusion and outlook.

2. Experimental Setup

2.1. UAV Description

The UAV chosen as a potential sensor carrier for a sonic anemometer is the Foxtech D130 X8 (Figure 1), a multi-copter with coaxial contra-rotating 28-inch (71 cm) propellers, arranged in four pairs that share the same axis of rotation. The propellers in each pair spin in a contra-rotating set-up driven by eight brushless electric motors (T-MOTOR U10II KV100). In its default flying configuration it is powered by two 6S LiPo batteries. The weight of the frame including the motors is approximately 9 k g . Together with the batteries, the system has, without scientific payload, a take-off weight (TOW) of about 15 k g and a maximum endurance in hovering mode of 45 min , depending on the atmospheric conditions. With a maximum TOW of 36 k g the system leaves a large margin for either additional scientific payload or added battery capacity to extend the flight time. More information on the UAV’s specifications can be found in Table 1, and a more detailed description of the system is given in Ghirardelli et al. [49].

2.2. Drone Mounting

Outdoor spaces offer enough room for the airflow from drone propellers to spread. However, these areas also bring unpredictable ambient air movement, making it challenging to identify the specific effects of the propellers on the airflow. Experimenting indoors, within a confined and largely controlled environment, also raises certain challenges. Indoor environments are subject to secondary air circulations created by the interaction of the generated flow within the boundaries of the limited indoor space (i.e., floor, walls, and ceiling), potentially leading to largely disturbed measurements of the targeted PIF. Despite these potential challenges we opted for an indoor setting, focusing on a design that minimizes the potential flow disturbances discussed above.
We found that the drone’s size compared to the hall’s ceiling height was too large for indoor hovering without significant interference between the PIF and the floor. By positioning the drone at a 90-degree angle to its normal flying direction (as shown in Figure 2), which aligns the downwash with the ground, we anticipated a notable reduction in disturbances caused by ground effects.
The mounting rack, sketched in Figure 2, consisted of a frame housing the ground control station (GCS), two power supplies units (PSU), and a horizontal extension arm holding the drone, placing its center 2.3   m away from the base. The rack was attached to a forklift, enabling both vertical and horizontal positioning upfront the experiment. The UAV was oriented centrally towards the large bay door in order to direct the downwash out of the hall (see the right panel in Figure 3 and the schematic view from above in Figure 5). Opening the door during the experiments allowed the far downwash to exit the enclosed hall volume, while the primary area of interest for our PIF measurements, a few rotor diameters beneath the drone, remained indoors, well shielded from external factors.
The fixed mounting of the drone during the experiments had additional benefits, including absolute stability for precise flow sensor positioning, continuous power from two SkyRC Technology eFUEL power supplies with a maximum output of 50 A at 25 V DC, and wired connection from a Panasonic TOUGHBOOK running Mission Planner GCS software to the ArduPilot-based flight controller to regulate the throttle settings for the motors.

2.3. Sonic Anemometer Measurement Rig

Five Campbell Scientific CSAT3 ultrasonic anemometers were used to measure the PIF. The anemometers were mounted to a mast with vertical spacing increments of 0.50   m , resulting in a measurement range from 0.60 m –2.70 m above the ground. The mast was attached to the undercarriage of a height-adjustable table, allowing for an additional 0.40   m of vertical adjustment, extending the maximum measuring height to 3.10   m above ground. The whole measurement rig was then attached to a Euro-pallet and could thus be moved manually with a hydraulic jacklift (see also right panel in Figure 4). Each of the sonic anemometers sampled the 3D wind vector and the sonic temperature at a frequency of 10 Hz , and was connected to a Campbell Scientific CR3000 data logger for data storage and synchronization.
To associate velocity measurements to spatial points, distance measurements to a reference surface, located at the respective wall of the hall perpendicular to the movement of the rig, were taken using the Benewake TF02 LIDAR range finder. Those data were also recorded by the Campbell Scientific CR3000 data logger and thus exactly synchronized with the sonic anemometer data. The height of the sonic anemometers was noted down manually for each series of 2 D scans (see Sect. 4) and for redundancy also measured by two barometers (one at a fixed level and one attached to the height adjustable measurement mast). Both the LIDAR range finder and the two barometers were logged at a sampling frequency of 1 Hz . The sensor assembly was placed on three Euro-pallets for mobility and the desired starting height. A schematic of the sensor setup is shown in Figure 4. Real-time data access and visualization were achieved using a Panasonic TOUGHBOOK running LoggerNet software and connected to the CR3000.

3. PIF Measurements

3.1. Measurement Pattern

Figure 5 depicts the experimental coordinate system and its location in the hall. The z-axis is chosen normal to the floor, increasing with height above the hall surface. The downwash propagates along the y-axis in the positive direction, whereas the x-axis is parallel to the rotor plane with the positive direction pointing away from Wall A. The measurement planes, designed to intersect the downwash across the x-z plane, were defined at intervals equivalent to the drone’s rotor diameter, D, of 71 cm. The first cross-section was placed at 2 D downstream from the rotor plane, the last one at 10 D. Additionally, measurement planes tangential to the downwash, spanning the y-z plane, were defined along the centerline of the downwash, i.e., x / D = 0 , and at x / D = ( 1 , 2 , 3 ) . All measurement planes cover the height interval between 0.60 m to 3.10 m above ground. The solid black lines in Figure 5 represent the chosen measurement pattern, each line corresponding to one x-z or y-z cross-section. At each cross section, the array of anemometers was manually moved back and forth with a velocity in the order of 0.05  m s−1 using a jacklift. For robust statistics, this measurement pattern was repeated 3 times for a given altitude setting of the height-adjustable table, thus sampling every point along the height of the 5 sonic anemometers 6 times. After that, the height-adjustable table was raised by 0.10   m , and the procedure was repeated until the whole cross-section was sampled in a vertical resolution of 0.10   m . Assuming quasi-stationary flow conditions, this measurement approach is expected to provide detailed insights into the average speed and structure of the drone downwash.
The entire measurement pattern was carried out for two different throttle settings of 35% and 45%, corresponding to a TOW of roughly 15 kg and 20 k g .

3.2. Data Processing

Upon the completion of data collection, two data sets were gathered. Data set 1 (DS1) contains all the velocity vectors u x , u y , and u z collected by the five sonic anemometers with a sampling frequency of 10 Hz . The vectors u x , u y , and u z describe the measured wind speed relative to the coordinate system used for the experiment, where u x refers to the wind speed along the x-axis, u y along the y-axis, and u z perpendicular to the ground. The second data set (DS2), important for the exact localization of the flow measurements, contains the distance measurement to the relevant reference point and was gathered at 1 Hz . The data was recorded over two days, the 24th and 25th of October 2022. In total, about 6 hours of measurements resulted in 219 760 relevant wind velocity vectors in 10 Hz resolution in DS1 for each of the 5 anemometers, corresponding to 21 976 spatial data points in DS2. To ensure the viability of the collected data points the resulting data sets were first cleaned for missing data and physical outliers. The lidar distance measurements in DS2 contained several unphysical jumps and spikes. Considering that the data points were collected in a continuous and monotonous movement of a few m s−1, all consecutive data points moving more than 0.3   m were dismissed and replaced by linearly interpolated values. For sequences of data points where the neighbors were also not deemed viable, a linear motion profile was fitted, based on the start and end time and position of that corresponding pass. In DS2, 12 % of data points were initially recorded as ’not a number’ (NaN). The distribution of these measurements throughout the time series was uniform, displaying no discernible pattern in their occurrences. The cause behind this substantial number of faulty measurements remains unknown and can only be subject to speculation. However, after applying linear interpolation, the occurrence of NaN measurements decreased to zero.
Following the correction of the spatial data set, DS1 and DS2 were merged based on their recorded common timestamps. DS2 was linearly interpolated from the 1 Hz original data onto the 10 Hz resolution of DS1, resulting in a positioning of each wind measurement. The wind velocity vectors parallel to the motion of the anemometer assembly were adjusted by subtracting the calculated velocity of the assembly movement. For further analysis and interpretation, the flow measurements were now binned over 0.10   m spatial intervals. This results in two multi-dimensional arrays: one for the y-z cross section and another for the x-z cross-section. Each two-dimensional array in the x-z plane spanned 3.5   m in the y direction and 2.5   m in the z-direction, while in the arrays in the y-z plane spanned 6 m in the x-direction and 3.5   m in the z direction (Figure 6). Each matrix within these arrays represents a measurement plane, with each cell storing all velocity vectors captured over the corresponding 0.10   m spatial interval (Figure 6). On average, each cell contained about 35 velocity vectors. Velocity data outliers in each cell were eliminated using statistical measures, employing a 2.5 median absolute deviation cutoff threshold individually for each spatial cell.
In the next step, the mean velocity and mean velocity variance were calculated for each cell. The cell-averaged mean PIF was estimated as the magnitude of the mean velocity vector of that cell, using
P I F i , j , k ¯ = u i , j , k ¯ 2 + v i , j , k ¯ 2 + w i , j , k ¯ 2 .
The PIF variability of each cell was estimated using turbulent kinetic energy (TKE) per unit mass calculated using the following equation:
e i , j , k = T K E i , j , k m i , j , k = 1 2 ( u i , j , k ′ 2 ¯ + v i , j , k ′ 2 ¯ + w i , j , k ′ 2 ¯ ) ,
applying cell-wise Reynolds decomposition,
χ i , j , k ′ = χ i , j , k − χ i , j , k ¯ ,
with χ ∈ u , v , w .

4. Results and Discussion

4.1. Mean PIF

For the x-z measurement planes, the perspective of the plotted results is looking from the UAV along the y-axis downstream through the relevant plane, as indicated by ”Perspective A” in Figure 5. Conversely, for the y-z measurement planes, the perspective is looking at the downwash along the x-axis from a location of positive x-value, as indicated by ”Perspective B” in Figure 5, i.e., representing cross sections of the main downwash at different distances from the center line of the drone.
Visualizing the matrices within these tensors as heat maps, we get a good visual representation of the PIF produced by the drone as it propagates downstream. In these visualizations, the PIF corresponds to the cell-averaged velocity magnitude (Eq. 1), where brighter colors denote higher values. The depicted results, showcased in Figure 7 and Figure 8, illustrate the progression of the PIF for a throttle setting of 35%. Figure 7 depicts the mean PIF measured across the x-z planes, employing ”Perspective A” in Figure 5. These planes are set at increments of one rotor diameter, beginning at a distance of two rotor diameters downstream of the mean rotor plane, y = 2 D , and extending to y = 9 D . Figure 8 depicts the y-z planes, parallel to the main downwash direction seen from ”Perspective B” in Figure 5. The first plane is at the center line (CL, at x = 0 ) of the UAV, and subsequent planes are measured at x = 1 D and x = 2 D from the center. Each colored cell in both figures corresponds to a 10 by 10 cm cell within their respective matrices.
The development of the PIF can be characterized as follows:
  • y = 2 D , 3 D : Distinct hotspots of high PIF are clearly visible, one for each rotor set. The peak velocities within the range of 6 m s−1 to 7 m s−1 are observed at the center of the four hotspots. Beyond a distance equivalent to 1.5 D rotor diameter in x- or z-direction from the drone’s center, there is no substantial observed flow disturbance.
  • y = 4 D , 5 D : The initially distinct hotspots begin to converge, forming a more consolidated and uniform downwash structure. The region with the highest velocities is centered behind the drone, with peak velocities of 5 m s−1 to 6 m s−1 for a throttle setting of 35%. The data indicate a more pronounced onset of downwash expansion at this stage. The figures also show that the air between the downwash and the floor, starts to speed up.
  • y = 6 D , 7 D : The hotspots corresponding to the four-rotor sets have now completely dissipated. The downwash center, marked by the region of greatest velocities, remains largely centered relative to the drone. Notably, the downwash starts to expand asymmetrically, indicating the potential effect of wall, floor and ceiling in the far-flow of the downwash.
  • y = 8 D , 9 D : At this stage, the downwash center has shifted towards the left by about 0.5 D . Peak velocities around 4 m s−1 were observed. The edges of the downwash lean more towards the bottom left.
It is worth noting here that both the acceleration of the air under the downwash at a y-position of 4 D to 5 D and the predominant expansion of the downwash in the negative x-direction, can likely be attributed to the specifics of the test setup and the available space. The asymmetric expansion of the downwash in the z-direction is probably a consequence of the drone not being mounted at a sufficient height, leading to air acceleration between the ground and the downwash. The expansion of the downwash in the negative y-direction is most likely a result of inflow through the bay door, driven by the need to maintain mass equilibrium as the drone displaces air out of the confined space through the bay door. Qualitative measurements of this compensation flow were conducted, indicating its presence, but not exactly quantifying its strength and extension.

4.2. PIF Variabillity

Observing the trends discerned from the results thus far, and keeping in mind that the primary area of interest for sensor placement of a boom-mounted wind sensor is within the first few D downstream of the rotor plane, the PIF effect on the flow will be analyzed by TKE cross-sections up to a downstream distance of y = 5 D .
Preprints 95947 g009
Figure 9. Total TKE across x-z planes at different distances Y below the drone at 35% throttle setting.
Figure 9. Total TKE across x-z planes at different distances Y below the drone at 35% throttle setting.
Preprints 95947 g009
The results demonstrate that the region with the most significant variability in PIF lies at the interface where the faster-moving downwash meets the stagnant surrounding air, i.e., the area with the strongest velocity shear. It is also important to highlight that the variations in PIF are predominantly concentrated within the downwash region. This underpins the idea that a substantial portion of PIF disturbance for a boom-mounted sonic anemometer can be effectively mitigated by siting the sensor below the rotor plane at a distance of 1 D to 2 D away from the drone’s center.

4.3. Comparison with CFD Simulations and Environmental Observations

One main motivation for the low-cost experiment presented in this study was to provide a first evaluation for the realism of CFD simulations for the modeling of PIF of our chosen multicopter setup. Corresponding simulations for this system, with a focus on wind sensor placement, have recently been performed by [49].

4.3.1. CFD Simulations

The simulations presented and discussed in the following have been performed by Ansys Fluent. It has to be noted that the actual geometry of the drone is simplified, e.g., by assuming an actuator disc instead of rotating propellers [50,51,52], and neglecting the body of the drone. Thus, the drone is represented by eight two-dimensional discs that create an instantaneous pressure jump in the flow. In the simulation, the domain is a cube with each side measuring 20 D. Within this cubic domain, the drone is modeled with its rotors plane oriented horizontally, ensuring that the downwash is directed vertically downward. The drone is positioned precisely in the horizontal center of the cube, equidistant at 10 D from each side wall. Vertically, it is placed slightly towards the top, 7 D from the ceiling and 13 D from the bottom, allowing for ample space for airflow and downwash development. To enhance the stability and convergence of the solution, the top surface of the cube is defined with an inlet velocity of 0.8  m s−1. Conversely, the opposite face at the bottom is designated as a pressure outlet. Results are rotated to align with the experimental measurements’ reference frame, and PIF values are refined by subtracting the background flow.
Details on the CFD setup, including the meshing procedure and parameters can be found in [49] and [53].
Figure 10 shows the resulting vertical velocities in the downwash modelled at different distances below the drone. The distances presented are chosen in accordance with the upper four panels in Figure 7. The modeled vertical velocities are, in general, higher (max. values of around 9 m s−1 for a distance of 2 D) than our measurements in the hall (max. around 7 m s−1). This might partially be explained by the simplified representation of the drone in the model simulations, in particular the intrinsic uncertainty in relating the chosen throttle setting for the drone in the experiment to an exact value of the instantaneous pressure jump of the actuator discs in the CFD simulations. In addition is the flow simulated without background turbulence, a feature that is typically present in natural flows. Despite this difference in the flow magnitude, compares the structure and dimension of the downwash between model and measurement very well. Both identify four clearly separated downwash areas at distances of 2 D and 3 D. The slightly elliptical shape is caused by the tilt angle of 8° for the propellers of our Foxtech drone, which is also implemented as a corresponding tilt of the actuator discs in the CFD. For distances beyond 4 D, the individual downwash centers have clearly merged into a single, large, and rather homogenous downwash zone.
The general overestimation in the downwash velocity described above can also be seen in the comparison of the cross sections along the y-z plane along the centre line of the drone (upper panels in Figure 8 for the measurements and Figure 11 for the CFD simulations). The consecutive cross section at 1 D distance shows a fast decrease of the downwash effect, that has nearly fully disappeared at 2 D. The width of the downwash area and its structure are rather similar. Some minor differences observed in the experiment can be associated to the interaction of the downwash with the floor of the hall at the lower boundary, and the observed counterflow entering the open bay door close to the ceiling, inhibiting the widening of the downwash at the top.

4.3.2. Outdoor Observations

Another possibility of a qualitative comparison and validation of our measurement results arises from a visualization of the PIF of a hovering multi-rotor drone over a calm water surface in low-wind conditions. Figure 12 shows three consecutive still pictures extracted from a video of a slowly descending drone above a lake. At heights above 15 D, the downwash at the surface is not strong enough to create any visual disturbance (picture to the left). Closer to the surface, visualized in the picture in the center for a height of 7 D, a clear single central downwash appears, surrounded by a ring of small surface waves indicating the flow divergence in the outflow region. At a level of around 4 D, the single downwash quickly transforms into four separated downwash zones created by the individual propellers (picture to the right). Although created by rather different drone systems (the Foxtech D130 octocopter in X8 configuration for the laboratory tests, and the DJI Matrice 300 RTK quadcopter hovering over the lake), the resulting flows show quite intriguing similarities that indicate a considerable degree of generality in the PIF created by a drone with four centers of rotation. The most prominent one is the transformation of the individual rotor flows to one combined downwash flow area, occurring in a distance of approximately 3 D to 5 D below the rotor plane, in good agreement with both our measurements and the CFD simulations.

5. Conclusions

Given the need to understand how drone systems affect the surrounding air to ensure optimal sensor placement, and the lack of an adequate wind tunnel facility at hand, this paper has shown that there are ways to obtain the desired measurements and information on a shoestring budget. The described approach of tilting the drone by 90 degree, and measuring the now horizontal downwash with a rig of five sonic anemometers in cross sections in various distances perpendicular and parallel to the downwash direction, has shown that it is possible to obtain reliable and realistic measurements of the PIF generated by the drone. Some distortions in the results were observed at locations further away from the rotor plane, which can be attributed to the specific constraints of the available experimental indoor space. Despite those unavoidable constraints and the implied potential interactions with ceiling, wall and floors, we are confident that the approach of conducting the experiment in a controlled environment poses more advantages than backdraws.
The data obtained from this experiment corroborated the results obtained from precursory CFD simulations Ghirardelli et al. [49]. The characteristics of the experimental results matched well with the characteristics of the simulations, even though the conditions of the two environments differed slightly. An observed overestimation of the CFD modelled downwash velocities in the order of 20% can be mainly attributed to the uncertain relationship between the throttle setting of the drone and the chosen pressure jump prescribed at the location of the actuator discs in the CFD simulations. With respect to our sensor placement considerations, this overestimation indicates that following the CFD simulations will provide a rather conservative approach of the flow disturbances created by the PIF.
The retrieved dataset can serve as an expansion of the existing foundation to be used for future analysis and comparison with other model simulations, e.g., in combination with a complementary experimental dataset of high-resolution, in-flight measurements of the drone downwash of the Foxtech D130 X8 drone by a short-range lidar WindScanner system described by Jin et al. [53].
The combination of CFD simulations with targeted measurements provides in our opinion a reliable and cost-efficient framework to address the topic of optimal sensor placement on multi-rotor drones. Simulations can be efficiently performed for a large number of drone configurations and environmental flight conditions, thus identifying potential sweet spots for the sensor mounting. Those areas have then to be investigated in more detail by corresponding high-quality measurements.

Author Contributions

Conceptualization, A.F., M.G., S.K. and J.R.; methodology, A.F., M.G., S.K., T.K., and J.R.; software, A.F., M.G., S.K., and E.C.; formal analysis, A.F, M.G. and S.K.; resources, M.G., S.K., and J.R.; data curation, A.F., M.G., and J.R.; writing—original draft preparation, A.F., M.G., S.K., E.C., and J.R.; writing—review and editing, A.A., M.G., S.K., E.C., T.K., and J.R.; visualization, A.F., M.G., S.K., E.C., and J.R.; supervision, M.G., S.K., and J.R.; project administration, J.R. and S.K.; funding acquisition, J.R. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

Parts of this research was funded by the European Union Horizon 2020 research and innovation program as part of the Marie Skłodowska-Curie Innovation Training Network Train2Wind (grant no. 861291; https://www.train2wind.eu/, last access: 13 October 2023).

Data Availability Statement

The data are available from the corresponding author upon request.

Acknowledgments

Many thanks to our institute engineers Anak Bhandari, Helge Bryhni, and Tor de Lange for their help and assistance in preparing the measurement setup and accepting the several day blockage of our storage hall at Marinoholmen for the experiment. The authors are also very greateful to Alizee Lehoux from Uppsala University for the provision of the video material of a hovering drone that formed the basis for Fig. 12 in the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
ABL atmospheric boundary layer
CFD computational fluid dynamic
GCS ground control station
PIF propeller induced flow
PSU power supply unit
RPAS remotely piloted aircraft system
TKE turbulent kinetic energy
TOW take-off weight
UAS uncrewed aerial system
UAV uncrewed aerial vehicle

References

  1. Stull, R.B. An Introduction to Boundary Layer Meteorology; Springer Netherlands: Dordrecht, 1988. [Google Scholar]
  2. Wyngaard, J.C. Turbulence in the Atmosphere; Cambridge University Press: Cambridge, 2010. [Google Scholar] [CrossRef]
  3. Jensen, N.; Hjort-Hansen, E. Dynamic Excitation of Structures by Wind—Turbulence and Response Measurements at the Sotra Bridge. Report No. STF71 A 1977, 78003. [Google Scholar]
  4. Kristensen, L.; Jensen, N. Lateral coherence in isotropic turbulence and in the natural wind. Boundary-Layer Meteorology 1979, 17, 353–373. [Google Scholar] [CrossRef]
  5. Kaimal, J.C.; Wyngaard, J.; Izumi, Y.; Coté, O. Spectral characteristics of surface-layer turbulence. Quarterly Journal of the Royal Meteorological Society 1972, 98, 563–589. [Google Scholar] [CrossRef]
  6. Jensen, N.O. Simultaneous measurements of turbulence over land and water. Boundary-Layer Meteorology 1978, 15, 95–108. [Google Scholar] [CrossRef]
  7. Mauder, M.; Foken, T.; Aubinet, M.; Ibrom, A. Eddy-Covariance Measurements. In Springer Handbook of Atmospheric Measurements; Foken, Ed.; Springer: Cham, 2021; chapter 55. [Google Scholar]
  8. Mahrt, L. Stably Stratified Atmospheric Boundary Layers. Annual Review of Fluid Mechanics 2014, 46, 23–45. [Google Scholar] [CrossRef]
  9. Kral, S.T.; Reuder, J.; Vihma, T.; Suomi, I.; Haualand, K.F.; Urbancic, G.H.; Greene, B.R.; Steeneveld, G.J.; Lorenz, T.; Maronga, B.; et al. The Innovative Strategies for Observations in the Arctic Atmospheric Boundary Layer Project (ISOBAR): Unique Finescale Observations under Stable and Very Stable Conditions. Bulletin of the American Meteorological Society 2021, 102, E218–E243. [Google Scholar] [CrossRef]
  10. Porté-Agel, F.; Bastankhah, M.; Shamsoddin, S. Wind-Turbine and Wind-Farm Flows: A Review. Boundary-Layer Meteorology 2020, 174. [Google Scholar] [CrossRef] [PubMed]
  11. Veers, P.; Dykes, K.; Lantz, E.; Barth, S.; Bottasso, C.L.; Carlson, O.; Clifton, A.; Green, J.; Green, P.; Holttinen, H.; et al. Grand challenges in the science of wind energy. Science 2019, 366, eaau2027. [Google Scholar] [CrossRef] [PubMed]
  12. Mikkelsen, T.; Sjöholm, M.; Angelou, N.; Mann, J. 3D WindScanner lidar measurements of wind and turbulence around wind turbines, buildings and bridges. In Proceedings of the IOP Conference Series: Materials Science and Engineering. IOP Publishing, 2017, Vol. 276, p.012004. [CrossRef]
  13. Cheynet, E.; Flügge, M.; Reuder, J.; Jakobsen, J.B.; Heggelund, Y.; Svardal, B.; Saavedra Garfias, P.; Obhrai, C.; Daniotti, N.; Berge, J.; et al. The COTUR project: remote sensing of offshore turbulence for wind energy application. Atmospheric Measurement Techniques 2021, 14, 6137–6157. [Google Scholar] [CrossRef]
  14. Ogawa, Y.; Ohara, T. Observation of the turbulent structure in the planetary boundary layer with a kytoon-mounted ultrasonic anemometer system. Boundary-Layer Meteorology 1982, 22, 123–131. [Google Scholar] [CrossRef]
  15. Hobby, M.J. Turbulence Measurements from a Tethered Balloon. PhD thesis, University of Leeds, 2013.
  16. Canut, G.; Couvreux, F.; Lothon, M.; Legain, D.; Piguet, B.; Lampert, A.; Maurel, W.; Moulin, E. Turbulence fluxes and variances measured with a sonic anemometer mounted on a tethered balloon. Atmospheric Measurement Techniques 2016, 9, 4375–4386. [Google Scholar] [CrossRef]
  17. Blanc, T.V.; Plant, W.J.; Keller, W.C. The Naval Research Laboratory’s Air-Sea Interaction Blimp Experiment. Bulletin of the American Meteorological Society 1989, 70, 354–365. [Google Scholar] [CrossRef]
  18. Nambiar, M.K.; Byerlay, R.A.E.; Nazem, A.; Nahian, M.R.; Moradi, M.; Aliabadi, A.A. A Tethered Air Blimp (TAB) for observing the microclimate over a complex terrain. Geoscientific Instrumentation, Methods and Data Systems 2020, 9, 193–211. [Google Scholar] [CrossRef]
  19. Elston, J.; Argrow, B.; Stachura, M.; Weibel, D.; Lawrence, D.; Pope, D. Overview of small fixed-wing unmanned aircraft for meteorological sampling. Journal of Atmospheric and Oceanic Technology 2015, 32, 97–115. [Google Scholar] [CrossRef]
  20. Pinto, J.O.; O’Sullivan, D.; Taylor, S.; Elston, J.; Baker, C.B.; Hotz, D.; Marshall, C.; Jacob, J.; Barfuss, K.; Piguet, B.; et al. The Status and Future of Small Uncrewed Aircraft Systems (UAS) in Operational Meteorology. Bulletin of the American Meteorological Society 2021, 102, E2121–E2136. [Google Scholar] [CrossRef]
  21. de Boer, G.; Argrow, B.; Cassano, J.; Cione, J.; Frew, E.; Lawrence, D.; Wick, G.; Wolff, C. Advancing Unmanned Aerial Capabilities for Atmospheric Research. Bulletin of the American Meteorological Society 2019, 100, ES105–ES108. [Google Scholar] [CrossRef]
  22. Båserud, L.; Reuder, J.; Jonassen, M.O.; Kral, S.T.; Paskyabi, M.B.; Lothon, M. Proof of concept for turbulence measurements with the RPAS SUMO during the BLLAST campaign. Atmospheric Measurement Techniques 2016, 9, 4901–4913. [Google Scholar] [CrossRef]
  23. Mansour, M.; Kocer, G.; Lenherr, C.; Chokani, N.; Abhari, R.S. Seven-Sensor Fast-Response Probe for Full-Scale Wind Turbine Flowfield Measurements. Journal of Engineering for Gas Turbines and Power 2011, 133, 081601. [Google Scholar] [CrossRef]
  24. Calmer, R.; Roberts, G.C.; Preissler, J.; Sanchez, K.J.; Derrien, S.; O. ; apos;Dowd, C. Vertical wind velocity measurements using a five-hole probe with remotely piloted aircraft to study aerosol–cloud interactions. Atmospheric Measurement Techniques 2018, 11, 2583–2599. [Google Scholar] [CrossRef]
  25. Alaoui-Sosse, S.; Durand, P.; Medina, P.; Pastor, P.; Lothon, M.; Cernov, I. OVLI-TA: An Unmanned Aerial System for Measuring Profiles and Turbulence in the Atmospheric Boundary Layer. Sensors 2019, 19, 581. [Google Scholar] [CrossRef]
  26. Witte, B.; Singler, R.; Bailey, S. Development of an Unmanned Aerial Vehicle for the Measurement of Turbulence in the Atmospheric Boundary Layer. Atmosphere 2017, 8, 195. [Google Scholar] [CrossRef]
  27. Wildmann, N.; Hofsäß, M.; Weimer, F.; Joos, A.; Bange, J. MASC – a small Remotely Piloted Aircraft (RPA) for wind energy research. Advances in Science and Research 2014, 11, 55–61. [Google Scholar] [CrossRef]
  28. Wildmann, N.; Ravi, S.; Bange, J. Towards higher accuracy and better frequency response with standard multi-hole probes in turbulence measurement with remotely piloted aircraft (RPA). Atmospheric Measurement Techniques 2014, 7, 1027–1041. [Google Scholar] [CrossRef]
  29. Lenschow, D.H.; Spyers-Duran, P. Bulletin No. 23. Technical Report 23, NCAR Research Aviation Facility, 1989.
  30. Crawford, T.L.; McMillen, R.T.; Dobosy, R.J.; MacPherson, I. Correcting airborne flux measurements for aircraft speed variation. Boundary-Layer Meteorology 1993, 66, 237–245. [Google Scholar] [CrossRef]
  31. Drüe, C.; Heinemann, G. A Review and Practical Guide to In-Flight Calibration for Aircraft Turbulence Sensors. Journal of Atmospheric and Oceanic Technology 2013, 30, 2820–2837. [Google Scholar] [CrossRef]
  32. Vellinga, O.S.; Dobosy, R.J.; Dumas, E.J.; Gioli, B.; Elbers, J.a.; Hutjes, R.W.a. Calibration and quality assurance of flux observations from a small research aircraft. Journal of Atmospheric and Oceanic Technology 2013, 30, 161–181. [Google Scholar] [CrossRef]
  33. Davenport, A.G. The response of slender, line-like structures to a gusty wind. Proceedings of the Institution of Civil Engineers 1962, 23, 389–408. [Google Scholar] [CrossRef]
  34. Shimura, T.; Inoue, M.; Tsujimoto, H.; Sasaki, K.; Iguchi, M. Estimation of Wind Vector Profile Using a Hexarotor Unmanned Aerial Vehicle and Its Application to Meteorological Observation up to 1000 m above Surface. Journal of Atmospheric and Oceanic Technology 2018, 35, 1621–1631. [Google Scholar] [CrossRef]
  35. Hofsäß, M.; Bergmann, D.; Denzel, J.; Cheng, P.W. Flying UltraSonic - A new way to measure the wind. Wind Energy Science Discussions, 2019; 1–21. [Google Scholar] [CrossRef]
  36. Thielicke, W.; Hübert, W.; Müller, U.; Eggert, M.; Wilhelm, P. Towards accurate and practical drone-based wind measurements with an ultrasonic anemometer. Atmospheric Measurement Techniques 2021, 14, 1303–1318. [Google Scholar] [CrossRef]
  37. Bailey, S.C.C.; Sama, M.P.; Canter, C.A.; Pampolini, L.F.; Lippay, Z.S.; Schuyler, T.J.; Hamilton, J.D.; MacPhee, S.B.; Rowe, I.S.; Sanders, C.D.; et al. University of Kentucky measurements of wind, temperature, pressure and humidity in support of LAPSE-RATE using multisite fixed-wing and rotorcraft unmanned aerial systems. Earth System Science Data 2020, 12, 1759–1773. [Google Scholar] [CrossRef]
  38. McConville, A.; Richardson, T. High-altitude vertical wind profile estimation using multirotor vehicles. Frontiers in Robotics and AI 2023, 10, 1–14. [Google Scholar] [CrossRef]
  39. Li, Z.; Pu, O.; Pan, Y.; Huang, B.; Z, Z.; H, W. A Study on Measuring the Wind Field in the Air Using a Multi-rotor UAV Mounted with an Anemometer. Boundary-Layer Meteorology 2023, 187, 1–27. [Google Scholar] [CrossRef]
  40. Prudden, S.; Fisher, A.; Mohamed, A.; Watkins, S. A flying anemometer quadrotor: Part 1. 7th International Micro Air Vehicle Conference and Competition - Past, Present and Future, 2016. [Google Scholar]
  41. Fabrizio Schiano, J.A.M. Towards Estimation and Correction of Wind Effects on a Quadrotor UAV. 2014. [CrossRef]
  42. Paz, C.; Suárez, E.; Gil, C.; Vence, J. Assessment of the methodology for the CFD simulation of the flight of a quadcopter UAV. Journal of Wind Engineering and Industrial Aerodynamics 2021, 218, 104776. [Google Scholar] [CrossRef]
  43. Deters, R.W.; Krishnan, G.K.A.; Selig, M.S. Reynolds number effects on the performance of small-scale propellers. 2014. [CrossRef]
  44. Kutty, H.A.; Rajendran, P. 3D CFD Simulation and Experimental Validation of Small APC Slow Flyer Propeller Blade. Aerospace 2017, 4. [Google Scholar] [CrossRef]
  45. Zheng, Y.; Yang, S.; Liu, X.; Wang, J.; Norton, T.; Chen, J.; Tan, Y. The computational fluid dynamic modeling of downwash flow field for a six-rotor UAV. Frontiers of Agricultural Science and Engineering 2018, 5, 159–167. [Google Scholar] [CrossRef]
  46. Guo, Q.; Zhu, Y.; Tang, Y.; Hou, C.; He, Y.; Zhuang, J.; Zheng, Y.; Luo, S. CFD simulation and experimental verification of the spatial and temporal distributions of the downwash airflow of a quad-rotor agricultural UAV in hover. Computers and Electronics in Agriculture 2020, 172, 105343. [Google Scholar] [CrossRef]
  47. Guillermo, P.P.H.; Daniel, A.M.V.; Eduardo, G.G.E. CFD Analysis of two and four blades for multirotor Unmanned Aerial Vehicle. Institute of Electrical and Electronics Engineers Inc., 2018. [CrossRef]
  48. Lei, Y.; Lin, R. Effect of wind disturbance on the aerodynamic performance of coaxial rotors during hovering. Measurement and Control (United Kingdom) 2019, 52, 665–674. [Google Scholar] [CrossRef]
  49. Ghirardelli, M.; Kral, S.T.; Müller, N.C.; Hann, R.; Cheynet, E.; Reuder, J. Flow Structure around a Multicopter Drone: A Computational Fluid Dynamics Analysis for Sensor Placement Considerations. Drones 2023, 7, 467. [Google Scholar] [CrossRef]
  50. Rankine, W.J.M. On the mechanical principles of the action of propellers. Transactions of the Institution of Naval Architects 1865, 6. [Google Scholar]
  51. Froude, R.E. On the part played in propulsion by differences of fluid pressure. Trans. Inst. Naval Architects 1889, 30, 390. [Google Scholar]
  52. Sayigh, A. Comprehensive Renewable Energy; Vol. 2, Elsevier, 2012.
  53. Jin, L.; Ghirardelli, M.; Mann, J.; Sjöholm, M.; Kral, S.T.; Reuder, J. Rotary-wing drone-induced flow – comparison of simulations with lidar measurements 2023. [CrossRef]
Preprints 95947 g001
Figure 1. Image of the rotary-wing UAV Foxtech D130 X8 and its overall dimensions.
Figure 1. Image of the rotary-wing UAV Foxtech D130 X8 and its overall dimensions.
Preprints 95947 g001
Preprints 95947 g002
Figure 2. Schematics of the mounting rack for the drone with the control station, the power supply unit (PSU), and the unmanned aerial vehicle (UAV).
Figure 2. Schematics of the mounting rack for the drone with the control station, the power supply unit (PSU), and the unmanned aerial vehicle (UAV).
Preprints 95947 g002
Preprints 95947 g003
Figure 3. Mounting rack for the drone assembled to a forklift and the sonic rig used during the experiment (left panel) and view from behind through the 90° tilted drone towards the bay door that was open during the measurements (right panel).
Figure 3. Mounting rack for the drone assembled to a forklift and the sonic rig used during the experiment (left panel) and view from behind through the 90° tilted drone towards the bay door that was open during the measurements (right panel).
Preprints 95947 g003
Preprints 95947 g004
Figure 4. Schematic of the ultrasonic anemometer assembly (left panel) and its realization in action during the measurements (right panel).
Figure 4. Schematic of the ultrasonic anemometer assembly (left panel) and its realization in action during the measurements (right panel).
Preprints 95947 g004
Preprints 95947 g005
Figure 5. Sketch of the measurement setup with perspective A: View across the downwash along the y-axis downstream for the x-z measurement planes; and perspective B: Side-view of the downwash along the x-axis for the y-z measurement planes.
Figure 5. Sketch of the measurement setup with perspective A: View across the downwash along the y-axis downstream for the x-z measurement planes; and perspective B: Side-view of the downwash along the x-axis for the y-z measurement planes.
Preprints 95947 g005
Preprints 95947 g006
Figure 6. Schematics of the tensor structure of the measurement volume for the cross sections in the x-z (left) and y-z planes (right).
Figure 6. Schematics of the tensor structure of the measurement volume for the cross sections in the x-z (left) and y-z planes (right).
Preprints 95947 g006
Preprints 95947 g007
Figure 7. Total PIF across x-z planes measured at different distances Y below the drone at 35% throttle setting. The corresponding distances along the y-axis are given in the lower left corner of each panel, The black circles mark the position and extent of the UAV frame (without propellers).
Figure 7. Total PIF across x-z planes measured at different distances Y below the drone at 35% throttle setting. The corresponding distances along the y-axis are given in the lower left corner of each panel, The black circles mark the position and extent of the UAV frame (without propellers).
Preprints 95947 g007
Preprints 95947 g008
Figure 8. Total PIF across y-z planes measured at different distances X from the vertical plane that cuts through the center of the drone at 35% throttle setting. The corresponding distances along the x-axis are given in the top left corner of each panel
Figure 8. Total PIF across y-z planes measured at different distances X from the vertical plane that cuts through the center of the drone at 35% throttle setting. The corresponding distances along the x-axis are given in the top left corner of each panel
Preprints 95947 g008
Preprints 95947 g010
Figure 10. Total PIF across x-z planes from CFD simulations at different distances Y below the drone.
Figure 10. Total PIF across x-z planes from CFD simulations at different distances Y below the drone.
Preprints 95947 g010
Preprints 95947 g011
Figure 11. Total PIF across y-z planes from CFD simulations at different distances X from the centerline (CL) of the drone.
Figure 11. Total PIF across y-z planes from CFD simulations at different distances X from the centerline (CL) of the drone.
Preprints 95947 g011
Preprints 95947 g012
Figure 12. Visualization of the downwash of a DJI Matrice 300 RTK quadrocopter drone hovering in different heights above a calm lake surface. For levels above 15 D (rotor diameters) the downwash don’t reach the surface (picture to the left). At 7D, a clear single central downwash is visible (picture in the centre) that transforms quickly to four individual downwash zones at around 4 D above the surface (picture to the right). The stills have been extracted from a video gratefully provided by Alizee Lehoux, Uppsala University.
Figure 12. Visualization of the downwash of a DJI Matrice 300 RTK quadrocopter drone hovering in different heights above a calm lake surface. For levels above 15 D (rotor diameters) the downwash don’t reach the surface (picture to the left). At 7D, a clear single central downwash is visible (picture in the centre) that transforms quickly to four individual downwash zones at around 4 D above the surface (picture to the right). The stills have been extracted from a video gratefully provided by Alizee Lehoux, Uppsala University.
Preprints 95947 g012
Table 1. Key technical specifications of the Foxtech D130 X8 UAV.
Table 1. Key technical specifications of the Foxtech D130 X8 UAV.
Dimensions
Width (tip to tip with 28 inch propellers) 1.88 m
Height 0.74 m
Diagonal wheelbase 1.2 m
Weight frame 9 k g
Weight with batteries 15 k g
Frame arm length 0.46 m
Propeller size 28 inch ( 0.71   m )
Propeller pitch 8°
Propulsion System and Autopilot
Speed Controller T-MOTOR Flame 80A ESC
MOTOR T-MOTOR U10II KV100
Propeller Foxtech Supreme 2880 Pro CF
Flight Controller Pixhawk Cube Orange
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated