有无化学反应的射流及尾流的稳定性(英文)

射流或尾流可视为无固壁边界而流场速度剖面具有拐点的流动。在无粘不可压流中,有拐点的速度剖面总是不稳定的,因此,射流及尾流总比固壁边界的边界层更不稳定。对于无粘可压缩流,情况有所不同。文献[1,2]曾指出,高M数流动条件下,射流或涡层可能是稳定的。最近,在对超音速燃烧问题的研究中,考察一下化学反应对可压缩射流或尾流稳定性的影响,也许是有意义的。本文评述了有无化学反应的射流或尾流的稳定性问题,发现化学反应对射流或尾流的稳定性有显著的影响。     

定常可压缩无粘流非线性方程解的完全边界积分表示式

 本文以ρV为因变量将可压缩等熵无粘流非线性方程的解表示成完全边界积分形式,从而将问题的计算维数降低一维。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

求解可压缩流动的同位网格SIMPLE方法研究

在Rhie-Chou动量插值的基础上,推导了同位网格可压缩SIMPLE算法.经过无粘流超音速凸包算例和激波/湍流边界层干扰算例计算发现,如果对流项采用高阶有界HLPA格式,密度插值采用一阶迎风和中心差分的混合格式,这种算法能够很好地模拟凸包超音流的流动现象,在采用了新型GAO-YONG湍流模型后也能够较好地模拟激波/湍流边界层干扰.      

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

基于PSE的单股剪切混合流稳定性分析

将抛物化稳定性方程(PSE)方法应用到可压缩单股剪切混合流的稳定性研究中.采用并发展了适用于自由剪切流的高精度数值方法,包括六阶紧致格式、坐标变换以及渐近边界条件等,对PSE进行有效求解.通过求解相似边界层方程得到更准确的剪切层内基本流;求解线性稳定性理论(LST)控制方程得到扰动的初始条件,并通过流向空间推进方法对扰...     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

自适应笛卡尔网格Ghost Cell方法研究

在笛卡尔网格中,利用Ghost Cell方法处理浸入边界,模拟二维无粘可压缩流。针对静止物体,比较了各种不同的Ghost Cell边界条件下的熵误差,总压误差,以及阻力系数。此外为提高激波分辨率,将Ghost Cell方法与基于叉树结构的自适应笛卡尔网格算法相结合,数值结果显示本文的方法是切实可行的。     

壁判据用于计算流体力学（CFD）可信度评估

本文把作者提出的近壁干扰剪切流动(ISF)全域理论与流体运动方程组及流体在壁面上无滑移条件相结合导出一组壁面判据.壁判据为计算流体力学(CFD)仿真的可信度评估提供了基于流体理论的一条直接验证途径.对不可压缩流动的十一个熟知的NS方程组精确解,包括二维驻点和斜入射再附点流,二维分离点和背风驻点流,轴对称驻点和背风驻点流,旋转圆盘附近的三维Von Karman 流、收缩和扩张渠道流和非定常斜入射三维驻点流;以及经典边界层及其无粘外流和相似性边界层及其粘外流的NS方程组解,提出可用于验证近壁流动计算的几个壁相关函数.证实它们准确满足所有壁判据;说明壁判据可用来检验NS方程组数值解近壁计算结果的计算精度并验证其可信度.     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

Investigation of the stability of inviscid compressible swirling flows

The stability of an inviscid compressible swirling flow between two concentric cylinders is analytically investigated. Two stability criteria are derived for compressible swirling flows under narrow gap approximation by an analytic method analogous to Ludwieg's method. Finally, the effect of compressibility is discussed.     

有限长圆管内刚体旋转流的直接数值模拟

通过对三维不可压缩粘性流体的直接数值模拟,研究了有限长圆管内刚体旋转流在初始扰动下,扰动的增长和动力学演化过程.首先,通过改变流动的Reynolds数和旋转角速度ω,分析了轴对称流动中Reynolds数和ω对流动失稳的影响.对于高Reynolds数轴对称流动,利用数值模拟得到的不稳定结果与无粘线性稳定性理论预测的结果一致.最后,为进一步研究流动不稳定发展的三维效应,又进行了完全三维流动的模拟.数值结果表明,三维流动中不稳定的主导模态为螺旋模态,模态在线性段指数增长,导致流动产生螺旋型的漩涡破裂,之后经过非线性段的增长后达到饱和,流动最终发展为轴对称泡型旋涡破裂的稳定状态.     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上，提出了一种探测并识别二维激波干扰模式的新算法，从3个层面详细介绍了该算法的实施流程。首先，采用基于当地流场参数设计的传统激波探测方法，辨识出激波附近的一系列网格单元；其次，通过经典的K-means聚类算法将这些激波单元划分成许多簇，并根据簇的相邻信息定义每个簇的类别；最后，设定相关准则对某些紧邻的簇进行合并，进而确定各个激波干扰点的位置，记录各条激波分支所对应的簇，采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明，该算法不受网格类型的限制，不仅可以保证最终拟合的激波线具有较高的位置精度，还可以清晰地识别出流场中多激波干扰的模式，同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

用隐式近似因子分解法计算喷管跨音速流场

本文用隐式近似因子分解法结合贴体曲线坐标计算了定常、无粘、跨音速喷管流场,选择了参数变化比较剧烈的小喉部曲率半径喷管,陡壁收敛段喷管和潜入喷管三种算例.计算结果表明,该方法收敛快,并且具有良好的精度.Courant数可以取至7.这种隐式法与显式法相比较,花费的机时少得多.本文的方法可望推广到求解N-S方程、PNS方程,和二相有粘喷管流动计算.