冰晶粒子运动过程中的相变特性

针对冰晶在温暖环境下运动时固液汽耦合的相变现象，基于欧拉法建立了冰晶粒子的运动相变模型和计算方法。计算了冰晶粒子在强迫对流环境下的融化相变过程，与实验结果对比验证了运动相变模型和计算方法的准确性。针对NACA0012翼型计算了冰晶绕流运动时的热力学特性，得到了冰晶粒子到达撞击表面时的融化状态与收集系数。研究了冰晶粒径大小、初始球形度、气流相对湿度和温度对运动相变的影响。结果表明：冰晶粒子运动相变模型可以有效地评估冰晶结冰风险，冰晶粒子的融化速率主要取决于粒子直径、球形度、气流温度、湿度等因素，环境温度为288 K时冰晶粒子的融化时间为27.5 s，而相同条件下环境温度为302 K时的融化时间仅有5.2 s。     

翼面结冰过程中的冰晶运动相变与黏附特性

冰晶的运动、相变、黏附现象广泛存在于航空发动机内部结冰过程中，开展相应的计算方法研究，系统分析冰晶运动相变与黏附特性是研究发动机内部结冰、保证飞行安全的迫切需求。针对冰晶运动相变及黏附特性计算，建立了拉格朗日框架下冰晶运动-传热传质耦合的数值计算方法，分析了基于Monte Carlo方法的冰晶撞击与黏附收集系数计算方法，基于NNWICE平台开发了相应的计算程序。以NACA0012翼型为对象，计算研究了扁椭球状、六角平板状和长椭球状3种非球状冰晶颗粒与球状冰晶的运动、传热传质过程，系统分析了来流温度等参数对冰晶黏附特性的影响，得到了冰晶形状对冰晶运动轨迹和融化过程的影响规律及来流总温与液态水含量/总水含量（LWC/TWC）对冰晶黏附特性的影响规律。相关工作为进一步分析混合相结冰数值模拟计算奠定了工作基础。     

气固两相喷嘴平面叶栅能量损失的试验研究

通过对气固两相压气机平面叶栅后,沿栅距和叶高分布的气流压力和速度的测量与分析,研究了气固两相流流经压气机平面叶栅时的速度系数和能量损失系数随粒子浓度和攻角改变而变化的规律。试验是在攻角i=-10°、0°气流速度c=10m/s及粒子质量浓度α=0％、2.0％和2.7％的工况下进行的。结果表明:气固两相流流经压气机平面叶栅时能量损失比单相纯气流的大;而且随着粒子浓度和攻角增大,两者的差别也愈大。     

砂尘对压气机叶片的侵蚀及性能影响

为了应对飞机在起飞和降落的过程中会面临多砂环境，尘土、砂粒等异物颗粒会随空气被吸入航空发动机内的情况，明 确压气机部件的磨损和砂粒的吸入对压气机特性产生的影响，以航空发动机吞砂作为研究背景，以跨声速轴流压气机转子NASA Rotor 37为研究对象，采用两相流动数值模拟的方法对砂粒的运行轨迹、侵蚀分布和特性变化开展仿真研究。获得了不同的吞砂 条件和转速工况对颗粒运动轨迹、侵蚀分布和侵蚀量的影响，并分析了砂粒在不同吞砂条件和不同转速下引起的压气机性能和内 部流场分布的变化规律。结果表明：大部分颗粒进入压气机后，随气流进入叶栅通道，在叶栅通道内与叶片压力面发生碰撞；吞砂 导致压气机的增压能力减弱，对于直径为20、150、300 μm的吞砂工况，其增压范围分别缩小了3.83%、2.13%和1.01%。     

进口畸变条件下轴流压气机转子内流的PIV研究

用50%堵塞比的插板畸变屏模拟低速轴流压气机进口流场畸变,用PIV(粒子图像速度场仪,ParticleImageVelocimetry)技术测量均匀和畸变进气条件下压气机转子叶片通道内的流场结构。实验表明,进气畸变条件下压气机流场在转子叶片通道内不同的周向位置会产生四种不同的流动结构状态,同时,进气畸变大大增加了漩涡扰动、流动亏损和转子叶排气流分离的强度,这样必然会对压气机的性能产生很大的影响。实验也验证了经典平行压气机理论的缺陷。      

航空发动机结冰试验中水滴特性参数分析

通过建立水滴在低温低马赫数气流中的传热传质以及运动数学模型，和以Fluent软件解析雾化喷嘴周围气流的总压分布和速度分布的基础上，分析了水滴直径、水滴速度、水滴表面温度和液态水含量等特性参数，在不同的来流速度、来流温度和水滴初始温度条件下的变化。结果表明，在来流速度、水滴初始速度以及直径大小相同条件下，水滴与气流的温差，是影响发动机进口截面液态水含量的主要因素。对结冰试验过程中模拟结冰云雾气象条件的参数选取进行了建议。     

轴流压气机叶尖射流扩稳数值计算模型

基于体积力方法将压气机叶片对气流的作用力和做功简化为源项,对环形无叶片通道求解带源项的Navier-Stokes(N-S)方程,实现了压气机全通道流场的计算,通过在转子上游机匣施加射流建立了轴流压气机叶尖射流扩稳数值模型.模型中的源项根据叶片排进口气流参数进行变化,通过判断施加的人工扰动变化确定压气机失稳边界.采用该模型对某跨声速单级压气机在有、无射流下的性能和稳定性进行了计算分析,无射流下计算的该压气机性能和稳定边界与试验结果都能较好地一致,其中稳定边界流量仅相差0.7%;有射流下射流对压气机具有明显的扩稳效果,将计算结果与国外公开文献的结果进行了对比,验证了模型的可靠性.      

大型结冰风洞云雾场适航应用符合性验证

结冰风洞云雾场符合性是大型结冰风洞适航应用的基础。为验证3 m×2 m结冰风洞云雾场符合性，发展了基于SAE ARP5905的结冰风洞云雾场符合性验证方法，针对主试验段构型，开展了结冰云雾场符合性验证试验，获得了试验段内液滴尺寸和液态水含量拟合关系，考察了喷嘴水压、液滴尺寸、试验段气流速度和喷嘴数量对试验段液态水含量的影响，形成了主试验段结冰云雾控制包线。结果表明：主试验段内液滴尺寸分布具有显著的单峰分布特征，体积中值直径（MVD）模拟范围近似在10~75 μm之间；试验段中心处液态水含量随着喷嘴水压和MVD的增大而增大，同时近似与试验段气流速度成反比，与喷嘴数量成正比；增大喷嘴水压和喷嘴数量会提高试验段内云雾液态水含量空间均匀性，但是增大气流速度却会减弱试验段内云雾空间均匀性；3 m×2 m结冰风洞主试验段结冰云雾控制包线可以覆盖大部分适航条例25部附录C结冰气象条件，但对低液态水含量结冰条件的模拟仍存在局限。     

进口畸变下缝式机匣处理改善轴流压气机性能的研究

为了研究进口畸变下缝式机匣处理改善轴流压气机性能的机理,以亚声速单级轴流压气机的孤立转子为研究对象,在53.4%设计转速下,采用10通道非定常的数值模拟方法,开展进口畸变下轴向倾斜缝机匣处理改善轴流压气机性能的机理研究。采取在进口径向延伸段中部沿周向设置栏杆的方式产生畸变。与进口均匀实体壁机匣相比,进口畸变实体壁机匣的失速裕度改进量和峰值效率改进量分别为-5.88%和-1.44%,而进口畸变带机匣处理的失速裕度改进量和峰值效率改进量分别为24.37%和-1.09%。进口均匀实体壁机匣时,压气机的失速类型为典型的突尖型失速,即由叶顶间隙泄漏流引起的堵塞所造成。进口畸变引发叶顶处吸力面的气流分离,且气流分离沿周向呈非轴对称分布,进口畸变后叶顶通道内的低速区在很大程度上还是由叶顶间隙泄漏流所造成,只不过相比于进口均匀,叶顶处吸力面的气流分离提前发生,并且反过来又加剧叶顶间隙泄漏流的负面影响,即进口畸变后压气机失速的主要诱因没有发生本质的变化。机匣处理上、下游的喷射、抽吸过程分别对叶顶通道上游、中部的低速区气流起到激励、吹除作用,其有效地抑制了由叶顶间隙泄漏流所造成的叶顶前缘溢流和叶顶通道内压力面附近的回流,进而提高了叶顶通道的流通能力。     

端弯叶片通道出口二次流的研究

实验并且采用考虑粘性的涡传输方程计算端弯压气机叶片通道出口二次流。给出了３种端弯叶片及对应的原始直叶片通道出口气流偏转角的实验值和计算预测值。研究表明：采用上述方程能比较准确地预测出口气流偏转角；叶片后缘减弯可减小流向二次涡强度。     

高主流湍流度下密度比对涡轮导叶全气膜冷却特性的影响

为了研究高主流湍流度下二次流密度比对涡轮导叶全气膜冷却特性的影响,使用热色液晶测量了在主流湍流度为15%,二次流密度比为1.0和1.5下三维涡轮导叶的气膜冷却效率和换热系数。二次流与主流质量流量比为7.0%和12.5%。结果表明:二次流密度比增大可以降低冷气射流的动量,小流量比工况下,在叶片前缘和压力面前半段,动量较低的二次流在高主流湍流度的影响下更易耗散,增大二次流密度比使冷却效率明显降低;大流量比工况下,二次流动量降低使气膜孔后区域冷气贴附性增强,气膜冷却效率和冷气覆盖效果均得到提升。小流量比工况下,二次流密度比增大对叶片表面换热的影响较小;大流量比工况下,二次流密度比增大使吸力面中弦区域和压力面后半段的平均换热系数比分别降低15%和25%。     

大型结冰风洞中冰晶热/力平衡特性数值研究

为明晰大型结冰风洞中冰晶热/力平衡特性,发展了基于欧拉法的冰晶运动和传热传质耦合计算方法,模拟了典型大型结冰风洞构型内冰晶运动和传热过程,从沉降收缩、动量平衡和热平衡三个方面,考察了颗粒形状和体积密度的影响.结果表明:颗粒形状和体积密度对小尺寸冰晶的热/力平衡特性无显著影响.降低颗粒球形度和体积密度会抑制大尺寸冰晶沉降...     

高马赫数涡轮发动机射流预冷特性研究

射流预冷涡轮基发动机在高空高马赫数工作时对冷却水和液氧具有迫切的需求。本文以气液相变冷却机制为切入点，开展高空模拟试验进气预冷段内水-液氧射流冷却的数值分析，考虑真实雾滴颗粒运动的热力现象，基于欧拉-拉格朗日多相流方法解析气液两相热质传输过程，分析水-液氧混合射流对高马赫数涡轮发动机预冷段内流动及换热特性的影响规律。结果表明，水-液氧射流雾化蒸发的效果具有即时性，基于水雾-水蒸汽比热大和汽化焓高的特点，水雾浓度对主流总温降和总压恢复占主导性；而液氧浓度有利于降低湿空气的热流密度。在射流浓度2%-8%时，预冷段总压降系数为0.84%-1.27%，总温降系数范围为2.15%-15.12%，即温降范围为12.92K-90.89K。为平衡高空高马赫数时冷却水和液氧的需求，需控制水-液氧的射流比例，液氧射流量建议小于60%的总射流浓度。在“40%水-60%液氧”的射流比例时预冷段内流动和传热特性达到局部最优。在发动机物理转速不变时，射流冷却后预冷段内湿空气来流质量流量增幅0.22%-9.39%，其中空气和水蒸气含量的贡献份额分别约为71.8%和28.2%。因此，射流预冷有利于涡轮发动机在高马赫数时具有更高的加速度。     

液氮抽真空制冷的理论及试验

在考虑液氮总质量变化、储罐罐体热容的前提下,基于化学势温差,传热系数和传质系数等参数,详细建立了液氮降压制备过冷液氮的数学模型。理论计算的液氮温度及质量变化和制取过冷液氮试验值吻合较好,结果表明:该数学模型能准确描述液氮相变过程,且获得了液氮抽真空制冷中相变传热的化学势温差为1.21 K、传热系数为2 942.37 W/(m2·K)及传质系数为0.004 4 m/s。      

三维机翼结冰模拟

针对机翼结冰问题,利用欧拉法求解空气水滴两相流流场,通过Fluent软件用户自定义函数功能编程求解三维机翼表面水滴局部收集系数以及结冰情况.在计算中.将水溢流方向沿机翼弦向和展向分解,以此对Messinger结冰热力学模型进行改进,使之更符合三维情况,并在此基础上进行三维传质传热分析;修正传统的二维表面对流换热系数求解...     

R141b在矩形微尺度通道中的两相流传热特性

设计搭建水力直径分别为1mm和0.5mm的矩形微尺度通道实验台,研究了以R141b型制冷剂作为工质的两相流沸腾传热特性。实验取热流密度为1～16kW/m2、质量流速为111.1～333.3kg/(m2·s)和质量干度为0～1,分析了三者对平均传热系数的影响,探究影响换热的主导因素。结果表明:热流密度较高时,平均传热系数随热流密度增加而减小,流动换热主要受到沸腾传热的影响;当质量流速较大且热流密度较低时,平均传热系数随热流密度增加而有所增长;热流密度较低时,平均传热系数随质量流速变化明显,热流密度升高到一定值后,质量流速对平均传热系数的影响很小;当质量流速处于111.1～333.3kg/(m2·s)时,平均传热系数随质量干度的增加而减小。      

基于多孔介质局部非热平衡方法的刷式密封耦合传热特性

建立了刷式密封流动传热与摩擦接触的多物理场耦合分析的计算模型.推导了描述刷式密封在局部非热平衡条件下耦合传热过程的双能量方程,该方程既考虑了刷丝与气流的对流换热,同时也考虑了刷丝束内导热过程的各向异性.采用所建立的耦合计算模型数值预测了刷式密封在不同压比和转速下的温度分布.结果表明:刷式密封的摩擦热效应使刷丝与转子相接触的附近区域形成局部高温,刷丝的高温区域主要集中在围栏高度以下;泄漏气流与刷丝间强烈的对流换热使刷丝的温度分布和气流的温度分布非常相近;刷丝的最高温度随压比增大而升高,但围栏高度以下区域的平均温度随压比增大而降低;刷丝的最高温度和围栏高度以下区域的平均温度均随转速升高而升高;摩擦热效应引起的泄漏气流的温升效应使刷式密封泄漏量随转速升高而减小.      

边界层吹吸气对高负荷扩压叶栅性能的影响

采用边界层流动控制能够有效抑制扩压叶栅的流动分离。以某大弯折角低稠度扩压叶栅为研究对象,利用数值模拟手段研究了原型、叶片表面边界层单独吹气以及吹吸气相结合等边界层控制手段下的流场和叶栅性能变化情况。结果表明,无论是单独吹气还是吹吸气相结合的边界层控制方法,都能有效控制扩压叶栅中的边界层分离,从而较大幅度地增大叶栅负荷,并降低气动损失;计算表明,吹气和吸气的效果不尽相同,且吹吸气口位置及吹吸气流量对边界层的流动亦有较为明显的影响。其中采用1.7%的吹气流量,结合1.38%的吸气量,可以使静压增压比提高15%以上,而损失系数降低至原型的20%以内。     

飞机结冰过程的液/固相变传热研究

基于液/固相变的基本理论, 对飞机结冰过程的液/固相变传热特性进行了研究, 建立了液膜及冰层生长的分析模型, 并采用该模型对来流条件各参数对冰层生长特征及其速率的影响进行了分析.研究表明, 在飞机结冰过程中, 冰层生长速率及不同冰型的形成不但受来流参数的影响, 而且与固、液相区的内部传热特性相关, 因而考虑各相区的传热过程是提高冰层生长预测准确性的重要途径.所建模型考虑了液/固相变的非线性特征, 对冰层生长速率及冰型特征有着较好的预测作用.      

单孔冲击式帽罩前缘流动换热特性数值分析

为了提高航空发动机帽罩冲击防冰结构的设计分析水平,对单孔冲击式帽罩前缘结构的流动换热特性进行数值研究,分析了不同冲击孔径与不同冲击雷诺数对帽罩前缘速度流场、换热系数与努塞尔数的分布规律。结果表明:在冲击雷诺数一定的条件下,冲击孔径越大,射流核心速度和前缘壁面附近的气流速度越小,前缘冲击区形成的涡流团越大,当孔径D=6 mm时,小孔径冲击下前缘区整体换热效果不如大孔径的,而在滞止区的换热效果则要优于大孔径的;当D>12 mm时,孔径大小对壁面换热基本没有影响;在冲击孔径相同时,增大冲击雷诺数使得冲击射流、前缘壁面附近及侧壁曲面通道内的气流流速增大,冲击区内的涡流团则逐渐减小;冲击雷诺数的增大也增强了前缘冲击区的换热特性。