Articles | Volume 18, issue 13
https://doi.org/10.5194/amt-18-3035-2025
https://doi.org/10.5194/amt-18-3035-2025
Research article
 | 
10 Jul 2025
Research article |  | 10 Jul 2025

A correction algorithm for rotor-induced airflow and flight attitude changes during three-dimensional wind speed measurements made from a rotary unoccupied aerial vehicle

Yanrong Yang, Yuheng Zhang, Tianran Han, Conghui Xie, Yayong Liu, Yufei Huang, Jietao Zhou, Haijiong Sun, Delong Zhao, Kui Zhang, and Shao-Meng Li
Abstract

A hexacopter unoccupied aerial vehicle (UAV) was fitted with a three-dimensional sonic anemometer to measure three-dimensional wind speed. To obtain accurate results for three-dimensional wind speeds, we developed an algorithm to correct biases caused by the rotor-induced airflow disturbance, UAV movement, and attitude changes in the three-dimensional wind measurements. The wind measurement platform was built based on a custom-designed integration kit that couples seamlessly to the UAV, equipped with a payload and the sonic anemometer. Based on an accurate digital model of the integrated UAV–payload–anemometer platform, computational fluid dynamics (CFD) simulations were performed to quantify the wind speed disturbances caused by the rotation of the UAV's rotor on the anemometer during the UAV's steady flight under headwind, tailwind, and crosswind conditions. Through analysis of the simulated data, regression equations were developed to predict the wind speed disturbance, and a correction algorithm for rotor disturbances, motions, and attitude changes was developed. To validate the correction algorithm, we conducted a comparison study in which the integrated UAV flew around a meteorological tower from which three-dimensional wind measurements were made at multiple altitudes. Comparison between the corrected UAV wind data and those from the meteorological tower demonstrated excellent agreement. The corrections result in significant reductions in wind speed bias caused mostly by the rotors, along with notable changes in the dominant wind direction and wind speed in the original data. The algorithm enables reliable and accurate wind speed measurements in the atmospheric boundary layer made from rotorcraft UAVs.

Share
1 Introduction

Wind measurement is crucial in various fields of research and application, including meteorology and environmental sciences. Accurate wind characteristics facilitate modeling of atmospheric transport patterns (Gryning et al., 1987; Stockie, 2011), remote sensing data verification (Drob et al., 2015), model input data assimilation (Gousseau et al., 2011; Vardoulakis et al., 2003), and digital modeling result optimization (Booij et al., 1999; van Hooff and Blocken, 2010). In particular, wind profile measurements near the surface can improve the understanding of atmospheric boundary layer (ABL) dynamics and micrometeorological turbulence at the surface (Seibert et al., 2000), allowing detailed understanding and model descriptions of energy and mass exchanges between air and surfaces and transport processes. The recent development of unoccupied aerial vehicles (UAVs) has provided an opportunity for the measurement of wind fields in three dimensions with high spatial resolutions (McGonigle et al., 2008; Martin et al., 2011; Kim and Kwon, 2019). Their small size, low flight altitude, high mobility, and ability to support sensing devices make UAVs ideal platforms from which to measure wind in the ABL (Thielicke et al., 2021; Shaw et al., 2021; Stewart et al., 2021). Multi-rotor UAVs allow flexible control of flight attitude and stationary hovering, and they can carry varying payloads depending on the number of rotors (Villa et al., 2016; Riddell, 2014; Bonin et al., 2013; Stewart et al., 2021), offering significant advantages in capturing high-resolution wind characteristics in low-altitude conditions (Anderson and Gaston, 2013; McGonigle et al., 2008).

UAVs are often employed to measure wind characteristics both directly and indirectly. Indirect measurement methods involve utilizing pre-installed sensors on the UAV (Elston et al., 2015), in conjunction with specialized flight patterns and wind retrieval algorithms (Bonin et al., 2013; Rautenberg et al., 2018; Gonzalez-Rocha et al., 2019) to achieve wind speed measurement. While these methods offer advantages of operational simplicity and cost-effectiveness, their core principle relies on inversely estimating wind speed through dynamic parameters such as thrust, attitude angles, and flight velocity (Crowe et al., 2020; Donnell et al., 2018; Sikkel et al., 2016; Simma et al., 2020). However, their accuracy is critically dependent on both the measurement precision of the inertial measurement unit (IMU) and the computational reliability of inversion algorithms. Specifically, inherent noise interference in IMU sensors (e.g., a gyroscope's angular rates can be severely affected by external disturbances of up to 0.5° s−1) (Hoang et al., 2021; Neumann and Bartholmai, 2015), combined with uncertainties in parameter configuration within inversion algorithms (the root mean square error (RMSE) of wind speed estimation is 1–1.4) (Bonin et al., 2013), can lead to significant deviations in wind speed estimations. Furthermore, these methods typically assume constant aerodynamic parameters for UAVs, an assumption that often fails to hold in practical complex wind field environments (Bonin et al., 2013).

In contrast, direct measurement methods entail installing additional wind sensors on the UAV to obtain real-time wind information in the field. Porous probes (Soddell et al., 2004; Spiess et al., 2007), pitot tubes (Niedzielski et al., 2017; Langelaan et al., 2011), and anemometers (Rogers and Finn, 2013; Nolan et al., 2018) are commonly used sensors. Sonic anemometers are a more prevalent choice for rotorcraft UAVs, capable of measuring wind speed by detecting changes in the speed of sound travel between different sensors (Thielicke et al., 2021). Recent experiments have demonstrated that under highly turbulent conditions, UAV equipped with properly installed sonic anemometers in wind tunnels can achieve wind speed measurements with RMSEs ranging from 4.3 % to 15.5 % compared to bistatic lidar (Thielicke et al., 2021). Due to the increasing use of rotorcraft UAVs for wind measurements, sonic anemometers are recognized as one of the most promising methods in terms of measurement accuracy and precision.

Sonic anemometers have been mounted onto rotary-wing UAVs for measuring wind speed to varying degrees of success. Typically, an anemometer is mounted at a position along the central axis above the UAV, with data adjusted for the additional wind speed signals induced by UAV motion and attitude changes. Nevertheless, the strong airflow perturbations caused by the rotating propellers can distort real wind flow patterns and significantly affect the accuracy of wind measurements (de Divitiis, 2003). However, these distortions were not considered in the adjustment algorithms. To address this issue, researchers have developed several new correction methods. The first method involves mounting the anemometer along the central axis high above the UAV where the rotor wash effects on the wind speed measurement are believed to be limited (Shimura et al., 2018; Barbieri et al., 2019). Johansen et al. (2015) concluded that, for anemometers at about 40 to 45 mm above the multi-rotor plane of small UAVs, the flow influences from rotors are negligible. However, this method may not be suitable for hexacopters and octocopters due to the high position required, which may raise safety and flight control concerns. The second method involves new corrections based on experiments in an indoor area to measure wind velocity signal bias caused by the rotors during flight and then subtracting the bias (Palomaki et al., 2017). Palomaki et al. (2017) quantified rotor-induced wind speed errors as 0.5 m s−1 compared to tower-mounted anemometers and subtracted these errors from the directly measured wind speed values in subsequent analyses (Palomaki et al., 2017). However, this method is limited by the size of the indoor area, inadequate for full simulations of real UAV rotor speed and attitude changes during flight, and insufficient for the development of a comprehensive correction scheme. Additionally, it does not take into account the detailed coupling of true winds with propeller downwash. The third method is similar to the second except the use of wind tunnels to establish a more accurate relationship between increased airspeed and UAV motion or attitude parameters (Thielicke et al., 2021; Neumann and Bartholmai, 2015). While effective in determining numerical relationships, the method is limited by the high cost of wind tunnel experiments (Dao et al., 2023) and, more importantly, by the additional errors introduced by reflected airflows from the wind tunnel walls and ground (Haleem, 2021; Pettersson and Rizzi, 2008), as well as having the same issues as full simulations of real UAV rotor speed and attitude changes during flight.

The flaws in these correction methods could be addressed using computational fluid dynamics (CFD) simulations to analyze the airflow generated by the UAV's propellers. As far as we know, CFD has been employed to analyze airflow patterns around drones but has not been utilized to correct wind measurements obtained from UAVs (Oktay and Eraslan, 2020; Hedworth et al., 2022). In this paper, we introduce a three-dimensional wind speed correction algorithm for sonic anemometer wind measurements taken from a rotary UAV. This algorithm considers the rotor-induced airflow of the UAV, based on CFD simulations, along with the UAV's motion and attitude changes during flight. The accuracy of the algorithm is confirmed by comparing the corrected wind speeds with those measured from a meteorological tower at multiple altitudes. These results could contribute to ongoing efforts aimed at enhancing the performance and reliability of UAV-based wind speed measurement techniques. Additionally, they pave the way for potential applications, such as quantifying pollutant emissions from industrial complexes (Han et al., 2024).

2 Method

2.1 Digital model representation and simulation tool

2.1.1 Digital model representation

A six-rotor UAV (KWT-X6L-15, ALLTECH, China), equipped with six 32 cm diameter propellers driven by M10 KV100 brushless DC motors, was the platform from which wind was measured. The UAV has a symmetrical motor wheelbase of 1765 mm with an unloaded takeoff weight of 22.5 kg and a maximum flight speed of 18 m s−1. It has a flight endurance >30 min while carrying its maximum payload of 15 kg.

A miniature three-dimensional ultrasonic anemometer (Trisonica-Mini Wind and Weather Sensor, Anemoment, USA) allowed the measurement of wind speed under 15 m s−1 with an accuracy of ±0.1 m s−1 and a resolution of 0.1 m s−1 and wind direction of 0–360° with an accuracy of ±0.1° and a resolution of 0.1°. It was set at 70 cm above the plane of the propellers of the UAV and mounted on a custom-designed carbon fiber tube and frame, which was further mounted onto a rectangular carbon fiber support base attached to the underbelly of the UAV body to minimize the effect of propellers-induced flow on the anemometer measurement. The xtytzt coordinate axes of the anemometer, with its center as the origin, were set to be parallel to the xyz axes of the aircraft body frame. The mounting of the three-dimensional anemometer is shown in Fig. 1a.

A base digital model of the UAV was provided by its manufacturer for the present CFD simulations. The digital model was further augmented with the accurate digital representation of the three-dimensional anemometer and its mounting frame. Furthermore, considering that the UAV wind measurements are usually tied to other air measurement applications, additional payload simultaneously attached to the UAV underbelly was necessary. Such a payload on the UAV needs to also be included in the digital model for the CFD simulation. In the present case, we added the digital model of a 6.37 kg air sampler developed in our group (Yang et al., 2024) to the UAV base digital model (Fig. 1b).

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f01

Figure 1The establishment of the coordinate system and numerical simulation model for the UAV wind measurement platform. (a) The UAV wind speed measurement platform. (b) The digital model of the UAV wind measurement platform in the 3D CFD model simulation domain.

Download

2.1.2 Simulation tool

The CFD simulations were conducted using SolidWorks Flow Simulation 2022, a pressure-based finite-volume solver employing a fully coupled turbulence modeling approach. It employs an adaptive Cartesian mesh approach for three-dimensional solid meshing, with the governing equations being the Navier–Stokes equations for simulating the interaction of fluids, and the turbulence model utilizing the standard kε two-equation model (Jonuskaite, 2017).

The selection of SolidWorks Flow Simulation was driven by its seamless integration with CAD geometries, which eliminated potential errors associated with STL file conversions for our complex multi-rotor UAV design. Additionally, its wall functions for boundary layers effectively resolve gradient variations in boundary layers around rotating blades, reducing trial and error related to near-wall settings. The built-in solver convergence adopts a phased approach to multiple variant scenarios, decreasing the need for re-runs caused by insufficient convergence and thereby conserving computational costs. Its unique turbulence model automatically determines flow regimes (laminar, transitional, and turbulent), ensuring shorter turbulence model setup times while maintaining enhanced model accuracy (Azmi et al., 2017; Ramya et al., 2015).

While ANSYS Fluent offers advanced transient turbulence models (e.g., DES/LES), its computational cost for equivalent spatial resolution was typically higher than that of SolidWorks (Afaq and Ahmad, 2023). Given our need to simulate over 100 operational scenarios, SolidWorks' balance of engineering accuracy and computational tractability was deemed optimal for deriving the empirical correction algorithm.

For CFD simulations, the complete digital model for the UAV and its payloads was set in the xsyszs simulation coordinate system in SolidWorks, on a 1:1 scale (Fig. 1b).

2.2 Simulation scenarios

For the UAV flight simulation, we considered over 100 flight envelope scenarios, including parameters such as UAV flight altitude, wind direction, and wind speed. Since the UAV's predominant flights are within the atmospheric boundary layer, characterized by significant variability in the wind speed and direction, a flight envelope for the UAV in the simulated environments was set up for the complete UAV digital model for flight altitudes of 30 and 1000 m. The lower height (30 m) represents the typical operational altitude for industrial UAV applications within the boundary layer, while the upper height (1000 m) corresponds to the altitude where atmospheric flow transitions to more stable, low-density, freestream conditions. These flight envelopes were designed for the UAV to be subject to headwind, tailwind, and crosswind relative to its flight direction. Under the constraint that the UAV can only operate under true wind speeds ≤18 m s−1 and assuming the applicability of the correction algorithm to most flight scenarios, CFD simulations were conducted for the UAV under these three wind directions. The simulations encompassed the following flight envelopes as listed in Table 1: the UAV flew at ground speeds of 18, 14, 10, and 8 m s−1, respectively, and adapted to wind speeds of 1.5, 3.3, 5.4, 7.9, 10.7, and 14 m s−1. It should be noted that the numerical simulations were conducted by converting wind speed and ground speed into airspeed through vector synthesis.

Table 1The simulation flight envelope scenarios for the UAV-based wind measurement platform.

Download Print Version | Download XLSX

2.3 Flight parameters

The movements of the UAV through air, including takeoff, ascent–descent, attitude changes, turning, and horizontal flights, are driven by the rotary propellers, whose power requirement is closely tied to the weights of the UAV and its payload as well as to the relative motions of the UAV in air. During a normal flight, the UAV adjusts its inclination angle and propeller speeds in order to achieve a set ground speed for flight. By analyzing the gravity G, pull T, and wind resistance D experienced by the UAV under flight conditions (Fig. 2), its inclination angle θ and propeller rotation speed M can be calculated according to Eqs. (1)–(5) (Quan, 2017).

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f02

Figure 2Schematic diagram of forces acting on a UAV.

Download

(1)tanθ×mg=D,(2)p×(sinθ×Sxoy+cosθ×Sxoz)=D,(3)0.5ρ(Vwind+VUAV)2=p,(4)cosθ×mg=T,(5)T=CT×ρ×M602×Dp4,

where θ is the inclination angle of the UAV; m is the combined weight of the UAV and the payloads (i.e., the air sampler and the anemometer along with its installation frame in the present case), calculated to be 28.869 kg; g is the gravitational constant at 9.8 m s−2; D is the wind resistance in newtons; Vwind is the wind speed in m s−1; VUAV is the ground speed of the UAV in m s−1; p is the wind pressure on the UAV in N m−2; Sxoy and Sxoz are the projected surfaces of the UAV in the horizontal direction and vertical directions, determined to be 0.296 and 0.229 m2, respectively; CT is the rotor pull coefficient with an experimentally determined value of 0.048542; Dp is the UAV propeller diameter at 0.8128 m; ρ is the air density in kg m−3; T is each rotor pull in newtons; and M is the rotation speed of the rotors in revolutions per minute (rpm).

The complete flight envelope was defined by combinations of critical parameters, including wind directions, wind speeds, airspeeds, ground speed, inclination, wind resistance, pull, and M. A series of CFD simulations were conducted to systematically evaluate the simulated wind field characteristics for each unique parameter set within this envelope.

2.4 Simulation parameters

The simulation parameters primarily include the computational domain and mesh, and fluid and environmental properties, as well as the rotating region. During the CFD flow simulations of the UAV using SolidWorks, the computational domain dimensions (3.3×3.3×3.3 m3) were determined by prioritizing the analysis of flow field distribution around the anemometer while balancing computational costs. The computational domain was divided into two parts with different spatial resolutions based on the grid sizes, considering the computational time and accuracy required for resolving the details of the digital UAV model. The first part was the global domain with a grid size of 0.23×0.23×0.23 m3, providing a lower spatial resolution. The second part was a nested subdomain within the global domain, specifically defined for the position and dimensions of the anemometer to simulate the measured velocities. The grid size for this nested subdomain was set at 0.0125×0.0125×0.0125 m3, providing a higher spatial resolution. The total number of grids in the computational domain was 1.11×108, and the specific grid configurations are shown in Fig. 3. The wall is set as an adiabatic wall, and its roughness is set to 0. The fluid was modeled as air with characteristics of turbulent and laminar flow. To isolate the rotor-induced flow dynamics from background atmospheric turbulence, a turbulence intensity of 0.1 % and a turbulence length scale of 0.012 m were set. This low turbulence intensity minimizes confounding effects from ambient atmospheric fluctuations, while the length scale corresponds to the anemometer frame width ( 0.01 m) to resolve rotor-generated eddies. These assumptions prioritize the systematic bias correction for rotor-induced airflow. The atmospheric pressure was adjusted to 1.01×105 and 9.00×104 Pa at altitudes of 30 and 1000 m, respectively, and the atmospheric temperature was assumed to be 25 °C at both altitudes. The relative humidity at different altitudes was determined based on the prescribed pressure and temperature corresponding to each altitude. The detailed configurations of these parameters are listed in Table 2.

The UAV's airspeed and aerodynamic angles were configured according to the different flight parameters described in Sect. 2.2 and 2.3. To represent the rotor digitally, six virtual cylinders of the same volume were used to encapsulate the six rotors, with their circumferences matching the rotating trajectory of the propeller tip. These virtual cylinders were treated as the rotational regions in the CFD simulation, with their rotation directions aligned with the actual rotation direction of the UAV's propellers. The rotation direction from rotor nos. 1 to 6 was alternately clockwise and counterclockwise, and the rotation speed for each flight condition was obtained from Eqs. (1)–(5).

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f03

Figure 3Grid configuration of the computational domain.

Download

Table 2Simulation parameter configuration.

Download Print Version | Download XLSX

To ensure relatively accurate simulations, two categories of flow field properties were specified as computational objectives prior to the start of the simulations, and the simulations were terminated upon convergence of the simulation results for all objectives. The first category comprised global-domain computational objectives, including average total pressure (PG), average velocity (VG), average vertical velocity (VGy), and average forward velocity (VGz), where the subscript G denotes the global domain. The second category consisted of subdomain computational objectives, which included the average velocity (Vs) and three-dimensional average speed components Vsx, Vsy, and Vsz at the anemometer position in the simulation coordinate system. It is noteworthy that the aforementioned average values refer to the spatial averages over the global domain or subdomain.

Upon simulation completion, these velocity components (Vsx, Vsy, Vsz) were further converted to velocity components at the anemometer sensor position (ux_sensor, uy_sensor, uz_sensor) according to the coordinate system shown in Fig. 1a and b and Eqs. (6)–(8) below. The converted velocities, ux_sensor, uy_sensor, uz_sensor, were subtracted from the airspeed (denoted as ux_air, uy_air, and uz_air) setting for each CFD simulation to estimate the false wind signals arising from the induced flow by the UAV rotors, expressed with Δux, Δuy, and Δuz, respectively, using Eqs. (9)–(11).

(6)ux_sensor=-Vsz(7)uy_sensor=Vsx(8)uz_sensor=-Vsy(9)Δux=ux_sensor-ux_air(10)Δuy=uy_sensor-uy_air(11)Δuz=uz_sensor-uz_air

In other words, the false wind signals Δux, Δuy, and Δuz are the terms that must be determined and corrected for in the wind measurements from the UAV.

3 Results and discussion

3.1 Example analysis of simulation results

According to Sect. 2.2, this study develops a series of simulation scenarios for the UAV digital model under various combinations of altitude (30 and 1000 m), wind direction (tailwind, headwind, and crosswind), ground speed (8 to 18 m s−1), and wind speed (1.5 to 14 m s−1). To demonstrate the flow field characteristics around the UAV under various scenarios, one UAV hovering scenario and six representative simulation scenarios were specifically selected for analysis as examples.

Figure 4 presents the cross-sectional view of the velocity flow field around the UAV in a hovering state under wind-free conditions. In this scenario, the surrounding velocity field is solely generated by the rotationally induced flow from the UAV's own rotors. The simulated 2–4 m s−1 airflow around the anemometer originates exclusively from rotor rotation, demonstrating that the rotor-induced flow during hovering inherently interferes with wind speed measurements by the anemometer.

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f04

Figure 4The velocity flow field distribution of the UAV's hovering state.

Download

The other six scenarios include UAV flight simulations at altitudes of 30 and 1000 m, with a ground speed of 8 m s−1 and a wind speed of 5.4 m s−1, under tailwind, headwind, and crosswind conditions. Figures 5 to 7 present cross-sectional views of the surrounding flow fields during UAV flight under these conditions. In the figures, color gradients represent the magnitude of the velocity, while arrows indicate both the direction and the magnitude of the velocity. Overall, under varying wind conditions, the direction and speed of airflow around the UAV show significant differences. While the airflow direction around the UAV remains relatively consistent at the both simulation altitudes, the airflow speed at 1000 m is slightly higher than at 30 m, particularly under tailwind conditions. Specifically, based on the ground speed, wind direction, and wind speed settings, the UAV's airspeed relative to the wind is 2.6, 13.4, and 5.4 m s−1 in tailwind, headwind, and crosswind scenarios, respectively.

As shown in Fig. 5a and b, in the tailwind scenario, the maximum downwash velocity occurs directly beneath the UAV rotors, with the airflow directed vertically downward. The next highest velocities are observed around the sides and above the rotors, where the airflow follows an inward and downward trend. The wind speed at the anemometer location is minimally influenced by the UAV rotors, meaning the measured wind speed represents the true airspeed. In the headwind scenario (Fig. 6a and b), the highest airflow velocity is detected near the area directly above the rotors, with the airflow also following an inward and downward pattern. The lowest velocity is found directly below the rotors, where the airflow moves upward and outward. At the anemometer's location, some interference from the UAV rotors is present, so the wind speed at this point is a combination of the true airspeed and the rotor-induced velocity. As exhibited in Fig. 7a and b, in the crosswind scenario with wind blowing from left to right, the airflow around the UAV resembles that in the headwind scenario, but the overall flow field is deflected to the right due to the crosswind, with relatively low airflow velocity. In the scenario with wind blowing from right to left, the flow field shifts to the left.

These simulation results show that the flow field around the UAV varies significantly depending on both the presence/absence of wind and its directional characteristics, and the anemometer experiences different levels of interference accordingly. Thus, accurately quantifying the interference of the UAV rotors on the anemometer is essential. However, in practical application scenarios, it is also necessary to comprehensively consider additional airflow disturbances induced by the UAV's own motion and attitude fluctuations and to develop corresponding dynamic compensation algorithms.

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f05

Figure 5Simulation flow field example results of the UAV wind measurement platform in the tailwind scenario. Panels (a) and (b) represent the longitudinal cross-sections of the simulation flow fields for the UAV at altitudes of 30 and 1000 m, with a ground speed of 8 m s−1 and a wind speed of 5.4 m s−1.

Download

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f06

Figure 6Simulation flow field example results of the UAV wind measurement platform in the headwind scenario. Panels (a) and (b) represent the longitudinal cross-sections of the simulation flow fields for the UAV at altitudes of 30 and 1000 m, with a ground speed of 8 m s−1 and a wind speed of 5.4 m s−1.

Download

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f07

Figure 7Simulation flow field example results of the UAV wind measurement platform in the crosswind scenario. Panels (a) and (b) represent the longitudinal cross-sections of the simulation flow fields for the UAV at altitudes of 30 and 1000 m, with a ground speed of 8 m s−1 and a wind speed of 5.4 m s−1.

Download

3.2 The effect of flight altitude on rotor interference with anemometer measurements

Through simulating the flight of the UAV in all simulation scenarios, the false signals produced by the UAV rotors on the anemometer at different altitudes and wind characteristics were captured. Initially, the influence of flight altitude on the false signals was examined.

The simulated flight data under tailwind and headwind conditions were integrated into a unified dataset since the UAV flight velocity vector is parallel to the tailwind and headwind velocity vectors during normal flight. The simulated false wind signals on the anemometer in the airframe x, y, and z directions, caused by the propeller-induced airflow under tailwind and headwind conditions, were represented by ΔuxT/HW, ΔuyT/HW, and ΔuzT/HW, respectively. For the tailwind and headwind datasets, according to the Wilcoxon nonparametric test for paired samples, as shown in Table 3, the differences in ΔuxT/HW, ΔuyT/HW, and ΔuzT/HW were not significant (p<0.05) at either the 30 m or the 1000 m altitudes. Therefore, in the presence of a tailwind or headwind, the interference from the UAV propeller-induced flow on the anemometer measurement can be considered independent of the flight altitude in this altitude range.

Similarly, the simulated false wind signals for the crosswind conditions on the anemometer in the x, y, and z directions were represented by ΔuxCW, ΔuyCW, and ΔuzCW. The Wilcoxon nonparametric test of paired samples was also applied (shown in Table 1) between the two altitudes. No significant differences were found for ΔuxCW and ΔuzCW between the two altitudes, but there was an obvious discrepancy for ΔuyCW (p=0.00) at the two altitudes. This indicates that, under crosswind conditions, the disturbances of the UAV propeller in the x and z directions of the anemometer are not altitude dependent but that, in the y direction, it is necessary to distinguish the altitude. This behavior can be attributed to differences in the interaction between the y-direction component and rotor rotational momentum caused by variations in atmospheric density at different altitudes under crosswind conditions.

Table 3Wilcoxon nonparametric tests for paired samples of false wind velocity signals between 30 and 1000 m flight altitudes.

Download Print Version | Download XLSX

3.3 Rotor interference with anemometer measurements

This study employs a regression fitting to explore the relationship between the false wind signals generated by the UAV rotors' airflow and the UAV's airspeed. Under tailwind and headwind conditions, the false wind signals (ΔuxT/HW, ΔuyT/HW, and ΔuzT/HW) on the anemometer resulting from the UAV rotor-induced flows at both flight altitudes were aggregated and fitted as dependent variables in a regression using ux_sensor as the independent variable. As shown in Fig. 8a, b, and c, good linear relationships were found between ΔuxT/HW, ΔuyT/HW, and ΔuzT/HW and the simulated velocity components in the x direction (ux_sensor), respectively. The specific relationship is described by Eqs. (12) to (14). Thus, using the UAV velocity components in the x direction, the false wind signals caused by the UAV propellers can be determined and removed from the raw measured wind velocity from the anemometer.

For crosswind conditions, regressions were fitted with false wind signals (ΔuxCW and ΔuzCW) as dependent variables and ux_sensor as the independent variable in the same way (see Fig. 9). A linear relationship was observed between the false wind signals in both x and z directions (ΔuxCWand ΔuzCW) and ux_sensor, with the specific expressions in Eqs. (15) and (16), respectively. As described in Sect. 3.2, ΔuyCW was sensitive to flight altitude under crosswind conditions; hence ΔuyCW at 30 and 1000 m altitude (Δuy(30)CW and Δuy(1000)CW) were regressed against uy_sensor for the two flight altitudes separately. The Δuy(30)CW exhibited a linear relationship with uy_sensor, as shown in Eq. (7). However, the correlation coefficient between Δuy(1000)CW and uy_sensor was found to be lower than 0.5, indicating that Δuy(1000)CW may be considered independent of uy_sensor. Therefore, the average value of Δuy(1000)CW (0.006 m s−1) was regarded as the Δuy(1000)CW at this flight altitude.

Despite the dependence of ΔuyCW on flight altitudes, Δuy(30)CW and Δuy(1000)CW are confined to a similar numeric range. Therefore, they may be roughly considered to represent Δuy for lower altitudes (e.g., 0 to 500 m) and higher altitudes (e.g., 500 to 1000 m), respectively.

Hence, for crosswind situations, the wind velocities in the x, y, and z directions measured by the anemometer are corrected by subtracting ΔuxCW, ΔuzCW, and Δuy(0–500)CW, which are estimated from ux_sensor/uy_sensor, or at relatively high flight altitudes using a constant value of 0.006 m s−1 for Δuy(501–1000)CW.

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f08

Figure 8Regression fit of artificial velocity ΔuxT/HW, ΔuyT/HW, and ΔuzT/HW with ux_sensor for tailwind and headwind flight conditions at two altitudes. In the figure, simulation data are marked with black dots, fitted curves are indicated with black lines, the 95 % confidence bands are identified as green shading, and the 95 % prediction bands are represented with dashed gray lines.

Download

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f09

Figure 9Regression fit of false wind velocity signals ΔuxCW, ΔuzCW, and Δuy(0–500)CW with ux_sensor/uy_sensor for crosswind flight conditions at two altitudes. The symbols in the figure are the same as in Fig. 6.

Download

(12)ΔuxT/HW=0.51+0.061×ux(13)ΔuyT/HW=-0.01+0.70×uy(14)ΔuzT/HW=1.22+0.17×ux(15)ΔuxCW=0.71+0.071×ux(16)ΔuzCW=0.84+0.13×ux(17)Δuy(0–500)CW=-0.0043+0.19×uy(h=0500m)(18)Δuy(501–1000)CW=0.006(h=5011000m)

In Eqs. (17) and (18), the variable h represents the flight altitude of the UAV.

3.4 The overall correction algorithm

3.4.1 Motion and attitude compensation correction of the UAV

In addition to the false wind signals caused by propeller rotations, additional false wind velocity signals from the anemometer can be attributed to UAV movement and attitude (pitch, roll, and yaw) changes during flight, and as such they also need correction. When the UAV moves horizontally and vertically relative to the ground, the velocity vector measured by the anemometer is a vector combination of the true wind velocity and the UAV's ground velocity. Consequently, the ground velocity of the UAV (vx and vz, with vy always being 0 due to no motion in the y direction) contributes false wind velocity components to measurements by the anemometer. Moreover, the UAV's flight attitude undergoes adjustments in the pitch, roll, and yaw Euler angles (θ, φ, and ψ, respectively) in order to compensate for aerodynamic resistance or adapt to flight plans. These adjustments lead to the anemometer measuring additional velocities resulting from the rotational rates of the attitude angles (μ(θ) and μ(φ), with μ(ψ) remaining zero due to the alignment of the rotational axis of ψ with the line connecting the UAV's center of gravity and the anemometer). Furthermore, the effect is further amplified by the distance (r) between the anemometer and the UAV's center of gravity. It is noteworthy that there is currently no reported correction algorithm for the influence of attitude angle variations on anemometer wind velocity measurements from UAVs. To obtain accurate wind information, after eliminating the aforementioned interferences, the wind velocities (ux,uy, and uz) observed by the anemometer in the airframe coordinate (x, y, and z directions) were transformed into the north–east–down (NED) ground coordinate using the direction cosine matrix (DCM) as given in Eq. (19).

(19)uNuEuD=DCM(θφ,ψ)uxuyuz+vx0-vz+-μ(θ)μ(φ)0,(20)DCMθ,φ,ψ=cosψ-sinψ0sinψcosψ0001cosθ0sinθ010-sinθ0cosθ1000cosφ-sinφ0sinφcosφ,

where DCM is defined by Eq. (20); uN, uE, and uD refer to corrected north, east, and down components of wind velocity in the ground coordinate; vx and vz are the motion velocities of the UAV in the x and z directions, respectively, which are directly provided by the GPS receiver output of the UAV or can be directly computed from the UAV longitude/latitude coordinate output; and μ(θ) and μ(φ) represent the product of the pitch rate ω(θ) and roll rate ω(φ), respectively, with the rotation radius r, which is the distance between the anemometer and the center of gravity of the UAV, as defined in Eqs. (21) and (22). Due to the alignment of the anemometer's z axis with that of the UAV, the variation in yaw ψ does not introduce false wind speed into signals from the anemometer in the airframe coordinate, resulting in μ(ψ) being equal to zero.

(21)μ(θ)=ω(θ)×r=d(θ)dt×r,(22)μ(φ)=ω(φ)×r=d(φ)dt×r,

where ω(θ) and ω(φ) are defined as the differentiation of θ and φ, respectively, with respect to time t.

3.4.2 Compensation correction for induced-flow disturbance by UAV rotors

Based on the statistical analyses of the fluid simulation results in Sect. 3.3, the regression relationships between the false wind velocity signals generated by the propeller rotation and the simulated wind components sensed by the anemometer are integrated into the motion and attitude correction algorithm of UAV given in Eq. (19). The updated wind velocity correction algorithm is given as Eq. (23), whose second and third vectors on the right side of Eq. (3) represent the contributions of the propeller-induced wind signals under tailwind/headwind and crosswind conditions, respectively, to ux, uy, and uz, with A and B defined in Eqs. (24) and (25) to quantify their magnitudes. Since the measured wind velocities ux and uy from the anemometer correspond to the simulated ux_sensor and uy_sensor, respectively, the regression relationships are modified by replacing ux and uy with ux_sensor and uy_sensor, respectively. This yields the estimations of the false wind velocity signals, Δux, Δuy, and Δuz, under different wind directions, in relation to ux and uy, as specified by Eqs. (12)–(18). Using Eq. (16), the actual wind velocity components, including north wind (uN), east wind (uE), and vertical wind (uD), are computed after correcting for the effects of the UAV's rotor propeller disturbance, motion, and attitude on the wind signal measurements from the anemometer.

(23)uNuEuD=DCMθ,φ,ψuxuyuz-A×ΔuxT/HWA×ΔuyT/HWA×ΔuzT/HW-B×ΔuxCWB×ΔuyCWB×ΔuzCW+vx0vz+-μ(θ)μ(φ)0(24)A=uxux2+uy2B=uyux2+uy2

3.5 Validation of the correction algorithm

A comparative experiment was designed to verify the effectiveness of the correction algorithm described in Eq. (23). The experiment primarily compares three different kinds of wind data: the first is the three-dimensional wind vector corrected only for UAV motion and attitude compensation (Eq. 19 and denoted as VO); the second includes additional corrections for UAV rotor interference, along with motion and attitude compensation (Eq. 23 and denoted as VR); and the third is the three-dimensional wind directly measured by the meteorological tower (denoted as VT). The comparison experiment was conducted with the UAV flying wind boxes around the 80 m meteorological tower within the Experimental Base of the Beijing Key Laboratory of Cloud, Precipitation and Atmospheric Water Resources. The meteorological tower was equipped with three-dimensional ultrasonic anemometers positioned at heights of 30, 50, and 70 m, with one anemometer in the north and one in the south (see Fig. 10). Experiments were conducted during the daytime on 19 July 2022, with neutral atmospheric stability to minimize thermal boundary layer effects on vertical wind variability.

The UAV flew around the tower in a box flight path at a horizontal distance of about 10 m away from the tower, at all three heights. During these flights, the UAV maintained a commanded horizontal speed of approximately 5 m s−1, a value selected as a compromise between achieving sufficient spatial sampling resolution and maintaining stable flight attitude control. A total of 30 independent wind-box flights were conducted, with each altitude (30, 50, and 70 m) sampled 10 times. Each flight lasted approximately 13 min, generating over 800 valid data points per altitude. Given the potential interference from near-surface vegetation with the 30 m anemometer on the tower, wind velocities acquired by the UAV at 50 and 70 m heights during steady flight intervals were analyzed herein. Using a 3σ threshold of the mean value of the entire dataset to exclude data outliers caused by sudden gusts or UAV maneuvers (such as turning), data during steady UAV flight periods were retained.

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f10

Figure 10Comparative experiment on wind measurements between the UAV and the meteorological tower.

Download

Figure 11 presents the VO, VR, and VT time-series data acquired during the dual-altitude flight tests of the UAV at 70 and 50 m, with the 70 m altitude test data collected prior to 15:05 UTC+8 and the 50 m altitude test data obtained after 15:10 UTC+8. The results in Fig. 11a demonstrate that at elevated wind speeds (>3 m s−1), the wind velocities of VR were substantially lower than those of VO. The root mean square relative errors between VR and VT and between VO and VT are 0.28 and 0.37, respectively, with the former being approximately 24 % smaller than the latter. This indicates that the correction effect of Eq. (23) is especially pronounced in strong wind conditions. In contrast, under gentle wind speeds (≤3 m s−1), VR exhibited greater consistency with VO, but there was still a significant downward revision in the average speed in VR. The average wind speeds of VO, VR, and VT were 2.4, 1.91, and 1.81 m s−1, respectively, with VR exhibiting a 22 % decrease compared to VO. The statistical analysis using the Wilcoxon signed-rank test confirmed a significant difference (p<0.01) in wind speed between VO and VT, whereas no significant differences (p>0.01) were found between VR and VT. This suggests that after compensating for UAV motion, attitude, and rotor interference, the wind speed measured by the UAV anemometer is more closely aligned with that measured directly by the meteorological tower. Moreover, under stronger winds, the wind direction values of VR, VO, and VT were relatively similar, yet at weaker winds, VR showed a small low bias of about 3.3 % (Fig. 11b). The mediocre performance of VR under low wind speeds may originate from the disruption of stable maneuverability in drone rotors caused by low wind speeds, which in turn leads to the failure of the correction algorithm based on CFD steady-state simulations.

Figure 12 presents the wind rose diagrams, offering a detailed overview of the wind speed and direction characteristics for VR, VO, and VT. Compared to the prevailing wind direction frequency (north wind, 39 %) of VT, the dominant wind direction frequency errors for VO and VR are 40 % and 5 %, respectively, demonstrating the superiority of VR in correcting the prevailing wind direction frequency. Meanwhile, deviations in the secondary components introduced by VR (e.g., northwester wind) indicate directions for subsequent model optimization. These analyses indicated that Eq. (23) can effectively correct wind measurement biases induced by UAV disturbances, motion, and attitude changes, particularly at higher wind speeds.

In addition, it should be emphasized that while this study primarily relied on meteorological tower data for algorithm validation, cross-validation through industrial emission scenarios has further confirmed the algorithm's robustness in complex flow fields (Han et al., 2024).

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f11

Figure 11Comparison of wind speed and wind direction time series for VR, VO, and VT. (a) Comparison of wind speed time series for VR, VO, and VT. (b) Comparison of wind direction time series for VR, VO, and VT. (Note that the meteorological tower measurements recorded wind data at 5 s intervals, while the UAV-based measured and corrected wind data were processed with a 10 s sliding average to suppress rotor-induced high-frequency noise, followed by 5 s non-overlapping averaging to align temporally with the tower's 5 s output interval.)

Download

https://amt.copernicus.org/articles/18/3035/2025/amt-18-3035-2025-f12

Figure 12Comparison of wind roses for VO, VR, and VT.

Download

3.6 Discussion on the limitations of the algorithm

The current development of algorithms based on idealized steady-state CFD simulations relies on two key assumptions: low environmental turbulence intensity (0.1 %) and turbulence length scales dominated by anemometer geometric parameters (0.012 m). While this idealized setup effectively isolates rotor-induced flow distortion, its turbulence characteristics fundamentally differ from natural atmospheric conditions. However, it is crucial to emphasize that the algorithm's applicability under turbulent conditions remains valid. This is because rotor-induced wind speed deviations exhibit systemic long-timescale characteristics, whereas atmospheric turbulence primarily affects measurement accuracy through random fluctuations in wind speed and direction with instantaneous nature. This temporal-scale distinction enables our correction algorithm to effectively eliminate systemic biases while minimizing the impact of transient turbulence effects. Nevertheless, it should be noted that under stable atmospheric conditions (low wind speeds) as discussed in Sect. 3.5 or extreme weather scenarios, such airflow environments may disrupt the stable maneuverability of UAV rotors or obscure the systemic drainage effects of rotors, potentially leading to a nonlinear degradation in algorithm accuracy.

In addition, another limitation of our study is the assumption of a smooth surface in CFD simulations, which does not fully capture the impact of surface roughness on wind speed variations near the ground. In reality, surface roughness elements (e.g., vegetation, buildings, or terrain irregularities) alter the wind profile, increasing turbulence and wind shear in the atmospheric surface layer. This effect is particularly relevant to UAV-based wind measurements at low altitudes.

To further enhance the correction algorithm's applicability under diverse environmental conditions, future research will focus on the following aspects: conducting sensitivity studies under different turbulence intensity conditions, implementing supplementary correction modules specifically targeting atmospheric turbulence, and incorporating surface roughness length parameters in future CFD simulations. Although atmospheric turbulence presents significant challenges for UAV-based wind measurements, the correction framework established in this study has demonstrated its effectiveness in improving measurement accuracy across diverse meteorological conditions, thereby laying a critical foundation for developing reliable UAV-based wind measurement systems.

4 Conclusions

The scenarios involving direct measurements of wind fields within the atmospheric boundary layer using multi-rotor UAVs have become progressively commonplace, heightening the significance of accurate wind assessment. However, the rotor propellers during UAV flight introduce additional induced flows at the anemometer location, leading to false wind speed signals. For the present UAV–anemometer–payload configuration, a CFD-based method was used to simulate the process of the UAV wind measurement platform during stable flights under headwind, tailwind, and crosswind conditions. The analyses of induced airflows surrounding the anemometer led to a predictive tool for disturbance airflows. Building upon the UAV motion and attitude correction algorithm, a correction algorithm was proposed for the combined false wind signals from UAV rotor propeller disturbance, motion, and attitude changes during UAV flights. Through comparison of the corrected wind speeds derived from measurements taken from the UAV platform and concurrent three-dimensional wind measurements from a nearby meteorological tower, the validity of the correction algorithm has been demonstrated. Although the algorithm still has certain limitations, it provides a feasible approach for the direct measurement of wind speed from multi-rotor UAV flights.

In conclusion, this study represents a significant advancement in three-dimensional wind speed measurement using UAV platforms, presenting a practical and effective method for direct and accurate wind measurement. This technological breakthrough not only creates a strong foundation for precise wind field measurements with UAVs but also provides potential avenues for the accurate quantification of gaseous pollutant emissions based on UAV. The outcomes of this work carry considerable scientific importance and offer valuable practical applications.

Data availability

Data are available from the corresponding author upon reasonable request.

Author contributions

YY conducted CFD simulations, performed data analysis, and drafted the initial manuscript. YZ contributed to the development of calibration algorithms and provided algorithm programming. YY and TH designed and conducted flight measurement validation experiments under the guidance of SML. DZ provided the site for the flight measurement validation experiments. CX, YL, and YH assisted in refining simulation scenarios and algorithms. JZ, HS, and KZ contributed to the improvement of data analysis methods. SML reviewed and edited the manuscript to ensure the accuracy and completeness of the research.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Financial support

This research has been supported by the National Natural Science Foundation of China Science Fund for Creative Research Groups (grant no. 22221004).

Review statement

This paper was edited by Huilin Chen and reviewed by three anonymous referees.

References

Afaq, M. and Ahmad, R.: Comparison Between ANSYS Fluent and Solidworks Internal Flow Simulation for Analysis of A Fuzzy Logic Controller-Based Heating/Cooling System in A Mobile Robot Design, in: 2023 International Conference on Robotics and Automation in Industry (ICRAI), IEEE, Peshawar, Pakistan, 3–5 March 2023, 1–6, https://doi.org/10.1109/ICRAI57502.2023.10089571, 2023. 

Anderson, K. and Gaston, K. J.: Lightweight unmanned aerial vehicles will revolutionize spatial ecology, Front. Ecol. Environ., 11, 138–146, https://doi.org/10.1890/120150, 2013. 

Azmi, M. F. M., Marzuki, M. A. B., and Bakar, M. A. A.: Vehicle aerodynamics analysis of a multi-purpose vehicle using CFD, ARPN Journal of Engineering and Applied Sciences, 12, 2345-2350, 2017. 

Barbieri, L., Kral, S. T., Bailey, S. C. C., Frazier, A. E., Jacob, J. D., Reuder, J., Brus, D., Chilson, P. B., Crick, C., Detweiler, C., Doddi, A., Elston, J., Foroutan, H., Gonzalez-Rocha, J., Greene, B. R., Guzman, M. I., Houston, A. L., Islam, A., Kemppinen, O., Lawrence, D., Pillar-Little, E. A., Ross, S. D., Sama, M. P., Schmale, D. G., Schuyler, T. J., Shankar, A., Smith, S. W., Waugh, S., Dixon, C., Borenstein, S., and de Boer, G.: Intercomparison of small unmanned aircraft system (sUAS) measurements for atmospheric science during the LAPSE-RATE campaign, Sensors, 19, 2179, https://doi.org/10.3390/s19092179, 2019. 

Bonin, T. A., Chilson, P. B., Zielke, B. S., Klein, P. M., and Leeman, J. R.: Comparison and application of wind retrieval algorithms for small unmanned aerial systems, Geosci. Instrum. Method. Data Syst., 2, 177–187, https://doi.org/10.5194/gi-2-177-2013, 2013. 

Booij, N., Ris, R. C., and Holthuijsen, L. H.: A third-generation wave model for coastal regions: 1. Model description and validation, J. Geophys. Res.-Oceans, 104, 7649–7666, https://doi.org/10.1029/98JC02622, 1999. 

Crowe, D., Pamula, R., Cheung, H. Y., and De Wekker, S. F.: Two supervised machine learning approaches for wind velocity estimation using multi-rotor copter attitude measurements, Sensors, 20, 5638, https://doi.org/10.3390/s20195638, 2020. 

Dao, M. H., Zhang, B., Xing, X., Lou, J., Tan, W. S., Cui, Y., and Khoo, B. C.: Wind tunnel and CFD studies of wind loadings on topsides of offshore structures, Ocean Eng., 285, 115310, https://doi.org/10.1016/j.oceaneng.2023.115310, 2023. 

de Divitiis, N.: Wind estimation on a lightweight vertical-takeoff- and-landing uninhabited vehicle, J. Aircraft, 40, 759–767, https://doi.org/10.2514/2.3155, 2003. 

Donnell, G. W., Feight, J. A., Lannan, N., and Jacob, J. D.: Wind characterization using onboard IMU of sUAS, in: 2018 Atmospheric flight mechanics conference, Atlanta, Georgia, 25–29 June 2018, 2986, https://doi.org/10.2514/6.2018-2986, 2018. 

Drob, D. P., Emmert, J. T., Meriwether, J. W., Makela, J. J., Doornbos, E., Conde, M., Hernandez, G., Noto, J., Zawdie, K. A., McDonald, S. E., Huba, J. D., and Klenzing, J. H.: An update to the Horizontal Wind Model (HWM): The quiet time thermosphere, Earth Space Sci., 2, 301–319, https://doi.org/10.1002/2014EA000089, 2015. 

Elston, J., Argrow, B., Stachura, M., Weibel, D., Lawrence, D., and Pope, D.: Overview of small fixed-wing unmanned aircraft for meteorological sampling, J. Atmos. Ocean. Tech., 32, 97–115, https://doi.org/10.1175/Jtech-D-13-00236.1, 2015. 

Gonzalez-Rocha, J., Woolsey, C. A., Sultan, C., and De Wekker, S. F. J.: Sensing wind from quadrotor motion, J. Guid. Control Dynam., 42, 836–852, https://doi.org/10.2514/1.G003542, 2019. 

Gousseau, P., Blocken, B., Stathopoulos, T., and van Heijst, G. J. F.: CFD simulation of near-field pollutant dispersion on a high-resolution grid: A case study by LES and RANS for a building group in downtown Montreal, Atmos. Environ., 45, 428–438, https://doi.org/10.1016/j.atmosenv.2010.09.065, 2011. 

Gryning, S. E., Holtslag, A. A. M., Irwin, J. S., and Sivertsen, B.: Applied dispersion modelling based on meteorological scaling parameters, Atmos. Environ., 21, 79–89, https://doi.org/10.1016/0004-6981(87)90273-3, 1987. 

Haleem, A.: Vehicular Aerodynamics Wind Tunnel Testing of Unmanned Aerial Multirotor Vehicles and Wall Interference Corrections, PhD Thesis, Toronto Metropolitan University, 2021. 

Han, T., Xie, C., Liu, Y., Yang, Y., Zhang, Y., Huang, Y., Gao, X., Zhang, X., Bao, F., and Li, S.-M.: Development of a continuous UAV-mounted air sampler and application to the quantification of CO2 and CH4 emissions from a major coking plant, Atmos. Meas. Tech., 17, 677–691, https://doi.org/10.5194/amt-17-677-2024, 2024. 

Hedworth, H., Page, J., Sohl, J., and Saad, T.: Investigating Errors Observed during UAV-Based Vertical Measurements Using Computational Fluid Dynamics, Drones, 6, 253, https://doi.org/10.3390/drones6090253, 2022. 

Hoang, M. L., Carratu, M., Paciello, V., and Pietrosanto, A.: Noise Attenuation on IMU Measurement For Drone Balance by Sensor Fusion, in: 2021 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Glasgow, UK, 17–20 May 2021, 1–6, https://doi.org/10.1109/i2mtc50364.2021.9460041, 2021. 

Johansen, T. A., Cristofaro, A., Sorensen, K., Hansen, J. M., and Fossen, T. I.: On estimation of wind velocity, angle-of-attack and sideslip angle of small UAVs using standard sensors, in: 2015 International Conference on Unmanned Aircraft Systems (ICUAS), Denver, CO, USA, 9–12 June 2015, 510–519, https://doi.org/10.1109/icuas.2015.7152330, 2015. 

Jonuskaite, A.: Flow simulation with SolidWorks, Bachelor's Thesis, Plastics Technology, Arcada University of Applied Sciences, Helsinki, Finland, 2017. 

Kim, M.-S. and Kwon, B. H.: Estimation of sensible heat flux and atmospheric boundary layer height using an unmanned aerial vehicle, Atmosphere, 10, 363, https://doi.org/10.3390/atmos10070363, 2019. 

Langelaan, J. W., Alley, N., and Neidhoefer, J.: Wind field estimation for small unmanned aerial vehicles, J. Guid. Control Dynam., 34, 1016–1030, https://doi.org/10.2514/1.52532, 2011. 

Martin, S., Bange, J., and Beyrich, F.: Meteorological profiling of the lower troposphere using the research UAV ”M2AV Carolo”, Atmos. Meas. Tech., 4, 705–716, https://doi.org/10.5194/amt-4-705-2011, 2011. 

McGonigle, A., Aiuppa, A., Giudice, G., Tamburello, G., Hodson, A., and Gurrieri, S.: Unmanned aerial vehicle measurements of volcanic carbon dioxide fluxes, Geophys. Res. Lett., 35, L06303, https://doi.org/10.1029/2007GL032508, 2008. 

Neumann, P. P. and Bartholmai, M.: Real-time wind estimation on a micro unmanned aerial vehicle using its inertial measurement unit, Sensor. Actuat. A-Phys., 235, 300–310, https://doi.org/10.1016/j.sna.2015.09.036, 2015. 

Niedzielski, T., Skjøth, C., Werner, M., Spallek, W., Witek, M., Sawiñski, T., Drzeniecka-Osiadacz, A., Korzystka-Muskała, M., Muskała, P., and Modzel, P.: Are estimates of wind characteristics based on measurements with Pitot tubes and GNSS receivers mounted on consumer-grade unmanned aerial vehicles applicable in meteorological studies?, Environ. Monit. Assess., 189, 1–18, https://doi.org/10.1007/s10661-017-6141-x, 2017. 

Nolan, P. J., Pinto, J., Gonzalez-Rocha, J., Jensen, A., Vezzi, C. N., Bailey, S. C. C., de Boer, G., Diehl, C., Laurence, R., Powers, C. W., Foroutan, H., Ross, S. D., and Schmale, D. G.: Coordinated unmanned aircraft system (UAS) and ground-based weather measurements to predict Lagrangian coherent structures (LCSs), Sensors, 18, 4448, https://doi.org/10.3390/s18124448, 2018. 

Oktay, T. and Eraslan, Y.: Computational fluid dynamics (Cfd) investigation of a quadrotor UAV propeller, International Conference on Energy, Environment and Storage of Energy, Kayseri, Turkey, 5–7 June 2020, 21–25, 2020. 

Palomaki, R. T., Rose, N. T., van den Bossche, M., Sherman, T. J., and De Wekker, S. F.: Wind estimation in the lower atmosphere using multirotor aircraft, J. Atmos. Ocean. Tech., 34, 1183–1191, https://doi.org/10.1175/JTECH-D-16-0177.1, 2017. 

Pettersson, K. and Rizzi, A.: Aerodynamic scaling to free flight conditions: Past and present, Prog. Aerosp. Sci., 44, 295–313, https://doi.org/10.1016/j.paerosci.2008.03.002, 2008. 

Quan, Q.: Introduction to Multicopter Design and Control, Springer, Singapore, 384 pp., https://doi.org/10.1007/978-981-10-3382-7, 2017. 

Ramya, P., Kumar, A. H., Moturi, J., and Ramanaiah, N.: Analysis of flow over passenger cars using computational fluid dynamics, Int. J. Eng. Trends Technol., 29, 170–176, https://doi.org/10.14445/22315381/IJETT-V29P232 2015. 

Rautenberg, A., Graf, M. S., Wildmann, N., Platis, A., and Bange, J.: Reviewing wind measurement approaches for fixed-wing unmanned aircraft, Atmosphere, 9, 422, https://doi.org/10.3390/atmos9110422, 2018. 

Riddell, K. D. A.: Design, testing and demonstration of a small unmanned aircraft system (SUAS) and payload for measuring wind speed and particulate matter in the atmospheric boundary layer, University of Lethbridge, Canada, 2014. 

Rogers, K. and Finn, A.: Three-dimensional UAV-based atmospheric tomography, J. Atmos. Ocean. Tech., 30, 336–344, https://doi.org/10.1175/JTECH-D-12-00036.1, 2013. 

Seibert, P., Beyrich, F., Gryning, S. E., Joffre, S., Rasmussen, A., and Tercier, P.: Review and intercomparison of operational methods for the determination of the mixing height, Atmos. Environ., 34, 1001–1027, https://doi.org/10.1016/S1352-2310(99)00349-0, 2000. 

Shaw, J. T., Shah, A. D., Yong, H., and Allen, G.: Methods for quantifying methane emissions using unmanned aerial vehicles: a review, Philos. T. Roy. Soc. A, 379, https://doi.org/10.1098/rsta.2020.0450, 2021. 

Shimura, T., Inoue, M., Tsujimoto, H., Sasaki, K., and 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, J. Atmos. Ocean. Tech., 35, 1621–1631, https://doi.org/10.1175/JTECH-D-17-0186.1, 2018. 

Sikkel, L., de Croon, G., De Wagter, C., and Chu, Q.: A novel online model-based wind estimation approach for quadrotor micro air vehicles using low cost MEMS IMUs, in: 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), IEEE, Daejeon, Korea (South), 9–14 October 2016, 2141–2146, 10.1109/IROS.2016.7759336, 2016. 

Simma, M., Mjøen, H., and Boström, T.: Measuring wind speed using the internal stabilization system of a quadrotor drone, Drones, 4, 23, https://doi.org/10.3390/drones4020023, 2020. 

Soddell, J. R., McGuffie, K., and Holland, G. J.: Intercomparison of atmospheric soundings from the Aerosonde and radiosonde, J. Appl. Meteorol., 43, 1260–1269, https://doi.org/10.1175/1520-0450(2004)043<1260:IOASFT>2.0.CO;2, 2004. 

Spiess, T., Bange, J., Buschmann, M., and Vorsmann, P.: First application of the meteorological Mini-UAV 'M2AV', Meteorol. Z., 16, 159–170, https://doi.org/10.1127/0941-2948/2007/0195, 2007. 

Stewart, M., Martin, S., and Barrera, N., and Barrera, N. (Eds.): Unmanned aerial vehicles: fundamentals, components, mechanics, and regulations, Unmanned Aerial Vehicles, Nova Science Publishers, Hauppauge, New York, USA, 71 pp., 2021. 

Stockie, J. M.: The mathematics of atmospheric dispersion modeling, SIAM Rev., 53, 349–372, https://doi.org/10.1137/10080991X, 2011. 

Thielicke, W., Hübert, W., Müller, U., Eggert, M., and Wilhelm, P.: Towards accurate and practical drone-based wind measurements with an ultrasonic anemometer, Atmos. Meas. Tech., 14, 1303–1318, https://doi.org/10.5194/amt-14-1303-2021, 2021. 

van Hooff, T. and Blocken, B.: Coupled urban wind flow and indoor natural ventilation modelling on a high-resolution grid: A case study for the Amsterdam ArenA stadium, Environ. Model. Softw., 25, 51–65, https://doi.org/10.1016/j.envsoft.2009.07.008, 2010. 

Vardoulakis, S., Fisher, B. E. A., Pericleous, K., and Gonzalez-Flesca, N.: Modelling air quality in street canyons: a review, Atmos. Environ., 37, 155–182, https://doi.org/10.1016/S1352-2310(02)00857-9, 2003.  

Villa, T. F., Gonzalez, F., Miljievic, B., Ristovski, Z. D., and Morawska, L.: An overview of small unmanned aerial vehicles for air quality measurements: Present applications and future prospectives, Sensors, 16, 1072, https://doi.org/10.3390/s16071072, 2016. 

Yang, Y., Zhou, J., Xie, C., Tian, W., Xue, M., Han, T., Chen, K., Zhang, Y., Liu, Y., and Huang, Y.: A New Methodology for High Spatiotemporal Resolution Measurements of Air Volatile Organic Compounds: From Sampling to Data Deconvolution, Environ. Sci. Technol., 58, 12488–12497, 2024. 

Download
Short summary
A wind speed correction algorithm for multirotor unoccupied aerial vehicles (UAVs) was developed using computational fluid dynamics (CFD). An integrated compensation algorithm was designed to account for the effects of UAV motion, attitude changes, and rotor-induced airflow on wind speed measurements. Comparative experimental results confirmed the effectiveness of the proposed compensation algorithm.
Share