粘性气体中粘性液体射流分裂与雾化机理研究

采用线性稳定性分析的方法对粘性气体中的粘性液体射流的分裂与雾化机理进行了分析,数值计算表明:液体射流分裂与雾化过程中存在一临界气体韦伯数We2c=1,We2＜We2c对应的是射流分裂过程,We2＞We2c对应的是射流雾化过程,射流分裂过程和雾化过程的机理有所不同.当We2＜1时,We2对射流分裂过程具有稳定性的作用;当We2＞1时,We2对射流雾化过程起着不稳定性的作用.液体Reynolds数Re1在整个射流过程中始终起着不稳定性的作用,气体Reynolds数Re2的作用却相反.气液密度比Q,即气动力在整个射流分裂与雾化过程中始终起着不稳定性的作用.     

双毛细管法高温气态环己烷黏度测量研究

为了探索双毛细管法在碳氢燃料高温气态黏度测量中的应用,本文采用双毛细管法测量了亚临界压力下环己烷的高温气态黏度。测量温度范围：301.9 K-567.1 K（间隔约10K）;压力范围：0.5MPa-4.0MPa（间隔1.0 MPa）。考虑螺旋管圈的二次流效应,使用基于雷诺数（Re）或迪恩数（De）的函数关系式对黏度测量结果进行离心修正。在测量参数条件下,气态环己烷的雷诺数Re ≤ 1984.5。研究表明：双毛细管法测量环己烷液态和气态黏度时,螺旋管圈的二次流效应具有显著差别;液态环己烷的离心修正系数小于 4.32%,气态环己烷的离心修正系数最高达到 49%;离心力二次流效应是影响环己烷气态黏度测量精度的最重要因素;离心修正关系式的准确度对环己烷气态黏度测量结果的不确定度产生显著影响。     

气液比对气液同轴离心式喷嘴缩进室内流动过程的影响

为了研究气液同轴离心式喷嘴缩进室内部非定常流动过程,采用Level Set和VOF相耦合的方法结合网格自适应技术对缩进长度为8 mm的液体中心式气液同轴离心式喷嘴流动过程进行了数值仿真研究,计算得到了较为精细的液膜一次破碎过程、流场结构和压力振荡特性。结果表明：液膜的破碎模式受气液比的影响较大,随气液比的增加,液膜破碎模式由曲张表面波主导的破碎变为穿孔破碎。此外,清晰获得了自激振荡过程,分析了缩进室内部压力场及速度场分布特征,研究发现随着气体压降增加,气体环缝出口会出现膨胀波和激波,形成一个“扇环形”的超声速流场区域；激波后气流分离,出现旋涡,形成局部高压区,旋涡中心随激波面周期性地上下移动,致使局部压力出现周期性振荡。     

幂律流体液膜破裂的线性稳定性分析

为对凝胶推进剂撞击式喷嘴产生的液膜的破裂行为进行研究,使用线性稳定性分析方法并结合幂定律本构方程对幂律流体液膜的稳定性进行了研究,并对流体物性参数、流变参数等对液膜流动稳定性的影响规律进行了分析.结果发现:稠度系数、流动指数和表面张力越大,液膜表面扰动波的增长率越小,波长越大.而较高的液膜速度、气液密度比和液膜厚度都使扰动波增长率变大,但波长随液膜速度和气液密度比的增大而减小,而液膜厚度的情况却相反.理论计算与实验值的对比表明,线性稳定性分析能较准确地预测液膜的破碎长度,但预测的表面波长却大于实验值.      

旋转状态涡轮叶片吸力面单孔气膜冷却实验

采用稳态液晶测温方法,系统研究了1.5级涡轮叶片吸力面在旋转状态下的气膜冷却特性.实验中,主流经加热压缩后冲击涡轮转动,基于动叶弦长的涡轮进口主流雷诺数为8×104.射流分别采用空气和二氧化碳,其对应射流-主流密度比分别为1.03和1.57.实验转速为630,700 r/min和737 r/min,对应旋转数分别为2.092,2.324和2.448.吹风比从0.3到3.0变化.结果表明,吸力面上,气膜冷却效率随吹风比的增大先上升后下降,存在一个最佳吹风比,使冷却效果最好;增大密度比有利于增加气膜覆盖面积;旋转降低了气膜冷却效率;气膜向低半径方向偏转,但并不十分明显.      

横流作用下气液两相环状流射流液膜的雾化特性实验研究

为了更好地认识横向气流对气液两相环状流射流雾化过程的影响,本文采用高速摄像等技术,针对气液两相环状流射流液膜在横向气流中的破碎与雾化特性开展了实验研究。研究发现,气液两相环状流射流在横流作用下能够实现稳定的雾化,射流液膜的破碎与雾化具有周期性和不连续性特征;查清了雾化过程中射流液膜存在三种不同的破碎模式：爆式破碎、分段式破碎和环膜破碎,并对每种破碎模式下的液膜破碎特征进行了研究;结合实验结果统计分析,获得了射流液膜不同破碎模式的动力学条件和变化规律;同时对环状流射流液膜不同破碎阶段的雾化液滴的粒径分布进行了统计分析,发现三种破碎模式下,液膜的爆式破碎产生的雾化液滴粒径更小,雾化效果较好,同时提高环状流表观气速和横流速度也能够促进射流液膜的雾化。     

旋转状态下孔排构型气膜冷却特性实验研究

 通过对旋转状态下带有孔排构型平板叶片模型表面的气膜冷却换热特性进行实验研究,得到了不同吹风比、旋转数、主流雷诺数和密度比情况下气膜孔下游流向位置上绝热效率的变化规律,同时提出了一个新的用于评价气膜冷却效果优劣的无量纲参数——有效覆盖率(EAR)。实验采用空气和二氧化碳作为冷却工质,利用热色液晶测温（TLC）技术对测试表面的二维温度场进行测量,并通过旋转拍照系统对液晶图片进行即时采集。研究结果表明,吹风比对旋转状态下二次流的流动形态起着决定性作用。随着吹风比的增大,二次流经历了附壁流动、分离流动和分离再附壁流动3种形态。当旋转数增大时,吸力面和压力面上的绝热效率呈现先上升后下降趋势,且压力面上整体绝热效率要高于吸力面对应值。主流雷诺数的增大导致壁面绝热效率出现降低趋势,但整体冷却效果变化不大;冷气对主流密度比的增大则有助于壁面冷却效果的改善。此外,EAR可以对不同工况下的壁面气膜冷却效果进行量化比较和评判,因而极具工程推广价值。     

低韦伯数非牛顿射流撞击破碎直接数值模拟

非牛顿射流的撞击破碎在液体火箭推进系统中被广泛用于燃料的喷注雾化,然而人们对其破碎物理机制却知之甚少。本文将采用基于液体体积法的直接数值模拟(DNS)工具,研究夹角为90°的2个等直径低韦伯数射流撞击现象,并分析二者形成的单一对角射流特征及其破碎机理。研究结果表明,撞击形成的单一对角射流直径较原射流直径大1.66倍,并在头部形成液滴诱导破碎的发生。除了头部破碎,在对角射流的发展过程中还观察到一类液柱破碎,表现为射流表面不稳定波不断发展形成新的弯曲波破碎,并产生卫星液滴及液滴的融合。伴随两股射流撞击的发生,气液两相交界面的面积也不断减小,同时,射流内部的黏性也不断变化,在本文的低雷诺数和低韦伯数条件下,流体内部黏性系数变化超过10%。     

吹风比对旋转涡轮叶片气膜冷却的影响

对1.5级涡轮叶片在旋转状态下不同吹风比时的气膜冷却特性进行了实验研究.实验中基于动叶弦长的涡轮进口主流雷诺数为1.6451×105,冷却工质采用二氧化碳,对应主流射流密度比为1.57,实验涡轮转速为475 r/min,对应旋转数为1.901,吹风比为0.5~2.0.采用稳态液晶方法测温.结果表明:①压力面上,随吹风比的增大,气膜冷却效率升高,气膜覆盖区域增大,气膜轨迹的偏转程度减弱;②吸力面上,随吹风比的增大,气膜冷却效率先上升后下降,气膜覆盖区域亦先增加后减少,气膜轨迹的偏转程度不明显;③射流流动的曲率半径影响气膜对壁面的附着.      

液体中心式同轴离心喷嘴液膜破碎特性仿真研究（液体火箭推进）

为了实现对不同工况下液体中心式同轴离心喷嘴液膜破碎特性的数值模拟研究,采用网格自适应加密技术、耦合的Level-set和Volume of Fluid（CLSVOF）方法对气液界面进行捕捉,利用改进延迟分离涡模拟（IDDES）方法模拟湍流。分析了液膜的破碎模式、喷雾锥角、破碎长度以及流场特性。通过观察分析得到：随着气液相互作用的增强,液膜破碎模式依次经历主导表面波发展导致的液膜破碎、Rayleigh-Taylor（R-T）和Kelvin-Helmholtz（K-H）不稳定性引起的液膜破碎,以及气动破碎模式。随着气液动量通量比（Momentum Flux Ratio,MFR）的增大,喷雾锥角和破碎长度逐渐减小且呈渐进趋势,发现无量纲喷雾锥角和破碎长度均与MFR_^-A成正比例关系。相同液膜破碎模式而不同工况时,主要流场特征一致。     

超声悬浮乙醇液滴气动破碎特性研究

采用超声波悬浮液滴的方式,选取直径分别为1.0mm、1.5mm和2.0mm的乙醇液滴为研究对象,在韦伯数We=15~90的范围内,对其在脉冲气流中的破碎特性进行了实验研究。利用高速摄像机拍摄破碎现象,分析了液滴的破碎模式与运动特性。研究结果表明,随着韦伯数的增加,液滴先后产生袋状破碎、袋状/蕊心破碎、羽状/液膜稀释破碎。三种不同的破碎模式本质上都是袋状破碎——边缘袋状破碎与核心袋状破碎,低We数下主要为核心袋状破碎,高We数下主要为边缘袋状破碎。We数越大,破碎形成的袋直径越小,数量越多,破碎越剧烈。初始直径增大会使液滴转变破碎模式的临界We数增大。脉冲气流中乙醇液滴的迎风面运动遵循匀加速规律,其无量纲位移以We-1/2的比例缩小后是无量纲时间的简单二次函数。液滴在横向的扩散过程分为两个阶段,一个是惯性控制阶段,另一个是毛细效应阶段。     

Flow and Heat Transfer Characteristics in Rotating Two-pass Channels Cooled by Superheated Steam

In a modern gas turbine,using superheated steam to cool the vane and blade for internal convection cooling is a promising alternative to traditional compressor air.However,further investigations of steam cooling need to be performed.In this paper,the three-dimensional flow and heat transfer characteristics of steam are numerically investigated in two-pass square channels with 45° ribbed walls under stationary and rotating conditions.The investigated rotation numbers are 0 and 0.24.The simulation is carried out by solving the Reynolds averaged Navier-Stokes equations employing the Reynolds stress turbulence model,especially considering two additional terms for Coriolis and rotational buoyancy forces caused by the rotating effect.For comparison,calculations for the air-cooled channels are done first at a Reynolds number of 25 000 and inlet coolant-to-wall density ratio of 0.13.The results are compared with the experiment data.Then the flow and heat transfer in steam-cooled channels are analyzed under the same operating conditions.The results indicate that the superheated steam has better heat transfer performance than air.Due to the combined effect of rotation,skewed ribs and 180° sharp turn,the secondary flow pattern in steam-cooled rotating two-pass channels is quite complex.This complex secondary flow pattern leads to strong anisotropic turbulence and high level of anisotropy of Reynolds stresses,which have a significant impact on the local heat transfer coefficient distributions.     

液晶瞬态技术测量带侧向流扰流柱通道端壁换热

采用热色液晶瞬态测量技术测量带侧向流扰流柱通道端壁全表面换热系数的分布,研究了侧流比及雷诺数对换热的影响,其中,侧流比为0.25~1.0,雷诺数为3×104~9×104。结果表明:(1)侧流比对扰流柱通道的流动形态及端壁换热有重要影响;(2)存在一个临界侧流比,在临界侧流比以下,流动形态沿主流方向呈错排流状态;在临界侧流比以上,流动形态沿侧流方向呈错排流状态;在临界侧流比附近,流动为顺排流动状态,方向在主流和侧流方向之间;(3)侧流比较小或较大时,扰流柱通道端壁换热较强;在临界侧流比附近,换热相对较弱。     

Atomization mechanism of a charged viscoelastic liquid sheet

In order to study atomization mechanism of a viscoelastic liquid sheet in an electric field,the spatial–temporal stability analysis of a viscoelastic liquid sheet injected into a dielectric stationary ambient gas in the presence of a vertical electric field is conducted. The dispersion relations of both sinuous and varicose disturbance modes are solved to explore the spatial–temporal instability of a charged viscoelastic sheet, by setting both the wave number and frequency complex. A parametric study is performed to test the influence of the dimensionless parameters on the absolute instability of the sheet. The results show that the increase of liquid Weber number and time constant ratio, or decrease of gas to liquid density ratio and Reynolds number, can damp the absolute instability. The effect of the liquid elasticity depends on the value of time constant ratio: when time constant ratio is small, the increase of liquid elasticity could amplify absolute growth rate, but the effect is weak when the elasticity number is relatively large; when time constant ratio is large, the increase of liquid elasticity cannot affect the absolute growth rate. Moreover, the variation of electrical Euler number can hardly influence the absolute instability of a charged viscoelastic sheet.     

小高径比扰流柱冷却通道的换热和流动特性

采用数值模拟的方法,对涡轮叶片尾缘处圆形小高径比扰流柱冷却通道的换热和流动特性进行了研究,分析进口雷诺数和扰流柱间距对冷却通道换热和流动特性的作用过程.结果表明:进口雷诺数的提高能够有效改善冷却通道端壁的换热性能,但这种改善能力随着进口雷诺数的提高而逐渐减弱,同时降低冷却通道的压力损失系数.在两种扰流柱间距中,流向间距是影响端壁换热性能的主要因素,随着流向间距的减小,冷却通道换热性能逐渐变好,压力损失系数降低;横向间距是影响冷却通道流动损失的主要因素,两者大小成反比关系.在通道计算中,扰流柱平均换热性能约是端壁平均换热性能的1.8倍,端壁换热权重约是换热面积比0.824倍,同时该权重几乎不受进口雷诺数的影响.      

小展弦比机翼加装格尼襟翼的低雷诺数试验

通过风洞试验研究了在低雷诺数下加装格尼襟翼的小展弦比机翼气动特性,机翼展弦比为1.67,格尼襟翼为1%~4%弦长高度,试验雷诺数分别为2.0×105和5.0×105.天平测力和表面测压的试验结果表明:低雷诺数下小展弦比机翼加装一定高度的格尼襟翼后,升力系数明显提高,加装1%弦长高度的格尼襟翼还能够提高机翼的升阻比.这是因为在试验雷诺数下,合适高度的襟翼在提高了机翼升力的同时并未显著增大机翼阻力.对比不同试验雷诺数下格尼襟翼的作用效果,表明格尼襟翼能够减少低雷诺数气流分离的不利影响,并且在较小的雷诺数下这种作用更加显著.关于格尼襟翼对低雷诺数层流分离现象的影响,还需要通过细致的流场显示技术进行研究.      

三棱柱的低雷诺数绕流研究

数值模拟了截面为等腰三角形的棱柱在不可压牛顿流体中的二维绕流问题,分析了不同的截面尺度比以及三角形顶点迎风和背风安置时的尾流稳定性与阻力系数.研究发现对于较大截面尺度比的棱柱,在相同雷诺数的情况下顶点迎风相较于顶点背风的情况更易失稳从而产生侧向力,但阻力系数相对较小,当雷诺数为160时,最多能够减少19.4%.进一步地对比研究了顶点迎风与顶点背风在不同截面尺度比时升、阻力系数及涡脱落频率随雷诺数的变化关系.      

带肋和双排出流孔通道的流动特性

研究了通道进口雷诺数和总出流比对带肋和双排出流孔通道流量系数和压力分布的影响。实验研究的通道入口雷诺数为3×104~1.5×105,通道总出流比为0.09~0.22。结果表明:通道总出流比较小时,流量系数沿流向减小。通道进口雷诺数增加,流量系数先增加,之后基本不变;通道总出流比较大时,流量系数基本不变;各工况下总压系数沿流向依次经历迅速减小、基本不变、继续减小的过程;沿流向各位置上的总压系数在通道进口雷诺数为6×104~9×104时最小;出流比增大,沿流向各位置上的总压系数随之增大。     

旋转冲压发动机高速动静混合气体轴承性能分析

为了满足旋转冲压发动机对高速支撑的要求,本文对动静压混合高速气体轴承进行理论分析与数值研究。首先通过旋转冲压发动机的工作条件确定了气体轴承的供气压力,对描述轴承内气体流动的雷诺方程采用牛顿迭代与有限差分法进行求解,获得不同偏心、不同转速下轴承内气体压力分布并分析动静压耦合机理。同时分析了不同供气孔排数对压力分布与承载能力的影响,给出了不同转速、不同供气孔排数下轴承所能支撑的最大转子重量,为下一步旋转冲压发动机转子系统设计奠定基础。     

再入飞行器和模型转换雷诺数计算

本文采用近似的激波关系、球锥体压力分布、平衡气体状态方程、边界层厚度计算公式及转捩准则等,按照流管法计算了4种再入飞行器的外流转捩雷诺数和转捩距离,并给出了高硅氧、碳石英、碳碳球锥模型的转捩雷诺数范围,对再入飞行器和模型的计算结果作了分析和比较。