流过泰氟隆烧蚀表面的化学非平衡粘性激波层数值研究

本文对有泰氟隆烧蚀的化学非平衡粘性激波层流场进行了数值模拟。应用空间推进与总体迭代相结合的求解方法，采用１９组元、２９种反应的空气—泰氟隆化学反应系统，对化学非平衡粘性激波层流场与泰氟隆壁面烧蚀传热过程进行耦合求解，研究了平衡催化壁和非催化壁条件对粘性激波层烧蚀流场的影响。     

高空滑移条件对化学非平衡粘性流动的影响

本文用数值方法计算了不同高度滑移边界条件对钝体绕流流场特性的影响。采用化学非平衡粘性激波层控制方程,边界条件含壁面滑移、激波滑移和无滑移条件。计算结果表明,考虑滑移条件后,激波层较厚,激波处的电子数密度为一有限值。文中给出了计算结果,并与其它文献的结果作了比较。     

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

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

11 组元化学反应气体粘性激波层钝体绕流数值计算

本文采用１１组元化学模型对双曲体粘性激波层化学非平衡绕流流场进行了数值计算,给出了压力、温度,Ｎ＋２、Ｏ＋２、Ｎ＋、Ｏ＋和ＮＯ＋摩尔浓度及ｅ－数密度在驻点的分布,并与７组元、５组元的计算结果作了比较。     

高超音速化学平衡流粘性激波层数值计算

本文给出了层流化学平衡粘性激波层方程对于细长解析体的数值解。这些求解结果是通过空间步进方法得到的,并且与参考文献里的计算结果作了比较,以估价该方法的精确度。结果表明,这种控制方程完全耦合求解的方法能够给出相当精确和稳定的结果,这些结果与参考文献的计算结果相当吻合。     

气体分子碰撞传能及反应的DSMC 仿真

研究了气体分子的非弹性碰撞机理,给出了分子内自由度激发、反应的几率选择模式及分子内自由度能量交换模式。以钝体驻点线流场为仿真实例,比较了ＤＳＭＣ仿真结果与粘性激波层方程计算结果。二者在较低高度时的一致表明ＤＳＭＣ仿真中采用的模式是合理的。     

激波对密度间断层干扰及界面效应的数值模拟

运动激波在密度间断层传播时 ,引起密度分界面的变形 ,相应激波形状发生改变。以高密度层为例 ,采用TVD格式模拟这一过程 ,并与半解析解比较分析 ,初步讨论了激波对高密度层干扰作用引起高密度层的上扬和高密度层内激波与壁面的反射现象     

两道同侧扫掠激波作用下形成的三维平板边界层

为探索一种减弱扫掠激波/边界层干扰强度的侧压方式,采用数值计算方法研究两道同侧扫掠激波作用下形成的三维平板边界层。对三维边界层的流动机理进行研究,并与相同气流偏转角的一道扫掠激波作用下形成的三维边界层进行定量比较。研究表明,存在两道同侧扫掠激波时,第一道扫掠激波使侧壁附近边界层变薄,进而使第二道扫掠激波与再附形成的边界层的相互作用减弱,从而总的效果是两道扫掠激波与边界层的相互作用强度比一道扫掠激波的弱;两道扫掠激波的气流总偏转角与一道扫掠激波的相同,相交后汇聚成的激波强度与一道扫掠激波的基本相同,之前被两道扫掠激波"扫"起来的边界层在汇聚激波作用下还呈现锥形流动,分离线和来流的夹角也与一道扫掠激波的基本相同。     

超声速扩压器中激波串结构的数值模拟

通过求解由BL湍流模型封闭的二维、轴对称及三维雷诺平均N S方程 ,数值模拟了等截面超声速扩压器中由激波 /附面层干扰诱导的复杂流场 ,比较了二维直管、圆截面直管及三维矩形截面直管中的流场特性、激波串长度及压强恢复程度。在来流马赫数为 3 0的二维直管计算中 ,采用四步Runge Kutta显式方法数值仿真了激波串自激振荡过程 ,并与实验结果作了对比分析     

基于激波场张量捕捉的超声速弹丸跟踪

超声速弹丸产生的激波蕴含丰富的运动状态信息。在分析弹丸激波形成机理和建立激波形成前后参数关系的基础上,建立空间张量的激波捕捉三维Euler模型和声信号融合算法,并进行空间谱辨识实验;为实现反演跟踪,结合外场试验背景,提出动态g-h-k滤波跟踪方法,验证激波声场反演跟踪的有效性,为探索弹丸激波场进行弹丸跟踪提供了新方法。     

等截面隔离段中激波串结构的数值模拟

采用二阶TVD格式的有限体积法偶合Baldwin Lomax代数湍流模型求解雷诺平均Navier stokes方程 ,数值模拟了二维矩形等截面隔离段中的流动现象 ,很好地模拟出由激波 /附面层干扰所形成的复杂的激波串流场结构。计算结果与国外有关文献的结果进行了比较。     

预处理方法在全速度流场数值模拟中的应用

基于预处理方法,采用有限差分方法,LU-SGS隐式求解了二维可压Navier-Stokes方程,数值模拟了无粘突包、粘性空腔、化学非平衡流场,并与试验结果进行了对比。计算结果表明,所采用的数值模拟方法能够求解低速及亚、跨、超全速度流场和化学非平衡流场,所采用的预处理方法在亚声速范围内保持数值精度的前提下,加速了收敛。     

全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢.原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩...     

流动相似律研究

有别于经典的量纲分析方法和对流动控制方程及其定解条件的无量纲化分析方法, 论文给出了一种分析流动相似律的新思路, 亦即将复杂的流动过程分解为若干个子过程, 然后分析每一个子过程中反映其本质性态的特征量。依据这一思路导出的高超声速热化学非平衡流动相似律与文献采用经典分析方法得到的流动相似律一致。     

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

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

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

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