欧拉方程多重网格方法及三维跨声速大迎角流动数值计算

本文给出一种粗、细网格时间步长可调整匹配至最优组合的多重网格方法以获得在最细网格上最大的加速收敛效益，二重网格中节省机时可达３５％。文中的数值方法可用于飞行器跨声速大迎角无粘有旋的涡干扰流场的模拟研究。     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Ｅｕｌｅｒ方程定常解质量的非结构网格自适应方法和模拟结果。流场算法为结点型有限体积法和Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式，网格用阵面推进法生成，生成初始网格时使用均匀的直角坐标背景网格。     

在非结构网格上求解非线性航空声学问题的高精度有限体积法

介绍了有限体积数值方法,该方法适用于在任意非结构网格上求解线性和非线性的航空声学问题.本方法基于角点-中心的多参数格式,可在笛卡尔网格上达到六阶精度,对于可能的不连续性采用了自适应耗散.通过一系列算例研究了该方法的特性,结果表明:在模拟谐振型管中的噪声抑制中,所提出的方法是很有效的.     

三角翼低速动态大攻角气动特性试验研究

对后掠角分别为X=60°、70°和80°尖前缘平板三角翼模型和一个前缘后掠角为（76° 40°）的双三角翼模型在低速风洞中作大攻角俯仰摩波运动，振幅am=30°、60°和90°，缩减频率K=0.01～0.12，基于根弦长雷诺数Re=2.76×10^5～8.23×10^5。进行了六分量动态气动载荷测量，动态流动显示和70°三角翼上翼面非定常压力测量，并分别与对应的静态试验结果比较。分析了运动参数包括缩减频率、振幅和Re数、后掠角对气流动态迟滞特性的影响。     

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

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

无波动,无自由参数的耗散差分格式（NND）在喷流计算中的应用

本文用文献[1]建立的无波动、无自由参数的耗散差分格式(NND格式),对无粘喷流流动进行了数值模拟,并就一典型的自由射流流动,与国外高精度、高分辨率的PLM差分格式(Piecewise-Linear Method)的计算结果和实验结果作了比较,得到了很一致的结果,且所用的计算时间少于PLM格式。     

高模碳纤维表面冷等离子体处理技术

高模碳纤维复合材料有较高的比强度、比模量，但纤维和基体树脂间的粘结性差，使得复合材料层间剪切强度很低，限制了复合材料优异性能的发挥。     

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

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

大曲率弯道和叶栅内的湍流数值模拟

本文在贴体曲线坐标系下应用非交错网格技术，发展了适于求解复杂流动的压力修正算法，采用速度的协变物理分量为动量方程的求解变量，推导得出一个适于任何曲线坐标系下控制单元体界面流通量的插值模式，可以有效地抑制不合理的波动压力场，计算表明，算法稳定收敛，计算值与试验值吻合良好，同时能预示通道内复杂的流动现象。