直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

瓦状特征型面塞式喷管三维计算与实验

为了进一步了解瓦状塞式喷管的性能,采用NND差分格式求解三维N S方程和空气冷流对6单元瓦状特征型面塞式喷管进行了数值模拟和实验研究。研究模型的内喷管面积比为4,总面积比为40,设计压强比为1047。计算得到了流场马赫数和塞锥表面压强分布、喷管推力系数效率,以及不同压强比下中心平面、过渡平面和边缘平面的塞锥表面压强变化规律。计算结果与实验数据吻合得较好,效率数值最大相差1%。实验塞式喷管最大的推力系数效率为0 995,同钟型喷管相比,具有很好的高度补偿能力:从地面到高空,效率在0 93～0 995之间变化。和以前简化型面的4单元瓦状塞式喷管相比,实验和数值模拟均说明塞锥特征型面的优化设计提高了喷管性能,更充分体现了塞式喷管的高度补偿特性,可以成为未来工程应用的选择方案。     

“瓦”状塞式喷管的数值模拟与实验

为了寻求高性能和更接近工程应用的发动机，提出了一种内喷管为轴对称喷管，塞锥为凹面的“瓦”状塞式喷管，分析了这种塞式喷管的优缺点，并针对一研究模型进行了数值模拟和实验比较，数值模拟采用NND格式求解曲线坐标下的三维平均雷诺的N－S方程，并用k-ε两方程湍流模型封闭方程组，实验研究采用酒精和氧气作为推进剂进行了热试车；研究模型的内喷管面积比为3．24，总膨胀比为22．15，设计压力比为220，结果显示“瓦”状塞锥改善了塞锥的流场，并且当压力比在16．8－220的范围内变化时，其相对理想喷管的喷管效率在0．90－0．96内变化，对发动机设计作进一步改进，其性能有望进一步提高。     

塞锥型面简化与截短对塞式喷管性能的影响

用NND差分格式求解N－S方程，对截短直锥塞式喷管和用特征线法设计的截短型面锥塞式喷管的流场和性能进行了比较。成功地捕获了流场激波、膨胀波及反射压缩波，计算结果与实验数据吻合较好，对不同型面塞式喷管在不同高度下的推力性能和塞锥截短对发动机性能的影响。研究表明，相同长度的截短直锥发动机总推力比型面锥塞式喷管发动机低1％－1.5%；相对30％截短型面塞锥发动机，50％截短型面塞锥发动机的总推力可提高约2.5%。     

求解可压缩N－S方程的无粘－粘性分数步法

用按物理过程作时间分裂的无粘－粘性分数步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

烧蚀反应动力学参数的确定

本文提出了一种理论模型与实验相结合确定碳基材料烧蚀反应动力学参数(活化能和指前因子)的方法。首先建立了碳基材料热化学烧蚀的湍流边界层模型,运用数值方法对边界层守恒方程组、表面烧蚀和固相热传导进行耦合求解,计算过程中不需要事先假设烧蚀的控制类型。其次,在电弧加热器中模拟火箭发动机喷喉的烧蚀环境进行烧蚀材料与氧化性气体反应的烧蚀实验。最后,运用所建立的模型和在电弧加热器上的实验结果进行预测比较,确定烧蚀反应的动力学参数。作为实例,确定了石墨与CO_2的烧蚀反应的动力学参数。     

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

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

喷管冷吹风实验的推力架设计

本文介绍了柔性推力架和刚性推力架在喷管风洞实验中的表现.讨论了: 为什么风洞实验中使用柔性推力架达不到推力测量的精度要求?为什么试车台上用柔性推力架又是可行的.通过分析,找到了问题的根源在于风洞实验中用的是气体介质,而试车台上用的介质是液体推进剂.