A method for predicting the deposition and drift of aerial spraying droplets of a strip-shaped composite crown layer

Abstract

The application discloses a strip-shaped composite canopy aerial application droplet deposition drift prediction method, relates to the technical field of plant protection spraying simulation, and obtains crop canopy and unmanned aerial vehicle parameters of a field experiment as first simulation parameter information; a rotor slip grid model, a dynamic parameter porous medium canopy model and a canopy droplet probability capture model are respectively constructed by using the first simulation parameter information; a three-dimensional model of a calculation domain is constructed in combination with the rotor slip grid model and the dynamic parameter porous medium canopy model, and grid division is carried out; second simulation parameter information of strip-shaped composite canopy crop simulation spraying is set; and the second simulation parameter information, the three-dimensional model of the calculation domain after grid division and the canopy droplet probability capture model are combined to carry out simulation statistics and analyze the downwash flow field distribution and the droplet deposition law under a strip-shaped composite irregular and heterogeneous canopy. The influence of the downwash air flow and the special canopy structure on the droplet migration and deposition is analyzed, so as to provide a reference for the strategy optimization of strip-shaped composite canopy unmanned aerial vehicle application.

Description

Technical Field

[0001] The present application relates to the technical field of plant protection spray simulation, and in particular to a method for predicting the deposition and drift of droplets during aerial spraying of strip-shaped composite canopies. Background Art

[0002] The soybean-corn strip intercropping system effectively increases soybean yields while ensuring stable corn production, making it a crucial measure for ensuring national grain and oil security. Large-scale implementation of this intercropping system urgently requires support from precision pesticide application strategies. Traditional field droplet deposition experiments typically use indirect measurement methods such as water-sensitive paper and polyester film sheets to measure droplet deposition distribution within the crop canopy. However, field experiments often face numerous limitations, including high costs and complex sources of error.

[0003] In recent years, CFD (Computational Fluid Dynamics) numerical simulation has become increasingly widely used in agricultural aerial spray research. This method uses Lagrangian particle tracking to simulate droplet migration trajectories. Chinese patent application number CN202211506013.7 discloses a CFD-based method for simulating inter-row wind-delivered spraying. This method uses a porous media model to represent the crop canopy structure and successfully simulates the deposition of droplets on the porous media surface under the influence of airflow. However, this method fails to reflect the dynamic changes in canopy resistance properties under the influence of auxiliary airflow. Furthermore, this method only deposits droplets on the porous media canopy surface, failing to simulate the deposition process within the canopy.

[0004] At present, there is an urgent need for a numerical simulation method of UAV droplet drift and deposition applicable to the irregular and heterogeneous strip-shaped composite canopy of soybean and corn. Summary of the Invention

[0005] In order to solve the above technical problems, this application proposes the following technical solutions:

[0006] In a first aspect, the present application provides a method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy, comprising:

[0007] Obtain crop canopy and drone parameters from the field experiment as the first simulation parameter information;

[0008] Using the first simulation parameter information, a rotor sliding grid model, a dynamic parameter porous media canopy model, and a canopy droplet probability capture model are constructed respectively;

[0009] Combining the rotor sliding grid model and the dynamic parameter porous media canopy model to construct a three-dimensional model of the computational domain and perform grid division;

[0010] Setting the second simulation parameter information of the strip composite canopy crop simulation spray;

[0011] Combining the parameter information, the three-dimensional model of the computational domain after grid division and the canopy droplet probability capture model, simulation statistics are performed to analyze the downwash flow field distribution and droplet deposition law under the strip-shaped composite irregular and heterogeneous canopy.

[0012] In a possible implementation, obtaining crop canopy and drone parameters from a field experiment as first simulation parameter information includes:

[0013] Through the aerial spraying experiment, the wind field data under the rotor under standard load was obtained, as well as the wind speed data of the airflow into and out of the two crop canopies at different flight altitude and flight speed combinations;

[0014] The three-dimensional size of the equivalent porous medium canopy model of the crop is determined by field measurements of the crop strip composite planting pattern, crop plant height, and plant leaf area parameters, and the canopy leaf area density is calculated:

[0015]

[0016] in: LAD is the canopy leaf area density, S leaf represents the total leaf area of ​​the plant, S Indicates the area occupied by the plant, h Indicates plant height;

[0017] Determine the spray atomization parameters based on the actual drone nozzle model and operating flow rate.

[0018] In a possible implementation, the use of the first simulation parameter information to respectively construct a rotor sliding grid model, a dynamic parameter porous media canopy model, and a canopy droplet probability capture model includes:

[0019] By combining the inverse modeling method and comparing the simulated flow field data at different speeds with the measured data, the standard rotor speed is determined and the rotor sliding grid model is constructed;

[0020] Determine the basic equation of the airflow velocity at the canopy outlet and determine the undetermined parameter function relationship of the basic equation of the airflow velocity at the canopy outlet, and obtain the dynamic parameter porous medium canopy model in combination with the porous medium model preset in the simulation software;

[0021] The droplet deposition behavior in the canopy space is divided into two processes: droplet blade contact and adhesion. The deposition probability formula of the droplet canopy in each grid unit is obtained, and a canopy droplet probability capture model coupled with dynamic parameters is constructed.

[0022] In one possible implementation, the inverse modeling method is used to determine the rotor standard speed by comparing the simulated flow field data at different speeds with the measured data, and to construct the rotor sliding grid model, including:

[0023] Use an optical scanner to collect rotor point cloud data, and import the obtained point cloud data into reverse modeling software to reconstruct the rotor three-dimensional model;

[0024] The rotor three-dimensional model is used to carry out a rotor downwash flow field simulation experiment through sliding grid technology. The rotor speed is determined by comparing the measured flow field data under the same operating conditions with the simulated flow field data under different speeds, and the rotor sliding grid model is obtained.

[0025] In one possible implementation, determining the basic equation for the airflow velocity at the canopy outlet and determining the undetermined parameter function relationship of the basic equation for the airflow velocity at the canopy outlet, and obtaining the dynamic parameter porous medium canopy model in combination with a porous medium model preset in the simulation software, includes:

[0026] Taking the canopy inlet airflow velocity, leaf area density, inertial drag coefficient and porosity as input features, the calculation equation for the canopy outlet airflow velocity is derived:

[0027]

[0028] in: V is the predicted value of the airflow velocity at the canopy outlet; is the airflow velocity at the canopy entrance; is the canopy inertial drag coefficient; is the canopy porosity; k is the number of iterations;

[0029] The residuals of wind speed data of two crop canopies were collected through computational simulation and field experiments. The functional relationship between the inertial drag coefficient and porosity change was constructed using the least squares method based on the trust region, and the undetermined parameters in the canopy outlet airflow velocity equation were identified. The objective function of the least squares method is as follows:

[0030]

[0031] The iterative calculation of the trust region algorithm is shown as follows:

[0032]

[0033]

[0034] in: is the parameter vector to be identified, is the optimal identification parameter vector obtained using the least squares method; is the simulated measurement value of the airflow velocity at the canopy outlet; is the parameter vector at the kth iteration when the least squares method is iteratively calculated; The equation for the airflow velocity at the canopy outlet is Jacobian matrix; m is the number of experimental data groups; is the mth error function in the parameter vector The gradient vector at ; and Respectively The first and second parameters in the vector;

[0035] The functional relationship between the above-mentioned inertial resistance coefficient and the canopy porosity change is loaded into the porous medium model preset in the simulation software to obtain the dynamic parameter porous medium canopy model.

[0036] In one possible implementation, the droplet deposition behavior in the canopy interior is divided into two processes: droplet contact and adhesion. The deposition probability formula of the droplet canopy in each grid unit is obtained, and a dynamic parameter-coupled canopy droplet probability capture model is constructed, including:

[0037] According to the actual working conditions in the field, the fixed parameters of blade unit diameter, leaf area density, droplet density and air dynamic viscosity are determined;

[0038] Obtain the droplet size, velocity inertia, and airflow velocity at the center of each grid cell, calculate the contact probability and adhesion probability between each grid cell in the canopy and the droplet in real time, and finally calculate the capture probability:

[0039]

[0040]

[0041]

[0042]

[0043] in: is the capture probability; is the contact probability between the droplet and the leaf, is the droplet velocity; C 1 is the model empirical coefficient; is the time variation; is the probability of droplets adhering to leaves; is the droplet density; d is the droplet diameter; u is the air velocity; is the dynamic viscosity of air; D leafIndicates the width of vegetation leaves;

[0044] A canopy droplet deposition model with differential distribution of internal spatial deposition probability was constructed based on the deposition probability calculation formula.

[0045] In a possible implementation, the constructing of a three-dimensional computational domain model by combining the rotor sliding grid model and the dynamic parameter porous media canopy model, and performing grid division, includes:

[0046] Ignoring the fuselage of the rotor sliding mesh model, a multi-rotor equivalent representation of the plant protection UAV is adopted;

[0047] Based on crop phenotypic parameters, canopy wind resistance properties and droplet capture probability, a zoned differentiated porous medium equivalent method was used to characterize the crop canopy.

[0048] A multi-level grid division strategy is adopted to appropriately encrypt the grids of the rotor, sprinkler, crops and their surrounding areas.

[0049] In a possible implementation, the second simulation parameter information for setting the strip-shaped composite canopy crop simulation spray includes:

[0050] The DPM under the Lagrangian framework is used to simulate the droplet motion trajectory and the spray atomization parameters are set according to the actual working conditions;

[0051] The boundary pressure around the crops and near the drone is set to allow the droplets to escape freely, and the ground boundary is set as a no-slip wall that can capture the droplets.

[0052] In one possible implementation, combining the second simulation parameter information, the gridded three-dimensional computational domain model, and the canopy droplet probability capture model to perform simulation statistics and analyze the downwash flow field distribution and droplet deposition pattern under the zonal composite irregular and heterogeneous canopy, includes:

[0053] Importing the second simulation parameter information, the three-dimensional model of the computational domain after grid division, and the canopy droplet probability capture model into the post-processing software CFD-post and introducing the droplet deposition amount and deposition density variables;

[0054] The CFD-post simulation is used to statistically analyze the distribution characteristics of the irregular heterogeneous canopy downwash flow field and the droplet drift and deposition characteristics.

[0055] One possible implementation method also includes: through simulation experiments and analysis under a variety of typical operating conditions, the downwash flow field distribution and droplet deposition patterns under the crop strip-shaped composite irregular and heterogeneous canopy are revealed, and the drone flight control strategy is optimized.

[0056] In this application example, the dynamic canopy drag coefficient and porosity changes under the interaction of downwash airflow and the crop canopy are simulated, and the droplet and airflow kinetic energy parameters of each grid cell are obtained to simulate the droplet deposition process within the canopy space with dynamic parameter coupling. The influence of downwash airflow and special canopy structure on droplet migration and deposition is analyzed, revealing the spatial distribution patterns of downwash flow and droplet deposition, providing a reference for optimizing drone spraying strategies in zonal composite canopies. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] Figure 1 A schematic flow chart of a method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy provided in an embodiment of the present application;

[0058] Figure 2 A schematic diagram of the process of constructing a rotor sliding grid model provided in an embodiment of the present application;

[0059] Figure 3 A schematic diagram of constructing a 3D rotor model provided in an embodiment of the present application;

[0060] Figure 4 A schematic diagram of the three-dimensional model of the computational domain and its grid division provided in an embodiment of the present application;

[0061] Figure 5 The downwash airflow distribution diagram provided in the embodiment of the present application;

[0062] Figure 6 Schematic diagram of airflow distribution at different levels of soybean and corn canopies provided in the embodiments of this application;

[0063] Figure 7 A side view of the droplet deposition distribution provided in an embodiment of the present application;

[0064] Figure 8 Schematic diagram of droplet deposition distribution at different levels of soybean and corn canopy provided in the examples of this application. DETAILED DESCRIPTION

[0065] The present solution will be described below with reference to the accompanying drawings and specific implementation methods.

[0066] See also Figure 1 The present embodiment provides a method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy, including:

[0067] S101: Acquire crop canopy and drone parameters of a field experiment as first simulation parameter information.

[0068] In this example, corn and soybean crops were used at maturity. The experimental drone was a DJ-T30 six-rotor drone equipped with 16 fan-shaped hydraulic atomizing nozzles (model SX11001), with a maximum spray flow rate of 7.2 L / min.

[0069] In order to obtain the first simulation parameter information, an aerial defense experiment was conducted. A hot wire anemometer (405i, TESTO) was used to collect wind field data below the rotor under a standard load, as well as wind speed data of airflow into and out of the two crop canopies at different flight altitude and speed combinations.

[0070] Through actual field measurements, it was found that the planting patterns of soybean and corn crops were 30 cm × 10 cm and 40 cm × 10 cm (row spacing × plant spacing), and the plant heights were 1.0 m and 2.9 m, respectively. The leaf area density of soybean and corn canopies was calculated by the following formula: 3.80 m -1 and 6.03 m -1 .

[0071]

[0072] in, LAD is the canopy leaf area density; S leaf represents the total leaf area of ​​the plant; S Indicates the area occupied by the plant; h Indicates plant height;

[0073] Based on the actual operating conditions of the drone and the nozzle model, the nozzle flow rate was determined to be 0.34 L / min. The spray droplet size range was 50-250 μm, and the dispersion coefficient under the Rosin-Rammler distribution was 2.23.

[0074] S102: Using the first simulation parameter information, respectively construct a rotor sliding grid model, a dynamic parameter porous medium canopy model, and a canopy droplet probability capture model.

[0075] The rotor sliding mesh model is constructed by using the inverse modeling method to build a three-dimensional rotor model, and the simulated flow field data at different speeds are compared with the measured data to determine the standard rotor speed. The rotor sliding mesh model is obtained by using the sliding mesh technology. The rotor sliding mesh model construction process is as follows: Figure 2 shown.

[0076] The OKIO 5M 3D scanner was used to obtain the point cloud data of the UAV rotor. After noise reduction, packaging, and repair in Geomagic Studio 12, the rotor 3D model was obtained. Figure 3 shown.

[0077] Field experiments measured the vertical component of wind speed at a horizontal section 65 cm below the UAV rotor at flight speeds of 3 m / s and 5 m / s. Simultaneously, simulation experiments were conducted based on the aforementioned 3D rotor model, obtaining eight sets of simulated wind speed data at rotor speeds of 160 rad / s, 180 rad / s, 200 rad / s, and 240 rad / s, under the same flight parameters as the measured values. Finally, by comparing the actual wind speed data at the same flight speed during actual operations, the rotor speed of 200 rad / s was determined to best match the actual wind field distribution. Rotation and translation of the rotor model were achieved using sliding mesh technology, ultimately determining the rotor sliding mesh model.

[0078] To construct a dynamic parameter porous media canopy model, first determine the basic equation for the canopy outlet airflow velocity. Then, use the trust region-based least squares method to identify the unknown parameters in the canopy outlet airflow velocity equation. Finally, load the functional relationship into the pre-set porous media model in the simulation software to obtain the dynamic parameter porous media canopy model.

[0079] Specifically, the canopy inlet airflow velocity, leaf area density, inertial drag coefficient, and porosity are used as input features to derive the canopy outlet airflow velocity calculation equation:

[0080]

[0081] in: V is the predicted value of the airflow velocity at the canopy outlet; is the airflow velocity at the canopy entrance; is the canopy inertial drag coefficient; is the canopy porosity; k is the number of iterations.

[0082] The residuals of wind speed data of two crop canopies were collected through computational simulation and field experiments. The functional relationship between the inertial drag coefficient and porosity change was constructed using the least squares method based on the trust region, and the undetermined parameters in the canopy outlet airflow velocity equation were identified. The objective function of the least squares method is as follows:

[0083]

[0084] The iterative calculation of the trust region algorithm is shown as follows:

[0085]

[0086]

[0087] in: is the parameter vector to be identified, is the optimal identification parameter vector obtained using the least squares method; is the simulated value of the canopy outlet airflow velocity; is the parameter vector at the kth iteration step for the least square method iterative calculation; is the Jacobian matrix of the equation of the canopy outlet airflow velocity m is the number of experimental data sets; is the gradient vector of the mth error function at the parameter vector ; and respectively represent the 1st and 2nd parameters in the vector ;

[0088] The function relationship between the above inertia resistance coefficient and the change of the canopy porosity is loaded into the preset porous medium model of the simulation software to obtain the dynamic parameter porous medium canopy model.

[0089] The canopy droplet deposition model is constructed, the droplet deposition behavior in the internal space of the canopy is divided into two processes of droplet blade contact and adhesion, the deposition probability formula of the droplets in the canopy in each grid unit is obtained, and the dynamic parameter coupled canopy droplet probability capture model is constructed.

[0090] Specifically, in the embodiment, the fixed parameters of the blade unit diameter, the leaf area density, the droplet density and the air dynamic viscosity are determined according to the actual working conditions in the field.

[0091] The droplet particle size, the droplet inertia and the airflow velocity of the grid center point of each grid unit are obtained, the contact probability and the adhesion probability between the canopy grid units and the droplets are calculated in real time, and finally the capture probability is calculated:

[0092]

[0093]

[0094]

[0095]

[0096] wherein: is the capture probability; is the contact probability of the droplet and the blade, is the droplet velocity; C 1 is an empirical coefficient of the model; is the time variation; is the adhesion probability of the droplet and the blade; is the droplet density; d is the droplet diameter; u is the air flow rate; is the air dynamic viscosity; D leafrepresents the leaf width of the vegetation.

[0097] A crown layer droplet deposition model of internal space deposition probability difference distribution is constructed based on the deposition probability calculation formula. In this embodiment, the crown layer droplet capture algorithm is written by using c or c++ language according to the deposition probability calculation formula, and the crown layer droplet capture algorithm is used to simulate, predict or optimize the calculation model of the deposition behavior of droplets in the crop canopy in the pesticide spraying process. The droplet capture algorithm is introduced into the discrete phase and porous medium module through the user-defined function module in the simulation software, the standard function of the simulation software is enhanced, and thus the crown layer droplet probability capture model of internal space deposition probability difference distribution is constructed.

[0098] A crown layer droplet deposition model of internal space deposition probability difference distribution is constructed based on the deposition probability calculation formula.

[0099] S103, a three-dimensional model of a calculation domain is constructed by combining the rotor slip grid model and the dynamic parameter porous medium canopy model, and grid division is performed.

[0100] Based on the rotor slip grid model constructed in S102, the fuselage is ignored, and a 6-rotor equivalent is used to represent the plant protection unmanned aerial vehicle, as shown in Figure 4 (a). Based on the crop phenotype parameters, the crown layer wind resistance properties and the droplet capture probability, the soybean and corn canopy are equivalently represented by a partitioned differential porous medium. The three-dimensional model is imported into Mesh, and a multi-level grid division strategy is used to appropriately encrypt the grid of the rotor, the spray head, the crop and the surrounding area, as shown in Figure 4 (b).

[0101] S104, the second simulation parameter information of the strip-shaped composite canopy simulation spraying is set.

[0102] DPM setting, DPM in the Lagrangian framework is used to simulate the droplet motion trajectory, and the spraying atomization parameters are set according to the actual working condition. The DPM particle injection source is arranged below the rotor model according to the actual spray head distribution pattern, and the injection source is controlled to follow the rotor area for synchronous horizontal displacement through the slip grid technology.

[0103] Boundary condition setting of calculation domain, the boundaries around the crops and near the unmanned aerial vehicle are set as pressure outlets with a pressure of 0, and the droplets are allowed to escape freely. The ground boundary is set as a non-slip wall that can capture droplets.

[0104] S105, the second simulation parameter information, the three-dimensional model of the calculation domain after grid division and the crown layer droplet probability capture model are combined to simulate, count and analyze the downwash flow field distribution and droplet deposition law under the strip-shaped composite irregular and heterogeneous canopy.

[0105] The droplet deposition visualization is performed by importing the second simulation parameter information, the meshed three-dimensional computational domain model, and the canopy droplet probability capture model into the post-processing software CFD-post, and introducing the droplet deposition amount and deposition density variables. The post-processing software CFD-post is then used to statistically analyze the distribution characteristics of the irregular heterogeneous canopy downwash flow field and the droplet drift and deposition characteristics.

[0106] The simulated UAV downwash flow pattern distribution and canopy airflow velocity distribution are as follows: Figure 5 and Figure 6 As shown in Figure 1, the downwash airflow from the drone has a strong penetration effect on the soybean canopy and the upper part of the corn canopy (2.0-2.9 m). However, the penetration ability of the middle and lower parts of the corn canopy (0-2.0 m) is significantly reduced, which will make it difficult for a large number of droplets to reach the bottom of the canopy, which is extremely detrimental to the uniformity of droplet deposition. The simulated drone spray droplet distribution and canopy deposition are shown in Figure 1. Figure 7 and Figure 8 As shown, droplet deposition in the corn canopy decreased significantly with decreasing sampling altitude, with droplet deposition in both the soybean and lower corn canopies (0-1.0 m) remaining below 2.0 μL / cm². In summary, using this combination of flight parameters for drone spraying resulted in poor droplet deposition distribution in both crop canopies, making it difficult to achieve the desired pest and disease control results.

[0107] To optimize drone spraying strategies, this example also conducted simulation experiments using various typical operational parameter combinations, comparing and analyzing the effects of varying flight altitudes and speeds on airflow and droplet penetration across the canopies of two crops. By revealing the downwash flow distribution and droplet deposition patterns under the irregular, heterogeneous, and strip-like canopies of soybeans and corn, this study provides a reference for the precise optimization of operational parameters such as flight altitude and speed during drone spraying operations.

[0108] In the embodiments of the present application, "at least one" refers to one or more, and "more" refers to two or more. "And / or" describes the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B can represent the existence of A alone, the existence of A and B at the same time, and the existence of B alone. Among them, A and B can be singular or plural. The character " / " generally indicates that the previous and next associated objects are in an "or" relationship. "At least one of the following" and similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one of a, b and c can be represented by: a, b, c, ab, ac, bc, or abc, where a, b, c can be single or multiple.

[0109] The above description is merely a specific embodiment of the present application. Any person skilled in the art may easily conceive of variations or substitutions within the technical scope disclosed in this application, and such variations or substitutions shall be within the scope of protection of this application. The scope of protection of this application shall be subject to the scope of protection of the claims.

Claims

1. A method for predicting droplet deposition and drift during aerial spraying of composite canopies, characterized in that: include: Obtain crop canopy and drone parameters from the field experiment as first simulation parameter information; The rotor sliding grid model, the dynamic parameter porous media canopy model, and the canopy droplet probability capture model are constructed using the first simulation parameter information, including: By combining the inverse modeling method and comparing the simulated flow field data at different speeds with the measured data, the standard rotor speed is determined and the rotor sliding grid model is constructed; Determine the basic equation of the airflow velocity at the canopy outlet and determine the undetermined parameter function relationship of the basic equation of the airflow velocity at the canopy outlet, and obtain the dynamic parameter porous medium canopy model in combination with the porous medium model preset in the simulation software; The droplet deposition behavior in the canopy space is divided into two processes: droplet contact and adhesion. The deposition probability formula of the droplet canopy in each grid unit is obtained, and a dynamic parameter-coupled canopy droplet probability capture model is constructed. Combining the rotor sliding grid model and the dynamic parameter porous media canopy model to construct a three-dimensional model of the computational domain and perform grid division; Setting the second simulation parameter information of the strip composite canopy crop simulation spray; Combined with the second simulation parameter information, the three-dimensional model of the computational domain after grid division and the canopy droplet probability capture model, simulation statistics are performed to analyze the downwash flow field distribution and droplet deposition law under the strip-shaped composite irregular and heterogeneous canopy.

2. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The step of obtaining crop canopy and drone parameters of the field experiment as first simulation parameter information includes: Through the aerial spraying experiment, the wind field data under the rotor under standard load was obtained, as well as the wind speed data of the airflow into and out of the two crop canopies at different flight altitude and flight speed combinations; The three-dimensional size of the equivalent porous medium canopy model of the crop is determined by field measurements of the crop strip composite planting pattern, crop plant height, and plant leaf area parameters, and the canopy leaf area density is calculated: in: LAD is the canopy leaf area density, S leaf represents the total leaf area of ​​the plant, S Indicates the area occupied by the plant, h Indicates plant height; Determine the spray atomization parameters based on the actual drone nozzle model and operating flow rate.

3. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The inverse modeling method is combined with comparing the simulated flow field data at different speeds with the measured data to determine the standard rotor speed and construct the rotor sliding grid model, including: Use an optical scanner to collect rotor point cloud data, and import the obtained point cloud data into reverse modeling software to reconstruct the rotor three-dimensional model; The rotor three-dimensional model is used to carry out a rotor downwash flow field simulation experiment through sliding grid technology. The rotor speed is determined by comparing the measured flow field data under the same operating conditions with the simulated flow field data under different speeds, and the rotor sliding grid model is obtained.

4. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The method of determining the basic equation of airflow velocity at the canopy outlet and determining the undetermined parameter function relationship of the basic equation of airflow velocity at the canopy outlet, and obtaining the dynamic parameter porous medium canopy model in combination with the porous medium model preset in the simulation software, includes: Taking the canopy inlet airflow velocity, leaf area density, inertial drag coefficient and porosity as input features, the calculation equation for the canopy outlet airflow velocity is derived: in: V is the predicted value of the airflow velocity at the canopy outlet; is the airflow velocity at the canopy entrance; is the canopy inertial drag coefficient; is the canopy porosity; k is the number of iterations; The residuals of wind speed data of two crop canopies were collected through computational simulation and field experiments. The functional relationship between the inertial drag coefficient and porosity change was constructed using the least squares method based on the trust region, and the undetermined parameters in the canopy outlet airflow velocity equation were identified. The objective function of the least squares method is as follows: The iterative calculation of the trust region algorithm is shown as follows: in: is the parameter vector to be identified, is the optimal identification parameter vector obtained using the least squares method; is the simulated measurement value of the airflow velocity at the canopy outlet; is the parameter vector at the kth iteration when the least squares method is iteratively calculated; The equation for the airflow velocity at the canopy outlet is Jacobian matrix; m is the number of experimental data groups; is the mth error function in the parameter vector The gradient vector at ; and Respectively The first and second parameters in the vector; The functional relationship between the above-mentioned inertial resistance coefficient and the canopy porosity change is loaded into the porous medium model preset in the simulation software to obtain the dynamic parameter porous medium canopy model.

5. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The droplet deposition behavior in the canopy interior is divided into two processes: contact and adhesion of droplets to leaves. The deposition probability formula of the droplet canopy in each grid unit is obtained to construct a canopy droplet probability capture model coupled with dynamic parameters, including: According to the actual working conditions in the field, the fixed parameters of blade unit diameter, leaf area density, droplet density and air dynamic viscosity are determined; Obtain the droplet size, velocity inertia, and airflow velocity at the center of each grid cell, calculate the contact probability and adhesion probability between each grid cell in the canopy and the droplet in real time, and finally calculate the capture probability: in: is the capture probability; is the contact probability between the droplet and the leaf, is the droplet velocity; C 1 is the model empirical coefficient; is the time variation; is the probability of droplets adhering to leaves; is the droplet density; d is the droplet diameter; u is the air velocity; is the dynamic viscosity of air; D leaf Indicates the width of vegetation leaves; A canopy droplet probability capture model with differential distribution of internal spatial deposition probability is constructed based on the deposition probability calculation formula.

6. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The step of constructing a three-dimensional computational domain model by combining the rotor sliding grid model and the dynamic parameter porous media canopy model and performing grid division includes: Ignoring the fuselage of the rotor sliding mesh model, a multi-rotor equivalent representation of the plant protection UAV is adopted; Based on crop phenotypic parameters, canopy wind resistance properties and droplet capture probability, a zoned differentiated porous medium equivalent method was used to characterize the crop canopy. A multi-level grid division strategy is adopted to appropriately encrypt the grids of the rotor, sprinkler, crops and their surrounding areas.

7. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The second simulation parameter information for setting the strip-shaped composite canopy crop simulation spray includes: The DPM under the Lagrangian framework is used to simulate the droplet motion trajectory and the spray atomization parameters are set according to the actual working conditions; The boundary pressure around the crops and near the drone is set to allow the droplets to escape freely, and the ground boundary is set as a no-slip wall that can capture the droplets.

8. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 1, characterized in that: The method of combining the second simulation parameter information, the three-dimensional model of the computational domain after grid division, and the canopy droplet probability capture model to perform simulation statistics and analyze the downwash flow field distribution and droplet deposition law under the strip-shaped composite irregular and heterogeneous canopy, includes: Importing the second simulation parameter information, the three-dimensional model of the computational domain after grid division, and the canopy droplet probability capture model into the post-processing software CFD-post and introducing the droplet deposition amount and deposition density variables; The CFD-post simulation is used to statistically analyze the distribution characteristics of the irregular heterogeneous canopy downwash flow field and the droplet drift and deposition characteristics.

9. The method for predicting droplet deposition and drift during aerial pesticide application in a strip-shaped composite canopy according to claim 8, characterized in that: Also includes: Through simulation experiments and analysis under various typical operating conditions, the downwash flow field distribution and droplet deposition rules under the crop strip-shaped composite irregular and heterogeneous canopy are revealed, and the drone flight control strategy is optimized.

Images