农药喷施是减少病虫害、提高作物产量的关键手段^[1]，农药喷施过程受到诸如喷头型号、药液喷施量、液滴频谱、喷雾器上喷头姿态等的影响^[2-3]，不同的作业方式和作业环境往往导致不同的作业效果，不适宜的作业环境与作业方式会使药液飘失并产生药害。作为减少农药用量、降低农药残留和污染的途径之一，农业航空施药技术近几年在我国迅速发展，并受到各个相关领域的广泛关注^[4]。低容量或超低容量喷雾是农业航空施药技术的主要喷雾方式^[5]，由于喷施雾滴粒径小，易产生飘移，探明环境风速对喷施雾滴雾化特性及飘移的影响规律是施药技术的关键。野外实地试验^[6-7]、风洞试验^[8]、仿真模拟^[9]是测量喷雾沉积飘移的常用手段。野外实地试验难以控制风速、可重复性差、成本高昂。相较于野外实地试验，风洞试验和仿真模拟能够模拟真实风速、风向等气象环境，准确控制相关参数，试验重复性好、可靠性高，已成为研究航空喷施雾滴雾化及飘移特性的主流方法^[10-12]。风洞试验可以根据试验需求对风速、风向等参数进行精准控制，从而模拟无人机作业时的真实飞行环境，国内外许多学者利用风洞模拟自然风，对不同参数下喷头喷施雾滴参数特性进行了相关试验。Teske等^[10]以旋转雾化器为研究对象，利用澳大利亚昆士兰大学的开路式风洞对雾滴雾化效果进行了研究，并建立了相关数据模型。Kirk^[13]通过改变不同管道压力、风速和喷头安装角度，对CP系列喷头进行了雾滴粒径拟合，预测不同参数下喷头雾滴粒径大小。Fritz等^[14]和Hoffmann等^[15]基于美国农业部农业航空研究中心的航空施药风洞，研究了高低速条件下的雾滴粒径分布规律；Martin等^[16]利用该风洞对静电喷雾在不同风速条件下的雾滴分布规律进行了研究。曾爱军等^[17]针对常用扇形喷头雾化特性进行测试，通过对飘移潜在指数的计算，发现雾滴粒径和风速是影响雾滴飘移的主要因素。张慧春等^[18]通过在昆士兰大学的风洞试验，研究了风速、喷头结构型号、药剂和采样距离对雾滴飘失的影响，建立了包含这4个因素的多元非线性雾滴飘移特性模型，提供了判断喷头雾谱等级的量化标准。唐青等^[19]基于北京农业智能装备技术研究中心自主设计制造的IEA-I型高速风洞，采用马尔文Spraytec雾滴粒度仪对标准扇形喷头和空气诱导喷头在高速气流条件下的雾化特性进行了研究，发现随着风速的增加，雾滴体积中径逐渐减少；管道压力的变化对标准扇形喷头的雾滴体积中径影响较大。

在仿真模拟方面，随着流体力学的不断发展和计算机软硬件性能的提升，相关学者基于仿真模拟方法、针对喷施雾滴雾化和分布特性也开展了一系列的研究工作。Dekeyser等^[2]通过计算流体力学 (Computational fluid dynamics，CFD)模拟技术，对几种常用喷头的羽流分布和附带空气流进行研究，证明了液体分布与产生的空气流动直接相关。Duga等^[20]考虑树木结构、冠层风流和喷雾器的运动，通过建立三维CFD模型，评估计算果园喷雾器喷雾沉积和飘移，并对不同喷嘴布置的苹果园进行了试验验证。Hong等^[21]开发了一个综合CFD模型，预测了空气辅助喷雾器在树冠内部和周围的风速分布，并与实际测试结果对比，结果表明该模型可以合理预测空气辅助喷雾器排出的空气分布。孙国祥等^[22]基于三维模拟离散相技术建立了不同喷雾高度和风速下的雾滴沉积量和沉积率预测模型。刘雪美等^[23]利用三维流场模拟多相流模型研究了在自然风影响、辅助气幕胁迫和自身重力作用下，连续相和雾滴粒子群离散相耦合的交互作用对雾滴沉积的影响。风洞试验和仿真模拟是研究喷头雾化及雾滴沉积飘移特性的重要手段，鲜见有人将风洞试验与仿真模拟结合起来研究航空喷施作业药液雾滴雾化和沉积过程。本文以航空喷施作业常用的Lechler系列LU 120-015和LU 120-03标准扇形压力式喷头为研究对象，通过仿真模拟将喷头喷施雾滴分布特性可视化，探究雾滴的沉积和飘移分布特性，结合实际的风洞试验验证模拟结果的准确性，并对其喷施过程的雾滴飘移特性进行分析，以期为航空植保喷施实际作业过程中作业参数的选择和关键部件改进提供数据支撑和理论指导。

1 模型构建及模拟 1.1 几何模型建立及网格划分

如图1所示，在Ansys15.0 Workbench下集成三维建模软件Geometry中建立模拟计算域模型，模拟计算域设为长20.0 m (x：−2.5~17.5 m)、入口处宽2.0 m (z：−1.0~1.0 m)、出口处宽2.2 m (z：−1.1~1.1 m)、高1.1 m (y:0~1.1 m)的箱体。为方便后期统计，在计算域底面以原点(0，0，0)为圆心建立一个狭长椭圆面A_face，用于后期统计雾滴准确沉积，沉积面在距喷雾口高度为l的平行面上，A_face是以喷嘴轴心为圆心、a为半长轴、b为半短轴的椭圆。计算公式如下：

$ a=l \times {\rm{tan}}\left( {{\theta _1} + {\theta _2}} \right){\rm{,}} $

(1)

$ b=l \times {\rm{tan}}\left( {{\rm{arctan}}\frac{{{l_{\rm{a}}}}}{{{l_{\rm{b}}}}} + {\theta _2}} \right){\rm{,}} $

(2)

$ {A_{{\rm{face}}}={\text{π}} ab}{\rm{,}} $

(3)

式中，a为准确沉积区域的半长轴、b为准确沉积区域的半短轴，m；A_face为雾滴准确沉积区域，m²；θ₁为喷射半角，(°)；θ₂为喷雾扩散角，(°)；l为沉积面高度，此处为0.6 m；l_a为喷口半短轴，此处为0.095 mm；l_b为虚拟原点到喷口距离，此处为1.200 mm。同时，在计算域内x=2~15 m处以1 m等间距建立多个面，用以切割上述雾滴沉积面，用Meshing进行四面体网格划分，为保证后期结果准确性，划分过程优先考虑网格质量，网络划分单元总数为300 163个。

1.2 数值计算模型 1.2.1 连续相模型选择

本文采用基于压力的瞬态模拟求解，考虑重力作用。为精确模拟和捕捉空气流场细节，根据喷雾气流特征，连续相湍流模型采用Launder和Spalding提出的标准k-ε模型^[24]，该模型在保证雷诺应力求解计算精度的条件下具有较好的稳定性和经济性，近壁面处理方法为标准壁面功能。

1.2.2 离散相模型选择

雾滴输送方程采用欧拉−拉格朗日方法^[25]求解，其离散相颗粒运动方程如下：

$ \frac{{{\rm{d}}{{\vec u}_p}}}{{{\rm{d}}t}}=\frac{{18\mu }}{{{\rho _{\rm{p}}}d}}\frac{{{{{C}}_{{D}}}{{Re}}}}{{24}}\left( {\vec u - {{\vec u}_{\rm{p}}}} \right) + \frac{{{{g}}\left( {{\rho _{\rm{p}}} - \rho } \right)}}{{{\rho _{\rm{p}}}}} + \frac{{\rho d}}{{2{\rho _{\rm{p}}}{\rm{d}}t}}\left( {\vec u - {{\vec u}_{\rm{p}}}} \right), $

(4)

式中， $ \overrightarrow{u} $ 为连续相速度，m/s； $ \overrightarrow{{u}_{\rm{p}}} $ 为颗粒速度，m/s；ρ_p为颗粒密度，kg/m³；ρ为空气密度，kg/m³；d为颗粒直径，m；g为重力加速度，m/s²；μ为黏性系数；Re为相对雷诺数；C_D为曳力系数。

喷嘴模型采用平面扇形喷嘴模型(Flat-fan atomizer model)，雾化破碎模型采用线性化不稳定液膜雾化模型。喷雾模拟考虑离散相雾滴碰撞和聚合，判断雾滴合并和反弹主要根据Rourke等^[26]得到的临界值，临界值是碰撞韦伯数和集合雾滴管与小液滴半径的函数，计算公式如下：

$ {b_{{\rm{crit}}}}=\left( {{r_1} + {r_2}} \right)\sqrt {{\rm{min}}\left( {1.0,\frac{{2.4f}}{W}} \right)}, $

(5)

式中，b_crit为判断雾滴碰撞合并或反弹的临界值，m；r₁、r₂为小液滴半径，m；f为r₁/r₂的函数；W为碰撞韦伯数，雾滴碰撞二次破碎模型选用Taylor破碎模型(TAB)。

1.2.3 模拟试验参数设置

1) 离散相喷射源参数：离散相材料为水，离散相释放位置坐标为(x，y，z)=(0，0.6，0)；喷嘴轴向矢量分量为(x，y，z)=(0，−1，0)；雾滴质量流率为0.02 kg/s；喷射半角为60°；喷雾扩散角为6°。启用非稳态粒子追踪，采用离散随机游走模型，粒子释放时间尺度常数为0.01 s。

2) 边界条件参数：对于流体边界条件(Boundary condition)：箱体域x=−2.5 m平面为速度入口边界(Velocity-inlet)，速度分别设置为0、1、2、3、4、5和6 m/s；x=17.5 m平面为自由出流边界(Outflow)；其他面均为壁面边界(Wall)。对于离散相的边界条件类型(DPM boundary condition type)，入口处和出口处为逃逸(Escape)，地面为形成液膜(Wall-film)，其余均为反弹(Reflect)，计算区域内部材料为空气。

3) 求解方法及模拟参数：压力−速度耦合方式为SIMPLE，瞬态公式为二阶隐式，亚松弛因子保持默认设置。求解开始前对全局进行标准初始化，求解时迭代时间步长为0.010 s。当研究风速对雾滴粒径影响时，模拟喷施时长0.005 s，喷施开始前开启风场，待模拟区域内流场稳定后开始喷雾，喷雾结束后延迟0.010 s统计流场域中的全部粒子粒径分布，目的是让喷施的离散相粒子可以在空中充分碰撞或聚合，此时离散相粒子距离喷头平均长度为0.35 m，达到常规的雾滴统计位置；当研究雾滴沉积飘移特性时，模拟喷雾时长5 s，同样，喷雾开始前开启风场，待模拟区域内流场稳定后开始喷雾，喷雾结束后风场保持开启一段时间，直到空间内雾滴被地面捕捉或离开计算域。

2 风洞试验 2.1 试验设备和材料

本试验所用风洞位于华南农业大学国家精准农业航空中心风洞实验室，为符合ISO国际标准(ISO 22856)的高低速复合风洞。试验过程中为消除因壁面湍流和雾滴反弹造成的试验误差，风洞底面铺设防飞溅人工草皮。风速测量设备为加野麦克斯公司生产的Kanomax 6036-BG Anemomaster带压力传感器式数字风速计。喷雾系统是由农业农村部南京农业机械化研究所设计的喷雾控制系统，该系统由延时继电器、储水箱、增压泵、卸压阀、压力表、喷施管道以及喷头(德国Lechler系列LU 120-015和LU 120-03标准扇形压力喷头)组成，通过调节减压阀出口压力，能够精确控制喷施压力，通过调节继电器模式，可以定时喷施。荧光度检测设备为天津港东科技发展有限股份公司生产的F-380荧光分光光度计及配套软件。

2.2 风洞试验设计

将喷头安装在风洞内距离风洞底面0.6 m的中心支架上，调整喷嘴方向，使喷嘴竖直朝下，喷施扇面椭圆长轴方向和来风方向垂直。试验介质选用可溶性荧光示踪剂(质量分数为0.5%的罗丹明−B)水溶液，采用直径为1 mm的聚乙烯线进行雾滴收集。在下风向距离喷头平面2 m的位置，由风洞地面向上0.1~0.8 m放置8根间距为0.1 m的收集线，用以检测穿过垂直平面空气的雾滴通量，分别命名为V1~V8；其中，设置距离地面为0.1 m收集线的目的是消除由于风的湍流和雾滴在地面的飞溅导致对收集线的污染。此外，沿水平方向在距离地面0.1 m高的位置以1.0 m的间隔距离放置13根收集线，以检测喷雾从2~15 m范围内的水平飘失，分别命名为H3~H15。根据植保无人机的航空植保喷施作业规范，环境侧风大于3级风时不能进行喷施作业，选择1、3、6 m/s的风速进行风洞试验，每次喷施时间为5 s，喷施结束后待附着在线上的雾滴充分晾干，戴上一次性橡胶手套收集聚乙烯线并单独放置在编好号的自封袋中，避光低温保存，及时处理，每组试验重复3次。风洞试验布置如图2所示。

2.3 评价指标及数据处理

通过Fluent自带离散相粒子统计功能对沉积高度为0.6 m的地面进行统计，计算获得雾滴飘移特性。通过对沉积在底面上不同单位距离长方形块内的离散相质量统计，可以获得不同侧风条件下离散相飘移情况。将收集的聚乙烯线用40 mL去离子水充分振荡洗涤，取样通过荧光分光光度计测量得到示踪剂的质量浓度(μg/L)，质量浓度乘以洗脱液体积(L)，得到采集线上沉积的示踪剂质量(μg)，通过竖直排列布置的8根水平收集线上获得的数据，可以拟合得到测试垂直面内雾滴飘移体积通量随高度分布的n次多项式曲线 $ \dot{v}\left(y\right) $ ，同时得到飘移量分布的特征高度(h，m)^[27]，其定义为：

$ h = \frac{{\displaystyle\int _0^{{{{h}}_{\rm{N}}}}{\dot{{v}}}\left( {{y}} \right){{y{\rm{d}}y}}}}{{\displaystyle\int _0^{{{{h}}_{\rm{N}}}}{\dot{{v}}}\left( {{y}} \right){{{\rm{d}}y}}}}{\rm{,}} $

(6)

式中， $ \dot{v}\left(y\right) $ 为通过测试截面内任意高度雾滴的体积通量，L/s；y为测试垂直面的单个计算高度，m；h_N表示测试垂直面的累积计算高度，m；拟合准确度均超过97%。

雾滴水平飘移率(R_h)以距离地面0.1 m布置的13根聚乙烯线收集的雾滴沉积量占喷头喷施雾滴的百分比来表示，计算公式如下：

$ {R_{\rm{h}}} = \frac{{{A_{\rm{d}}}}}{{{T_{\rm{a}}}}} \times 100{\text{%}} {\rm{,}} $

(7)

$ {A_{\rm{d}}} = \sum\limits_{{{i}} = 1}^{{n}} {{{{d}}_{{i}}}} \left( {\frac{{{s}}}{{{w}}}} \right){\rm{,}} $

(8)

$ {T_{\rm{a}}} = vc{\rm{,}} $

(9)

式中，A_d为喷雾过程中，雾滴在收集线上的沉积总量，μg；n为收集线的数量；d_i为第i根收集线上示踪剂的沉积，μg；s为收集线间的距离，m；w为收集线的直径，本试验中的w为0.001 m；T_a为喷施示踪剂的质量，μg；v为喷雾体积，L；c为示踪剂质量浓度，μg/L。

3 结果与分析 3.1 数值模拟结果

无侧向风条件下，喷头的雾滴在风洞内部的沉积分布模拟情况如图3所示。由图3可以看出，无侧向风时，喷嘴离散相雾滴以喷射半角为60°的扇形面分布，并以6°的扩散角向外扩散，除空间中的极少数无序离散相粒子外，绝大部分离散相粒子均在式(1)~(3)所求得的椭圆范围内。

图4为无侧风条件下雾滴在椭圆面内准确沉积比例(R_a)和水平飘移率(R_h)的沉积分布，不同颜色代表不同的离散相沉积质量浓度，蓝色向红色渐变代表浓度逐渐加大。由图4可知，在无侧向风时，计算机模拟所得离散相雾滴在底面上的沉积质量浓度云图与通过式(1)~(3)计算所得区域基本一致，均沉积在以喷嘴为中心的椭圆区域。由于喷嘴正下方中间区域直接喷射且喷口中间较宽，所以离散相粒子能大量沉积；底面中心向两侧的沉积量逐渐下降，主要是由于喷嘴喷口两侧较窄，且喷施雾滴在z轴方向的运动势能随喷施时间的增加而减小，因此，雾滴粒子在喷施范围内雾量沿z方向总体上呈先增大后减小的正态分布。另外，由沉积比例图可以看出，风洞壁两侧处也出现较大浓度的雾滴沉积，这主要是由于单喷头喷幅两侧的部分雾滴在运动势能减小的情况下向两侧发生了扩散，并与风洞侧壁碰撞反弹至底部从而产生一定的雾滴沉积。

为更好地反映出不同侧风风速参数下的离散相雾滴分布特性，选择适当的视角以反映离散相雾滴沉积质量浓度区间。图5为侧风风速为0~6 m/s时雾滴在计算区域内R_a和R_h的沉积分布图。由图5可以看出，随着侧风风速的增加，雾滴水平飘移越来越明显，高浓度沉积区域逐渐偏离喷雾中心轴线(x=0 m)朝下风向偏移，沉积质量浓度大于0.008 kg/m³的x轴区域边界由5 m逐渐推移至13 m。对于沉积在底面上的雾滴，x轴向沉积质量浓度分布变化明显，z轴向雾滴沉积质量浓度分布接近对称，从z轴中心轴线向两侧呈现出由低到高再到低的“M”型分布。对于下风向的离散相沉积质量浓度，底面与两壁交接处有2段沉积质量浓度很高的窄条，在距两侧壁面0.1 m处存在一块沉积质量浓度较低的区域，且随着风速增大，区域沿x轴方向逐渐拉长。

图6为在x=1 m平面内气流的z方向分速度云图和风洞内x轴方向的流线图。由图6的云图可见，雾滴粒子由蓝向红运动的方向表示z方向分速度由负变为正；雾滴粒子在x=1 m处的分速度近似对称分布，表明离散相粒子在x=1平面及其附近空间内存在着涡旋运动。图6的流线图可以看出，从入口处进入的平行气流通过扇形喷雾面后形成2个关于xy平面对称的涡旋，涡旋的影响范围超出喷雾的初始位置，导致雾滴在一定空间范围内出现上扬；在此处两侧下风向的小雾滴很容易被卷到喷施扇面后方，导致小雾滴均沉积在靠近风洞中心平面处，此时离散相呈现出非流体特性，说明离散相雾滴和连续相空气相互耦合且彼此影响。

为进一步表明侧向风速对雾滴飘移的影响，用狭长椭圆内雾滴沉积质量占总的离散相沉积质量的比例来表示某一风速下的R_a，用下风向2 m后地面和出风口捕获雾滴沉积量占总的沉积质量的比例来表示某一风速下的R_h，根据计算机模拟针对风洞内底面上不同区域的雾滴沉积质量的统计结果，计算雾滴在不同风速下的R_a和R_h，结果如图7所示。

由图7可以看出，随着侧向风速的增加，模拟离散相的R_a呈指数下降，由14.11%降到0.66%；R_h呈线性增加，由14.25%增加到60.58%。可见，随着侧向风速的增大，离散相雾滴飘移程度越来越严重。为进一步分析沉积结果和侧向风速的相关性，对R_a、R_h与风速(v)分别进行回归分析，回归方程分别为：R_a=0.1476e^−0.529v (R²=0.995)和R_h=0.0796v+0.1456(R²=0.995)，对应的显著性水平P值分别为0.0063和0.0036，均达极显著相关(P<0.01)，表明随着侧向风速的增加，离散相R_a呈指数式降低、R_h呈线性增加；侧向风速的增加不利于喷施雾滴在目标区域内的准确沉积，容易造成飘移。

3.2 风洞试验结果及数据对比分析

对LU 120-015和LU 120-03喷头在不同风速下的雾滴飘移进行风洞试验，喷施压力为0.3 MPa，侧向风速分别为1、3和6 m/s，结果如图8所示。由图8可知，在侧风风速为1 m/s时，喷头在下风向2 m处的竖直平面内雾滴通量由下至上逐渐减小；在侧风风速为3和6 m/s时，雾滴通量先增大后减小。可见，当侧向风速较小时，在喷头压力作用下，雾滴上扬程度低，主要向底部运动；当侧向风速增加时，较大的侧向风速会导致喷出的雾滴上扬，此时雾滴在竖直平面内的分布呈先增加后减少的趋势。

对于下风向2 m及大于2 m的雾滴沉积飘移量，在不同侧风风速条件下，随着收集线离喷头水平距离的增加，雾滴沉积量均呈现逐渐下降的趋势，但下降速率有一定的差异。当侧向风速为3 m/s时，随采样距离的增加，雾滴沉积量下降速率大于侧向风速为6 m/s时。相对于大雾滴，小雾滴沉降速度慢，随着风速逐渐增加，侧向风对雾滴水平飘移距离的影响也逐渐增大，小雾滴更容易受到侧向风的影响，且风速越大，雾滴的飘移距离越远。当侧向风速为3 m/s时，飘移的雾滴及数量非常小，沉积量下降速率较快；侧向风速为6 m/s时的沉积情况相反，且同一位置采集到的雾滴沉积量大于3 m/s时的雾滴沉积量。可见，侧向风速较大时雾滴飘移的数量及粒径增加，从而导致在同一位置处采集到的雾滴沉积量较大。通过计算得到3种侧风风速下雾滴飘移量分布的特征高度(h)分别为0.175、0.200和0.245 m，对应的水平飘移率分别为0.4%、48.1%和75.1%，这与雾滴飘移量随着风速的增大而逐渐增加的规律相吻合。

通过统计模拟结果中竖直单位距离内的离散相粒子沉积质量，获得模拟结果的飘移量数据。将风洞内水平布置的V1(H2)、H3~H15聚乙烯线上的罗丹明−B沉积质量与收集距离进行线性拟合，反推得出在一定距离下的雾滴飘移量数据。通过SPSS软件进行相关性分析可知，在1、3和6 m/s侧风风速下，风洞实测与数值模拟的雾滴飘移量的相关系数分别为0.905、0.995和0.978。对计算机模拟的雾滴水平飘移率(R_ht)和风洞试验中雾滴水平飘移率(R_hs)进行相关分析，其表达式为R_ht=1.888R_hs−0.3533 (R²=0.963)，对应的显著性水平P值为0.038 (P<0.05)，表明两者具有较强的相关性，从而验证了模拟结果的准确性。

4 结论

本文以航空植保无人机平面扇形喷头为研究对象，采用CFD离散相模型的粒子跟踪技术方法模拟雾滴在低风速条件下的沉积飘移特性，并结合相近情况下的风洞试验对扇形喷头的雾滴沉积飘移特性进行了测试和验证。研究结果表明：

1)通过对模型结果的统计和分析，随着侧向风速的增加，模拟离散相粒子的准确沉积率(R_a)由14.11%呈指数下降到0.66%，水平飘移率(R_h)由14.25%呈线性增加到60.58%，可见，随着侧向风速的增大，离散相雾滴飘移程度越严重，雾滴水平飘移越明显。

2)在不同的侧向风速下，雾滴在竖直平面内的分布有差异。计算得出1、3和6 m/s风速下雾滴的飘移量分布的特征高度(h)分别为0.175、0.200和0.245 m，水平飘移率分别为0.4%、48.1%和75.1%，可见，随侧向风速的增加，雾滴飘移程度急剧增加。

3)通过对比计算机模拟与风洞测试结果可知，设置合适的喷施参数，计算机模拟可以获得与风洞测试试验有一定相关性的雾滴飘移结果；对比计算机模拟及风洞测试试验的水平飘移率发现，两者具有显著的相关性(R²=0.963，P<0.05)，说明数值模拟和风洞测试试验结合对于雾滴飘移具有较好的预测效果。