三阶HWCNS的构造及其在高超声速流动中的应用

对网格质量要求高、计算稳定性差和计算效率低是制约高阶精度格式应用于高超声速复杂流动模拟的重要因素。针对这些问题,发展了三阶精度的混合节点半节点加权紧致非线性格式(HWCNS3),改进其光滑测试因子和非线性权得到了HWCNS3-OP,并给出了它们的频谱特性。利用Lax和Osher-Shu算例测试了格式对间断和高频波的捕捉能力;通过钝锥和航天飞机的高超声速绕流算例,考察了HWCNS3-OP在真实流动模拟中热流和气动力的预测精度及其计算效率。研究结果表明:HWCNS3-OP具有较高的分辨率和良好的间断捕捉能力,高频波捕捉能力相对HWCNS3提高了约3倍,相对守恒律的单调迎风中心格式(MUSCL)提高了约4倍;HWCNS3-OP计算稳定性较好,计算效率相对五阶HWCNS提高了2~3倍,HWCNS3-OP是一种较适合高超声速复杂流动模拟的高阶精度格式。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

二维跨音速流动计算的显式多步有限面积方法

本文按Jameson等的思路,采用新的人工粘性项系数,建立显式多步有限面积方法,给出绕NACA0012翼型和RAE2822翼型跨音速流动的计算结果。     

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

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

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

钝锥超声速粘性绕流的隐式推进求解方法研究

本文提出了一种新的推进求解方法。通过对模型方程的稳定性分析和钝锥超声速绕流的数值模拟,表明该方法既克服了文献中推进格式步长有下界和稳定性差的不足,也克服了亚声速区压力不好处理的困难。文中建立的一阶、二阶隐式推进格式,稳定性好,推进步长无下界,计算结果的精度较高。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

旋涡内流计算

本文采用涡量-向量势表示的不可压缩三维流动的基本方程及其边界条件。针对小曲率弯管流,忽略主流方向的扩散项而使基本方程抛物化并利用流量守恒条件确定轴向压力分布。为避免横向速度向量散度不为零所带来的麻烦,本文提出两种计算方案。计算结果证明:用于求解三维抛物化方程组的这两种方案都是方便可行的。     

吸力面附面层抽吸在三维高负荷扩压叶栅中的作用机制

对某矩形高负荷扩压叶栅在不同弦向位置开设全叶高抽吸槽的5组方案进行了数值研究,分析了抽吸槽弦向位置等参数对抽吸量分布规律的影响;通过叶栅实验探究了局部展向抽吸方案的效果.数值仿真的计算域包含吸附叶片内部的真空腔,边界条件按照实验条件设置.研究发现:全叶高抽吸方案的抽吸量沿展向大致呈C型分布;叶高中部和端部的主要抽吸效果都体现在叶高中部流场,端部的抽吸量对叶栅角区的回流有一定的抑制效果.抽吸量沿展向的分布规律受叶栅流道和叶片内腔流场的共同作用,因此应根据三维高负荷扩压叶栅流场的具体特性对吸力面抽吸槽/孔进行细化设计.      

粘性三维流动计算的任意拉格朗日—欧拉法

推导了旋转叶轮内流三维粘性流动计算的ALE格式。以网格的虚拟运动计算对流作用;采用粘性子循环时步、拉格朗日子循环时步加速收敛;改进上风格式为部分施主单元格式来减小扩散同时实现良好的计算性能;由通量限制控制耗散截断误差;计算中引入短波长阻尼以消除不规则性;以亚网格(SGS)湍流模式模拟湍流运动。      

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

基于BVD原理的高保真空间重构方法

简要综述了一类基于单元边界变差最小化(Boundary Variation Diminishing,BVD)原理,设计双曲守恒律高保真数值格式的空间重构方法。BVD原理要求尽量减少通过重构得到的网格边界两侧物理量之间的差,从而能够有效地控制黎曼求解器中的数值黏性。BVD方法针对数值解的空间分布特征,选择多个函数作为空间重构的候补函数,并根据BVD判定准则从候补函数中选取最合适的函数进行空间重构。BVD判据不需要根据求解对象进行经验参数(阈值)的调整。选用适当的候补函数和BVD准则,可以完全避免现有算法中为抑制数值振荡而必须采用的非线性限制。BVD格式能在抑制数值振荡的同时,有效地控制数值耗散,可以对光滑解与间断解都获得高保真的计算结果。本文概述了BVD方法的基本思想、设计相关格式的基本思路,以及一些具有很强实用价值的BVD格式。并通过单相和两相可压缩流动的一些典型算例验证BVD格式的特点和优势。     

11组元轴对称热化学非平衡流场的数值研究

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,采用双温度、十一组元反应气体模型,利用隐式ＮＮＤ有限差分格式的时间预处理技术数值求解了ＦＩＲＥＩＩ飞船热化学非平衡流场。文中同时给出了七组元模型的计算结果,并进行了结果比较。     

绕振动旋成体的跨声速流数值分析

本文采用Ｅｕｌｅｒ方程和ＴＶＤ格式数值求解了小幅振动旋成体的非定常跨声速绕流问题，，假设在绕小幅振旋成体的非定常气流中，物理量的变化为平均定常扰动物理量的高阶小量，然后对Ｅｕｌｅｒ方程进行摄动展开，导出前两阶（零阶及一阶）摄动量的控制方程和相庆的边界条件，并进行数值求解。文中仅给出了振动旋成体的非定常压力分布。     

定常吹/吸气控制凸包分离的数值研究

针对凸包流动分离主动控制开展了参数化数值研究。凸包模型为NASA设计用于流动分离主动控制的内外流通用模型。重点研究运行参数变化时定常吹/吸气影响。首先采用试验数据标定了所开发程序,然后进行了一系列参数化数值模拟。结果显示:定常吸气作用明显不同于定常吹气。对于所给定几何模型,不管吸气流量增加到多大,吸气都无法消除分离,而吹气可以做到。详细展示流动结构并讨论了流动机理。      

特征无反射边界条件应用于二维圆柱绕流模拟

采用特征无反射边界条件和吸收层技术相结合的边界处理方法,研究了二维圆柱绕流模拟的边界反射问题。近年来在计算气动声学中提出了将大涡模拟和声学传播方程结合的混合方法求解声场,基于计算量的考虑,须在较小的计算区域内进行大涡模拟,这就需要在计算区域较小时,边界的反射尽可能小。作为大涡模拟研究的技术储备,将该方法应用于二维圆柱绕流的模拟中。数值结果表明这种方法在计算区域较小时,能起到很好的无反射效果,计算结果与前人在计算区域取得很大时的结果比较吻合。     

新型低能耗飞机液压系统低温试验方案研究

由于机载液压系统产品在整个飞行包线内会经历低温环境的考验,为了充分验证飞机用液压系统产品在低温环境下性能的可靠性,本文针对液压系统开展低温试验的多种设计方案进行对比,并对其中的三种低温试验方案进行分析,明确各自适应的试验类型,最后通过对比各种方案的优缺点,提出了一种可以满足长时间、低功耗、大流量的液压系统低温试验的全新方案。该方案通过实测满足了试验要求,并达到了预期高效节能的效果。     

计算空气动力学的流场和图形显示技术

本文阐述了计算空气动力学流场和图形显示的重要性,它的基本方法和主要内容;介绍了在VAX计算机上研制的计算空气动力学流场和图形显示软件的主要功能。