求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

钝体头部边界层特征谱分析和转捩途径预测

利用线性稳定性理论,分别对钝头体高超声速绕流边界层中的亚声速区域、跨声速区域和超声速区域的特征谱进行了系统分析.结果表明,根据线性稳定性理论,所分析的各速度范围流场,对所有单频波的扰动都是稳定的.但随速度的增加,稳定波的衰减率越来越低;随激波到物面距离的增加,s族扰动波相速度的实部向无量纲速度1靠拢,说明头部主激波的作用在一定意义上与平面槽道流动的上壁面相似;研究工作还表明,流场中存在着低衰减率扰动波,且无论那个速度范围内的流场,其最低衰减率扰动波的相速度实部都在0.7左右.根据钝头体特征谱分析,本文可能发生的"bypass"型转捩机理进行了理论分析,提出了几条可能的"bypass"型转捩途径.     

火箭燃气喷流的数值和实验研究

通过数值计算和实验,研究了２１和９火箭发动机燃气喷流流场。数值计算采用激波高分辨率、二阶精度ＴＶＤ格式,得到包括射流边界、桶形激波、反射激波、Ｍａｃｈ盘和滑移线的清晰的流谱。实验测试应用Ｆ－Ｐ激光干涉仪和自动测压系统,得到流场Ｆ－Ｐ干涉照片和压强分布。计算结果与实验结果符合较好     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

块结构的MUSCL算法求解三维进气道流场

本文将ＭＵＳＣＬ格式扩展为三维算法,并用于湍流流场的数值模拟。通过求解超声双斜楔压缩面流动,得到了含有激波、嵌入激波、剪切流和湍流边界层等流动特征的复杂流动结构,数值解与实验数据吻合得很好。文中将此算法与块结构网格技术相结合,使之能求解具有三维复杂几何形状的流场。并应用该算法模拟了Ｆ－１６战斗机进气道在５种典型飞行状态下的三维湍流流场,计算结果清楚地揭示了这些情况下的三维流动现象,得出了进气道的性     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

带座舱飞船高超声速再入非平衡流场的数值研究

本文针对带座舱飞船高超声速再入大气层过程中存在的严重气动加热现象 ,利用简单隐式TVD差分格式和激波捕捉法 ,数值求解三维化学非平衡Navier Stokes方程 ,其中化学模型是 7组分 1 5个化学反应的空气化学模型 ,对带座舱飞船再入高度为 70km和速度为6 7km/s的化学非平衡流场进行了数值模拟 ,给出了带座舱飞船再入迎角为 0°,1 0°和 2 0°情况下的壁面热流、表面压力和电子数密度等参数分布     

双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

低信噪比下分组交织器识别

针对现有的分组交织器识别算法计算复杂高且容错性差缺点,从分组交织后的同步码分布规律出发,提出了一种新的识别算法。首先,利用数据矩阵统计特性,给出了在任意矩阵列数下,同步码和随机业务数据位置上的概率密度分布函数,基于最小错误判决准则,设定了同步码检测门限,同时基于3倍标准差准则,求解出稳健的交织周期识别门限;其次,分析了数据矩阵中每一行与每一列累积量之间的对应关系,提出了一种快速交织周期遍历方法,使得数据矩阵的构建次数大大减少;最后,总结了4个分组交织后同步码分布规律,通过遍历同步码序列,利用同步码之间的位置关系,实现交织同步位置、分组交织列与交织行参数快速识别。仿真结果表明:所提算法具有较强的低信噪比容错性,在信噪比为-6 dB条件下,参数识别率能够达到98%以上,同时与现有的算法相比,其性能提升4~10 dB且计算效率明显提高。      

一种改进的MUSCL格式

给出了一种改进的二阶精度的ＭＵＳＣＬ格式，并用来求解Ｅｕｌｅｒ方程，改进后的ＭＵＳＣＬ格式在线化插值步中，用Ｖａｎ Ｌｅｅｒ的插值斜率对特片变量进行插值，并对与熵波、剪切波对应的特片变量使用了小耗散的插值斜率，在迎风步采用Ｒｏｅ的近似黎曼解求得迎风通量。比较一维和二维的计算结果均可以看出，改进后的ＭＵＳＣＬ格式在消除数值振荡的同时，明显提高了流场间断的分辨率。