喷雾模型在cfd软件中的应用.docxVIP

喷雾模型在cfd软件中的应用 喷雾的子模型 如果在汽油车的工作过程中进行数值模拟，则燃料喷雾器的精确模拟对于预测燃烧和排放非常重要。对喷雾的模拟是通过各具特色的子模型来考虑的。气液动量交换、湍流扩散、液滴蒸发、喷雾破碎、液滴碰撞和喷雾碰壁等子模型组成了一个复杂的喷雾模型。本文主要讨论喷雾破碎模型的研究现状及其在CFD软件中的应用。 1 喷雾破碎模型 内燃机喷雾破碎机理复杂,大量的实验研究表明,在高韦伯数的喷嘴附近区域,喷雾破碎不仅受到气液两相间的相互作用的影响,而且受到喷嘴内部流动现象如湍流和空穴等影响;而在距喷嘴较远的喷雾下游区,喷雾破碎将完全受空气动力的作用,与喷嘴形式无关。由O’Rourke和Amsden提出的经典TAB破碎模型将喷雾液滴的振动和变形与弹簧质量系统相类比,建立液滴变形方程,通过变形量判断液滴破碎形式,适宜于低韦伯数的喷雾模拟,因为对于特别高的韦伯数,喷雾液滴散落分布,用弹簧质量系统类比是不适合的。由Reitz和Diwakar提出的GM模型将喷雾液滴破碎分为袋形破碎和剥裂破碎两种形式,通过临界韦伯数进行判断。这2种模型从描述单个液滴破碎的角度出发,建立喷雾破碎模型,可用于柴油机和直喷汽油机喷雾破碎的模拟,具有应用简单的优点。然而,由于这2种模型不能直接反映喷嘴内部流动对喷雾破碎的影响,同时射流破碎毕竟不同于液滴破碎,因此,在反映喷雾破碎机理上有着先天的不足。目前,表面波不稳定性理论是研究射流碎裂机理中最成功的理论,在内燃机喷雾模型中的应用也越来越多。WAVE模型认为KH不稳定波的增长导致了液滴从射流表面剥落下来,模型适宜于高相对速度和高环境密度的液滴破碎过程。在WAVE模型的基础上,发展了KH-RT模型。Huh-Gosman模型认为喷嘴内部湍流对喷雾破碎过程至关重要要。LISA模型认为KH不稳定波在液膜表面产生,其中短波控制表面波的增长,用于对产生空心锥喷雾的喷嘴附近雾场进行模拟。同时,考虑喷嘴内空穴流动对喷雾破碎影响的模型也越来越多,孔式喷嘴流动模型是其中的一种。本文以孔式喷嘴流动模型、WAVE模型、KH-RT模型、Huh-Gosman模型和LISA模型在内燃机CFD软件中的应用为出发点,对其进行讨论。 2 空穴产生的原因 大多数CFD软件将燃油喷雾视为由若干离散的具有代表性的统计样本——液滴组(particle)组成,在喷嘴出口处需要给定每个代表液滴的初始位置、大小、速度和温度。由于缺乏喷嘴附近的喷雾破碎的实验数据,传统地采用喷雾下游处实验数据作为喷嘴的入口条件,是一种不得已而为之的方法。对于孔式喷嘴喷雾,最常用的理论计算方法是采用离散的blob液滴模型,假定初始油滴的半径和速度分别为喷孔出口的几何半径和平均喷射速度。实际上,blob模型忽略了喷嘴内部流动对喷孔出口处流动状态产生的影响。在喷孔内可能发生的流动有包括层流和湍流在内的单相流动、空穴流动和介质倒流。现代柴油机中的高压喷射使得在喷孔中最容易发生空穴流动,如图1所示。图中R为喷孔入口半径,D为喷孔出口直径,L为喷孔长度。喷嘴流动模型在blob液滴模型的基础上,提出了有效喷射速度和有效喷孔面积的概念。 通过计算喷孔内入口收缩流动区域(C点)的静压力pc来判断流动是否处于空穴状态 pc=p1-ρl2U2cUc=UCcCc=1√1C2ct-11.4RD(1) p1=p2+ρl2(UCd)2(2) 其中Cd=1√Κinlet+fL/D+1 式中 ρl——燃油密度 Cd——流量系数 p1、p2——入口处(1点)、出口处(2点)压力 Uc——收缩处(C点)的流动速度 Cc——收缩系数,按Nurick拟合的式(1)计算 Cct——理论常数,约等于0.611 U——名义平均流速,通过式(2)计算 f——壁面摩擦因数 Kinlet——入口处(1点)的损失系数 当pc大于液体的饱和蒸汽压pvapor时,喷嘴内流动处于单相流状态,则设定液滴出口速度等于U,初始液滴索特平均直径D32等于喷孔出口直径;当pc小于pvapor,则设定流动处于全空穴状态,需要重新计算一个新的入口处压力p1和流量系数Cd p1=pvapor+ρl2U2cCd=Cc√p1-pvaporp1-p2 通过求解收缩口到出口之间的质量连续和动量守恒方程求得有效喷射速度Ueff和有效喷射面积Aeff为 Ueff=Uc-p2-pvaporρlUAeff=AUUeff 式中 A——喷孔截面积 因此,喷射液滴SMD为 D32=√4Aeff/π 喷嘴的几何形状影响着空穴的产生,R/D比值越低、喷嘴内入口越尖,空穴在孔内发生的位置将越早。越尖的喷嘴形成的有效速度越高,有效直径越小。L/D比值只影响空穴发生的时间,不影响空穴发生的程度。KIVA、FIRE和FLUENT均提供了空穴喷嘴流动模型,用于对高压柴油机喷雾液