导弹发射管内燃气流场的实验分析

针对密封盖未及时开启造成发射装置的损坏问题进行了实验分析。实验测量了对发射管壁静压和密封盖受到的燃气射流冲击压力，利用子波分析方法求取了实验数据的基于子波系数的频带能量谱，并对含有射流噪声的测量数据进行了子波滤波，得到了各点压力低频成分的变化情况，同时得到了密封盖所受冲击冲量与激波振荡幅值的关系。研究表明密封盖不及时开启会在弹管间隙中形成高幅值的起始冲击波和燃气回流激波振荡。这种破坏性的流动结构会导致发射装置的损坏。     

液体火箭发动机推力室粘性流场数值模拟和实验验证

运用耦合点隐式方法的MacComak两步差分格式结合κ-ε湍流模型求解肼类燃料发动机推力室中的粘性化学反应流动。化学反应采用17种组分，12个主要反应的有限速率的化学反应模型，得到了流动参数在推力室中的分布情况。其结果同理论分析的相一致，最后用实验验证了用数值模拟方法的正确性，研究结果对发动机的设计、性能和理论研究具有重要的实际意义和参考价值。     

武器分离及舱门开启过程数值模拟研究

基于非结构重叠网格技术,发展了一套能较好模拟空腔流动和武器分离问题的数值计算方法.通过模拟典型空腔流动和WPFS模型分离轨迹,并与试验结果比较,说明本文方法可较好模拟空腔复杂流动和武器分离轨迹.采用发展的数值方法对某内埋弹舱舱门开启过程和内埋武器分离轨迹进行了数值模拟,研究了舱门运动的非定常效应对舱门气动特性的影响.研究表明:一定条件下,舱门运动产生的非定常效应对舱门气动特性影响较大,舱门和舱门控制系统设计时必须考虑；在本文给定条件下,弹舱内的内、外侧导弹均能安全分离.通过研究,为内埋弹舱的舱门设计和内埋武器分离安全评估提供了技术手段和依据.     

GRP格式推广用于多维复杂流动的数值计算

采用算子分裂技术，将原始求解一线流动问题的高精度、高分辨率GRP格式推广用于包含各种不同间断解及其相互作用的多维复杂流动计算。以轴对称管端出口处激波绕射传播为例，作了数值模拟，得到与实验结果吻合很好的流场精细结构。     

微波等离子推力器喷管流动数值模拟

为了准确把握微波等离子推力器喷管流动的机理与特点,采用冻结流和非平衡流两种模型对其进行了对比数值模拟分析。非平衡流模型考虑了流动过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分、4反应的有限速率化学反应模型,采用二阶精度NND格式数值求解耦合化学反应源项的N S方程组。数值模拟的结果揭示了喷管内的流场结构,反映了喷管内的离解电离状况,得到了推力器的推力和比冲。分析表明,数值方法有效,计算结果合理,具有工程应用价值,能为相关研究提供参考。     

驻点壁面催化速率常数确定的研究

以平衡流动作为热环境估算的依据，提出了用数值求解非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和实验测量热流值确定模型表面材料催化速率常数的方法。用５组分１７个化学反应Ｄｕｎｎ－Ｋａｎｇ空气化学模型和轴对称热化学非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对激波管中球头和平头圆柱模型绕流流场进行了数值模拟，给出了驻点热流随催化速率常数变化的分布，并根据激波管实验测量的热流值确定了表面材料Ｐｔ、ＳｉＯ２、Ｎｉ和某种     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

GPR格式推广用于多维复杂流动的数值计算

采用算子分裂技术，将原始求解一维流动问题的高精度、高分辨率ＧＲＰ格式推广用于包含各种不同间断解及其相互作用的多维复杂流动计算。以轴对称管端出口处激波绕射传播为例，作了数值模拟，得到了与实验结果吻合很好的流场精细结构。     

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

叶根间隙泄漏流对跨声压气机转子性能的影响

为研究压气机转子上游叶根轴向间隙泄漏流对转子性能的影响程度以及泄漏流与转子流场之间的相互作用机制,对带有叶根泄漏流的跨声压气机转子流场进行了数值模拟。计算时分别采用了两种模化途径,即耦合求解转子上游叶根间隙内空腔与转子通道内的流场,以及求解将叶根间隙模化为进口边界的转子流场。结果表明:叶根泄漏流可使叶根附近流动堵塞增大,造成转子气动性能的恶化。另外,叶根泄漏流的周向速度分量越大,泄漏流对转子流场的不利影响也越大。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

轴流涡轮基元级静叶稠度特性的数值研究

针对不同静叶稠度(静叶叶片数)的轴流涡轮基元级进行了非定常数值模拟,研究了静叶稠度对涡轮基元级流动状态和损失情况的影响.结果表明静叶稠度的改变对涡轮基元级流动状态和损失的影响直接与动叶稠度相关.静叶稠度的改变影响通过涡轮基元级的流量,在动叶稠度不变时,会引起涡轮基元级反力度的改变.静叶稠度增加到一定程度时,会使气流在静子中的膨胀加速过于剧烈而产生激波损失及激波与边界层干涉带来的边界层分离损失.静叶稠度减少到一定程度时,会使转子中的流动状态极大恶化,进口极大的负攻角致使动叶压力面发生大范围的分离.存在一个最佳的静叶稠度,使涡轮基元级的损失最小.      

RP-3航空煤油模拟替代燃料的化学反应详细机理

在化学激波管中对RP-3航空煤油的着火特性进行了实验测量,获得了多工况下RP-3航空煤油的着火延迟时间.根据RP-3航空煤油的化学组成及物理特性,提出了由体积分数分别为0.65,0.1,0.25的正癸烷、甲苯与丙基环己烷3种组分组成的模拟替代燃料,并形成了该模拟替代燃料的化学反应详细机理.采用该化学反应详细机理对该模拟替代燃料在化学激波管中多工况下的着火特性进行了数值计算,并与实验数据进行了对比分析.结果表明:在不同压力与当量比下,RP-3航空煤油着火延迟时间的对数与着火温度的倒数呈线性关系,并且随着火温度与压力的升高以及当量比的降低,着火延迟时间逐渐缩短;同时,在各工况下采用该化学反应机理计算得到的该模拟替代燃料着火延迟时间与RP-3航空煤油着火延迟时间的实验值吻合良好.      

不同飞行状态下固体火箭发动机尾喷焰数值研究

为了研究飞行状态对固体火箭发动机尾喷焰的影响,建立了含化学反应项和组分输运项的N-S(Navier-Stokes)控制方程,并采用MUSCL(monotonic upstream-centered scheme for conversation laws) Roe格式的有限体积法进行求解.在此基础上,结合热力计算结果对某固体火箭发动机在不同飞行状态下的化学反应流和单一组分流进行仿真计算.结果表明:初始倾角和第一马赫波节长度随飞行高度的增加而增大,随来流马赫数的升高而减小.同单一组分流相比,化学反应流对第一马赫波节的影响较初始倾角显著.射流轴线温度因飞行状态的不同有较大涨落差异.初始倾角、第一马赫波节由来流总压决定,射流轴线温度由激波结构、复燃效应、来流条件等因素共同决定.      

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

燃料射流隔板对受限超声速反应混合层流场和燃烧特性影响

为了研究燃料射流隔板壁面扰动对于受限超声速反应混合层流场和燃烧特性的影响,基于抽象物理模型激波和膨胀波理论推演,得到了隔板诱发壁面扰动下的受限超声速反应混合层流场结构,并采用数值模拟方法进行了验证。数值结果表明,隔板诱发壁面扰动下的受限超声速反应混合层流场主要由冷态主导,且相对无隔板扰动多了回流区、激波、膨胀波、波与反应混合层相互作用等复杂现象。在此基础上,对于不同隔板厚度计算的结果表明,随着隔板厚度的增加,隔板下游的回流区增大,上下膨胀产生的压力不平衡加剧,反应混合层会产生偏斜,回流区厚度和偏斜距离与隔板厚度成正比。此外,第一道反射激波存在降低点火延迟的作用,且存在一个隔板厚度阈值,阈值以下随着厚度增加点火延迟随之降低。同时,后续的多道反射激波导致混合层发展的局部起伏和局部燃烧增强。     

共用尾喷管多管脉冲爆震发动机数值模拟研究

采用基元反应模型和显式迎风TVD差分格式数值模拟方法,研究共用尾喷管多管脉冲爆震发动机(Pulse Detonation Engine,简称PDE)各爆震管之间相互影响.采用点隐算法解决化学反应引起的刚性问题,模拟了共用尾喷管内的流场,以及爆震波退化为激波遇到尾喷管收敛斜面反射及反射波向爆震管上游传播的过程.研究结果为多管PDE的共用尾喷管设计提供了理论基础.