碳基材料氧化烧蚀的快/慢反应动力学问题

在前期单/双平台问题研究的基础上[1]，就学术上和工程应用中关注的碳在空气中燃烧的"快反应"和"慢反应"问题进行了深入分析，研究发现被广泛使用50余年的"慢反应"并不存在，而被弃置的"快反应"则真实存在，且具有重要的应用价值。采用"快反应"动力学数据，同时考虑CO、CO₂两个烧蚀产物，得到的无量纲烧蚀速率随温度的变化曲线存在两个平台，其中，温度稍低情况下出现的第一平台对应的主要烧蚀产物为CO₂，温度稍高情况下出现的第二平台对应的主要烧蚀产物为CO，且第一平台对应的无量纲烧蚀速率恰好是第二平台的1/2。过去常被忽略的CO₂扮演了重要角色，由它产生的第一平台，将以往文献中看似完全独立、毫无关系的"快反应"和"慢反应"曲线建立了联系。理论分析表明:第一平台之前的快速上升段属于从速率控制区到扩散控制区的过渡区，第一平台及其以后的区域都属于扩散控制区（包括两个平台之间的连接线），它是由反应生成物CO与CO₂的分压比δ从0到∞的变化引起的，与表面化学反应动力学条件完全无关。由"双平台"理论得到的从低温到高温、由速率控制区经由过渡区到达扩散控制区的整条烧蚀速率曲线，与实验结果完全吻合。     

三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态

通过引入依赖于密度的物面法向速度变换wr=-1/ρh1h2∫ozh1h2(e)ρ/(e)tdz,描述物理速度空间(ua,va,wa=u,v,w+wr)具有无滑移壁面条件的三维可压缩非定常连续方程可转换成变换速度空间(u,v,w)内具有无滑移条件定常连续方程.因此,采用定常壁面分离的分析方法和结论,再通过变换和研究wr的贡献,给出了三维可压缩非定常壁面分离的判则以及分离线附近的流动形态.研究指出,二维和三维情况下,都出现伴有壁外附着的壁面分离情况.数值模拟证实了理论和结论.     

一种新的间断侦测器及其在DGM中的应用

根据单元交界面左右变量的差别,提出了一种新的间断侦测器构造方法。该间断侦测器的构造原理简单,编程实现容易。针对一维和二维Euler方程,我们将此间断侦测器用于间断Galerkin格式的数值计算中。数值实验表明本文构造的间断侦测器能够准确捕捉到激波的位置,从而只在间断区域引入限制器,在减少计算量、保证光滑区计算精度的同时,提高了激波等强间断的分辨率,能够明显地改善流动精细结构的模拟精度。     

细长圆锥低超声速绕流尾迹流时空演化特性的研究

本文用数值方法研究了在大的迎角范围内，采用局部抛物化流动模型所得到的有关圆锥有迎角绕流尾迹流中周期解的谐分岔、亚谐分岔以及环上Ｈｏｐｆ分岔这三类不稳定性的时空结构演化与非线性相互作用。分析得到了两种不可通约频率（ｆ１，ｆ２）的强非线性相互作用以及由此而产生的频率拓宽所形成的ｍｆ１±ｎｆ２多级频带结构。研究表明：当第三个不可约频率被激发，不论时间，还是空间方向的发展都会存在混沌。在近尾区，流态随攻角增大的演化历程上，通向混沌的道路和准周期性道路相接近，但包含有带倍周期性的亚谐分岔。在固定迎角的情况下，由近尾到远尾区，通向混沌的道路也接近于准周期道路，但带有阵发性。     

高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

细长圆锥低超声速绕流层迹流时空演化特性的研究

本文用数值方法研究了在大的迎角范围，采用局部抛物化流动模型所得到的有关圆锥有迎角绕流尾迹流中周期的谐分岔，亚谐分岔以及环上Ｈｏｐｆ分岔这三类不稳定性的时空结构演化与非线性相互作用。     

跨声速三角翼旋涡破裂过程的数值模拟

本文利用张涵信提出的时，空二阶精度的隐式混合通量分裂的ＮＮＤ格式，通过求解非定常ＮＳ方程模拟了马赫数Ｍ为０．８、迎角α为２１．８°时某外形三角翼背风区旋涡破裂的发展，演变过程，并对旋涡区域的截面流态和涡轴上奇点附近的空间流动特征进行了分析。     

细长锥体有攻角绕流的非对称分岔及其涡结构的数值模拟

本文用ＮＮＤ格式计算了具有复杂涡结构的细长锥体有攻角绕流问题，着重研究了非对称分离的分岔现象以及分岔前后涡结构的演变，数值计算结果表明，大约在ａ／（２β）＝１．０时，Ｎ－Ｓ方程出现分岔解，同时还给出了分岔前后涡结构及物理最分布的特征。     

粒子侵蚀的相似规律及实验模拟问题

本文探讨了粒子侵蚀的相似规律,给出了粒子侵蚀实验的模拟参数,讨论了相似规律在实验模拟中的应用及实验模拟问题。     

超声速三维粘性分离流动的反扩散数值解法

本文将文献[1,2]的时间分裂的反扩散方法推广到三维粘性分离流动。计算结果表明,这种算法是可行的。