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

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

用多重网格方法计算旋翼跨声速无粘流场

发展了一种加快悬停旋翼无粘流场计算收敛速度的多重网格方法。由于悬停旋翼流场中存在不可压区域，同时旋翼尾涡系统的发展需要较长的时间，使得旋翼流场的收敛速度远低于固定翼流场，因此研究旋翼流场的多重网格算法具有重要意义。空间离散格式采用了中心有限体积方法，时间推进应用了五步龙格-库塔法。采用3层网格的V循环，对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明：尽管多重网格方法对旋翼流场的加速收敛作用不如对固定翼流场的加速收敛效果，但是多重网格方法仍然可以显著地加快旋翼流场收敛。     

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

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

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

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

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

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

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

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

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

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

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

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

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。