采用弯扭导向器改型涡轮级提高发动机性能的三维粘性流场的数值研究

对某型涡喷发动机涡轮级导向器进行了全三维弯扭设计以后,在保证涡轮前总温、总压和燃气流量的前提下,合理调整尾喷管面积,提高涡轮后背压,保证在设计转速下涡轮输出功与原型涡轮级相同,从而保证了压气机的设计点和喘振裕度,同时提高了发动机推力。带来这样效果的原因是弯扭导向器使涡轮级的效率提高,导致涡轮出口总压的提高。本文对改型和原型涡轮级进行了全三维粘性数值模拟和气动分析。      

气动静/动稳定性的相关性及粘性滞后时间的研究

风洞试验证实,物体的气动静稳定性和气动动稳定性之间,存在着反向演变现象。研究表明,它主要是绕流的粘性效应诱导的气动量和它们对物体姿态运动的时间滞后引起的,并且,该气动量中的静稳定性粘性诱导量和动稳定性诱导量之间是负相关的。本文给出了相关分析以及粘性诱导量之间的负相关模式,并用作者得到的风洞动态试验结果证明了上述负相关模式的实在性。此外,还根据上述风洞试验结果,计算了粘性动态相关因子Ｋｄ与粘性滞后的     

超声速粘性剪切流空间稳定性对称紧致差分数值分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Ｍ数、Ｒｅ数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快／慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。     

曲率和挠率对圆截面螺旋管道粘性流动的三阶作用

利用Ｇｅｒｍａｎｏ的曲线正交坐标系下的完全Ｎ－Ｓ方程,采用摄动方法求解了圆截面螺旋管道内的粘性流动,得到完全三阶摄动解,研究了曲率和挠率对二次流动和轴向速度的一阶、二阶和三阶作用。在曲率一阶项的作用下,截面上形成方面相反的两个二次涡,轴向速度最大值移向外侧;曲率的二阶项使二次涡的中心移向外侧;在曲率和挠率的高阶项（二阶和三阶）的共同作用下,靠近下半截面的涡增加,同时,二次涡的中心和轴向速度最大值移     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

超燃冲压发动机缩比燃烧室流场数值模拟

为了采用二维N-S方程来研究带有垂直喷射的三维问题,即拓展二维模型的适用范围,提出了一种针对喷射的简化模型,即采用源项加质来近似模拟喷射的质量添加。为了验证这种简化的可行性,针对文献中的实验条件进行了对比计算。结果表明,数值模拟与实验点符合得较好。运用以上简化模型,对放置在自由射流实验台上的缩比超燃冲压发动机燃烧室流场,采用五组元单步反应模型进行了数值模拟,得到了各主要气动参数及组分质量分数的分布,表明,根据数值模拟所得到的壁面静压值与实验壁面压力符合得较好。      

粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。     

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。