低雷诺数下小型无人机翼型气动特性分析

以小型无人机翼型研究为背景,利用基于线性稳定性理论的eN方法对对小型无人机常用的翼型CLARKY在雷诺数Re=1.0×10~5、5×10~5、1.0×10~6,迎角由-5°~20°时的气动性能进行了计算和对比分析。随着雷诺数的增大,翼型上表面的转捩位置不断向前缘移动,气流分离则由完全分离逐渐转变为层流分离泡结构,使得翼型的最大升力系数和临界迎角增大,阻力减小,最大升阻比显著增大,有利迎角逐渐减小,翼型CLARKY的气动特性逐渐得到改善。     

基于γ-Re_(θt)转捩模型的低雷诺数翼型数值分析

通过综合运用Michel转捩判据和γ-Re_(θt)转捩模型的方法实现对低雷诺数翼型的气动分析。首先通过Michel转捩判据和层流分离位置分别估算转捩动量厚度雷诺数,然后分别用Langtry和Tomac提出的经验关系式计算转捩动量厚度雷诺数和转捩区长度,通过UDF(User-Defined Functions)将不同组的3个经验关联值写入Fluent软件对E387翼型(雷诺数为3×10~5)进行数值分析,并与试验值进行了对比。结果表明:基于Michel转捩判据比基于层流分离更能准确预测转捩位置,基于层流分离预测的转捩位置过于靠前,但是2种方法的气动力计算结果与试验值都比较吻合,Langtry和Tomac的经验关系式计算结果之间差别不大,Tomac方法能够计算出分离泡,且与Selig油流试验吻合较好。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟,研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型,基于结构化拼接网格,对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性,雷诺数对翼型气动特性具有显著影响。随雷诺数降低,时均分离泡外形增大、位置后移,平均阻力系数增大,特别是在较低雷诺数下,翼型升阻力系数随时间出现振荡现象。进一步研究表明,造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低,流动黏性增大,导致分离剪切层速度梯度减小,流动发生转捩及再附位置后移,直至翼型表面不再发生转捩和再附。     

低雷诺数下翼型不同分离流态的大涡模拟

采用高精度大涡模拟方法,对5°来流迎角、马赫数0.4、三个不同雷诺数(55 000、100 000和150 000)的NACA 0025翼型进行仿真,研究低雷诺数条件下翼型的气动特性。通过对比分析3种工况的计算结果,发现翼型绕流存在两种不同的分离流态:Re=55 000和100 000时,翼型上表面出现大尺度的开式分离,形成宽的尾迹区;Re=150 000时,上表面边界层分离后再附到翼型表面,形成时均化的闭式分离泡,尾迹宽度明显减小。无论哪种流态,Kelvin-Helmholtz(K-H)不稳定性均对层流分离诱导转捩过程起重要作用,雷诺数的增加导致转捩过程加速,时均化的分离区也从开式变为闭式。     

低雷诺数下层流分离的等离子体控制

为有效控制层流分离特性,消除或减弱低雷诺数时小迎角下的升力非线性现象,改善翼型升力特性,并通过翼型的上表面转捩带与油流显示测量对等离子体激励控制机理进行阐述,对厚度为16%椭圆翼型低雷诺数下的气动特性进行了风洞试验研究。在此基础上,在上表面前缘10%弦长处布置激励器,通过压力分布测量观察等离子体激励对层流分离的影响。试验结果表明：当翼型上表面仅发生层流分离时,等离子体激励和转捩带的作用类似,可以有效延迟或者消除后缘层流分离,从而增加升力;当翼型上表面出现层流分离气泡并发生再附现象时,等离子体可以有效减小或者消除层流分离泡的范围,从而减小升力;通过控制层流分离,占空循环等离子体激励可以实现对低雷诺数小迎角下的升力的线性控制。     

湍流度对低雷诺数翼型气动特性的影响研究

为研究湍流度对低雷诺数翼型气动特性的影响,采用经过风洞试验验证的基于γ-(Reθt)转捩模型的RANS数值模拟方法,针对典型低雷诺数翼型E387,选取不同雷诺数/湍流度状态开展了对比分析.研究结果表明,湍流度对翼型气动特性的影响十分显著且规律较为复杂.高湍流度可以在一定状态下使升阻比大幅提升,显著改善低雷诺数翼型的气动...     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

低雷诺数下翼型层流分离泡及吹吸气控制数值研究

在低雷诺数下Eppler387翼型表面会出现层流分离泡现象,因此本文使用Fluent求解器开展吹/吸气控制翼型表面层流分离泡的数值研究,主要探究了射流位置、射流角度、射流速度比三个控制参数对层流分离泡控制的影响规律.研究结果表明,采用与SSTκ-ω湍流模型耦合的γ～(Reθt)转捩模型可以准确预测层流分离泡的位置;吹/吸气可以有效抑制低雷诺数下层流分离泡的发展,明显提高低雷诺数下翼型升阻比.固定射流位置,较大吸气速度比和较小吹气速度比可分别获得较好的流动控制效果,且吸气控制比吹气控制对层流分离泡的抑制作用更加有效.     

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

基于IDDES框架的γ-Reθ转捩模型

针对高超声速流动的γ-Re_θ转捩模型在模拟强制转捩时存在捕获边界层内扰动不足的缺点,将RANS（Reynolds Averaged Navier-Stokes）框架改为IDDES（Improved Delayed Detached Eddy Simulation）框架,既能像基于RANS框架的转捩模型一样模拟复杂构型的自然转捩,又能发挥IDDES能够捕获更多脉动信息的优点,较为准确地模拟粗糙颗粒诱导强制转捩。通过对一系列简单构型的自然转捩及来流马赫数为6条件下平板上单个粗糙颗粒诱导强制转捩的模拟表明,模型既能体现γ-Re_θ转捩模型的优点,在自然转捩模拟中具有较强的鲁棒性,能够反映雷诺数等因素对转捩位置的影响规律;也能体现IDDES方法的优点,能够捕捉粗糙颗粒诱导的扰动及涡结构,从而较为准确地刻画出强制转捩的整个流程。     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

基于表面粗糙度的超高负荷低压涡轮叶片附面层控制

实验研究了表面粗糙度对PACKD-A低压涡轮叶型损失及附面层特性的影响.实验分别在定常来流与非定常条件下进行,非定常条件下的上游尾迹通过运动的圆棒来模拟.粗糙叶片通过在光洁叶片表面切槽,埋入砂纸制作成型.叶型损失与吸力面载荷使用气动探针与壁面静压孔结合压差传感器来测量,附面层流场使用热线探针来测量.结果表明:覆盖5.2%吸力面弦长(起始于44.3%吸力面弦长,终止于49.5%吸力面弦长),粗糙高度为8.82μm的控制方案在非定常条件下效果最佳,该方案可在整个考察雷诺数范围内(3×104~12×104)降低叶型损失;覆盖19.5%吸力面弦长(起始于30%吸力面弦长,终止于49.5%吸力面弦长),粗糙高度为20.91μm的控制方案在定常条件下效果最佳,该方案可在低雷诺数范围内(小于8×104)降低叶型损失、扩大叶型正常工作雷诺数范围,但在高雷诺数下(大于8×104),却带来了一定的额外损失.      

The Solar Quiet Chromosphere—Corona Transition Region

The chromosphere-corona transition region of the Sun enjoys both simplicities of character and complexities of character which result from its very thin geometrical extent. The simplicities derive from the reasonably clear view of the energy balance (both observationally and theoretically), while the complexities derive from both the proximity of the not-so-clearly viewed regions below and above, and the almost certain convolutions and perhaps discontinuities in the three dimensional geometry of the transition sheet. While observational resolution and spectral information has improved greatly in recent years, the problems associated with a single vantage point, the Earth and its environs, have not gone away. To understand the transition region we must resolve structures radially and temporally as well as in the plane of the sky. This revised version was published online in June 2006 with corrections to the Cover Date.     

粗糙壁面诱导的流动转捩数值模拟方法

为预测粗糙壁面诱导的流动转捩,对Langtry提出的γ-Re_θ转捩模型增加粗糙壁面的比耗散率和涡黏性的边界条件,并对模型中的经验关联函数——转捩动量厚度雷诺数进行修正,引入等效沙粒表面粗糙度作为变量,使模型具有一定的预测粗糙壁面诱导的流动转捩的能力.对粗糙平板自然转捩算例和变压力梯度平板绕流算例进行数值模拟,计算结果与风洞实验数据符合的较好.主要结论如下:表面粗糙度一般会增加边界层内的传热系数和阻力系数,同时使层流到湍流的转捩位置提前;自然转捩的转捩位置受表面粗糙度影响较大,与光滑壁面相比,平板算例中0.15mm的表面粗糙度使转捩位置提前40%;分离诱导转捩的转捩位置受表面粗糙度影响较弱,随着表面粗糙度逐渐增加,转捩位置和分离泡位置略有后移,且分离泡强度逐渐减弱,分离泡之后的阻力系数增加.      

后掠翼边界层横流不稳定转捩预测模型

通过对经典Falkner-Skan-Cooke三维边界层相似解的理论分析和数值求解,结合二维边界层转捩判据的思想,采用由试验数据标定的C1准则关系式求解横流不稳定转捩位移厚度雷诺数,建立了针对固定前缘后掠角机翼的横流转捩判据,并且通过方程求解和数据拟合得到了该转捩判据的数学结果.应用该模型对30°前缘后掠角的ONERA-M6机翼和45°前缘后掠角的NLF(2)-0415无限展长机翼进行了横流不稳定转捩数值模拟.模拟结果显示:改进后的转捩模型预测所得到的转捩位置精度较高,均与后掠翼横流试验数据吻合较好,从而证明了构建的横流不稳定转捩判据的合理性和实用性.      

间断Galerkin方法计算叶型转捩流动

使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.      

横向矩形微槽对高超边界层失稳的控制作用

模拟了马赫数为6的空间发展平板边界层,通过在平板表面添加二维横向微槽研究了其对基本流及第二模态扰动波的影响。结果表明:所研究的二维微槽构形与常规多孔涂层相比具有较大的尺寸(100 μm以上),微槽对基本流的流向速度影响很小,空间采样点上的频谱并未出现新的不稳定模态。微槽有利于减小摩擦阻力,槽内流动导致的压差阻力比摩擦阻力低1~2个数量级,总阻力随着开槽率的增加而减小,当开槽率为06时,减阻效果超过40%。在入口以两种方式添加了第二模态扰动波,包括单个扰动和多个叠加扰动,下游的演化结果说明横向微槽能够在一个宽带频率范围内对第二模态的增长起到明显的抑制作用,且控制效果随着开槽率的增大而增强。