不同雾化方式喷雾场的数值分析与试验验证.pdf

年 月 第 卷 第 期 刘秀娟不同雾化方式喷雾场的数值分析与试验验证刘秀娟，张慧春，（南京林业大学机械电子工程学院，江苏 南京 ；南京林业大学林业资源高效加工利用协同创新中心，江苏 南京 ）收稿日期：；修回日期：；网络出版时间：网络出版地址：基金项目：国家自然科学基金资助项目（）；江苏省现代农机装备与技术示范推广项目（）；江苏省重点研发计划项目（）；江苏省 工程项目（苏人 ）第一作者简介：刘秀娟（），女，吉林吉林人，讲师（通信作者，），主要从事流体机械及计算流体动力学研究第二作者简介：张慧春（），女，内蒙古呼和浩特人，教授（），主要从事农林信息技术与装备、智能植保机械和表型技术与平台研究摘要：为了研究植保雾化喷头的雾化效果，基于 中拉格朗日粒子跟踪气液两相流模型，分别结合压力雾化和空气雾化模型对不同风速下的雾滴粒径、运动速度、质量浓度分布规律进行了对比计算 计算结果表明：压力雾化模型中，雾滴粒径小于 的占比小于 ，粒径 的细雾滴占比约为 ，近喷雾口处雾滴速度更高，约是空气雾化模型的 倍；雾滴质量浓度分布更均匀；当风速达到 时，细雾滴分布比例增加，雾滴会出现二次破碎 空气雾化模型，中等雾滴粒径分布在 ，占比约为 ，雾滴速度衰减慢，喷雾口附近雾滴不易逆风运动；雾滴沉积范围小，喷幅较窄，大部分雾滴沉降在以喷头为中心的 范围内 此外，完成了低速风洞试验，收集雾滴沉积并进行卡方检验，结果表明 个雾化模型对于喷雾场中液滴沉积值均可信，但压力雾化模型计算值与试验值偏离度较小关键词：雾化方式；雾滴粒径；质量浓度分布；沉积；卡方检验中图分类号：文献标志码：文章编号：（）：刘秀娟，张慧春 不同雾化方式喷雾场的数值分析与试验验证 排灌机械工程学报，（）：，（），（）：（），（，；，）：，排灌机械工程学报第 卷 ，：；植物生长过程中的病害、虫害、杂草的控制与农药应用技术直接相关，而喷头的雾化方式直接影响着农药的应用效益 近些年病虫害精准防治技术越来越受到重视 ，为了研究喷头雾化效果，需要进行全面的实地现场测试 ，通常试验成本高昂，可重复性差，可操作性不强 因此，计算流体动力学（，）技术作为一种有效和成熟的工具，以其数值模拟的特点，解决了因试验条件限制难以解决的问题，得到了广泛应用目前，国内外一些学者已做了一些喷雾机工作时雾滴沉积分布或飘移 模拟研究 等 将 种不同的苹果树和梨树冠型集成到几何模型中，设计了 种空气辅助果园喷雾机，模拟果园气流对喷雾分布的影响 王景旭等 基于 模拟探讨了在温室环境中用气流辅助方式喷施农药时，施药对象周围的流场对雾滴飞行轨迹及雾滴附着行为产生的影响 等 将粒子束按直径离散成 的雾滴，以 分布规律从移动喷雾机口注入，作为 模型的雾滴源，研究苹果园的树冠上雾滴沉积、脱靶沉积和飘移 傅泽田等 建立了 型果园风送式喷雾机气流场速度分布仿真模型，并设计试验对模拟结果进行了验证为了研究植保雾化喷头的雾化效果，文中基于 技术，采用拉格朗日粒子跟踪气液两相流模型，分别结合压力雾化、空气雾化模型数值模拟压力雾化、空气雾化喷头，对比分析 种雾化的喷雾场特征 选用 公司的型号 压力喷头、空气雾化喷头进行试验验证，定量分析种雾化方式在受到不同横向气流干扰时，对雾滴粒径、沉积等性能的影响，以期为喷雾系统的前期设计提供参考 数值模型与计算方法 计算方法喷雾是气 液两相流相互作用的一个过程，根据喷雾流动特征，连续相模拟采用 湍流模型，其输运方程为（）（）()，（）（）（）()（），（）式中：为连续相密度，；为湍流动能，；为方向上的速度，；为湍流动能耗散率，；，分别为湍流动能、湍流动能耗散率的有效普朗特数的倒数，；为偏湍流黏度张量，；为平均速度梯度引起的湍流动能项，；为由于浮力影响引起的湍流动能项，；为可压缩湍流脉动对总的耗散率的影响，；，为自定义项，；为湍流动能耗散率引起的附加项，；，为经验常数，分别为 ，计算中考虑连续相与离散相雾滴间的双向耦合作用 在笛卡尔坐标系中离散相雾滴运动方程为 （）（），（）（），（）（），（），（）第 期刘秀娟，等不同雾化方式喷雾场的数值分析与试验验证 ，（）式中：为雾滴速度，；为气流速度，；为黏性系数，；为雾滴密度，；为雾滴直径，；为曳力系数；为雷诺数；为重力加速度，；为雾滴受到 方向的附加力，；为雾滴单位质量曳力，；，为量纲一常数常用植保雾化喷头的工作原理是单向压力雾化和低压空气雾化 压力雾化喷头在喷头内部将压力能转化为液体柱的动能，使其高速喷入静止或较低速的空气环境中，气液相对运动剧烈，液体表面张力及黏性减弱，在气动力作用下雾化成大量的细小液滴 空气雾化喷头在喷头单侧或双侧设置喷射泵或文丘里喉道，将雾化空气引入混合室，在混合室内，液气速度差引起气体与液体混合物破碎分解成雾滴，再经喷口喷出 文中在数值计算中，喷雾模型选用压力雾化模型和空气雾化模型 离散相材料为水，计算中考虑雾滴的碰撞和聚合，雾滴二次破碎选用泰勒比拟破碎模型；在计算雾滴运动轨迹时，采用离散相随机轨道模型描述气流的瞬时速度对雾滴运动轨迹的影响 几何模型建立及参数设置在前处理软件 下建立几何模型如图 所示，模型计算域长 宽 高为 ，网格划分总数为 喷雾口距进风口 ，距底面 ，喷雾方向与气流方向垂直，与重力方向一致 在进风方向上，距喷雾口 、底面 位置，轴方向自下向上每间隔 ，建立 个雾滴采样圆柱体；轴方向，距离喷雾口 ，和 建立 个雾滴采样圆柱体，圆柱体长 、直径 图 仿真几何模型 计算域左侧面边界条件为速度入口（进风口），速度值分别为 ，；计算域右侧面与大气相接，边界条件设为压力出口；其余面边界设为壁面喷雾参数为上游喷雾压力 ，喷雾半角 ，喷雾流量率 ，喷孔直径 ；喷射时间 ，喷雾扩散角 ，液膜常数 ，索带常数 试验设计 雾滴沉积试验为了验证 数值模拟结果的可靠性，设计了收集雾滴沉积试验装置，如图 所示 沉积试验中收集杆的位置与模拟计算采样圆柱体相同 试验选用压 力 雾 化 喷 头 和 空 气 雾 化 喷 头 ，在小型低速风洞进行喷雾试验 试验中喷头被固定在距进风口 、距地面 位置应用长 、直径 的碳纤维杆收集喷雾沉积量 采用质量浓度为 的荧光素钠溶液作为示踪剂 喷雾结束后，把收集杆装入管袋中，加入 清水充分振荡洗涤后，用荧光分光光度计（型号为 ）对荧光素钠质量进行测量，每组试验 次，取其平均值作为沉积量进行定量分析图 雾滴沉积试验布局 对雾滴沉积计算值和试验值进行归一化处理，即，（），（）式中：，为计算域 采样圆柱体上雾滴沉积量，；，为计算域内雾滴总质量，；，为试验收集杆沉积溶液中荧光素钠质量，；，为 喷施溶液中荧光素钠的总质量，将计算值和试验值进行卡方检验分析，即（，），（）式中：为试验次数如果由卡方检验式（）计算的卡方小于临界值 ，说明液滴沉积归一化计算值分布在显著水平 时与归一化试验值存在统计学关联 雾滴粒径试验为了研究风速对雾滴粒径的影响，设计图 所 排灌机械工程学报第 卷示试验装置 在小型低速风洞内，对喷头 和 进行喷雾压力 ，风速 ，和 的喷雾试验 试验中用相纸卡收集带有颜色的雾滴，在距离喷头 的下风处、距风洞底面高度 ，设置 个相纸卡采样点图 雾滴采集试验布局 试验结束后，使用 将相纸卡扫描成图片，应用雾滴分析软件 对图片进行处理，测出雾滴斑点面积，由式（）计算雾滴实际粒径 ，（）式中：为雾滴实际粒径，；为雾滴沉积在相纸卡上斑点的面积，试验结果与讨论 数值计算结果分析 雾滴粒径及运动速度图 为风速 时 种雾化模型下雾滴的粒径 及运动速度 由图 可见，大雾滴沉积在喷头附近，较小液滴悬浮在雾滴群上层，被风带到高处并顺风运动得更远 初始时，压力雾化模型雾滴粒径更小，顺风运动的雾滴随风飘移距离远达喷头 甚至更远 由图 可见，压力雾化模型雾滴近喷雾口处速度更高，约是空气雾化模型的 倍 空气雾化模型雾滴初始时与气流相对速度越大，受到外部气流阻力越大，因此在图 （标志区）显示喷雾口附近的雾滴不易逆风运动 空气雾化模型大部分雾滴沉降在距离喷头 范围内 说明液滴粒径和横向气流存在互作效应，共同影响着雾滴的运动图 不同雾化模型下的雾滴粒径及运动速度图 雾滴沉积质量浓度对比在风速 、距离喷雾口高度 ，和 时，种雾化模型雾滴质量浓度 分布，如图 所示图 不同雾化模型下距喷头不同高度时的雾滴质量浓度分布 第 期刘秀娟，等不同雾化方式喷雾场的数值分析与试验验证从图 中可以看出，雾化方式对雾滴质量浓度的空间分布存在较大的差异 压力雾化模型雾滴释放后质量浓度分布更均匀，质量浓度先增加后降低，且随着距离喷雾口水平距离增加，雾滴数量逐渐减少，雾滴更易受到风的影响，顺风飘移得更远、沉积覆盖表面积更大 空气雾化模型雾滴释放后质量浓度迅速达到峰值，雾滴主要沉积在喷雾口附近，外围质量浓度更高，且由于喷雾口处雾滴密集、粒径较粗，喷雾的后方受到了风的防护，横风难以穿透高度密集的液滴，导致雾滴沉积范围小，质量浓度迅速衰减，喷幅较窄 雾滴粒径分布对 种雾化方式的雾滴粒径数目分布 （某一尺寸液滴数目占液滴总数的数目分布）进行对比分析，如图 所示，为风速图 不同风速下 种雾化模型的雾滴粒径数目分布 由图 可以看出，压力雾化模型中，粒径小于 的雾滴占比小于 ，粒径 的细雾滴占比约为 空气雾化模型，中等雾滴粒径分布在 ，占比 左右 当风速由 增加到 时，种雾化方式对雾滴粒径影响均不明显 但是当横风速达到 时，雾滴粒径分布曲线有向左偏移的趋势，压力雾化模型细雾滴分布比例降低，极细雾滴分布比例增加 说明当风速达到 时，横风加剧了雾化场中空气与雾滴相互作用，导致部分粒径雾滴二次破碎，小雾滴比例增加 试验结果的对比验证 雾滴粒径对比不同风速 下，种雾化方式雾滴粒径 （雾滴体积分布占全部雾滴体积分布 的雾滴粒径）的计算值和试验值间的差异，结果见表 由表可见，风速的 和 时，风速对雾滴粒径影响不明显，但计算值比试验值结果偏大 风速 时，试验值 粒径比风速 和 时明显增大，这是因为风速为 时，雾滴顺风速度相对于降落速度大，撞击到收集雾滴的相纸后，部分雾滴反弹减速后再次沉积在相纸上，和原雾滴痕迹相连、合并，导致在相纸上留下的雾滴斑点变大，软件 数据处理时，不能准确测量雾滴粒径，导致雾滴粒径偏大表 风速对雾滴粒径的影响 （）压力雾化空气雾化试验值计算值试验值计算值 雾滴沉积通过对 种雾化方式计算值和试验值雾滴沉积对比，得到图 所示结果 图中 为距离喷雾口高度，为沿风向与喷雾口的距离，为雾滴沉积归一化值 从图中可以看出，风速对雾滴沉积的影响显著，随风速增加，喷雾场中同一位置处雾滴的沉积明显增大 风速 时，在距喷雾口高度上，压力雾化模型雾滴沉积计算值与 喷头试验值有较 好 的 一 致 性 空 气 雾 化 模 型 计 算 值 与 喷头试验值的变化趋势相同，但计算值偏大，部分原因可能是空气雾化模型模拟空气雾化喷头时气液相对速度不能准确预估对雾化模型数值结果进行可靠性评估，由式（）得到 ；查找卡方临界值表，当自由度为 ，概率为 的 （），故卡方检验表明，雾滴沉积归一化计算值与试验值存在统计学关联，说明文中所采用的数值计算喷雾方式对于流场中雾滴沉积结果的分析是可信的；，说明压力雾化模型比空气雾化模型计算值与试验值偏离度小 排灌机械工程学报第 卷图 雾滴沉积计算值与试验值的对比 结论雾化方式对雾滴粒径、运动速度、质量浓度等性能参数有显著的影响）数值计算结果表明：压力雾化下，雾滴初始沉降速度较大；细雾滴 粒径占比较大，当风速达到 时，细雾滴分布比例降低，极细雾滴分布比例增加 空气雾化下，中等雾滴粒径分布在 ，雾滴初始速度较低，不易逆风运动）在相同的工况参数下，种雾化方式中，计算值与试验结果在不同程度上存在定量一致性 压力雾化下，在风速 时，距离喷雾口不同高度，沉积计算值与试验值有较好的一致性 空气雾化下计算值比试验值偏差大一些，部分原因可能是空气雾化模型模拟空气雾化喷头时气液相对速度不能准确预估）风速由 增加到 时，对雾化模型与试验喷头 雾滴粒径影响不明显 当风速增加到 时，雾化模型 减小，试验喷头 结果显著增大参考文献（），：康灿，高可可，张永超 基于剪切效应的气动雾化喷嘴试验研究 排灌机械工程学报，（）：，（）：（），：郑加强，徐幼林 环境友好型农药喷施机械研究进展与展望 农业机械学报，（）：，（）：（）董福龙，周宏平，茹煜，等 对冲喷头喷雾场液滴分布特性试验 农业机械学报，（）：，（）：，（），第 期刘秀娟，等不同雾化方式喷雾场的数值分析与试验验证 ，（）：陆大勇，贾和坤，叶泽 高压条件下椭圆孔喷雾破碎雾化分析 江苏大学学报（自然科学版），（）：，（），（）：（），：王景旭，祁力钧，夏前锦 靶标周围流场对风送喷雾雾滴沉积影响的 模拟及验证 农业工程学报，（）：，（）：（），：，：傅泽田，王俊，祁力钧，等 果园风送式喷雾机气流速度场模拟及试验验证 农业工程学报，（）：，（）：（）曹建明 喷雾学 北京：机械工业出版社，：（责任编辑张文涛）檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨（上接第 页），（）：（），（）：杨治平，陈明昌，张强，等 不同施氮措施对保护地黄瓜养分利用效率及土壤氮素淋失影响 水土保持学报，（）：，（）：（）石小虎，曹红霞，杜太生，等 膜下沟灌水氮耦合对温室番茄根系分布和水分利用效率的影响 西北农林科技大学学报（自然科学版），（）：，（），（）：，（）雷宏军，肖哲元，肖让，等 水、肥、气耦合滴灌对温室番茄生长和品质的影响 干旱地区农业研究，（）：，（）：（）（责任编辑张文涛）