结冰风洞中液滴相变效应数值模拟

为明晰结冰风洞中液滴相变效应,发展了基于Euler法的气液两相传质传热耦合流动计算方法,探索了3 m×2 m结冰风洞主试验段构型相变效应对液滴传热过程的影响,开展了参数影响研究,评估了试验段内液滴过冷状态。结果表明：构型内液滴经历了先蒸发后凝结两个阶段,蒸发效应促进了准一维传热阶段液滴温度的下降趋势,使液滴温度趋于湿球温度,而凝结效应则抑制了三维收缩阶段液滴温度的下降趋势,进而增大了试验段液气温差,影响液滴过冷状态;增大初始相对湿度和试验段气流速度,会导致蒸发效应减弱而凝结效应增强,进而增强了相对湿度和气流速度对液滴过冷状态的影响程度,与此相反,初始液滴尺寸的增加,则会导致蒸发效应和凝结效应均减弱,进而减弱了液滴尺寸对液滴过冷状态的影响程度;在典型工况下,小尺寸液滴（液滴直径范围为40~100 μm）在高风速时（试验段气流速度大于164 m/s）将偏离过冷状态（液气温差超过3℃）。     

结冰风洞中过冷大水滴云雾演化特性数值研究

为明晰结冰风洞中过冷大水滴（SLD）云雾演化特性，发展了基于欧拉法的SLD液滴运动、传热和传质耦合计算方法，针对3 m×2 m结冰风洞主试验段水平收缩构型，分析了SLD云雾沉降收缩特性、动量平衡特性和热平衡特性，探索了液滴变形破碎的影响，评估了构型出口处SLD液滴动量平衡和热平衡状态。研究结果表明：直径超过250 μm的SLD液滴在构型内会发生显著形变，液滴尺寸越大则变形程度越强，尤其在160 m/s工况下，当液滴直径超过750 μm后，SLD液滴甚至会发生破碎；液滴变形破碎效应会增大液滴加速度和液滴温度下降率，促使SLD液滴趋近动量平衡和热平衡状态；SLD云雾（最大液滴直径小于1 000 μm）在构型出口处会出现显著的粒径浓度分层、动量分层和热分层现象，其中直径小于100 μm的小尺寸液滴速度快、温度低且不断凝结，趋于平衡态，但直径超过500 μm的大尺寸SLD液滴速度慢、温度高且不断蒸发，则显著偏离平衡态；增大试验段气流速度尽管会减弱SLD云雾粒径浓度分层程度，但会增强动量分层和热分层程度，尤其在160 m/s工况下，SLD云雾会均匀分布在构型出口中心区域内（-0.75 m < Y < 0.75 m且-0.5 m < Z < 0.5 m），与其平衡态间的最大速度差和温度差将分别超过18 m/s和20℃。     

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

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

大型结冰风洞气流场适航符合性验证

大型结冰风洞气流场适航符合性是大型结冰风洞适航应用的先决条件。为验证3 m×2 m结冰风洞气流场适航符合性,首先建立了结冰风洞气流场适航符合性验证方法,然后针对主试验段构型,开展了气流场适航符合性验证试验,考察了试验段气流速度和喷嘴干空气射流对流场特征参数(气流速度、气流偏角和气流湍流度)的影响,最后评估了试验段内气流场品质,获得了结冰风洞气流场控制包线。结果表明:喷雾耙结构会影响试验段内气流速度和气流偏角空间分布形态,进而导致了非均匀峰值区的形成;增大试验段气流速度会改善气流场品质,但喷嘴干空气射流会显著恶化试验段气流速度低于60 m/s的气流场品质;3 m×2 m结冰风洞主试验段气流场品质在主要试验速度范围内均满足适航审定要求。     

基于粒子轨迹的结冰风洞收缩段优化设计数值模拟

为了研究收缩段型面和液滴轨迹特性对结冰风洞流场品质与结冰试验效果的重要影响,获得收缩段的优化设计参数,采用数值模拟方法,使用二维简化模型对3m×2m结冰风洞收缩段的流场特性进行了研究,对几种典型收缩曲线的流场进行了计算,对比了收缩曲线对沿轴向压力分布、出口截面速度分布、不同粒径云雾粒子的极限释放距离及试验段入口云雾粒子分布的影响,结果表明:三次曲线与x_m为0.4的五次曲线不能满足流场合格指标,其他曲线流场特性差别不大,维氏曲线具有最小极限释放距离,在收缩段出口,三次曲线和五次曲线粒子包络面积最大,综合比较,x_m为0.5的五次曲线能有效兼顾流场均匀性、极限释放距离和收缩段出口粒子包络面积比.      

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

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

结冰风洞水滴直径标定方法研究

过冷水滴粒径大小是重要的结冰云雾参数,获知结冰风洞中的水滴直径,是得到定量结冰风洞实验结果的基础。对于结冰风洞内水滴直径单一或者分布比较集中的情况,提出了一种采用数值计算和结冰风洞实验相结合的手段标定水滴直径的方法。该方法首先采用拉格朗日法数值计算水滴运动轨迹,得到撞击极限随水滴直径变化的关系曲线,在此基础上,进行结冰风洞实验,测量实验得到的水滴撞击极限,通过在撞击极限与水滴直径关系曲线上进行插值,进而得到实验水滴直径大小。采用该方法对0.3m×0.2m结冰风洞内的水滴直径进行了标定,分别计算和测量了25m/s和35m/s两种速度条件下的水滴撞击极限,得到的水滴直径值相差不超过1μm,初步说明该方法的合理性。同时,对于结冰风洞内水滴粒径多尺寸分布的情况,还提出了相应的标定其容积平均直径MVD的方法,该方法在计算水滴收集率的基础上,通过测量驻点处的结冰厚度,实现对MVD的测量。采用本文提出的两种方法进行结冰风洞水滴粒径标定,只需要一般的长度测量工具即可进行,操作方便,成本低廉,克服了常规的水滴直径测量或标定需要专门设备的不足。     

超声速气流中简化微小液滴横向喷射的汽化过程探索

为研究超声速气流中简化液滴的汽化过程问题,本文分析了两相流计算中已有的两相传热模型,并对简化液滴绕流开展数值计算.在来流Ma≤0.6的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果相一致,而在来流Ma≥0.9的条件下,数值计算得到的简化液滴-气流之间的传热速率与已有模型得到的结果存在很大偏差.由此建立了考虑简化液滴与气流相对超声速相互作用的两相传热模型.进一步,采用Charles B.Henderson阻力系数关系式与新建立的传热模型,对不同直径简化液滴的运动与汽化开展工程计算.在来流2.7Ma的二维平板超声速流场中选取一个截面,作为气相流场,结果显示,(1)简化液滴与主气流存在相对超声速作用.当简化液滴直径dk≤0.12mm时,作用区域约为0.1m～0.4m,当dk>0.12mm时,作用区域明显增大,(2)简化液滴的穿透尺度不超过0.011m/m(深度/长度),时间尺度约为0.28ms～3ms,(3)简化液滴完成汽化的空间尺度约为0.1m(dk>=0.03mm)、0.45m(dk>=0.05mm)与1.24m(dk>=0.075mm),而当dk>0.09mm时,简化液滴完成汽化的空间尺度则大于1.9m.使用考虑简化液滴与气流相对超声速相互作用的两相传热模型与使用传统的传热模型对简化液滴的运动轨迹没有影响,而对简化液滴的汽化过程有较大影响.     

高性能亚音速平面叶栅风洞设计

论述了用于测量压气机和涡轮叶片在亚声速状态下的二维气动特性的平面叶栅风洞设计问题。风洞采取暂冲式直流型式,由压缩机和储气罐提供高压气流,稳定段前采用大开角扩散段,并在扩散段进口段设置整流锥,锥后加装两层球面整流网,收缩段与稳定段连接处加设唇口,唇口与收缩段型面曲线均为双三次曲线,为防止试验段由于内外压力差导致漏气所造成的压力损失,在试验段外部加装了驻室结构,出口气流通过消声塔降低噪音再排入大气。风洞整体参数居于国内领先水平。     

过冷大水滴条件下结冰相似准则

虑到中国现有结冰风洞在过冷大水滴(SLD)结冰条件模拟能力方面的局限,发展SLD结冰风洞及其相应的缩比方法迫在眉睫。研究了大水滴运动机理,并依据绕流流场相似、水滴运动和撞击特性相似、结冰过程热力学特性相似原理,提炼出SLD结冰条件下的6种相似缩比准则,利用结冰数值模拟软件FENSAP-ICE对所提相似缩比准则进行筛选与比较,从而基于结冰相似百分比量化判断准则找出给定条件下的最优准则法。结果表明:考虑SLD运动特性的撞击参数与ONERA法相耦合的方法以及撞击参数与Ruff法相耦合的方法的计算结果与实际成冰构型吻合良好；此外,撞击参数与ONERA法相耦合的方法在计算直径为171 μm的SLD时优于其他。