双重加权实质无波动激波捕捉格式的改进和应用

本文重新研究了双重加权实质无波动(DWENO)激波捕捉格式的构造过程，对模板光滑系数的定义，以及两重权值的计算方法进行了改进，并通过一维Riemann问题测试了这些改进，确定了参数的选择范围。最后，本文将改进后的DWENO(r=3)格式进一步应用于一维、二维典型问题的数值模拟，并比较了DWENO格式和WENO格式的精度和分辨率。     

飞行器亚跨超声速流气动力并行计算研究

隐格式并行求解是计算流体力学并行算法研究的难点,本文针对分布式存储并行处理系统,在异构编程PVM环境下,采用"局部隐"方法,使用发送/接收数据结构和自动剖分等技术,有效地实现了对飞行器亚跨超声速流气动力数值模拟串行程序的并行化改造,保证了并行计算结果与串行结果的一致性.该方法具有简单实用、并行效率高以及并行程序可移植性好和应用前景广泛等特点.     

Jameson/TVD格式应用旋翼翼型绕流分析

以一种常见的旋翼翼型为例 ,探讨了两种旋翼桨叶翼型绕流数值模拟方法 ,即 Jameson有限体积中心格式和有限体积通量修正 Jameson/ TVD混合格式。这两种格式的方法都采用有限体积法进行离散 ,并用龙格—库塔时间推进法求解。通过算例分析 ,比较了这两种方法的特点 ,为进一步有效地依靠求解湍流 Navier-Stokes方程来模拟完整的旋翼流场、旋翼阻力与功率奠定了基础     

高超声速飞行器外形热流密度分布计算的高精度方法研究

为提高高超声速条件下飞行器表面热环境的计算精度,发展了高精度算法WCNS-E-5,并开展了热流密度分布的数值计算研究.该算法采用的是五阶精度加权显式非线性格式,结合四阶精度的二阶偏导数差分近似、四阶精度的边界格式,同时还对网格偏导数及温度梯度也离散为四阶精度.通过求解三维Navier-Stokes方程,数值研究了高超声速飞行器外形如钝锥和双椭球体的物面热流密度分布.模拟结果表明,这种保证全流场高精度的WCNS-E-5得到的流场图像清晰、真实.与通常的二阶算法MUSCL对比,WCNS-E-5具有更高的流线分辨率,获得的热流密度更接近于实验测得分布.     

非定常喷流和尾腔流动的高分辨率数值模拟

导弹在运载器尾腔内的点火分离，涉及高温、高压气体在喷管和运载器尾腔内的三维非定常流动。采用分步显式和隐式ＮＮＤ格式对轴对称和三维多喷管非定常亚、跨、超混合喷流流场进行数值模拟，单喷管导弹的喷管流动和尾腔流动使用整体轴对称流动模型；对于多喷管导弹，喷管流动使用轴对称流动模型，尾腔流动为三维流动模型，喷管与尾腔流动间由边界参数进行耦合。算例表明，该方法能较好地捕获喷管和尾腔流动中的三维运动激波，计算结     

CPL编组18RMK语句浅析

民航总局于2006年3月10日发布了新的行业标准MT/T4007-2006《民用航空飞行动态固定电报格式》,并于2006年6月1日起实施,用以取代先前使用的旧版本MT/T4007-1998。其中,主要变化之一是取消了原来国内AFTN电报所使用的返航报RTN、备降报ALN,取而代之的是现行飞行变更报CPL。     

有关有限差分高精度格式两个应用问题的讨论

激波装配法结合有限差分高精度格式是发展全场一致高精度算法的一种途径。在对流场内的间断全部进行装配后,对高精度格式的应用研究成为发展本算法的主要研究问题。本文将有限差分的高精度格式应用于贴体坐标系时发现,对均匀流场,高精度格式因不满足几何守恒律而产生的数值误差比一阶迎风格式大,初步分析认为是由于高精度格式所用的模板比一阶格式更宽,涉及的网格点数更多,从而引入了更多的误差。而作者提出的基于离散等价方程的相容性算法可消除这一误差。此外,本文在利用激波捕捉法求解正方形均匀网格上的正激波运动问题时发现因通量分裂格式的使用,在激波处会产生随着特征线传播的非物理波动,这一波动在激波与网格不完全匹配时表现为多维波动相互干扰的虚假"数值湍流"现象,高精度格式的高分辨率特性使得这一现象更加明显。这是因为激波捕捉法假设激波为空间连续函数,用于包含激波的流场时必然得到数值解表示的过渡区,导致可信度评估困难,使用激波装配法可以避免这一问题。     

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

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

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

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

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

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