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

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

裙后向台阶超声速绕流特性数值研究

本文研究了在中部带有裙体后向台阶的轴对称体的超声速绕流特性。采用空间推进法数值模拟超声速无粘流场，对于轴对称后向台阶的分离区采用Ｃｈａｐｍａｎ－Ｋｏｒｓｔ理论模型进行处理，计算结果给出了与实验值吻合的物面压力分布以及超声速流场的流动图画。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

三角翼大迎角绕流的数值模拟

采用FINETM/HEXA非结构网格流场计算软件包对三角翼飞行器的低速大迎角绕流进行了数值模拟.结果表明,计算值与实验数据符合较好,采用Agglomeration 方法能大幅度提高收敛速度,而网格自适应方法用总体不大的网格单元数可较精细地模拟流场.本文还分析了前缘分离涡破裂前后的流动现象和旋涡横截面流线图谱的变化规律.     

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

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

超声速绕凹角流动的压力分裂形式简化N－S方程数值解

本文从压力分裂形式的简化Ｎ－Ｓ方程出发，用Ｒｕｂｉｎ格式对超声速粘性流动进行了数值求解，给出了超声速绕凹角附着流及带有小分离气泡流动的结果，与其它数值方法所得结果相比，符合较好。     

超声速绕凹角流动的压力分裂形式简化N—S方程数值解

本文从压力分裂形式的简化Ｎ－Ｓ方程出发，用Ｒｕｂｉｎ格式对超声速粘性流动进行了数值求解，给出了超声速绕凹角附着流及带有小分离气泡流动的结果，与其它数值方法所得结果相比，符合较好。     

平板绕流的抛物化Navier－Stokes方程分析

对于绕有限长平板不可压粘性层流问题,循抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＰＮＳ）方程,进行ＰＬＫ坐标摄动分析,推得与Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＮＳ）方程结果完全一致的摩擦阻力规律。该规律在１５＜Ｒｅ＜１０＿４范围内成立,证实ＰＮＳ方程适应中等Ｒｅｙｎｏｌｄｓ数流动。     

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

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

小Knudsen数稀薄气体绕流流场分析

本文从Ｂｏｌｔｚｍａｎｎ方程５矩和１３矩方程出发，应用三层模型分析了滑流国稀薄气体绕流流场的性态，提出了一个计算这种流场的完整模型。理论分析表明，摩阻系数的小Ｋｕｎｄｓｅｎ数效应是高除小量，可忽略不计；但必须考虑热传导系数的小Ｋｎｕｄｓｅｎ数效应。基于奇异摄动法分析，数值求解了绕二维平板流动流场，壁面摩阻系数和热流率的计算结果与现有的实验和理论结果相一致。     

双燃式冲压发动机中超燃燃烧室冷态流场数值模拟

采用Ｍａｃ Ｃｏｒｍ ａｃｋ 矢通分裂格式, 求解雷诺平均的Ｎ－Ｓ方程组, 采用Ｂａｌｄｗ ｉｎ－Ｌｏ－ｍ ａｘ 代数湍流模型和混合长度模型, 模拟超燃燃烧室内两股平行超声速来流的冷态流场, 取得了与试验相一致的结果, 计算表明回流区和射流流动参数对超燃燃烧室流场内波系的变化有着较大的影响。     

蜂窝夹层结构脱粘缺陷电磁锤敲击检测模型及持力时间分析

为了定量描述蜂窝夹层结构敲击检测中敲击头持力时间与脱粘缺陷直径和面板特性的依赖关系,建立了敲击过程的力学模型,推导和分析了持力时间随缺陷直径和面板拉力的变化情况,并与试验数据做了拟合和比较。结果表明,在拟合的面板拉力取值下,模型计算值与试验值吻合较好。模型可用来计算持力时间的数值和分析持力时间的变化趋势。对于同一面板的脱粘缺陷,持力时间随缺陷直径的增大而增大。对于同种材料面板的相同直径的脱粘缺陷,持力时间随面板厚度的增大而减小。对于同厚度面板的相同直径的脱粘缺陷,复合材料碳环氧树脂面板对应的持力时间短于铝面板。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

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

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

基于非线性耦合本构关系的改进边界条件

非线性耦合本构关系（NCCR）模型是在Eu的广义流体动力学方程（GHE）基础上,通过绝热假设、Eu封闭和Myong简化推导出的关于非守恒量（黏性应力与热流）的非线性代数方程,有效拓展了线性的纳维-斯托克斯-傅里叶（NSF）本构模型在非平衡流动中的模拟能力,为快速准确模拟连续与稀薄耦合流动问题提供了强有力的理论工具。针对该模型开展滑移边界条件研究,结合努森层内物理量非线性分布的特点,提出一套在物面处与模型精度相一致的非线性修正滑移边界条件。在有限体积框架下,采用AUSMPW⁺格式和LU-SGS方法以及NCCR的完整耦合求解算法,对不同稀薄程度的高超声速单原子氩气圆柱绕流和平板绕流问题进行数值模拟。研究结果表明,基于NCCR模型的修正边界条件准确刻画出物面努森层内流动的非线性特点,有效提高了固壁滑移边界的精度。采用非线性修正边界的NCCR模型准确预测了连续流、滑移流和过渡流域的物面压力、摩阻与热流系数。