电子膨胀阀射流噪声数值模拟
2014-06-05王志毅任夫磊程德威
黄 皓,王志毅,任夫磊,程德威
(1.浙江理工大学,a.机械与自动控制学院;b.建筑工程学院,杭州310018;2.奥克斯电气有限公司,浙江宁波315191)
电子膨胀阀射流噪声数值模拟
黄 皓1a,王志毅1b,任夫磊1a,程德威2
(1.浙江理工大学,a.机械与自动控制学院;b.建筑工程学院,杭州310018;2.奥克斯电气有限公司,浙江宁波315191)
应用Realizable k-ε湍流模型和Mixture两相流模型,对电子膨胀阀小开度下的二维流场进行了数值模拟,根据计算结果,使用宽频噪声计算方法得到了膨胀阀内流体噪声分布。对电子膨胀阀的改进模型进行了数值模拟,并对流体的速度场及其噪声辐射结果进行了分析,结果表明,阀口长度加长和阀口改为喇叭口结构可以减小射流噪声。
电子膨胀阀;射流噪声;FLUENT;数值模拟
0 引 言
电子膨胀阀是变频空调中必不可少的节流元件,通过与压缩机变容量技术的有效结合,可以有效提高变频空调系统的舒适性与节能性[1-2]。然而,当系统在特定工况运行时,电子膨胀阀会产生较大的噪声,有时甚至在室内都能听到。目前,对于电子膨胀阀噪声的研究很少,研究者主要是膨胀阀生产厂家以及空调制造商。珠海格力电器股份有限公司对冷媒流经电子膨胀阀时产生的节流噪声进行了研究,指出电子膨胀阀内气泡破裂噪声一般发生在1 000~2 000 Hz的中频区和5 000~6 000 Hz的高频区,并针对各种噪音提出了相应的控制措施[3]。广东美的制冷设备有限公司根据毛细管稳定性和准确性好的优点,在电子膨胀阀装机时,在其前后增加一定长度的毛细管,降低了电子膨胀阀阀口上下游冷媒的压力差,减小了冷媒的湍流强度,从而减小湍流噪声[4]。但是这些改进措施均基于特定的空调系统进行,没有提升电子膨胀阀自身的工作性能。电子膨胀阀产生的噪声主要为流体噪声,因此,从改进流道入手来解决电子膨胀阀的噪声问题是企业掌握核心竞争力的必然选择。
本文针对电子膨胀阀节流噪声问题,通过FLUENT软件开展流场和噪声分布的数值模拟,分析了膨胀阀的流道几何结构对流场和噪声的影响,对电子膨胀阀降噪具有一定的工程指导意义。
1 数学模型
由于电子膨胀阀内流体存在相变过程,本文采用Mixture多相流模型,同时开启空化模型。Mixture模型及空化模型的控制方程为:
a)连续性方程
b)动量方程
式中,ρ为流体介质密度,vi为介质在i方向上的速度,Xi为相应的空间i方向,Q为质量的产生率,Tij为广义应力张量,Fi为其他源项。而是引起射流噪声的主要原因,Tij包括三部分组成:
其中右端第一项为转移张量,第二项为粘滞应力张量,第三项为热传导引起的应力张量。在射流中,力的变化成为偶极子辐射源,应力变化则成为四极子辐射源。在亚声速和跨声速的自由喷射流中,粘滞应力与热传导的影响可以忽略不计,故稳态喷射流中唯一的噪声源来自于转移动量的梯度张量,这种声源一般为四极子声源[5-6]。因此,可以使用FLUENT噪声模型中的宽频噪声模型对电子膨胀阀射流噪声进行数值分析[7]。
c)紊流动能k及其耗散率ε的输运方程,
d)第二相的体积分数方程
关于流体噪声的理论研究,Lighthill最早给出了自由射流中湍流噪声总功率与喷射速度及出口直径的关系式:
其中,W 为湍流总声功率(W),ρ为流体密度(kg/m3),D为喷口孔径(m),v为流体速度(m/s),ρ0和c0分别为声传播介质密度(kg/m3)和声速(m/ s),K为Lighthill系数,实验值约为0.3×10-4~1.8×10-4。对于自由射流以及固体边界不起主要作用的喷射流,射流噪声受喷射速度的影响最大,阀口喷射速度越大,产生的湍流噪声也就越大。
Proudman[8]在Lighthill的噪声理论基础上,推导出一个各向同性湍流引起的声功率公式。Lilley[9]通过计算Proudman的最初理论中所忽略的变化较小的时间差分项推导出另一公式,得到单位体积的声功率为:
声功率级计算公式为:
式中,u为湍流速度(m/s),l为特征长度(m),a0为当地声速(m/s),α为一常数,Pref为参考声功率(W)。
2 电子膨胀阀流场数值模拟
2.1 几何建模及网格无关性分析
电子膨胀阀三维剖面图如图1所示,当流体从弯管流入、直管流出时为正向流动;当从直管流入、弯管流出时为反向流动,本文所有模拟均基于反向流动。为了缩短计算时间又能够体现流场特征,将此模型进行了简化,建立二维模型用于流场及噪声分析[10]。
在开度为200脉冲工况下,针对网格无关性验证,分别做了1万、2万和6万3种网格数量的模型,如图2所示。从图2中可见,阀口附近网格经过了加密处理,最小尺寸均小于10-5m。通过在流道出口处设置一个质量流量监测点,得到换算为圆管道后的流量分别为0.033 0、0.033 9、0.035 4 kg/s,其相对误差值均在5%以内,表明网格密度对计算结果的精度影响不大,本研究选择2万的计算网格。数值模拟所需的边界条件设置为给定进、出口压力条件,图形左侧显示为黄色的边界线设置为轴对称边界条件,其他均为无滑移壁面条件。
图1 电子膨胀阀三维剖面
图2 3种不同网格密度对比
2.2 流场模拟
数值模拟首先采用稳态求解器进行计算,然后以收敛的稳态解作为瞬态计算的初始值再次计算。计算过程采用的两相流模型为Mixture模型、湍流模型为Realizable k-ε模型、空化模型为Schnerr-Sauer模型。设定边界条件为进口压力2.5 MPa,出口压力0.8 MPa。流体物性参数为R410A两相参数,2.5 MPa条件下对应的R410A饱和液相密度为967.5 kg/m3、动力粘度为9.72×10-5kg/m·s,饱和汽相密度为105.67 kg/m3、粘度为1.4×10-5kg/m·s,定义液相为主相,汽相为第二相。求解方法设定为SIMPLE算法,松弛因子均设置在0.3~0.5之间。
根据计算结果,得到流体速度分布如图3所示。从图3可以看出喷口处速度最大,在喷口下游阀针壁面带拐角处速度梯度较其他地方大,此处湍流强度最高。从与其对应的湍流云图图4可以看出,此处即为射流的过渡区,其辐射出的噪声强度亦为最大,达到46.5 dB,如图5所示。将噪声计算结果与文献[11-12]对比后发现,模拟结果基本满足国标要求,因此,可认为该计算模型适用于电子膨胀阀的噪声模拟。
图3 电子膨胀阀流体等速线
图4 电子膨胀阀原始结构湍流动能分布
图5 电子膨胀阀流体噪声云图
2.3 流道模型改进及对比分析
电子膨胀阀产生射流噪声的主要原因是阀口前后压差大,过渡区内流体速度梯度大、湍流强度高,而湍流强度又与喷管形状和流体参数有关[13]。张勃等[14]通过实验,研究了不同宽高比的圆转矩形收敛喷管的射流宽、窄两个对称面上的湍流强度,分析结果表明,当宽高比大于8时,在射流中心线上距离喷口同一距离处,宽高比越大,湍流强度越大。这是因为喷口高度的减小,使得其对出口截面附近涡系发展的促进作用减弱。因此,本文尝试通过加长上游阀口长度,来促进涡系发展以减弱湍流强度。改进前、后模型的参数对比如表1所示。
表1 不同结构模型参数对比
将原始结构电子膨胀阀的数学计算模型,应用到改进后的电子膨胀阀上,得到了改进后的电子膨胀阀的模拟结果。阀口长度由2 mm加长为3 mm的湍流动能分布如图6所示,上游管道连接处由直角变为喇叭口的湍流动能分布(如图7所示)。与图4进行对比后发现,3种不同结构对应的混合区与过渡区的湍流分布均相似,改进后模型的湍流强度最大值均有不同程度的减弱,完全湍流区的湍流强度变大且占据更大的空间。当喷口加长后,出口处的涡系更容易控制在射流轴线范围内并向下游传播,射流结构中3个分区的连接更加紧密,湍流变化率趋于平缓,因此产生的湍流噪声也相应减小,从图8和图9可以明显发现这一变化。阀口加长后射流声功率减小了0.3 dB,喇叭口结构对应的声功率减小了1.2 dB,由此可以推断,加长阀口长度和阀口改为喇叭口结构可以减小射流噪声。
图6 加长阀口结构湍流动能分布
图7 喇叭口结构湍流动能分布
图8 加长阀口结构声功率分布
图9 喇叭口结构声功率分布
根据模拟的结果,得到3种模型喷口处的速度变化曲线,如图10所示。从图10中可以看出,两种改进后的电子膨胀阀,出口速度均比原始结构大,这是因为上游流道的改变为流体流动提供了较长时间的动力,所以,流体动量增加,速度变大。根据Lighthill理论,膨胀阀出口速度越大,射流噪声强度越大,噪声声压级越大,这与模拟结果正好相反。产生这种现象的原因,是Lighthill理论只适用于自由射流和固体边界不起主要作用的半受限射流,它无法对电子膨胀阀的射流噪声进行解释,这是由于电子膨胀阀的管道口径较小,内部流动很容易受到壁面的影响。因此,在进行理论计算时,必须将边界层的作用考虑在内。
图10 不同结构阀口速度变化对比
3 结 论
本文采用Realizable k-ε模型和Mixture两相流模型对小开度下3种不同结构的电子膨胀阀的流场进行了数值模拟,得到结论如下:
a)电子膨胀阀射流流场的湍动能分布与声功率级分布极为相似,验证了射流噪声强度主要取决于速度梯度这一理论。
b)Lighthill理论无法解释流场噪声分布和阀口处的速度分布对比结果,这说明对自由射流和固体边界不起主要作用的喷射流均适用的Lighthill理论不适用于口径较小的电子膨胀阀噪声计算。
c)通过宽频噪声模型计算得到的声功率分布图发现,加长阀口长度和阀口改为喇叭口结构的电子膨胀阀产生的射流噪声均比原始结构要小,即阀口长度加长和阀口改为喇叭口结构可以减小射流噪声。
[1]商萍君,易佳婷.电子膨胀阀的优势和发展趋势[J].制冷与空调,2007,21(2):72-76.
[2]闫长林.电子膨胀阀在家用变频空调器中的应用[J].家电科技,2013(3):52-54.
[3]何 珍.浅谈电子膨胀阀节流噪声的产生与控制[J].科技创业家,2013(4):79-81.
[4]林 竹,吴 空.变频电子膨胀阀空调冷媒流动异音分析研究[J].家电科技,2012(11):78-79.
[5]韩阳泉.喷注噪声的数值模拟[D].哈尔滨:哈尔滨工程大学,2007:45-48.
[6]马大猷.现代声学理论基础[M].北京:科学出版社,2004:94-95.
[7]安世亚太流体产品业务部.ANSYS Fluent航空气动噪声解决方案[J].航空制造技术,2012(9):103-104.
[8]Proudman I.The generation of noise by isotropic turbulence[J].Proc Roy Soc,1952,A214:119.
[9]Lilley G M.The radiated noise from isotropic turbulence revisited[R]//NASA Contract Report.Hampton:NASA Langley Research Center,1993:75.
[10]王志刚,王冬青,柳成文,等.电子膨胀阀内制冷剂两相流动过程的数值计算与实验验证[J].流体机械,1999,27(8):51-55.
[11]中华人民共和国国家经济贸易委员会.JB/T 10386-2002,家用和类似用途空调电子膨胀阀[S].2003.
[12]中华人民共和国工业和信息化部.JB/T 10212-2011,制冷空调用直动式电子膨胀阀[S].2011.
[13]Shih T H,Chen J Y,Lumley J L.Second order modeling of boundary-free turbulent shear flows[J].AIAA Journal,1992,30(6):1553-1560.
[14]张 勃,吉洪湖,曹广州,等.宽高比对矩形喷管射流湍流强度影响试验[J].航空动力学报,2010,25(10):2244-2248.
NumericaI SimuIation of Jet Noise of EIectronic Expansion VaIve
HUANG Hao1a,WANG Zhi-yi1b,REN Fu-lei1a,CHENG De-wei2
(1.Zhejiang Sci-Tech University,a.School of Mechanical Engineering&Automation;b.School of Civil Engineering and Architecture,Hangzhou 310018,China;2.AUX Electrical Appliance Co.,Ltd.,Ningbo 315191,China)
Numerical simulation of two-dimensional flow field was carried out for electronic expansion valve with Realizable k-εturbulent model and mixture two-phase flow model.Based on the calculation result,flow noise distribution in the expansion valve was gained with broadband noise calculation method. Numerical simulation of the improved model of the electronic expansion valve was conducted.Besides,velocity field of the flow and noise radiation results were analyzed.The results show longer valve port and trumpet structure of the valve port can reduce jet noise.
electronic expansion valve;jet noise;FLUENT;numerical simulation
O358
A
(责任编辑:康 锋)
1673-3851(2014)05-0512-05
2014-03-01
国家科技支撑计划项目(2012BAF01B05)
黄 皓(1988-),男,河南信阳人,硕士研究生,研究方向为流体机械理论及应用。
王志毅,E-mail:zywang@zstu.edu.cn