本文在前人研究的基础上，研究了甩油盘内外部的非均匀流动现象，使用数值模拟及实验的方法，对甩油盘的流动非均匀性进行了研究，同时初步尝试进行理论上的解释。
1 数值计算方法 本文使用ANSYS FLUENT商用仿真软件对直连式甩油盘的全盘内部流域及盘外雾化区域使用雷诺平均Navier-Stokes(RANS)方程的方法进行求解。该程序使用旋转坐标系的方式耦合离心力与表面张力，分别使用VOF连续相和DPM离散相的方法对甩油盘内部流动及甩油盘外雾化场进行分析，使结果更加精确。
本文仿真工作分为2个部分：①甩油盘内部的流场；②单一甩油盘孔及全盘中的流体离开甩油盘后由于高速运动与静止空气互相作用的流场。介质设置为水，2个部分分别有不同的流场域及不同的3D模型。
1.1 甩油盘内部流动的数值计算方法 内部流动采用瞬态计算，计算方法如表 1所示。
表 1 内部流动数值计算方法 Table 1 Numerical simulation methods for internal flow
| 计算方法类型 | 计算方法具体内容 |
| 计算域 | 全盘 |
| 网格划分方式 | 非结构化网格 |
| 湍流模型 | k-episilon |
| 两相流模型 | VOF |
| 进口边界条件 | 流量进口 |
| 出口边界条件 | 压力出口 |
| 甩油盘壁面边界条件 | 无滑移旋转壁面 |
表选项
1.1.1 计算域 本文采用计算完整甩油盘的方法对某型甩油盘进行了数值模拟。所仿真的甩油盘是某具体型号发动机上的甩油盘，总转速为22 000 r/min，下文使用百分比表示转速。其结构如图 1和图 2所示。
 |
| 图 1 甩油盘剖视图 Fig. 1 Cross-section of fuel slinger |
| 图选项 |
 |
| 图 2 甩油盘出口孔位置示意图 Fig. 2 Schematic of fuel slinger exit hole location |
| 图选项 |
图 1和图 2给出了甩油盘的具体结构，其组合件由连接管、甩油盘和后盖板焊接而成，在仿真中将它们共同形成的内腔视作第1部分的仿真域。甩油盘有18个孔，分为2排。离连接管较近的一排孔间距为20°或者40°，而距离连接管较远的一排孔间距为60°。
本文甩油盘的计算域具有几何对称性，但是并不选取其中的一部分进行计算，而是将甩油盘焊接组装件内部的整个空腔作为计算域。这是因为考虑到甩油盘内的液膜十分薄，即使小扰动也可能有十分重要的影响，甚至可能会有全盘范围的不稳定流动现象发生。
1.1.2 湍流模型及两相流模型 根据Wardle等[9]的研究，气液两相在几何对称的计算域内捕捉不均匀流动使用k-ε湍流模型、LES和DES(一种LES和RANS的混合模型)都能获得较为贴近的结果。但是考虑到后两者的计算对资源的消耗太过巨大，且本文使用瞬态方法计算两相流，本身计算难度较大，所以使用k-ε湍流模型仿真甩油盘内部流动。
对于捕捉甩油盘内部气液交界面的分布情况，考虑到甩油盘中液体多为连续的流体，所以采用欧拉方法的VOF法追踪单个网格中流体占网格的总体积比函数F来进行气液交界面的捕捉。若F=1，则该网格充满流体，若F=0，则该网格没有流体，若0 < F < 1，则流体占据部分网格。根据张健等[10]的总结，二维网格中VOF法的示意图见图 3，计算公式如下:
 | (1) |
 |
| 图 3 VOF法示意图 Fig. 3 Schematic of VOF method |
| 图选项 |
其输运方程(守恒形式)为
 | (2) |
采用如图 3所示的网格，阴影为液体，采用一阶差分对时间进行处理，可将式(2)转化为
 | (3) |
式中：
 | (4) |
1.1.3 网格划分和边界条件 由于甩油盘形状复杂(见图 4)，并使用标准壁面函数和无滑移壁面条件，第1层的网格高度为0.1mm，保证y+为50，使得FLUENT中的k-ε模型有效。相邻网格的膨胀比为1.1，网格总数为4×106。
 |
| 图 4 内部流场的网格划分 Fig. 4 Mesh of inner flow field |
| 图选项 |
连接管的进口边界条件为质量流量进口，具体流量为该小型航空发动机各个工况下的流量。出口为压力出口，压力设置为一个大气压。其余各甩油盘面均设置为旋转壁面。整个甩油盘内部的计算域按甩油盘中心轴根据每个工况给定旋转速度，壁面设置为与计算域旋转速度相同。
1.2 甩油盘外部雾化区域的数值计算方法 甩油盘外部单孔雾化区域采用稳态加DPM瞬态的方法计算，计算方法如表 2所示。
表 2 外部雾化区域的数值计算方法 Table 2 Simulation methods for outer atomization field
| 计算方法类型 | 计算方法具体内容 |
| 计算域 | 甩油盘外流场 |
| 网格划分方式 | 结构化网格 |
| 湍流模型 | k-episilon |
| 两相流模型 | DPM |
| 进口边界条件 | 速度进口 |
| 射流进口边界条件 | 压力雾化喷嘴 |
| 出口边界条件 | 速度出口 |
| 甩油盘壁面边界条件 | 无滑移旋转壁面 |
表选项
1.2.1 计算域 甩油盘的孔有28.5°的倾斜度，甩油盘甩出的水会集中成一根水柱，可认为该水柱是一个与其直径相同但扩散角为零喷雾张角也为零的压力雾化喷嘴喷出的水柱。
本文中单孔的液柱喷出后流场的计算域为一个内半径为70mm、外半径为140mm、扇形角为90°、厚度为60mm的扇形体，如图 5所示。其设计是根据发动机的尺寸决定的，在此尺寸下液体不会沿径向飞出燃烧室。
 |
| 图 5 单孔外流场的仿真模型 Fig. 5 Model to simulate outer flow field (one hole) |
| 图选项 |
根据甩油盘内部流动的结果，甩油盘外部流动也应该是不均匀的，所以在此对甩油盘的全盘外部流动进行仿真模拟，观察不均匀性对甩油盘的雾化场造成的影响。同时判定单孔计算的结果是否能够在一定意义上代表甩油盘全盘的雾化场情况。
使用直径远大于甩油盘燃烧室的环形仿真区域对甩油盘外的雾化场进行计算。同时，此仿真区域的直径大于实验的甩油盘保护罩，确保甩油盘保护罩内的雾化场能够得到有效模拟，该区域如图 6所示。此区域外直径为1400mm，内直径为甩油盘的外径，厚度为120mm，认为甩油盘孔在厚度方向上居中。
 |
| 图 6 甩油盘外流场全盘模拟区域 Fig. 6 Field to simulate outer flow field of fuel slinger (whole slinger) |
| 图选项 |
1.2.2 湍流模型及两相流模型 根据黄勇等[11]的研究，液柱在喷出甩油盘后，油雾将会变得十分细碎，占网格的比例越来越小，甚至小到被系统忽略的程度。所以，该情况下再使用VOF法就会导致较大的误差，解决方法是使用DPM法对喷雾进行模拟[12]。在前人研究的基础上[12-14]，本文选用k-episilon湍流模型。
1.2.3 网格划分和边界条件 单孔外部流动使用结构化网格，网格总数为102万。第1层网格高度为0.067mm，膨胀比为1.2，保证y+为50左右，使得使用壁面函数的k-episilon湍流模型有效。
全盘的外流场被划分为388万个结构化网格。第1层网格高度为0.067mm，膨胀比为1.2，保证y+为50左右，使得使用壁面函数的k-episilon湍流模型有效。
数值模拟中，压力雾化喷嘴与甩油盘孔位置相同。雾化空间模型中，除底部为甩油盘外壁，设定为壁面边界条件，其他面均设定为速度进口边界条件，进口速度相对于雾化区域为零。雾化区域以与甩油盘相反方向、大小相同的转速旋转。
2 数值计算结果 2.1 内流场仿真结果 图 7为70%转速时(总进流量与转速对应)甩油盘的总出流量与时间的关系。图中，红线表示设计总出流量，负值代表出流。可见，甩油盘的总出流量与时间相关，在0.54~0.64s内，最小总出流量为0.03kg/s，最大总出流量为0.122kg/s，差值为0.092kg/s，占最大总出流量的75.4%，具有强烈时间相关性。选取总出流量稳定在设计出流量时的结果对各孔的流量进行时间平均。
 |
| 图 7 70%转速下甩油盘总出流量与时间的关系 Fig. 7 Total outmass-flow rate versus time under 70% rotation rate of fuel slinger |
| 图选项 |
在进行时间平均后，取各转速孔流量的最大值、最小值、中位数和平均值，可得表 3。其中，与孔平均流量差值的表达式为
 | (5) |
表 3 不同转速下孔流量统计 Table 3 Statistics of mass flow rates of holes under different rotation rates
| 转速/% | 孔流量状态 | 孔流量/(kg·s-1) | 与孔平均流量差值/% |
| 60 | 最大值 | 6.72×10-3 | 88.24 |
| 中位数 | 3.72×10-3 | 4.20 | |
| 平均值 | 3.57×10-3 | 0 | |
| 最小值 | 4.68×10-4 | 86.89 | |
| 70 | 最大值 | 7.60×10-3 | 112.89 |
| 中位数 | 5.44×10-3 | 52.38 | |
| 平均值 | 3.57×10-3 | 0 | |
| 最小值 | 1.20×10-3 | 66.39 | |
| 80 | 最大值 | 2.35×10-2 | 263.78 |
| 中位数 | 5.76×10-3 | 10.84 | |
| 平均值 | 6.46×10-3 | 0 | |
| 最小值 | 4.02×10-3 | 37.77 | |
| 90 | 最大值 | 1.91×10-2 | 120.05 |
| 中位数 | 7.75×10-3 | 10.71 | |
| 平均值 | 8.68×10-3 | 0 | |
| 最小值 | 1.36×10-4 | 98.43 | |
| 100 | 最大值 | 2.00×10-2 | 33.33 |
| 中位数 | 1.70×10-2 | 13.33 | |
| 平均值 | 1.50×10-2 | 0 | |
| 最小值 | 1.00×10-2 | 33.33 |
表选项
式中：
转速为80%时，孔流量最大值与其对应转速下孔流量平均值的差值最高，达到了263.78%。孔流量最大值、最小值较平均值均至少相差33.33%，普遍接近或达到100%。
各孔流量情况如图 8所示。除100%转速外，各转速都存在流量向某几个孔集中的现象。转速为60%时，流量集中在孔1~孔9；转速为70%时，流量集中在孔8~孔18；转速为80%时，流量集中在孔17；转速为90%时，流量集中在孔8~孔15;转速为100%时，孔的流量呈现近正弦曲线型，同样也不遵循该甩油盘两排孔不同流量但同一排孔流量相同的设计初衷。说明甩油盘的内部流动具有不均匀性。
 |
| 图 8 流量与出口孔关系 Fig. 8 Mass flow rates versus exit hole |
| 图选项 |
2.2 内流场仿真结果原理 甩油盘中，油膜厚度最大不超过0.8mm，接近喷油孔时，油膜厚度为3~5μm[6]，如此薄的油膜符合Huppert[15]的情况。在重力场的作用下，液膜沿斜面向下运动，逐渐分离，形成如手指状的结构。原因是黏性力和表面张力的共同作用。同时，干燥的区域会在下游形成[16]。在液膜不够厚时，其能维持原状，这是因为它能在σls(液体与固体之间的表面张力)、σsg(气体与固体之间的表面张力)及σ(气体与液体之间的表面张力)作用下保持平衡[17]。
甩油盘中，离心力取代了重力的作用，也能够导致如图 9所示(越红代表水的体积分数越高)的螺旋手指状结构，这是甩油盘内流量不均的成因。
 |
| 图 9 甩油盘中形成的液面分离现象 Fig. 9 Water film split up inside fuel slinger |
| 图选项 |
在图 9中划分剖视图可得图 10，其与图 11中Huppert[15]通过实验验证的理论模型滞止点分离图十分近似，仿真结果符合其实验及理论分析。
 |
| 图 10 甩油盘A-A剖视图 Fig. 10 A-A cross-section of fuel slinger |
| 图选项 |
 |
| 图 11 Huppert的液面滞止点理论图[15] Fig. 11 Huppert's theory about liquid film stagnation point[15] |
| 图选项 |
2.3 单孔外流场计算结果 图 12为不同转速下孔流量最大值、中位数、平均值和最小值的SMD值。可见，转速是影响SMD最重要的变量。但同一转速下各孔不同的流量确实对SMD造成了影响。100%转速的SMD为42μm左右，60%转速的SMD为90μm左右，两者相差约100%。同时，同一转速下，70%转速时不同孔流量SMD差值达到最大，为11.17μm，相比于转速对SMD的影响小。
 |
| 图 12 SMD与孔流量统计值的关系 Fig. 12 SMD versus mass flow rate statistic value of holes |
| 图选项 |
2.4 全盘外流场计算结果 图 13和图 14分别为90%转速时使用由内流场仿真得到的孔流量及使用孔流量平均值仿真得到的结果。其他转速下液雾分布图差距不大，所以只取这2张图进行讨论。虽然在2.3节中可以得到一个结论：对于SMD的大小而言，流量的变化影响并不是很大。但是对于整个雾化场的分布而言，孔流量不均匀的影响是很大的。图 13在90%转速下，左下角的几个孔由于流量不大，外流场的雾化出现了空白区域，这可能对之后的燃烧影响很大。
 |
| 图 13 90%转速下使用内部流场得到的各孔流量仿真雾化场图 Fig. 13 Outer flow field simulated using different hole mass flow rates from inner flow field under 90% rotation rate |
| 图选项 |
 |
| 图 14 90%转速下使用孔流量平均值得到的各孔流量仿真雾化场图 Fig. 14 Outer flow field simulated using average hole mass flow rate under 90% rotation rate |
| 图选项 |
3 实验系统 实验系统由测试系统(高速摄像机、PDPA)、辅助系统(供水系统、控制系统、通风系统等)和旋转系统(高速电机及其附件)组成，如图 15所示。由高速电机带转的甩油盘工作转速是连续可调的。使用计算机控制工况参数，采集、记录及处理实验数据，监控实验过程等。图 16为实际实验系统图。
 |
| 图 15 实验系统示意图 Fig. 15 Schematic of experimental system |
| 图选项 |
 |
| 图 16 实际实验系统图 Fig. 16 Photo of experimental system |
| 图选项 |
试验台运行的原理为：经离心泵加压，供水系统的储水箱中的水经过供水管路被输送至高速电机的尾部。高速电机尾部与供水管路末端使用高速旋转接头连接，由高速旋转接头供水管路中的水进入高速电机空心的旋转主轴之内。到达高速电机旋转主轴前端的水，再通过特制的连接管从甩油盘中心进入到甩油盘内部。工作时，高速电机的旋转主轴和甩油盘同步旋转，在离心力的作用下进入到甩油盘内腔的水高速离开甩油盘各喷油孔，并形成液雾雾化。控制系统实现连续调节实验工况参数。使用调节阀和控制水泵可将流量连续调节，控制风机变频箱可将风机功率连续调节，控制高速电机可将电机主轴(即甩油盘转速)连续调节。控制系统还具有远程控制PDPA和高速摄像机的功能。实验时，改变供水流量和电主轴转速，应用PDPA和高速摄像机能够测试甩油盘在不同工况下的雾化性能，可通过结果分析流量和转速等对于甩油盘雾化性能的影响规律。
试验台及高速摄像机的性能见表 4和表 5。
表 4 试验台的技术指标 Table 4 Technique parameters of experiment rig
| 指标 | 内容 |
| 实验介质 | 水 |
| 介质温度 | 常温 |
| 供水压力 | 0 ~0.6MPa |
| 介质清洁度 | 不低于6级(GJB 420B—2015) |
| 流量调节精度 | ±0.2% |
| 供水流量 | 0~1200kg/h(连续可调) |
| 甩油盘转速范围 | 0~22000kg/min(连续可调) |
| 最大振动限定值 | ≤1g |
| 转速调节精度 | ±100r/min |
| 加速时间 | 从10500~22000r/min≤10s |
| 连续工作时间 | ≥1h |
表选项
表 5 高速摄像机性能 Table 5 Functionality of high-speed photography
| 参数 | 数值 |
| 最大分辨率/(像素×像素) | 1280×1024 |
| 最大拍摄频率(@分辨率1280×16)/(帧·s-1) | 62000 |
| 最大拍摄频率(@最大分辨率)/(帧·s-1) | 2550 |
| 工作温度/℃ | -40~50 |
| 传感器类型 | CMOS |
| 传感器规格/in | 1 |
| 注：1in=2.54cm。 | |
表选项
4 实验结果 如图 17所示，转速为1000r/min，水的质量流量为935.8kg/h，雷诺数Re=2.1×104工况下，甩油盘各孔流量较为均匀。图 17能够很好地反映甩油盘内部的加工状况，作为对比实验排除因为甩油盘内部加工误差导致各孔流量不均的状况。甩油盘内部不存在巨大的、能够几乎封死一个孔的加工误差，每个孔的流量十分均匀，流道通畅。推测在这样小转速的情况下，孔的流量能够均匀是因为甩油盘内部离心力的作用较小，而表面张力作用较大，液膜依旧能够保持比较好的环形，这样各孔都能持续地接触到液膜，进而流量较为均匀。
 |
| 图 17 甩油盘低转速下雾化形态 Fig. 17 Slinger atomization field under low rotation rate |
| 图选项 |
图 18和图 19为同一个实验工况，转速为6000r/min，水的质量流量为929kg/h，Re=13.1×104。这是2个不一样的孔，当它们转到同一个位置时，流量是不一致的。图 18中几乎看不到孔的流量，整个孔仿佛被堵住了一般，几乎完全不出水。图 19中可以观察到一条明显白色的液柱离开了甩油盘孔。这证明了2.1节中的结论，即甩油盘的内部流动具有不均匀性，其会影响甩油盘的喷雾均匀性。
 |
| 图 18 各孔流量对比(1) Fig. 18 Comparison of hole's mass-flow rate (1) |
| 图选项 |
 |
| 图 19 各孔流量对比(2) Fig. 19 Comparison of hole's mass-flow rate (2) |
| 图选项 |
如图 20所示，转速为6000r/min，水的质量流量为929kg/h，Re=13.1×104，很明显地展示出了在同一时间各孔有不均匀流量的情况。从左向右第2个红色箭头处流量十分大，而另外3个箭头所标的孔处流量几乎可以忽略不计。如果仔细观察全盘的流量，除了这4个非常明显的流量不一致的孔，可以发现其他孔流量也略有区别。
 |
| 图 20 各孔流量对比(3) Fig. 20 Comparison of hole's mass-flow rate (3) |
| 图选项 |
5 结论 本文针对工程应用的直连式甩油盘雾化装置的雾化过程开展了三维数值模拟及高速摄影实验研究。仿真中，内部流场为非稳态模拟，外部流场为稳态模拟。在本文几何及环境条件下，得到初步的分析结果。结果表明：
1) 甩油盘在最大转速为22000r/min，实际转速为60%~100%的条件下，液体流动不稳定，出现流量随时间变化的现象，总出流量为时间的函数，总出流量差值可达75.4%。
2) 甩油盘各孔出流量不一致，相差可达263.78%，这可能是由薄液膜在离心力和表面张力的作用下产生的分离现象所导致。
3) 转速相比各孔不同流量SMD影响更大。100%转速的SMD为42μm左右，60%转速的SMD为90μm左右，相差接近100%转速SMD的100%。同时，同一转速下，70%转速时不同孔流量SMD差值达到最大，为11.17μm，相比于转速对SMD的影响小。
4) 甩油盘雾化场的空间均匀性受到内部流动不均匀性的影响较大，液雾分布相当不均匀。
5) 高速摄影实验结果支持了以上结论，即孔流量确实是不均匀的。
参考文献
| [1] | MORISHITA T.A development of the fuel atomizing device utilizing high rotational speed[C]//ASME 1981 International Gas Turbine Conference and Products Show.New York: ASME, 1981: V001T04A007. |
| [2] | MAZALLON J, DAI Z, FAETH G M. Primary breakup of nonturbulent round liquid jets in gas crossflows[J]. Atomization and Sprays, 1999, 9(3): 291-312. DOI:10.1615/AtomizSpr.v9.i3.40 |
| [3] | DAHM W J A, PATEL P R, LERG B H.Visualization and fundamental analysis of liquid atomization by fuel slingers in small gas turbine engines: AIAA-2002-3183[R]. Reston: AIAA, 2002. |
| [4] | SALLAM K A, AALBURG C, FAETH G M. Breakup of round nonturbulent liquid jets in gaseous crossflow[J]. AIAA Journal, 2004, 42(12): 2529-2540. DOI:10.2514/1.3749 |
| [5] | 贾永忠. 甩油盘的雾化分析[J]. 推进技术, 1995, 16(4): 45-49. JIA Y Z. An analysis on atomization of fuel slinger[J]. Journal of Propulsion Technology, 1995, 16(4): 45-49. (in Chinese) |
| [6] | 宋双文, 罗辉.离心甩油盘特性浅析[C]//中国航空学会(动力)第三届小发动机, 中国航空学会(直升机)第三届动力与传动学术讨论会, 2000: 125-133. SONG S W, LUO H.The characteristics of the slinger injector[C]//The 3rd Small Engine Academic Symposium of Chinese Society of Aeronautics and Astronautics Chinese Society of Aeronautics and Astronautics, 2000: 125-133(in Chinese). |
| [7] | 熊纯, 宋双文, 陈剑, 等. 一种离心甩油盘雾化性能的试验[J]. 航空动力学报, 2012, 27(7): 1517-1522. XIONG C, SONG S W, CHEN J, et al. Experimental study of atomization performance on a fuel slinger[J]. Journal of Aerospace Power, 2012, 27(7): 1517-1522. (in Chinese) |
| [8] | 蒋雪辉. 气液2相离心甩油盘设计[J]. 航空发动机, 2016, 42(2): 7-10. JIANG X H. Design of a gas-liquid two-phase slinger[J]. Aeroengine, 2016, 42(2): 7-10. (in Chinese) |
| [9] | WARDLE K E, ALLEN T R, ANDERSON M H, et al. Free surface flow in the mixing zone of an annular centrifugal contactor[J]. AIChE Journal, 2010, 54(1): 74-85. |
| [10] | 张健, 方杰, 范波芹. VOF方法理论与应用综述[J]. 水利水电科技进展, 2005, 25(2): 67-70. ZHANG J, FANG J, FAN B Q. Advances in research of VOF method[J]. Advances in Science & Technology of Water Resources, 2005, 25(2): 67-70. DOI:10.3880/j.issn.1006-7647.2005.02.019 (in Chinese) |
| [11] | 黄勇, 林宇震, 樊未军, 等. 燃烧与燃烧室[M]. 北京: 北京航空航天大学出版社, 2009. HUANG Y, LIN Y Z, FAN W J, et al. Combustion and combustor[M]. Beijing: Beihang University Press, 2009. (in Chinese) |
| [12] | 王文波, 岳明, 杜春雨, 等. 喷嘴损伤对环形回流燃烧室性能的影响[J]. 航空动力学报, 2017, 32(6): 1335-1341. WANG W B, YUE M, DU C Y, et al. Effect of nozzle damage on performance of annular backflow combustor[J]. Journal of Aerospace Power, 2017, 32(6): 1335-1341. (in Chinese) |
| [13] | SHIH T H, LIOU W W, SHABBIR A, et al. A new k-epsilon eddy viscosity model for high Reynolds number turbulent flows:Model development and validation[J]. Computers & Fluids, 1995, 24(3): 227-238. |
| [14] | YAKHOT V, ORSZAG S A, THANGAM S, et al. Development of turbulence models for shear flows by a double expansion technique[J]. Physics of Fluids A, 1992, 4(7): 1510-1520. DOI:10.1063/1.858424 |
| [15] | HUPPERT H E. Flow and instability of a viscous current down a slope[J]. Nature, 1982, 300(5891): 427-429. DOI:10.1038/300427a0 |
| [16] | HARTLEY D E, MURGATROYD W. Criteria for the break-up of thin liquid layers flowing isothermally over solid surfaces[J]. International Journal of Heat and Mass Transfer, 1964, 7(9): 1003-1015. DOI:10.1016/0017-9310(64)90042-0 |
| [17] | RUCKENSTEIN E. On the break-up of thin liquid layers flowing along a surface[J]. International Journal of Heat and Mass Transfer, 1971, 14(1): 165-169. DOI:10.1016/0017-9310(71)90150-5 |
