高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

轴流压缩机粘性尾迹噪声研究

本文研究了微弯机翼在正弦阵风作用下的脉动升力,并导出计算公式,分析了动、静叶栅几何及气动参数对转子脉动升力及噪声级的影响。     

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

跨声速粘性流绕振荡翼型的非定常计算

本文采用时间和穴是二阶精度的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式及ｑ－ω二方程微分模式，结合网格自适应技术，数值模拟Ｎ－Ｓ方程，计算了跨声速下的翼型非定常运动，包括俯仰，浮沉和前后平移振荡。结果表明，压力分布和气动系数与实验基本符合，微分模式和自适应网格能够显著提高激波和边界层计算精度。     

超声速绕凹角流动的压力分裂形式简化N—S方程数值解

本文从压力分裂形式的简化Ｎ－Ｓ方程出发，用Ｒｕｂｉｎ格式对超声速粘性流动进行了数值求解，给出了超声速绕凹角附着流及带有小分离气泡流动的结果，与其它数值方法所得结果相比，符合较好。     

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

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

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

硅基复合材料粘性系数的参数辨识

硅基复合材料是目前国内比较成熟、综合性能较好的烧蚀型热防护材料之一,高温粘性系数是这类材料烧蚀性能的重要表征参数.本文利用系统识别原理,给出这类复合材料粘性系数的参数辨识方法及其应用实例.     

粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。