在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

多段翼型失速分离特性的计算

本文用高阶奇点分布板块法(Panel Method》及边界层理论,对多段翼型的失速性能进行计算。采用分离尾迹的模型,按尾迹是流线用迭代法求出分离尾迹边界(此为内迭代)。在尾迹上配置奇点,然后把尾迹考虑在内进行粘/位流迭代(此为外迭代)。收敛后可得多段翼型空气动力特性,且可计算到失速以后的情况。算例有单段翼型GA(W)-1,NACA 4412,和F111飞机的三段翼剖面,所得结果与实验值能很好地符合。     

基于Euler方程的直升机旋翼翼型反设计方法

直升机旋翼流场的求解和贴体网格的生成是进行旋翼翼型反设计的关键,本文采用Euler方程求解直升机旋翼翼型流场,并通过求解泊松方程生成围绕翼型的网格.在此基础上,结合线化二维亚音速和超音速流的小扰动理论,采用MGM方程作为反设计方程,基于有限差分方法建立一套直升机旋翼翼型反设计的迭代方法.在给定的目标压力分布条件下,与CFD方法相结合,反设计得出满足要求的翼型,并与目标翼型吻合良好,表明本文建立的翼型反设计方法能够有效地用于直升机旋翼翼型设计.     

翼上螺旋桨构型耦合气动特性及翼型优化设计

针对半环形式翼上螺旋桨构型，研究了螺旋桨-机翼耦合流场特性，并以短距起降（STOL）状态最优升阻特性为目标对机翼翼型进行全局优化。首先，针对螺旋桨-气动面耦合构型，通过动量源法与真实桨叶模型CFD的计算对比，分析动量源法用于该构型设计分析的可行性。其次，为得到有利于桨-翼耦合特征的新翼型，建立了翼上螺旋桨构型自由型面变形（FFD）参数化模型，采用遗传算法对翼上螺旋桨构型机翼翼型进行全局寻优设计，分析了优化翼型参数及流场变化规律。最后，将优化翼型用于三维半环形机翼，分析其流场特性与二维计算结果的异同，验证二维翼型优化的有效性。结果表明：真实桨叶多重参考系（MRF）方法不能准确计算翼上螺旋桨构型下的流场结构，而动量源法计算结果与真实桨叶滑移网格非定常方法较为吻合；采用二维动量源CFD方法进行翼型的遗传算法优化是有效的，受半涵道的保护，二维优化翼型的优势在三维构型中得到了有效继承；翼上螺旋桨构型的翼型优化应当着重关注翼面曲率变化，在本文计算状态下，通过增加桨盘附近翼面曲率、保持附着流动来加强Coanda效应，有效实现了气动增升，优化后机翼升力提高了22.51%，显著减弱桨盘后高压区并产生二次吸力峰值，同时保持了机翼负阻力特性。     

基于特征正交分解的跨声速流场重构和翼型反设计方法研究

在二维流场的重构问题中应用特征正交分解(Proper Orthogonal Decomposition,POD)数据处理方法。利用CFD技术计算得到的流场快照对气动力模型进行降阶,然后利用基于POD的降阶模型(Reduced Order Model,ROM)对所需的流场参数进行重构,在快照范围内可以得到高精度的结果,且具有一定的外插能力。在翼型反设计问题中该方法仍然是成功的,通过修正快照向量,利用基于POD降阶模型的数据重构方法,由已知的翼型表面压力分布通过反设计就能够高效精确地得到对应的最优翼型形状,这极大地简化了翼型反设计问题。本文分别在跨声速范围对RAE 2822翼型的流场重构和Korn翼型及NACA 63212翼型的反设计进行了验证,证明了基于POD的流场重构和翼型反设计方法的有效性和高效性。     

翼型失速的非线性动力学特性及其控制

以NACA0012翼型为研究对象,应用CBS( Characteristic-Based Split)有限元方法对翼型绕流流场进行数值模拟,着重对翼型失速现象的动力学机理和特征进行研究.首先给出了静止翼型在不同迎角下流场的瞬时流线图以及静态翼型的迎角-升力系数关系曲线,并且分析出此曲线具有滞后、突跳等典型的非线性动力学...     

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

具有自由流线边界二元叶栅的通用解法

本文提出一种满足自由流线边界条件二元叶栅的新解法——两个相关流场交叉迭代的计算方法,统一地论述了普通叶栅的上下游区、分流叶片叶栅和串列叶栅的二元亚音速流场计算问题,计算程序可以通用于这三种情况。 在分流叶片叶栅和串列叶栅的计算中提出了一种流场最小干扰的物理条件,从而可以计算出由分流叶片或分流面所形成的两个叶道之间的流量分配比。 介绍了数值计算结果,并与试验数据作了比较,这些计算与比较初步验证了本文方法,以及所建议的计算流量分配比物理条件的合理性,并对分流叶片叶轮的气动设计提出了一些有益的设想。     

利用流线法决定绕扁长椭球流动分离的位置

本文按照流线法把所测量的扁长椭球的压力分布用于决定涡层型三维流动分离,按边界层边缘处流线会聚的位置近似地决定分离线的位置。作为算例,计算了扁长椭球体在α=30°状态下的绕流,并在水洞中作了流场显示,数值结果和实验结果的对比表明,二者十分一致。分析表明,涡层型分离不仅与壁面极限流线的收敛有关,而且还与边界层边缘处流线会聚有关,前者的出现并不是引起涡层型流动分离的充要条件。     

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

基于天线安装的菱形翼无人机翼型优化设计

以菱形翼布局高空长航时(HALE)传感器无人机研究为背景,分析了翼型内部安装平面相控阵雷达天线的基本原理。通过推导雷达探测性能的估算方程,建立了在翼型中安装平面相控阵雷达天线的数学模型。采用数值模拟的方法,对典型雷诺数下的RAE2822翼型的气动特性进行了研究并与试验数据进行了对比,验证了使用有限体积法k-ωSST(Shear Stress Transport)湍流模型求解雷诺平均Navier-Stokes方程在这一状态下的可靠性和适用性;同时对不同厚度的NACA 64A系列翼型的流场结构和流动机理进行了分析,在此基础上提出了一种满足平面相控阵雷达天线安装情况下翼型的优化设计思想。优化结果验证了这一设计思想的可行性。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

超声速流场的流线-特征线坐标变换与应用

为了求解二维平面弯曲激波波后流场,讨论并发展了一系列基于流线-特征线坐标系变换的流场代数计算方法。该系列方法依据二维平面流场中气流角、静压沿特征线近似不变的特点,可快速求解平面弯曲激波波后流场。其中,等气流角近似法适合模拟流线轨迹,等静压近似法适合求解波后流场参数。在此基础上,又提出了一种改进的等气流角-等静压混合方法用于计算弯曲激波波后流场。和特征线法对比,等气流角-等静压混合法计算得到的流场基本特征与特征线法得到的结果相同,在来流马赫数Ma=6和Ma=4情况下误差分别仅为0.5%和0.15%,证实了该方法在求解二维平面弯曲激波流场中的适用性。     

多段翼型增升小翼纳维-斯托克斯方程分析

装有扁平板增升小翼的多段翼型流场数值分析已经完成。研究的小翼装在接近主翼段后缘，高度范围是0．25％～1．25％参考翼型弦长。二维数值模拟在嵌入式结构拓扑网格上采用不可压缩纳维-斯托克斯（N-S）方程求解分析。新的网格细化提高了网格边界重叠部位的计算精度。在固定雷诺数条件下，对装有窄缝襟翼的二段翼型研究了各种小翼的影响。计算和实验结果都表明，相对未装小翼的基本翼型，主翼段凹陷处装有小翼提高了最大升力和升阻比。计算得出的流线分布揭示，正是小翼的存在引起了襟翼流场分离的减弱。也研究了三段翼型在一个雷诺数范围内的流场情况。对优化的襟翼来说，计算和实验所得的雷诺数影响结果一致。当襟翼从最有利的位置移动时，数值分析表明，小翼有助于翼型重新优化，大概能起到所优化襟翼效益的1％量级的作用。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

一种新的边界处理方法在笛卡尔网格中的应用

采用自适应笛卡尔网格方法求解Euler方程,固壁边界条件通过一种离散强迫浸入边界方法引入,边界切割网格与流场网格保持一致,从而有效解决了小切割网格单元的时间步长限制问题.针对自适应笛卡尔网格的特点提出了新的虚拟点反射方法和插值方法,解决了数值吸力峰值问题.对RAE2822翼型、双NACA0012翼型和Suddhoo三元翼型的流动状况进行了无黏数值模拟,并与现有的理论解、单域以及分区结构网格解进行了对比.结果表明:该离散强迫浸入边界方法有效地引进了固壁边界条件,结合自适应笛卡尔网格技术能够准确模拟复杂几何外形下的流动状况.      

旋翼桨叶翼型流场和跨音速翼涡干扰噪声研究

基于桨涡干扰流动的噪声场分析了NACA0012翼型的远场气动噪声.采用Euler方程计算了NACA0012、NACA23012、OA207三种翼型的绕流场,选取Kirchhoff面选为包围整个翼型的一层网格,通过得到NACA0012翼型的Kirchhoff面上在桨尖涡干扰流动下的近场解,根据Kirchhoff方法求解远场声压.样例计算分析体现出:在计算桨涡干扰噪声时,旋翼CFD技术与Kirchhoff方法相结合的良好效果.     

翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。