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

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

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

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

Godunov方法在复杂流场数值模拟中的应用

本文根据Ｖａｎ Ｌｅｅｒ提出的方法给出了求解三维非定常Ｅｕｌｅｒ方程和三维定常Ｅｕｌｅｒ方程的二阶Ｇｏｄｕｎｏｖ方法。方法要点为：在各体积离散单元中,用流动变量的分段线性分布代替一阶方法中的分段常量分布,在确定线性分布的斜率中引入单调性限制条件;在离散单元界面构造Ｒｉｅｍａｎｎ问题,由Ｒｉｅｍａｎｎ间断解确定单元侧面的通量;用预测－修正两步法对离散方程推进求解,用上述方法结合单区网格或多区重迭网格     

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

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

用全位势方程计算细长翼-身组合体的可压缩绕流

用Ｈ－Ｏ网格、守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕大后掠细长翼－身组合体的可压缩流动。对ＡＧＡＲＤ－Ｂ翼－身组合体的计算表明,计算的自由流Ｍａｃｈ数可从亚音速直到低超音速。算出的升力与力矩特性与实验数据符合较好     

多部件组合体分块网格生成技术及应用

给出多部件组合体分块网格生在的一种方案及实践,包括交界面网格生成技术、空间网格生成方法及质量控制。以一个平尾和机翼具有一定高度差的机身－机翼－平尾－立尾组合体外形为例,生成与公共交界面或边界面正交或接近正交的分块网格。用于Ｅｕｌｅｒ方程跨声速来流几个迎角的计算,结果与实验符合较好。     

超声速前机身及翼—身组合体的全速势方程数值计算

通过数值求解全速热方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某些一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

翼-身-尾组合体的网格生成及 Euler 方程数值模拟

采用代数方法快速生成翼-身-尾组合体三维O-C型贴体网格。该代数网格生成方法使用双边界法,可保证物面附近近似正交,真实地模拟了机身头部及与机翼、尾翼交界处的形状,提高了代数方法生成网格的质量。并用中心格式有限体积法求解跨声速Euler方程,以某翼-身-尾组合体模型为例,计算结果与实验符合良好。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

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

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

全机复杂流场网格生成及N-S方程数值模拟

全机复杂流场网格生成对计算结果起着举足轻重的作用。本文首先通过合理的拓扑结构,分块对接的方法对复杂流场进行了网格划分,然后采用N-S方程有限体积法及分区对接法对翼—身—尾全机流场进行数值模拟。计算结果验证了网格生成的有效性和计算结果的合理性。     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

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

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

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

通过人工压缩性方法求解了三维定常不可压Ｅｕｌｅｒ方程和雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因子分解格式及其对角化形式，采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。计算模型为７０°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Ｅｕｌｅｒ方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

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

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

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

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

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

大展弦比飞机翼身Euler绕流计算

给出了用Ｅｕｌｅｒ方程对大展弦比翼身组合体绕流的数值模拟，将Ｋｒｏｌｌ的求解欧拉方程的有限体积法，推广到求解机翼机身流动问题，基本思想是使用保角变换、代数变换和椭圆方程变换混合方法生成机翼机身贴体网格。     

超声速前机身及翼－身组合体的全速势方程数值计算

通过数值求解全速势方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。