论文部分内容阅读
目前,流体力学中非线性偏微分方程的求解有各式各样的数值解法。基于单元的数值方法在求解流体力学问题时均需要利用单元进行计算区域的离散,这不便于求解和应用。无单元Galerkin(Element-Free Galerkin,EFG)法作为一种有巨大研究及应用价值的数值方法,其采用的是基于离散点的近似,不需要单元信息,不仅便于处理复杂边界,而且计算准确,因此不少学者将其应用于流体力学问题的求解。近年来EFG法在流体力学领域的应用仍然处于发展阶段,而且主要用于求解二维流场问题,但是很多实际流动问题往往不能简化为二维流动问题进行分析,故本文以低雷诺数下三维不可压缩粘性流体流动为研究对象,以EFG法和计算流体力学为理论基础,编写计算程序对其流动进行数值研究,全文的主要研究工作和所得结论如下: (1)推导出三维Stokes流动的EFG法离散控制方程,以推导的离散方程和EFG法理论为依据编写算法程序对具有精确解的三维Stokes方程进行了求解,验证了所编程序和算法是有效可靠的。将其应用于立体空腔Stokes流动问题的求解时,重点分析了不同速度比和高宽比对流动的影响,数值结果表明随着速度比S值的增加,其内部流态由四个临界值分成五个不同的阶段;随着高宽比A的变化,其内部漩涡也经历了由一个到两个,随后又只有一个,接着发展为多个漩涡的过程。 (2)考虑惯性项的作用,建立了三维定常Navier-Stokes(简称N-S)方程的EFG法离散表达式,编写程序数值研究了180°圆形截面弯管的流动特性,重点探究雷诺数和曲率直径比对其流场和漩涡生成的作用。结果表明,在曲率直径比相同时,入口雷诺数的取值影响着弯管内迪恩漩涡的结构和强度;在雷诺数相同时,理尼漩涡在小的曲率直径比弯管中出现的地方靠近弯曲段前沿,而在大的曲率直径比弯管中则往后一些。另外,弯管内速度和压力的分布受曲率直径比变化的影响更为明显。 (3)考虑到与时间有关的流动情况,采用θ加权法进行时间项的离散,完成三维非定常N-S方程的EFG法离散格式的推导,编写求解程序对并列圆柱在低雷诺数下的绕流展开分析。结果显示,并列圆柱的三维流动特性明显,不同雷诺数下并列圆柱中的尾流形态存在差异。 本文运用EFG法对三维定常Stokes流动、考虑惯性项的弯管定常流动以及考虑时间项的并列圆柱非定常流动展开了数值模拟,计算结果准确反映了其三维流场情况,给工程应用提供了一些有效的指导。