汽车空气滤清器进气阻力分析

2016-05-30 03:27袁志群黄飞健刘金武许西安
关键词:数值计算

袁志群, 黄飞健, 刘金武, 许西安

(1. 厦门理工学院 机械与汽车工程学院, 福建 厦门 361024;2. 中南大学 交通运输工程学院, 湖南 长沙 410075)



汽车空气滤清器进气阻力分析

袁志群1,2, 黄飞健1, 刘金武1, 许西安1

(1. 厦门理工学院 机械与汽车工程学院, 福建 厦门 361024;2. 中南大学 交通运输工程学院, 湖南 长沙 410075)

摘要:为了降低汽车空气滤清器进气阻力,应用计算流体动力学方法建立壳体和滤芯的耦合数值计算模型.利用正交分析方法研究进出气口结构参数对进气阻力的影响规律,总结进气阻力随流量变化的关系式,通过滤芯的流量-压降实验验证数值计算方法的可靠性.结果表明:进气阻力随着流量的增加呈抛物线趋势增加;在额定流量工况下,进气阻力随着出气口面积的增加而减小,负载时最大降幅达到37.2%;进出气口圆角半径越大,进气阻力越小,负载时最大降幅分别达到16.6%和36.2%.

关键词:空气滤清器; 进气阻力; 结构参数; 计算流体动力学; 数值计算

进气系统性能是发动机性能的重要保障,结构设计缺陷将导致进气系统阻力增大,影响发动机燃油经济性和动力性.此外,进气阻力增加会导致内部真空度增加,对各接口密封性提出更高要求.空气滤清器是发动机进气系统的重要组成部分.以往进行空气滤清器结构设计时主要依靠经验公式及工程师的经验,通过试制的方法进行.因此,结构设计初期不能准确地预测产品的性能,而需要等到试制品完成后,通过实验的方法分析产品性能,耗时耗力.通过计算流体动力学的方法,能适时准确判断产品性能,并且获得许多实验方法难以实现的数据,找出结构设计中存在的问题,及时进行结构改进,从而降低成本,缩短产品研发周期[1-6].文献[7]研究进出气管管径大小和插入长度对进气性能和消声性能的影响.文献[8-9]对空气滤清器流动阻力特性进行实验和仿真研究.从现有研究成果来看,更多的工作主要是对不同车型的空气滤清器进行改进,不具有通用性,无法全面反映空气滤清器的结构参数对进气阻力的影响规律.基于此,本文对乘用车常用的褶状干式空气滤清器进气阻力进行分析与研究.

1数值计算模型的建立

为了总结空气滤清器结构参数对进气阻力影响的普遍性规律,建立乘用车常用的褶状干式空气滤清器几何模型.该模型由上、下壳体及滤芯等3部分组成.利用UG建立内流场1∶1几何模型, 保留结构全部特征.空气滤清器数值分析模型与空气实际流动区域具有一致性,共分为上、下壳体自由流动区域和中间滤芯多孔介质流动区域.

(a) 计算域      (b) 面网格 图1 空气滤清器网格划分Fig.1 Air filter mesh

在前处理软件ICEM-CFD中,采用OCTREE方法对空气流动区域进行离散,如图1所示.在空气滤清器壳体表面拉伸出与其平行的三棱柱网格,以满足壁面函数的需求,精确模拟壳体表面的附面层.在流动变化剧烈的区域进行网格加密,通过采用不同的网格数目验证了网格的无关性.数值分析模型网格达90万.测试表明,划分的网格单元数足以消除网格密度的影响.

为了简化计算,对数值分析模型作如下3点假设:1) 空气滤清器内部流动速度远小于音速,空气密度变化基本为零,假设为不可压缩流体;2) 在入口流量一定的情况下,空气滤清器内部流动作定常处理;3) 滤芯内部各物性参数各向同性.

空气滤清器内部流场基本控制方程为三维不可压缩雷诺时均方程,湍流模型采用标准k-ε方程,壳体近壁区域采用壁面函数进行修正.计算采用二阶迎风格式,方程求解采用SIMPLE算法.入口采用速度入口边界条件(根据厂家提供的额定空气流量和入口面积转化成速度数据),出口采用压力边界条件,壳体壁面采用无滑移壁面边界条件,滤芯采用多孔介质模型.流体流经多孔介质时遵循Darcy定律,多孔介质流动控制方程是在雷诺时均方程中增加一个负的动量源项,由黏性损失项和惯性损失项两部分组成.对动量源项在滤芯厚度方向进行积分,得到多孔介质模型计算经验公式为

表1 滤芯压降计算值和实验值对比

式中:p为压力;Δx为滤芯厚度;v为渗滤速度;α-1为粘性阻力系数;c2为惯性阻力系数.

通过试验台架测得空气滤清器滤芯的压降(Δp)-流量(Q)变化关系,结合多孔介质模型计算公式,得到计算所需参数.通过数值分析的方法计算滤芯的压降,与实验数据进行对比,误差控制在6%以内,结果如表1所示.表1中:η为误差.表1验证了文中数值计算方法的准确性.

滤芯试验在空气滤清器性能试验台架上完成,其测试原理如图2所示.图2中:左端连接灰尘喷射器,用于测验滤芯的粉尘过滤性能;右端连接流量控制阀和抽气机,用于测验滤芯在不同流量下的性能.

图2 滤芯压降测试原理Fig.2 Testing principle of pressure dropin the filter element

2分析与讨论

2.1空气滤清器进气阻力分析

汽车空气滤清器流动性能的重要指标是进气阻力,其大小直接影响进气效率.进气阻力[10]定义为:当空气流过总成时,其上、下游压力损失.首先,应确定上、下游压力降,然后,根据上、下游由于管道截面积变化引起的流速变化所产生的动压加以修正.进气阻力的定义公式为

图3 空气滤清器进气阻力与流量关系Fig.3 Relationship between air intakeresistance and flow rate in the air filter

空气滤清器的进气阻力试验需要以空气滤清器的额定空气体积流量为标准,按额定空气体积流量的40%,60%,80%,100%,120%分别进行测取.空气滤清器在空载和负载情况下进气阻力与流量(Q)的变化关系,如图3所示.由图3可知:随着进气量的增加,进气阻力呈抛物线关系增长;在负载工况下,增加趋势更加明显,壳体本身所引起的进气阻力占主要部分.因此,合理的结构设计对于降低进气阻力至关重要.

应用最小二乘法对数据进行拟合,可得进气阻力与流量的关系式:1) 负载时,Δpr=0.032 6Q2+1.355 1Q;2) 空载时,Δpr=0.027 2Q2+0.041 7Q.

入口截面压力的分布云图,如图4所示.由图4可知:空气滤清器内部压力呈现递减趋势.空气从进气口进入上壳体以后,在正对入口截面的上壳体壁面上,气流受到阻滞,速度下降,压力上升.因此,出现局部的压力升高.结合图3数据分析可知:气流在空气滤清器内部流动过程中,压力损失主要来自于壳体本身,滤芯所占比重较小,并且随着流量的增加,趋势越明显.

空气滤清器的三维流场,以及入口截面速度云图与流线图,分别如图5,6所示.由图5,6可知:空气滤清器内部流场主要存在2个涡系,分布在上壳体中,分别在平行和垂直于滤芯平面的截面上;涡系的存在使气流不能顺畅地通过滤芯,并且滤芯表面的速度分布不均匀;进气管与壳体过渡部分为直角,气流从进气管进入后,存在局部的分离涡系;下壳体气流平均速度偏小,不能及时从出气口流出.这些因素都降低了进气效率,造成进气阻力偏大.

图4 入口截面压力分布云图 图5 空气滤清器三维流场Fig.4 Pressure distribution of inlet section Fig.5 3D airflow of air filter

(a) 速度云图 (b) 流线图图6 入口截面速度云图与流线图Fig.6 Velocity contour and flow chart of inlet section

图7 空气滤清器结构参数示意图Fig.7 Schematic diagram of the structuralparameter of the air filter

从以上分析可知,空气滤清器上壳体涡系较为复杂,下壳体平均速度偏小,因此,进气效率低,进气阻力偏大.

2.2空气滤清器结构参数对进气阻力的影响

空气滤清器结构参数示意图,如图7所示.以下将主要从对上下壳体流动特性影响较大的进出气口结构参数入手,分析进出气口结构参数,如进出气口面积比(S1/S2)、进气口圆角半径(R1)和出气口圆角半径(R2)对进气阻力的影响规律,结果如图8所示.

2.2.1进出气口面积比保持进气口截面积S1不变的情况下,改变出气口S2的大小,分析空气滤清器在额定空气流量工况下进气阻力随进出气口面积比变化规律,如图8(a)所示.由图8(a)可知:随着进出气口面积比的增加,空气滤清器在空载和负载情况下进气阻力基本呈线性趋势增加.这是因为进出气口面积比越大,出气口面积相对较小,气流从出气口流出不顺畅.因此,在结构设计时应该尽可能加大出气口面积,提高空气滤清器进气效率.

由图8(a)还可知:进出气口面积比变化对滤芯的阻力基本无影响.此外,进出气口面积比由1.5减小为0.8,相应的出气口面积增加89.8%,则空气滤清器空载时进气阻力降低49.6%,负载时进气阻力降低37.2%.

(a) 进出气口面积比        (b) 进气口圆角半径        (c) 出气口圆角半径图8 空气滤清器结构参数对进气阻力影响Fig.8 Influence of structural parameters of air filter on air intake resistance

2.2.2进气口圆角半径保持进出气口面积比为定值1.5,改变进气口与壳体的圆角半径R1,分析空气滤清器在额定空气流量工况下进气阻力随进气口圆角半径的变化规律,如图8(b)所示.由图8(b)可知:随着进气口圆角半径的增加,空气滤清器在空载和负载情况下,进气阻力都随着原角半径的增加而降低;增加圆角半径可以避免气流从进气口流入壳体后出现局部的气流分离,提高了进气效率;当圆角半径达到一定值时,进气阻力不再减小,圆角半径的变化对滤芯阻力基本没有影响.此外,当进气口圆角半径由0增加到10 mm时,空气滤清器空载时,进气阻力降低20.4%;负载时,进气阻力降低16.6%.

2.2.3出气口圆角半径保持进出气口面积比为定值1.5,改变出气口与壳体的圆角半径R2,分析空气滤清器在额定空气流量工况下,进气阻力随出气口圆角半径的变化规律,如图8(c)所示.由图8(c)可知:随着出气口圆角半径的增加,空气滤清器在空载和负载情况下, 进气阻力都随着原角半径的增加而降低.增加圆角半径相当于增加出气口截面积,气流从出气口流出更加顺畅,提高了进气效率.当圆角半径达到一定值时,进气阻力不再减小,圆角半径的变化对滤芯阻力基本无影响.此外,出气口圆角半径由0增加到15 mm时,空气滤清器空载时,进气阻力降低42.9%;负载时,进气阻力降低36.2%.

表2 正交分析方案和进气阻力结果

2.3进出气口结构参数正交分析法

应用正交分析方法分析进出气口面积比(S1/S2)、进气口圆角半径(R1)、出气口圆角半径(R2)等3个影响因素的主次关系以及相互作用规律.选用L9(33)三因素三水平正交分析方案,水平选取遵循等间距原则,计算得到9种方案进气阻力(Δpr),如表2所示.

按极差分析法得到正交试验结果分析,如表3所示.表3中:Yj,k表示j影响因子的k水平所对应的进气阻力平均值,由Yj,k大小可以判断j影响因素的优水平;Rj为j影响因素水平变动时进气阻力的变动幅度,Rj越大,说明该因素对进气阻力影响越大. 由表3及前述分析可知: 组合优化方案比单因素优化方案效果更好.

表3 正交分析表

本次试验对进气阻力大小影响的主次因素依次是:进出气口面积比、进气口圆角半径、出气口圆角半径.最优组合是进出气口面积比(S1/S2)为0.8,进气口圆角半径(R1)为15 mm,出气口圆角半径(R2)为15 mm.

3结论

通过对乘用车常用的褶状干式空气滤清器进气阻力进行分析与研究,得出以下5点结论.

1) 空气滤清器进气阻力随着流量的增加呈抛物线趋势增加.利用最小二乘法拟合出空气滤清器进气阻力与流量的关系式,对结构和尺寸相似的空气滤清器初期结构设计,预测进气阻力大小具有一定指导意义.

2) 在流量一定的情况下,空气滤清器进气阻力主要来自于壳体自身,滤芯所占比重较小,随着流量增加,这种趋势越明显.因此,合理的结构设计至关重要,其中,进出气口结构参数影响最大.

3) 进气阻力随着出气口面积和出气口圆角的增加而降低.在结构设计时,可优先考虑增加出气口圆角半径,当圆角半径达到约15 mm时,进气阻力基本保持不变,空载和负载时,降低幅度分别达到42.9%和36.2%.

4) 进气阻力随着进气口圆角半径增加而降低.当圆角半径达到10 mm左右时,进气阻力基本保持不变,空载和负载时,降低幅度分别达到20.4%和16.6%.

5) 组合优化方案优于任意单因素优化方案,结构设计时可结合实际情况选择进出气口面积比及进出气口圆角半径.对进气阻力影响的主次因素依次为:进出气口面积比、进气口圆角、出气口圆角.

参考文献:

[1]REBAL M,PRAT M,MEIRELES M,et al.A semi-analytical model for gas flow in pleat filters[J].Chemical Engineering Science,2010,65(1):2835-2846.

[2]NAGARAJAN G,KUMAR S,CHOWDHURY D.CFD analysis of air filters for an off-highway vehicle[C]∥SAE Motosports Engineering Conference and Exhibition.Detroit:[s.n.],2007:559-567.

[3]刘联鋆,郝志勇,刘迟.空气滤清器流动阻力与噪声特性的仿真和优化[J].汽车工程,2011,33(12):1092-1097.

[4]LIU Chi,HAO Zhiyong,CHEN Xinrui.Optimal design of acoustic performance for automotive air-cleaner[J].Applied Acoustics,2010,71(5):431-438.

[5]关英俊,张力锋,李想,等.汽车空气滤清器的模态分析及结构优化[J].计算机仿真,2015,32(4):124-128.

[6]何志霞,蒋兆晨,王硕,等,空气滤清器的空气动力学仿真及优化[J].中南大学学报(自然科学版),2012,43(3):1179-1184.

[7]王鹏宇,康宁.某型号空滤器进气系统的气动和消声性能研究[J].汽车科技,2013(4):1-7.

[8]唐刚志,张力,邓涛,等.发动机用平板型空滤器流动阻力特性分析和改进[J].重庆大学学报,2015,38(4):18-23.

[9]赵树恩,蒋艾伶.某发动机空滤器的CFD优化设计[J].机械设计与制造,2015(8):163-166.

[10]重庆汽车研究所.汽车用空气滤清器试验方法:QC/T 32-2006[S].北京:中国计划出版社,2006:1-30.

(责任编辑: 钱筠 英文审校: 杨建红)

Analysis on Air Intake Resistance of Automotive Air Filter

YUAN Zhiqun1,2, HUANG Feijian1, LIU Jinwu1, XU Xi′an1

(1. School of Mechanical and Automotive Engineering, Xiamen University of Technology, Xiamen 361024, China;2. School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China)

Abstract:To reduce the air intake resistance of automotive air filter, a coupled numerical model of shell and filter element was established by using computational fluid dynamics technique. The influence of the inlet and outlet structure parameters on air intake resistance was analyzed by orthogonal analyses method, and the relationship between the air intake resistance and the mass flow rate was summarized. The reliability of the numerical calculation method is verified by the experiment of flow rate versus pressure drop in the filter element. The results showed that the air intake resistance increased at parabola trend with the increasing mass flow rate. Under rated flow conditions, the air intake resistance decreased with the increasing outlet area and the maximum drop was 37.2% under load conditions. The air intake resistance decreased with the increasing inlet and outlet fillet radius, the maximum drop was 16.6% and 36.2% under load conditions.

Keywords:air filter; air intake resistance; structure parameters; computational fluid dynamics; numerical analysis

中图分类号:TK 402

文献标志码:A

基金项目:国家自然科学基金资助项目(51405517, 51405411); 福建省教育厅科技项目(JB13151); 福建省厦门市重大科技计划项目(3502Z20100080)

通信作者:袁志群(1983-),男,讲师,博士研究生,主要从事车辆空气动力学的研究.E-mail:yzqhnu@163.com.

收稿日期:2015-11-23

doi:10.11830/ISSN.1000-5013.2016.03.0268

文章编号:1000-5013(2016)03-0268-05

猜你喜欢
数值计算
期权定价理论方法研究
“流动与传热数值计算基础”教学方式思考
阵列感应测井在直井和斜井中的对比
浅谈MATLAB在数学建模中的应用
数值模拟两层有限大小地层中多电极电流场分布
平衡流量计流动特性数值计算分析
MATLAB软件可视化效果和数值计算在高等数学学习中的应用
方形散流器喉部对送风气流均匀性的影响分析