叶尖射流对风力机叶尖流场影响的数值研究

为了设计出更加适合非并网系统的风力机,采用在叶尖加入射流的方法来改变叶尖流场分布.在风力机叶片顶端沿弦长布置3个喷口,采用CFD数值模拟方法,通过改变风力机转速获得原型和带喷口的风力机模型的气动特性以及流场分布.发现在转速低于1 200r/min时,带有安装在不同位置的喷口的风力机功率增长率几乎都为零,射流在这一转速范围内对风力机的气动性能几乎没有影响.而转速高于1 200r/min时,随着转速的增大,喷口位于叶尖中部的风力机的功率增长率快速地增大,射流影响了75%以上叶高的表面的压力分布,在大转速下吸力面低压区范围较大,其叶尖涡涡量低于其他方案中,并且在下游扩散得比其他方案快,改善了风力机下游流场,提高了风力机效率.喷口布置在叶尖前缘时其叶尖涡的局部涡量较原型叶片稍大,降低了风力机功率的输出.喷口布置在叶尖尾缘时基本和原型叶片相同.该结论为设计适用于非并网系统的定桨距变转速风力机提供了基础.     

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

高负荷低压涡轮边界层转捩的实验研究

为了研究高负荷低压涡轮边界层的非定常转捩过程及雷诺数对尾迹边界层相互作用机制的影响,采用热线和表面热膜测试技术在高负荷叶栅实验台上对吸力面边界层的非定常时空演化进行了实验测试。结果表明:高负荷叶型的吸力面边界层在没有来流尾迹作用的情况下会不可避免地发生分离。分离剪切层中的K-H无粘不稳定机制主导了失稳转捩过程;来流尾迹的逆射流会与边界层相互作用,产生卷起涡结构,从而促进分离剪切层的转捩过程。边界层的分离被转捩产生的湍流条带及其后的寂静区抑制;随着雷诺数降低,分离点向上游移动11.4%吸力面弧长,分离泡的流向和法向范围扩大,尾迹引起的卷起涡结构尺度也随之变大。雷诺数通过改变分离剪切层的状态对转捩机制产生影响。     

亚声速等膨胀平面射流的初始流场结构

基于大涡模拟方法与5阶精度混合TCD/WENO格式,对等膨胀射流在马赫数为Ma=0.6条件下的初始流场结构进行了数值模拟。数值结果描述了射流等膨胀过程中,射流初始流场的形成、变化及发展过程。并发现了射流初始流场中涡核为射流失稳的起始点,研究了射流混合层与涡核剪切层上涡的产生、发展与破碎以及它们卷吸环境气体并与射流气体混合的整个过程,揭示了射流初始过程中涡核与涡环直径的变化趋势,得到了射流轴向速度与压力的分布。数值结果可为相关发动机的设计提供重要参考。     

吸力面上气膜冷却对涡轮叶栅流场影响的实验研究

利用气动探针测量和墨迹显示方法,对不同实验方案下,带吸力面气膜冷却的某型涡轮导向器叶栅流场结构进行了实验研究。结果表明,冷气射流与燃气主流的掺混以及卵型涡的形成,使得吸力面根部出现了与通道涡旋向相反的涡系;卵型涡始终以一定形式存在于叶片表面,直到叶栅出口与尾迹相互作用后才达到均匀状态;冷气射流很难进入到通道涡分离线与端壁所形成的三角形区域中,通道涡分离线明显向端壁方向下移。      

平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

间隙高度对自发射流抑制叶尖泄漏的影响

通过数值求解三维定常黏性雷诺时均N-S方程,获得了单孔叶尖自发射流条件下不同叶顶间隙的叶栅流场,对比分析了间隙高度对自发射流与叶尖泄漏流相互作用特性、叶尖泄漏流量以及叶片载荷的影响.结果表明:当叶顶间隙高度为1mm(t/H=0.5%)时,自发射流对泄漏流有明显的阻挡作用,泄漏流量比减少0.06%,同时叶片载荷增加1.39%.当叶顶间隙高度增大到4mm(t/H=2%)时,自发射流的阻挡作用及对叶片载荷的增加作用基本消失;减小间隙高度可以有效提高自发射流的控制效果,同时降低因分离造成的流动损失;自发射流的存在显著改变了间隙流场分布及叶尖吸力面附近静压系数分布,计算发现当泄漏流绕自发射流流过时,下游流场出现类似卡门涡街的涡分布现象.      

连续旋转爆轰发动机冷流场的混合特性研究

连续旋转爆轰发动机（CRDE）中燃料和氧化剂的快速掺混是实现爆轰波成功起爆和稳定传播的重要前提,然而目前国际上关于这方面的研究还相对较少。本文采用大涡模拟（LES）方法,对非预混CRDE中燃料和氧化剂的混合过程及其主要机理开展深入研究。研究结果表明,非预混CRDE流场中存在欠膨胀特征、大尺度涡结构,以及回流区等复杂的流动现象,其中由于Kelvin-Helmholtz（K-H）不稳定性产生的大尺度湍流涡结构是促进氢/氧混合的主要机制。此外,本文还考察了氧气喷注位置对非预混CRDE的流场结构和混合特征的影响,发现氧气喷注位置会影响射流剪切层形态、涡尺度,以及回流区分布等,进而影响氢气和氧气射流的混合过程和混合程度。与其他进气位置相比,氧气在靠近内壁面喷注时更有利于氢/氧的快速掺混。     

直孔射流对压气机叶栅流场影响的实验研究

为了探究直孔射流对压气机叶栅的影响,通过实验方法,结合流场显示技术和流场测试技术,对无控叶栅和直孔射流方案下的压气机平面叶栅在正攻角下的流场结构和气动性能进行了分析。结果表明:无控叶栅中吸力面存在三个螺旋点,而不同射流方案下螺旋点的数量和位置变化明显;无控叶栅端壁存在一个从吸力面起始的分离区,布置射流孔后,在射流孔前发展出马蹄涡,马蹄涡的两个分支的发展情况及其对流场影响随着不同射流方案呈现出不同的特点;射流孔的位置对控制效果有明显的影响,最佳方案减小了3.2%的总压损失,增大了1.86%的通流流量;在最佳方案下,吸力面螺旋点数量减少到了1个,端壁没有明显的尾迹出现,出口处高损失区的欠偏转和端区的过偏转均有所减弱。     

双孔射流对圆柱表面水滴撞击特性影响的数值分析

对于低导热系数复合材料制作的发动机进气部件防冰问题,冲击-热气膜复合防冰技术应用前景较好,但是射流对水滴撞击影响的还缺少研究。本文以三维开双孔圆柱为对象,采用基于欧拉-欧拉框架的三维水滴撞击数值模拟方法,分析了双孔射流对水滴撞击特性的影响。研究发现,由于开孔射流与主流之间的强耦合,射流对来流水滴存在吹袭和扩散等不同作用机制,影响规律与射流强度和物面位置密切相关。针对典型水滴直径(20μm)和来流速度(30m/s)的工况,取沿展向六个截面及孔间中心截面的水滴运动轨迹和涡量场进行了分析。在小吹风比(0.5和1.0)时,射流影响仅限于展向近孔截面,孔间水滴撞击系数相比于无气膜工况最大降幅为7.51%。吹风比为1.5时,射流吹袭和扩散作用更强烈,孔间水滴沿展向被吹向外侧,造成近孔截面水滴量减少而远孔截面增加,孔间截面水滴撞击系数比无气膜工况最大降幅为8.95%,水滴撞击系数由展向外侧向孔间的最大降幅约17.65%。研究初步揭示了射流气膜对圆柱表面来流水滴撞击特性的影响规律和形成原因,为后续冲击-热气膜复合防冰系统的设计提供参考。     

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

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

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

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

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

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

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.     

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

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

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

实验研究了表面粗糙度对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),却带来了一定的额外损失.      

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

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

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

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

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

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

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

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