基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

旋翼桨叶翼型流场和跨音速翼涡干扰噪声研究

基于桨涡干扰流动的噪声场分析了NACA0012翼型的远场气动噪声.采用Euler方程计算了NACA0012、NACA23012、OA207三种翼型的绕流场,选取Kirchhoff面选为包围整个翼型的一层网格,通过得到NACA0012翼型的Kirchhoff面上在桨尖涡干扰流动下的近场解,根据Kirchhoff方法求解远场声压.样例计算分析体现出:在计算桨涡干扰噪声时,旋翼CFD技术与Kirchhoff方法相结合的良好效果.     

燃油直接喷射入横向气流中雾化机理数值研究

采用商业软件Fluent对燃油直接喷射入横向气流中的破碎雾化过程进行数值模拟,选用欧拉-欧拉方法中的VOF（volume of fluid）耦合level-set方法模拟气液两相流动.数值模拟了横向气流中燃油直接喷射射流破碎雾化过程及其相应流场信息,分析了流场结构对液柱破碎雾化过程的影响,获得了横向气流中射流破碎点位置以及射流穿透深度,并与试验数据进行了对比.数值计算结果表明:①采用VOF耦合level-set算法能够较好地模拟射流喷入横向气流中燃油破碎雾化过程;②由于反向旋转的涡对的存在,使得液膜和液滴沿着展向（y向）输运,对于液滴的二次雾化起到促进作用;③射流液柱总是在喷嘴下游约8倍直径处开始破碎,且破碎点位置取决于喷嘴几何特性.      

基于平均扩展速率的裂纹扩展模型

不同于基于线弹性断裂力学（LEFM）的经典裂纹扩展模型,基于平均扩展速率的裂纹扩展模型通过形式简单的经验公式来描述随机载荷谱下单位飞行小时的平均裂纹扩展,该模型将随机谱当作"当量常幅谱"处理,模型中仅有两个特征参数,反映了所有复杂的载荷顺序效应及其他影响。基于该模型,进一步提出速率类比法则,可基于试验谱下扩展速率对未试验谱下扩展速率进行合理预测。首先介绍了基于Frost&Dugdale公式和基于Paris公式的平均扩展速率模型的内涵及参数推导方法,然后通过典型机身框关键部位模拟件在多种随机谱下的疲劳试验数据对平均扩展速率模型和速率类比法则进行试验验证。经对比分析,推荐使用基于Paris公式的平均扩展速率模型。经验证,采用基于Paris公式的模型和速率类比法则,可实现随机载荷谱下裂纹扩展速率和寿命的高精度预测。同时,速率类比法则的准确度与理论参数推导时选取的经典LEFM模型有关。平均扩展速率模型和速率类比法则为外场飞机单机寿命监控提供了一种可靠的疲劳寿命预测方法。     

基于Euler方程的直升机旋翼翼型反设计方法

直升机旋翼流场的求解和贴体网格的生成是进行旋翼翼型反设计的关键,本文采用Euler方程求解直升机旋翼翼型流场,并通过求解泊松方程生成围绕翼型的网格.在此基础上,结合线化二维亚音速和超音速流的小扰动理论,采用MGM方程作为反设计方程,基于有限差分方法建立一套直升机旋翼翼型反设计的迭代方法.在给定的目标压力分布条件下,与CFD方法相结合,反设计得出满足要求的翼型,并与目标翼型吻合良好,表明本文建立的翼型反设计方法能够有效地用于直升机旋翼翼型设计.     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

基于RBF插值边界条件的笛卡尔网格方法

介绍了一种基于RBF(Radial Basis Function,径向基函数)插值处理物面边界条件的笛卡尔网格方法.整个流场采用有限体积方法计算,但是对于物面附近切割的网格计算域不完整,因此通过RBF插值来得到物体内部实单元上的值,从而使得流场计算顺利的进行推进.相比于其他插值方法来说,RBF插值特别适合于多维离散点的插值,不会因为插值点的位置关系而产生病态的系数矩阵,而且使用比较少的网格点就可以得到同最小二乘方法相当的精度.将这种方法应用于圆柱和翼型绕流中,得到了满意的结果,证明了方法的有效性,还对边界条件的各种不同形式对结果的影响作了对比.     

固体火箭发动机的一种点火药量估算方法

针对传统经验公式估算点火药量存在的不足,从固相点火理论出发,推导得到固体发动机点火药量的计算公式.对不同配置的固体火箭发动机进行点火药量验算的结果表明,该公式估算的点火药量与发动机真实试验药量基本一致,估算精度明显优于传统经验公式,且公式通用性较好.      

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。