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

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

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

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

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

利用微气囊扰动进行流动控制的数值模拟研究

探讨一种可以用来模拟“三角翼上分布微气囊，从而控制流动获得滚转力矩”的数值方法。为微型飞机的气动设计提供一种工具。研究的内容包括：考虑微气囊的三角翼网格生成、流场NS方程计算、微气囊不同布局对流动的扰动等。     

HTPB复合固体推进剂含损伤和老化本构研究

通过热氧老化实验、声发射损伤检测实验和蠕变柔量测试实验确定HTPB复合固体推进剂非线性本构关系，并与Schapery非线性本构方程进行比较，结果表明，含损伤老化模型较好地解释了HTPB复合固体推进剂的非线性特征，可以应用分析HTPB复合固体推进剂的力学特性。     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

基于欧拉方程的跨音速外形优化设计方法

将流场分析和优化方法耦合起来，形成一种基于飞机轴向截面积分布的设计方法，用以减小跨音速和超音速飞机的波阻，在本文方法中，流场用欧拉方程进行解算，以阻力参数构成目标函数，以多个截面面积（或当量旋成体半径）为约束，选择合适的优化变量，采用传统的梯度法进行优化设计，文中利用非精确搜索措施，避免了用传统一维搜索所需的大量机时，本算法所需机时少，优化效果明显，可用于飞机的初步设计。     

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

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

一种三维弹性机翼气动优化新方法

基于控制理论方法和有限元分析软件,通过将网格生成、流场计算、有限元分析、粘性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,研究与发展了计及静气动弹性影响的三维机翼气动优化设计方法,其中,雷诺平均Navier-Stokes方程为主控方程,气动载荷和结构静弹性变形量由气动/结构方程的耦合迭代求解得到,目标函数梯度信息由共轭方程数值求解得到。典型大展弦比机翼气动减阻设计结果表明：研究及发展的计及静气动弹性影响的三维机翼气动优化设计方法是有效的,能够有效考虑静气动弹性的影响。     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.