基于阻力特性的布风器结构优化
冯国增，周志平，雷淑雅，刁海兵，沈艳琳，张东辉
(江苏科技大学 能源与动力学院，江苏 镇江 212003)
摘 要: 本文针对目前船舶中常用的布风器形式，通过调整3 种结构参数（即A，B，C），详细了解其对流动阻力系数的具体影响，并数值分析得出结构参数A 影响最大。本研究采用Ansys Fluent 软件，通过对布风器内的流动进行大量的模拟分析，提出更为优化的结构参数设计，数值计算出改进后布风器的总阻力系数比原始布风器减少20%，说明改进后布风器送风的均匀性趋于良好。本文的工作为布风器更为合理的设计打下基础，也为后期的降噪设计研究抛砖引玉。
关键词 ：布风器；阻力系数；送风均匀性；降噪
0 引 言
布风器是空调通风系统的重要组成部分，它在空调送风管路中起着终端散流作用，布风器有很多分类和分类方法，按其结构分类可分为喇叭式、细孔式、转球式和百叶箱式等[1-2]，按其安装位置分类可分为顶置式和壁置式，按其工作性能分类可分为直布式和诱导式，直布式的工作原理是将空调器的空气直接喷入舱室，只在室内引起一次空气流动，原理上这种工作方式会产生很大的噪声，为了保证其静音效果，在其近风口处加入大型消音箱来减小噪声，而这样做的缺点是体积过大，仅适用于大型货舱或客船的居住舱室[3-4]。诱导式布风器的原理是高速喷出一次风，引起一部分室内空气（二次风）进入布风器与一次风进行混合，其优点是能使舱内气体的速度和温度分布均匀[5-6]。从相关企业的市场情况看，目前应用于船舶舱室的先进布风器，还是处于引进仿制阶段，价位也较高，难以适应船舶舱室变风量送风的要求，布风器的设计改进是一个有待解决的课题。
从目前的研究现状看，有关布风器的流动模拟目前已有一些先期的研究。丁亮等[7]对一种顶式布风器的内部结构优化进行了模拟分析，已改善现有布风器静压箱体积小、消音棉厚度薄的设计缺点，并将其应用于海洋工程领域。孙丽颖等[8]系统研究了4 种不同布风器对船舶舱室气流分布影响的综合评价，为提高船舶居住舱室的舒适性提供了理论基础和技术参考。为进一步满足船舶舱室的空调温度灵活调节要求，邵宜南[9]提出应用于船舶的电加热布风器，布风器内加装电加热器的单风管空调系统不占用更多空间，又能保持船舶舱室的送风量和新鲜空气量，满足不同舱室需求，目前已被船东和船员广泛接受。
用来检验农村居民点规模在局部地区是否存在统计上显著的高值和低值，用于可视化表达空间“热点”区和“冷点”区，用于农村居民点规模分异研究[24]，其计算式为：
本文针对目前船舶中常用的布风器形式，通过结构参数的调整，了解其对流动阻力系数的具体影响，为布风器的优化设计打下基础，为布风器的降噪设计打下基础。本研究采用Ansys Fluent 软件，对布风器内的流动进行数值模拟，计算在不同工况下的阻力系数，详细了解不同结构参数对布风器阻力特性的影响。以期获得更为合理的设计方案，也为后续布风器的降噪设计打下基础。
1 物理模型与边界条件
布风器结构如图1 所示，布风器内部空气流动属于不可压缩粘性流体流动。模拟计算过程中采用的湍流模型选择Realizable k-ε 模型，对时间的离散采用2 阶隐格式。计算中采用的流体为空气（Pr=5.42），采用四面体非结构化网格对已建立的3 维几何模型进行网格划分，如图2所示。该网格共划分出3 个计算区域：进风段采用六面体网格；箱体、出风口均采用四面体网格，网格总数22 万。进风口选择为压力入口边界条件，出风口选择流动出口边界条件，布风器内部壁面设为壁面边界条件，流动介质为空气，最大风量为250 m3/h。针对此结构型式的布风器，模拟中对比了不同结构参数对布风器流动特性之影响，对每一结构形式分别计算了不同入口速度下的总阻力系数
为了有效地管理好高邮湖，充分发挥水利部门的统筹管理职能，切实保护好国家湖泊资源，要修订相关法律法规，完善湖泊各项规章制度，如高邮湖水环境影响评价制度、水行政执法巡查制度、湖泊联席会议制度、排污口设置审批制度、湖泊岸线规划制度、取水许可制度等。同时借鉴国外的经验，建立和完善目前的水功能区划制度、农业及养殖业水资源保护制度、污染事故紧急处理及预警制度、水资源补偿费征收制度等。
可知，以原有布风器内部结构参数A 取228 mm时为基准，当A 值增大时，计算出局部阻力系数值偏高，表明气体流动损失大；当A 值逐渐减小时，局部阻力系数大幅降低，并最终趋于平稳状态，整个布风器内部空气流动性能得到大幅改善。因此，影响因子A 对该布风器的气体流动阻力系数影响颇大，且当A 取值为215 mm 时，局部阻力系数能降低18%。
图 1 原始布风器三维结构图
Fig. 1 The original three-dimensional structure of the air distributor
图 2 布风器内部的计算区域
Fig. 2 Calculation area inside the air diffuser
速度场流线图，它们主要是为了能更加直观地表达各处速度大小的分布，流线越密集的地方速度就越大。从图中还可以直观地看出速度流线的走向，原始布风器不合理的内部结构导致出来的速度流线纠缠在一起，在局部形成旋涡。
对于上市公司而言，企业的管理相对复杂，企业的日常经营和管理很大程度上都依赖公司的管理层。管理者的积极性和管理效率对于企业未来的发展至关重要。但是纵观现阶段企业管理层的领导来看，对于领导者还缺乏有效的约束和激励措施。管理层的领导者对于上市公司内部治理结构来说，地位还比较尴尬，其既不是企业股权的所有者，也不是企业的投资者。其薪资水平大多都是来源于企业的绩效考核和评价。这种单一化的考核模式对于管理层领导者来说是不够公平的。纵观现阶段上市公司的发展来看，对于管理层领导者既没有系统的激励措施，同时也缺乏完善的约束机制。使得对于中低层领导者的管理还不够健全。不利于企业治理结构的改革。
为保证模拟结果的准确性，首先对布风器进行网格验证实验，如图3 所示。当网格数小于22 万，总阻力系数的数值较为波动，往后在网格数的不断增加的情况下，布风器的阻力系数值趋于平稳。所以模拟网
格数选择20 万左右。
图 3 原始布风器网格验证示意图
Fig. 3 Schematic of the original air distributor grid
2 模拟结果与讨论
2.1 布风器流动性能的影响因子模拟分析
对布风器3 个方向进行调整分别是布风器内部挡板与壁面之间的距离A、调整X 方向2 块挡板的长度B、调整静压箱的宽度C，通过计算流体动力学方法对这3 个影响因子进行模拟计算，再根据计算结果对布风器内部结构优化，3 个影响因子A，B，C 如图4 所示。
重视污水再利用，把污水处理设施建设作为城市基础设施建设的重要内容。对城市管道实行雨污分流，同时安排好污水回用设施的建设，使经过达标处理的污水能循环使用。建立污水集中处理良性运行机制，积极推进城市污水处理厂企业化经营步伐，同时应积极拓展融资渠道，鼓励和吸引社会资金和外资投向城市污水处理和回用设施项目的建设和运营，进而实现污水资源化。
图 4 布风器影响因子示意图
Fig. 4 Air distributor impact factor diagram
2.1.1 影响因子A 的计算结果
适当调整布风器内部挡板与壁面之间的距离（即A=212，213，214，215，216，222，228，234，240，246 mm），建立三维几何模型，采用四面体非结构化网格进行网格划分，计算结果如图5 所示。
模拟中通过调整布风器内一些关键结构参数，比如挡板与壁面间的距离（即A），x 方向挡板的长度（即B），静压箱y 向的宽度（即C）等，以达到减小布风器流动阻力的目的。模拟中进行了6 种方案的比较：1）挡板与壁面间的距离；2）x 方向消音棉的长度所示；3）增加Y 方向外框的长度；4）同时改进方案1 和方案2；5）同时改进方案2 和方案3；6）同时改进方案1、方案2 和方案3。根据模拟结果，提出优化的布风器结构设计方案，并对其进行了降噪测试。
图 5 局部阻力系数随影响因子A 变化图
Fig. 5 Local resistance coefficient with the impact of factor A changes
Realizable κ-ε 模型的湍动能及其耗散率输运方程为：
2.1.2 影响因子B 的计算结果
适当同时调整X 方向2 块挡板的长度（即B=202，203，204，205，206，207，208，213，218，223，228，233，238，243 mm），建立三维几何模型，采用四面体非结构化网格进行网格划分，计算结果如图6 所示。
适当同时调整X 方向两块挡板的长度（即C=390，400，410，420，430 mm），建立三维几何模型，采用四面体非结构化网格进行网格划分，计算结果如图7 所示。
图 6 局部阻力系数随影响因子B 变化图
Fig. 6 Local resistance coefficient with the
impact of factor B changes
由图8 可知，原布风器进出口风压差较大，在额定风量下，达到105 Pa 左右。由于原始布风器不合理的内部结构，导致流体刚进入静压箱时压力损失严重，而且在出风口处会形成旋涡；改进后的布风器内部流动压力损失下降了很多，流线分布也更均匀。从图9 可以看出，原布风器在入口处的速度逐渐降低，当进入静压箱后截面上出现一大块速度骤增区域，而流体在出来的时候速度又骤降许多，可见这样不合理的内部结构导致内部速度场分布不均匀且速度损失大。从图10 典型流线图也可以清楚了解到，原布风器出口的流线中心区域流速较小，四周区域流速较大，使得送风均匀性不佳。
2.1.3 影响因子C 的计算结果
但是，张恨水先生的文学功底确实为我们敬佩的。从头到尾，文中一直贯穿着的古诗、古情，让我在读的时候不是想尽快知道结果，而是慢慢享受这种韵味儿。这就是所谓的“大家”吧。
引理 1 对于样本模型Yi=m(Xi)+εi,i=1,2,···,n.^mH(x)为m(x)具有带宽矩阵H的局部线性估计,并满足文献[13]中的正则条件.设x为一个非边界点,则在给定X1,X2,···,Xn 下^mH 的偏倚为
图 7 局部阻力系数随影响因子C 变化图
Fig. 7 Local resistance coefficient with the impact
of factor C change chart
适当调整静压箱的宽度，通过数值模拟计算得出增加20 mm 时，整个气体流动性能得到进一步改善，但由于设备安装空间受限，所以在此不做特别详细的研究。
2.2 小结
基于当前实用型船用布风器，本节先分析其内部结构构造情况，再分别调整3 种不同的影响因子（即A，B，C），然后采用计算流体动力学方法进行数值模拟，最后整理所有数据结果，分析得出：影响因子A 对该布风器内气体流动阻力影响至关重要，影响因子B 则仅次之，过度减小结构参数B 会影响布风器降噪效果，影响因子C 则是在受限要求内做稍微的调整，使得布风器内空气流动性能更佳。
3 原始与优化后布风器的流动性能对比分析
首先对原始布风器进行数值模拟，在额定风量下入口速度为8.85 m/s，压力云图、速度云图以及速度流线图如图8～图11 所示，由此可以很好了解布风器内部的压力和流速分布。
可知，不断减小布风器内部结构参数B 的值，局部阻力系数的降低速率由快到慢，直至趋于稳定；当影响因子B 取值为206 mm 时，改进后布风器的局部阻力系数降至最低，是原有布风器的85%。由此说明，优化结构参数B 对改善布风器气体流动性能也起着不可或缺的作用。
图 8 原始布风器的压力场
Fig. 8 Original air blower pressure field
图 9 原始布风器的速度场
Fig. 9 Velocity field of the original air distributor
图 10 原始布风器速度场流线图（正视方向）
Fig. 10 Original air distributor velocity field flow chart （frontal direction）
图 11 原始布风器速度场流线图（俯视方向）
Fig. 11 Original cloth velocity field flow diagram（looking down）
式中：
为由于平均速度梯度引起的湍动能产生；
为用于浮力影响引起的湍动能产生；
为可压速湍η流脉动膨胀对总的耗散率的影响； ；
和
为常数；
别为湍动能及其耗散率的湍流普朗特数。在Fluent 中，作为默认值常数，
进一步对原始布风器的4 种工况（即进口流量100 m3/h，150 m3/h，200 m3/h，250 m3/h）进行分析研究，发现大部分工况压力分布极不均匀，在出口截面处的速度大小差异太大，经过研究分析，得出挡板与壁面间的距离大小对整个布风器的流动性能影响最大。而且，布风器内的噪声与流动阻力存在着相互关联，一般而言，流动阻力越大，噪声也越大。所以这里从减小布风器内的流动阻力角度对其内部结构进行优化。经大量的模拟分析，最终采用以下综合的优化方案：“回”字型布风器的消音棉沿Y 向内缩共13 mm，沿X 方向上下均向左缩减22 mm，另外适当扩大外框Y 向距离变为430 mm。
图12～图15 是布风器改进方案6 的压力云图、速度云图以及速度流线分布，入口流速为8.85 m/s。由图14 可知，改进后的布风器流线分布更加平顺，走向基本一致，没有出现速度骤增骤减的情况，而且相对原布风器来讲，流线分布较为混乱的区域也更少，也就是布风器内的旋涡区域得以缩小。当入口流速为8.85 m/s时，改进方案6 的布风器压力损失大约为75 Pa，较原始布风器下降很多。这表明改进后的效果显著良好。
高校图书馆不仅肩负起校园建设的诸多使命，而且全面推动了校园文化建设的精神与培育和传播，在培养高校人才以及推进学术发展与研究方面具有突出的贡献。
图 12 改进方案6 布风器的压力场
Fig. 12 Improvement scheme six pressure field of the air diffuser
图 13 改进方案6 布风器的速度场
Fig. 13 Improvement program six velocity field of the air diffuser
图 14 改进方案6 布风器速度场流线图（正视方向）
Fig. 14 Improvement scheme six blower velocity field flow chart （frontal direction）
图 15 改进方案6 布风器速度场流线图（俯视方向）
Fig. 15 Improvement scheme six blower velocity field
flow chart （looking down）
图 16 不同改进方案下的阻力系数变化
Fig. 16 Changes in the resistance coefficient under different improvements
图 17 不同改进方案下的静压降
Fig. 17 Static pressure drop with different improvements
从图16 和图17 可以看出，在同一工况下，原始布风器阻力系数最大，这意味着气流在布风器内的能量损失也是最大；通过调整布风器内部结构参数，可使布风器的流动性能得以改善。出口截面的压力和速度分布均匀性比原始布风器更好。
肿瘤标志物是一种能够反映肿瘤存在情况的物质，其是通过患者的肿瘤细胞分泌情况在人体的细胞、血液以及体液中的表现来反映肿瘤的情况。此种标志物不但可以存在与良性肿瘤中，并且在恶性肿瘤中也同样存在，在肿瘤早期，这种肿瘤标志物就会出现，临床上通过对标志物水平的检测，来对肿瘤患者疾病的发展情况进行判断 [3]。
表1～表3 为原布风器和不同改进方案的静压损失和总阻力系数对比，在不同的改进方案中，方案6 的总阻力系数比原始布风器降低20%左右，效果最为显著。显然，改型后的布风器减阻效果很明显，更具有实用价值。
布风器阻力系数和降噪设计之间存在紧密联系，气动噪声的主要来源就是流体的旋涡，一般而言，布风器内旋涡强度越大，噪声就越大。下一步研究的主要方向是将阻力特性、流动旋涡特性与噪声相关联。
表 1 原始布风器的压损和总阻力系数
Tab. 1 The original air pressure loss and the overall drag coefficient
表 2 改进方案2 的静压损失和总阻力系数
Tab. 2 Improved scheme 2 static pressure loss and total drag coefficient
表 3 改进方案6 的静压损失和总阻力系数
Tab. 3 Static pressure loss and total resistance coefficient for improved scheme six
4 结 语
通过系统的模拟分析，研究发现：所提出来的6 种方案均能有效改进此“回”字型布风器。改进方案1～改进方案3 均是从单方面因素进行改进的，其中改进方案2 的效果更为明显，即消音棉的长度适当减小，这意味着这一因素是影响布风器综合性能的重要影响因子。综合比较，改进方案6 的方案局部阻力系数可降至1.60，达到预期目标，布风器的结构优化方案优先选择此改进方案。这些工作也为后期的降噪设计打下基础。
参考文献：
[1]JG/T20-1999 空气分布器性能试验方法[S], 1999.
[2]刘鑫, 韩星星, 郭宝坤, 等. 船用布风器出口流道的结构优化[J]. 制冷与空调, 2017, 17(10): 24-28+79.
[3]郭宝坤, 李慧子, 刘鑫, . 船用布风器冬季工况射流流场模拟及试验研究[J]. 制冷与空调, 2016, 16(12): 26-30.
[4]王在忠. 某船空调噪声分析与降噪改进研究[J]. 中国修船,2014, 27(4): 24-28.
[5]权崇仁, 王洋, 于立庆, 等. 舰船典型舱室气流组织数值模拟[J].中国舰船研究, 2015, 10(6): 107-113.
[6]秦晓. 舰船住舱气流组织数值研究[D]. 大连: 大连海事大学,2015.
[7]丁亮, 魏成革, 高指林, 等. 一种顶式布风器的优化设计[J].中国造船, 2012, 53(S2): 180-183.
[8]李以通. 布风器形式对船舶居住舱室气流分布的影响研究[D].哈尓滨: 哈尓滨工业大学, 2013.
[9]邵宜南, 边建松. 电加热布风器[J]. 机电设备, 2003(3): 23-31.
Research on the fluid characteristics of wind distributor
FENG Guo-zeng, ZHOU Zhi-ping, LEI Shu-ya, DIAO Hai-bing, SHEN Yan-lin, ZHANG Dong-hui
(Jiangsu University of Science and Technology, School of Power and Energy, Zhenjiang 212003, China)
Abstract : In this paper, a detailed understanding of the specific influence of the three kinds of structural parameters (A,B, C) on the flow resistance coefficient is given in view of the current form of the air distributor currently used in ships. Numerical analysis shows that the structural parameter A has the greatest impact. In this study, Ansys Fluent software was used to conduct a lot of simulation analysis of the flow in the air blower, and a more optimized design of the structure parameters was proposed. The total drag coefficient of the air blower was calculated to be 20% less than that of the original air blower.Indicating that the improved air distribution of the air blower tends to be good. The work of this paper laid a foundation for a more reasonable design of the air diffuser, and also made a valuable contribution to the later research of noise reduction design.
Key words: wind distributor；pressure drop factor；air supply uniformity；noise reduction
中图分类号 ：F24
文献标识码： A
文章编号： 1672 - 7649(2019)07 - 0071 - 06 doi：10.3404/j.issn.1672 - 7649.2019.07.014
收稿日期: 2017 - 12 - 29
基金项目: 江苏省产学研合作计划资助项目(BY2016073-10)
作者简介: 冯国增(1971 - )，男，副教授，研究方向为传热传质学。
标签：布风器论文; 阻力系数论文; 送风均匀性论文; 降噪论文; 江苏科技大学能源与动力学院论文;