液氧/甲烷发动机变截面冷却通道传热数值研究

为提高液体火箭发动机推力室再生冷却通道的冷却效率,对液氧/甲烷发动机推力室变截面冷却通道的耦合传热进行数值模拟,探究了冷却通道的高宽比对跨临界甲烷的湍流流动和对流传热的影响。燃气-冷却通道-冷却剂的三维耦合计算采用一种改进的迭代耦合方法。研究结果表明:在冷却通道横截面积不变时,增大冷却通道高宽比可以降低喉部燃气侧壁面最高温度。冷却通道的高宽比越大,冷却剂压力损失越大。但过大的高宽比会导致压力损失急剧增大,且进一步降低喉部壁面最高温度的效果不明显。燃气侧壁面温度在变截面冷却通道的突扩突缩处出现局部下降,且下降幅度会随着高宽比的减小而增加。大高宽比冷却通道中,喉部侧壁面附近发生传热恶化的范围有限,主要在肋侧壁面附近的下半部分。研究结果为推力室变截面再生冷却通道的设计提供了参考。     

非线性数值模拟柱面气膜密封动态特性

针对气膜密封研究中密封气膜动态特性分析难点问题,采用非线性数值模拟法对柱面气膜密封系统的动态特性进行研究[D].建立了柱面气膜密封系统的动力学模型,给出了基于有限元数值方法的瞬态气体雷诺方程（含时间项）与系统动力学方程的耦合求解流程和方法,数值仿真给出了柱面密封环的中心轨迹图、确定了一定转速下的临界无量纲质量、不同密封环无量纲质量时等量位移扰动下系统的响应图,分析了密封环无量纲质量对系统稳定性的影响.结果表明:光滑柱面气膜密封系统存在临界无量纲质量,临界时系统受扰后处于半频涡动状态,随着密封环无量纲质量的增加,受扰后系统收敛时间变长,稳定性变差.      

航空拖曳诱饵系统的动态特性研究简

 为准确预测航空拖曳诱饵系统能否干扰成功,建立了系统的物理数学模型并对其动态特性进行了仿真研究。应用集中质量法,将柔性拖曳缆绳离散为一系列由阻尼弹簧连接的节点,建立了缆绳的动态模型;对诱饵进行受力分析,建立了诱饵的六自由度模型;提出了缆绳与诱饵的耦合条件,使模型更加精确。分别对诱饵释放过程中,以及释放完成后载机机动时系统的动态特性进行了仿真研究,给出并分析了缆绳的形状、张力和诱饵的姿态角等参数的变化规律。结果表明：为避免出现“鱼钩”现象,应尽可能减小释放诱饵的初速度与载机空速的夹角;应按梯形速度释放诱饵,以使缆绳中拉力的最大值较小。释放完成后,应控制载机最大飞行速度,以避免缆绳进入载机的高温尾喷流区;载机作盘旋时,缆绳在载机的圆形轨迹之外,且载机飞行速率一定时,角速度越大,缆绳向外趋势越大,越有利于避开载机的尾喷流区。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

超燃冲压发动机内外流场三维燃烧数值模拟

以某型煤油燃料超燃冲压发动机整机模型为研究对象,对发动机外流场的流动与发动机内的流动燃烧过程开展了耦合三维数值模拟研究.燃烧室壁面静压计算结果与实验数据吻合良好,平均相对误差为4.9%,验证了数值方法的有效性.结果表明:内流场的流动燃烧状况既受到外流场捕获空气流量的限制,又会反过来影响着外流场的激波结构;泄压孔对保持发动机中各部件流量匹配起到关键作用;由于进气道内激波与边界层的相互作用以及泄压孔诱发的斜激波的影响,流场在第1处燃油喷射孔附近呈现出燃烧流场不对称特征.      

气动阀门自激振动机理及动态稳定性

基于弹簧振子的扰动响应特性,提出了气动阀门气固耦合自激振动的稳定性机理.以单向阀为例,根据小扰动原理,构建系统的气固耦合动力学模型及稳定性模型,通过求解线性方程组的特征根,得到了单向阀系统的稳定工作区间及参数影响特性,稳定性分析模型的有效性由单向阀气动试验验证.对于单向阀,存在压力-流量的临界稳定性曲线.工作压力一定时,当单向阀工作流量小于相应压力下的临界稳定流量,单向阀处于不稳定状态,一个微小的扰动都将会导致阀芯周期性振荡;反之,单向阀工作稳定,阀芯将处于稳定开度.增加单向阀阀芯阻尼、弹簧刚度及减小阀门进口直径在一定程度上将有助于提高单向阀的工作稳定性.     

叶顶间隙对高压涡轮动叶应力影响的数值研究

针对叶顶间隙直接影响高压涡轮动叶气动特性及结构可靠性的问题,对不同叶顶间隙下的涡轮动叶进行了气热耦合数值计算.通过对流体域与固体域在物理边界上能量的传递进行耦合计算,揭示了泄漏涡流运动对叶片所受温度、压力分布的影响规律,分析了叶片热应力、气动力变化对结构强度的影响.结果表明:随着叶顶间隙高度的增加,泄漏流速增加,热应力变化明显,叶顶间隙高度在0.3%~3.0%变化过程中,叶片最大等效应力增大2.6%.所得结果可为叶顶间隙控制提供理论依据.     

固体火箭发动机转角涡脱落分类

应用三维大涡模拟(LES)方法对C1x模拟试验器进行了数值模拟,以研究转角涡脱落引起的声涡耦合现象.结果表明:与二维和轴对称LES相比,三维LES方法可以准确捕获转角涡脱落运动的三维特性及锁频现象.根据转角涡脱落是否在腔内形成声反馈循环过程可将其分为两类问题,第1类形成声反馈循环涡;第2类形成自耗散涡.燃烧室的台阶扩张比、台阶前后长径比是影响两类转角涡脱落现象的重要结构参数.     

机身减速板流动特性研究（英文）

现役高机动战斗机普遍采用机身减速板来减小飞行速度和转弯半径并提高机动能力。采用物面测压及空间流场测量相结合的实验方法,在机身减速板开度60°,机身迎角0°～70°条件下,研究了机身减速板铰链力矩随迎角的变化规律,分析了减速板迎风侧和背风侧的流动结构。研究结果表明：减速板铰链力矩按迎角可分为3个区域：常值区（α=0°～16°）,减速板铰链力矩基本不变,因为减速板迎风侧正压力逐渐减小,而背风侧负压力逐渐增加,两种相反的变化趋势相互抵消。非线性增长区（α=16°～32°）,减速板铰链力矩显著增加,因为减速板铰链力矩主要贡献区为背风侧,该迎角区内减速板背风侧存在一对不断增强的旋涡,背风侧负压力显著增加。在非线性衰减区（α=32°～70°）,减速板铰链力矩在迎角32°～36°范围内急剧减小,因为在迎角36°减速板背风侧旋涡流动变为速度较低的再附流动;减速板铰链力矩在迎角36°～44°范围内逐渐增加,因为该迎角区作用于减速板迎风侧的机身涡不断增强,导致减速板迎风侧正压力显著增加;减速板铰链力矩在迎角44°～70°范围内逐渐减小,因为该迎角区作用于减速板迎风侧的机身涡不断减弱直至破裂,导致减速板迎风侧正压力逐渐减小。     

电热除冰的热力耦合特性及其对冰层的影响研究

采用计算和实验相结合的方法,初步研究了电热除冰过程中的热力耦合特性及其对冰层的影响。在电加热条件下,耦合外部气动力载荷的作用,采用有限元方法计算了不同热流密度下表面冰层和蒙皮间界面法向和切向应力的分布,比较了加热/不加热条件下界面法向和切向应力分布的差别,研究了冰层最大主应力随热流密度的变化规律。研究发现,电加热条件下,在表面冰层融化前,热力耦合特性将造成冰层内部应力的显著增加,从而造成冰层局部区域的破裂,加速冰层的破坏。同时,设计的原理性实验结果验证了热力耦合特性对冰层的破坏影响。研究结果对于电热除冰理论和除冰技术的发展有现实意义。