一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

自身场磁等离子体推力器数值仿真流场特性分析

磁等离子体推力器以其推力大、比冲高等特点,成为未来深空探测、星际航行任务首选的电推力器类型,研究推力器内部等离子体流场特性,有利于解释推力器出现的物理现象。针对特定自身场磁等离子体推力器,建立磁流体模型,使用TVD Lax-Friedrich格式以及ADI方法对推力器内部及羽流进行数值求解,得到等离子体流场及等离子体参数分布情况,仿真结果显示在高电流工况下,等离子体羽流更加集中,轴向加速效果更加显著,但高电流模式下阴极温度较高,不利于阴极的使用寿命。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。     

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式,具有十分优异的分辨率特性。但在求解多维方程组时,HCW-R格式需要求解块状三对角方程组,因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P,构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式,其可沿迎风方向、由边界处向内推进求解,避免了处理三对角或块状三对角方程组,从而其计算代价与显式格式无异。虽然就分辨率而言,HCW-E稍逊于HCW-R,但前者的计算效率要显著高于后者。因此,当花费相同的计算代价,HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况,开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质,应用合适的高阶精度格式对各项分别进行离散求解,并采用经典几何声学射线法计算声爆传播路径,实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟,验证了本文高阶离散方法的可靠性。进一步的结果分析表明,采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性,在同等网格数目上可获得更高分辨率的计算结果,且其网格收敛性远优于传统二阶精度格式,同时在计算效率上也有一定优势;在Burgers方程的各效应项中,非线性项的影响更为明显,因而采用高阶离散格式时其优势更为突出;热黏性吸收项对数值结果影响较小,采用高阶离散格式对计算精度的提升效果也不明显,实际计算中仍可以采用传统离散格式,甚至忽略该项的贡献。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

一种改进的紧致WENO混合格式

给出了一种用于激波捕捉计算的守恒型混合紧致WENO格式.混合格式可视为五阶WENO格式和五阶守恒紧致格式的加权平均.在格式构造中,本文采用了间断分辨率更好的WENO权因子计算方法以及新的子格式权因子计算方法.对于Euler方程,仅对混合格式中WENO格式部分做特征投影处理,不仅获得了特征型WENO格式良好的间断分辨率,同时又保证了计算效率.数值实验验证了本文混合格式的高时间效率和高流场分辨率.     

超燃进气道分离流动的凸起物控制机理研究

双模态超燃冲压发动机由于压力扰动可能发生不起动现象,造成推力严重下降,对飞行稳定性与飞行安全具有很强的破坏性。不起动初始阶段主要受到激波与边界层相互作用引起的流动分离影响,采用5阶特征型WE-NO格式与3阶TVD型Runge-Kutta格式的高精度数值方法,求解三维Navier-Stokes方程,研究了利用凸起物作为涡发生器的被动控制方法,及其对高超声速流动分离现象的控制效果。结果表明高精度数值格式能够捕捉到复杂精细的流动分离结构,总结了摩阻、压力等在分离再附位置的变化规律;发现凸起物通过诱导形成局部流向涡进行流动控制,能够改变压力分布,减弱分离强度,影响分离结构。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...