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

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

非平衡流计算方法及其模拟激波诱导振荡燃烧

给出一种数值模拟高超声速化学非平衡流动的计算方法，采用流动项与化学反应生成源项解耦处理，即可以有效地解决非平衡流动计算中遇到的刚性问题，还具有公式推导简单和应用方便的优点。从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超音速H2／Air预混合气体钝体绕流场，较好分辨出流场中激波诱导振荡燃烧非定常流动现象，表明计算方法具有较高的精度。对计算方法的效率和网格数的影响也进行研究。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

非结构网格热化学非平衡Navier-Stokes方程组求解技术研究

采用非结构混合网格数值求解三维热化学非平衡Navier-Stokes方程组获得高超声速流场各组分气动热力学参数,并计算物面热流值。利用有限速率化学反应模型和二温度气体模型分别描述化学非平衡和热力学非平衡过程,通过Park振动——离解模型刻画化学非平衡和热力学非平衡耦合效应。无黏通量项采用AUSM格式离散;对无黏通量项和源项采用隐式时间推进格式,增强数值求解过程稳定性。数值算例表明,为得到正确的物面热流值,非结构混合网格边界层y+需减小至0.5,且网格增长率不超过1.4;采用元素摩尔数守恒关系式计算化学反应源项,可有效降低数值残差积累对计算小组分气体的不利影响,提高各气体组分计算准确度。     

三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

塞式喷管流场数值模拟

采用ＮＮＤ２Ｍ差分格式,从考虑к－ε湍流模型的理想气体三维薄层近似Ｎ－Ｓ方程出发,对塞式喷管内外流场进行了数值模拟。计算结果与国外文献十分一致。通过流场分析,对塞式喷管流动机理有了清楚的认识。计算表明ＮＮＤ２Ｍ差分格式及相应计算软件适合于模拟复杂喷流流场,可用于塞式喷管发动机的研究。     

COIL化学非平衡流动的一种解耦方法研究

提出了一种计算氧碘化学激光器化学非平衡流动的解耦方法.求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分混合反应流场.空间离散格式为AUSM up格式,用四步龙格-库塔方法作显式时间推进.假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型.使用一种"点隐式"方法来处理化学源项的刚性问题.计算结果表明本文方法可以有效解决计算中的刚性问题,与现有隐式算法相比减少了大量的矩阵运算,公式推导简单,易于实现和改进.     

C-H-O-N体系液体火箭发动机喷管流场数值模拟

采用弱耦合点隐式方法的MacCormack格式对C－H－O－N体系的液体火箭发动机喷管粘性流场进行了数值模拟，采用了考虑化学反应和不考虑化学反应的两种模型，其中化学动力学模型为C－H－O－N体系的12组分、14反应有限速率的化学反应模型。数值模拟得到了流场参数在喷管中的分布，结果表明，数值模拟结果与理论分析一致，为发动机的设计及试验提供了参考。     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

塞式喷管流场计算

采用NND差分格式，对塞式喷管轴对称流场进行了数值模拟。通过流场分析，可以对塞式喷管火箭发动机有清楚的认识。计算表明，NND格式有很好的分辨率；无粘和有粘的情况流场结构和性能相差悬珠；塞式喷管的自适应能力很好。     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

喷管型面对固体火箭发动机性能的影响

用时间相关法数值计算跨音速喷管流动时,其控制方程分别采用守恒型和非守恒型方程.通过对JPL喷管进行计算发现:在相同的精度下,非守恒型方程的收敛速度比守恒型方程的收敛速度快一倍多.采用修正阻尼法(DASHCO法)进行数值计算,其收敛速度均提高六倍以上.最后,采用直接优化法研究了喷管收敛段型面及喉部型面对发动机性能的影响.     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

Research on vacuum plume and its effects

In vacuum environment, the exhaust flow of attitude control thrusters would expand freely and produce the plume, which possibly causes undesirable contamination, aerodynamic force and heating effects to the spacecraft. Plume work station (PWS) is developed by Beihang University (BUAA) for numerically simulating the vacuum plume and its effects. An approach which combines the direct simulation Monte Carlo (DSMC) method and difference solution of Navier-Stokes (N-S) equations is applied. The internal flows in nozzles are simulated by solving the NS equations. The flow parameters at nozzle exit are used as the inlet boundary condition for the DSMC calculation. Experimental studies are carried out in a supersonic low density wind tunnel which could simulate the 60-80 km altitude environment to investigate the plume and its effects. To demonstrate the capability of PWS, numerical simulations are performed for the vacuum plume of several typical attitude control thrusters. The research results are of great help for the engineering design.     

固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

计算非对称突然膨胀槽道流动的加速多重扫描耦合推进方法

本文提出加速多重扫描耦合推进方法,并用于非对称突然膨胀槽道流动计算。数值结果与实验数据吻合得很好,从而表明这个方法是有效的,精确的。依据对数值结果的分析,作者给出能够描写存在大范围分离的不可压缩简化Navier-Stokes方程组的最优形式,并且对这组方程的椭圆型数学性质作了进一步的分析。作者还对扫描过程的收敛性,推进过程的稳定性进行了理论分析,此外还讨论了不同精确度的差分格式对解的影响。