带副翼三维机翼绕流的Euler方程解

提出了一种求解带副翼偏转三维机翼绕的Ｅｕｌｅｒ方程计算方法。采用对接分区网格民分区求解算法的结合,有铲地求解绕引外形的复杂流动,提出了一种满足通量守恒的内边界耦吉选用ＶａｎＬｅｅｒ分裂格式离散玩粘通量项,并构造了一种Ｌｉｍｉｔｅｒ函数以保证ＴＶＤ性质。数值算例表明本方法是求解带操纵纵面偏转的机翼绕流的有效方法。     

来流条件对叶片表面附面层转捩影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明:来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。      

Global Response of Martian Plasma Environment to an Interplanetary Structure: From Ena and Plasma Observations at Mars

As a part of the global plasma environment study of Mars and its response to the solar wind, we have analyzed a peculiar case of the subsolar energetic neutral atom (ENA) jet observed on June 7, 2004 by the Neutral Particle Detector (NPD) on board the Mars Express satellite. The “subsolar ENA jet” is generated by the interaction between the solar wind and the Martian exosphere, and is one of the most intense sources of ENA flux observed in the vicinity of Mars. On June 7, 2004 (orbit 485 of Mars Express), the NPD observed a very intense subsolar ENA jet, which then abruptly decreased within ∼10 sec followed by quasi-periodic (∼1 min) flux variations. Simultaneously, the plasma sensors detected a solar wind structure, which was most likely an interplanetary shock surface. The abrupt decrease of the ENA flux and the quasi-periodic flux variations can be understood in the framework of the global response of the Martian plasma obstacle to the interplanetary shock. The generation region of the subsolar ENA jet was pushed towards the planet by the interplanetary shock; and therefore, Mars Express went out of the ENA jet region. Associated global vibrations of the Martian plasma obstacle may have been the cause of the quasi-periodic flux variations of the ENA flux at the spacecraft location.     

飞行过载对燃烧室化学反应流场影响

为了研究飞行过载对固体火箭发动机燃烧室化学反应流场影响,以Liang模型模化铝滴燃烧,以有限化学反应速率模型模化湍流燃烧,对过载条件下发动机内流场进行了数值分析,数值结果与试验结果取得了趋势上的一致。研究表明,文中采用的数值计算方法可有效重现发动机热结构故障点;飞行过载改变了流场温度、粒子浓度、化学反应速率等参数分布;过载条件下燃烧室绝热结构表面铝滴积聚及剧烈的化学放热反应是导致其异常烧蚀的原因之一,铝滴局部积聚燃烧会导致温度场畸变;热结构设计必须与流动结构匹配。     

基于MEMS IMU的舵面运动测量仪

针对飞行器升降舵、副翼角运动的测试需求,提出了单轴捷联姿态系统的原理及算法,介绍了基于MEMS陀螺仪和加速度计的舵面运动测量仪的构成,并对测量仪的工作过程进行了仿真.仿真表明,测量仪的性能满足舵面运动测量的需要.     

直升机桨叶刚柔耦合特性及计算方法分析

为了研究桨叶刚柔耦合特性,采用中等变形梁模型处理桨叶弹性变形,分别采用小转角和有限转角两种方法处理桨叶绕铰的刚性运动。以桨叶扬起下坠碰撞问题为研究对象。计算分析表明:(1)弹性桨叶存在较明显的刚性运动和弹性变形耦合;(2)挥舞角较小时,两种处理方法计算结果一致,均与试验数据吻合较好;(3)挥舞角较大时,刚柔耦合增强,两种处理方法的计算值差别较为明显;(4)桨叶弯曲刚度对桨叶刚柔耦合特性有较明显影响,刚度越小,耦合越强;(5)桨叶如果比较柔软,小转角处理方法计算值相对有限转角处理方法,相位滞后、幅值增加均较为明显,计算桨叶动响应时需引起注意。      

旋转对内冷通道气膜孔流量系数的影响

为了获得涡轮叶片旋转时对气膜孔流量系数的影响,采用非结构化网格及标准k-ε紊流模型,求解三维N-S方程,对带90°肋和气膜孔出流的旋转矩形通道内的三维流场进行了数值模拟,气动参数变化范围是:通道入口雷诺数Re=60000,150000,罗斯贝数Ro=0,0.11,0.22,气膜孔总出流比为0.22和0.09。Ro≠0时,旋转效应对气膜孔流量系数有明显影响,通道顺时针旋转时,科氏力由上壁面指向带气膜孔的下壁面,引起流量系数增加;通道逆时针旋转时,情况相反。计算结果还表明,在通道内同一径向位置处的两个气膜孔的流量系数是不同的,通道顺时针旋转时,进入左侧孔的流体和气膜孔轴线的夹角小于进入右侧孔的流体和气膜孔轴线的夹角,导致左侧孔的流量系数大于右侧孔的流量系数;通道逆时针旋转时,情况相反。Ro=0时的计算结果与实验数据符合很好。     

顶部喷气对高速轴流压气机性能及流场的影响

喷气是提高压气机稳定工作裕度的有效措施之一。采用数值模拟方法研究了顶部喷气对某轴流压气机气动性能的影响。数值计算所获得实壁机匣总性能与其试验结果相吻合。喷气数值计算结果表明，转子上游约30％轴向弦长处的喷气推迟了失速的发生。通过详细地分析顶部喷气对压气机顶部区域流场结构的影响，揭示了其扩稳机理。     

主燃孔对旋流杯下游流场的影响

采用PDA（Phase Doppler Analyzer）测量系统对带双轴向旋流杯、主燃孔和冷却气流的模型燃烧室的主燃区流场（以下简称真实结构）进行了实验研究,测量了燃烧室内流场的特性,并与不带主燃孔和冷却气流的头部旋流杯（以下简称纯头部结构）的气流场进行了对比.发现真实结构燃烧室内的回流区不再是轴对称,主燃孔和二次气流的存在使回流区明显被压扁（最扁处的长轴约为短轴的1.6倍）;真实结构的回流区长度明显缩短,其长度L与旋流杯出口直径D的比值（L/D）约为1.3,而纯头部结构的L/D约为2.5.     

含损伤粘弹性本构及其在有限元分析中的实现

为研究固体发动机药柱的损伤粘弹性行为,针对一种耦合损伤的三维非线性粘弹性本构模型,采用积分算法对其进行了数值离散,导出了本构方程的增量有限元形式.针对该增量型本构方程对应的非线性方程组,建立了采用Newton-Raphson迭代法求解的应力更新方法,并基于Abaqus软件和Fortran语言开发了相应材料本构子程序,最后应用算例对所设计数值方法及程序进行了验证与分析.结果表明,本文有限元解与对应解析解吻合良好,保持了很高的精度,所设计方法及程序有效,可应用于实际型号发动机药柱损伤结构分析.