论文部分内容阅读
诞生于21世纪的涡流冷壁推力室概念(Vortex Combustion Cold-Wall Chamber,VCCWC)是涡流技术工程化应用的典型代表,该推力技术的核心内容为限制在推力室范围内、具有特殊速度分布形式的涡流,该涡流存在轴向速度反向,在推力室壁面与燃烧反应区之间形成外涡冷流,该外涡冷流对推力室壁面起到了隔离热量传递、吸收热量传递的双重作用。本文采用数值计算与理论推导相结合的研究方法,对VCCWC内的核心流动:双向涡流进行系统深入的研究,以理论模型结果与实验结果为参考,揭示受限空间内双向涡流流动及其雷诺应力特征,进而提出其外涡冷流所具有的隔离、吸收热量传递的双重作用。以实验数据和解析模型作为参考,以RSM模型计算结果作为对照,研究了基于各向同性涡粘性假设湍流模型中的旋转修正方法在双向涡流数值计算中的适用性。三种旋转修正方法中,附加旋转修正方法性能最差,附加旋转修正项对双向涡流几乎没有作用。漩涡修正方法在修正因子取默认值0.07时也没有太大改善。曲率修正方法在切向速度预测上相较于前两种方法有了非常显著的改善,但是在轴向速度分布预测上,曲率修正方法依然无法反映出双向涡流最本质特征:壁面附近的轴向速度反向。通过三种旋转修正方法与实验及理论解的比较可知,只有切向速度分布与轴向速度分布同时符合双向涡流理论,才能确认双向涡流数值计算结果的正确性。从Boussinesq各向同性假设出发,以实验数据和理论模型作为参考,研究了双向涡流雷诺应力方向性特征。根据计算结果可知,二次压力-应变RSM模型在预测双向涡流方面存在一定优势,其准确程度较线性RSM模型有了一定提高。双向涡流场雷诺应力随涡雷诺数增大而增大,雷诺应力方向性均显著,在径向雷诺应力峰值的无量纲半径位置处,径向雷诺应力与切向雷诺应力差距达到了切向雷诺应力的150%左右。RNG模型所得雷诺应力与RSM模型存在较大差别,各向同性假设使得RNG模型所得雷诺应力趋于均匀。雷诺应力预测错误最严重的是切向雷诺应力,其分布趋势与两种RSM模型相去甚远。径向雷诺应力与轴向雷诺应力虽然给出了正确的雷诺应力分布趋势,但是数值大小上与准确值相去甚远。未来对涡粘性湍流模型进行考虑双向涡流特征的修正时,需要人为引入雷诺应力方向性特征。以实际涡流冷壁推力室中的收敛喷管为研究对象,数值模拟研究了收敛段入口尺寸所形成的肩径比对推力室内压力场及速度场的影响。肩径比远大于无量纲涡幔半径时,收敛段入口处与肩部涡之间形成的出口涡破坏了流动的顺畅性,大的肩径比通常无法达到推力室喷管的收缩比要求。肩径比与无量纲涡幔半径相当时,形成与理论模型非常接近的理想双向涡流,此时推力室效率达到最高。当肩径比小于无量纲涡幔半径时,收敛段内会形成虚拟收敛段,该虚拟收敛段不依赖于肩径比,入射口氧化剂来流中的一部分会直接从喷管流出而并不参加内涡燃烧反应,因而推力室效率显著降低。旋流条件与双向涡流特性耦合研究中,通过数值计算发现,切向速度及最大切向速度均随入射倾角的增大而减小,随入射压降增大而增大,随长径比的增大而减小。不同入射倾角、不同入射压降、不同长径比下整个发动机内最大切向速度的无量纲径向位置均恒定在0.19附近。对于不同长径比工况,长径比为1.0时,最大切向速度从发动机顶端到入射口附近逐渐增大,长径比为1.5时,最大切向速度从发动机顶端到入射口附近先增大后减小。切向速度及最大切向速度的轴向衰减率维持在3%以内。无量纲涡幔半径从发动机顶端开始线性增大到入射口附近,变化范围为0.71至0.82。入射倾角相同时,随入射压降的增大燃烧室长径比对最大切向速度大小的影响将随之增大通过理论推导方法,对涡流冷壁推力室内的传热过程进行了理论分析与建模过程,建立了热流量平衡关系式。通过理论分析指出,外涡氧化剂冷流存在两个作用:外涡氧化剂冷流的存在,导致内涡向壁面传热的热阻增大了 14倍以上;由于所吸收的一部分热量用于自身温度的升高,外涡氧化剂冷流对壁面事实上起到了冷却作用,导致外涡向壁面的传热量只占内涡总传热量的60%以下,因此可知外涡氧化剂冷流起到了阻碍热量传播与冷却推力室壁面的双重作用。