半柔壁喷管气动设计关键控制参数研究

采用数值模拟的手段,对半柔壁喷管流场进行模拟.主要目的是通过研究流场的均匀性,确定喷管各气动设计控制参数对设计结果的影响,为半柔壁喷管的设计优化提供指导.本文着重研究了不同喷管最大膨胀角θA,喷管半消波区控制参数θB,喷管半消波区分布控制参数m,边界层修正角,收缩段型面形式等因素对流场的影响.对得到的结果进行了分析,初步确定了进行半柔壁喷管设计的控制参数,为下一步工作打下了基础.     

固体火箭发动机碳基材料喷管热化学烧蚀特性

为了准确预示固体火箭发动机碳基材料喷管的烧蚀率,依据热化学烧蚀理论,建立了喷管传热烧蚀的二维轴对称气-固-热耦合计算模型,计算通过FLUENT壁面化学反应模型完成,无需事先假设烧蚀控制机制。针对70-lb BATES发动机喷管进行了烧蚀计算,研究了推进剂配方、氧化性组分、燃烧室压强对喷管烧蚀的影响。结果表明:烧蚀率计算值与试验测试值吻合较好;烧蚀率分布遵循喷管内壁热流密度分布规律,在喉部上游入口处达到峰值;烧蚀率随推进剂Al含量增加而降低,随燃烧室压强升高而近似正比例增大;H2O是决定烧蚀的主要氧化性组分。     

固体火箭发动机喷管喉部凝相颗粒粒度分布实验

设计了一种新的收集固体火箭发动机喷管凝相颗粒的实验装置,针对典型的HTPB复合推进剂,开展了喷管喉部凝相颗粒的收集实验和粒度分析,研究了燃烧室压强和收敛角度对喷管喉部颗粒粒度分布的影响规律。研究结果表明,喷管喉部的凝相颗粒在0.27～50μm之间都有颗粒存在,凝相颗粒主要集中在0.3～15μm之间,粒径大于15μm的颗粒较少;燃烧室压强对颗粒粒径有较大影响,随着燃烧室压强的升高,凝相颗粒粒径变小,粒度分布更为集中;燃烧室压强相同的条件下,收敛角度对喷管喉部的凝相颗粒粒度分布影响较小。     

多喷管射流气动声学特性的数值研究

采用大涡模拟结合动力Smagorinsky亚格子模型对四喷管超声速高温射流进行了三维非稳态数值模拟。基于非稳态计算的结果,利用FW-H面积分方法对四喷管射流远场的声学特性进行计算,并将结果与单喷管得到的结果进行对比。结果表明,多喷管射流形成了复杂的流场,其产生的气动声场并不仅仅是单喷管简单的叠加,由于四股射流之间的相互影响使得声场也产生了新的特点,不但表现在总声压级数值上,还对改变了射流噪声辐射的指向性,最大声压级的位置由原来的50°左右变化到30°。     

合排气二元收敛喷管与发动机流量匹配设计方法

发展了一种基于喷管三维黏性数值分析的内外涵气流质量流量匹配设计方法.该方法主要是通过建立以喷管待定几何参数为未知量的非线性无量纲内外涵气流质量流量匹配方程组,然后采用逐步线性迭代法结合喷管三维黏性数值分析对该非线性方程组进行求解,得到满足发动机内外涵流量匹配的喷管设计方案.通过数值算例验证分析,结果表明：①该方法可以快...     

摆动喷管结构分析的半解析法

针对摆动喷管结构轴对称的特点，采用半解析法，选用精度较高的八结点等参数环元，对缠绕复合材料旋转体在非轴对称载荷作用下建立了有限元方程，在微机上采用ＦＯＲＴＲＡＮ语言编写了计算程序。通过算例，验证了此种方法的正确性和可行性，应用本文所提供的计算方法和程序，对某型固体火箭发动机缠绕复合材料摆动喷管进行了静力分析，从而给出了喷管静力分析的一种工程计算方法。     

喷口辐射热泄漏对延伸段壁温分布的影响

引入喷口的热辐射泄漏损失，导出辐射冷却推力室壁温的计算方程。分析表明，考虑与不考虑喷口热损失二者相比较，前者壁温约下降７０℃。本结果可为喷管延伸段材料的选取提供更为科学的理论依据。     

真空小喷管羽流场的Monte Carlo直接模拟

利用直接模拟的Ｍｏｎｔｅ Ｃａｒｌｏ方法数值模拟了真空喷管羽流场，获得了符合规律的数值模拟结果。由于将随机取样频率法用于轴对称问题中，解决了计算机内存和速度上的困难，同时也证明了Ｒ－Ｓ－Ｆ方法的可靠性。     

两种椭圆型方程求源项方法在喷管内流场网格生成中的应用

用椭圆型方程生成二维网格，发展两种Hilgenstock求源项方法，实现对网格与边界间距及与边界正交的控制。     

椭圆截面喷管排气辐射热流数值分析

文中对直升机红外抑制器椭圆截面喷管排气的辐射热流进行数值分析,并预估排气温度分布、速度和燃气成分分布对辐射热流的影响。为了计算辐射热流和考虑尾焰对周围环境的传热,采用形式简单并可用通用方程计算程序计算的六热流辐射模型。由于排气流场内温度和密度变化较大,故紊流模型采用考虑密度变化的修正双方程k-ε模型。本文中所用的比热C_(pm)是随气流温度和成分而变化,温度用迭代法求得。喷管出口截面为椭圆形,排气为三元自由射流,控制方程为三元抛物形偏微分方程,可采用前进积分求解,这样既可满足计算精度要求,又可节省机时。并根据三元喷流周向边界条件周期性变化的特点,采用CTDMA解法。计算结果与试验结果基本一致,说明本计算方法是可行的。