基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

某小型涡喷发动机二维数值仿真

整机二维仿真模型由带有粘性项的二维欧拉方程,燃烧模型和损失模型构成.利用任意曲线坐标系来适应复杂几何边界计算,采用隐式高阶精度Godunov格式求解非定常欧拉方程,能够自动捕捉激波,数值稳定性高.计算的总性能参数与设计值对比显示:单部件的仿真精度高于整机,涡轮仿真精度高于压气机.子午平面上的计算结果分析表明:离心压气机出口处的分离流动、燃烧室掺混孔的简化对总性能参数计算结果影响较大.     

直升机旋翼非定常气动力及结构动力响应计算

以二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为模型方程,采用Ｊａｍｅｓｏｎ有限体积并引入变系数隐式残值光顺及固接网格计算了直升机旋翼不同展向位置处翼剖面的非定常气动力,进而采用插值的方法模拟直升机在悬停及前飞状态下,旋翼的桨叶所受的非常气动力载荷,并将之与桨叶的结构振动方程相耦合,求解其时域内的结构响应,非定常下气动力计算结果与实验结果做了对比,两者比较吻合,旋翼结构响应计算结果合理。     

跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

两种跨声速气动弹性问题分析研究

在非结构运动网格基础上,采用中心有限体积法进行空间离散和双时间方法进行时间推进求解非定常欧拉方程.通过与气动力方程的联立求解,在时域内用四步龙格-库塔方法求解结构运动方程.分析和研究了二维嗡鸣和三维机翼颤振这两种跨声速非线性气动弹性问题.二维嗡鸣问题的研究考虑了翼面-舵面系统的缝隙间网格运动、缝隙对嗡鸣的影响和扰流片对嗡鸣的抑制.耦合多自由度Lagrange结构运动方程数值模拟了三维机翼的颤振问题.通过跨声速标模算例AGARD445.6机翼的颤振计算,计算的颤振临界速度与实验值有5%左右的误差,验证本方法的正确性.由于本方法是在对外形具有良好普适性的非结构动网格基础上完成的,具有良好的工程实践价值.     

大小叶盘结构连续参数模型和振动模态

给出了用于研究大小叶片整体叶盘结构固有振动特性分析的连续参数模型。模型用 Timoshinko梁模拟叶片 ,用平板模拟轮盘 ,并在叶片和轮盘间引入弹簧元件模拟叶片和轮盘间的耦合关系 ,简化了两者间的边界条件和连续性条件。采用 Gram-Schmidt方法生成正交多项式作为李兹容许函数进行模型的离散化 ,简化了所得到的频率方程求解过程 ,有利于求解过程的数值稳定性 ,并可方便地通过对方程系数性质的分析了解结构振动模态的性质。利用上述的分析原理和求解方法 ,进行了简化的大小叶片整体叶盘结构的计算 ,计算结果与有限元计算结果具有很好的一致性      

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

翼型近尾迹流中大尺度涡结构的理论模型

根据流动稳定性理论相干结构的相关理论 ,提出了一种二维波模型来描述翼型近尾迹流动中的大尺度涡结构 ,通过采用对近尾迹流动区域分区处理的方法计算了近尾迹区域的大尺度涡结构。计算所得的流线图 ,涡量分布以及雷诺应力分布具有与实验结果相似的分布规律。说明该二维模型能够较好的描述近尾迹流动中这种大尺度涡结构的运动学和动力学特性      

超音进气道的亚音扩压器流场的数值计算和实验研究

本文对二维、亚音扩压管道中的粘性流动进行了数值计算,并对同样的扩压器模型进行了风洞吹风实验。 根据要计算的二维、亚音、壁面扩张角不太大的管道中的粘流,我们对N-S方程进行了抛物化假设,化成可沿流程推进求解的抛物型方程,再用有限差分进行计算。 经与实验结果对比,证明本文采用的结合零方程紊流模型,在几何流线坐标系下求解抛物型方程的方法对解决二维、亚音扩压管中分离区以前的流动问题,可得到满意的结果。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

飞行器刚体运动与气动弹性相互作用问题研究

通过联立求解空气动力学基本方程、飞行动力学运动方程和弹性结构振动方程,在时间域内模拟和分析了大展弦比飞机纵向动力学稳定性问题。结合动网格技术,气动力计算采用基于欧拉方程的计算流体力学方法,结构变形和飞行姿态位置变化统一为模态表示方法,通过松耦合将飞行器姿态稳定性和结构变形稳定性施行了模拟。以某大展弦比机翼飞机为算例,研究了其刚体运动和机翼的弹性振动的相互影响。结果表明：对于具有大展弦比机翼的飞机,其机翼的低阶弹性模态易与飞机飞行中本身的刚体模态发生耦合,从而导致飞机机翼的气动弹性发散以及飞机本身刚体运动稳定性的改变。对于这类飞机,在其气动弹性和飞行稳定性的分析和设计中必须充分考虑到两种运动的相互影响。     

进气畸变对轴流压气机性能影响的准三维计算

为了提高分析进气畸变对轴流压气机气动性能和稳定性影响的能力,发展了一种沿叶高将压气机分割为若干个子压气机、采用二维可压缩非定常的Euler方程描述各子压气机回转面上无叶区的流动,并利用二维可压缩"激盘"取代各子压气机的转子和静子叶排的预测进气畸变影响的准三维模型。稳定边界的确定采用Koch的最大失速静压升系数法,并引用了临界畸变角的概念。计算分析了总温畸变、总压畸变以及总温 总压组合畸变沿多级轴流压气机的传递和对压气机性能和稳定性的影响。计算结果与现有的实验数据十分一致。     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

多级压气机非定常流动的数值模拟

进行了多级压气机动静叶片排干扰的三维非定常无粘流动数值模拟研究，采用高分辨率格式求解非定常Euler方程组。对某以级压气机的非定常无粘流场的计算结果表明，即使在不考虑粘性影响的情况下，压气机动静叶排干扰的非定常效应仍然是比较显著的，各叶片排所受到的非定常影响是不相同的，计算得到的叶片表面压力周期性波动的结果可为叶片强度计算提供较为准确的数据。     

叶尖泄漏流对压气机稳定性影响的数值分析

为研究叶尖泄漏流对稳定性的影响,发展了一种叶尖泄漏涡模型,并且在由课题组开发的TUSIAC(three dimensional and unsteady stall inception analysis code)程序中实现。该程序将转/静子叶排模化为三维激盘,并在无叶区求解三维非定常欧拉方程,黏性的影响通过特性曲线体现,因此对计算资源的要求较低。数值模拟的结果表明:在叶尖间隙的影响下,压气机的性能和稳定裕度均有所下降。叶尖间隙不改变失速先兆的类型,但使得失速团旋转速度加快,周向尺寸减小。所发展的叶尖涡模型能够预测叶尖间隙尺寸对压气机稳定边界与失速起始过程的影响,从而为在设计初期考虑压气机的气动稳定性、并且优化压气机设计,提供了一种较为实用的方法。      

二维非定常欧拉方程的摄动解在超声速矩形翼上的应用

本文利用超声速线化理论给出的矩形翼载荷系数解析解作为权函数,对沉浮和俯仰谐振的二维平板翼压强系数进行加权修正,获得矩形翼谐振气动力的近似解,计算结果与线化理论数值结果符合得很好,本方法是超声速和高超声速矩形翼非定常气动力的一种统一和快速的计算方法。     

基于非结构运动对接网格的旋翼前飞流场数值模拟

发展了一种基于线性弹簧理论的非结构运动对接网格方法,通过求解非定常Euler方程模拟了旋翼的前飞流场。利用运动区网格的转动模拟旋翼的旋转运动,并运用线性弹簧动态网格技术处理旋翼的挥舞和变距运动,方程求解中采用一种精度较高的插值和重构方法。算例表明,本文的方法能够有效地模拟旋翼的前飞流场。     

跨音速颤振计算方法研究

采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。