基于涡流发生器的风洞侧壁附面层控制

为降低风洞侧壁附面层对半模型数据的影响,在前期数值模拟的基础上,研制了一种适用于2.4m跨声速风洞半模型试验段侧壁的梯形涡流发生器,并进行了试验验证.结果表明:加装涡流发生器效果明显,亚声速范围内能够使附面层厚度降低20%~30%,对主气流均匀性影响可忽略;加装后半模型零升阻力系数降低,升力线斜率增大,压力中心向机身移动,体现了明显的附面层减薄效果,证明所研制的涡流发生器可应用于亚声速半模型试验中.      

一类新的速度入流条件在管道流模拟的应用

基于经典边界层理论,发展了一类简化的、按边界层厚度描述的速度入流边界条件。基于入流条件,对低速后向台阶流动和超声速压缩拐角流动进行了数值模拟。通过细致比较空间流场、边界层内速度型、表面摩擦阻力和压力分布等特性,对该类型边界条件进行了计算确认。结果表明,所发展的简化边界层厚度速度入流条件提法正确、合理,具有描述方法简单、鲁棒性好、适合工程计算的优点,可显著简化不同流速的管道流动的数值模拟。     

吸附式低反动度轴流压气机气动设计方法在多级轴流压气机中的应用

随着对吸附式低反动度轴流压气机内部流动细节的逐步深入,进一步完善了多级吸附式低反动度轴流压气机气动设计思想,完成了3级吸附式低反动度轴流压气机气动设计.三维黏性的数值计算结果表明:在第1级动叶入口叶尖切线速度为370m/s的前提下,通过只在首级静叶和末级静叶中进行附面层抽吸,实现了总压比为6.1的3级吸附式低反动度轴流压气机气动设计.附面层总抽吸流量占入口流量的11.3%.在不考虑附面层抽吸对流动效率影响的前提下,3级效率达到88.1%.      

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

尾缘锯齿降低叶栅噪声的数值模拟

采用大涡模拟与声类比的方法研究了尾缘锯齿对涡轮叶栅噪声的影响.设计了两种不同的尾缘锯齿,对比了Re=3.3×105(基于叶片弦长与叶栅出口速度)下两种不同结构锯齿尾缘叶栅与直尾缘叶栅的声功率.结果表明:尾缘锯齿可以降低叶片吸力面边界层分离噪声约5dB,降低尾缘涡脱落噪声约10dB.进一步的研究表明,尾缘锯齿可以降低叶片尾缘附近表面的压力脉动幅值约50%,将展向相关尺度较大的涡破碎成展向相关尺度较小的涡,并消除尾缘脱落涡,这三者的综合作用使噪声得到降低.      

剪切层与边界层组合流动的线性不稳定性分析

对不可压缩剪切层与边界层的组合流动完成了线性稳定性研究.组合流动的数学模型为Blasius边界层相似解与双曲正切函数的叠加,采用整体数值方法求解组合流动的稳定性方程,并验证了程序的准确性及网格无关性.研究给出组合流动的不稳定模态的辨识,即边界层模态和剪切层模态.在此基础上研究了剪切层对边界层模态不稳定性的影响以及壁面对剪切层模态的影响.由于剪切层的存在,使边界层模态中性曲线向左上方平移,临界雷诺数减小.此外,边界层模态不稳定性得到增强或抑制的影响,取决于扰动频率以及剪切层速度比的变化.组合流动中壁面边界层促使剪切层不稳定性减弱,主要表现在低频区域;而在高频区域,剪切层不稳定性几乎不受壁面边界层的影响.      

基于数值模拟的流场附面层边缘识别方法

附面层边缘通常取在速度达到主流速度0.99倍的位置,而复杂流场中主流流动往往并不均匀,给附面层边缘的准确识别造成了困难。为解决此问题,提出了用"参考主流"代替实际主流识别附面层边缘的方法:通过零剪切力滑移壁面边界条件下数值模拟得到不受附面层干扰的参考主流,在根据附面层定义确定附面层边缘时以该参考主流中的速度代替实际的主流速度。通过斜楔压缩和弯曲压缩两个超声速压缩流场对该识别方法进行了验证,所得到的斜楔压缩出口截面上附面层厚度与采用实际主流速度判断得到的厚度相对误差仅4.1%。根据该方法的识别结果对弯曲压缩型面设计进行附面层修正后,弯曲激波高度与无黏设计值之间的误差从修正前的2.0%降低至0.3%,压缩面末端压力的相对误差从修正前的6.6%降低至2.3%。该方法避免了指定主流速度的主观性,识别结果较为准确。     

基于大涡模拟的离心式喷嘴雾化过程研究

为实现离心式喷嘴雾化过程的精确数值仿真,探究喷嘴内部流动特性与外部液膜破碎形式,采用基于大涡模拟的仿真方法,对一种典型的四进口离心式喷嘴进行研究,仿真结果揭示了喷嘴内部相界面的振荡现象与外部液膜的破碎细节,并通过耦合流体体积法（VOF）与离散相模型（DPM）,获得液滴粒径的空间分布特征。研究结果表明：在液体填充过程中,喷嘴内的气液相界面存在波动与褶皱,形状并不稳定,内部的空气芯直径呈现正弦模式的振荡变化,喷嘴出口液膜厚度沿周向分布不均,这些因素导致出口附近的液膜表面出现扰动。在不同的进口条件下,不稳定性导致液膜表面上的扰动波形式不同。进口压力为0.3MPa时,液膜破碎由开尔文-亥姆霍兹（K-H）不稳定性产生的轴向正弦波所导致,产生沿周向分布的环形液带;在0.7MPa下,液膜表面开始出现由瑞利-泰勒（R-T）不稳定性引发的周向扰动波;随着压力增加至1.1MPa,液膜的破碎则由R-T不稳定性主导,产生沿轴向分布的液带结构,随后在气动力与表面张力的作用下破碎成液滴。二次雾化破碎后,喷嘴外部截面内的粒径呈“单谷”分布,液滴平均粒径计算结果与实验的最大相对误差为5.1%,与实验数据吻合度较高。     

高能合成煤油GN-1性能研究 --“液体火箭推进”专刊

结合火箭发动机工程应用环境,研究了GN-1煤油高温高压热物性、安全特性、传热结焦性能以及点火延迟性能,并和现役火箭煤油作了对比分析。采用商业仪器对GN-1煤油在最高200℃、最高压力25MPa范围内的密度、黏度、定压比热容、导热系数、表面张力进行了实验测量。在实验数据的基础上,依据基团贡献法、基于比容平移法则的P-R方程、摩擦理论、广义对应态原理分别对GN-1煤油在最高350℃、最高压力60MPa范围内的密度、黏度、定压比热容、导热系数、表面张力进行了理论计算,建立了GN-1煤油的密度方程、黏度方程、导热系数方程,并将方程的计算值与实验值进行了对比,计算偏差较小。对GN-1煤油和火箭煤油的安全性能进行对比研究,GN-1煤油的闪点为40℃,自燃温度为305℃,高于火箭煤油(225℃);燃点为47℃,低于火箭煤油(82℃),GN-1煤油的爆炸极限范围为0.44%-2.9%（40℃）。GN-1煤油和火箭煤油的急性经口毒性LD50＞5000mg/kg,均属于第五级化合物（实际无毒）。在入口压力10 MPa,流速10 m·s-1,内壁温480℃条件下,GN-1煤油的传热系数比火箭煤油提高14.4%,建立了传热准则方程。GN-1煤油出口油温220℃时试验段平均结焦速率是出口油温150℃时的4.43倍,GN-1煤油不锈钢材质管路中试验段平均结焦为高合管材质管路中的22.3%。在970K-1105K温度范围内,GN-1煤油的点火延迟时间为320μs-471μs,是火箭煤油的55.6%-69.3%。相关研究可对发动机可靠设计及应用提供重要参考。     

基于传热反问题方法的N_2O/C_2H_4预混推进剂燃烧室热流测量研究

测量液体火箭发动机的热载荷是获取燃烧室内部信息的重要方法。为了获取N_2O/C_2H_4预混推进剂燃烧室内壁的热载荷,建立了液体火箭发动机的热流计算的反问题方法,该方法基于对燃烧室壁面温度场的直接求解,通过对轴向多个位置测量温度的反演计算得到燃烧室内壁热流和温度。研究表明:应用文中建立的传热反问题方法能够较为准确地获得热流随时间及空间的分布;热电偶的位置对计算准确性有明显的影响,与理论深度偏差在0.2mm以内的随机深度偏差可导致超过4%热流反演误差;N_2O/C_2H_4预混推进剂燃烧室热流及温度沿轴向逐渐降低,表明燃烧室内的反应释热过程主要在燃烧室头部附近发生。