环流喷管的轴对称湍流流场计算

对火箭发动机环流喷管进行了二维数值仿真,用数值方法求解Navier－Stokes方程,数值格式为显格式时间推进格式,湍流模型选择可靠性较好的两层代数湍流模型。由于所研究的流场边界形状复杂,计算中采用了分区的方法。通过多个算例的对比,确定了使环形喷管推力效率较高的几何形状。画出了环流喷管内部的计算网格图、速度矢量图、等马赫线图、等压线图及喷管上壁面及轴线的压力分布图。     

四机并联发动机低频动态特性分析

为了获得发动机结构的低频动态特性,采用有限元方法进行了液体火箭一级四机并联发动机(含机架)的计算模态分析。由于发动机为复杂的非对称结构,建模时按质量和刚度等效原则将推力室的夹层结构简化成单层板,其余结构大都按空间梁简化,通过改变材料密度保证了模型质量与发动机稳态工作时的质量一致。简化后的计算规模适中,既反映了结构的主要特征又不浪费计算机资源和时间。计算结果表明了四机并联发动机系统中有多阶轴向运动较大的模态值得注意,需避免它们和发动机的低频脉动及箭体的低频振动耦合。     

固体火箭喷管颗粒尺寸分级两相跨音速流场计算

本文中对固体火箭喷管颗粒尺寸分级的两相跨音速流场作了计算.气相控制方程采用隐式近似因子分解法求解,尺寸分级的颗粒控制方程采用特征线求解,然后,二者进行充分的耦合,可以获得固体火箭发动机含有任意颗粒质量分数和不同颗粒尺寸时轴对称喷管跨音速流场的参数分布.文中讨论了不同颗粒半径和质量分数对流场的影响,对单一颗粒尺寸和颗粒尺寸分级的参数进行了比较.两相耦合计算的迭代收敛速度取决于气相,本文中气相方程求解的格式除部分边界外是隐式的,CFL数可取至6左右,收敛速度快.特别是对颗粒尺寸分级的计算,得益更大,其得益的倍数为颗粒的分级数.     

微波与等离子体之间的相互作用

根据微波和等离子的基本理论以及它们之间相互作用的原理,建立了由微波能、谐振腔内电场、工质流场、等离子体平衡场等组成的轴对称偏微分方程组。通过对模型的计算,获得了放电管内等离子体区的形状、电场、电子浓度和温度场,以及气体温度场。     

涡轮泵转子的临界转速研究(Ⅰ)均匀支承转子临界转速的传递矩阵法

根据转子振动时的受力和位移状况,分别导出轴、盘和支承的传递矩阵,再根据转子的结构布局导出频率和振幅等的计算方程,可求得转子的各阶临界转速及振幅。传递矩阵法对转子的不同布局方案运用十分方便灵活。还给出了减振阻尼、支承弹性的组合刚度计算方法。计算表明,支承的弹性、阻尼和质量对转子的振动特性（临界转速和振幅）有重大影响,转子与支承的匹配对转子的振动特性能起调节控制作用。     

微波等离子推力器同轴谐振腔内的电磁场特性

为了研究同轴谐振腔内导体的头部形状、位置对腔内电磁场分布的影响规律,用时域有限差分法(FDID)对微波等离子推力器(MPT)TEM谐振模式的同轴谐振腔中的电磁场进行了数值模拟,并结合实验现象对数值模拟的结果进行了分析讨论。结果表明,内导体头部形状的改变,改变了同轴腔内电磁场的分布规律和电场强度最大集中区域的位置。内导体位置的改变,只改变电磁场强度的大小,且内导体愈接近同轴腔的底面(即：放电间隙δ愈小),腔内电场强度愈强。同时,锥形头部的内导体,更有利于MPT的启动和稳定工作,为采用同轴谐振腔的MPT工程样机的研制提供了依据。     

小型微波等离子推力器研究

 介绍约100W 功率微波等离子推力器的系统结构, 谐振模式, 电磁场分布, 电磁场、等离子体和流场的耦合, 喷管中的等离子体复合流动, 真空羽流, 原理实验样机性能及应用前景。     

100W微波等离子推力器的调试与地面实验

微波等离子推力器(MPT)谐振腔只有在谐振状态下，微波能量才能被高效地用于加热工质产生推力。采用微波无源器件回波损耗的测试方法对100W MPT谐振进行精确地调谐，分析微波能量的吸收效率及谐振频率带宽，研究腔体尺寸、微波耦合探针位置在同一谐振频率条件下的匹配性。在地面条件下，对100W MPT系统进行了实验，利用推力器腔体内的高电场强度可直接使系统推力器腔体中工质电离，其最高能量吸收效率可达98．3％。     

预燃室三维湍流和燃烧过程的数值模拟(Ⅱ)数值模拟结果及分析

用三维湍流N－S方程和单步快速不可逆化学反应描述液氧－煤油液体火箭发动机预燃室内的三维湍流和燃烧过程。采用同位网格和SIMPLE算法求解控制方程，得到了喷注单元和预燃室内参数的详细分布。结果表明，预燃室结构设计合理，其出口处燃气浓度、温度分布均匀，质量加权平均温度与实际温度拉近，同时表明，预燃室头部的喷注单元和液氧二次喷注孔的结构排列，喷注单元的流动和燃烧状况，液氧二次喷注孔的入口参数等，对预燃室出口燃气温度等参数分布的均匀性影响很大。     

涡轮泵转子的临界转速研究(Ⅱ)──非均匀支承转子临界转速的传递矩阵法

当支承的组合刚性沿圆周不均匀时,可将其刚性分解到两个互相垂直的主平面中,在主平面中具有最大和最小刚性。同时,把转子振动时的受力和位移分解到主平面中,用八阶矩阵分别建立轴、支承、盘的传递矩阵。然后根据具体的转子布局方案可导出临界转速和振幅计算方程。研究表明,非均匀支承转子与均匀支承转子相比,频率谱有质量差别,前者可能出现的临界转速数增加一倍,在刚性小的方向首先发生振动,转子振动的轨迹为椭圆。