跨音速翼型和机翼的反设计计算方法

给出一种跨音速翼型和机翼的反设计计算方法。对所应用的积分方程反方法引入人工粘性项;采用Riegels因子法消除前缘奇性;对强激波问题采用光滑-松弛过程;并将方程中的系数积分成解析形式;对二维翼型反设计计算还提出了一种封闭形式的正则化条件。算例结果表明,该方法对跨音速翼型和机翼设计是一种有效的工具。     

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

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

楔—后掠柱激波干扰流场和热流的工程算法

本文针对航天飞机机翼前缘的简化模型,即楔-后掠柱模型,在大后掠、高超声速条件下,提供了一种计算柱体前缘激波干扰流场和热流的工程计算方法。在中国科技大学流体力学实验室的激波风洞中取得了激波干扰的流场显示结果。计算的结果与国内外实验数据的比较表明,这是一种有效的方法。     

超声速和高超声速机翼俯仰导数—当地流活塞理论解法

本文给出计算超声速和高超声速尖前缘三维机翼俯仰导数的一个近似解析方法。该方法是根据当地流活塞理论和片条理论导得的。在激波附着于机翼前缘情况下,本方法可用于计算任意迎角,颇为一般的平面形状、尖前缘剖面机翼的俯仰导数。中小平均迎角的数值结果与有限的试验数据和其它理论作了比较,示例表明本方法的结果令人满意。     

典型低音爆构型的近场音爆计算研究

音爆的预测是超声速民用飞机设计的关键技术。采用多块结构化网格有限体积法求解Euler/N-S控制方程,对SeebALR旋成体和三角翼两个典型低音爆构型国际标准算例的近场音爆进行了计算,与风洞试验数据进行了对比分析。结果表明,对于上述标准算例的简单绕流场进行近场音爆计算时,流场黏性效应不是必须考虑的因素;随着方位角的变化,三角翼尾部激波从单个激波演变为多激波;来流攻角的变化不影响头部激波强度,机翼前缘激波强度随攻角的增加而逐渐增强。     

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

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

一种近无激波机翼修型的工程设计方法

本文利用现有解全速势方程的一种分析方法和“虚拟气体”概念,提出一种近无激波机翼修型的工程设计方法。避免原“虚拟气体”设计方法中在超音速域内的推进,而采用部分穿透速度边界条件的办法来获得修型的效果。以实现接近于无激波机翼的绕流流场。算例结果表明这种工程设计方法使用灵活,效果良好。     

前缘剖面形状对80°／60°后掠双三角翼上涡流运动影响的数值模拟研究

采用上风高分辨率格式———通量差分分裂格式离散求解三维可压非定常薄层N S方程 ,数值模拟前缘剖面形状对 80°/60°后掠双三角翼上涡流运动的影响。计算模型包括尖前缘、菱形前缘、圆形前缘等三种不同前缘形状的机翼。计算结果表明 :三种不同剖面形状的前缘可以诱导产生不同的前缘分离 ,形成的各前缘剪切层的特点也不同。尖前缘机翼的边条翼涡和外翼涡合并点最靠前 ,其次为菱形前缘和圆形前缘的机翼。在大迎角情况下 ,三种机翼上的内翼涡发生破裂。圆形前缘机翼上的内翼涡破裂点比其它两种情况下破裂点的位置较靠前。涡的合并对二次涡的结构和特点也有显著的影响。     

超声速边界层内小激波的数值模拟

本文通过有限差分计算求解了NS方程,对超声速二维平板边界层的非线性稳定性问题进行了模拟.流场的来流马赫数为4.5.平板从尖前缘开始,流场在前缘处产生一道前缘激波.数值计算模拟了边界层对成倍频关系的双频扰动波的吸收问题.通过计算发现,在一定频率下边界层内可出现小激波.计算给出了小激波的发展过程.     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

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

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

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

跨声速钝体绕流数值模拟研究

通过数值求解ＮＳ方程，进行了三方面工作（１）格式和算法对返回舱跨声速流场特性的影响研究；（２）网格形态对捕捉声速弱激波的有效性研究；（３）模拟了无支杆和有支杆返回舱跨声大攻角绕流流场，给出了后球体上的跨声速激波和后支杆上的第二道激波。     

跨音速三维机翼及翼身组合体非守恒型全位势流AF3算法

 本文对跨音速非守恒型全位势流,用Baker的快速有限差分法,改变其网络拓扑,研制出分析机翼及翼身组合体跨音速绕流的实用程序。其中包括:对“C-H”型贴体坐标网格提出AF3算法的具体实施方案取得了快速收敛效果;采用张量形任意曲线坐标上的非守恒全位势方程,使得程序中跨音速求解部分既能用于分析单独机翼又能用于翼身组合体跨音速绕流计算。     

AF-2格式的加速收敛参数研究

一、引言 Ballhaus等人提出的求解跨音速定常流动的隐式近似因式分解法（AF-2）具有收敛快速,占计算机存储少等优点。文献[1,2]分别对于完全小扰动方程和全位势方程的计算给出了很好的例证。文献[3]表明上述结论对于横向小扰动位势方程也是成立的。AF-2格式的快速收敛主要是通过加速收敛参数的循环取值来实现的。因此,合理正确地选取加速收敛参数序列对AF-2格式至关重要。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

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

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