非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。     

COIL化学非平衡流动的一种解耦方法研究

提出了一种计算氧碘化学激光器化学非平衡流动的解耦方法.求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分混合反应流场.空间离散格式为AUSM up格式,用四步龙格-库塔方法作显式时间推进.假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型.使用一种"点隐式"方法来处理化学源项的刚性问题.计算结果表明本文方法可以有效解决计算中的刚性问题,与现有隐式算法相比减少了大量的矩阵运算,公式推导简单,易于实现和改进.     

二维热化学非平衡流动的非结构网格DSMC方法及其应用

研究了二维热化学非平衡流动非结构网格DSMC方法实现的过程。提出了一种新型的高效搜索算法,该算法不仅可以跟踪模拟分子在网格之间的迁移,而且在搜索过程中可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射非确定论判据。为加快流场的时间发展历程,设计了适合DSMC方法的动态局部时间步长技术,使其可以应用到定常和非定常流场的计算。利用Fortran90的动态分配内存技术编制了计算程序。最后对过渡流域高超声速圆柱绕流进行了数值试验,计算结果初步验证了该算法的可行性。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

化学非平衡高超声速电离流动数值研究

本文对非平衡化学区应高超声速绕流零攻角轴对称钝头体电离空气粘性流动进行了研究，通过带有化学源项的Ｎａｖｉｃｒ－Ｓｔｏｋｅｓ方程对流动作数值模拟，总的连续方程由单个组元守恒方程代替。高温电离空气中存在７个主要的组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ，和ｅ￣－。化学动力学模型由７个组元和这些组元之间６个可能的反应来描述。发展了一种数值方法以求解这组完全耦合的与时间有关的偏微分方程组。利用一类新的有限差分方法来求解这组以守恒形式表示的控制方程。这是一种以ＴＶＤ格式来离散粘性流通量的方法。在数值实验中采用了ＴＶＤＹｅｅ－Ｈａｒｔｅｎ格式，对全部算例都是用激波捕捉法。关于化学组元，我们试图处理完全催化壁、非催化壁和有限催化壁等三种不同的表面边界层条件。计算是对ＲＡＭＣ－Ⅱ试验条件进行的，得出了７１ｋｍ高空情形的典型流场数值解。对表面压力，Ｓｔａｎｔｏｎ数和热流率给出了化学反应气体结果与完全气体结果的比较。计算的流场电子密度与ＲＡＭＣ－Ⅱ实验数据进行了比较，计算结果与实验数据符合是相当好的。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

再入飞行器化学非平衡流动的隐式向量算法研究

在YH－2计算机上对再入飞行器化学非平衡流流场的NND隐式程序进行向量化计算，在对程序作了重构优化后，向量加速比达3．4。这是对化学非平衡流N－S方程隐式程序进行向量化计算的成功尝试，也为进一步实现程序的并行执行奠定了良好的基础。     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

非平衡流解耦方法及其计算激波诱导燃烧的应用验证

介绍并发展了一种改进的化学非平衡流动解耦方法.将流动控制方程与化学反应生成源项解耦处理,组分对流项的矢通量分裂基于流动方程密度对流项的分裂形式给出.对流项采用五阶WENO格式求解,化学反应源项的刚性采用简化的隐式方法进行处理.将该方法应用于预混H₂/O₂和H₂/Air激波诱导燃烧计算,得到的定常及非定常燃烧过程与实验观测相符,振荡频率与实验测量较吻合,结果表明该方法应用于多组分多步反应流计算是可行的.      

高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

普朗特-迈耶流中湿空气非平衡凝结影响的研究

普朗特-迈耶流是飞行器设计过程中一种重要的基本流动.本文发展了求解湿空气非平衡凝结流动的数值方法,对凸角分别为90°和40°的湿空气来流条件下的普朗特-迈耶流动进行了数值模拟,并且与试验结果进行了对比.90°凸角时,膨胀扇内静压分布的计算值与试验结果吻合良好.40°凸角时,膨胀扇内密度的计算值与试验值吻合良好;随着来流相对湿度的增加,凝结起始位置向上游移动而最大成核率减小.膨胀扇内的湿空气非平衡凝结显著改变了压力、马赫数、温度等参数的分布.     

高超声速非平衡流场多个振动温度模型的数值研究

利用数值求解热化学非平衡Navier-Stokes方程的方法,建立了考虑一个和多个振动温度的高超声速热化学非平衡流场的CFD计算程序,并对半球模型和球锥模型进行了数值计算,研究结果表明:(1)各个分子组分的振动温度分布于考虑一个振动温度模型的振动温度周围,即一个振动温度模型的振动温度是多个振动温度模型的各个分子组分振动温度的一种平均值;(2)在模型身部,振动温度的峰值高于平动温度的峰值;在模型近尾出现振动温度冻结并高于平动温度。     

再入体实验模型热化学非平衡绕流流场的数值分析

利用数值求解一个和多个振动温度热化学非平衡Navier-Stokes方程的CFD计算程序,对爆轰风洞球锥试验模型的热化学非平衡绕流流场进行了数值模拟,分析了分子组分振动温度在全流场(头身部流场和底部流场)中的分布规律.研究结果表明:(1)计算的压力、电子数密度以及流场的光辐射数据与试验数据符合较好;(2)在再入体身部,大多数分子的振动温度峰值大大高于平动温度的峰值;在再入体近尾出现振动温度冻结并大大高于平动温度的现象;CO2的两个振动态的振动温度分布非常接近平动温度分布.     

超声速空间发展燃烧反应剪切层放热效应分析

为探索化学反应对混合燃烧的影响机理,本文从流体力学基本方程出发,采用直接数值模拟方法研究了二维7组元/8反应化学非平衡剪切层的流动,讨论了数值上考虑化学反应放热对流动混合的影响,并分析其机理.研究发现:在相同的基频扰动促发下,考虑化学反应时混合层厚度变薄,且混合层的厚度会随着燃料组分的增大先变小后增大,涡量厚度发展至饱和的位置逐渐靠后;反应放热和组分含量影响湍流切应力和湍动能的变化,涡量厚度随着湍动能以及湍流切应力规律性变化,随着燃料组分的增加,湍动能、湍流切应力及混合层厚度非单调变化,存在极值点.     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。