五阶HWCNS格式在低速复杂流场中的应用

采用五阶精度显式HWCNS格式求解雷诺平均NS方程，利用多块对接结构网格技术，对30P-30N多段翼型进行网格收敛性研究，在不考虑转捩的情况下采用SA一方程湍流模型研究HWCNS格式与二阶精度MUSCL格式对该翼型压力分布和典型站位速度型的影响，并与实验结果进行对比分析。采用HWCNS格式和SA一方程湍流模型模拟梯形翼高升力构型低速复杂流场，通过对总体气动特性和压力分布的分析，探讨五阶精度显式HWCNS格式在低速复杂外形流动中的应用能力。     

五阶HWCNS在低速复杂流场中的应用

采用五阶精度显式混合加权紧致非线性格式求解雷诺平均NS方程;利用多块对接结构网格技术,对30P-30N多段翼型进行网格收敛性研究。在不考虑转捩的情况下,采用SA一方程湍流模型研究混合加权紧致非线性格式与二阶精度MUSCL格式对该翼型压力分布和典型站位速度型的影响,并与实验结果进行对比分析。采用混合加权紧致非线性格式和SA一方程湍流模型模拟梯形翼高升力构型低速复杂流场,通过对总体气动特性和压力分布的分析,探讨五阶精度显式混合加权紧致非线性格式在低速复杂外形流动中的应用能力。结果表明,对30P-30N三段翼型,采用全湍流模拟方法可以得到较好的压力分布;对梯形翼高升力构型,在附着流和边界层小分离情况下混合加权紧致非线性格式有较好的模拟能力。     

WCNS格式在热化学非平衡流数值模拟中的应用研究

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。首次将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算。进行了自由流速度5.0~5.7km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到了正确的流场参数分布,壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

加权紧致非线性格式在热化学非平衡流数值模拟中的应用

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算,进行自由流速度5.0~5.7 km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到正确的流场参数分布、壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明,WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

固定舵二维弹道修正弹网格模型研究北大核心

为确定固定舵二维弹道修正弹的网格划分方法,分别建立非结构网格和结构网格模型。应用Fluent软件对所得网格模型进行数值计算,获得修正弹在不同马赫数下的气动参数。对所得气动参数与风洞试验结果进行对比分析,确定不同马赫数下更准确的网格模型。分析结果表明:2种网格模型均可应用于修正弹的气动数值计算,但在准确性方面,结构网格模型更适用于跨声速条件下的数值计算,非结构网格模型更适用于高声速条件下的数值计算。在不同条件下,将2种模型结合起来进行数值计算,会得到更加准确的气动参数。     

外挂式导弹机弹分离气动干扰特性研究

以经典的机翼/挂架/导弹组合模型为例,采用重叠网格软件系统和计算流体力学技术,从机弹干扰工况简化和气动干扰特性影响因素分析两个方面对外挂式导弹机弹分离气动干扰特性进行研究。导弹分离轨迹参数和气动干扰系数的数值预示结果与捕获轨迹试验结果吻合,表明该计算方法能有效预测机弹分离轨迹和分析导弹与载机间复杂气动干扰现象。根据计算流体力学结果,从马赫数、机翼攻角、导弹攻角等方面,给出导弹在不同分离工况下的气动干扰规律,并采用增量系数法对缺失工况进行一阶外插处理的气动干扰数据外推方法,可应用于机载外挂空基武器的机弹分离轨迹预示和气动干扰特性设计中,具有重要的工程应用价值。     

固定舵二维弹道修正弹网格模型研究

为确定固定舵二维弹道修正弹的网格划分方法,分别建立非结构网格和结构网格模型。应用Fluent软件对所得网格模型进行数值计算,获得修正弹在不同马赫数下的气动参数。对所得气动参数与风洞试验结果进行对比分析,确定不同马赫数下更准确的网格模型。分析结果表明:2种网格模型均可应用于修正弹的气动数值计算,但在准确性方面,结构网格模型更适用于跨声速条件下的数值计算,非结构网格模型更适用于高声速条件下的数值计算。在不同条件下,将2种模型结合起来进行数值计算,会得到更加准确的气动参数。     

用多重网格TVD方法求解三维高超音速粘性流场

本文采用时间相关法求解了三维Ｎ－Ｓ方程。在计算中，对于压力和对流项采用Ｒｏｅ的通量差分分裂技术和Ｏｓｈｅｒ－Ｃｈａｋｒａｖａｒｔｈｙ的ＴＶＤ格式，对于粘性项，采用了中心差分格式。通过数值模拟得到了钝头双锥体高超音速层流流场在不同攻角下的数值解。计算的重点放在壁面热流率的预测。在攻角为４°和１２°的状态下，Ｓｔａｎｔｏｎ数的计算值与实验值进行了比较，二者吻合较好。为了加快收敛速率，计算时采用了多重网格方法。     

高超声速飞行器前缘流-热-固一体化计算

针对高超声速流动气动加热与结构传热的复杂耦合问题,探索和研究基于有限体积法的高超声速流-热-固一体化求解方法,将流场与结构温度场进行统一建模与数值模拟。该方法避开了传统气动加热和结构传热耦合求解方法在时间域内进行流场与结构温度场耦合交替迭代计算所带来的大量数据交换与计算,将流场与结构温度场作为一个物理场,采用统一的控制方程进行求解。采用典型高超声速绕流二维圆管稳态或非稳态流-热-固耦合算例对该一体化方法进行验证,稳态时圆管驻点温度最高达到648 K,非稳态时的热流密度和结构温度与参考文献和实验值吻合较好,由此证明了该方法的可靠性和正确性。与耦合计算方法的对比分析结果表明:该一体化求解方法所得计算结果更接近实验值,并且计算量和网格依赖性都相对较小,具有更好的稳定性和计算精度,能为高超声速飞行器一体化热防护设计提供有效的理论和技术支撑。     

子母弹结构特征对分离特性影响分析

超声速下子母弹的分离过程是一种包括多体运动边界与复杂气动流场干扰的多体分离问题。采用嵌套网格技术,对流体力学方程组及六自由度刚体动力学方程组进行了耦合求解,重点研究了子母弹的结构特征(包括子弹的数目、有无弹翼及排列方式等)对分离过程中流场特性的影响。结果表明:嵌套网格技术能较好地模拟超声速子母弹分离过程中的复杂扰动流场;子弹体数目越多,弹体间的气动耦合效应越强;尾翼的存在能够极大提高弹体的静稳定性;相较于"圆周型"分布,"中心圆周型"分布趋于加快弹体间的分离。     

再入飞行器高超声速三维粘性流场的数值模拟

用数值方法研究大钝头倒锥体再入飞行器高超车速完全气体有攻角绕流的流场特性。控制方程为完全气体三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,数值方法采用激波捕捉的ＮＮＤ隐式算法,计算区域为含底部的全流场。文中给出了流场特性和气功力系数,并与实验结果进行了对比。     

AUSMPW+格式在高超声速热化学非平衡流数值模拟中的应用

将AUSMPW+格式应用到高超声速热化学非平衡流场的数值模拟中。为提高精度,采用了三阶MUSCL插值方法。与LU SGS方法结合,提高了单步计算效率和收敛性。采用热化学非平衡十一组元气体模型求解了非定常轴对称Navier Stokes方程组,得到了收敛结果。数值模拟结果与文献结果进行了对比,并在弹道靶中进行了钢质圆球的实验验证。计算结果与文献、实验的对比说明,AUSMPW+格式可以在热化学非平衡流的数值模拟中精确地捕捉到强弓形激波,得到准确的空气动力系数。     

高超音速粘性流场数值计算中差分格式的对比分析

本文通过三维高超音速粘性流场的求解对三种常用的差分格式： Ｈａｒｔｅｎ 的 Ｔ Ｖ Ｄ 格式、修正后的 Ｏｓｈｅｒ Ｃｈａｋｒａｖａｒｔｈｙ 的 Ｔ Ｖ Ｄ 格式以及 Ｊａｎｅｓｏｎ 的当地极值递减 （ Ｌ Ｅ Ｄ） 格式对激波的分辨率、边界层计算的准确性和格式的计算效率等三个方面进行了较细致的对比分析， 为三维高超音速复杂流场的数值模拟打下基础。     

交叉激波湍流边界层干扰问题数值模拟

基于BANS方程,采用有限体积法离散计算区域,对流项使用Harten TVD格式进行流场求解,湍流部分应用Menter SST两方程湍流模型,建立了隐式LU-SGS方法的三维CFD计算程序,数值模拟了来流马赫数4.96条件下的18°×18°对称双楔产生的交叉激波/湍流边界层干扰问题,计算了不同剖面沿展向压强分布以及流场的详细结构,并将计算结果与相关文献的实验与计算结果进行了比较.     

三维高超音速干扰流场数值模拟与实验验证

采用Osher Chakraverthy的TVD格式、Baldwin Lomax湍流模型和LU SSOR隐式方法求解了完全NS方程 ,数值模拟了三维高超音速绕流与横向喷流干扰流场。并在高超音速炮风洞中开展了喷流实验研究 ,对该喷流流场的数值计算结果进行实验验证     

二元混压式高超声速进气道流场数值模拟

使用有限差分法,TVD二阶迎风格式,Baldwin-Lomax湍流模型,对二元混压式高超声速进气道设计状态和非设计状态的粘性流场和性能进行了数值模拟,并分析了出口反压及隔离段长度对进气道流场的影响。结果表明隔离段长高比对进气道出口参数影响较大;隔离段出口反压不能太高,否则会引起进气道不起动。     

高超声速流动计算中LU-SGS隐式算法的应用研究

在高超声速条件下,对原始LU-SGS格式及其改进方法的收敛速度做了深入地比较分析,目的是进一步更好地将LU-SGS算法用于工程上复杂外形的计算模拟当中。二维圆柱,三维钝锥及空天飞机算例的结果表明：(i)对于高超声速粘性流动的计算,粘性项应进行隐式处理;(ii)BLU-SGS方法给出的内迭代方式的收敛性优于DP-LUR方法所给出的内迭代方式;(iii)LU-SGS算法中雅克比系数矩阵的计算方式对计算量及收敛性影响较大,若采用精确的矩阵形式则在流动无分离情况下能取得快速收敛的效果,而在含有流动分离的情况因受稳定性的影响精确的矩阵形式的收敛表现不及对角近似形式。     

高超声速飞行器滑翔制导方法综述

阐明高超声速飞行器滑翔制导的基本问题,分析滑翔制导过程面临的复杂多约束、机动任务要求、参数扰动等研究难点;分别就国内外标准轨迹制导方法和预测-校正制导方法相关研究现状展开综述,指出了这两类方法中存在的问题。在此基础上,提出高超声速飞行器滑翔制导研究中亟待解决的关键问题,并指出未来滑翔制导方法的研究热点。     

双原子分子气体的定常隐式全流域多尺度算法

巨大的计算资源需求极大地阻碍了统一气体动理学格式的应用。采用宏观预估技术,基于Boltzmann-Rykov模型方程发展全流域适用的保守恒定常隐式算法,协同求解宏观方程和微观方程以加速收敛。在单元界面,通过模型方程特征差分解构造简单高效的多尺度数值通量,并结合非均匀非结构速度空间和速度空间自适应技术进一步降低计算需求、提升计算效率。超声速和高超声速平板绕流和圆球绕流的数值结果验证了算法的准确性与高效性。结果表明,算法能够准确求解二维和三维双原子气体多尺度流动问题,且相比于显式离散统一气体动理学格式可加速一个量级。     

碳-酚醛材料烧蚀热解对再入流场特性影响的数值计算

数值模拟存在碳-酚醛材料烧蚀的高超声速再入流场,分析烧蚀和热解对流场热化学参数、电子数密度分布等的影响。采用19组元双温度的热化学模型,耦合热化学非平衡流Navier-Stokes方程组和烧蚀壁面边界条件,进行定常烧蚀流场求解;通过对比无烧蚀、非催化和辐射平衡壁温条件下的流场分析烧蚀的影响;讨论了壁面处碳-酚醛材料热解产物化学组成的确定方法,研究了不同热解率的影响。以RAM-C球锥的两个典型飞行条件(速度7.65km/s、高度61km和71km)为代表的研究表明:最主要的烧蚀热解产物是CO、H2、H,烧蚀产物和烧蚀的影响均局限于边界层内;烧蚀使原子和离子组元含量下降,当离子组元含量峰值出现在边界层内时,烧蚀使电子数密度峰值下降;随热解率增加烧蚀影响程度增强,烧蚀在后身区影响范围大于头部区,随飞行高度增加烧蚀影响范围扩大。