连续变弯度翼型动态气动特性数值模拟

针对翼型后缘连续变弯度运动中后缘边界精确数值模拟的问题,提出了基于二维多项式的时空曲面拟合方法,实现了对后缘边界时空位置的精确模拟。在此基础上,基于OpenFOAM发展了翼型后缘连续变弯度与大幅度俯仰运动耦合的运动边界数值模拟,并计算了翼型耦合运动的气动力,讨论了后缘线性/非线性变形对翼型大迎角动态气动特性的影响规律。结果表明:后缘运动对翼型俯仰运动的升阻特性影响显著,特别在翼型大幅度俯仰时后缘非线性变形对升阻特性改善效果比线性变形大6%~10%。同时还研究了翼型俯仰与后缘变形运动相位差对气动特性的影响。特别地,当相位差为180°时,后缘运动使动态失速时的最大升力提高50.3%,平均升力提高34.6%;当2种运动相位差为0°时,后缘运动使动态失速时的最大阻力降低39.7%,平均阻力降低30.2%,最大升阻比提高22.3%,平均升阻比提高16.8%;同时,翼型在动态俯仰过程中出现负阻力现象,对产生负阻力的原因进行了分析。这些结果可用于指导连续变弯度后缘控制律的设计。      

基于k-ωSST模型的同心筒发射装置流场数值模拟

围绕同心筒发射装置流场数值模拟及结构优化展开研究,流场数值模拟基于轴对称Navier-Stokes(N-S)方程、k-ω shear stress transfer(SST)湍流模型,并运用域动分层动网格方法,以同心筒初始方案为参照,研究了同心筒发射过程的动态流场,在此基础上,分别对内筒收敛段的收敛角度以及筒口扩张段的扩张角度进行了优化研究,研究结果表明:内筒收敛段使得燃气排导更为通畅,筒内的热环境有所改善,并且内筒收敛段起到了一定的防喷作用,收敛角越大,筒底热环境越好;筒口扩张段对发射过程动态流场也有明显影响,收敛角为15°,扩张角为20°时,导弹从开始发射到完全出筒过程中热环境最为良好。      

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

环月飞行器返回轨道设计与优化

针对传统月地返回轨道设计方法在从初始轨道直接转移的局限性，提出一种基于Lambert问题的月地转移轨道设计方法。该方法重新选取了轨道约束参数，并在笛卡尔坐标系下建立了返回轨道动力学模型，并基于该模型全面分析了返回轨道关于约束参数的变化特性。此外，针对带有约束条件的返回轨道优化设计问题，提出了快速收敛的函数构造法，并给出了有效的目标函数形式。以某一在轨环月飞行器的轨道参数作为初始轨道进行返回轨道计算，结果验证了所提出返回轨道设计方法以及函数构造法的有效性。     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型，对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟，对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价，清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构，即直角模式、对角模式和螺旋模式。随着截面尺寸的减小，爆轰波胞格图样与波阵面结构随之产生适应性的变形，乃至出现胞格数目的减少和部分三波线的消失；最终，直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式；临界尺寸下，得到了转换过程中新的过渡结构和压力脉冲振幅的变化，并以此初步讨论了直角和对角传播模式的稳定性。     

格序群的极子群在凸l－子群格中生成的闭子格

设Ｇ是一个格序群，Ｐ（Ｇ）是Ｇ的极子群集合，Ｃ（Ｇ）是Ｇ的凸ｌ－子群集合，Ω（Ｇ）是Ｐ（Ｇ）在Ｃ（Ｇ）中生成的闭子格，本文通过对Ω（Ｇ）的研究得到了一系列重要结果，对有限值格序群给予了一系列新的刻画，证明了ｇ是基元当且仅当ｇ－是ｇ的唯一的Ω－值等结论，实际上对Ｐ．Ｃｏｎｒａｄ提出的问题＂Ｄｅｓｃｒｉｂｅｔｈｅｅｌｅｍｅｎｔｓｉｎ＂进行了富有意义的探索。     

基于SOVA的并行级联空时格码的译码算法

并行级联空时格码调制(Turbo—STCM)将并行迭代译码结构(Turbo)应用于空时编码，可以极大地提高空时格码的编码增益。本基于一种改进的(Soft—output Viterbi a1gorithm，SOVA)译码算法，提出了针对Turbo—STCM的并行迭代译码方案。仿真结果表明，在保持分集增益的前提下，Turbo—STCM和原始的空时格码调制相比，编码增益有了极大的提高．同时采用该SOVA算法译码复杂度低，速度快，更利于实际系统的开发实现．     

飞机型架参数化设计系统的开发与应用

本文介绍了AutoCAD二次开发飞机型架参数化设计系统的结构、功能、应用。说明对AutoCAD进行二次开发使得工装型架设计绘图效率大为提高。     

基于全线程树数据结构的笛卡尔网格高效生成技术

笛卡尔网格方法具有易于自适应、自动化程度高、网格质量好等优势。由于其非贴体特性,在处理复杂构型或者复杂流动问题时候往往存在网格量过大的问题,不容忽视。笛卡尔网格生成的效率直接影响了整个计算周期的长短,有必要发展高效的生成技术。对于笛卡尔网格而言,决定网格生成效率的关键在于网格数据结构,其直接影响计算量和存储量。针对三维构型进行笛卡尔网格生成,发展了邻居查询更快捷、内存利用率更高的全线程树数据结构,并在本文的方法框架下进行了适应性应用和改进。同时,为了高效地判断网格单元类型,构建了物面单元的快速检索方式,并引入了染色方法,进一步提高网格生成效率。还提出了一种鲁棒的奇异性检测算法,保证网格单元类型判断的鲁棒性。在流场解自适应方面,采用的是速度散度和旋度相结合的三维判据,以保证对于多种流动特征的捕捉能力。通过圆球、导弹、翼身组合体、机翼-副翼等算例进行了考核验证,经对比,网格自适应位置与理论解吻合较好,且网格单元生成耗时短、平均耗时受物面网格分布影响小,证明了方法的可靠和高效性。     

气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.