燃烧与流场在线测量诊断方法研究进展

介绍了上海理工大学颗粒与两相流测量研究所近年来在燃烧和流场在线测量诊断技术方面的研究进展,包括：改进基于光纤光谱仪的火焰辐射温度测量技术,提高火焰温度的测量精度与范围并反演出火焰辐射率及液体燃料组分;辐射与吸收光谱相结合测量脉冲燃烧火焰温度和水蒸气浓度;基于激光诱导荧光法测量喷雾液滴温度;强烟尘水滴干扰条件下基于差分吸收光谱法（DOAS）测量排放污染气体（SO2和 NOx ）浓度;基于饱和蒸汽原理的光谱法研制汞连续排放监测标准系统;光谱法与图像法结合测量高温物体温度场;基于单帧单曝光图像法测量多相流速度场与粒度分布;基于可调谐半导体激光吸收光谱法（TDLAS）同步测量液态纯水水膜蒸发过程中液态参数（水膜厚度及温度）与气态参数（水膜上方水蒸气的温度）等。     

气膜孔形状对排孔下游换热的影响

实验研究了气膜孔几何形状、二次流雷诺数及吹风比对孔排下游的局部换热系数的影响。所用孔形是簸箕形孔、圆锥形孔及圆柱形孔;实验的参数范围为二次流孔径雷诺数Re=10000～25000,二次流吹风比M=0.3～2.0,在上述范围选取了26个工况分别对3种孔形进行了实验。      

矩形横截面螺旋管湍流流动与传热特性的数值研究

 采用了剪切应力输运(SST)k-ω两方程湍流模型并考虑近壁低雷诺数的修正对矩形横截面螺旋管内冷却水流动和传热特性进行了数值研究。数值分析了在不同入口雷诺数、曲率半径以及扭距条件下,螺旋管内的温度、速度场以及流线的变化,讨论了螺旋管内、外壁面对流传热系数的差异及产生机理,同时与直通管道传热性能进行了比较。研究发现由于离心力的作用,螺旋管内存在显著的二次流动,管内、外侧壁面对流传热存在差异。旋转一周后,螺旋管即进入了流动稳定状态,入口雷诺数可以显著提升螺旋管整体的对流换热效率,扭矩和曲率对内外壁面传热效果的影响不大,而窄高型的横截面构型可以显著改善螺旋管的传热效果。研究结果对应用矩形横截面螺旋管的冷却设计提供参考。     

叶片前缘气膜冷却效率的实验研究

采用放大的半圆柱状表面模拟涡轮叶片前缘的形状,对叶片前缘多排圆柱形孔的气膜冷却效率进行了实验研究。测出了孔排区及其下游的局部冷却效率。研究了主流雷诺数及平均吹风比对冷却效率及二次流(冷气)出流轨迹的影响。实验参数范围是:主流雷诺数Re=42000～127000,平均吹风比M=0.8～2.0,测量分8个工况进行。      

壁面约束下椭球颗粒曳力特性研究

颗粒两相流系统广泛存在于自然界和工程应用中,实际颗粒往往并非球形且受其附近壁面的影响。基于欧拉-拉格朗日方法的颗粒两相流数值模拟已成为非球形颗粒两相流系统设计的重要工具,然而其准确性受限于非球形颗粒曳力模型。因此,深入研究壁面约束下非球形颗粒曳力特性并发展近壁修正预测模型具有重要意义。本文采用直接数值模拟方法对有限雷诺数（0.5≤Re≤10）下,不同颗粒长轴取向与来流之间的夹角以及颗粒-壁面间距的单椭球颗粒曳力特性开展详细数值研究,分析比较不同工况下曳力系数的变化,探明雷诺数、攻角以及壁面间距对椭球颗粒曳力系数的影响规律,并评估现有自由流中椭球颗粒曳力关系式的适用性,最后发展了考虑近壁修正的椭球颗粒曳力预测模型,为优化非球形颗粒两相流数值模拟工具提供重要支撑。     

双参数传热实验的液晶瞬态测量不确定度分析

以气膜冷却为例,利用多元隐函数求导的方法推导了双参数传热实验中采用液晶进行瞬态测量时的双参数不确定度计算关系式,采用不确定度传递系数研究了双参数不确定度与温度测量不确定度的关系,得到了不确定度传递系数随主流及二次流无量纲温度的变化曲线,分析了主流及二次流无量纲温度参数对双参数测量不确定度的影响.结果表明,提高主流温度或者增大两次测量中二次流的温差可以减小实验双参数测量不确定度.      

涡轮叶栅端壁表面换热系数的测量

在大尺寸低速开式叶栅传热风洞中对一种高压涡轮导向叶片构成的直叶栅通道端壁表面的换热特性进行了试验研究。在雷诺数为 750 0 0～ 30 0 0 0 0的范围内,采用热膜法测量了端壁换热系数的分布,研究了Re数对叶栅端壁换热的影响规律;同时给出了端壁平均 St数与来流 Re数之间的准则关系式     

来流湍流度及端壁效应对涡轮叶片上对流换热的影响

利用大尺寸低速开式叶栅风洞对涡轮叶片表面的对流换热进行了实验测量,叶片高度方向布置了3个测量位置,距端壁距离分别为5mm,30mm及150mm.重点对比研究了高来流湍流度、端壁效应及来流雷诺数对涡叶片中部及根部区的过渡起始点、过渡区长度及换热的影响。实验参数范围是:来流湍流度Tu=0.75～13.5%,来流雷诺数Re=60000～240000.      

流体喉部推力调节特性实验

采用空气与水作为二次流工质,进行流体喉部的冷流实验,研究了固体火箭发动机流体喉部的推力调节特性.分析了不同二次流工质、注射方式,注射流量下的推力响应时间、扼流性能、推力偏角和推力效率.实验结果表明:注射液态二次流推力响应时间更短;扼流性能、推力偏角与二次流的注射位置及注射角度有关,且随流量比的增大而增大;相同的流量比下,气态二次流的推力性能要比液态二次流的效果更好,但提供相同的流量比,液态二次流需要压比更小,且流量比的调节范围更大.      

凝结相变对低温风洞雷诺数试验的影响

为了研究凝结相变对低温风洞雷诺数试验能力的影响,基于Fluent软件建立了气-液两相凝结流动模型,相变模型采用考虑非等温效应修正的经典成核理论和Gyarmathy液滴生长理论;针对不同来流压力和不同试验雷诺数工况条件,对氮气绕流NACA 0012翼型进行了数值模拟。模拟结果表明:随着来流温度的降低,翼型附近区域气体膨胀越过气态饱和线并达到过冷状态,进一步降低来流温度则会在宏观层面上观测到凝结相变对当地流场的改变,与无凝结相变的流场相比,释放的潜热加热气流导致马赫数降低及压力系数的改变;在确保不破坏翼型气动性能试验的前提下,充分利用气体过冷区域来降低来流压力以此减少驱动功率和液氮喷入量是可行的,或者保持来流压力不变提升低温风洞的试验雷诺数。     

流量突变条件下煤油流动与对流传热特性实验

对流量突变条件下RP-3航空煤油的圆管流动与传热过程进行了实验研究。燃料流动与传热实验是在超临界压力为3MPa、燃料温度为300～650K、雷诺数为3000～60000范围内进行的。实验同步测量了圆管煤油温度、壁温、油压的沿程分布,并应用基于动量与能量守恒的非定常控制体分析方法,获得了航空煤油壁面摩阻系数与传热努塞尔数的分布。实验结果表明:当煤油温度未超过临界值（约650K）、处于过压液态时,该实验条件下的入口流量的突增或者突减过程对摩阻系数与努塞尔数的影响很小。      

液体火箭发动机推力室再生冷却流动与传热计算研究

为了研究液体火箭发动机推力室再生冷却流动与传热的快速仿真方法,建立了推力室再生冷却的准二维模型,对航天飞机主发动机开展了再生冷却流动与传热计算仿真研究,对比分析了再生冷却准二维模型和三维模型的仿真计算结果。研究表明,两种计算模型均可较好地预测推力室燃气及再生冷却剂的流动和传热。三维模型计算精度高,但计算用时较长。计算得到的航天飞机主发动机的燃气侧壁面最高热流密度为162.2MW/m2,最高壁温为1159.7K,冷却剂温升为244.0K,压降为8.5MPa。准二维模型计算结果精度略有降低,但计算时间较三维模型减小了90%。四个参数与三维模型计算结果的差异分别为0.3%,4.4%,8.6%和4.5%,在可接受范围内。本文的准二维模型计算时间短,适用于液体火箭发动机再生冷却结构的方案筛选和优化设计,三维模型计算精度高,适用于设计完成后的校核计算。     

转子叶片径向受限的“冲击-气膜出流”冷却结构流动与换热

为揭示转子叶片径向受限的"冲击-气膜出流"冷却结构流动换热规律,以某型双层壁叶片肋化分割形成的冷却单元为研究对象,通过数值模拟的方式,对冲击雷诺数Rej,旋转数Ro,无因次温比(Tw-Tf)/Tw等参数变化下流场和换热特性变化规律展开研究。结果表明:在哥氏力和离心力作用下,受限空间内存在射流偏转、径向二次流动以及二次冲击等现象;流动的径向受限可抑制射流偏转,强化冲击换热;相同的旋转数Ro下,逆转向冲击(叶背区)换热努赛尔数Nu比顺转向冲击(叶盆区)高8%。在研究的参数范围内,数值模拟和试验结果说明径向受限周向出流结构能有效的抑制旋转对换热的削弱。     

高湍流度时全气膜涡轮叶片表面冷却和换热特性的实验研究

为获得高主流湍流度时全气膜涡轮叶片表面的冷却和换热特性,在跨声速风洞中实验研究了质量流量比(MFR)和主流雷诺数(Re)对叶片表面气膜冷却效率和换热系数比的影响。在叶片前缘布置了5排圆形孔,在吸力面和压力面分别布置了3排和6排圆形孔,实验结果由嵌入在叶片中截面的热电偶测得。实验中基于弦长的主流雷诺数的范围为3.0×105～9.0×105,叶栅出口马赫数Ma为0.8, MFR的范围是5.5%～12.5%,主流湍流度Tu为14.7%。实验结果表明:主流雷诺数升高显著增强了叶片表面的换热,使层流边界层到湍流边界层的转捩位置提前。对于吸力面S/C0.2的区域(S/C为当地弧长与弦长之比),气膜冷却效率受MFR影响明显,当MFR大于7.7%时提高MFR会导致气膜冷却效率降低;该区域的换热系数比在中低雷诺数时受MFR影响较小,在高雷诺数时随MFR升高而升高。压力面S/C-0.7区域的气膜冷却效率随MFR升高而升高,-0.7S/C-0.4区域的气膜冷却效率受MFR影响较小,对于整个压力面而言,MFR升高提高了叶片表面的换热系数。相对于叶片其它区域,压力面后半段区域和吸力面的气膜冷却效率受雷诺数影响较大。     

离心叶轮内三维湍流流场的LDV测量

应用激光多普勒测速仪和旋转编码器实现了对后弯闭式离心叶轮内三维湍流流场的测量。较详细地描述了设计工况下回转流面上主流速度的分布特点 ;绘出了径向流面上速度矢量分布图 ,结果表明在旋转流场中 ,流体粘性对气流方向具有不可忽视的作用 ;实验测量还得到了在曲率和旋转因素的综合作用下 ,流道中部和出口处产生的由强到弱的二次流涡 ,涡面垂直于主流方向     

低含气率影响下圆柱动静压轴承特性分析

为研究低含气率条件下,温度和气泡含量对润滑油黏度的影响,采用涡轮油和氮气制备油气两相流,在不同温度下采集含气率和动力黏度,拟合得到两相流的黏度计算模型.针对带有深浅腔的圆柱动静压轴承,基于拟合的黏度模型,结合有限元法和有限差分法联立求解油膜Reynolds方程、能量方程及两相流等效密度模型,分析了不同偏心率下承载力、摩...     

流动偏转器对机翼失速特性的控制

介绍了一种控制机翼失速特性的新技术--流动偏转器.通过对二维翼型的数值模拟,研究了流动偏转器高度对控制效果的影响.采用PIV测量作为流动显示手段,验证了较低雷诺数(Re=6.32×10<'5>)时流动偏转器对机翼流动分离的控制效果.基于数值模拟结果和止交设计方法,对流动偏转器在更高雷诺数(Re=1.76×10<'6>)下进行了风洞测力实验研究.研究结果表明,流动偏转器可以有效控制机翼失速特性,能够抑制机翼大攻角下的流动分离,推迟失速攻角和增加升力.对测力实验结果的止交分析还给出了以16°到30°攻角范围内平均气动力最佳为目标的最优水平组合.     

节流制冷对低温贮箱主流液体控温影响分析

为研究低温流体节流特性及冷量引入对贮箱主流体控温影响,首先对低温流体节流干度和体积含气率进行了分析,结果表明质量占比较小的气相占据了大部分空间体积,对流动速率及换热产生较大影响;建立了节流制冷性能测试平台,采用液氮工质开展了节流前压力为0.3～0.37 MPa工况下的减压降温试验,节流前后降温达到了11.3～14.2 K;在集成节流阀孔的热力学排气系统（TVS）系统中,通过节流制冷使贮箱流体产生了平均6.5 K的温降,将贮箱压力控制在150～160 kPa范围;冷量的引入使主流区液体温度持续波浪式降低,气液界面热分层处的降温效果更加明显。     

旋转通道内流阻特性实验的实现与验证

针对涡轮叶片的蛇形内冷通道内流阻特性的研究,在北京航空航天大学航空发动机气动热力国家级重点实验室的旋转涡轮叶片内冷通道换热实验台上构建了旋转工况的测压系统.该测压系统具有高精度、多路选通、高压高旋转数等特点.在通道进口雷诺数从20000~70000,旋转数从0~1.025的范围内,实验研究了旋转状态下,冷态与热态流场下方形截面光滑U形通道流阻系数.实验结果与国外同类实验对比验证了构建的实验系统的可靠性和优越性.实验结果表明:低雷诺数下静止工况的流阻随雷诺数增大而增大,并在雷诺数增大到一定值后转而减小.冷态下流阻随旋转数增大而增大,低旋转数下旋转对热态流阻影响并不显著,高旋转数下热态流阻随旋转数增大而显著增大.      

Accurately predicting hypersonic transitional flow on cone via a symmetry approach

A new algebraic transition model is proposed based on a Structural Ensemble Dynamics (SED) theory of wall turbulence, for accurately predicting the hypersonic flow heat transfer on cone. The model defines the eddy viscosity in terms of a two-dimensional multi-regime distribution of a Stress Length (SL) function, and hence is named as SED-SL. This paper presents clear evidence of precise predictions of transition onset location and peak heat flux of a wide range of hypersonic Transitional Boundary Layers (TrBL) around straight cone at zero incidence, to an unprecedented accuracy as validated by over 70 measurements for varying five crucial influential factors (Mach number, temperature ratio, cone half angle, nose Reynolds number and noise level). The results demonstrate the universality of the postulated multi-regime similarity structure, in characterizing not only the spatial non-uniform distribution of the eddy viscosity in hypersonic TrBL on cone, but also the dependence of the transition onset location on the five influential factors. The latter yields a novel correlation formula for transition center Reynolds number which takes similar functional form as the SL function within the symmetry approach. It is concluded that the SED-SL model simulates TrBL around cone with uniformly high accuracy, and then points out to an optimistic alternative way to construct hypersonic transition model.