WSPM-System: Providing Real Data of Rotor Speed and Pitch Angle for Numerical Simulation of Downwash Airflow from a Multirotor UAV Sprayer

: The accurate setting of input parameters in the numerical simulation of downwash airflow from a UAV sprayer is important for acceptable simulation results. To provide real data of simulation parameters (rotor speed and pitch angle) for the numerical simulation of downwash airflow, a wireless simulation parameter measurement system (WSPM-System) was designed and tested in this study. The system consists of hardware and software designed based on Arduino and LabVIEW, respectively. Wireless communication was realized by nRF24L01. The lattice Boltzmann method (LBM) was applied for the numerical simulation of downwash airflow. The results showed that the valid communication distance of the WSPM-System was 100 m, with a packet loss rate of less than 1%. While hovering, the rotor speed dropped by about 30% when the load of the UAV sprayer changed from 16 kg to 4 kg, which resulted in the maximum vertical downward velocity ( V VD ) on the horizontal detection surface dropping by about 23%. Under forward flight, the rotor speed in the front ( n 1 , n 6 ) and rear ( n 3 , n 4 ) of the UAV sprayer, respectively, showed a negative linear correlation and positive linear correlation with flight speed (R 2 > 0.95). Meanwhile, the rotor speed in the middle ( n 2 , n 5 ) was consistent with the rotor speed while hovering under the same load; the pitch angle showed a positive linear correlation with flight speed (R 2 > 0.94). A correlation analysis of measured and simulated values of the V VD revealed that the numerical simulation of downwash airflow with the parameters provided by the WSPM-System was reliable (R 2 = 0.91). This study confirmed that the input value of the rotor speed in the fluid software needed to be determined according to the application parameters of the UAV sprayer, thus providing a feasible method and system for obtaining real simulation parameters. of the UAV sprayer’s downwash airflow and pesticide applications. This study has contributed to improving the reliability of this system’s numerical simulation results.


Introduction
The chemical control of pests and diseases is an integral activity in agricultural production systems, and it plays an important role in increasing the quality and yield of crops [1,2]. In recent years, a new spray application technology, based on an unmanned aerial vehicle (UAV) sprayer, has been rapidly developed, particularly in East Asian countries (such as China, South Korea, and Japan) [3], and it achieved good application results in field crops such as rice paddies, corn, wheat, and cotton [4][5][6][7]. In addition, the spraying of orchards with multirotor UAV sprayers has also been reported [8][9][10][11]. Additionally, the data (coverage, residue, and penetration) collected from one of the most difficult pest control systems (almond orchard) in US agriculture support reported successes in China [12]. Compared with a single-rotor helicopter (fuel-powered or electric-powered), electric multirotor UAV sprayers with a smaller payload (10-16 kg) are highly efficient, highly flexible, and affordable, and they are rapidly being commercialized for crop protection in small, fragmented, and complex terrains, as well as in irregular farms, thus increasing their worldwide interest [12,13]. Additionally, the reduced exposure of workers to pesticides makes the UAV sprayer more attractive [10]. At present, more than 90% of agricultural UAV sprayers in China used for crop protection are electric multirotor UAV sprayers [14]. In addition, the use of UAV sprayers has been extended to fertilizer applications [15], the release of predatory mites and other natural enemies [16], and for sowing rice seeds [17]. Moreover, during the COVID-19 pandemic, the UAV sprayers have been used for spraying disinfectants [18].
Although UAV sprayers have been adopted in commercial spray scenarios and have emerged as one of the main pesticide application approaches in China [19], there is still a lack of environmental risk assessment studies regarding off-target losses. The off-target losses of pesticides, including spray drift, result in heavy contamination of the ecosystem, including water resources and soil, as well as human health threats [20][21][22][23]. Based on European Regulation 1107/2009 [24], the objective of protecting human and animal health and the environment should take priority over the objective of improving crop production. European Directive 2009/128/CE [25] pointed out that the impact of pesticide uses on human health and the environment should be assessed by conducting an environmental risk assessment. However, in addition to field tests [14,26], a more efficient way to achieve this goal is to simulate the UAV sprayer's downwash airflow and pesticide applications by using more efficient and low-cost numerical simulation technologies. This means higher requirements for the accuracy of the numerical simulation of downwash airflow using fluid software, such as Fluent, XFlow, etc. Hence, the input values of simulation parameters (rotor speed and pitch angle) in fluid software are very important.
In recent years, the downwash airflow from UAV sprayers has become a topic of interest. The use of numerical simulation technology to study downwash airflow is not only important for understanding spray drift but can also provide a rich database for the environmental risk assessment of drifted pesticides. Different sources of airflow can lead to different spatiotemporal distribution characteristics for airflow; as such, insights into the spatiotemporal distribution of airflow are an important prerequisite in the analysis of droplet deposition dynamics. Regarding traditional air-assisted sprayers, Qi et al. [27] found that the airflow from an air-assisted sprayer could effectively improve the distribution of droplets in tree canopies. Gu et al. [28] studied the unimpeded air velocity profiles of an air-assisted five-port sprayer. Their results showed that the airflow distribution pattern was different with various application parameters. Hong et al. [29] developed an integrated CFD model to predict air velocity distributions inside and around tree canopies blown by an air-assisted pesticide sprayer. Similar to the horizontal and/or vertical airflow from air-assisted sprayers, the downwash airflow of the UAV sprayer plays an important role in the transportation of droplets from their source (nozzles) to the target [30]. Strong airflow can increase on-target deposition by reducing droplet flight time and reducing the effect of meteorological conditions on the spray [29]. Inevitably, downwash airflow can increase the risk of droplets' off-target losses [26,30], and changes the shape of plants [31][32][33].
Studies of downwash airflow from multirotor UAV sprayers are essential for increasing our understanding of the performance characteristics of UAV sprayers. Researchers used an anemometer to measure downwash airflow in a real application of a UAV sprayer [34][35][36] and carried out a large number of field spray tests [4][5][6][7][8][9][10][11][12]37]. These data not only reflect the final effect of downwash airflow on droplet deposition and drift more realistically, but also provide a reference for the reasonable selection of UAV sprayer types and their application parameters. However, due to uncontrollable environmental conditions (ambient wind speed magnitude and direction, etc.) and limited test costs, test equipment, and test methods, it is not yet possible to obtain the complete spatiotemporal distribution of downwash airflow and local vortex patterns of UAV sprayers through field tests [19]. Numerical simulation technology can effectively solve the above problems and complement field measurement. To this effect, this approach has become a mainstream research tool for the downwash airflow from UAV sprayers.
At present, researchers have simulated the downwash airflow from single-rotor [38][39][40], quadrotor [41][42][43], six-rotor [30,[44][45][46][47], and eight-rotor [48] UAV sprayers. These research contents were mainly focused on the spatiotemporal distribution of the downwash airflow and its influence on droplet movement under various application parameters (flight speed, flight altitude, load, etc.) of specific UAV sprayers. However, more importantly, these published studies seldom pay attention to the accuracy of the input values of the simulation parameters in fluid software. On the one hand, it is important to acquire acceptable simulation results; on the other hand, multirotor UAV sprayers are equipped with fixed pitch propellers that are only operated by varying rotor speed. For a specific multirotor UAV sprayer, the rotor speed and pitch angle are different under various application parameters, resulting in different downwash airflows as determined by flight principles [49]. Generally, a stronger downwash airflow can enhance spray penetration in canopies and reduce spray drift from natural wind [26]. However, current studies on the numerical simulation of downwash airflow have not paid enough attention to the above problem. This is not conducive to obtaining reliable simulation results.
To date, simulation parameters (rotor speed and pitch angle) have been obtained through the following paths: (1) calculating them according to the dynamic model; (2) Output data from flight records; (3) direct measurement with the measuring system. Path 1 is mainly used in the theoretical design stage of multirotor UAVs, which may not reflect the data of UAV sprayers in real applications; some parameters involved in the dynamic model are still determined through complex experiments [50,51]. Path 2 remains the most viable path for obtaining simulation parameters, but it is limited by the fact that the flight record of UAV sprayers in the market does not contain data on rotor speed and pitch angle, and the flight controller of commercial UAV sprayers is not open source. Therefore, it is quite difficult to output the pitch angle value from the attitude sensor provided by the UAV sprayer. More importantly, various sensors equipped with UAV sprayers cannot output rotor speed values. However, rotor speed is an important input parameter in the numerical simulation of downwash airflow using fluid software, such as Fluent, XFlow, etc. Considered comprehensively, Path 3 not only obtains the actual data of rotor speed and pitch angle in UAV sprayer applications, but it also does not require complicated flight control and attitude calculation expertise; thus, it is the most feasible choice at present. However, most of the existing measurement equipment comprises platforms in the laboratory to measure the dynamic parameters of one rotor or a mini-aircraft [52][53][54][55]; this is not suitable for monitoring the simulation parameters of the multirotor UAV sprayer in real-time.
Currently, only a few published studies on the numerical simulation of downwash airflow have concerned themselves with changes in rotor speed and pitch angle of multirotor UAV sprayers under varied application parameters. This work aims to monitor the real simulation parameters of multirotor UAV sprayers in an application through a specially designed wireless simulation parameter measurement system (WSPM-System). The system consists of hardware and software designed based on Arduino and LabVIEW, respectively. Wireless communication was realized by nRF24L01. Meanwhile, to evaluate system performance and collect simulation parameter data, a system test experiment and a simulation parameters measurement experiment were carried out. Eventually, the data obtained from WSPM-System were successfully applied in the numerical simulation of downwash airflow according to the lattice Boltzmann method (LBM).

Multirotor UAV Sprayer
Experimental tests were performed using an electric six-rotor UAV sprayer (DF-16L, Henan Difengde Aviation Technology Co., Ltd., Zhoukou, China), which is commonly used in China (Figure 1). A tank with a 16 L capacity was equipped on the UAV sprayer with four replaceable nozzles. When spraying different crops, the application parameters were different according to the manufacturer's recommendation. In addition, the system was fitted with a GPS-guided autonomous flight device to remotely control the UAV sprayer in each spray application. The specifications and parameters of the UAV sprayer are listed in Table 1.

Wireless Measurement System
The hardware and software of the WSPM-System were designed with a keen consideration of real application scenarios. The hardware was for data collection and transmission, whereas the software was for data reception, display, and storage. In this work, the hardware design was according to the Arduino (www.arduino.cc, accessed on 21 December 2020) open-source platform with low cost and a short development cycle, and used sensors and the wireless module compatible with the open-source philosophy. The Lab-VIEW (www.ni.com, accessed on 01 September 2021) platform was adapted to develop the software. The system framework is as shown in Figure 2.

Hardware Design
Considering the real-time monitoring of rotor speed and pitch angle data of the UAV sprayer in real applications, the hardware needed to be installed on the UAV sprayer. On this note, it was vital for the size and weight of the hardware to be as small as possible in addition to having an anti-interference capability in the dynamic environment. The design had six identical Hall sensors that were for monitoring the rotor speed of the six-rotor UAV sprayer. Each Hall sensor individually monitored the change in one rotor speed. Compared to a photosensitive sensor, the Hall sensor can effectively avoid the impact of dust and light and is more suitable for the application environment of a UAV sprayer. An attitude sensor (JY61, WitMotion Shenzhen Co., Ltd., Shenzhen, China) was selected to monitor the change in pitch angle, which was equipped with the MPU6050. The MPU6050 devices combine a 3-axis gyroscope (±250/500/1000/2000°/s) and a 3-axis accelerometer (±2/4/8/16 g) on the same silicon die. The attitude sensor JY61 integrates an attitude solver, the Kalman filter algorithm, to acquire accuracies of 0.05° and 0.1° under static and dynamic conditions, respectively (www.wit-motion.com, accessed on 19 March 2021).
In addition to sensors for monitoring data, the hardware of the WSPM-System also included a wireless module, an Arduino board, and a lithium battery. The nRF24L01 of a low-power wireless transceiver chip was selected as the wireless module, which worked in the 2.4 GHz universal ISM band. This module contains the transmitter used by the MCU (TTL-nRF24L01) and the receiver used by the computer (USB-nRF24L01). An Arduino Mega 2560 (Rev3) was used as a datalogger due to its memory capacity and the number of IO ports. It is a microcontroller board based on the ATmega2560 and has 54 digital input/output pins, 16 analog inputs, and contains everything needed to support the microcontroller (www.arduino.cc, accessed on 21 December 2020). More importantly, it has six digital pins for interrupt, which was required for simultaneous monitoring of the six rotors in real-time. To connect multiple components on an Arduino Mega 2560 board, an expansion board was selected to expand the number of Vin and GND pins. In addition, a 9 V/11,800 mAh lithium battery was used to directly power the system through the direct current (DC) power jack. The basic parameters and appearance of the hardware can be seen in Table 2.
The code was developed in the open-source Integrated Development Environment (IDE) of Arduino. The function format of external interrupt was 'attachInterrupt (interrupt, function, mode)'. As shown in Equation (1), it evaluates the rotor speed by calculating the time required to complete one lap. The RISING was selected to trigger when the pin went from low to high. Additionally, it contained the Wire library required for the I 2 C communication protocol used by the MPU6050; this open-source library can be found in the GitHub repository (www.github.com, accessed on 15 March 2021). Although the roll angle, yaw angle, and pitch angle data can be obtained from the JY61, only the pitch angle data were needed as output data in this study. All the data of rotor speed (1−6) and pitch angle were sent to the serial port where the nRF24L01 was connected according to the preset sampling interval, using the Arduino serial commands where n is the rotor speed, rpm; microst is the trigger moment of the t-th, microsecond; microst+1 is the trigger moment of the (t + 1)-th, microsecond. Component 3144 MPU6050 nRF24L01+ ATmega2560 Data detected a n θ Input (Output) voltage (V) 7-12 (5 or 3.3) (5, 9, or 12) 11,800 mAh Appearance a n (rotor speed, rpm); θ (pitch angle, °). The range of pitch angle is (−90°, +90°), and the value is positive in the forward flight of the UAV sprayer.
The hardware integration and installation effects are shown in Figure 3. The Arduino board was packaged with a 3D printed shell, and the interface of six Hall sensors and one JY61 were reserved. Then, it was mounted on the landing gear of the UAV sprayer together with the lithium battery. According to the rotor numbers in Figure 3f, the Hall sensors were installed in sequence. Since the rotor was fixed on the brushless motor to rotate synchronously, the magnet was fixed on the brushless motor. To ensure measurement accuracy, the horizontal mounting distance of the Hall element to the magnet fixed on the brushless motor was <5.0 mm [30]. In addition, since only the pitch angle data were collected from the attitude sensor, JY61, the installation of the JY61 did not require its coordinate system to coincide with the fuselage coordinate system of the UAV sprayer. However, the XOY plane of the JY61 coordinate system should be parallel to the X'O'Y' plane of the fuselage coordinate system as much as possible. At the same time, the directions of X and X', Y and Y' axes should be consistent, respectively. The total weight of the hardware was about 0.95 kg, accounting for 2.6% and 2.3% of the standard take-off weight and maximum take-off weight, respectively; this did not affect the flight performance of the UAV sprayer. In addition, this total weight was counted in the load of the simulation parameters measurement experiment.

Software Design
The graphical user interface (GUI) of the software was designed based on LabVIEW, and the VISA tool was used to realize serial communication to receive data. VISA is a protocol to provide communication between LabVIEW and any other device. The GUI and program block diagram of the WSPM-System can be seen in Figure 4. The data from the hardware were inputted through the serial wireless module at a 115,200 baud rate. The seven parameters (rotor speed 1-6 and pitch angle) were displayed on the front panel of LabVIEW in the form of a linear graph and table. Meanwhile, the data were written into a file (Excel.CSV) in real-time. The detailed application is divided into three main stages: hardware preparation, connection to GUI, and data collection ( Figure 5). In the first stage, the code runs just after powering up the Arduino; after that, the system initialization is completed through establishing a communication between the sensors and the devices. Finally, a time-based sample of the sensors is performed, and then data transmission is completed according to the preset sampling interval in a loop mode. In the second stage, the VISA tool completes initialization after starting the GUI of the software. Then, we selected the correct COM device and baud rate, and then clicked 'OPEN' to complete the connection and started receiving data continuously. In the third stage, the data are written into the linear graph, table, and Excel.csv, respectively. We then clicked 'CLOSE' to stop the data writing, and then clicked 'OPEN' to restart the data writing and generate a new Excel.csv or exit the system. The data recorded in Excel.csv are then ready for data analysis.

Wireless Simulation Parameter Measurement System Test Experiment
Firstly, the communication performance of the system was tested. Considering the actual application demand of the system, wireless communication was used to transmit data; the variables tested mainly included sampling interval and communication distance. The evaluation index of the test was packet loss rate. Richer monitoring data can be obtained with a higher sampling frequency, and the application of the system within the valid communication distance can ensure that the data are collected in real-time and completely by the software. To this effect, the tests were carried out under communication distances of 80, 100, and 120 m. There were no obvious obstacles at the tests site. Meanwhile, another variable that needed to be controlled was the sampling interval. The study considered intervals of 0.1, 0.5, and 1.0 s.
Secondly, the stability of the data collected by the system was tested under good communication performance. Since six identical Hall sensors were used, only the rotor speed data monitored by one of the Hall sensors were tested. The rotor speed data collected from the WSPM-System were compared with the reference data collected from the RC41 tachometer (±0.5‰ to ±1.5‰ rpm, Jintan Runchen Electric Appliance Factory, Changzhou City, Jiangsu Province, China). Considering that the accuracy of the system may be affected by rotor speed size, the tests were carried out under five gradients. In addition, to verify whether the attitude sensor was accurately positioned, the pitch angle data monitored by the attitude sensor were tested under three gradients. The UAV sprayer was placed on a flat ground, inclined 0°, 5°, and 10°, respectively. A laser rangefinder (±0.3°, MILESEEY ® , Shenzhen, China) was used to calibrate these angles.

Simulation Parameters Measurement Experiment
To collect simulation parameters of the UAV sprayer under various application parameters, the simulation parameters measurement experiment was carried out based on the specially designed WSPM During the experiments, meteorological data were monitored in real-time using a Kestrel ® 4500 Weather & Environmental Meter (Nielsen-Kellerman Co., Boothwyn, PA, USA). Weather parameters including temperature, relative humidity, ambient wind speed magnitude, and direction were recorded every two seconds. To avoid the impact of downwash airflow from the UAV sprayer, the weather meter was about 40 m from the nearest flight route. The maximum ambient wind speed throughout the day occurred in the afternoon (after 12:00) with a southerly wind (160°-220°), less than 3 m s −1 , and there was almost no ambient wind in the morning (7:00-11:00). To minimize the impact of ambient wind on the flight performance of the UAV sprayer, the measurement window with almost no ambient wind was selected for each treatment according to the wind speed displayed by the meteorological meter in real-time. Once the ambient wind speed exceeded 1 m s −1 , the experiment was stopped immediately, and the next measurement window waited. As shown in Figure 6b, a laptop with the software of the WSPM-System was placed on the western side of the experimental area 10 m away from the UAV sprayer flight route (from A to B) to collect data. The forward flight distance was maintained at about 150 m. This allowed the system to not only collect rich data but also to collect data that were transmitted stably within the valid communication distance.
Based on the flight principle of multirotor UAV [49], different loads need different rotor speeds to provide lift, which leads to different downwash airflow. In theory, different flight speeds correspond to different pitch angles, and the pitch angle is 0° in the hovering state. Meanwhile, a change in flight altitude can change the spatiotemporal distribution of downwash airflow, and the lift of the rotor may change due to the ground effect [56], which leads to change in rotor speed. Therefore, the simulation parameters measurement experiment was conducted in three steps. Step 1-The simulation parameters of the UAV sprayer while hovering were measured. Since the downwash airflow had a significant backward tilt when the UAV sprayer was flying forwards [30], the hovering state was more susceptible to ground effects under the same load and flight altitude. Within the tank capacity of the UAV sprayer, four load gradients were selected (LG = 4, 8, 12, and 16 kg). Considering the normal operation height recommended by the manufacturer, three flight altitude gradients were selected (FA = 2, 3, and 4 m) for each load gradient, which resulted in a total of twelve treatments. Considering the proportion of the hardware weight of the WSPM-System, the evaluated weights of water added to the tank were 3.05, 7.05, 11.05, and 15.05 kg. The measurement position was at point A (Figure 6b), and, for each treatment, after the UAV sprayer was hovering steadily, the data were collected continuously, and the time was about one minute.
Step 2-The simulation parameters of the UAV sprayer in forward flight were measured. Based on the tests and analysis obtained in the first step, the flight altitude was fixed at 3 m in this step. Similar to the treatments in the first step, four load gradients of 4, 8, 12, and 16 kg were selected. According to the normal flight speeds recommended by the manufacturer, seven flight speed gradients (FS = 1, 2, 3, 4, 5, 6, and 7 m s −1 ) were selected for each load gradient, which resulted in a total of twenty-eight treatments. As shown in Figure 6b, for each treatment, after the hovering of the UAV sprayer stabilized at point A, the software of the WSPM-System was opened to collect data and forward flight after about ten seconds. The mean value of the pitch angle collected within the ten seconds was used to correct the pitch angle in the forward flight. After about 150 m forward flight, the software of the WSPM-System was closed at point B.
Step 3-The vertical downward velocity (VVD) of downwash airflow from the UAV sprayer in hovering was measured. Previous studies [30,44] indicate that the VVD is the main body of the airflow velocity in the downwash airflow, which plays a key role in droplets' delivery to the target during the spraying and can improve the deposition of droplets in the crop canopy. As shown in Figure 6c, airflow velocity was measured using six hot wire anemometers (Testo 405i, Titisee-Neustadt, Germany). The six hot wire probes (#1, #2, #3, #4, #5, and #6) were evenly distributed on circles of a radius of 0.82 m and located directly below rotors 1, 2, 3, 4, 5, and 6. To avoid the effect of the ground spreading airflow, the hot wire probe was 1 m away above the ground [30]. The anemometer recorded data every two seconds and transmitted them via Bluetooth to an application on a smartphone. For each case, after the UAV sprayer was hovering steadily at point A, the Smart Probe APP was opened to collect data and then closed after about one minute. Subsequently, the calculated average velocity was recorded as VVD of the marked point [30]. Finally, the data of VVD were used in a numerical approach for validation.

Numerical Simulation of Downwash Airflow
The simulation parameter data from the WSPM-System were used in the numerical simulation of downwash airflow from the UAV sprayer. A total of six cases considering application parameters (flight altitude and load) in hovering were selected to evaluate the distribution characteristics of downwash airflow. Initial conditions are listed in Table 3. Table 3. Application parameters of the UAV sprayer in hovering and size of the computational domain.

Case Application Parameters Computational Domain (m) FA (m)
LG (kg) Length Width Height  1  2  16  12  12  4  2  3  16  12  12  5  3  4  16  12  12  6  4  3  4  12  12  5  5  3  8  12  12  5  6  3  12  12  12  5 A full-scale three-dimensional (3D) model of the UAV sprayer was created for numerical simulation. As shown in Figure 7, the fuselage was modeled in Creo Parametric 2.0 (PTC Inc.) software; a 3D optical scanner (Zhanhao Industrial Design Co., Ltd., Shenzhen, China) was used to scan the rotor to generate point cloud data, and then the rotor's 3D model was reconstructed by Geomagic Studio (Geomagic Inc., Research Triangle Park, North Carolina, US) software. The computational domain was set up as a virtual wind tunnel forming a cuboid with length × width × height (12 m × 12 m × (FA + 2) m). The lattice Boltzmann method (LBM) was used to simulate the downwash airflow from a mesoscopic point of view. This approach not only avoids the numerical dissipation problem that may be caused by manually dividing the grid but also has a strong dynamic adaptive refinement ability for the prediction of turbulent wakes. More detailed information about the turbulence model, boundary conditions, and spatial discretization used in this study is shown in a previous study [30]. At t = 1.3 s, the airflow in the computational domain fully dispersed in all cases studied, and this study adopted this value as the simulation time.

Optimized Extraction of Simulation Parameter Data
The simulation parameter data from the WSPM-System needed to be optimized and extracted to reduce accidental errors and ensure reliability. For each treatment, seven parameters with ni (i = 1, 2, 3, 4, 5, 6) and θ were included in the original data set. At first, the data of each parameter was arranged in ascending order; then 10% of the data in both the front and rear sides was deleted. The remaining 80% of the data were used to calculate the mean value of each parameter. The calculating equations are as follows: where ni is the mean value of i-th rotor speed (rpm); θ is the mean value of pitch angle (°);  (2) and (3), the acquired values were subsequently used for the data analysis and numerical simulation.

Analysis Method of Collected Data Stability
The coefficient of variation (CV) was used as an indicator for the stability of the collected data in this paper. A lower CV value indicates better stability of collected data. For comparing the variability among several data groups without the same unit of measurement or the mean value, the CV was applied more meaningfully and widely than the standard deviation [57]. This is because the CV is the relative variation and a dimensionless measure. The calculating formula is as follows: where CV is the coefficient of variation (%); M is the value of ni or θ (rpm or °); SD is the standard deviation value of rotor speed or pitch angle (rpm or °).

Regression Analysis
Unary linear regression analysis for each treatment in the simulation parameters measurement experiment was performed to produce the regression equations with FS as an independent variable and ni, and θ as dependent variables, respectively. For the simple linear regression model, it is assumed that the error has constant variance. However, when fitting the experimental data, the instrument error may need to be considered (which reflects the accuracy and precision of a measuring instrument). Therefore, the assumption of constant variance in the error is violated. Thus, it is necessary to assume εm to be normally distributed with nonconstant variance, and the error acts as variance, which can be used as a weight in fitting (www.originlab.com/doc/Origin-Help, accessed on 15 April 2021). The equations are as follows: where β1 is the slope and β0 is the intercept; ε is a random error term; N is the total number of data points; xm and ym are the independent variable and dependent variable of m-th data point, respectively; wm is the weighting factor of m-th data point; wm The regression analysis was performed using the OriginPro 2018C (OriginLab Corporation, Northampton, MA, USA) software package.

Numerical Approach to Validation
The lattice Boltzmann method (LBM) has been successfully used in previous studies [30,[38][39][40][41]. Nevertheless, taking into account the difference between the UAV sprayer type and the measurement scenario, Step 3 in the simulation parameters measurement experiment was carried out to verify the reliability of the LBM. Initially, for each case, six sensors were established in the computational domain to monitor and output the simulated values of VVD over time. The position of the sensor corresponded to the hot wire probe in Step 3. After the downwash airflow stabilized (from 1.0 to 1.3 s), the simulated values were extracted and averaged, which were the final simulated values of the marked point. Then, a comparison between the measured and simulated values of the mean VVD was carried out through correlation analysis.

Performance Test of the WSPM-System
As shown in Table 4, data loss was not found within a communication distance of 80 m. When the communication distance was 100 m, the packet loss rate was 1%; when the communication distance was 120 m, the packet loss rate increased, but not above 4%. The test results further showed that the sampling interval had no great effect on data transmission. In addition, the communication performance tests were also carried out at 140 m, but the results are excluded in Table 4. This was because the communication distance may have reached the limit of the nRF24L01, and there was almost no data transmission. Based on the above results and analysis, although this system can be used for data collection at 120 m, the valid communication distance adopted in this study was 100 m to ensure a good transmission effect in a dynamic environment. To obtain the richest data within a limited valid communication distance, a sampling interval of 0.1 s was considered in the simulation parameters measurement experiment. The results of the stability of the collected data are shown in Table 5. Comparing the results of rotor speed measured under five gradients, there was no significant difference in the mean value (M), the standard deviation (SD), and the coefficient of variation (CV) between WSPM-System and RC41. This indicated that the data collection of the system designed in this study was reliable. Although within the test range the SD value increased to a certain extent with the measured rotor speed, the increase was not more than ±23 rpm. The minimum and maximum values of CV from the WSPM-System and RC41 appear in Text.1 (0.46% and 0.33%, respectively) and Text.2 (1.42% and 1.05%, respectively). Comparing the results of pitch angle measured under three gradients (0°, 5°, and 10°), the M values calculated from the data collected from WSPM-System were 0.36°, 4.97°, and 10.14°, respectively. All SD values were 0.01°; however, they decreased successively to 1.82%, 0.15%, and 0.10%, respectively. Considering that the angle measurement accuracy of the laser range finder used was not high, this implied that the JY61 installation position was relatively correct. Therefore, in the simulation parameters measurement experiment, the pitch angle was also corrected. CV 0.46% 0.33% * 1.42% 1.05% * 0.87% 0.80% * 0.58% 0.51% * 0.51% 0.58% * 1.82% 0.15% 0.10% Note: M (Mean value); SD (Standard deviation); CV (Coefficient of variation); '*' means the reference data from RC41.

Simulation Parameters in Hovering
For Step 1, as can be seen in Figure 8, the six rotor speeds were basically same for each treatment, with SD values not more than ±27 rpm. Meanwhile, the six rotor speeds had the same distribution pattern under the three flight altitude gradients with the same load. This was translated to mean that the flight altitude had no effect on the rotor speed within the measurement range, i.e., no ground effect was found. Note that, due to limited experimental scenarios, these results apply to similar crop and ground conditions (the average height of wheat about 7 cm) as well as this specific type of UAV sprayer. However, it was observed that the rotor speed changed with the load. When the loads were 4, 8, 12, and 16 kg, the rotor speeds were maintained about at 2097, 2404, 2706, and 3001 rpm, respectively. As the load changed from 16 kg to 4 kg, the rotor speed dropped about 30%. This observation was in agreement with the analysis results based on its flight principle. Therefore, these results suggested that the input value of the rotor speed in the fluid software needed to be determined according to the application parameters of the UAV sprayer. Figure 9 shows the measured pitch angle values in the hovering UAV sprayer as a box chart. It can be seen that the total mean value of the pitch angle was −1.41°. This means that even if the UAV sprayer were hovering, there may be differences between the measured pitch angle and the theoretical value (0°). Therefore, to obtain a more accurate pitch angle value of the UAV sprayer at various flight speeds, not only does the JY61 attitude sensor need to be installed in the correct position, but the measurement of the pitch angle in hovering is also necessary for its correction in forward flight.

Simulation Parameters in Forward Flight
For Step 2, the mean values of measured rotor speed in forward flight of the UAV sprayer are presented in Figure 10. All the SD values were not more than ±30 rpm. As can be seen in Figure 10, under four load gradients, six rotor speeds (n1, n6; n3, n4; n2, n5) showed a similar trend with an increase in FS. Specifically, the rotor speeds in the front (n1 and n6) decreased with an increase in FS and the rotor speeds in the rear (n3 and n4) increased with an increase in FS. It was noted that the rotor speeds in the middle (n2 and n5) did not change with an increase in FS. Combining these values with the regression equations analysis, it can be found that R 2 values more than 0.95 (maximum, 0.998; minimum, 0.955) indicate a negative linear correlation between the rotor speeds in the front (n1, n6) and FS; there was a positive linear correlation between the rotor speeds in the rear (n3, n4) and FS. On the contrary, R 2 values less than 0.12 (maximum, 0.118; minimum, 0.003) indicated that there was no correlation between the rotor speeds in the middle (n2, n5) and FS. More importantly, n1, n6, and n3, n4 decreased or increased at almost the same rate. Based on these data, as the FS increased from 1 m s −1 to 7 m s −1 , the rotor speeds in the front and rear changed about 10% (n1 and n6 decreased, n3 and n4 increased).  Figure 11 shows the measured pitch angle values in the forward flight of the UAV sprayer, which had been corrected using the collected data in hovering. It can be seen that, for all load gradients, the pitch angle increased with an increase in FS. Combining the regression equations analysis, it can be found that R 2 values greater than 0.94 (maximum, 0.945; minimum, 0.992) indicate a positive linear correlation between the pitch angle and FS. On the whole, the pitch angle increased from 0.79° to 5.80° when the FS was increased from 1 m s −1 to 7 m s −1 ; these values were calculated based on the mean values of the four load treatments.  Figure 12 shows the vorticity distribution diagram of downwash airflow while hovering from the front view. The axial view and the top view are given in Figure A1. The results of Step 1 revealed that the change in flight altitude did not affect rotor speed. However, the impact of flight altitude on the downwash airflow distribution can be observed in Figure 12a-c. Specifically, as the flight altitude increased, the vorticity under the fuselage decreased, which implied that the vortex motion of the downwash airflow with a lower flight altitude was stronger at the same load. In addition, the range of ground spreading airflow also decreased to a certain extent with increased flight altitude ( Figure  A1a-c). More importantly, the change in flight altitude led directly to the change in downwash airflow velocity. Figure 13 presents the distribution of VVD on the horizontal detection surface; the distance between the detection surface and the ground is 1 m. The negative value of VVD indicates that the velocity direction is vertical upwards. As shown in Figure 13a-c, as the flight altitude increased, the distribution of larger VVD was transferred from the surrounding area to the central area. For FA = 2 m, the direction of VVD in the central area was vertical upwards with a maximum value of 7.9 m s −1 (Figure 13a). This phenomenon did not appear in other cases.

Simulation Results of Downwash Airflow
Additionally, the results of Step 1 revealed that the load affected the rotor speed. From Figure 12d, e, f, and b, it was observed that, as the load increased, the vorticity under the fuselage increased, meaning that the vortex motion of the downwash airflow with a larger load was stronger at the same flight altitude. In addition, the range of ground spreading airflow also increased to a certain extent with the increase in load ( Figure A1d, e, f, and b). For different loads, the distribution of VVD on the horizontal detection surface is shown in Figure 13d, e, f, and b, respectively. It was observed that VVD increased with the load. It is worth noting that, when the loads were 4, 8, 12, and 16 kg, the maximum VVD values were 13.8, 16.2, 18.3, and 17.9 m s −1 , respectively. The maximum value of VVD on the horizontal detection surface dropped about 23% when the load was changed from 16 kg to 4 kg. For LG = 16 kg, although the maximum value of VVD decreased compared with the other three load gradients, the VVD values greater than 15 m s −1 had a larger distribution range. Moreover, as shown in Figure 14, an R 2 value of 0.91 indicated a good correlation between the measured and simulated values of the VVD.

Discussion
The use of UAV sprayers for spraying pesticides has rapidly developed in China and has been widely used in commercial applications in recent years. According to preliminary statistics from the National Agro-Tech Extension and Service Center, as of mid-November 2020, the number of UAV sprayers in use for agriculture was about 80,000, and the operating area was nearly 53 million hectares [58]. Given its advantages, the environmental risks caused by off-target losses of pesticides due to factors such as downwash airflow should be concerns for farmers, researchers, manufacturers, regulators, and decision makers. However, very few published studies on the numerical simulation of downwash airflow consider the value of rotor speed and pitch angle as inputs in fluid software based on the application parameters of UAV sprayers. The WPSM system designed in this study can contribute to providing accurate input parameters for the future numerical simulation of downwash airflow from UAV sprayers, and it is expected to contribute to indepth understandings of downwash airflow characteristics in real spray applications.
In experiments, not only was the performance of the WSPM-System tested, but so was the simulation parameter data of the multirotor UAV sprayer when hovering and in forward flight. According to the result obtained from the communication performance test, the valid communication distance of the WSPM-System was 100 m with a packet loss rate of less than 1%, which meets the basic requirements of the simulation parameters measurement experiment set out in this paper. In addition, the test results of the stability of the collected data indicated that data collected using the WSPM-System designed in this study was reliable. For rotor speed, there was no significant difference in the M, SD, and CV values between the WSPM-System and RC41 (Table 5). For each treatment, the pitch angle in hovering was used for pitch angle correction in forward flight. The objective of this step was to ensure that the pitch angle data of the UAV sprayer could be accurately monitored by JY61.
When the UAV sprayer was hovering, the impact of flight altitude on rotor speed was not found (Figure 8). In other words, the ground effect did not affect the rotor lift. On the contrary, load had a significant effect on rotor speed, which dropped about 30% when the load was changed from 16 kg to 4 kg. Since the three flight altitude gradients selected in this study were the conventional application altitudes of the UAV sprayer, the impact of flight altitude on the next measurements of stimulation parameters with a similar scenario (field crops) can be ignored. However, it is worth noting that the application technique of UAV sprayers in field crops has become increasingly mature, more and more manufacturers and service organizations have paid attention to their application in economic crops of higher additional value, especially orchards [26]. On the one hand, there is a huge difference in the crop canopy between field crops and orchards; on the other hand, tree canopies have an important effect on the distribution characteristics of downwash airflow [46,47]. Therefore, when the UAV sprayer is used for spraying in orchards, it is still unclear whether the flight altitude can be ignored in the measurement of simulation parameters.
When the UAV sprayer was in forward flight, an impact of flight speed on the rotor speed and pitch angle was found (Figures 9 and 10). More importantly, the rotor speed in a symmetrical position showed the same law (n1and n6; n3 and n4; n2 and n5). This means that, in more subsequent simulation parameters measurement experiments, for quadrotor, six-rotor, and eight-rotor UAV sprayers with a symmetrical layout, only half of the rotors need to be monitored. Although the applicability of the WSPM-System to other types of multi-rotor UAV sprayers has been expanded, it is necessary to note that all rotors need to be monitored when an experiment is carried out in a dynamic environment. Under the impact of a crosswind or a complicated flight attitude, the rotor speed may change according to different laws. In addition, the pitch angle showed a positive linear correlation with FS (R 2 > 0.94), which increased from 0.79° to 5.80° when the FS increased from 1 m s −1 to 7 m s −1 . The flight speed of 1-7 m s −1 corresponded to the pitch angle of 1-7°, which was used by Wen et al. [41] in their simulation of a quadrotor UAV sprayer for pesticide applications. This difference can be attributed to the different types of UAV sprayer.
The data of the rotor speed obtained from the WSPM-System were successfully applied in the numerical simulation of downwash airflow. The numerical simulation results revealed the impact of flight altitude and load on downwash airflow velocity. The correlation analysis of measured and simulated values of the VVD revealed that the numerical simulation of downwash airflow was reliable (R 2 = 0.91). The distribution of larger VVD transferred from the surrounding area to the central area as the flight altitude increased. Only FA = 2 m, the vertical upward of VVD, appeared in the central area, and the maximum value was 7.9 m s −1 . This observation can be explained by the evolution of the downwash airflow based on a previous study [30]. The evolution of the downwash airflow from the top to the ground produced a spiral wake vortex, a turbulent wake, a spreading airflow, and a vortex ring. Induced by the wingtip vortex, a spiral wake vortex appeared below the fuselage. As the airflow developed downward, the six spiral wake vortices approached each other so that serious interference occurred. This destroyed the spiral wake vortex structure, causing it to lose its coherence and diffuse into a turbulent wake. When FA = 2 m, the horizontal detection surface was located in the area where the six spiral wake vortices approached each other and caused interference, while the horizontal detection surface in other cases was located in the turbulent wake area. This also revealed that the spiral wake vortex was only distributed in a fixed area below the rotor, and was not affected by flight altitude.
The rotor speed dropped by about 30% as the load changed from 16 to 4 kg, which resulted in the maximum value of VVD on the horizontal detection surface dropping by about 23%. This means that, in a single spray application process of the UAV sprayer, the downwash airflow velocity gradually decreased as the liquid in the tank was released continuously; however, it was not clear whether the droplet deposition could be changed, making this phenomenon worthy of further exploration. For LG = 16 kg, the maximum value of VVD was smaller than for LG = 12 kg, which may be caused by the turbulence characteristics of downwash airflow. Moreover, the range of ground spreading airflow decreased or increased to a certain extent with the increase in flight altitude or load, respectively. This was caused by the difference in the lateral velocity of the ground spreading airflow. Simultaneously, this difference was also the reason for the different canopyairflow vortex structures caused by the downwash airflow in the rice canopy [31].
For future developments, we identified the limitations of the WSPM-System that should be improved or supplemented. For instance, although the valid communication distance of 100 m was suitable for the simulation parameters experiment in this paper, the demand for a longer valid communication distance can be realized by replacing the wireless communication accessories. Secondly, the WSPM-System's lithium battery (0.53 kg) accounts for 55.79% of the total hardware weight (0.95 kg). To expand the applicability of the WSPM-System, further work can be performed to highlight ways to reduce the weight of the lithium battery as much as possible under the premise of satisfying system performance through power consumption tests. In addition, we identified works that were worthy of being carried out based on the WSPM-System. The real data from WSPM-System will be used in the numerical simulation of downwash airflow from unmanned aerial vehicle (UAV) sprayers in an orchard to further explore the interaction effects between the downwash airflow and tree canopies.

Conclusions
In this work, to provide real data of simulation parameters (rotor speed and pitch angle) for the numerical simulation of downwash airflow, a WSPM-System was designed; the performance test results show that the system met the study requirements. The data of simulation parameters obtained from the WSPM-System revealed the effect of the UAV sprayer' application parameters on rotor speed and pitch angle. On the one hand, the rotor speed dropped by about 30% when the load was changed from 16 kg to 4 kg while hovering. On the other hand, when the flight speed was increased from 1 m to 7 m s −1 , the rotor speeds in the front and rear changed by about 10% (n1 and n6 decreased, n3 and n4 increased), and the pitch angle increased from 0.79° to 5.80°. Eventually, the data of the rotor speed obtained from the WSPM-System were successfully applied in the numerical simulation of downwash airflow. Additionally, the validated numerical simulation results show that, as the load changed from 16 kg to 4 kg, the maximum value of VVD on the horizontal detection surface dropped by about 23% due to the rotor speed dropping. This result means that the input value of the rotor speed in the fluid software needs to be determined according to the application parameters of the UAV sprayer. Based on the results in this study, the WSPM-System can be used to provide real data of simulation parameters for numerical simulations of the UAV sprayer's downwash airflow and pesticide applications. This study has contributed to improving the reliability of this system's numerical simulation results.