变容积密集烤房的CFD分析与试验研究

叶大鹏 沈碧河 张炳辉 黄俊炜 谢立敏

摘要：為保障密集烤房装烟密度，降低烘烤的能源消耗，研发一套适用于密集烤房的变容积系统。在完成变容积装置的设计后，基于CFD方法模拟分析装置与烟叶的不同距离对烤房内部气体分布均匀性的影响。通过烘烤试验获取实际烘烤数据，对模拟值加以验证。试验结果表明：当隔板与烟叶距离为0cm时，流速不均匀系数Kv为0.40，温度不均匀系数Kt为0.41，距离为10cm时，Kv=0.41，Kt=0.43；距离为20cm时，Kv=0.42，Kt=0.49。装烟区9个测量点的温度模拟值与实测值基本吻合，误差在6%以内。变容积烤房在装烟量为一半时，相比常规烤房的燃料消耗可节约13.4%。研究结果表明：当隔板与烟叶距离为0cm时烤房内部的气体分布最均匀；CFD模型与数值模拟结果具有可靠性；变容积装置具有较好的保温效果，可保证装烟密度，降低烤烟能耗。

关键词：密集烤房；变容积；CFD方法；数值模拟

中图分类号：S229+.2  文献标识码：A  文章编号：2095-5553 （2024） 03-0090-07

CFD analysis and experimental research on variable volume bulk curing barn

Ye Dapeng1， 2， Shen Bihe1， Zhang Binghui3， Huang Junwei1， Xie Limin1， 2

（1. College of Mechanical and Electrical Engineering， Fujian Agricultural and Forestry University， Fuzhou，350100， China; 2. Key Laboratory of Agricultural Information Perception Technology in Fujian Province， Fuzhou，350100， China; 3. Tobacco Science Research Institute of Fujian Tobacco Monopoly Bureau， Fuzhou， 350003， China）

Abstract：

In order to ensure the density of smoke in the bulk curing barn and reduce the energy consumption of baking， a variable volume system suitable for the bulk curing barn was developed. After the design of the variable volume device was completed， the influence of different distances between the device and tobacco leaves on the uniformity of gas distribution inside the bulk cuing barn was simulated and analyzed based on the CFD method. The actual roasting data was obtained through the roasting test， and the simulated values were verified. The test results showed that when the distance between the separator and the tobacco leaf was 0 cm， the non-uniform coefficient of flow velocity Kv was 0.40， the non-uniform coefficient of temperature Kt was 0.41， Kv=0.41， Kt=0.43 when the distanc was 10cm， Kv=0.42， Kt=0.49 when the distanc was 20cm. The simulated temperature values of 9 measurement points in the cigarette loading area were basically consistent with the measured values， with an error of within 6%. When the capacity of the variable volume curing barn was half， the fuel consumption can be saved by 13.4% compared to the conventional curing room. The research results indicated that when the distance between the partition and the tobacco leaves was 0cm， the gas distribution inside the curing room was the most uniform. The CFD model and numerical simulation results had a certain degree of reliability. The variable volume device has good insulation effects， which can ensure the density of tobacco loading and reduce energy consumption of flue-cured tobacco.

Keywords：bulk curing barn; variable volume; CFD method; numerical simulation

0 引言

密集烤房是专门用于卷烟烟叶烘烤的设备，具有容量大、装烟密度高等特点［1］。烤烟时密集烤房的装烟密度对烤烟的经济效益有直接关系，装烟密度过高或过低对于烤烟的经济效益都有不良的影响。烤房的能源总消耗与用电量随着装烟密度的提高而增加，但分摊到单位重量干烟上的烘烤成本则相对降低。因此为了节省能耗，提高经济效益，密集烤房的装烟密度不宜过低；但过高的装烟密度又将对烤后的烟叶等级造成影响，导致均价降低，降低经济收益［2］。另外，不合适的装烟密度会对烟叶内在化学成分协调性和评吸质量带来不利影响［3］。相关研究表明，装烟密度为70kg/m3时下烤后烟叶整体品质最高，经济效益最好［4， 5］。

目前的卷烟生产过程中，鲜烟叶采用逐片采摘烘烤的方式。由于天气或水肥差异等各方面因素的影响，烟叶成熟速度不同，在一个采收周期内收获的烟叶装到烤房后会有所剩余或不足，大多数烟农装烟时只能改变装烟密度将烤房装满，容易导致装烟密度过高或过低，无法达到最优装烟密度，从而增加能耗，影响烟叶品质，降低烤烟经济效益。因此，寻求一种密集烤房容积调节的方法，可根据装载的烟叶数量调节烤房的实际烘烤区域，使烟叶保持在最佳装烟密度下烘烤显得极为必需。

近些年，CFD方法在研究各种农业设施室内的气体流场分布中被广泛应用［68］。在各种畜禽舍的通风换气与湿热传递分布等问题的研究中，多运用CFD方法进行分析，如：运用计算流体动力学方法对猪舍冬季气流环境进行研究［9］、对保育猪舍垂直通风进行CFD模拟及挡风板优化设计［10］。一些屋舍结构的优化问题也可通过CFD方法得到较好的解决［11， 12］，对密集烤房的诸如通风孔尺寸、位置之类的结构优化也多使用CFD模拟仿真方法［13， 14］，但有关烤房的容积调节研究却鲜有报道。

基于以上讨论，本文从划分烤房区域，隔离装烟区，减小加热空间的角度研究设计了一套适用于密集烤房的容积调节系统。使每次烤烟总能以最佳装烟密度进行烘烤，通过烤烟对比实验确定隔板厚度，利用CFD三维数值模拟方法分析隔板与烟叶的不同距離对装烟区温湿度分布均匀性的影响，探究隔板最佳安装位置，提高烤烟均匀性，并降低能源消耗。

1 材料与方法

1.1 新型密集烤房基本结构与工作原理

目前常用的密集烤房根据热空气流动方向的不同可分为气流上升式与气流下降式，二者除风机的位置与朝向不同外结构基本相同［15］。实验所用烤房是一座位于福建省三明市永安市的气流上升式密集烤房，烤房可分为燃烧室和烘烤室，烘烤室内部尺寸为8.0m×2.7m×3.5m（长×宽×高）。烘烤室的天花板、左右墙面以及大门门板均由5cm厚的聚氨酯夹芯板构成，内部设有烟架，用于承载烟杆。烤房工作时燃烧室内的气体被加热炉加热，由风机吹向出风口，出风口出来的热空气往大门方向移动的同时慢慢上升，与温度较低的烟叶进行热交换后到达天花板，通过回风口排出。烤房结构与热风流动路线如图1所示。

在图1所示的烤房中搭建变容积装置，该变容积装置主要由保温面板与铝合金骨架构成，六块面板拼接后形成隔板，搭接在烟架上，将烤房内部分隔为装烟区与无烟区，如图2所示。

为保证装置的密闭性，隔板的保温材料尺寸略大于烤房壁面之间的距离，面板互相挤压减少空隙，因此，材料需要具备一定的可压缩性和压缩复原性；此外，所选取的保温材料需有较好的隔热能力以减少装烟区的能量损耗。常见的保温材料及其物理性质如表1所示。

气凝胶毡保温性能最好但是复原性差且密度过大；聚氨酯保温性能次之，但易出现鼓包与开裂；酚醛脆性较大，而聚苯乙烯则吸水性高，均无法满足隔板的要求。珍珠棉保温性能较好，压缩性与压缩复原性良好，密度较小但容易弯曲，可以通过骨架加固，使隔板烘烤时不会在风压的作用下发生变形。综上，最终选定EPE高密度珍珠棉作为变容积装置的保温材料。

变容积装置的隔板厚度对装置的保温性能影响巨大，隔板厚度越大，不同面板之间相互挤压后的接触面积越大，密闭性保温性越好，但装置的造价也越高。通过测定不同厚度下装烟区与无烟区的温度差确定面板厚度。

综合考虑每千克干烟的平均烘烤成本与实际生产情况，研究的装烟量为最佳装烟密度条件下装满烤房容积一半，即在靠近燃烧室的一侧，每层装4m长的烟叶。

1.2 烤房建模与网格划分

根据现有的烤房结构，等比例建立烘烤室的几何模型。忽略烟架，观察窗等所占比重较小的结构。模型由变容积装置分隔成两部分，右侧装烟区内的三层烟叶简化为三个多孔介质区域，模型结构简单，在ICEM CFD 19.2中划分六面体结构化网格并进行局部加密，共产生347734个节点，344667个单元，如图3所示。

1.3 控制方程与模型选择

烟叶烘烤过程中烤房内部气体保持流动，整个过程遵循质量、动量和能量三大守恒定律［18］；

质量守恒方程（连续性方程）如式（1）所示。

аρаt+а（ρu）аx+а（ρv）аy+а（ρw）аz=Sm（1）

式中：u、v、w——速度矢量在x、y、z方向分量，m/s；ρ——烤房内流体密度，kg/m3；t——时间，h。动量守恒方程（N-S方程）如式（2）所示。

ааt（ρui）+ааxi（ρuiuj）=－аpаxi+аτijаxj+ρgi+Fi（2）

式中：p——静压，Pa；τij——应力张量；gi、Fi——i方向上的重力体积力和外部体积力，N。

能量守恒方程如式（3）所示。

а（ρT）аt+аSymbolQC@（ρuT）=аSymbolQC@kcpgradT+ST（3）

式中：cp——比热容，J/ （kg·K）；T——温度，K；k——流体的传热系数；ST——黏性耗散项。

烟叶烘烤的空氣流动属于湍流，选择标准的k-ε模型进行求解［16］，求解过程中将挂烟区域视为孔隙率为0.87%［17］的连续多孔介质。

1.4 边界条件与数值计算

运用Ansys Fluent软件的2022R1版本进行数值求解。将烘烤室的热空气入口设置为速度入口，平均风速3.5m/s，温度42℃，湍流模型采用标准的k-ε模型，入口处设置湍流参数，湍流强度I=3.7%，水力直径D=0.7m。出口设置为常压状态的压力出口。烤房壁面与天花板设置为绝热壁面，不考虑其与外界环境的换热，即以压力基进行瞬态求解，压力速度耦合方式为SIMPLE，压力和动量采用二阶迎风离散格式。

2 试验与分析

为验证变容积装置对烤烟燃料消耗的影响，选取同一烤房群内的两座相邻烤房进行对比试验。一座为安装了变容积装置的试验烤房，一座为未装变容积装置的对比烤房。在两座烤房内部搭建温湿度监测系统，检测并记录实时烘烤数据。在烤房两侧壁面与中间烟架三个平面上共布置18个温湿度探头，探头型号为RS485型pvc探头。在两座烤房以相同装烟量，相同开烤时间进行烘烤，通过RS-XZJ-100-Y型监控主机每隔10min记录一次烘烤数据并发送至云平台。同时记录两座烤房的燃料消耗情况。探头的安装位置如图4所示。

烘烤试验开始的时间为2023年7月4日，于福建省三明市永安市小陶镇进行，烘烤时间为250h。

2.1 隔板厚度的确定

变容积装置隔板的密闭性和保温性对于烤房的能源消耗有很大影响。隔板的保温性越好，装烟区的热量往无烟区的传递越弱，无烟区的温度也就越低。无烟区与外界温差越低，烤房的热量损耗越小，节能效果越明显。在装烟量为烤房容积一半的情况下测试隔板厚度为5cm，10cm时的变容积装置保温性以确定最佳隔板厚度。

选取以隔板为对称平面的四个对称观测点P6、P7、P14和P16所采集的数据进行对比。其中，P6、P14位于无烟区内，分别距离隔板约90cm与200cm；P7、P16位于有烟区内，距离隔板约90cm与200cm。烘烤过程中，厚度为5cm的隔板隔离下装烟区与无烟区的温度变化与温差如图5所示。

厚度为10cm的隔板隔离下装烟区与无烟区的温度变化与温差如图6所示。

由图5和图6可以看出，随着烘烤时间的推移，位于有烟区的P7与P16位置的温度逐渐上升，在烤房的温度调控系统作用下，温度曲线与烤烟烘烤工艺曲线吻合。位于无烟区的P6与P14位置的温度曲线出现周期性波动，温度下降处的记录时间为夜晚，应是夜晚温度降低所致。P6与P7、P14与P16的温度差值也出现周期变化，总体呈现上升趋势，表明随着烘烤温度的上升，隔板两侧的温度差呈增大趋势。P6的温度始终小于P7，P14始终小于P16，图5与图6中的监测点温度差值十分接近，表明10cm的隔板保温性能与5cm相近。但10cm隔板的厚度增加一倍，设备成本提升，综合考虑设备造价与节能效果，变容积装置的隔板厚度定为5cm。

2.2 隔板最佳位置确定

2.2.1 不同隔板位置的装烟室气体流速场分布

在模拟仿真的后处理阶段，选取烤房内距离地面0.35m、1.75m和2.85m的三个水平截面Z1、Z2、Z3与一个垂直地面与隔板的中心截面Y，展示烤房的气体流速场分布情况。分别创建隔板与烟叶距离为0cm、10cm与20cm的不同模型，三个模型在Z1、Z2、Z3截面上的气体流速分布如图7所示。

图7（a）所示的隔板与烟叶距离为0cm时，空气从入口进入后向上扩散，Z1截面上方均为烟叶区域，截面上各处气体扩散速率相近，气流撞到隔板后回流，隔板附近风速出现周期性波动。Z2截面在隔板一侧风速高于燃烧室一侧，高速区较分散。Z3气流在往出口运动时速度逐渐增大。

图7（b）所示的隔板与烟叶距离为10cm时，在Z2截面的气流高速区集中的隔板一侧，相比图7（a）中的Z2截面更集中，应是隔板与烟叶区域存在间隙，进入烘烤室的气流到达隔板后通过间隙上升至天花板，经过Z2截面时往燃烧室运动时速度下降明显。由于气体在隔板处快速上升，使得Z3截面的气流往出口运动的气体流速增加，平均流速大于图7（a）。

图7（c）所示的隔板与烟叶距离为20cm时，隔板与烟叶区域间隙增大，间隙处的气体流速增大，Z2截面高速区更集中，Z3截面平均流速大于图7（b）。

2.2.2 不同隔板位置的装烟室温度场分布

三个模型的不同温度场分布如图8所示。

由图8（a）可以看出，隔板与烟叶距离为0时，热空气由进风口进入后，在到达隔板的过程中，从烟叶区域的最低端向上扩散。在同一水平面各处扩散速度相近，使得同一水平高度的烟叶区域受热情况相近。图8（b）中Z1、Z2截面温度几乎一致。Z3截面由于出口变窄，在出口附近的角落气流难以到达，产生低温区，但该区域无烟叶，影响较小。随着时间的推移烟叶慢慢达到与进口空气相近的温度。

图8（c）所示隔板与烟叶距离为10cm时，热空气在到达隔板时，由于隔板与烟叶存在空隙，入口的空气进入烘烤室后很快便从空隙处排向出口，温度传递呈“C”型向内发展，烟叶区域慢慢达到均匀的温度，但比入口处温度略低。图8（d）中Z1截面温度均匀，Z2截面在隔板附近的温度高于燃烧室一侧，Z3截面在装烟区的中心区域温度较低。

图8（e）所示隔板与烟叶距离为20cm时，隔板与烟叶空隙增大，热空气进入烘烤室后更快从出口排出，“C”型温度传递更明显，由于热空气与烟叶接触时间变短，出口下方的烟叶升温更慢，出现低温区。图8（f）中Z1与Z3截面温度均匀，Z2截面从隔板往燃烧室温度递减，两端温度差异较大。

2.2.3 最佳隔板位置

引入不均匀系数K［19］描述烤房内部各个物理场的不均匀程度，评价隔板位置对烤房环境的影响。速度场不均匀系数Kv、温度场不均匀系数Kt的公式如式（4）所示。

Kv，t=1m∑mj=11n∑ni=1（Ti－T-）2（4）

式中：j——水平截面数；i——该水平面上的测量点数；T-——该平面各个测量点的温度平均值;m——截面总数;n——同一截面的测量点总数。K值越小表明均匀程度越好。

在三个Z截面的装烟区分别取装烟区的对角线交点和两条对角线上的四个四等分点计算不均匀系数，不同隔板位置的Kv与Kt见图9所示。

由图9可见，随着隔板位置与烟叶的距离增大，气体流速的不均匀系数Kv与温度不均匀系数Kt均增大。隔板位置的变化与气体流速的不均匀系数曲线较为平缓，位置变化与温度不均匀系数曲线陡峭。综上，隔板距离为0cm时烘烤室的气体流速场与温度场均匀程度最好。

2.3 模型可靠性验证

记录试验烤房中装烟区的9个监测点的烘烤温度數据，以烘烤时间t=68h时的温度实测数据与仿真模拟值进行对比，结果如表2所示。

2.4 燃料消耗分析

在试验烤房与对比烤房内装烟194杆，占据烤房装烟容积的一半，装烟密度为70kg/m3。以颗粒状生物质燃料作为能量来源。记录两座烤房每日的生物质燃料消耗，如图10所示。

烘烤完毕后，两座烤房的燃料料斗内均无燃料剩余。由图10可以看出，在烘烤过程中，试验烤房每日的燃料消耗量均小于对比烤房，并且在烘烤后期，试验烤房的燃料消耗明显减少。整个烘烤实验中，对比烤房共计消耗生物质燃料820kg，加装了变容积装置的试验烤房共计消耗生物质燃料710kg。相较于对比烤房，试验烤房共节约燃料110kg，约节省13.4%。

变容积装置将烘烤室分隔为装烟区与无烟区，缩小了烤房需要加热的空间。在烤房没有满烤时，从燃烧室进入烘烤室的热空气在遇到保温隔板后更快上升，烟叶更快升温，减少热空气在无烟区运动过程中的能量损耗。另外，变容积装置使无烟区的温度与外界温差减少，降低了烤房无烟区与外界的热交换产生的热量损失，达到节约能源的作用。

3 结论

1） 本文设计了一套变容积装置，并建立了变容积烤房的CFD仿真模型，分析该装置在新型密集烤房内的安装位置对烤烟室气流分布均匀性的影响。模拟结果显示，当隔板与烟叶距离为0cm时，流速不均匀系数Kv为0.40，温度不均匀系数Kt为0.41，10cm时Kv=0.41，Kt=0.43，20cm时Kv=0.42，Kt=0.49，隔板安装在紧靠烟叶的位置装烟区流场分布最佳。

2） 对比烤烟室内流场分布的模拟和实际烘烤数据测，装烟区9个测量点的温度实测值与模拟值基本吻合，相对误差＜6%，表明所建立的CFD模型及模拟具有可靠性。

3） 变容积装置的隔板厚度为5cm，材料为EPE保温棉可保障装置的密闭性与保温性。烘烤试验结果表明，烤房装烟量为50%时，加装了变容积装置的烤房比常规烤房节约13.4%的燃料消耗。

参 考 文 献

［1］刘光辉， 聂荣邦. 我国烤房及烘烤技术研究进展［J］. 作物研究， 2011， 25（1）： 76-80.

［2］黄建， 张海伟， 卢瑞杰. 密集烤房装烟密度对烟叶质量及烘烤效益的影响［J］. 湖北农业科学， 2014， 53（7）： 1595-1598.Huang Jian， Zhang Haiwei， Lu Ruijie. Effects of loading densities on leaf qualities and benefits of bulk cured tobacco ［J］. Hubei Agricultural Science， 2014， 53（7）： 1595-1598.

［3］张滔. 密集烤房装烟密度对烟叶质量的影响［D］. 长沙： 湖南农业大学， 2013.Zhang Tao. The effect of tobacco density in dense curing on tobacco leaf quality ［D］. Changsha： Hunan Agricultural University， 2013.

［4］贺非， 皮本阳， 孙敬钊， 等. 装烟密度对烤后烟叶品质的影响［J］. 农业工程技术， 2020， 40（26）： 32-33.

［5］詹军， 李伟， 武圣江， 等. 密集烤房装烟密度对烤烟上部烟叶质量的影响［J］. 作物研究， 2011（5）： 473-476.Zhang Jun， Li Wei， Wu Shengjiang， et al. Effects of loading density of bulk curing barn on quality of upper leaves in flue-cured tobacco ［J］. Crop Research， 2011（5）： 473-476.

［6］吴照学， 梁伟， 鲍恩财， 等. 基于CFD的发酵床网养鸭舍冬季环境模拟与优化［J］. 中国农机化学报， 2023， 44（12）： 33-39.Wu Zhaoxue， liang Wei， Bao Encai， et al， Simulation and optimization of winter environment for duck house with fermentation bed net based on CFD ［J］. Journal of Chinese Agricultural Mechanization， 2023， 44（12）： 33-39.

［7］汪开英， 李开泰， 李王林娟， 等. 保育舍冬季湿热环境与颗粒物CFD模拟研究［J］. 农业机械学报， 2017， 48（9）： 270-278.Wang Kaiying， Li Kaitai， Li Wanglinjuan， et al. CFD simulation of indoor hygrothermal environment and particle matter of weaned pig building ［J］. Transactions of the Chinese Society for Agricultural Machinery. 2017， 48（9）： 270-278.

［8］吴志东， 房俊龙， 巴文革， 等. 冬季采暖保育猪舍送排风管道组合换气系统设计与评价［J］. 农业工程学报， 2021， 37（10）： 152-158.Wu Zhidong， Fang Junlong， Ba Wenge， et al. Design and evaluation of combined ventilation system with air supply and exhaust ducts for a nursery heating piggery in winter ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2021， 37（10）： 152-158.

［9］周麗娜， 徐述胤， 马精阳， 等. 基于计算流体力学的猪舍冬季气流环境模拟与优化研究［J］. 黑龙江畜牧兽医， 2022（7）： 49-54， 59， 127-128.Zhou Lina， Xu Shuyin， Ma Jingyan， et al. Simulation and optimization of winter airflow environment in piggery based on CFD ［J］. Heilongjiang Animal Science and Veterinary Medicine， 2022（7）： 49-54， 59， 127-128.

［10］房俊龙， 吴爽， 吴志东， 等. 保育猪舍垂直通风CFD模拟及挡风板优化设计［J］. 东北农业大学学报， 2022， 53（5）： 59-68.Fang Junlong， Wu Shuang， Wu Zhidong， et al. CFD simulation of vertical ventilation and optimal design of wind shield in nursery pig house ［J］. Journal of Northeast Agricultural University， 2022， 53（5）： 59-68.

［11］李国建， 杨林， 许金钗， 等. 基于CFD的笼养鸡舍冬季通风方案模拟与验证研究［J］. 黑龙江畜牧兽医， 2021（21）： 57-62， 151-153.Li Guojian， Yang Lin， Xu Jinchai， et al .Simulation and verification studies of ventilation schemes in stacked-cage hen house based on CFD in winter ［J］. Heilongjiang Animal Science and Veterinary Medicine， 2021（21）： 57-62， 151-153.

［12］尚峰军， 赵洁. 采用CFD技术分析和评估蛋鸡舍温热环境质量［J］. 中国家禽， 2021， 43（2）： 61-68.Shang Fengjun， Zhao Jie. Using CFD technology to analyze and evaluate the quality of the warm environment of layer house ［J］. Chinese Poultry， 2021， 43（2）： 61-68.

［13］Martínez-Martínez V， Baladrón C， Gomez-Gil J， et al. Temperature and relative humidity estimation and prediction in the tobacco drying process using artificial neural networks ［J］. Sensors， 2012， 12， 14004-14021.

［14］Bai Zhipeng， Guo Duoduo， Li Shoucang， et al. Analysis of temperature and humidity field in a new bulk tobacco curing barn based on CFD ［J］. Sensors， 2017， 17（2）： 279-279.

［15］詹军， 樊军辉， 宋朝鹏， 等. 密集烤房研究进展与展望［J］. 南方农业学报， 2011， 42（11）： 1406-1411.Zhang Jun， Fan Junhui， Song Chaopeng， et al. Comparative analysis of bulk curing barn and traditional flue-curing barn ［J］. Journal of Southern Agriculture， 2011， 42（11）： 1406-1411.

［16］刘兆宇， 罗会龙， 陈颐， 等. 烟叶烘烤干筋期密集烤房热湿环境数值分析［J］. 江苏农业科学， 2021， 49（1）：160-165.

［17］孙力.基于水分迁移模型变孔隙率密集烤房烟叶干燥数值模拟及实验分析［D］. 昆明： 昆明理工大学， 2020.Sun Li. Numerical simulation and experimental analysis of tobacco leaf drying in the closed bulk curing barn， based on variable porosity and moisture migration model ［D］. Kunming： Kunming University of Science and Technology， 2020.

［18］王福军. 计算流体动力学分析—CFD软件原理与应用［M］. 北京： 清华大学出版社， 2004.

［19］邓书辉， 施正香， 李保明， 等. 低屋面横向通风牛舍空气流场CFD模拟［J］. 农业工程学报， 2014， 30（6）： 139-146.Deng Shuhui， Shi Zheng Xiang， Li Baoming， et al. CFD simulation of airflow distribution in low profile cross ventilated dairy cattle barn ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2014， 30（6）： 139-146.

基金项目：福建省农业信息感知技术重点实验室项目（2021ZDSYS0101）；中国烟草总公司福建省公司科技计划项目（2022350000240084）

第一作者：叶大鹏，男，1971年生，福建宁德人，博士，教授；研究方向为农业智能机械。E-mail： ydp@fafu.edu.cn

通讯作者：谢立敏，女，1985年生，福建南平人，博士，讲师；研究方向为机器人运动控制。E-mail： lucy_min@163.com