翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

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

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

民机翼身整流影响分析与设计研究

针对某民机翼身组合体进行了翼身整流设计,以减弱机翼机身附面层之间的相互干扰,消除翼身结合处气流的分离,改善该民机的阻力性能。将翼身整流的多种设计方案(不同类型,不同安装位置,不同的几何尺寸等)配置到某民机翼身组合体模型上,用CFD方法对其流场进行粘性绕流数值模拟并对结果进行优化分析,得到了针对该民机减阻效果较好的翼身整流设计方案。     

带襟、副翼三维机翼粘性绕流计算

运用嵌套网格方法严格考虑机翼切口和副翼、襟翼端面几何形状生成了合理的计算网格 ,并应用雷诺平均 Navier-Stokes方程和 Johnson-King湍流模型计算了带副翼三维机翼和带双襟翼三维机翼的粘性绕流。计算实践表明 ,用该方法可以很好地预计压力分布和最大升力。     

翼-身组合体跨音速守恒型全位势流AF2高效算法的应用

 本文用Joukowski交换将圆截面机身的每一横截面逐一转绘为垂直切缝,因而将翼-身组合体转绘为某单独机翼。我们用作者采用过的数值保角转绘法对转换的单独机翼生成O型贴体坐标网格,然后用逆转绘构造翼-身组合体绕流的O型贴体坐标网格。跨音速计算采用守恒型金位势方程,精确边界条件和AF2高效有限差分迭代算法。本文指出,只需少量修改,就可将单独翼分析程序扩展为可同时适用于翼-身组合体。本文考虑的组合体其机身为有限或无限长,圆形截面,机翼前后机身轴线允许弯曲;机翼有任意平面形状,但翼尖弦长不能太短。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

用全位势方程计算侧滑机翼的跨声速绕流特性

本文对大后掠小展弦比细长机翼，由垂直于机翼纵轴的每一横流截面生成Ｏ型网格，从而形成对机翼流场的Ｈ－Ｏ型网格，用守恒型全位势方程，差分方法和隐式近似因式分解ＡＦ２迭代算法计算侧滑机翼的流场。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

基于非结构重叠网格的二维外挂物投放模拟

本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。     

亚声速全机纵横向气动特性计算方法

本文介绍了一种涡格方法,它适用于亚声速飞机的纵横向气动特性计算。该方法中,升力面用布置在中弧面的平均平面上的单涡格层,或用布置在上、下表面的平均平面上的双涡格层模拟,机身则由布置在表面上的四边形涡环或由布置在一组同轴棱柱面上的涡栅来模拟。地面效应用布置镜像涡系的办法处理。文中讨论了机翼尾涡偏斜对平尾、立尾的影响。 该方法计算的结果与有关的理论和实验进行了比较。这些比较表明,此方法计算精度满足工程设计的要求。