弹丸超声速绕流的Navier-Stokes模拟

本文轴对称Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组出发,采用高分辨率对称ＴＶＤ格式数据模拟了包括底部回流区在内的弹丸超声速绕流流场。     

高超声速钝体绕流的奇异摄动理论

对于二维和轴对称的物体的无粘高超声速绕流,当比热比Υ→1和来流Mach数M→∞时,激波将贴近物面,Chernyi等人应用薄激波层理论给出了尖头物体绕流问题的渐近展开解,它属于正则摄动理论。但对于钝头体,Chernyi等人的理论遇到了困难,因为该理论中对各物理量的量阶估计在物面附近是不适用的。本文在对物面附近各物理量的量阶重新进行估计的基础上,以Chernyi等人的渐近展开式作为外解,并在物面附近给出了各物理量的新的内解的渐近展开式,从而用奇异摄动理论中匹配的渐近展开法,得到了钝头体绕流问题的解。文中还给出了渐近展开式中首项的解析表达式结果。     

超声速平板绕流DSMC/EPSM混合算法研究

本文将平衡粒模拟方法（EPSM）与蒙特卡罗直接模拟方法（DSMC）结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于4时超声速平绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性。将DSMC方法与混合算法的计算效率进行了比较。     

超声速高超声速球锥绕流的边界层稳定性特点初探

本文对一个球锥组合体,在超声速来流及零迎角条件下,对其边界层的稳定性特点进行了分析并讨论了转捩预测问题。结果表明转捩总是最先出现在来流马赫数不太大,或来流雷诺数较高并且具有较大的周向速度或周向波数的情况下。     

具有气体动力和物理化学非均匀性的超声速物体绕流

本文研究的问题是超声速气流流过流体时产生的气体动力和物理化学非均匀性及其对物体空气动力特性的影响。     

细长圆锥低超声速绕流层迹流时空演化特性的研究

本文用数值方法研究了在大的迎角范围,采用局部抛物化流动模型所得到的有关圆锥有迎角绕流尾迹流中周期的谐分岔,亚谐分岔以及环上Ｈｏｐｆ分岔这三类不稳定性的时空结构演化与非线性相互作用。     

特型钝体绕流动态特性的低速实验研究

在低速自由射流风洞上 ,利用热线及相关设备对由半球头部、颈部及旋成体构成的特型钝体周围非定常流动的基本特性、联合特性及湍流度分布进行了测试与分析 ,实验工况为迎角α=- 1 5°～ 1 5°、侧滑角 β =0°～ 1 0°。结果表明 ,在上述工况范围内 ,该钝体周围流场中速度脉动的能量分布平坦 ,属宽频带随机信号 ,流场中没有发生明显的流动分离 ,流场动态品质良好 ;凹陷区内气流的脉动以不同的速度向下游空间传播 ;在上述α、β变化范围内 ,钝体颈部凹陷区的最大湍流度高达 1 0 .9～ 2 4.8%。     

用空间推进法求解钝锥三维粘性绕流

本文将文献〔１〕中提出的方法推广应用于三维可压缩Ｎ－Ｓ方程，方中给出了钝锥超声速三维粘性绕流的结果。得到了背风面分离的流动图画，计算结果与实验相吻合。     

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。     

带有分离的二维可压粘性喷管流动计算

采用二维、时变、可压纳维-斯托克斯(N-S)方程计算喷管流动。流动可为无粘、层流或湍流。湍流模型可采用代数混合长模型,一方程和二方程模型。喷管形状为二维(平面或轴对称)。计算网格内点采用MacCormack显式格式。采用局部细网格,快速解法以及时间间隔光滑方法,以加速高雷诺数下收敛。边界点采用参考平面特征线格式计算,粘性项作为源项处理。给出了初步算例结果,与其它计算和实验数据比较,符合良好。     

用隐式多重网格法计算三维粘性流动

求解高雷诺湍流流动时,边界层法向网格间距较流向和展向相比非常小,因此边界层中存在高度伸缩的网络,这将大大降低多重网格的求解效率。通过谐调的处理多重网格过程的各个细节,既提高了解的精度,又克服了网格展弦比的影响,使得计算效率提高了6～7倍。     

钝锥有攻角分离流动的数值模拟及其分析

本文定性分析了开式分离的性状,并对钝锥有攻角超声速绕流的开式分离作了数值模拟。分析指出,开式分离可能存在两种形态,第一种分离线的起点为正常点,第二种分离线的起始为鞍、结点(包括螺旋点)的组合。对于第一种形态,分离线的起点是横向分离的起始点,除分离线外,分离面上的流线不是从分离线的起点发出的。对文中计算的情况,流动属第一种开式分离。计算证实了定性分析的结论。计算和分析均指出,对第一种开式分离,在分离的起始区域,分离流面尚未卷曲,但在下游,则变成卷曲面。文中还研究了围绕物体的流管在分离诱导下的变形情况。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。