高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

冷热两股气流在圆管内掺混的速度场、温度场与辐射通量场的数值方法

用数值方法求解冷热两股轴对称平行射流在圆管内掺混的流场、温度场和辐射通量场。计算采用一种通用性的数值方法,其特点是强烈地以物理本质为基础,而不仅仅是单纯的数学推导。这样做,能方便地分析和解释计算结果,使之更真实地接近于实际。 通过计算,分析了当引射系数改变时,掺混的流场,温度场和辐射通量场的变化规律,其结果与实验数据比较一致。     

三维激波／湍流边界层干扰产生的起始分离的预测

本文对由平板上直立尖劈产生的激波/湍流边界层干扰引起的起始分离进行了研究,提出了一种预测起始分离条件的方法。边界层内的分析采用了Johnston的三角形模型,边界层外主流分析则利用了Prandtl-Meyer函数。预测方法的结果表明与文献中的实验数据符合良好,而且比McCabe理论和Korkegi准则及Lu的半经验关系有更好的物理基础。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

高超音速后掠激波与边界层干扰流场特性

应用表面油流和液晶温度显示、表面热流率和压力测量四种测试技术,在Ma_∞=7.8、单位长度雷诺数尺Re_∞=3.5×10~7/m条件下,研究了30°迎角尖前缘翼面诱导激波和湍流边界层干扰流场特性。结果表明:表面流动是准锥型,在强相互作用下,流动发生了二次分离和再附,导致表面油流线的流向角、温度、热流率和压力分布均出现明显的凹坑。     

平板湍流摩阻系数的显式表达式

 湍流摩阻系数的计算是工程技术中经常遇到的问题。计算时所依据的公式多为半经验半理论的结果,如Karman公式,Van Driest公式和Prandtl公式等。这些公式是以隐函数给出的,给计算带来不便。但是不难发现,这类公式经过适当的变换,可以化为统一形式的超越代数方程,这种方程没有有限形式的解析解,本文给出了一种近似解析解。根据这种方法除了导出上述公式的显式表达式外,还对White公式作了修正。当然解超越代数方程最根本的方法是迭代法,本文对此也作了简单讨论。     

高超音速湍流分离激波结构不稳定性的实验研究

用铂膜电阻温度计测量了前向台阶诱导激波与湍流边界层相互作用流场中的表面热流率脉动。试验条件是:自由流马赫数为7.8,单位长度雷诺数为3.5×10~7米~(-1)。给出相互作用区平均热流率和脉动热流率分布。结果表明:在激波诱导的高超音速湍流分离流中,激波结构是不稳定的,产生一个间歇区域。在间歇区中,表面平均热流率由未扰动湍流边界层的热流信号和低频高幅热流脉动迭加而成,出现一个极大值。相互作用愈强,愈大,间歇区域愈长。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。