超声速前机身及翼—身组合体的全速势方程数值计算

通过数值求解全速热方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某些一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

热空腔-喷气流的组合辐射

讨论了热空腔(喷气发动机的加力筒)和喷气流的组合辐射问题, 推演了组合辐射的辐射传输微分方程, 导出了有边界辐射条件下的解。讨论了由多个温度不同的壁面组成且腔内充满高温燃气的热空腔辐射传输方程解的具体形式, 导出了它向后半球的光谱辐射强度的表达式;同时还推导了热空腔-喷气流组合辐射的光谱辐射强度的表达式。最后给出算例, 算例结果表明, 热空腔-喷气流组合辐射的空间分布具有与实测分布较一致的梨形包络, 其光谱分布也与实测的光谱分布的一致性较好。      

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

给出了建立在全三维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冰翼型的气动特性进行了计算与分析.