计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

跨声速伽辽金有限元法所用上风技术的改进

伽辽金有限元法是解跨声速绕流的有效方法,对超临界流动,为了自动捕捉到激波,通常采用了上风技术,但导致了气流通过激波时质量守恒条件得不到满足。本文通过Bateman变分原理和加罚泛函找到了强制满足质量守恒约束条件的附加项,大大改善了用伽辽金有限元法解超临界绕流的计算结果。     

翼身组合体的地面效应

本文用涡格镜像法计算小迎角时翼身组合体的地面效应。在机身上分布离散的马蹄涡系(涡格)。在机身内部分布机翼的马蹄涡系对机身(圆)的镜像涡系。在机身表面分布空间源(汇)。若机身迎角不为零时,还要沿机身轴线分布二维偶极子。用这种翼身组合体的对地镜像来模拟地面的作用。本文主要研究小迎角的情况。故翼身组合体及其镜像均与地面平行。在TQ—16电子计算机上,用ALGOL—60语言编成程序,计算了翼身组合体在接近于地面时纵向空气动力特性,其结果为:升力曲线斜率随高弦比的减少而增加、低头力矩随之增加,压力中心也随着变化。     

航天器高超声速气动力的数值计算

本文把航天器的外形用数以百计的平面四边形元素来近似。利用高超声速绕流时面积元素之间的干扰可忽略的假设,分别计算出各面积元素的气动力,积分后得出航天器的气动力,算例的结果是令人满意的。本文也描述了航天器表面网格的生成技术和计算机绘制航天器在不同视角下外形图的原理。本文方法可供航天器初步气动设计时用。     

翼-身组合体跨音速绕流全位势积分方程数值计算

 用积分方程方法求解 Prandtl- Glauert算子表示的全位势方程 ,并计算了翼 -身组合体跨音速绕流。用 Murman- Cole差分格式计算空间场源强度 ,以捕捉激波。计算结果与相应的实验结果符合良好     

亚音速侧滑薄翼气动特性的数值计算方法

本文利用变换坐标的方法,把处理来流平行于机翼对称面的涡格法推广到有侧滑的情况去.文中计算了亚音速流中侧滑薄翼的气动特性。机翼可以是任意平面形状和具有上反角。文中给出了一些算例,并把其结果与其他方法和实验数据比较,结果是令人满意的。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

计算亚音速定常流中偏转副翼时机翼—机身组合体气动特性的涡格镜像法

本文用涡格镜像数值解法计算了亚音速定常流中偏转副翼(或襟翼)时翼身组合体的气动特性。方法适用于简单副翼和差动副翼。压缩性的影响本文不采用P—G转换,而是直接使用可压缩空间后掠马蹄涡和源(汇)。计算方案已编成Algol—60语言源程序,计算了2个型号偏副翼时气动特性并和实验结果进行比较。