用高精度紧致差分格式分块耦合求解二维粘性不可压缩复杂流场

本文提出了张量形式的粘性不可压缩流场方程的高精度分块耦合求解方法。在分块算法中用重叠一层的方法来实现子块间的数据交换，同时克服了块剖分所造成的数值奇性。利用本文所提出的算法成功地模拟了椭圆柱绕流问题的时间发展过程，如尾流，旋涡的产生，脱落过程以及Karman涡街的形成等。     

一种用于捕捉间断多维任意单元的有限元法

本文根据差分格式构造过程中对三阶、四阶人工粘性及熵增条件的分析，构造出适合于非结构单元的有限元法。计算实例表明，这种方法能较好地捕捉间断，激波前后的伪振荡基本消除。     

一种新型颤振激励系统的特性研究

在非结构网格的基础上，以欧拉方程为流场的主控方程，采用格心格式空间离散的有限体积法和四步Runge—Kutta显式时间推进的方法．对一种新型颤振激励系统的流场进行了数值研究。研究了该系统的工作原理和激励力特性。研究结果表明，该系统能很好地激励飞机结构。     

用非结构网格和欧拉方程计算多段翼型绕流

用Ｄｅｌａｕｎａｙ法划分三角形构造非结构网格时，如果布点不当，会产生边界穿透，使网格生成失败，本文分析了几种布点情形，提出不容易产生边界穿透的布点方式。这一原理已用于生成绕多段翼型的非结构网格。     

二维网格生成技术及应用

论述了采用求解偏微分方法与代数方法相结合,生成二维网格的技术及其在内流场和外流场中的应用.这种方法便于计算程序的编制,具有较好的通用性.     

一种制造网格资源集成与发现方法

提出了一种基于功能分类的资源集成和层次分布式与对等网(P2P)混合模式的资源发现方法.以某军机主机厂的部分资源为例开发了资源模型和内部资源发现原型系统,并通过一个基于P2P的制造网格仿真系统验证了外部资源发现算法的有效性.     

一种新的计算非稳态流动的自适应无网格方法

基于带权点填充布点方法,结合指定流场参数梯度的分布,提出了一种新的无网格自适应方法计算非稳态问题。该方法的自适应探测器 E 由节点的权重和参数梯度组成。使用自适应探测器识别出加密或者粗化的区域,形成加密空腔和粗化空腔,然后重新填充布点。新节点的权重由自适应探测器 E 结合该处的梯度大小通过预测-校正迭代算法计算得到。另外,由于参数梯度在激波附近波动范围比较大,所以新节点的权重存在最大值和最小值。首先对预先设置压强梯度的流场区域进行自适应布点,并与传统自适应布点结果进行对比,点云分布图显示所提自适应方法重新生成的点云结构疏密分布更加合理。进一步,将此自适应方法应用于 Riemann 问题和激波碰撞圆柱问题,计算结果表明该方法可以在节点数目较低的情况下显著提高运动激波的分辨率。在激波碰撞圆柱问题中,比较了使用自适应算法和非自适应算法得到相当的结果所使用时间,前者是后者的20.6%,因此该自适应算法在计算效率方面也具有较大的优势。     

关于 CFD 高精度保真的数值模拟研究

本文由四部分组成。第一部分研究了利用高阶格式、网格、物理模型和高性能计算机求解 NS 方程时,计算数据的高精度保真性的要求以及高精度保真的指标。第二部分研究了利用大规模计算开展大涡模拟和直接数值模拟存在的问题。第三部分分析了在高雷诺数下 NS 方程计算方法和网格的关系。最后提出了建立计算数值验证、确认的新方法。     

喷管结构形式对两级PDE性能的影响分析

为了研究喷管结构形式对两级 PDE 性能的影响及作用规律,以氢气和空气混合物为例,采用 FLUENT 软件中的 k-ε湍流模型,使用非平衡壁面函数、PISO 算法及基于梯度的动态自适应网格加密方法,对收敛、扩张、收扩等不同喷管结构形式的凹面腔内环形向心射流聚心碰撞产生激波会聚起爆爆震波的过程进行了数值模拟,并对两级脉冲爆震发动机产生的推力和冲量进行了计算。通过对不同时间点凹面腔产生的推力和冲量、喷管产生的推力和冲量、两级 PDE 总的推力和冲量三者之间的对比和分析,发现凹面腔产生的推力和冲量是两级 PDE 总的推力和冲量的主要来源,不同喷管结构形式对两级 PDE 推力和冲量性能的影响有较大的差别。以冲量为例,收敛喷管可以提高凹面腔的冲量,但是由于自身同时产生过多负的喷管冲量,因而总冲量较小;扩张喷管虽然对凹面腔冲量的提高作用不大,但是由于其自身可以提供正的喷管冲量,所以总冲量较大;收扩喷管的收敛段也可以较小幅度的提高凹面腔冲量,但同时产生的负喷管冲量会和扩张段产生的正喷管冲量发生抵消,使总的喷管冲量比较低,发动机总冲量也较小。综合考虑各方面因素,扩张喷管性能最优。可见在两级 PDE 喷管结构形式的选择上,要同时考虑凹面腔冲量和喷管冲量,权衡两者,才能选择出合适的喷管。     

倾转旋翼典型飞行状态气动特性的CFD分析

针对倾转旋翼特有的悬停、巡航、过渡3种飞行模式,基于CFD方法分别对倾转旋翼的上述3种典型飞行状态进行了数值模拟.首先,基于运动嵌套网格思想建立了适合倾转旋翼运动模拟的多层运动嵌套网格系统.然后,基于RANS(Reynolds averaged Navier-Stokes)方程建立了适用于非定常流场分析的CFD模拟方法.为了提高整个流场的计算效率,时间推进采用高效的隐式LU-SGS(lower upper symmetric Gauss Seidel)方法,并运用SPMD(single program multiple data)并行计算模式.最后,采用所建立的网格生成方法和流场求解器,进行了不同飞行状态下的倾转旋翼气动特性数值分析,验证了建立CFD方法的有效性并得出了多种飞行状态下的气动特性.结果表明:悬停状态下倾转旋翼相对于直升机旋翼有更加明显的下洗流影响区域.巡航状态两旋翼中心间距为2.5R时,倾转旋翼之间存在明显的气动干扰.过渡状态下,倾转旋翼提供的拉力在倾转角为45°前下降较快,且在倾转角较小时,倾转旋翼桨叶表面拉力分布仍能表现出直升机旋翼的特性.