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

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

机翼平面形状参数摄动面元法

在内部Ｄｉｒｉｃｈｉｅｔ边界条件的基础上,藉助于对奇异积分的解析微分,可用解析法计算偶强密度对机翼平面形状乡数的敏感性偏导致（建立一般敏感性方程所需）。机翼平面形状参数摄动后其偶强分布即可藉线性外推得出。从而,压力分布、升力和俯仰力矩系数即可快速确定。振动面元法与相应的低阶面无法计算结果表明,该方法有良好的准确度且使计算工作量大大减少（摄动外推所需ＣＰＵ时间比相应的低阶面元法少两个量级）。     

计算Cauchy和Mangler积分主值的新方法

对于Cauchy积分主值曾提出过一些数值积分法,主要是抽去奇异性方法和内插型或Gauss型方法。对于前者,要求积分核在积分区间上有连续的一阶导数;对于后者,实际上只能求出节点(某类函数零点)处的积分主值,且当积分核不是多项式时     

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

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

亚音速升力面气动敏感性导数计算

具有任意曲线前缘的亚音速升力面的气动敏感性导数由核函数法给出。用自适应积分法计算弦向积分,用Multhopp法结合抽去奇点,计算Mangler积分主值。将积分核展成Chebyshef多项式的渐近展开式以保证结果的收敛性。最后将广义力系数及其敏感性导数表示成简单形式,对椭圆、矩形和后掠机翼作了计算,所得结果在升力面理论精度范围内与直接由核函数法得到的结果一致;而且所得到的偏导数可在飞机设计中分析综合用于多学科优化。