关于三维Navier－Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵［Ａ］和［Ｂ］，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程计算（即Ｍ∞＝３．０，α＝０°和设计工况Ｍ∞＝２．６５，α＝０°），获得了满意的全场结果；对于Ｍ∞＝２．６５的设计工况，同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的Ｎ－Ｓ求解能在普通微机上进行。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

三维湍流高速进气道内外流场的高效高分辨率解

本文提出一种高效、高分辨率求解三维复杂流场的隐式算法，并成功地用于求解高速进气道内外流场的三维、可压缩、雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，获得了三个工况五种流场的三维数值解；与有关实验数据比较，结果令人满意。这种算法由Ｊａｍｅｓｏｎ－Ｔｕｒｋｅｌ的ＬＵ分解法和张涵信ＮＮＤ格式的数值通量耦合而成，湍流采用ＢａｌｄｗｉｎＬｏｍａｘ模型。大量数值计算表明，该算法能捕捉到复杂的激波系，有较好的激波     

边界层分离数值解

用N-S方程直接求解分离流的工作近几年已有不少。但这种计算比较费机时,需要存储最也大。在实际应用中,边界层理论仍然是研究分离流的一种好的近似方法。本文对二维、定常、不可压层流的分离流动进行了数值分析。主要内容有:(1)从典型的以流函数表示的边界层方程出发,对分离流进行数值计算。Keller和Cebeci曾用此方程求解过附体流动的反问题。在具有相似性解的条件下,这个方程简化为著名的Falkner-Skan方程,在初始剖面,这个方程简化为Blasius方程。(2)给     

旋涡内流计算

本文采用涡量-向量势表示的不可压缩三维流动的基本方程及其边界条件。针对小曲率弯管流,忽略主流方向的扩散项而使基本方程抛物化并利用流量守恒条件确定轴向压力分布。为避免横向速度向量散度不为零所带来的麻烦,本文提出两种计算方案。计算结果证明:用于求解三维抛物化方程组的这两种方案都是方便可行的。     

关于底部分离流的研究

底部分离流是气体动力学中的重要研究课题之一,底部流在工程应用上主要是计算底部压力及底部热流,提供流场结构的定性以至定量的描述,由于底部流场本身的复杂性,实验研究,特别是在高超声速范围内的实验数据还不够多,理论分析方法也不够完善,尤其是底部传热的研究远远不及底部压力的研究,这是因为在底部流场结构不太清楚的情况下,受热过程不容易得到合理的处理,一些半经验的方法均有很大的局限性,     

关于三维Navier—Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵[A]和[B]，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Navier-Stokes方程计算（即M∞=3.0,a=0°和设计工况M∞=2.65,a=0°），获得了满意的全场结果；对于M∞=2.65的设计工况。同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的N-S求解能在普通微机上进行。