固冲发动机补燃室非稳态沉积过程分析

采用随机轨道模型,基于动网格技术,分析了固冲发动机补燃室非稳态沉积过程中燃气流动边界形状变化对燃气流动及沉积分布的影响,得到了补燃室内壁面的非稳态沉积规律,对直连试验中出现的进气道出口拐角烧蚀沟槽现象进行了理论分析。研究发现,不同粒径的粒子在补燃室内的分布特性差异较大;凸起的沉积块可改变燃气流向,使进气道出口拐角处的高温区变大,易形成沟槽;点火初期,掺混区进气道之间的沉积速率较大。     

螺旋桨滑流对增升装置气动特性影响研究

采用准定常N—S方法,通过单独螺旋桨模型计算结果与非定常、试验结果的对比,验证了计算螺旋桨滑流的可行性。重点对带两台螺旋桨发动机的增升装置流场特性进行了计算分析,研究了螺旋桨滑流对飞机流场和气动性能的影响。研究结果表明：螺旋桨滑流明显改变了机翼表面的压力分布,使飞机升、阻力系数增大。为螺旋桨飞机设计提供了一种快速、有效的数值模拟手段,在工程领域具有一定的应用前景。     

二维Delaunay网格的一个约束边恢复算法

针对二维约束Delaunay网格生成中约束边恢复问题,提出了一个基于平面扫描策略的约束边恢复算法,证明其收敛性。给出了算法的计算复杂度,说明了算法能有效减少＂相交测试＂的次数。实现算法并用复杂算例验证了所得结论。     

三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。     

潜载导弹近筒口点火数值仿真

针对无横流影响的潜载导弹近筒口点火瞬态流场变化情况,建立多相流数学模型,考虑空穴模型,基于动网格技术,实现了潜载导弹弹射出筒、近筒口点火方式发射过程尾空泡的演变过程以及弹体、固壁的压强变化。仿真结果表明,筒口效应导致固壁各处的压强阶跃;近筒口点火后,尾空泡不断的颈缩、膨胀导致压强脉动,但不会产生断裂和回击现象。在筒口气团闭合后点火对筒口产生更大的压强脉动,有可能对固壁产生破坏。     

基于非结构网格的离子推力器羽流数值模拟方法研究

为适应复杂的航天器外形,文章采用基于非结构网格的数值模拟方法研究离子推力器真空羽流的基本特征。首先,介绍了非结构网格划分的方法,并分析了建模区域生成各种拓扑关系的特点;然后,提出了粒子穿越网格的具体算法,包括有利于减少计算量的结合矩阵重标号的带宽压缩算法;最后,开展了模拟结果和试验结果的对比研究,验证了以上数值模拟方法的可行性。     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

机翼-机身-短舱-挂架外形气动优化设计方法

采用基于径向基函数的无限插值方法进行了复杂外形的动网格生成,并针对其存在的可能导致壁面附近出现负体积与网格质量下降的问题,提出对不与物面直接相连的边的位移采用径向基函数插值,对与物面相连的边的位移进行线性插值获得位移量的方法解决该问题.采用无限插值方法建立了一种简单有效的物面相贯线处理方法.利用离散共轭方法计算目标函数梯度,对DLR-F6 机翼-机身-短舱-挂架外形的机翼和短舱进行了几何外形参数化与气动外形优化设计.结果表明:考虑短舱安装方式的优化设计较不考虑短舱安装方式的优化设计可降低大约0.0001的阻力系数.通过全机优化设计,将全机阻力系数降低了0.00153.      

利用多状态试车数据修正发动机部件特性

针对经典的航空发动机部件特性修正因子法的不足,提出一种于多状态试车数据的修正方法.利用发动机部件匹配模型方程以及试车参数计算值和实测值的差值组成目标方程组,采用改进的免疫粒子群算法优化求解获得部件特性修正系数.以多个不同状态下的修正系数为基础,利用滑动最小二乘方法拟合修正系数曲面,进而得到修正的部件特性图.数值试验表明:综合利用不同状态下发动机试车数据进行部件特性修正,克服了传统方法的不足.提高了修正部件特性在整个工作范围的精度.      

确定时间序列均值和方差函数的移动多点平均方法

为了提高有限样本或小样本情况下时间序列均值和方差函数的确定精度,以保证时序建模、分析和预测精度,提出一种确定序列趋势项的移动多点平均方法.该方法能够得到时间序列均值中非周期部分,结合样本周期图法得到的周期项,可综合得到其均值函数,并可进一步得到时序的标准差函数.Monte Carlo模拟对比结果表明:该方法不仅能够保证时间序列段内分析精度,而且能够有效提高时间序列的预报精度.