二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

双立尾对战斗机隐身特性的数值模拟

为了降低战斗机的雷达散射截面（RCS,Radar Cross Section）特性,研究了双立尾外倾对战斗机RCS特性的影响.对战斗机的三维数字样机进行网格划分,基于物理光学法数值模拟出双立尾向两侧对等偏转对三翼面战斗机的头向、侧向和尾向的RCS特性的影响,并对战斗机的缩比实体样机进行电磁测试,证明了方法的可行性和数值模拟的准确性.数值模拟结果表明,当双立尾从0°向外偏转10°时,三翼面战斗机侧向RCS值缩减至原来的 9.8%,而头向和尾向的RCS值变化幅度较小.      

三维椭圆形方程网格生成中的新源项修正法

构造一种新的源项修正方法用于三维椭圆型方程网格生成，可以实现对网格与边界正交及网格与边界间距的双重直接控制。     

鸭翼布局中双立尾对全机气动及流场特性影响

在战斗机先进气动布局研究中,双立尾位置的选择始终是一个十分重要的问题.不适当的双立尾位置会给飞机纵横向气动特性带来严重的影响.对一种鸭翼布局的飞机模型,按3种不同的双立尾配置进行了气动力测量、流态显示,然后用PIV(Particle Image Velocimetry)进行了不同迎角下的流场测量.结果表明:双立尾处于飞机内侧后置内移位置其最大升力系数具有最大值.破裂过程及流场特性同无双立尾时的情况十分相似,进而说明双立尾同机翼涡的干扰主要是促进了涡的提早破裂,从而恶化了全机气动特性.      

非结构重叠网格显式装配算法

针对重叠网格中洞映射法占用过多物理内存的问题,发展了一种改进型洞映射法;基于相邻单元搜索法,发展了一种基于相邻阵面的贡献单元搜索法;通过将割补法与隐式切割技术相结合提出了一种非结构重叠网格显式装配算法。该算法首先生成一套包围物面的笛卡儿网格,其次存储所有与物面边界相交的笛卡儿网格信息,最后根据所存储笛卡儿网格与所需判断的网格单元的相对位置来判断其是否为洞内单元。在成功判断出所有洞内单元后,以当前洞边界为初始阵面推进,同时以各个网格单元的物面距离为判别标准对重叠区域进行优化,生成最终插值边界。所提算法优化了传统非结构重叠网格装配过程,具有物理内存占用低,贡献单元搜索次数少以及计算效率高等特点。通过2个典型复杂流动算例验证了所提算法的准确性与适用性。      

多段翼型贴体网格的生成方法

本文给出了一种多段翼型绕流的分区贴体网格生成方法。采用泊松方程的解来给出贴体网格的坐标,泊松方程中的源项通过靠近边界的网格间距和正交性来决定,本文对分区网格的区域边界的选择作了讨论,采用生成的分区贴体网格计算了绕多段翼型位流的解,其结果与实验所得是很一致的。     

三维椭圆型方程网格生成中的新源项修正法

构造一种新的源项修正方法用于三维椭圆型方程网格生成，可以实现对网格与边界正交及网格与边界间距的双重直接控制．     

2D约束Delaunay剖分生成表面模型的表面网格

提出一种生成表面模型的表面非结构化网格的方法.该方法仅假定表面模型的每个面片在取定适当的投影平面时可看成是一个单值函数.对表面模型的每个面片,首先在其相应的投影面上进行二维约束Delaunay剖分,然后对网格顶点进行插值.由于采用联动剖分的策略使得单独生成的各个表面网格在共同边界处匹配.与通常的生成表面网格的网格前沿法相比,本方法无需预先对域的边界进行离散,边界的离散体现在约束Delaunay剖分中恢复限定线段的边界细分过程中,减少了用户干预.通过合理地指定2D约束Delaunay三角化时网格单元尺寸分布函数,可以有效控制表面网格对表面模型的逼近精度和实现自适应的表面网格.实验结果证明了本算法的有效性.     

带副翼的翼身组合体绕流的 Euler方程解

提出了一种求解带副翼偏转的翼身组合体绕流的Euler 方程计算方法.将对接分区网格与分区求解算法相结合,有效地求解了绕此外形的复杂流动.提出了一种满足通量守恒的内边界耦合条件.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种限制器(Limiter)函数以保证TVD性质.数值算例表明本文方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

四元数在战斗机飞行仿真中的应用

阐述了在战斗机的飞行仿真中采用四元数的原因和四元数的几何意义，并按照中国战斗机飞行仿真中采用的坐标系，介绍了三个问题，四元数运动学方程；四凶数与坐标变换矩阵的关系；用四元数表示的欧拉角；与此同时给出了一组数字仿真曲线，最后讨论了四元数法优于欧拉法之处。     

推力矢量和二元喷管

推力矢量是第四代战斗机获得过失速机动性的必要条件,而二元喷管又是推力矢量控制中性能较为优越的一类喷管。本文较全面地阐述了目前世界上推力矢量和二元喷管的发展和现况、二元喷管的构造和性能、二元喷管对改进飞机升阻特性和机动的作用。推力矢量能显著提高飞机在空战中的作战效率以及推力矢量控制所需的效能也在文中作了介绍。     

矢量喷流下平尾偏转对飞机气动性能的影响

通过系列参数试验,研究有、无矢量喷流作用下飞机平尾偏转对飞机气动特性、操纵特性以及飞机绕流/矢量喷流之间干扰特性的影响.实验结果表明:在亚音速中小迎角下,有、无矢量喷流时平尾偏转并不影响飞机的纵向静稳定性,只是改变飞机的零升力矩系数m_z0和零升力迎角α₀,但平尾偏转可以较大幅度地改善前体机翼的绕流形态,减缓机翼涡及边条涡的破裂.此外当飞机处于失速状态时,矢量喷流对于涡破裂失速流动表现出强烈的干扰作用,对飞机的气动特性产生了较大的影响,其干扰区域不仅局限于飞机后体,而且还延伸至前机体,该有利气动干扰量可以达到飞机气动力的10%以上.     

非对称空战体系对抗研究

对战斗机分代及其空战能力评估进行了研究,提出了空战体系对抗中,交战双方武器装备非对称度的概念和计算方法.基于体系对抗物理描述和数学模型,发展了非对称战役优势参数的概念.通过对典型的对称、非对称和极限非对称战役想定的作战模拟和分析,探讨了战斗机武器装备的数量与质量关系问题,装备质量太差就很难用数量去弥补.计算结果表明考虑非对称度影响的分析更接近当今战争的实际.      

一种战斗机外形涡流场Euler方程数值模拟

 用三步显式格式时间推进求解物理空间曲线网格上有限体积离散的Euler方程,数值模拟一种战斗机外形的涡流场.计算的机翼表面压力分布与实验符合较好,用Euler方程捕获到机翼前缘分离涡(主涡)及翼面上的二次分离涡.Euler方程解中出现的二次涡可能是逆压梯度、尖锐边缘和人工粘性共同作用的结果,它的流谱与实验定性符合.     

交叉航路航空器碰撞风险评估

在与经典的Reich模型进行比较后,对一种基于事件的交叉航路航空器碰撞风险模型进行了改进,补充考虑了处于平飞阶段的航空器因在垂直方向上相对运动而导致的碰撞风险.分析了其水平碰撞概率较Reich模型所得结果保守的原因,采用蒙特卡洛仿真的方法产生了两架航空器之间的相对位置误差数据,针对多个概率分布模型进行了拟合,比较了不同概率分布模型对航空器碰撞风险的影响.结果表明:代表严重误差的尾部数据直接决定了航空器碰撞风险,但因为其数据较少,所以除广义帕累托分布外,其他分布都不能很好地对尾部数据进行拟合.