液液同轴离心式喷嘴可以看做是两个离心式喷嘴组合而成,如图1所示^[3]。离心式喷嘴内流动的相关研究工作一般都是针对单个离心式喷嘴进行的,因此对于单离心式喷嘴,其内部流动特性研究非常深入,发表的相关文献很多^[4-12],有学者对其进行了深入细致的总结^{[1, 13]},本文不再赘述。对于拥有环形流道的外路离心式喷嘴的内部流动,学界关注较少,受喷嘴空间和试验观测手段限制,内部流动特性的研究以仿真居多。Vasque^[14]基于OpenFOAM和层流VOF模型分别开展了收口型的内喷嘴以及外喷嘴三维仿真,给定恒定的入口和出口压力,将仿真与实验结果进行对比,发现数值仿真能很好地预测内、外喷雾锥角,然而得到的流量系数与实验值存在一定误差。Ronceros Rivas^[7]等人则用VOF方法和RNG k-ε湍流模型对收口型内喷嘴和敞口型外喷嘴组成的液液喷嘴开展不同工况的数值仿真,结果表明仿真流场和试验吻合较好,而敞口型喷嘴和收口型喷嘴因为旋流室结构差异,导致内流场压力、速度上存在较大差异。汪凤山^[15]使用大涡模拟(LES)研究液液喷嘴流量系数、喷雾锥角以及外部液膜的雾化过程,给定入口和出口压力边界条件,发现旋转方向相同的内外两路同时喷射时,两路相互影响,雾化锥角较两路均有所增加,随着喷嘴等直段长度的增加,内、外喷嘴的液膜速度和流量系数均逐渐减小,喷嘴流量系数、喷雾锥角和出口液膜速度随着喷注压降的增加而增加。宋大亮^[16]采用VOF方法,将内喷嘴替换为堵头,改变堵头长度及直径。研究内喷嘴对收口型外喷嘴流动特性的影响,仿真得到的锥角和压降与试验吻合较好,显示了仿真的准确性,通过仿真得出内喷嘴对于外喷嘴流场有干扰,内喷嘴直径过大情况下,外喷嘴液膜流动并不会出现气核,应控制内喷嘴直径不大于外喷嘴出口气涡直径,这样设计出的喷嘴更优。吴宝元^[17]则研究了喷嘴几何特性和结构参数对敞口型喷嘴流量系数的影响,并对基于动量原理的流量系数进行了修正,提高了其精度。

对于液液同轴离心式喷嘴,根据喷嘴缩进长度的不同,整个喷嘴内外流动可以分为外混流动、临界流动和内混流动3种,如图2所示^[18]。喷嘴内部流动主要研究临界流动和内混流动两种状态下液膜在喷嘴内部的流动特性。内喷嘴缩进的液液同轴离心式喷嘴内部流动特性主要包括内喷嘴液膜撞击外喷嘴液膜、喷嘴内部液膜表面波的发展、液膜厚度的变化以及燃料与氧化剂的混合特性等。也有部分学者研究了内外喷嘴液膜互击、液膜厚度、

液膜表面波波长和频率等方面。随着喷嘴缩进长度的增加,内喷嘴产生的液膜会撞击外喷嘴液膜。Kim^[18]指出喷嘴缩进主要影响喷嘴内部液膜撞击点位置,有缩进时液膜在喷嘴内部混合,当内喷嘴液膜撞击外喷嘴壁面时会产生一种“撞击波”,由于黏性的存在,波的振幅随着缩进长度的增加呈指数递减,缩进越大,喷嘴出口这种波的振幅就越小。关于缩进室内部流场,受试验测量手段限制,目前试验研究较少,而仿真研究更多关注内外喷嘴单独喷注时的流场状况。Kim^[19]使用电导率法测量了8个不同缩进长度的组合喷雾出口液膜厚度,其中内喷嘴为收口型,外喷嘴为敞口型,他指出缩进会增加喷嘴的内混程度,缩进室内部液膜厚度存在低频振荡,这种低频振荡与质量流量的时空振荡以及缩进室内的混合过程有关,燃气发生器和燃烧室内的低频振荡和液膜之间的相互作用以及推进剂供应系统的振荡有关。反向旋转的液膜与同向旋转的液膜相比,内部混合增强(缩进增加)都使得液膜表面波频率增加,但是液膜表面波振幅却表现出相反的趋势。

液液同轴离心式喷嘴的外部流动主要包括内外液膜的相互作用(聚合与分离)、液膜的形态与破碎机理等。

理论分析是研究液膜破碎机理的重要手段,但由于其复杂性,当前的研究较少。丁佳伟^[20]在其博士论文中基于线性稳定性理论建立了描述液液同轴圆环旋转射流表面波稳定性的色散方程,实现了时间模式下该色散方程的数值求解,结果表明近反对称波形的扰动对液膜破碎过程起到支配作用,此外,作者还研究了流速比、气液密度比、表面张力以及液膜曲率等对射流稳定性的影响规律,为三维扰动下的液液同轴圆环旋转射流表面波稳定性研究奠定了基础。

不少学者研究了不同条件下的同轴射流喷雾形态,其中Sivakumar做了开创性的工作。研究者们使用高速相机拍摄得到瞬态喷雾图像,典型的内外喷嘴液膜聚合和分离的喷雾如图3所示。Ramamurthi^[21]采用试验研究了不同压降下的喷雾形态和破碎特性,发现郁金香形的喷雾形态对环境和喷注工况微小改变非常敏感,当喷嘴出口的离心力超过表面张力两个数量级,喷雾从郁金香形转变为分散的锥形。Sivakumar^[22]根据内外喷嘴液膜韦伯数的不同,将两液膜的流动状态分为聚合状态(状态Ⅰ)、分离状态(状态Ⅱ)、滞后状态(状态Ⅲ)以及自激振荡状态(状态Ⅳ),如图4所示,对于状态Ⅳ,Sivakumar指出该状态存在于外喷嘴韦伯数小于50的工况条件下,且振荡频率随内喷嘴韦伯数的增加而增大。

Kim^[18]通过改变缩进长度从而导致3种不同的液膜混合方式:外混、临界混合和内混,他研究了不同缩进状态下工况参数对液膜流动形态的影响,如图5所示。其中,对于无缩进的外混流动,可以明显看到同轴喷雾从分离到聚合的状态。而对于有缩进的内混流动,同样的韦伯数下,液膜始终处于聚合状态。

值得注意的是,液液喷嘴的自激振荡现象鲜有学者研究,而一般认为自激振荡与燃烧不稳定关系密切,因此值得深入探讨。对于内外液膜的聚合和分离过程,Sivakumar^[23-24]指出即使没有发生聚合,内、外锥形液膜之间也存在明显的相互作用,推测两个液膜之间的压力波动是造成这种现象的主要原因,液膜的聚合是由于内、外液膜间负压区的存在而形成,这一点在其他学者的仿真研究中得到证实。此外,液膜的聚合与分离过程存在滞后性,该滞后现象会强烈影响燃烧效率。

喷雾锥角是液液同轴离心式喷嘴外部最直观的流动特性,锥形液膜的聚合和分离对喷嘴整体的喷雾锥角具有重要影响。内外两个喷雾相互作用形成聚合喷雾,聚合喷雾的喷雾锥角介于内喷嘴和外喷嘴的喷雾锥角之间,锥角的大小与内外路射流的动量比有关^[25]。而国内外众多学者的仿真和理论分析表明,理论与仿真得到的喷嘴整体的喷雾锥角,与实验结果吻合较好^{[13, 25-28]},研究表明液体物性和几何参数对于喷雾锥角有重要影响。当忽略了黏性影响,理论喷雾锥角值与试验值往往有较大误差。

同离心式喷嘴一样,液液同轴离心式喷嘴的喷雾锥角受工况和结构参数的影响显著。工况参数方面,Rashid^[29]发现喷注压降相同时,外喷嘴喷雾锥角比内喷嘴的喷雾锥角更大。徐顺^[30]、王尧^[31-32]发现内外喷嘴单独工作时,喷雾锥角随喷注压降增加而增加; 内外喷嘴同时工作时,喷雾锥角随内喷嘴喷注压降的增加而减小,随外喷嘴喷注压降的增加而增加,如图6所示。Eberhart^[33]用冷流实验研究了不同反压对喷雾锥角的影响,结果表明,喷雾锥角随室压增加而减小,当燃烧室压力增大时,测量得到的喷雾锥角明显偏离理论值。结构参数方面,Kim^[18]分析了喷嘴缩进对喷雾锥角的影响,发现随着缩进长度的增大,喷雾锥角减小(图7)。Alves^[34]研究了切向孔数量对液液喷嘴的雾化和混合特性的影响。发现切向孔增多,喷雾分布更加均匀,而混合效率和喷雾锥角基本不变。

液膜破碎是液液同轴离心式喷嘴外部流动特性的重要方面,而破碎长度是液膜破碎最重要的特征参数。在工况参数的影响方面,Rashid^[29]研究了喷注压降对液液同轴离心式喷嘴破碎长度的影响,指出随着喷注压力的增加,液膜破碎长度变小。Eberhart^[33]用冷流实验研究反压影响,结果表明液膜破碎长度与室压成反比。Ahn^[35]发现对于同轴喷嘴,反压越大,破碎长度越小,喷注压降越大,破碎长度也越小。Kim^[18,36-37]指出液膜撞击点离喷嘴出口越远,破碎长度也就越大。对于内混状态下的喷嘴,因为黏性的存在,使得撞击波的能量随着缩进的增加逐渐减小,造成撞击波的振幅减小,液膜受到的扰动减小(图8)。Kim还发现随着内喷嘴韦伯数增加,无论是否存在缩进,内喷嘴破碎长度均减小,而随着外喷嘴韦伯数增加,破碎长度先迅速增加,到一个最大值后又逐渐减小。丁佳伟^[38]研究了喷嘴直径及长度等结构参数对液液同轴离心式喷嘴雾化特性的影响,发现增大外喷嘴直径,减小外喷嘴及内喷嘴长度可以减小射流破碎长度。

雾化特征包括SMD和速度等,均是评价喷嘴性能的重要标准。

SMD(式中简记为S)常被用来反映喷雾场的平均粒径,其定义为:S=∑D³_iN_i/(∑D²_iN_i),式中:D_i为液滴直径; N_i为直径为D_i的液滴数量。内外锥形液膜的聚合和分离会显著影响喷嘴的SMD。Sivakumar^[39]发现液膜聚合与分离过程将导致喷雾参数的急剧变化。聚合后的液膜将比分离状态下平均粒径大40%～50%,内部流速恒定时,组合喷雾的SMD随外部液膜喷射速度的增加先增大后逐渐降低。SMD受喷嘴结构、工况参数影响很大。Kim^[37]使用不同流体分析喷嘴缩进对喷雾特性的影响,如图9所示。指出随着缩进长度的增大,SMD先增加后减小,不同介质得到的最大SMD差异明显。丁佳伟^[38,40]使用VOF耦合大涡模拟对液膜破碎过程进行数值仿真研究,仿真结果表明,增大外喷嘴直径,减小外喷嘴及内喷嘴长度有利于增强雾化,而减小内喷嘴直径SMD增加。同轴射流的平均直径在内喷嘴射流和外喷嘴射流之间。和单独的离心式喷嘴相比,同轴射流液滴尺寸分布更广泛。

由于内外锥形液膜的聚合和分离与喷嘴的喷注压降等工况参数密切相关,因此工况参数对喷嘴雾化特性的影响规律尤为重要。Soltani^[41]研究了不同质量流量下喷嘴的喷雾特性,指出内喷嘴流量对SMD具有显著影响,随着喷雾远离喷口,最大SMD向外周移动并不断增大。丁佳伟^[42]研究了喷注压降对喷雾的影响,发现表面波频率随着喷注压降增加而增加。单独的内喷嘴表面波频率比同轴喷雾更高,说明更不稳定,因此更容易破碎,而单独的外喷嘴喷雾是最稳定的,更进一步指出SMD随着内喷嘴压降增加而减小。内喷嘴对同轴喷雾SMD的影响更大。对于内喷嘴压降较低时的同轴喷雾,液膜厚度是SMD的主导因素。郭志辉^[43]对比了脉冲和稳态喷雾,发现脉冲的存在导致雾化恶化,其喷雾中心区域直径分布不均,SMD在某个值附近脉动。脉冲时间减小,粒径的不均匀性增加(图 10)。Eberhart^[33]用冷流实验研究了不同反压对喷雾特性的影响。结果表明,当燃烧室压力增大时,一次雾化区域内的液滴直径与室压成反比,如图 11所示。 也有一些学者使用理论公式、仿真等方法对SMD进行研究,Han^[44]将Lefebvre和Jone的SMD经验公式结合得到的SMD与测量值很接近。Ghorbanian^[45]采用广义回归神经网络来预测喷嘴任何轴向或径向位置处的SMD,结果表明,神经网络预测的结果与实验数据吻合较好。

速度分布也是雾化特性的重要方面,Soltani^[41]研究了不同质量流量下喷嘴的喷雾特性,内喷嘴流量对速度具有显著影响,随着喷雾远离喷口,液滴最大速度减小。脉动速度的分布特点是相同的,但脉冲喷雾中心区域液滴脉动速度大于平均速度。Eberhart^[33]用冷流实验研究了不同反压对喷雾特性的影响。结果表明,当燃烧室压力增大时,一次雾化区域内液滴的平均速度与室压成反比。Ghorbanian^[46]采用广义回归神经网络来预测喷嘴任何轴向或径向位置处的液滴速度,结果表明,神经网络预测的结果与实验数据吻合较好。郭志辉^[43]发现稳态喷雾切向速度小,轴向为典型的双峰分布,如图 12所示。

混合特性也是喷雾特性的重要方面,对燃烧过程有重要影响。Radke^[47]向水里面添加碘化钾配成溶液,用CT进行喷雾质量分布研究,发现喷雾锥明显不均匀,由几个质量浓度随喷注压力变化的小块组成。Kim^[48]研究了缩进对混合和雾化的影响,如图 13所示,发现在轴向位置超过3 cm后,内混的混合效率比临界混合和外混的效率高10%,而外混有更好的雾化,因此在混合效率和雾化性能上存在一个折中方案。他认为实际用于火箭发动机的喷嘴,混合效率比雾化质量重要,因为如果雾化不好,大液滴可以在高温和高压环境下燃烧掉,因此内混是最优的,但也有学者认为临界混合是最优的。Han^[44]发现质量分布和混合效率对缩进比SMD更敏感,缩进5.25 mm时,有最大混合效率; 同时也指出存在一个最佳的喷嘴缩进,它就是使混合效率最高的时候。Woodward^[49]指出采用喷嘴缩进形成的缩进室作为混合室是解决液氧性能不好的合理方式。大缩进的喷注器由于混合良好使得喷注器性能显著提高并且没有造成结构损失和较大的压力损失。Seol^[50]研究了缩进对雾化和混合的影响。混合效率达到最大时的缩进对SMD的影响不大,然而,质量分布和混合效率改变了很多。因此可以推论,最佳缩进数存在,就是当混合效率达到最大的时候。除缩进外,Alves^[34],Long^[51]研究了切向孔个数对推进剂混合比的影响,发现切向孔数量增加,提高了喷雾周向质量分布均匀性,而混合效率没有随着切向通道数量的增加呈现出显著的变化。