低速大迎角三角翼旋涡流动的有粘与无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Euler方程和雷诺平均Nayierstokes方程，应用Beam－Warming近似因子分解格式及其对角化形式，采用BaldwinLomax代数湍流模型。计算模型为70°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Euler方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

大迎角体涡,翼涡干扰的Euler方程数值模拟

本文用Ｅｕｌｅｒ方程的简化区算法计算翼－身组合体绕流。当大迎角时，计算能自动捕捉到大后掠前缘的分离涡，机射的分离涡由强加Ｋｕｔｔａ条件生成，形成体涡与翼涡的相互干拢。本文对加长前机身的ＡＧＡＲＤ－Ｂ翼－射组合体模型，给出机射有分离和无分离时，横流截面的速度矢量分布、法向力和压力中心，分析体涡的影响。     

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

翼身组合体复杂涡系相互干扰的数值模拟

本文应用有限体积法求解Ｅｕｌｅｒ方程，数值模拟翼身组合大攻角、跨声速绕流。成功地模拟了机身体涡、机翼脱体涡、翼尖侧缘涡的形成、发展及其相互间的干扰。     

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

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

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

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

考虑湍流影响的前缘分离涡运动与破碎特性的研究

本文以Ｎ－Ｓ方程和Ｈａｌｌ涡核模型假设为基础，导出了描述湍流的涡核运动方程。利用差分计算方法，对三角翼前缘分离涡运动及其破碎特性进行了数值求解，分析了涡核流场的结构和各物理量的变化特性，反映了三角翼前缘分离涡运动及其破碎特性的实质，说明了湍流涡核方程能更有地模拟三角翼前缘分离涡运动和准确地确定涡核破碎位置，并得出了攻钐雷诺数对破碎位置的影响，计算结果与实验数据十分吻合。     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

跨音速机翼结构动力响应的计算研究

通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Ｅｕｌｅｒ方程,并用Ｊａｍｅｓｏｎ提出的有限体积方案,四步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进求解;结构模型采用二自由度ＰＡＰ模型。经分析,计算结果合理。     

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

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

三维水下爆炸气泡的数值模拟

在理想不可压势流假定下,采用高阶有限元离散边界积分方程求解水流场。假定泡内气体遵循等熵关系。用Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方法模拟气泡的演化。计算结果与球泡“精确解”及实验值比较,表明数学模型合理且三维数值方法精度高。计算给出气泡迁移、射流等重要动力学行为和物体受力分析     

可压缩涡环和平面激波相互作用的数值研究

本文从非定常的Ｅｕｌｅｒ方程出发,进行了可压缩涡环和平面激波相互作用的数值研究。     

跨声速翼型有粘／无粘强干扰计算

本文介绍一种二元跨声速激波-边界层强干扰的计算方法。边界层计算采用湍流边界层积分反方法,它借助Whitfield和Swafford提供的既适合附着流,也适合分离流的速度剖面表达式。跨声速无粘流用全速势方程模拟。通过边界上排溢速度来考虑粘性的影响,用有粘/无粘迭代得到粘性流解。本方法计算的结果与其它方法以及实验的结果进行了比较,证明该方法可以在工程上推广使用。     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件,采用数值求解Ｅｕｌｅｒ方程的方法模拟模型在风洞中的绕流场,然后将洞壁边界条件用远场边界条件替换,用同样方法算出模型在自由流中的绕流场,从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性,采用重叠网格技术,对飞机大迎角标模和ＳＢ－０３模型进行了亚声速大迎角的洞壁干扰计算与修正。