柱面气相散心爆轰波胞格演化规律研究

本文采用二阶精度NND格式,应用改进的二阶段化学反应模型,通过求解二维Euler方程对柱面气相散心爆轰波胞格演化过程进行了数值模拟.计算结果表明在传播过程中,空间尺度的扩张导致了散心爆轰波后气流的自然膨胀,使得多波结构的爆轰阵面呈现出显著的胞格自组织特性.根据计算结果与理论分析,本文归纳了五种胞格演化模式,分别命名为内凹波阵面会聚、波阵面扭结、褶皱波面失稳、胞格自合并和三波点滑移,并定义了各种模式的物理特征,分析了其相关的演化机制和规律.     

高阶精度线性耗散紧致格式的渐近稳定性

文献〔２〕提出的高阶精度线性耗散紧致格式,（ＤＣＳ）能较好地抑制高波数振荡。本文进一步研究表明：内点耗散参数α能够控制耗散紧致格式ＤＣＳ３和ＤＣＳ５的渐近稳定性,边界格式参数α１和α对以上格式的渐近稳定性的影响是次要的,并给出了 些参数的取值范围。通过平面可压缩Ｃｏｕｅｔｔｅ流的稳定性数值模拟表明,当格式ＤＣＳ５满足渐近稳定性要求且具有耗散时,计算结果与线性稳定性分析的结果十分吻合,否则,会出现较     

用BGK格式对平面可压缩自由混合剪切流动的数值模拟

本文介绍了ＢＧＫ－Ｂｏｌｔｚｍａｎｎ格式的构造与用有限体积法计算网络内宏观量的方法。文中还介绍了可压缩平面混合剪切层计算域在不同出口马赫数下的边界条件。然后将上述方法应用于四种不同情况的可压缩平面混合剪切层的数值模拟。还详细介绍了计算结果,并与其他方法所得的的结果以及实验数据进行比较。表明本文的结果较好,方法是有效的。     

残量分布格式初步研究与应用

介绍了一种新的流场计算数值格式——残量分布格式。它是一种多维迎风格式,适于三角形网格和四面体网格。给出了该格式的标量方程和Euler方程。算例显示:该格式能处理稳态与非稳态问题,计算精度高,有较强的捕捉激波能力,且结构紧凑,易实现编码。     

两种湍流模型及其可压缩修正在高超压缩拐角中的评估(英文)

采用在高阶精度差分格式对SA湍流模型和SST湍流模型及相应的可压缩修正模型在高超拐角流中进行了评估。可压缩修正方式考虑了密度梯度、压力膨胀和湍流马赫数等方法。为了减小数值误差与模型误差之间的混淆,控制方程的对流项采用了5阶精度的加权紧致非线性格式(WCNS-E-5),粘性项采用了一种半结点/结点交错的4阶中心格式。通过对马赫数为9.22的15度拐角和34度拐角湍流的模拟,考察了原始湍流模型及其修正模型的效果。计算表明:原始SST模型对高超拐角湍流的预测比原始SA模型准确,这种准确主要体现在对分离区预测、再附点附近压力和热流峰值预测上。通过混合采用Catris和Shur等的方法对 SA模型进行可压缩修正可以大大改进模拟效果。在SST模型的可压缩修正方法中,Catris的修正方法最好;考虑压力膨胀修正后得到的分离区远远偏大,比本原始SST模型更差;考虑湍流马赫数的修正方法得到的分离区偏小。本文还给了部分低阶格式的计算结果,高阶格式与低阶格式相比,对分离区大小、再附点附近的压力和热流峰值等的预测准度有所改进。     

基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

基于GPU和隐式格式的CFD并行计算方法

从图形处理器(GPU)架构特点出发,提出了基于数据并行的隐式计算流体力学(CFD)求解方法,空间离散格式采用迎风Roe格式,计算网格适用于结构和非结构网格。采用统一计算设备架构(CUDA)技术实现了GPU上的隐式CFD并行计算。分别在Intel Core2 Quad3.0GHzCPU和NVIDIAGTX280 GPU上进行了计算,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的计算加速比达到28倍,采用隐式格式计算加速比达到了28.7倍,同时计算加速比随计算规模的增加而增加。计算结果和实验结果较为吻合。     

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

用四阶Runge-Kuta法和TVD格式求解两维跨超音流场

利用坐标转换给出两维任意曲线坐标系下雅可比矩阵特征值、特征向量及Harten型TVD格式的详细表达式, 并结合四阶Runge-Kuta法发展了一种新的、有效的求解超音、跨音流场的算法。数值计算表明, 本文的方法可有效地求解跨音流场。