复杂组合体亚音速气动力计算改进的面元法

发展一种低阶面元法计算复杂组合体的亚音速气动力。本方法采用两种形式的速度势的内部第一类边界条件，通过对速度势的数值微分得到速度，进而计算气动力参数。本方法克服了早期低阶面元法对面元剖分敏感和高阶面元法的运算量大的缺点，兼有高阶面元法的精度高和早期低阶面元法的运算量小的优点。用本文给出的方法，计算了多个模型（包括型号）的气动特性，计算结果和实验结果有很好的一致性，表明本方法灵活、稳定、准确和有效。     

多段翼型分析和反设计位流计算

本文采用第二代低阶面元法给出了多段翼型分析和反设计位流计算方法。按照曲率变化情况划分板块可很好地提高计算效率。     

计及边界层汇流效应的多段翼型失速特性计算

 本文给出计及边界层汇流效应的多段翼型失速特性的解法。用高阶奇点分布面元法（Panel Method）求位流解,然后进行各翼段粘性尾迹形状迭代,并解出正常边界层和汇流边界层特性。当翼段上有后缘分离时,还要确定分离尾迹的形状,用位移厚度当量源（汇）模拟粘性效应。进行粘/位流迭代直至收敛。上述方法在超过多段翼型失速迎角时仍然有效。计算结果与实验数据比较,符合良好。     

亚声速机翼快速摄动面元法

在内部Dirichlet问题提法的面元法基础上，用解析法求几何外形摄动导数（偶强密度对机翼外形坐标的敏感性偏导数），从而快速确定机翼表面压力分布、升力系数和俯仰力矩系数。由于本文以偶强密度常值分布的低阶面元法为基础，故比以偶强密度二次分布的高阶面元法为基础的摄动面元法在最费机时的偏导数矩阵计算上要快一个量级，而由于以符合流场特性的物理内插代替加权几何内插，两者在确定物面压力分布时准确度却基本一致，此外，本方法对内存要求较低，可处理较相应高阶方法更多的面元数。     

Dirichlet边界条件低阶面无法的改进

对基于内部Ｄｉｒｉｃｈｌｅｔ边界条件的低阶面元法的速度分布提出一个改进算法。由极限过程结合高阶奇异积分的解析计算直接计算物面速度分布取代传统的偶强梯度法,以避免偶极强度边缘外推的困难。改进算法改善了ＶＳＡＥＲＯ分析中诱导阻力系数估值过低的问题。对后掠翼、前掠翼和翼身组合体等算例的计算表明,计算结果的准确度与高阶面元法相近。此外,直接积分法还可正确地计算外部流场速度分布。这是偶强梯度法不能完成的。     

M=1音速流翼型设计计算

来流马赫数为1.0的翼型绕流型态不同于其它跨音速时的流动型态。在数学上,它的远场边值呈奇异性。本文利用Chaplygin转换推导出精确的速度面流动方程,又根据小扰动理论,采用Tricomi方程的奇点近似解析解建立了对应在速度面上来流奇点的边界条件。在用有限差分法进行数值计算时,创建一种新的差分格式,并在不同速度区采用不同座标系统进行线松驰迭代计算,计算表明,其结果是稳定的和收敛的。本文用Spline方法确定出五个仟意的速度面上的边界曲线,经计算,在获得五个音速流流场的同时,也得到五个封闭的音速流翼型。     

多段翼型前缘缝翼位置的优化设计研究

本文讨论优化前缘缝翼位置的位流设计方法。优化变量为缝翼相对于主翼的缝隙(Gap)、覆盖量(Ouerlap)和偏角δ_s,目标函数为主翼上的压力峰值。应用高阶面元法计算多段翼型压强分布。用Powell优化法使主翼上压力峰值减至最小,以延迟多段翼型的失速,增大最大升力系数。本方法已用于计算NACA64A010两段和四段翼型以及Foster三段翼型,所得结果与实验数据和位流/边界层耦合设计法的结果有很好的一致性。     

多段翼型失速分离特性的计算

本文用高阶奇点分布板块法(Panel Method》及边界层理论,对多段翼型的失速性能进行计算。采用分离尾迹的模型,按尾迹是流线用迭代法求出分离尾迹边界(此为内迭代)。在尾迹上配置奇点,然后把尾迹考虑在内进行粘/位流迭代(此为外迭代)。收敛后可得多段翼型空气动力特性,且可计算到失速以后的情况。算例有单段翼型GA(W)-1,NACA 4412,和F111飞机的三段翼剖面,所得结果与实验值能很好地符合。     

用粘流—无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进，改进的半反－局部联立耦合方法。考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好，算例表明，翼型高升力状态计算必须包括尾流的作用，也应用对势流计算压力与实际压力之间     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

旋成面叶栅气动反命题的新解法(Ⅱ)

本文利用映象平面ξ-η的特点, 将旋成面叶栅气动正命题的微分-积分解法与叶型修正的“喷气模型”结合, 建立了一种适用于旋成面叶栅气动反命题的新解法。叶片表面给定的气动参数分布作为反命题的边界条件, 在求解过程中直接得到满足。通过对一例实际透平叶型的计算验证了本文解法的可靠性和实用性。      

跨声速翼型多目标优化设计方法

本文通过采用多目标分级优化方法并基于求解全速势方程的跨声速粘流翼型计算方法，研究发展了一种多参数、多约束和多目标的跨声速翼型数值优化设计方法，应用该方法可以从普通低速翼型和超临界民办型出发通过多目标优化后得到在跨声速区的多个马赫数下阻力系数最小化的翼型几何外形，设计实践表明，该方法具有收敛快，调用目标函数次数少等优点。     

流-势函数法的人工密度格式在叶型设计中的应用

本文采用流-势函数法进行跨音速叶型设计,用人工密度格式取代类型相关格式。与相关格式相比,人工密度格式具有全场算法统一、收敛过程稳定等优点。本文提供的设计算例表明,该方法能有效地进行叶型的初始工程设计。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。