计算三维振动机翼上跨-亚音速压力分布的局部线化面元法

本文以小扰动位势流理论为基础,采用局部线化假设,建立了三维振动机翼的跨音速下洗积分方程;给出了跨-亚音速流核函数的计算方法;用推广的偶极子栅格法计算了矩形谐振动机翼上的非定常压力分布。结果表明,本方法可以较准确地计算非线性的跨-亚音速非定常流动,关于跨-超音速流动的数值解法,尚需作进一步的研究。     

非定常气动力计算中的偶极子格网法

从亚音速、跨音速和超音速三个方面叙述了非定常气动力计算中的偶极子格网法，且给出了其基本解法和计算实例。     

叶轮机平面叶栅的气动设计及绕流计算

通过全速势方程的数值积分，利用ＡＦ２格式，给出了一种用于叶轮机中叶栅型面的气动设计和绕流计算的便捷，对于稳态亚音速，跨音速及超音速流动均适用。与目前普及很广的。在超音区应用人工压缩性的方案相比，该算法有明显的优势。     

叶轮机平面叶栅的气动设计及统流计算

通过全速势方程的数值积分，利用ＡＦ２格式，给出了一种用于叶轮机中叶栅型面的气动设计和绕流计算的便捷方法，对于稳态亚音速、跨音速及超音速流动均适用。与目前普及很广的、在超音区应用人工压缩性的方案相比，该算法有明显的优势。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

计算超音速后掠翼定常和非定常载荷系数分布的核函数法

 本文所提供利用核函效法计算超音速后掠翼定常和非定常载荷系数的方法和公式,不仅适用于直前缘情况,而且也适用于具有弯折前缘的情况,稍加修改也可适用于翼身组合体情况,从而推广了超音速核函数法的应用范围。     

飞行器亚、超音速气动载荷数值计算

应用Greem函数法数值模拟飞行器的气动载荷。给出了求解亚、超音速定常和非定常载荷的统一方法。数值模拟中釆用双曲四边形代替四边形元素;采用联合流场概念改进诱阻计算;并用计算机绘图原理核查输入的几何参数是否正确。以矩形翼及双三角翼航天飞机(089B模型)为例进行了数值模拟,结果与试验数据符合较好。     

一种高速风洞三维洞壁干扰壁压信息修正法

给出以壁压信息法为基础的高速风洞三维模型试验洞壁干扰修正方法。该方法要求洞壁附近为亚音速流动,但允许模型附近出现超音速区及激波。方法可用于各种通气壁或实壁风洞的亚跨音速试验。最后对四个模型在三种试验段中的二十多种试验进行了修正,其结果和NASA非线性修正方法的结果吻合得很好。     

亚声速谐振升力面的点偶极子法

本文介绍计算亚声速谐振薄翼非定常气动载荷的点偶极子法(DPM)。网格划分方式与偶极子网格法(DLM)相类似。此法将网格上的载荷分布集中置于网格中央剖面的1/4弦点上,控制点取在3/4弦点上,无需作数值积分计算。很容易用于计算形状复杂机翼的非定常气动载荷。 文中给出了矩形翼和后掠翼定常和非定常气动载荷的数值结果,其精度数值与DLM的相当,计算效率提高一倍。     

机翼跨音速无粘绕流数值计算

本文介绍机翼跨音速定常无粘绕流的一种数值计算方法。选用精确速势方程作为问题的数学模型。在直角坐标系中,经过适当的坐标变换,先把后掠翼变成矩形翼,再把无限的物理空间变成有限的计算空间。在计算空间中,用有限差分格式(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用线松弛迭代求解。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

叶轮机械颤振稳定性工程预测方法在跨声风扇中的进一步探讨

通过求解雷诺平均Navier-Stokes(N-S)方程并采用影响系数法,对三个基本正交模态的三维振荡跨声风扇叶片绕流问题进行了研究,基于刚体运动假设和模态叠加法,得到了可用于叶轮机械设计阶段颤振稳定性评估的稳定性参数图.结果表明,跨声风扇内部,激波对非定常气动力的分布具有主导作用;在一定情况下振型对颤振稳定性有重要影响,应将其作为颤振稳定性设计的重要参数之一;跨声风扇和亚声低压涡轮的稳定性参数图对比分析表明,两者稳定性参数分布形式具有相似性.      

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

叶轮机械颤振稳定性工程预测方法在跨音风扇中的进一步探讨

通过求解雷诺平均NS方程并采用影响系数法,对三个基本正交模态的三维振荡跨音风扇叶片绕流问题进行了研究,并基于刚体运动假设和模态叠加法,得到了可用于叶轮机械设计阶段颤振稳定性评估的稳定性参数图。结果表明,跨音风扇内部,激波对非定常气动力的分布具有主导作用;通过对稳定性参数图的分析表明,在一定情况下振型对颤振稳定性有重要影响,应将其作为颤振稳定性设计的重要参数之一;跨音风扇和亚音低压涡轮的稳定性参数图对比分析表明,二者稳定性参数分布形式具有相似性。     

地面颤振模拟试验中的非定常气动力模拟

地面颤振模拟试验作为一种颤振研究的新方法,可以有效地弥补传统气动弹性试验的不足。对地面颤振模拟试验的主要难点,即非定常分布式气动力集成减缩加载的方法开展研究:基于亚声速偶极子格网法和活塞理论建立了亚声速以及超声速翼面的非定常气动力模型,通过曲面样条插值以及有理函数拟合获得了试验时域减缩气动力;提出以颤振关键模态的振型为优化目标,使用遗传算法搜寻气动力最优减缩位置的优化方法;建立了闭环系统的时域状态空间模型,使用颤振时域仿真结果与频域理论结果进行对比,对比发现二者误差可控制在3%以内。研究结果表明,该文提出的非定常气动力模拟方法可以很好地表征翼面非定常气动力分布特性,可以作为地面颤振模拟试验研究可靠的理论基础。     

高背压进气道中内外流耦合作用的大涡模拟

采用大涡模拟研究了出口堵塞比为50.8%的轴对称进气道流动,重点考察了内外流耦合作用下流动的非定常特性。采用国家数值风洞（NNW）工程仿真软件进行数值模拟,得到的壁面平均压力、瞬时压力分布与试验数据符合良好。分析表明：为匹配出口背压,进气道在喉道区域形成激波串结构,使内流道流场分为上游超声速区、中部激波串区以及下游亚声速区;在激波串区,剧烈的逆压梯度产生了分离激波、激波串、分离区及分离剪切层等复杂结构;伴随着激波串运动和边界层大尺度分离,进气道壁面压力出现宽频脉动特征。脉动压力的时空分布表明：内流道脉动压力以扰动波的形式传播,为此建立的声反馈模型能较好地预测亚声速区的主导频率。相关性分析表明：激波串运动受上下游流动耦合作用,其中,频率为St=0.7的运动主要受上游流动影响,频率为St=0.9的运动主要受下游压力扰动波影响。