用虚拟区域法结合虚拟边界条件求解对称及管壁流动问题

分析了用基于Lagrange乘子的虚拟区域法数值求解时,翼型后缘存在的奇异问题.为了保持直角坐标网格结构简单及可利用快速算法的特点,提出引入虚拟边界条件控制虚拟流动的方法处理奇异问题.用虚拟区域法数值模拟了翼型不可压对称绕流;并把该方法推广应用到管壁流动的数值模拟,给出了喷管内流场的算例.数值结果表明,虚拟边界条件的引入有助于控制Lagrange乘子的变化,从而改善数值模拟的精度.     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

介质阻挡放电等离子体对翼型流动分离控制的实验研究

在低速开口风洞中进行了等离子体激励器对NACA0015翼型流动分离控制的实验研究.采用PIV技术,对翼型绕流流场进行了测量,显示了施加等离子体激励后流场的变化.通过五分量天平对升力和阻力的测量,研究了激励电压和激励频率对翼型流动分离控制的规律.研究表明,低风速下在翼型前缘施加等离子体激励,能够有效地控制翼型流动分离,在来流为20m/s时,最大升力系数增加11%,失速迎角增加6°;在给定的流动状态下,激励电压和激励频率存在一个阈值,不同迎角下该阈值不同,迎角越大,分离越严重,对激励强度的要求也越高.     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围.在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术.通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式.研究表明:在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值.     

翼型动态失速气动力二次峰值数值模拟研究

为揭示翼型动态失速状态下气动力二次峰值的发生机理,基于运动嵌套网格技术、有限体积方法、LU?SGS隐式格式和Roe?MUSCL格式建立了俯仰振荡翼型非定常流场的数值模拟方法.首先,基于所建立的数值方法对NACA0012翼型在深度动态失速状态下的气动特性进行模拟,计算结果与试验数据吻合良好,验证了数值模拟方法的准确性.然...     

用于弹性旋翼流场模拟的网格变形方法研究

建立了一个基于改进"Ball-Vertex"的具有拓扑无关性的网格变形方法,用来精确模拟桨叶发生弹性变形后贴体网格的运动以及旋翼近场气动力环境,为更好地分析弹性变形对桨叶贴体网格的影响提供了一套解决方案.该网格变形方法从基本的弹簧模拟思路出发,用一网状弹簧系统来模拟桨叶贴体网格的变形和运动.在网格结点周围引入"Ball-Vertex"概念,通过添加相关的冗余约束,该方法保证网格结点在变形过程中不会越过与该点所对应的对角面,从而避免了畸形网格单元的产生,可有效地用于较大变形情况下的流场数值计算.使用所建立的方法,首先对绕二维翼型的结构网格和非结构网格进行了网格变形模拟,以表明方法的有效性;然后,将该方法应用于振荡翼型的流场数值求解,进一步验证了所建立的方法;最后,结合该网格变形方法和旋翼流场求解方法对弹性旋翼流场进行了数值计算.在此基础上,得出了一些有意义的结论.     

直升机旋翼翼型及桨叶气动外形反设计分析

建立了一个基于余量修正思想的直升机旋翼反设计计算方法,用于直升机桨叶的气动外形设计研究.使用Poisson方程为控制方程生成围绕桨叶的贴体网格;在悬停状态下建立了以Euler方程为主控方程的旋翼流场求解方法,并采用了嵌套网格方法进行数值计算;在流场计算及网格生成基础上,采用MGM方程作为翼型反设计方程,建立了一套直升机旋翼翼型及桨叶气动外形的反设计方法.应用该方法,分别对二维翼型以及悬停状态下的旋翼桨叶进行了反设计分析.反设计结果表明,在给定的目标压力分布条件下,使用本文方法分别获得了满足要求的二维翼型及直升机桨叶外形,并与目标压力吻合良好.     

旋翼跨音速非定常黏性绕流的高效CFD模拟方法

为提高旋翼跨音速黏性绕流CFD模拟的效率,建立了一套基于隐式LU-SGS算法和OpenMP并行策略的旋翼非定常流场高效数值求解方法。首先,基于二维剖面翼型的Poisson方程求解和网格插值、翻折方法生成绕桨叶的O-C-O型贴体正交网格,并采用高效的"扰动衍射"(Disturbance diffraction method,DDM)挖洞方法、并行化″Inverse map″(IM)的贡献单元搜索方法,构建了旋翼运动嵌套网格系统。在此基础上,以耦合S-A湍流模型的非定常RANS方程为主控方程,对流通量采用高精度的Roe-MUSCL格式进行离散,时间推进采用隐式LU-SGS双时间方法,同时采用OpenMP并行策略加速旋翼流场求解。最后,运用所建立的方法分别对悬停ONERA 7A旋翼、前飞Caradonna-Tung(C-T)旋翼无升力状态及SA349/2旋翼有升力状态的气动特性及涡尾迹特征进行了模拟,通过与试验值对比验证了文中方法在旋翼非定常流场CFD模拟中的有效性和高效性。     

一种绕曲面物体的可压缩湍流混合格式(英文)

提出了一种中心-迎风型混合格式。在该混合格式中,中心差分格式和Roe通量差分裂格式进行混合,它们之间的切换通过一个二进制开关函数实现。为了验证该混合格式在计算绕曲面物体可压缩湍流问题时的可靠性,尤其是带激波的流动问题,采用分离涡模拟方法计算了3个典型的问题。研究结果表明,当前数值结果与已有的实验数据相符较好,这说明该混合格式可以用来研究带激波和湍流的曲面物体可压缩绕流问题。     

矩形渡槽槽体结构绕流流场数值模拟研究

基于二维非定常不可压缩雷诺平均N-S方程,采用RNG(renormalization group)k-ε两方程湍流模型,对矩形渡槽在空槽和满槽情况下的二维绕流特性进行了分析.计算网格为四边形结构网格,采用有限体积法对微分方程进行离散,应用SIMPLE算法解决压强-速度耦合关系.引入空腔流动研究成果对空槽和满槽风载体形系数基本相同的原因进行了理论分析,揭示了空槽时槽体内复杂的旋涡流动结构.     

极低湍流度下翼型边界层及近场尾流的雷诺应力测量试验研究

在西北工业大学低湍流度风洞中，用热线风速仪在对称翼型边界层及近场尾流中，于２０多处剖面上移测速度型、湍流强度、雷诺正应力及切应力等湍流参量。着重讨论了各雷诺应力分量的分布。结果表明，在边界屋和尾流宽度方向及流动方向上均有明显的规律性，其有关特征值明显低于常规风洞湍流度下的结果。     

应用后缘装置的跨声速层流翼型多岛优化设计

高雷诺数状态下,自然层流技术(Natural laminar flow,NLF)是减小机翼表面湍流摩擦阻力的有效方法。然而由于层流翼面上大范围顺压梯度的存在使得后缘处的恢复压差更大,产生更强的激波。因此在减小摩擦阻力的同时又增加了激波阻力。本文采用后缘装置(Trailing edge device,TED)来控制翼型后缘处的激波强度,基于线性稳定性理论(Linear stability theory,LST)的eN方法对流动进行转捩判断,进而应用多岛并行多目标进化算法(multi-objective evolutionary algorithm,MOEA)以获得大范围层流区域和弱化激波强度为目标对翼型进行优化设计。优化结果表明合作均衡策略耦合进化算法可以快速地捕捉到该多目标问题的Pareto阵面解,阵面上翼型的波阻力和摩擦阻力性能较初始翼型大大改善。同时,采用后缘装置控制激波强度时,无论在设计点还是偏离设计点时,优化后翼型均具有良好的升阻力特性和鲁棒性。     

基于γ-Reθ转捩模型的高超声速复杂构型转捩模拟

在大规模并行可压缩Navier-Stokes求解器AHL3D框架上，搭建了γ-Re_θ转捩模型。在该模型基础上，通过高超声速二维进气道构型算例，加入了可压缩性修正，使其能够模拟可压缩性对转捩位置的影响，同时通过修改分离诱导转捩关键参数，增加了模型对粗糙颗粒诱导强制转捩的敏感性，最后对耦合关系式也进行了适当的修改。为了验证修改后的模型对高超声速自然转捩和强制转捩的预测能力，对Ma7.4 Ames全尺寸模型和单个粗糙颗粒诱导Ma6的平板转捩进行了模拟。结果表明:与原始模型相比，修改后的模型的转捩位置被大大推迟，并且粗糙颗粒诱导转捩的作用被加强，与实验结果吻合良好。采用此模型对X-51A的20%缩比进气道模型在普渡大学Ma6静音风洞中的试验状态进行了模拟，模型不仅能够反映来流湍流度对转捩的影响，也能反映转捩带对转捩的促进作用。结果显示修改的转捩模型在高超声速复杂构型的转捩预测及研究中具有很好的应用潜力。     

基于合成射流的旋翼翼型动态失速控制研究

针对直升机旋翼工作环境下来流速度和迎角(Angle of attack,AoA)耦合引起的动态失速问题,建立了基于合成射流的旋翼动态失速控制的数值模拟方法。采用运动嵌套网格方法,通过对翼型的平移和旋转实现变来流速度-变迎角的耦合。以积分形式的雷诺平均N-S方程为主控方程,空间离散使用Roe格式,时间离散为隐式LU-SGS方法,以OA209翼型为研究对象,在翼型上表面放置合成射流激振器,开展了射流位置、动量系数、无量纲频率以及偏角等参数对轻度失速、深度失速下翼型动态失速控制的研究。研究发现,轻度失速下,射流位置靠近气流分离点时(20%c附近,c为翼型弦长),对逆压梯度引起的轻度失速控制效果最佳。深度失速下气流分离点虽在5%c之前,但射流位于前缘分离泡后端(10%c附近)时控制效果较好。大迎角需要较大的动量系数才能有效控制。射流频率对涡结构的尺寸和数量会产生一定影响,能改变气动特性波动幅度。较小的射流偏角对轻度失速的控制更有效,而深度失速则需要较大的偏角。     

Body-Fitted Momentum Source Method for Predicting Rotor Aerodynamics Characteristics in Hover

To gain high efficiency for the simulation of the aerodynamic characteristics of the rotor in hover,body-fitted momentum source(BFMS)method is proposed. In this method,the actual blade geometry is represented by the single layer of volume grid surrounding the blade. Aiming at correctly simulating the aerodynamic characteristics of the discrete cells along the chordwise of blade airfoil section,a new distributed force model is proposed. For comparison,the RANS method with S-A turbulence model and the steady rotor momentum source(SRMS)method based on embedded grid systems are established,respectively. And the grid connecting methodology is improved to embed the blade into the background grids for the three methods. Then,simulations are performed for the hovering Caradonna-Tung rotor by these methods,and the calculated results are compared with the available experimental data. Moreover,the pressure distributions along the blade are compared with the conventional momentum source methods. It is demonstrated that the BFMS method can be employed as an effective approach to predict rotor aerodynamic characteristics with a low computational resource and reasonable accuracy.     

地面效应下的不同翼型亚声速气动特性分析

采用数值方法研究了亚声速地面效应条件下不同翼型的气动特性,进一步以Ma=0.5来流工况为例,研究了翼型参数和飞行高度对气动特性的影响。计算结果表明在Ma为0.5、迎角为6°的地效情况下,翼型弯度减小,更容易在翼型前缘产生激波阻力;翼型下翼面后缘弯度增大使得后缘压力更高,升力系数和低头力矩相应增大;随着飞行高度的减小,地效作用加强,翼型下翼面压力增大,下翼面的升力增量大于上翼面吸力损失,机翼升力系数和升阻比增加越来越显著。     

低雷诺数效应对0.6m连续式风洞性能影响

风洞试验中通常采用降低运行总压的方法来扩大风洞雷诺数模拟的下边界。同常压试验相比，在低雷诺数条件下，风洞流场是否存在明显变化，风洞流场品质是否满足指标要求，直接影响风洞试验数据的精准度。为了研究低雷诺数效应对0.6m连续式跨声速风洞性能（包括轴流式压缩机性能、总压及马赫数控制精度、流场均匀性）的影响，调试人员在0.6m连续式跨声速风洞中开展了大量相关试验，本文在对试验结果进行整理、分析的基础上，给出了Re对风洞性能的影响。结果表明:（1）Re对压缩机性能、总压控制精度、马赫数控制精度、流场均匀性都有明显影响，当Re_c<5×105（c=0.1√A）时，雷诺数效应明显，且Re越小，影响越大。（2）0.6m连续式跨声速风洞能够真实、准确反映Re对风洞流场性能和测力试验数据的影响规律，是开展高空低雷诺数飞行器、翼型、发动机等性能研究的理想地面模拟试验平台。     

后退式微型后缘装置对翼型气动特性影响的实验研究

介绍了装有后退式微型后缘装置(Rearward Mini-TED)的 NACA23012翼型在低雷诺数条件下的表面压力分布、气动力和 PIV 速度场的风洞实验结果,并与 NACA23012原型翼的对应测量结果进行了对比分析,以探讨 Mini-TED 装置对翼型流场、气动特性产生的影响。本实验风速为15m/s,以弦长为特征量的雷诺数为 Re ≈1.3×105,翼型表面压力分布采用测压孔和压力传感器测量,通过积分获得翼型升力和压差阻力,并利用尾耙测量翼型受到的总阻力。结果表明,后退式 Mini-TED 翼型改变了翼型周围的流场速度分布和尾流流动结构,导致上翼面吸力和下翼面的压力升高,使翼型升力增加,但压差阻力也增加。同时发现后退式 Mini-TED 翼型使前驻点位置后移,加快了上翼面的流动速度,后缘分离受到抑制。     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。