壁面边界对撞击合成动量角的影响研究

为了研究动量比对针栓式喷注器撞击雾化合成动量角的影响规律,基于三相的与水平集耦合流体容积法(CLSVOF)对贴壁液膜/自由液膜撞击的合成动量角进行了数值模拟,并与自由液膜/自由液膜撞击的合成动量角进行了对比,结合试验结果及理论预估结果详细考察了有无壁面条件下合成动量角与动量比的关系,深入分析了造成两者之间差异的根源,揭示了壁面边界对合成动量角产生影响的作用机制。结果表明:CLSVOF方法计算的合成动量角与试验结果一致,最大相对误差约为10%,大多数工况点的相对误差小于5%;有壁面边界的贴壁液膜/自由液膜撞击合成动量角显著大于无壁面边界的自由液膜/自由液膜撞击的;仅一路流体贴壁的撞击合成动量角与常用的入口动量比理论预测值最大相差20°以上,而两路均无贴壁或者均贴壁的撞击合成动量角与理论预测值吻合很好。这一显著差异的根源在于两者撞击形成的高压区分布显著不同。壁面的存在使得撞击点附近形成的高压区对两路流体作用不对称,有壁面边界时壁面承受高压迫使其对贴壁流体有强的作用力,导致垂直于壁面方向动量不守恒,根据入口动量比预测的理论不再适用。

针栓式喷注器液膜下漏率预估模型

为了对针栓式喷注器液膜下漏率进行准确预测,基于针栓式喷注单元喷雾场结构分析,结合理论推导、数值仿真及试验研究3种方法建立了液膜液束各自变形的相对变形量模型;在考虑液膜液束变形的基础上,引入相互影响系数表征多喷注单元间相互影响,建立了实际阻塞率和实际下漏率模型。通过数值仿真及试验结果的多参数充分验证,结果表明:理论预估模型与数值仿真及试验结果一致性较好。液膜液束相互作用下,液膜绕液束流动和液束根部横截面前后缘位置移动不同步导致的展向变宽分别是液膜和液束发生变形的主因,且有效动量比越大,液膜相对变形越大,液束相对变形越小。对于一定阻塞率的几何结构,结果表明:下漏率随着有效动量比的增大而增大,增大趋势呈先快后缓,且实际下漏率均小于几何下漏率,这是由膜束变形导致的实际阻塞率比几何阻塞率更大造成的。另外,发现液膜下漏率仅与表征流场结构(有效动量比)及几何结构的无量纲参数(液膜厚度与液束直径之比和阻塞率)有关,与喷射速度的绝对值无关,并给出了模型中的常系数供工程设计预估使用,对从设计初期就考虑针栓头的热防护问题具有重要的指导意义。

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构,基于自适应网格加密技术和分三相计算的PLIC VOF(Piecewise Linear Interface Calculation Volume of Fluid)方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析,通过对两路推进剂分别进行界面追踪,获得了膜束撞击雾化混合过程的详细结构特征,与高速摄影试验结果定性定量对比均吻合较好,验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析,结果表明:膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形,膜束同时发生弯曲变形和横截面变形。另外,膜束撞击同时受到正压和剪切应力作用,导致了一系列复杂涡流现象,使得相互作用增强,雾化混合均增强,这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后,还发现膜束撞击喷雾场液滴分布呈现分区结构特征,分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区,实际中应兼顾雾化特性和混合特性,选取中等动量比膜束撞击,这可为针栓式喷注器的理论研究和工程设计提供重要参考。