NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

跨音速翼型和机翼的反设计计算方法

给出一种跨音速翼型和机翼的反设计计算方法。对所应用的积分方程反方法引入人工粘性项;采用Riegels因子法消除前缘奇性;对强激波问题采用光滑-松弛过程;并将方程中的系数积分成解析形式;对二维翼型反设计计算还提出了一种封闭形式的正则化条件。算例结果表明,该方法对跨音速翼型和机翼设计是一种有效的工具。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

欧拉方程多段翼型绕流计算的嵌套网格方法

本文采用嵌合体埋入技术，将主段翼型和襟翼分成两子域，对多段翼型的复杂构型进行了分区迭代求解，在有限体体积离散的基础上，求解欧拉方程计算各子域流场的解，通过各流场的耦合迭代到多段翼型绕流的解，典型的带３０％襟翼ＧＡ（Ｗ）－１翼型算例表明，该计算方法稳定，能较地模拟多段翼型中主段翼型与襟翼交接处的复杂流动现象计算结果与实验值有较好的一致性。     

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

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

大展弦比机翼模型设计对翼型流场气动特性的影响

采用SST两方程湍流模型,通过求解非定常Navier-Stokes方程,模拟了大展弦比机翼风洞模型振动条件下的翼型流场,总结了翼型不同振动状况下的流场和气动力特点,分析了模型设计中的不同振动情况对风洞试验结果的影响。研究结果表明：在大展弦比机翼风洞模型的设计中,将翼型的重心设计在机翼的弹性轴之后,对风洞试验的精度较为有利。此结论对大展弦比机翼的风洞实验模型设计有指导意义。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

三维翼跨音速守恒型全位势流高效差分算法的应用

本文对跨音速守恒型全位势流,用Holst的高效有限差分算法,研制出分析三维翼跨音速流的实用程序。与Holst的TWING程序不同,本文用快速保角转绘生成二维正交O型贴体坐标网格,并以此为基础构造三维贴体坐标网格,因而节省网格生成的机时。本文不仅将方法用到有中等后掠的厚翼,还提供大后掠大根梢比薄翼的计算经验。对根梢比很大,前后缘后掠角相差非常大的机翼,指出本算法的困难,并分析了困难原因。