机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

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

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

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

支撑翼布局客机总体参数对结构重量的影响

针对支撑翼布局客机在总体设计阶段对结构重量分析评估的需求，建立结构重量与支撑翼总体参数之间的定量描述关系。利用有限元分析求解工程梁模型的方法建立考虑静强度和颤振约束的支撑翼结构分析模型，采用全速势方程加边界层修正的方法建立气动载荷计算模型，在此基础上建立支撑翼总体参数与主承力结构尺寸的优化流程，实现支撑翼结构概念设计优化及重量分析。以支撑翼客机为对象的算例研究表明，应用这一优化流程能够分析总体参数对结构重量的影响规律，并实现定量描述，可以为支撑翼布局客机概念方案设计提供技术支持和数据参考。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

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

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

三维复杂无粘流场的自适应八叉树结构直角网格算法模拟

运用自适应八叉树结构直角网格对三维复杂流场进行了划分 ,并采用点面判别法进行外法向量的判断 ,使与之关联的与物面相交网格的体积计算更为简便准确。在此基础上 ,采用以中心差分为基础的 Jameson的有限体积法进行 Euler方程求解 ,在通量计算过程中 ,采用以面为基础的算法 ,减少了计算工作量。对多个外挂、翼身组合体及全机等复杂外形绕流流场进行了模拟 ,均取得了较好的效果     

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

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

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

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

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

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

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

复合材料超临界机翼跨音速颤振实验与数值计算研究

颤振课题是飞机设计过程中常常遇到的一个关键技术问题。以支线客机ARJ21超临界机翼颤振特性研究为背景,在俄罗斯TsAGI的T-106风洞中完成了该复合材料机翼跨音速颤振实验,基于N-S方程和无限插值方法（TFI）生成三维贴体运动网格对超临界机翼跨音速颤振进行了并行计算。结果表明：复合材料的超临界机翼在跨音速区域具有跨音速颤振＂凹坑＂现象;与风洞实验结果相比,有较好的一致性,为使用超临界机翼的运输类飞机跨音速颤振特性预计提供了一定的参考。     

An inverse code based on a residual-correction method for wing design

A fast, efficient and reliable code for wing design has been developed at INTA coupling a residual-correction method by Bauer, McFadden and Garabedian to an inviscid solver for wing analysis in transonic flow. Smoothing procedures, including Bézier cubic splines, are used to avoid irregularities of the wing surface, as well as the twist distribution. A non-conservative full potential finite differences code is used as the flow solver. The redesign of two wings in transonic flow are presented with promising results. The computing time required to achieve convergence up to a good level of accuracy is less than 1 hour in a Workstation of 300 Mflops. The code is actually a powerful tool for preliminary wing design. Future work will be addressed to the addition of viscous effects and to studies on the optimization of target pressure distributions.     

Numerical simulation of viscous flows with adaptively refined Cartesian grid

Quadtree-based Cartesian grid was automatically generated from specified geometry.Adaptive refinements were performed according to geometric parameters and solution of flow field.An altered CCST (curvature corrected symmetry technique) approach was proposed to apply solid wall boundary conditions.Driven flows in a square cavity and flows around NACA0012 airfoil were simulated and compared with the result of published structured grid and stretched Cartesian grid.The results show that solid wall boundary condition are accurately applied by current altered CCST approach,while incompressible/compressible subsonic,transonic and supersonic viscous flows are adequately simulated with adaptively refined Cartesian grid.      

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

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

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

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

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

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

Aerodynamic benefits of adaptive wing technology

Especially in the case of large transonic transport aircraft, flight conditions change considerably during a typical mission. This is accounted for by multiple but fixed design points which compromise aircraft performance. Employing adaptive wing technology where the wing geometry, or other means of flow control, adjusts the flow development to the changing freestream and load conditions allows us to explore fully the flow potential at each point of the flight envelope. Various means of flow control by geometric adaptation and by direct boundary layer control have been investigated within the German national program ADIF and the EU-project EUROSHOCK II and their potential explored. Here, corresponding results are presented and discussed, indicating the applicability and benefits of the adaptive control methods considered. It is also demonstrated that, generally, flow control must be adaptive to work in real aeronautical conditions since these conditions change within the mission flight envelope.