叶轮机械失速颤振的全三维数值分析方法

给出了建立在全三维N-S方程求解基础上的叶轮机械失速颤振分析的数值方法。包括用定常N-S方程求解三维粘性失速流和用三维N-S振幅方程求解三维粘性失速振荡流。为了验证方法的准确性, 用压气机环形叶栅的失速流实验和振荡叶片振荡流场的实验作了对比, 同时应用本方法分析了某30万千瓦汽轮机末级叶片的失速颤振。      

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

前缘机动襟翼气动力设计计算研究

本文应用欧拉方程对翼身组合体加装前缘机动襟翼构型进行数值计算，通过跨音速气动特性的计算和分析，得出前缘机动襟翼对翼身组合体气动特性的影响规律，为某型机加装机动襟翼的改进设计提供数值计算依据。     

编队飞行卫星轨道的初步设计研究

使用两体相对运动的Hill方程进行小卫星编队飞行队形的初步设计，给出了编队飞行中各绕飞小卫星的轨道要素确定过程。以空间圆形编队为例，在未考虑摄动的情况下，通过仿真验证了设计的有效性。     

三维复杂外形横截面拓扑结构随迎角变化研究

通过数值求解抛物化的椭圆型方程生成复杂外形三维网格;采用交替方向隐式分解的隐式NND格式求解全N S方程模拟"类升力体"外形在高超声速下的大迎角流动;采用张涵信发展的流动拓扑结构理论分析给出了"类升力体"外形垂直于体轴的横截面流线随迎角变化的拓扑结构;此外,通过计算发现:大于20°迎角后,在部分横截面背风对称线上出现结构不稳定的鞍点相连现象。     

非结构网格中LU-SGS隐式算法的非平衡性影响

在非结构网格中用LU-SGS隐式算法求解欧拉方程时,两个近似分解因子的项数可能相等,为平衡;也可能不相等,为非平衡.采用对网格重新编号的方法很难达到我们希望的平衡.本文对非平衡性的影响进行了探讨.对二维问题,四边形单元的非结构网格,设计出平衡与非平衡的编号方式.先对标量模型方程分析LU-SGS隐式算法的增长因子,然后通过数值试验来验证这种非平衡性的影响.结果表明,尽管非平衡时也能达到收敛,但平衡却远优于非平衡的情况.     

考虑传质传热效应的翼型积冰计算

本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.     

多级轴流压气机非设计性能的数值预估

采用LU-SGS隐式解法以及叶排间掺混面模型数值模拟压气机三维粘性流动,研究了多级轴流压气机的设计与非设计状态性能的预估问题。对一台五级高压压气机进行了详细的流场数值计算,包括了设计转速的100%,90%及80%等多个转速及各种工作状态。结果表明,本文所建立的数值方法与计算机程序对多级轴流压气机各种非设计状态的计算都具有收敛容易和快速的优点,计算得到的压气机总性能和各分级性能与设计参数进行了对比,证明该方法和程序的预估结果是比较准确的。根据设计转速下的计算结果,预估了压气机设计点的喘振裕度。此外,还对不同状态下叶片进口气流角的变化进行了分析。      

基于B-P本构模型的涡轮盘应力分析

针对Bodner-Partom统一粘塑性本构方程,采用Euler前差显式积分方法,编制了用户子程序,把该本构方程结合进了大型通用有限元软件ANSYS,进行了一维和二维的算例验证,并针对某涡轮盘进行了三维应力分析。该程序可考虑材料变形的率相关、循环加载等情况,适用于航空发动机高温部件非线性结构分析。      

雷达方程中部分参数数学推导

用Blake算法计算雷达作用距离时,有部分参数是通过查阅相关的图表来求得的,如可见度因子D0、大气损耗La和天空噪声温度Ta'的值需要通过表来确定,因此在实际运用中发现其存在误差.提取出这3个参数的数学计算模型,并给出算例计算值,通过算例计算值与图表中所查找的值进行对比,精确度都在满意的范围之内.