多孔表面抑制第二模态失稳的最优开孔率和孔半径分析

多孔介质、微槽道和超声波吸声材料等可用于抑制高超声速边界层第二模态扰动波。通过声波在无穷长小管道中的传播模型给定多孔介质表面的扰动边界条件,采用时间模式的线性稳定性理论分析了多孔介质表面对边界层稳定性的影响。对马赫数6热壁平板边界层的考察表明,多孔介质表面不仅可以大大推迟第二模态扰动波的中性点,还可以大大抑制最不稳定第二模态扰动波的幅值增长率。为了找到最优控制参数,在较广的参数范围内考察了多孔介质表面的开孔率和孔半径对第二模态扰动波幅值增长率的影响,得到了最优开孔率和孔半径沿流向的分布。最后,还考察了基本流的当地非平行性(法向速度)对控制效果和最优控制参数的影响。     

热化学非平衡高超声速平板边界层线性稳定性分析

高超声速边界层流动转捩是近期空气动力学研究的热点问题。对于环境扰动较小的自然转捩过程,稳定性分析已被证明是研究扰动演化的重要手段。另一方面,高超声速边界层内的温度会随着马赫数的升高而快速上升,极高的温度会引起所谓的高温真实气体效应,使得量热完全气体假设失效,从而对边界层稳定性和转捩产生影响。本文针对高温热化学非平衡气体,利用空气5组分模型开展了平板边界层的线性稳定性分析,重点研究了热化学过程对模态稳定性的影响,并探究了边界层离散谱模态的演化和同步过程。研究表明,对于由第二模态主导的高超声速二维边界层:(1)扰动相比基本流更趋向于热化学冻结态;(2)扰动方程中新出现的非平衡源项的扰动项对稳定性影响很小,非平衡过程主要是通过改变基本流剖面来间接影响稳定性;(3)声速是影响第二及更高模态的重要参数,热化学平衡态假设引起的声速计算式的变化能够解释边界层温度和厚度降低时第二模态频率反而降低的非常规趋势。     

超声波吸收层在天地往返运载器热防护中的应用与研究

超声波吸收层能够抑制高超声速边界层中的第二模态波,延迟边界层转捩,从而降低天地往返运载器热防护系统的重量、成本和复杂度。为了研究超声波吸收层的吸波特性,建立了线性稳定性理论模型。首先,基于四阶龙格-库塔法,采用打靶方式,获得了高超声速边界层基本流。然后,将超声波吸收层的导纳引入稳定性方程边界条件,利用切比雪夫谱方法,求解时间模式下第二模态波的增长率。对于槽道构型超声波吸收层,模型预测结果与已有数据进行了对比,二者符合良好。在几何参数相同的情况下,圆孔构型超声波吸收层对应的第二模态波增长率整体上低于槽道构型。当微腔深度为0.3,微腔数量为6,两种构型超声波吸收层对应的第二模态波增长率具有最小值。第二模态波增长率的最小值出现在微腔数量最少的超声波吸收层上,这是由第二模态波在超声波吸收层顶部散射产生了新的高频扰动波所造成。     

利用横向槽改善气膜冷却效率的实验

为了获得横向槽对气膜冷却效率的影响规律,分别对气膜孔出口全开槽、上游开槽和不开槽三种气膜冷却结构形式的平板实验模型进行了实验研究,获得了带横向槽气膜冷却结构的冷却效率随吹风比的变化规律和开槽形式对冷却效率的影响规律,以揭示横向槽改善气膜冷却效率的机理.研究表明:在横向槽的作用下,冷却射流向主流的垂直扩散受到抑制,二次流在流出横向槽后能更好的贴覆绝热壁面,提高了气膜冷却效率,并且在气膜孔出口处开槽,对提高气膜孔下游的平均冷却效率有较大的作用,上游开槽结构的下游平均冷却效率高于全开槽和不开槽结构的平均冷却效率.      

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

多孔壁面对高速边界层最优增长条带二次失稳的影响规律

边界层转捩是高超声速飞行器设计中的关键基础理论问题。当环境扰动强度较高时，将在模态扰动失稳区上游发生由最优增长条带二次失稳触发的亚临界转捩。为评估多孔壁面在亚临界转捩中的控制效果，以超/高超声速平板边界层流动为研究对象，建立了基于伴随抛物化稳定性方程的优化系统与求解方法。以最优扰动非线性演化形成的三维条带边界层为新的基本流动开展全局稳定性分析，研究表明：多孔壁面对第一模态频率范围内的二次失稳扰动为促进作用，对第二模态频率范围内的二次失稳扰动起抑制作用，并且转折频率接近局部快/慢模态的同步频率，对于工程应用中多孔涂层的布置方案具有一定的指导意义。     

吹吸扰动对壁湍流边界层摩擦阻力的影响

对风洞中平板湍流边界层施加局部周期性吹吸扰动,并用IFA-300热线风速仪测量了其近壁区域不同法向位置瞬时速度的时间序列信号.利用湍流边界层对数律平均速度剖面,计算湍流的壁面摩擦阻力.对平板湍流脉动速度信号用子波分析进行多尺度分解,通过条件采样以及相位平均的方法获得不同尺度下的相干结构波形和能量分布,来研究吹吸扰动对平板湍流边界层的摩擦阻力的影响.结果表明,扰动使得窄缝附近流动阻力增加,而下游区域流动阻力减小.      

激波对气膜干扰作用的抑制机理研究

为了研究激波对气膜冷却效果的破坏机理，并消除这种影响，以平板壁面为基础设计了一种带有卸压槽的壁面结构。通过数值计算研究了主流马赫数为3.2，冷流马赫数分别为1.0、0.6和0.4三种工况下开槽壁面对激波破坏的抑制作用。结果表明，在有激波入射的条件下，开槽壁面比平板壁面具有更好的流场结构，可使激波导致的近壁气膜的分离区最多减小至原来的三分之一，并有效减弱气膜入射后在肩部产生的反向涡旋对，这很好地抑制了气膜的卷吸和与主流的掺混。计算显示开槽壁面最大能够使壁面冷却效率提高6%，且这种作用效果与通过卸压槽的气流流量大小有关。此外仿真结果表明，在相同条件下波前卸压较波后卸压效果更好。通过合理安排卸压槽位置及槽面宽度，可以将总压损失控制在合理范围内。     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

超声速边界层转捩实验可控扰动的引入方法

通过数值模拟的方法,研究了在二维平板超声速边界层中,由壁面处安装的激振薄膜引入的扰动的演化过程.主要得到了如下结论:由激振薄膜引入的扰动能很快演化成标准的T-S波,但扰动的各物理量的演化为T-S波的速度并不相同,即各物理量并不同时形成T-S波,速度最先形成,而温度最晚.在相同的条件下,由不同激振形式或不同宽度的薄膜引入的扰动的演化速度基本相同.     

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

对原始的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),却带来了一定的额外损失.      

固定转捩在改善振荡来流下低雷诺数翼型气动性能中的应用

低雷诺数下,翼型表面易于出现的分离转捩现象会降低翼型的气动性能,采用数值计算方法探讨了固定转捩在改善低雷诺数下翼型气动性能中的应用。定常来流下,雷诺数为4×104时,在翼型E387上表面自由转捩位置之前的一定位置处使翼型固定转捩,计算结果表明分离区减小了,升力系数和升阻比明显提高。振荡来流下,雷诺数在8×104的50%范围内变化时,分别计算了NACA4412在0°迎角和E387在3°迎角时自由转捩及在60%弦长处固定转捩的升阻特性,结果显示自由转捩的翼型当雷诺数减小时,性能急剧恶化,而采用固定转捩的翼型受其影响要小得多,具有更加稳定的气动性能。     

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.     

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

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

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

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

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

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

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

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