气体涡轮燃烧室流动与燃烧三维数值模拟

用三维湍流N -S方程和单步快速不可逆化学反应描述涡轮燃烧室内的流动与燃烧过程。采用非交错网格中的SIMPLE算法求解控制方程并对涡轮燃烧室内的流动与燃烧过程进行三维数值模拟 ,得到了有代表性的圆筒形燃烧室内燃气参数的详细分布情况。模拟结果表明 :对圆筒形涡轮燃烧室内的流动和燃烧而言 ,预计结果与测量结果趋势相同。本文是研究涡轮燃烧室的初步尝试     

CS-01高空台推力测量和校准装置研制

推力是涡喷涡扇发动机高空模拟试验的重要测量参数之一。CS－01高空台原推力测量采用杠杆、砝码、测力秤系统，现采用原位校准，液压加载、应变式推力传感器来测量和校准，从而消除了“高度差”的影响，降低了稳态推力测量系统的测量不确定度。     

LRE模型燃烧室燃烧过程数值模拟

建立了任意斜交曲线坐标系下液体火箭发动机（ＬＲＥ）内部工作过程的气液两相湍流化学反应流模型。应用该模型对气氢液氧模型燃烧室内部燃烧过程进行了数值模拟，得到了全流场的速度分布图形、燃气组分等值线、流场等温线、流场等压线、流场等马赫数线、壁面温度与径向平均温度曲线和燃烧效率曲线。计算结果表明：用数值模拟方法分析ＬＲＥ内部工作过程是可行的。     

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

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

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

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

环型超声速空气引射器零二次流流场数值研究

采用二维轴对称雷诺平均方程和标准κ-ε双方程湍流模型，数值研究了环型超声速引散器零二次流的流场结构及盲腔压强的变化，空间上采用二阶迎风格式进行耦合求解，时间上采用显示的Runge-Kutta方法进行迭代推进。结果表明，引射器几何参数不变的情况下，启动后的盲腔压强与引射气流总压之比为一常数；喷管马赫数不变情况下，喷管出口面积与混合室入口面积比越小，盲腔压强越低，扩压器性能越好，启动要求的总压越低，对超声速空气引射器的设计具有指导意义。     

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

激光微推进器用于微小卫星的姿轨控仿真

为研究激光微推进器用于微小卫星在轨空间飞行任务的可行性,利用已建立的激光推进微小卫星仿真平台,对微小卫星在轨姿态调整、轨道维持和轨道转移等飞行任务进行了仿真和分析。仿真结果表明,星载激光微推进器可满足在轨微小卫星的高控制精度姿态调整、近地观测轨道维持以及较长时间要求的轨道转移等空间任务对推进系统的要求。     

激光冲击残余应力场的有限元分析

建立模拟激光冲击残余应力场的非线性弹塑性有限元模型,预测AISI 304奥氏体不锈钢厚、薄靶材经激光冲击后表面以及内部的残余应力场分布.模拟过程中,考虑材料高应变率(106s-1)下的力学性能和激光诱导冲击波的精确加载,分析激光功率密度、激光光斑尺寸、激光脉宽、激光冲击次数、激光单面与双面冲击等激光冲击参数以及初始残余应力等因素对不锈钢靶材残余应力分布的影响.数值模拟结果与X射线衍射法测量值吻合较好.在有初始残余拉应力(250MPa)存在的情况下,激光冲击仍能使304奥氏体不锈钢靶材形成残余压应力层,说明激光冲击工艺可提高奥氏体不锈钢焊接构件的抗应力腐蚀开裂能力.