径向简谐激励下航空发动机机匣的振动分析

根据航空发动机机匣薄壁结构的几何特征,将机匣简化为一端简支一端固支的圆柱壳结构,建立了系统的动力学方程,并采用多尺度法对动力学方程进行了求解,进而研究了径向简谐激励下机匣动力学方程的零解及非零解的稳定性,并讨论了系统单模态运动发展到复合模态运动以及能量在模态之间的转移过程,最后采用数值方法研究了径向简谐激励下机匣的振动特性.结果表明:随着径向简谐激励的增加,系统表现为周期运动、倍周期运动及混沌运动交替出现.      

稀薄流到连续流的气体运动论统一算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数的各流域均适用的气体分子速度分布函数简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而在速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。发展可用于速度空间宏观取矩的离散速度数值积分方法，获取物理空间各点的流动参数，由此发展一套能有效模拟各流域三维绕流问题的气体运动论统一算法。研究气体运动论数值算法所适合的并行方案，基于统一算法的HPF并行实现，建立一套能有效模拟不同流域复杂外形体绕流的HPF并行算法软件。通过对不同Knudsen数的一维、二维、三维气体绕流问题进行数值计算表明，计算结果与有关实验数据及其它途径得到的研究结果吻合较好，证实了本文发展的统一算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。     

径向气体轴承-柔性转子耦合系统动力学研究

采用轨迹法对径向气体轴承(纯动压径向气体轴承、动静混合径向气体轴承)-柔性转子耦合系统的非线性动力学特性进行研究.建立柔性转子系统的多自由度模型,通过变方向隐式(ADI)方法实现瞬态气体润滑雷诺方程(含时间项)与转子动力学方程的耦合求解,通过数值仿真获得了系统在不同偏心质量作用下的非线性气膜力、轨迹图、相图、频谱图、能量谱图、分岔图、振型图及重力作用下轴承的平衡位置分布图.针对所采用的轨迹法的特点,研究了相应的非线性动力学参数获取方法.结果表明,轨迹法能够很好的描述径向气体轴承转子动力学特性及其气膜涡动现象,为径向气体轴承-柔性转子系统设计奠定基础.      

中心进气转静盘腔一维流动模型的改进

采用数值模拟方法,对动盘面边界层夹带流量进行分析,以改进转静盘腔一维流动模型的计算精度.结果表明:边界层速度分布的1/7幂律假设低估了自由盘边界层夹带流量,从而影响了一维流动模型对源区边界的计算精度;传统的关联式未能有效地归纳不同工况下转静盘腔内动盘面边界层夹带流量,从而影响了一维流动模型对核心区气体转速的计算.通过归纳数值模拟结果中的自由盘和转静盘腔内动盘面边界层夹带流量,对传统一维流动模型中源区边界及核心区气体转速的计算提出修正,并对修正后的一维流动模型进行实验验证,表明该修正明显改进了一维模型对腔内径向压差的计算精度.      

激光与空气作用过程中辐射能量损失的研究

激光与空气相互过程中,与辐射相关的各种过程起着特别重要的作用。研究了高温气体辐射对气体流体力学过程的影响。辐射能量密度和压力对气体状态和运动的影响完全可以忽略不计,而流场的辐射能量损失不能忽略。利用光线追踪法,运用辐射输运方程,根据光线在网格中的传播路径,计算空气吸收的激光能量;根据网格厚度,确定辐射能量损失。数值求解含能量源项的流体控制方程,对激光击穿空气的等离子流场演化过程进行了模拟,得到了辐射能量损失随入射激光能量的变化规律。     

SPH方法对气液两相流自由界面运动的追踪模拟

气液两相流动是自然界中常见的流动现象,对其进行数值模拟时要求必须能够准确跟踪界面运动变形。本文利用光滑粒子流体动力学方法,结合微可压缩模型(SCM),引入界面控制方法XSPH速度修正以及Van Derwaals状态方程修正,对典型的气液两相流动如二维溃坝、气泡上浮等问题进行了数值模拟,分析了空气和水相互作用机理以及界面运动规律,并同实验结果和其他数值计算结果进行了比较。结果表明,该方法在模拟多相介质界面运动问题准确有效,可用于处理更为复杂的多相流动工程问题。     

激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

非均匀来流稳态畸变流场的数值计算

本把TEACH计算机程序推广到三维，用定常，绝热，不可压的时均三维N－S方程和连续方程，并运用紊流动能和紊流耗散率的微分输运方程封闭雷诺方程，模拟计算了非均匀来流稳定畸变流场。数值计算中除了对通用变量进行欠松驰，还对修正压力通过惯量进行了欠松驰，惯量取为压力修正方程中的系数|b|，保证了数值计算的稳定性和收敛性。在进口马赫数小于0．35时，流场畸变指数的计算值和试验值符合很好。     

基于Newton-Raphson算法的径向轴承瞬态热流体动力润滑的研究

建立了径向轴承在载荷和速度突然变化时的三维数学模型，模型中考虑了轴瓦的热变形，在油膜和轴瓦交界面采用热流连续的理想边界条件，数值模拟轴承的瞬态温度场，并对轴承的瞬态性能进行分析。在每一瞬时，用Newton—Raphson算法同时求解Reynolds方程、膜厚方程和轴颈运动方程获得轴承油膜的压力分布和轴颈中心的运动速度，然后数值积分压力分布得到轴承的油膜力，差分运动速度得到轴颈中心位置和运动加速度。用一有效的有限差分法同时求解油膜和轴瓦的温度控制方程。最后将Reynolds方程和能量方程通过节点压力和温度相耦合获得轴承的瞬态三维温度场。结果表明本所介绍的方法收敛快，大大节约计算时间。     

航空发动机高压旋转涡轮盘腔流动与换热

在同模型实验参数范围内，且数值解和实验解基本吻合的前提下，对发动机原型参数进行数值模拟，求解了旋转系下流体流动与换热微分方程组，通过对气体域和固体域进行共扼计算，并将动量方程与能量方程福合求解。结果表明利用低参数下实验数据支持的数值计算程序对真实参数的模拟是合理的，对于复杂的非线性现象，数值外推法比传统实验外推更加准确。     

总压形式的能量方程与热阻问题

本文根据热力学第二定律的熵方程和稳定流的能量方程导出了用总压表示的一维定常流的能量方程。这是不同于稳定流能量方程和伯努利方程的另一种形式的能量方程。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

高亚声速翼型柔壁自适应壁风洞实验技术

本文通过解跨声速小振动方程求解假想外场，给出进行自适应壁风洞实验的迭代方案。结果表明：二者具有较好的一致性。本文研究表明：两种方案均有很好的应用前景。     

柔性自适应壁风洞的翼型实验技术

给出一种新的二维跨音速柔性自适应壁风洞实验迭代方案,计算方法和实验验证结果。根据实验时实验段上下壁和模型上下表面实测压强分布对风洞内、外流场进行非线性数值模拟;计算流线化壁面的形状,进行自适应实验。用该迭代方案,在堵塞比ε=8％,实验段高与翼型弦长比值H/c=1.5情况下,对NACA-0012翼型进行了高亚音速验证性实验。实验结果与国外大风洞无干扰实验结果吻合很好。实验时迭代次数仅需1～2次。实验结果展示了自适应壁风洞实验技术用于翼型跨音速实验的前景。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

二维超音速激波风洞及其应用

简述了二维超音速激波风洞的设计要点和性能,并给出在该风洞中低凸台诱导激波和湍流边界层相互作用的实验结果。实践表明：在被驱动段和喷管之间安装一个前缘光滑的矩形剖面短管道;并将二维喷管精加工成型,就能获得实验所需的均匀超音速气流。     

跨声速开槽壁超临界翼型洞壁干扰的N－S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一步，其结果就符合得更好，这为下一步用Ｎ－Ｓ方程进行洞壁干扰修正打下了基础。     

跨音速自适应风洞的任意初始壁方案

提出一种新的二维柔壁自适应风洞进行跨音速三维模型实验时,基于任意初始壁状态的洞壁调整方案,讨论了在较高跨音速实验Ma数的实验条件下,采用洞壁预调效应处理的洞壁自适应方案的实验结果。     

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

稀薄流到连续流的气体运动论统一数值算法初步研究

从非线性模型Ｂｏｌｔｚｍａｎｎ方程出发,引入简化速度分布函数、使用离散速度坐标法对速度空间进行离散、降维,去掉分布函数对速度分量的连续依赖性;采用时间分裂法,将简化速度分布函数松驰变化方程分解为源项碰撞变化方程、对流运动方程,进行耦合计算,应用ＮＮＤ耗散差分方法直接模拟气体分子速度分布函数;发展离散速度数值积分法,通过宏观取矩获取物理空间各点的流动参数,从而建立一套能有效模拟各流域气动问题的简化的