当前箭弹正向超远程、高精度的方向发展,但火箭推力偏心却不同程度地影响着火箭弹的射程和精度。为增加火箭弹射程和精度,国外在研究火箭喷管推力偏心特性、模拟喷管内流场方面进行了大量的研究工作,包括三文流[4,5,6]、多相流、非平衡化学反应流以及燃气舵矢量喷管流等。对于喷管流场的模拟方法已经采用了Euler方程和雷诺平均N-S方程。这些研究的目标主要是提高燃气能量的利用率,侧重于提高主推力和推力系数,使其与喷管设计的目标相吻合。9990
对于在非对称扰动因素作用下形成的非对称流场,国内外相关的研究报道并不多见。国外在80年代中期以前主要用一文流附加小扰动来处理喷管内非对称流场;用三文特征线法模拟喷管扩张段的三文流场;用点泉流加小扰动在球面坐标系下求解锥形喷管内非对称流场等;用TVD法模拟喷管扩张段的三文流场。80年代后期以来没有发现国外在喷管非对称流场数值模拟方面的相关报道,也未发现将数值模拟方法与微推偏喷管设计理论相联系的研究报道。
关于喷管流场数值模拟格式,751七十年代主要用MacCormack格式和Beam-Warming格式求解Euler方程和N-S方程。用于喷管内流场数值模拟的格式还有Lax-Werdroff格式、LuSSor方法、Newton迭代法、TVD格式[7]和特征线法[4,5,8,9,10]。自从80年代起,A.Harten提出TVD格式,TVD格式被广泛采用。而特征线法作为一种成熟的数值模拟格式,在求解双曲型偏微分方程时,它具有透彻地了解流动特性的优点,也具有精度高、收敛速度快、程序编程较简便的优点。
喷管非对称内流场的试验也已做过了不少[1,3,7]。用冷气替代燃气测量喷管壁面上的分布压力,然后积分分布压力求侧向力和侧向力矩的大小;地面发动机静止试验采用57-2航箭发动机在751分力试验台上进行静止试验;火箭模型地面飞行试验利用已测试过侧向力和力矩的发动机A和B。B用微推偏喷管,A用普通喷管,对准立靶,发射A和B,B的概率误差比A小得多。该研究机构已具备测试喷管非对称内流场产生的侧向推力和主推力的条件。其测试设备主要采用751分力试验台、示波器及微机采样系统组成的试验装置,这些可以为研究不同类型的火箭发动机推力偏心特性提供保障。
在喷管内流场数值模拟领域国内己开展过大量的研究工作,而有关非对称流场数值模拟研究主要集中在南京理工大学火箭实验室里。90年代以前用一文小扰动法研究喷管非对称流场。由于计算方法的限制,计算精度只适用于喷管扩张半角较小的锥形喷管;90年代中期以来相继发展了高分辨率的TVD格式[7,11]和高精度双特征线法[9,10,12,13]模拟喷管内三文非对称流场。目前,已经研究了几何偏心、总温总压分布不均匀等扰动因素引起的推力偏心特性。这些研究主要是以特征线法和三文有限体积TVD格式模拟火箭喷管内三文定常非对称流场。在诸如喷管优化设计等方面需要反复多次计算喷管内流场的过程中,特征线法显现出独特的优越性,而且特征线数值计算法的程序通用性较好。
南京理工大学张庆民、武晓松[14] 等人应用巴特勒提出的五面体双特征曲线网格方法对火箭喷管扩张段内的非定常非对称流场进行了数值模拟,获得了锥形喷管侧向推力随时间的变化曲线。分析表明,该文计算的侧向推力结果与定常状态下的数据是吻合。证明五面体双特征曲线网格方法可用来研究非定常流动条件下的推力偏心特性。
武晓松、王政时[10] 等人研究了当气流偏离抛物型喷管超音速段对称轴时所引起的喷管推力偏心,用三文特征线方法数值模拟了喷管内的非对称流场,得到了在不同入口扩张半角时,喷管侧向推力分量沿轴向的变化及其零点位置在轴线上的分布。结果表明,侧向力第一零点和第二零点的位置随入口扩张半角的变化都是非线性的。 火箭喷管推力偏心特性国内外研究现状:http://www.751com.cn/yanjiu/lunwen_8886.html