线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

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

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

机翼跨声速绕流的压强极小积分有限元算法及应用

本文从全速位方程出发。利用压强极小积分,给出了机翼跨声速绕流的有限元算法,并编制出相应的FORTRAN语言程序。文中给出了适于有限元法的网格生成技术,采用了线松弛迭代法和逐步加密法进行数值求解,提供了M6机翼和一个三角翼的计算结果。     

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

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

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

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

民机的一种新型布局形式——翼身融合体飞机

 讨论了一种区别于机翼与圆柱机身的传统民机外形——翼身融合体（BWB）外形。从给定容积下减少浸润面积和摩擦阻力的考虑发展出了翼身融合体的概念。在满足系列设计约束条件下讨论了适于大型（800座）客机的翼身融合体的技术和商业上的可能性。概念设计研究表明,翼身融合体外形可获得比常规外形更好的性能。巡航升阻比可从传统外形的19提高至23。480座的BWB外形与同量级的传统外形（A380-700）相比,最大起飞重量可减少18%,每座位燃油消耗可减少32%,具有良好的环保性。通过气动载荷分布规律的反设计,翼型修型和三维优化设计等气动设计方法可以有效减阻和提高BWB的气动性能。由于BWB的高度融合性,讨论了多学科优化设计方法的重要性和有效性。BWB概念能够发展成BWB系列机,进一步提高巡航马赫数也是可能的。     

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

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

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

控制超临界翼型边界层分离的微型涡流发生器数值模拟

本文基于任意曲线坐标系下的雷诺平均Navier-Stokes方程和Spalart-Allmaras一方程湍流模型,采用对接拼接网格技术和多重网格加速收敛技术,对安装有叶片式微型涡流发生器的超临界机翼翼身组合体进行了数值模拟,研究了微型涡流发生器的高度和弦向安装位置对超临界机翼附面层流动控制的机理以及对超临界机翼气动性能的影响规律。     

动力增升对翼身融合飞机沉降特性的改善

翼身融合体飞机起降阶段的“沉降问题”是影响该类布局飞行安全的典型问题,为此提出了一种将动力风扇竖直嵌入机翼中的动力增升方案。分析了动力增升系统对飞行的气动影响,提出了动力增升系统的匹配设计和评价方法,探讨了动力增升系统布局的主要约束。在综合考虑动力增升系统推质比、俯仰力矩配平以及操稳特性等限制因素基础上,基于示例翼身融合体无人机,匹配设计了动力增升系统及其布局方案,针对过渡飞行状态设计了姿态稳定控制律。飞行试验表明:使用动力增升能缓解进近阶段升降舵操纵所引起的下沉率突增与高度沉降现象,同时可将验证机的起、降速度降低20%以上,且该方案易于在小型翼身融合体飞机上实现。      

一种新的变前掠翼无人机气动布局

研究了一种新的变前掠翼无人机气动布局概念,在低、亚、跨和超声速状态下可通过改变机翼的前掠角来获取最佳的气动性能。根据设计指标和翼身融合技术初步设计了其几何外形,并采用三维Navier-Stokes方程数值模拟和对比分析了5种不同任务构型的气动特性。结果表明：①在Ma=0.6巡航时,平直翼加挂副油箱构型最大升阻比为14.55,而三角翼构型仅为8.29;②在Ma=0.4机动时,45°前掠翼构型失速迎角达到38°且具有最大的升力系数2.455,较平直翼构型提高了4.9%;③在Ma=1.5高速突防时,三角翼零升阻力系数最小,比平直翼加挂副油箱构型减小了14.4%,最大升阻比提高了34.6%;④所有计算状态下俯仰力矩特性均良好。研究结果验证了变前掠翼无人机气动布局新概念的合理性和先进性,可为高性能无人机的设计提供参考。     

变前掠翼气动布局纵向性能研究

发展和改进了一种变前掠翼气动布局及其设计思想,设计了集战斗和轰炸于一体的通用任务模式,并根据设计指标和翼身融合技术初步设计了变翼几何外形。采用N-S方程计算了该几何外形不同构型的纵向气动性能,分析了机翼前掠角对不同任务模式下纵向气动性能的影响。结果表明,变前掠翼气动布局可通过改变机翼前掠角来获取最佳纵向气动性能。     

A fast and automatic full-potential finite volume solver on Cartesian grids for unconventional configurations

To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algo rithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual (GMRES) algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution con verging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body (BWB).The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.     

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

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

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

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

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

无尾布局纵向操纵的嵌入式舵面概念研究

基于提高W型无尾翼身融合布局大迎角俯仰操纵能力的设想,探索一种机身下表面嵌入式俯仰操纵舵面的新概念。采用CFD方法,研究纵向两个位置、两种弦向长度组成的四种机身下表面嵌入式舵面构型对无尾气动布局性能的影响及其流动机理,分析了嵌入式舵面提供纵向操纵控制的作用原理。研究结果表明,下表面嵌入式舵面可作为无尾布局大迎角俯仰操纵的有效补充措施,为无尾布局大迎角俯仰操纵提供了新的技术途径,具有进一步深入研究的价值。     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

单级入轨运载器气动布局研究

对单级入轨运载器气动布局进行了研究,设计了一种翼身融合体气动布局,并进行了布局优化。计算结果表明该布局在Ma=6的时候,最大升阻比Kmax可以达到3.755,从高速到低速都能维持在4左右。通过与已有乘波体方案的对比,表明新设计的翼身融合体是一种很有潜力的单级入轨运载器气动布局可选方案。