超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

JJ-5飞机失速/尾旋特性数值计算和图像显示研究

本文从飞机六自由度全量方程出发,计算、模拟了JJ-5飞机失速/尾旋特性和各种改出尾旋的方法,对JJ-5飞机特技机动飞行及其相应的尾旋特性进行了计算分析;并把计算结果和试飞结果进行了比较,结果表明数值计算与试飞基本上符合,结合计算机实时成像技术,对相应的计算结果进行了动态实时显示,显示结果直观、逼真,可望运用于有关尾旋飞行训练中.     

高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究

采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.     

真实飞行器无粘绕流直角切割网格数值模拟

本文采用直角切割网格，利用MGAERO无粘绕流解算器对三角翼真实飞行器带或不带有外挂物气动特性进行数值计算。计算马赫数为0．6～1．2。研究了网格密度对飞行器气动力特性的影响，同时与风洞实验结果进行了比较，符合程度令人满意。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

轴对称欠膨胀羽流流场的分离系数矩阵解法和实验研究

本文对轴对称超音速羽流流场进行了数值计算和实验研究。数值计算方法是空间步进的分离系数矩阵差分法,并用粘性和无粘分裂法处理粘性问题。同时对同轴超音速喷流进行了流场观察和测量,所得纹影照片清晰地显示了尾喷流场的详细结构。计算结果在近场与文献[1]一致,在远场与本文的实验结果基本一致。计算实践表明,与中心差分格式相比,分离系数矩阵方法的稳定性好,精度高,捕捉的激波耗散小,且无数值振荡。     

亚超音速全机纵向气动系数数值计算方法的改善与发展

本报告是在亚超音速位流解基础上改进发展的数值计算方法。用面分布有限基本解计算作用在物面上的压力,用数值积分方法计算作用在物体上的力与力矩。本文采用联合流场、有效剖面吸力概念以及顺流向片条作可压缩紊流附面层的积分解,来改进阻力特性的数值计算。 数值方法已编制成计算程序。 本文提供了用这个计算程序计算的若干算例,并与实验结果作了对比。比较表明,计算结果与实验值有良好的一致性。     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

可压缩流动转捩的迭代求解

本文叙述了高速可压缩边办流动稳定性和转捩点预报的数值计算方法，特别是Ｒａｙｌｅｉｇｈ反迭代法与边界层渐近匹配方法的配合，有效地提高了计算精度，节省了存储和计算时间，本文给出的评定界层第一模式的稳定性分析和转捩区数值预报结果，与已有的实验结果一致。     

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点.本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、PreV和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性.     

煤油超声速燃烧的数值分析

利用计算流体力学软件对煤油在所设计的双模态超声速燃烧室内的超声速喷雾燃烧进行了数值模拟。采用离散液滴模型、概率密度函数紊流扩散燃烧模型和紊流k ω模型计算了在飞行马赫数为5,煤油与空气的当量比为0 551时的情况。数值结果得到的壁面静压分布与实验数据基本一致;计算所得到的总压损失系数是0 696,非常接近实验测量值0 707,但计算得到的燃烧室燃烧效率远比实验值高。     

CFM56-5B型发动机风扇叶片鸟撞有限元仿真

针对CFM56-5B型发动机风扇叶片受鸟击问题,目前常采用实验和数值仿真两种方法,前者因成本条件客观问题难以实现,因此本文采用数值仿真来深入研究这一难题。本文运用CATIA三维软件进行建模,最大程度与实际风扇叶片相符,鸟体则采用目前应用广泛的椭球体流线型模型,运用LS-DYNA隐式求解方法对叶片进行离心应力初始化,建立了鸟撞叶片的数值计算模型,分析了鸟撞叶片的动态响应,通过大量仿真,得出鸟撞叶片的机械损伤程度与鸟的速度、受撞击位置之间的关系,为发动机风扇叶片的抗鸟撞设计提供帮助。     

复杂回旋流场的数值模拟和实验研究

本文通过四角喷射炉膛内回旋流场的数值计算和实验测试的对比研究，比较了几种计算方法的预估结果与实验结果的符合程序。结果表明：为克服当流动方向与网格线交角过大时所产和珠过分数值扩散，引入流向平衡粘性方法与κ－ε湍流模型相结合可获得与实验结果相一致的计算结果，而采用混合分格式加κ－ε模型所得到的计算结果比流向平衡粘性方法加等效常湍流粘性系数这样的简化处理的计算结果还要差。     

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

椭圆翼型高速气动特性实验研究

针对旋转机翼飞机主机翼翼型设计的特殊要求,在中国空气动力研究与发展中心FL-21风洞对16％相对厚度椭圆翼型马赫数0．4下的高速气动特性进行了实验研究,实验采用表面测压和尾排型阻测量技术。实验结果与数值计算结果的对比分析表明,数值计算具有较好的精确度,椭圆翼型的升力线斜率较小,压心位置靠前,力矩系数随迎角线性增大,最小阻力系数迎角并不像传统对称翼型那样为零度。     

二维非常化学非平衡流动高精度数值模拟及其与实验的比较

本文成功地将二阶迎风ＴＶＤ数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中，气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向。文中详细给出了气体力学方程和组元质量守恒方程组全耦合求解过程，最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论，本文所得结果与文献［４］提供的实验结果和文献［１］提供的数值结果相符合。     

S弯管中可压三元紊流控制方程和数值计算

本文通过对瞬态N-S方程的质量加权平均处理,对不可压缩流的k-ε两方程紊流模型做了可压缩性的修正。用有限差分法进行了数值计算。计算结果与实验数据吻合很好,这预示着数值方法可能会发展成为研究S弯管中高亚音速三元流动的一种手段。     

轴对称喷管内流场的加权余数有限元计算

本文采用非定常技术,用加权余数有限元方法给出喷管内流场的数值计算方法,对收敛一扩散喷管的实例计算表明:本文计算的壁面和中心线M数分布曲线与实验结果是一致的。