基于变域变分有限元的翼型反设计

反问题作为一种自由(未知)边界问题近年来得到了广泛的研究。变域变分通过把自由边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而自由边界求解和流场分析可以完全耦合进行。本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型，计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果。文中还采用了两种计算网格(H型和O型)来设计翼型，比较了它们的精度，计算表明本方法完全能用于翼型的反设计。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

跨声速粘性流绕振荡翼型的非定常计算

本文采用时间和穴是二阶精度的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式及ｑ－ω二方程微分模式，结合网格自适应技术，数值模拟Ｎ－Ｓ方程，计算了跨声速下的翼型非定常运动，包括俯仰，浮沉和前后平移振荡。结果表明，压力分布和气动系数与实验基本符合，微分模式和自适应网格能够显著提高激波和边界层计算精度。     

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

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

翼型动失速和非定常载荷的综合翼型数据方法

本文给出受迫谐振翼型的动失速工程估算方法。本方法基于凤洞试验，综合分析翼型动、稳态特性之间的差别与动、稳态条件之间的关系，建立一套估算动、稳态特性之间差别的经验公式，修正稳态特性，得到相应的动态特性。用本方法计算了三种翼型不同动态条件（包括深失速和后掠）的动失速特性，并与测量和文献结果进行了比较，结果符合得相当好。     

基于多项式混沌法的翼型不确定性分析及梯度优化设计

耦合伴随方法和非嵌入式多项式混沌法，发展了高效、可靠的不确定性梯度优化设计方法。利用伴随方程法求解目标函数对不确定性变量的导数，发展了一种梯度增强型多项式混沌法。通过亚声速和跨声速下等多种算例可以证明该方法可以提高不确定性分析的效率和精度。同时，利用基于方差分解的全局敏感性分析方法对不确定性变量的敏感性进行了量化。建立了多项式混沌耦合伴随方程的统计矩梯度求解方法，并结合梯度增强型多项式混沌法搭建不确定性梯度优化设计系统。基于该优化设计系统对二维低亚声速和跨声速翼型开展确定性及不确定性优化设计研究。优化结果显示，相比于确定性优化设计，不确定性优化设计通过合理权衡确定性性能和不确定性性能，可提高抵抗马赫数和迎角不确定性扰动的能力，同时优化性能均值和标准差。其中阻力系数均值最大可降低17%，阻力系数标准差最大可降低80%。而确定性优化设计可能导致性能鲁棒性的降低。     

Lorentz力的大小和作用位置对翼型绕流及升力和阻力的影响

分布在弱电介质溶液中的电磁力(Lorentz力),可以有效地控制边界层的流动,且Lorentz力的大小、方向以及作用位置可根据实际需要来调整。在转动水槽中,通过示踪粒子流场显示以及升力和阻力的测试,对Lorentz力对翼型绕流以及升力和阻力的影响进行了实验研究,并利用拟压缩方法数值求解贴体坐标中的控制方程,对实验中的相应问题进行数值研究。基于实验和计算结果,对电磁力控制时的Lorentz力的大小以及作用位置对翼型绕流和升、阻力的影响及其内在原因进行了讨论。     

基于实用参数控制的翼型参数化方法

在翼型的优化设计中,通常要对二元翼型进行参数化方法表示,其目标是利用较少的参数得到足够大的设计空间。分析了Hicks-Henne型函数叠加法和PARSEC几何参数法的特点和不足,提出了利用厚度、弯度等实用参数进行翼型定义的方法,该方法可以通过厚度和弯度参数更加直观地理解翼型轮廓和翼型基本气动性能。通过分析3种方法对5种翼型的最小二乘拟合结果得出结论:实用参数法的误差分布比较平均,且能够比另外两种方法更精确地表示翼型。     

高阶面元法计算多段翼型之位流

本文给出了一种求解多段翼型位流的高阶面元法。面元和奇点的分布分别采用了二阶和三阶的形式。由于高阶面元法较为精确地反映了绕流体的几何形状和奇点表面分布,本文给出了比低阶面元法更为精确的解,并且消除了低阶面元法计算中经常出现的翼型后缘处解的波动。