超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

二维凹槽超声速湍流流动数值模拟

本文采用时间、空间均为二阶精度的NND差分格式,以及Menterk ωSST二方程湍流模型,数值模拟了二维开式凹槽超声速粘性流动。用组分扩散方程计算了凹槽内粒子驻留时间。本文给出了M∞=3,长深比L/D=3,后壁倾斜角为90°和30°两种几何形状凹槽的计算结果,与相应的实验符合甚好。计算亦表明,该模型较高地估计了涡粘性,从而使流动更稳定。     

超声速开式空腔流动的数值模拟

研究了不同长深比、不同雷诺数条件下,超声速开式空腔流动特性的变化。算法采用空间上5阶的WENO格式,时间上具有TVD性质的3阶Runge-Kutta方法,湍流模型选用B－L模型。数值模拟获得了不同的时刻的注场图像,压力波、涡以及空腔结构之间的相互耦合作用,使得空腔内存在典型的周期性压力脉动,但不同长深比空腔的流动性存在较大差别。     

瞬时起动的空腔初期流动

采用有限体积方法求解非定常可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，研究了受限空腔在射流驱动下瞬时起动的空腔初期流动问题，数值计算结果给出了时间序列的速度矢量图和空腔出口流量随时间的变化，展示了射流、分离、剪切层和旋涡的相互作用及其演变过程。根据流动结构特征，将空腔的初期流动过程分为三个阶段：初始涡结构形成阶段、三涡结构形成阶段和两涡结构形成阶段。     

空腔流动特性风洞试验初步研究

通过测量动、静态压力，对两个不同比例的空腔平板模型进行了测压数据相关性试验研究，同时开展了空腔流动模拟技术，即相似准则研究。确定了在跨超声速条件下不同比例空腔模型测压数据的相依关系。结果表明：不同比例模型在Ma≤0．8，Ma≥2．0时，静态压力测量数据具有相关性；0．95≤Ma≤1．2时，静态压力测量数据相关性差．但呈现统一的规律，动态压力频率特性在全部试验Ma下具有相关性。     

二维超声速混合层流动稳定性的数值分析与并行计算

本文用DNS方法对两股空气的二维超声速混合流动(Mc=0.5)的稳定性问题进行了计算模拟.通过对采集的流动参数在时空两个方面进行的分析,发现流动失稳后出现了间歇结构.有限差分计算采用了并行计算技术,在微机并行环境下可高效率地完成大规模计算任务.     

部分平均Navier-Stokes方法模拟超声速斜面空腔流动

采用部分平均Navier-Stokes (PANS)方法对超声速斜面空腔流动进行了数值模拟研究,评估其对于超声速非定常湍流模拟的性能,并与雷诺平均N-S(RANS),脱体涡模拟(DES)的计算结果及实验数据进行对比.研究表明:模化湍动能比例全场可变的PANS方法预测的速度、壁面压力和摩擦力系数分布与DES的结果非常接近,都与实验值吻合得较好,优于RANS的计算结果;非定常雷诺平均N-S(URANS)计算的自由剪切层近乎为二维稳态的,而DES和PANS方法可以求解出更为丰富的流动结构;模化湍动能比例全场统一的PANS方法虽然可以求解出相比RANS更多尺度的流动结构,但在近壁区不能回归到RANS模型,预测的湍流边界层的速度型偏离对数律,后续的流动计算也偏离实验数据.      

带凹腔的超声速燃烧室燃烧流场数值模拟

对带凹腔结构的燃烧室二维甲烷燃烧流场进行了数值模拟.采用迎风3阶精度MUSCL格式求解二维含组分守恒N-S方程,湍流模型采用剪切修正的RNG k-e湍流模型,分别分析了凹腔不同的长深比和导流槽结构对燃料燃烧的影响;对喷甲烷燃烧工况进行了计算研究.结果表明:凹腔可以提高燃烧效率,却使总压恢复系数降低;凹腔的长深比越高,燃烧效率越高,总压恢复系数越低;在总压恢复系数较高的情况下,采用导流槽可进一步提高燃烧效率.     

三维绕空腔流的有限元数值模拟及流场拓扑结构

本文用一种高分辨率、无波动有限元法对三维空腔绕流进行了数值模拟。从计算结果中可以清晰地看到激波的移动，剪切层的波动，旋涡的变化和三维分离，其流场分辨率是令人满意的。     

凹槽非定常超音速流动的研究

用SST的两方程紊流模型,采用时间和空间上二阶的方法,求解雷诺平均N-S方程,得到开放型凹槽的流动特性的变化。数值模拟了飞行马赫数为5时的测试结果,得到不同时刻流场参数的变化。凹槽内流体质量的增加和减少,有助于油气的混合和燃烧,对于维持燃料的超音速燃烧是非常关键的,对超燃冲压发动机燃烧室的设计具有重要的价值。     

细长锥体有攻角绕流的非对称分岔及其涡结构的数值模拟

本文用ＮＮＤ格式计算了具有复杂涡结构的细长锥体有攻角绕流问题，着重研究了非对称分离的分岔现象以及分岔前后涡结构的演变，数值计算结果表明，大约在ａ／（２β）＝１．０时，Ｎ－Ｓ方程出现分岔解，同时还给出了分岔前后涡结构及物理最分布的特征。     

带双凹槽燃烧室中超声速流场的数值模拟

针对新近提出的双凹槽和预燃室结构的燃烧室概念,采用混合通量分离法及NND格式,对其中的冷态流场进行了数值模拟.结果表明:在没有喷流条件下,该结构是典型的开式流动,有利于混合,在工程应用范围内,可以用二维结果粗略地估算三维情况;但在有喷情况下,开式结构被打破,必须对三维流动作真实模拟.     

雾化燃料在超声速气流中横向喷射混合的数值模拟

本文用双流体模型对液态雾化燃料在直简形燃烧室中横向喷射问题进行了数值研究,气相用二阶迎负ＴＶＤ格式求解三维全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,液相用二阶预估、校正ＮＮＤ格式求解三维Ｅｕｌｅｒ方程。描述相间相互作用的常微分方程用预估、校正的二阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求解。分别考察了喷口处气相压力、相间速度滑移、液滴直径和燃烧室入口处预置的超声速流向涡对雾化燃料混合效果的影响。结果表明：气相较液相的扩     

外压式超声速进气道流场的数值模拟

采用CFD方法对某型外压式超声速进气道流场进行了二维数值计算,叙述了网格生成、边界条件确定方法,控制方程求解收敛技巧以及计算结果的分析.数值计算结果与理论相符较好.     

前掠翼布局中鸭翼气动影响的数值模拟

 采用三维Navier-Stokes方程和剪切应力输运（SST）湍流模型,就鸭翼不同位置和形状对前掠翼鸭式布局气动性能的影响进行数值模拟,并针对风洞试验方法难于分部件研究机翼、鸭翼以及机身各自气动特性的缺点,对布局升阻特性按部件分解研究并分析流动机理。研究结果表明：前掠翼鸭式布局气动性能（特别是在大迎角情况下）与鸭翼位置及其形状紧密相关,高位近距后掠式鸭翼可以与机翼产生更为有利的相互干扰,与无鸭翼布局相比最大升力系数提高约28.3%,最大升阻比提高约15.4%,大大地提高了前掠翼布局的纵向气动性能。该研究结果可为先进前掠翼布局飞机的预研和发展提供一定的理论参考。     

超声速自由剪切层流动的数值模拟和理论分析

本文从二维完全Ｎ－Ｓ方程出发，采用ＮＮＤ２Ｍ差分格式，计算了超声速自由剪切层流场，并且研究了上游扰动在剪切层内的传播规律。在数值模拟的基础上对流动涡结构进行了分析。     

当量比对单凹腔超燃燃烧室流场影响的三维数值模拟

 采用有限体积形式的上下三角分解（LU）隐式格式和k-g湍流模型,求解了多组分的N-S方程组。以北京航空航天大学宇航学院直连式超燃实验台的煤油燃烧实验为算例,开展了当量比φ＝0.32~0.97下5个工况的煤油超声速燃烧流动的三维并行数值模拟。计算结果表明：燃烧主要发生在煤油喷嘴所在的下壁面附近,流动出现亚声、超声的分层现象,高当量比工况下燃烧放热形成高的反压时,凹腔所在的燃烧室段出现激波串流动结构。计算得到的壁面压力分布曲线与实验结果符合良好。