电大尺寸散射体的RCS计算方法研究

为了更加有效的求解电大目标散射计算问题,在避免谐振区效应情况下获得复杂目标的隐身雷达散射截面(RCS)特性,引用混合场积分方程(CFIE),在快速算法的求解迭代过程中采用共轭梯度算法(CG)的收敛技术,能够稳定的求解电大尺寸的RCS。计算和分析了金属球双站RCS和金属立方体双站RCS,并和精确计算、相关参考文献进行了比较,证明了方法的正确性,计算结果稳定,在工程实际上有较大的应用价值。     

基于欧拉方程的跨声速翼型设计

给出了一种基于欧拉方程的跨声速翼型设计方法。方法以Ｔａｋａｎａｓｈｉ提出的“正反迭代余量修正”设计原理为基础，在气动力分析模块中，以欧拉方程为基本控制方程，并采用Ｗａｌｚ边界层方法对其进行粘性修正；反设计模块采用经过改进的二维翼型设计方法。方法的程序经过几个设计实例主题证明是有效而实用的。     

关于三维Navier－Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵［Ａ］和［Ｂ］，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程计算（即Ｍ∞＝３．０，α＝０°和设计工况Ｍ∞＝２．６５，α＝０°），获得了满意的全场结果；对于Ｍ∞＝２．６５的设计工况，同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的Ｎ－Ｓ求解能在普通微机上进行。     

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

高速风洞支架干扰数值修正研究

提出了一种计算高速风洞支架系统对飞行器模型纵向气动力干扰量的数值计算方法 ,从跨声速全位势积分方程出发 ,编制了适用于飞行器全机模型及其带支架情况下的跨声速绕流计算程序。通过对双垂尾模型和GBM 0 3模型两个算例的计算 ,讨论了尾支撑位置及其几何外形参数对模型气动力的影响 ,并对GBM 0 3模型带短支杆情况下的纵向实验结果进行了修正。表明该方法对于分析研究风洞模型支架干扰问题并进行支架干扰修正是可行的、有效的 ,可以作为选择尾支撑位置及其几何外形参数和对跨声速风洞纵向实验结果进行支架干扰修正的工具。