弯头钝锥的超音速无粘绕流数值计算方法

本文用差分法给出了弯头钝锥的超音速无粘绕流值计算方法。为解决物体体轴弯曲给推进计算带来的困难,本文通过Euler方程自变量的转换,将方程的求解面变为可以同推进方向倾斜的平面,即对Euler方程进行“斜推”或“平推”求解。由数值算例表明,本方法能得到较好的结果。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

小钝头锥无粘平衡、非平衡绕流流场的数值分析

本文给出钝头细长锥无粘、平衡和非平衡高超声速绕流流场的数值分析。对存在化学反应的空气,考虑了13种气体组分和22个化学反应,包括N_2、O_2和NO的振动松弛与离介过程之间的耦合反应。头部亚-跨声速区和后身高超声速区分别用直线法和有限差分法进行计算。在来流条件M_∞=20～26,Re_∞=(p_∞u_∞α)/μ_∞=1.5×10~3～5.5×10~6(α为头部半径)的参数范围内,给出了绕流流场和等离子体鞘层诸特性,如激波形状、电子密度和电子碰撞频率等。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

弯轴物体无粘超音速流场的数值计算

本文应用作者发展起来的“分离奇性”差分方法,计算了绕弯头锥的三维无粘超音速流场。计算的物型是具有过渡段的弯头锥模型。本文给出了三种物型的计算结果。     

二维非均匀射流遇障碍时波系的变化

本文研究了非均匀射流在二维管道内的钝头体绕流问题。所采用的Godunov格式适合于不规则区域上双曲守恒系统的流场计算,数值计算的结果给出该流场的主要特征和激波结构,对于认识这一物理现象和了解波的运动规律是有帮助的。     

高速列车绕流的数值模拟

给出了用低速Ｎ－Ｓ方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Ｔｒａｉｎ ３Ｄ，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

高超声速飞行器表面气动热和粘性摩擦力计算

本文给出了高超声速飞行器表面摩阻和传热系数（斯坦顿数）的计算结果，采用两种方法平面切面法亦即二维边界层近似法和工程方法计算了飞行器高超声速绕流的粘性效应，并对两种方法的计算结果作了仔细的比较，由文可见，对于在稠密大气层内，沿轨道运行头速度恒定的高超声速有翼飞行器，能够用本文所采用的两种方法计算其表面摩阻和热截荷，此二法可成功地应用于绕复杂形状物体的流动参数计算。     

绕流问题的区域收缩法及其数值解

本文提出了一种适于绕流问题的数值方法,它考虑了流场中不同区域内不同的物理现象,给出了外部边界条件的良好逼近。使用这种方法可以节省大量的计算机存贮单元和计算时间,也可以推广到三维绕流问题的数值计算。     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流,本文研究了垂直于物面的横截面上的流态结果。利用解析方法,给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律,从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场,计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

平面激波绕矩形障碍物后的绕射与反射

本文通过激波管二维光测模型实验和采用矩形网格的FLIC方法进行数值计算两种手段针对平面激波绕矩形障碍物后的绕射与反射的问题进行了分析和研究。文中给出了绕射波的波形以及平面激波绕过矩形障碍物后与地面相互作用而产生的从规则反射到马赫反射的波系图案。本文考虑了入射激波马赫数M_i及不同矩形障碍物高度H的影响,并将二维光测模型实验同数值计算的结果进行了比较。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

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

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

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

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

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。     

高超声速复杂外形热流的NND有限元数值计算方法研究

有限元法是一种加权积分的方法,这个特点和权函数的特性使得有限元法在计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文研究了高超声速复杂外形热流的NND有限元数值计算方法,并对高超声速球钝锥外形和带翼飞行器外形的绕流进行了数值模拟,算例结果令人满意.